Note: This discussion is about an older version of the COMSOL Multiphysics® software. The information provided may be out of date.
Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.
Large deformation of a 2-D square
Posted 17.09.2013, 18:53 GMT-4 Structural Mechanics Version 4.3b 2 Replies
Please login with a confirmed email address before reporting spam
Hello All
I am trying to apply a large shear strain on a square domain. The geometry is very simple, just a square where on the upper and lower edges, opposite velocities are applied. In order to prevent any rotation, I assume that the upper edge is fixed vertically, and I turn on the nonlinear geometry option.
I get divergence problems and to solve it I have tried every possible method (adaptive mesh, automatic refinement and even mesh displacement physics), nothing works. The process stops at 13% progress.
Thank you,
Hossein
I am trying to apply a large shear strain on a square domain. The geometry is very simple, just a square where on the upper and lower edges, opposite velocities are applied. In order to prevent any rotation, I assume that the upper edge is fixed vertically, and I turn on the nonlinear geometry option.
I get divergence problems and to solve it I have tried every possible method (adaptive mesh, automatic refinement and even mesh displacement physics), nothing works. The process stops at 13% progress.
Thank you,
Hossein
Attachments:
2 Replies Last Post 18.09.2013, 13:28 GMT-4
Please login with a confirmed email address before reporting spam
Posted:
1 decade ago
Hello Hossein,
The cause of your problem has a general interest, since it often turns up in support cases. To identify it, plot the smallest principal strain in the last converged time or parameter step. It is of the order -0.4.
When using geometric nonlinearity with a Linear Elastic material, the linearity actually says that the 2:nd Piola-Kirchhoff stress is proportional to the Green-Lagrange strain. (This is sometimes called a St. Venant-Kichhoff material). If you would model uniaxial tension with the same data, the force vs. displacement relation is far from linear at larger strains. It can actually be shown that the force is not unique at high compressive strains, and the limit is approximately -0.4.
Conclusion: If there are high compressive strains in your model, a Linear Elastic material will break down when the smallest principal strain is of the order of -0.4. It is necessary to use another material model.
Regards,
Henrik
The cause of your problem has a general interest, since it often turns up in support cases. To identify it, plot the smallest principal strain in the last converged time or parameter step. It is of the order -0.4.
When using geometric nonlinearity with a Linear Elastic material, the linearity actually says that the 2:nd Piola-Kirchhoff stress is proportional to the Green-Lagrange strain. (This is sometimes called a St. Venant-Kichhoff material). If you would model uniaxial tension with the same data, the force vs. displacement relation is far from linear at larger strains. It can actually be shown that the force is not unique at high compressive strains, and the limit is approximately -0.4.
Conclusion: If there are high compressive strains in your model, a Linear Elastic material will break down when the smallest principal strain is of the order of -0.4. It is necessary to use another material model.
Regards,
Henrik
Please login with a confirmed email address before reporting spam
Posted:
1 decade ago
Dear Henrik
Thank you for your kind reply. I had thought of every possible issue except the material. The point is that now by hyper elastic material (Neo-Hookean, with Lambda=0 , Mu=10^6) still I see the divergence when the principal strain is around -0.45. Do you have any suggestion on a more proper material choice?
Thank you,
Hossein
Thank you for your kind reply. I had thought of every possible issue except the material. The point is that now by hyper elastic material (Neo-Hookean, with Lambda=0 , Mu=10^6) still I see the divergence when the principal strain is around -0.45. Do you have any suggestion on a more proper material choice?
Thank you,
Hossein