CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)
-   -   Solids4Foam anisotropic solver (https://www.cfd-online.com/Forums/openfoam/238063-solids4foam-anisotropic-solver.html)

angupraveen August 20, 2021 08:49

Solids4Foam anisotropic solver
 
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 August 26, 2021 12:23

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


All times are GMT -4. The time now is 19:51.