Convergence problem using contact
- Comment
-
Topic Author
- Offline
- Junior Member
-
Less
More
- Posts: 25
- Thank you received: 0
11 years 1 month ago #7460
by Comment
Convergence problem using contact was created by Comment
Hello everyone,
I need help with my model of a compressive test. I currently work for a company that is interested to use Code_Aster instead of Abaqus. In order to do this, the first step is to compare the results of simple models. I have done a compressive test model with plasticity and then i have imposed a contact between the pieces, the model don't converge and after reading all the doc about this subject, i still can't identify the problem.

Is someone can help me and take a look ?
Thanks
I need help with my model of a compressive test. I currently work for a company that is interested to use Code_Aster instead of Abaqus. In order to do this, the first step is to compare the results of simple models. I have done a compressive test model with plasticity and then i have imposed a contact between the pieces, the model don't converge and after reading all the doc about this subject, i still can't identify the problem.


Is someone can help me and take a look ?
Thanks
- kwou
-
- Offline
- Moderator
-
11 years 1 month ago - 11 years 1 month ago #7464
by kwou
Interest: structural mechanics, solar energy (picture at 'my location' shows too little pv panels)
--
kind regards - kees
Replied by kwou on topic Re: Convergence problem using contact
hoi Dcomment88,
your problem is the MODI_MAILAGE command:
I didnot look into detail into the problem, but I guess you have to use ORIE_PEAU_3D=... because your problem, as stated, is 3D (mess file line 473 and further that reoccurs at the end of the file, line 635).
your problem is the MODI_MAILAGE command:
mesh2 = MODI_MAILLAGE(reuse=mesh2,
MAILLAGE=mesh2,
ORIE_PEAU_2D=_F(GROUP_MA=('ground', 'masterU', 'masterD', 'slaveU', 'slaveD', 'load', 'sym_y'),),
INFO=1,)
vous avez utilise le mot clé ORIE_PEAU_2d alors que le problème est 3d. utilisez ORIE_PEAU_3d
I didnot look into detail into the problem, but I guess you have to use ORIE_PEAU_3D=... because your problem, as stated, is 3D (mess file line 473 and further that reoccurs at the end of the file, line 635).
Interest: structural mechanics, solar energy (picture at 'my location' shows too little pv panels)
--
kind regards - kees
Last edit: 11 years 1 month ago by kwou.
- Claus
-
- Offline
- Moderator
-
Less
More
- Posts: 670
- Thank you received: 34
11 years 1 month ago #7465
by Claus
Code_Aster release : STA11.4 on OpenSUSE 12.3 64 bits - EDF/Intel version
Replied by Claus on topic Re: Convergence problem using contact
Uhm, assuming you've attached the right files, the error is right there in plain sight: vous avez utilise le mot clé ORIE_PEAU_2d alors que le problème est 3d. utilisez ORIE_PEAU_3d
Use ORIE_PEAU_3d indstead of 2d
/C
Use ORIE_PEAU_3d indstead of 2d

/C
Code_Aster release : STA11.4 on OpenSUSE 12.3 64 bits - EDF/Intel version
- Comment
-
Topic Author
- Offline
- Junior Member
-
Less
More
- Posts: 25
- Thank you received: 0
11 years 1 month ago - 11 years 1 month ago #7466
by Comment
Replied by Comment on topic Re: Convergence problem using contact
Hello,
Thanks for your quick answers
So yes initially i used ORIE_PEAU_3d but the model was not running well, so i changed it just to see what happened and i forgot to put back the right expression.
My problem is more about the settings of DEFI_CONTACT (CONTINUE VS. DISCRET and NEWTON VS. POINT_FIXE for the algorithms), how you choose them?
For this simulation i used a continuous method with POINT_FIXE for the GEOM and CONT algorithms but is there a better way to have a good convergence?
Thanks for your quick answers

