[cedwilli]

First, I would like to express my thanks to everyone who has contributed to solids4Foam, and especially the main developers.

I have been working on a thermal stress model where the material properties are not linear enough to use a single values for Young's modulus and have been trying to see if I can use a piece-wise linear function of stress vs. strain to define the mechanical properties for the thermal mechanical stress model.

I was able to install solids4Foam on an Ubuntu 20.04 LTS computer with Foam-Extend-4.1. I copied and successfully ran the thermoelasticity, hotCylinder tutorial, after adding a keyword entry for "sigma0" to the system/fvSchemes dictionary.

Next, to find out what "mechanicalLaw" options are available, I edited the constant/mechanicalProperties dictionary:

planeStress     no;

mechanical
(
    steel
    {
        type            thermoMechanicalLaw;
        rho             rho [1 -3 0 0 0 0 0] 7854;
        alpha           alpha [0 0 0 -1 0 0 0] 1e-05;
        T0              T0 [0 0 0 1 0 0 0] 0;
        mechanicalLaw
        {
//            type            linearElastic;
            type            bananas;
            E               E [1 -1 -2 0 0 0 0] 2e+11;
            nu              nu [0 0 0 0 0 0 0] 0.3;
        }
    }
);

Rerunning solids4Foam resulted in an error for type "bananas" and listed the following options:

--> FOAM FATAL IO ERROR: 
Unknown mechanicalLaw type bananas

Valid linearGeometry mechanicalLaws are : 

9
(
linearElasticMohrCoulombPlastic
linearElastic
linearElasticMisesPlastic
thermoMechanicalLaw
anisotropicBiotElastic
linearElasticFromFile
viscousHookeanElastic
orthotropicLinearElastic
poroMechanicalLaw
)

Then after substituting "linearElasticMisesPlastic" in place of "bananas" in the "mechanicalProperties" dictionary and copying the format of entries from the elastoplasticity/perforatedPlate tutorial's mechanicalProperties dictionary:

planeStress     no;

mechanical
(
    steel
    {
        type            thermoMechanicalLaw;
        rho             rho [1 -3 0 0 0 0 0] 7854;
        alpha           alpha [0 0 0 -1 0 0 0] 1e-05;
        T0              T0 [0 0 0 1 0 0 0] 0;
/*
        mechanicalLaw
        {
            type            linearElastic;
            E               E [1 -1 -2 0 0 0 0] 2e+11;
            nu              nu [0 0 0 0 0 0 0] 0.3;
        }
*/
        mechanicalLaw
        {
            type            linearElasticMisesPlastic;
            rho             rho [1 -3 0 0 0 0 0] 7854;
            E               E [1 -1 -2 0 0 0 0] 2e+11;
            nu              nu [0 0 0 0 0 0 0] 0.3;
            "file|fileName" "$FOAM_CASE/constant/plasticStrainVsYieldStress";
            outOfBounds     clamp;
            solvePressureEqn no;
        }
    }
);

Rerunning solids4Foam then results in a "function not implemented for this mechanical law" error.

Creating the mechanicalModel
Selecting mechanical law thermoMechanicalLaw
Selecting mechanical law linearElasticMisesPlastic
    Plasticity is nonlinear
Creating analyticalHotCylinder function object
Time = 1000

Evolving thermal solid solver
Solving coupled energy and displacements equation for T and D
setCellDisplacements: reading cellDisplacements


--> FOAM FATAL ERROR: 
Not implemented

    From function Not yet implemented for this mechanical law
    in file lnInclude/mechanicalLaw.H at line 506.

FOAM aborting

Is it possible to run the hotCylinder tutorial with a mechanical law for a material other than "linearElastic?" If so, am I proceeding correctly?