Henrik Sönnerlind
COMSOL Employee
Please login with a confirmed email address before reporting spam
Posted:
3 years ago
Oct 14, 2021, 3:42 p.m. EDT
The basis of the geometrically nonlinear formulation is what often is called a Total Lagrangian formulation. The strain measure that is used is the Green-Lagrange strain.
The inner part of the expression you have posted, is the tensor multiplication of the linear constitutive tensor and the elastic part of the Green-Lagrange strain. In the absence of inelastic strains, that would be the final 2nd Piola-Kirchhoff stress.
With inelastic deformations, things get trickier. By multiplicatively removing the inelastic part of the deformation, the elastic problem can be seen as computed in an intermediate configuration. To get from the undeformed to the final state, we do two things: First the inelastic stretching, followed by the elastic stretching. So the stresses computed by the expression above are of 2nd Piola-Kirchhoff type, but not in the correct configuration. Hence the extra pull-back operation.
If you are going to use large inelastic strains, caused by intercalation, I would suggest that you use geometric nonlinearity, and then add an External Strain node. Set Strain input to Deformation gradient, inverse. Then enter a diagonal tensor, where each element is (current_volume/original_volume)^(-1/3)
-------------------
Henrik Sönnerlind
COMSOL
The basis of the geometrically nonlinear formulation is what often is called a *Total Lagrangian* formulation. The strain measure that is used is the Green-Lagrange strain.
The inner part of the expression you have posted, C:\varepsilon_{el} is the tensor multiplication of the linear constitutive tensor and the elastic part of the Green-Lagrange strain. In the absence of inelastic strains, that would be the final 2nd Piola-Kirchhoff stress.
With inelastic deformations, things get trickier. By multiplicatively removing the inelastic part of the deformation, the elastic problem can be seen as computed in an intermediate configuration. To get from the undeformed to the final state, we do two things: First the inelastic stretching, followed by the elastic stretching. So the stresses computed by the expression above are of 2nd Piola-Kirchhoff type, but not in the correct configuration. Hence the extra pull-back operation.
If you are going to use large inelastic strains, caused by intercalation, I would suggest that you use geometric nonlinearity, and then add an *External Strain* node. Set *Strain input* to *Deformation gradient, inverse*. Then enter a diagonal tensor, where each element is (current_volume/original_volume)^(-1/3)
Please login with a confirmed email address before reporting spam
Posted:
3 years ago
Oct 19, 2021, 10:40 a.m. EDT
Thank you Henrik! Now everything makes sense to me.
Thank you Henrik! Now everything makes sense to me.