Hello all,

I am looking to implement a plane-stress isoparametric user-element routine that I intend to use with non-linear material behaviour. For linear elastic cases, the out-of-plane strain (epsilon_33) can be obtained in closed form using Hooke's law.

epsilon_ij = ( ((1+nu)*sigma_ij) / E ) - (nu*tr(sigma)*delta_ij / E), where tr(.) is the trace and delta denotes the kronecker operator.

In the above, as sigma_13 and sigma_23 are zero, we can clearly show that epsilon_13 and epsilon_23 are zero and can therefore be neglected.

My question is, can we also adopt the same idea for large strain case with non-linear but isotropic material behaviour?

i.e., can we also approximate F_13, F_23, F_31 and F_32 (components of def_grad) as unity assuming stress measure to be 1st Piola-Kirchoff for example in a hyperelastic material model?

or

the fact that epsilon_13 and epsilon_23 are zero for Hooke's law is just a result of linear elastic assumption?

I hope my question is clear and of course I can provide more information if needed.

Thank you for your time

Shree

Similar questions and discussions