So yes initially i used ORIE_PEAU_3d but the model was not running well, so i changed it just to see what happened and i forgot to put back the right expression.
My problem is more about the settings of DEFI_CONTACT (CONTINUE VS. DISCRET and NEWTON VS. POINT_FIXE for the algorithms), how you choose them?
For this simulation i used a continuous method with POINT_FIXE for the GEOM and CONT algorithms but is there a better way to have a good convergence?
Last edit: 11 years 1 month ago by Comment.
- Joël Cugnoni
-
- Offline
- Moderator
-
11 years 1 month ago - 11 years 1 month ago #7467
by Joël Cugnoni
Joël Cugnoni - a.k.a admin
www.caelinux.com
Replied by Joël Cugnoni on topic Re: Convergence problem using contact
Hi Dany,
To help you we need more details: please run a simulation with orie_peau_3d and post the message files for discrete and continue methods.
Non convergence can be caused by many different sources (integration of material law, friction stick slip, contact state oscillation, instabilities, too many constraints, mesh distortions leading to singular tangent matrix.)
Without message file it is impossible to tell what could go wrong.
As an advice, I found the following combination to be usually most robust:
discrete contact with Algo_CONT=LAGRANGIEN , REAC_GEOM=CONTROLE, ALGO_FROT=PENALISATION with E_T=100000
and for non Lin solver:
PREDICTION=ELASTIQUE, MATRICE=TANGENTE, PAS_ELAS_MINI =1E-3, Solveur MUMPS, MATRICE_DISTRIBUEE=oui, PCENT_PIVOT=70 (set mpi_nb_cpu=4 in asyl)
Another route to follow with continue method is to use penalty (start with penalty stiffnesses corresponding to Youngs modulus) for both friction and contact and point fixe for all algorithms.
Good luck and keep us posted.
Joel
To help you we need more details: please run a simulation with orie_peau_3d and post the message files for discrete and continue methods.
Non convergence can be caused by many different sources (integration of material law, friction stick slip, contact state oscillation, instabilities, too many constraints, mesh distortions leading to singular tangent matrix.)
Without message file it is impossible to tell what could go wrong.
As an advice, I found the following combination to be usually most robust:
discrete contact with Algo_CONT=LAGRANGIEN , REAC_GEOM=CONTROLE, ALGO_FROT=PENALISATION with E_T=100000
and for non Lin solver:
PREDICTION=ELASTIQUE, MATRICE=TANGENTE, PAS_ELAS_MINI =1E-3, Solveur MUMPS, MATRICE_DISTRIBUEE=oui, PCENT_PIVOT=70 (set mpi_nb_cpu=4 in asyl)
Another route to follow with continue method is to use penalty (start with penalty stiffnesses corresponding to Youngs modulus) for both friction and contact and point fixe for all algorithms.
Good luck and keep us posted.
Joel
Joël Cugnoni - a.k.a admin
www.caelinux.com
Last edit: 11 years 1 month ago by Joël Cugnoni.
- Comment
-
Topic Author
- Offline
- Junior Member
-
Less
More
- Posts: 25
- Thank you received: 0
11 years 1 month ago #7468
by Comment
Replied by Comment on topic Re: Convergence problem using contact
Hi Joel,
First thanks for your advices, my model runs perfectly with CONTINUE using a penalty on the contact. After that, i have introduced a friction (also with a penalty) but the simulation takes so much time, i started at 2 pm and at 5, it was still not finish....is it normal that the calculation takes so long?
If tomorrow my model has converged, i will try the DISCRETE approach to compare the robustness and the performance of the both method.
First thanks for your advices, my model runs perfectly with CONTINUE using a penalty on the contact. After that, i have introduced a friction (also with a penalty) but the simulation takes so much time, i started at 2 pm and at 5, it was still not finish....is it normal that the calculation takes so long?
If tomorrow my model has converged, i will try the DISCRETE approach to compare the robustness and the performance of the both method.
Moderators: catux
Time to create page: 0.145 seconds