Note: This discussion is about an older version of the COMSOL Multiphysics® software. The information provided may be out of date.
Discussion Locked This discussion was locked by a forum moderator.
How to deal with the nonconvergence at interface in solid mechanics
Posted 2015年5月27日 GMT+8 06:05 Mesh, Studies & Solvers, Structural Mechanics Version 4.4 8 Replies
Please login with a confirmed email address before reporting spam
Recently I encountered a problem. When I used solid mechanics module to simulate the thermal expansion deformation of a particle-matrix system(a sphere particle included in a matrix, axi-symmetric model). I cannot converge at some time. As fig. t430 and t440 shown, when time is 430s, the result is right while when time is 440s, the calculation cannot converge. However, if I only simulate the expansion of this particle, it can work well. So I want to make sure if it is the interface causes the problem. Is there any ways to deal with it? Thanks!
Best,
Rong
Please login with a confirmed email address before reporting spam
I also found when I even added a homogeneous temperature field on the center part, it cannot converge at some time. I think it is a very simple problem with large deformation, and in ABAQUS it can work well. But I don't know why I cannot handle with it in COMSOL. Attached is the pressure at different steps. Thanks!
Best,
Rong
Attachments:
Please login with a confirmed email address before reporting spam
Attachments:
Please login with a confirmed email address before reporting spam
I uploaded this COMSOL model. Can anyone help me check this model? It is a very simple model. I just added a thermal stain and want to get the stress distribution. Thanks!
Best,
Rong
Attachments:
Please login with a confirmed email address before reporting spam
For one, your model is not set up as axisymmetric.
Another problem is that the coefficient of thermal expansion for material 1 is set to a wild value (1[1/K]! ).
I hope this helps.
Best regards,
Jeff
Please login with a confirmed email address before reporting spam
Thanks so much! This model is just a simplified model. Even I build an axis-symmetric model (attached is that axis-symmetric model). It cannot converge yet. So I think it is not because of the axis-symmetric frame. Is it because the difference of Young's Module between two materials is two large?
According to the thermal expansion coefficient, finally I need a 50% thermal strain, so for simplification I set CTE equals to 1 and the temperature equals 50 0.5. Of cause I can define CTE as a smaller value while the temperature gradient would increase. So I think the result is the same. Thanks!
Best,
Rong
Attachments:
Please login with a confirmed email address before reporting spam
If you will have 50% strain, linear elasticity is not a suitable model. It will break down at large compressive strains.
There are at least two issues here:
1. Such a large thermal strain requires a multiplicative decomposition. This is used in the hyperelastic materials.
2. Interpretation of linear elasticity is different for different choices of stress and strain measures. See also www.comsol.com/blogs/why-all-these-stresses-and-strains . Two different softwares having different formulations will give vastly different results at 50% strain and linear elasticity. This is what you experienced in your first post.
Regards,
Henrik
Please login with a confirmed email address before reporting spam
Hi Henrik,
Thanks for your reply. Indeed, I have open the large deformation option including large plastic strain, so I think in my model the geometric nonlinear problem has been considered.
You said a large thermal strain requires a multiplicative decomposition. I think what you mean is the elastic deformation induced by thermal strain requires a multiplicative decomposition. However, is it right the multiplicative decomposition has been done when I used large plastic strain model? Now the elastic strain is the green-Lagrangian strain. But the model still cannot work. Do you have some advice for dealing with such a large thermal strain? Thank!
Best,
Rong
Please login with a confirmed email address before reporting spam
Hi Henrik,
I suddenly realized that only the center part of my model have considered the plasticity as well as multiplicative decomposition you mentioned while outer shell is pure elasticity. So maybe it is the point in your reply! I will modify my model soon to see if it is right. Thanks again.
Best,
Rong
Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.
Suggested Content
- KNOWLEDGE BASE Manually Setting the Scaling of Variables
- FORUM How to calculate the volume of deformation in Solid Mechanics interface
- KNOWLEDGE BASE What to do when a linear stationary model is not solving
- KNOWLEDGE BASE Structural Contact Modeling Guidelines
- BLOG How to Evaluate Stresses in COMSOL Multiphysics®
