[angupraveen]

Dear all,

Currently, I am trying to check the solution from the solids4Foam solver with the anisotropicBiotFoam written by another researcher Tian Tang. I would like to know, how I can reproduce these results using solids4Foam?

As I am relatively very new to solids4Foam, I would like to know what are the input files that have to be changed in order to replicate a similar result using s4F?

1. In anisotropic, we have different young’s moduli. How to add this in s4F?

2. How to know which solver to use?

The link to Tian Tang's Bitbucket is: https://bitbucket.org/tiantang/minig...utorials/soil/

Any information from you would be really helpful

Regards

Angu

[bigphil]

Hi Angu,

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

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

where constant/mechanicalProperties is given as:

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:

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:

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