Dear Samyak, It is definitely not true that the hydrostatic term -1/3*Tr(sigma) should be zero in an incompressible material. If this is the case you are simply calculating the deviatoric part of sigma instead of sigma. In order to get sigma you need to add the term p*Id(3). Actually the hydrostatic term you are looking for is simply equal to p.
Best regards Kostas On Thu, Feb 9, 2017 at 4:20 AM, samyak jain <[email protected]> wrote: > Dear getfem-users, > > I am currently trying to solve a contact problem between a hyperelastic > rubber and a rigid bosy and I need to calculate the pressure values on > either on the rubber. > > I am using Incompressible Mooney-Rivlin Hyperelastic law and if my model I > am adding also adding finite strain incompressibility brick. > > Now when I calculate Cauchy Stress from second piola kirchhoff stress, I > am getting the Hydrostatic term (-1/3Tr(sigma)) of the cauchy stress tensor > as zero which is what it should be as the material is incompressible. > > So, is there a way is getfem to calculate the hydrostatic pressure term > for such incompressible materials.I believe treating the material as > nearly incompressible (Poisson's ratio 0.499) is one way to solve it but I > don't know how it works or if it is implemented in the model. > > Could you guys please provide any help or suggestion to calculate the > hydrostatic pressure for such a case. > > Thanks a lot. > > Yours sincerely > Samyak > > _______________________________________________ > Getfem-users mailing list > [email protected] > https://mail.gna.org/listinfo/getfem-users > >
_______________________________________________ Getfem-users mailing list [email protected] https://mail.gna.org/listinfo/getfem-users
