CFD Online Discussion Forums - Solids4Foam anisotropic solver

Hi Angu,

solids4foam currently contains some examples of anisotropic mechanical law, e.g. see the following tutorial:

Code:

soldis4foam-release/tutorials/solids/linearElasticity/plateHole/plateHoleOrthotropic/constant/mechanicalProperties

where constant/mechanicalProperties is given as:

Code:

mechanical

(

    carbonFibre

    {

        type            orthotropicLinearElastic;

        rho             rho [1 -3 0 0 0 0 0] 7854;

        E1              E1 [1 -1 -2 0 0 0 0] 200e+9;

        E2              E2 [1 -1 -2 0 0 0 0] 50e+9;

        E3              E3 [1 -1 -2 0 0 0 0] 50e+9;

        nu12            nu12 [0 0 0 0 0 0 0] 0.3;

        nu23            nu23 [0 0 0 0 0 0 0] 0.3;

        nu31            nu31 [0 0 0 0 0 0 0] 0.075;

        G12             G12 [1 -1 -2 0 0 0 0] 20e+9;

        G23             G23 [1 -1 -2 0 0 0 0] 20e+9;

        G31             G31 [1 -1 -2 0 0 0 0] 20e+9;

        // Local material directions can be specified here or as fields                                                                                                                                      

        // If not specified, then the Cartesian unit vectors are assumed                                                                                                                                     

        //materialDirectionX (1 0 0);                                                                                                                                                                        

        //materialDirectionY (0 1 0);                                                                                                                                                                        

        //materialDirectionZ (0 0 1);                                                                                                                                                                        

    }

);

You can see that 3 Young's module are given, 3 shear modulii and 3 Poisson's ratios.

However, in your case, you want to reproduce the results from Tang's minigeotechfoam repo. I had not ported the anisotropicBiotFoam mechanical law (it is embedded in the solver in the calculateSigmaEpsilonOrtho.H header file), but, it is straight-forward to port so I have just pushed it to the development branch (commit e297b1e2). See the new tutorial:

Code:

solids4foam-release/tutorials/solids/poroelasticity/rodAndSeabed

The stress calculation is now consistent with anisotropicBiotFoam, however, I have yet to port the pressure equation formulation with anisotropic permeability. And the other required change is to port the mapped traction boundary condition.

Have a look at the pushed code and case and let me know if it works for you (git checkout development && git pull && ./Allwmake). If you are happy to, you can modify the pEqn in:

Code:

solids4foam/src/solids4FoamModels/solidModels/poroLinGeomSolid/poroLinGeomSolid.C

Alternatively, I will do it when I get a chance (let me know if above works OK first).

Philip