×

Notice

The forum is in read only mode.

Convergence problem using contact

  • Comment
  • Comment's Avatar Topic Author
  • Offline
  • Junior Member
  • Junior Member
More
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
Attachments:
More
11 years 1 month ago - 11 years 1 month ago #7464 by kwou
Replied by kwou on topic Re: Convergence problem using contact
hoi Dcomment88,

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.
More
11 years 1 month ago #7465 by Claus
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

Code_Aster release : STA11.4 on OpenSUSE 12.3 64 bits - EDF/Intel version
  • Comment
  • Comment's Avatar Topic Author
  • Offline
  • Junior Member
  • Junior Member
More
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?
Last edit: 11 years 1 month ago by Comment.
More
11 years 1 month ago - 11 years 1 month ago #7467 by Joël Cugnoni
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

Joël Cugnoni - a.k.a admin
www.caelinux.com
Last edit: 11 years 1 month ago by Joël Cugnoni.
  • Comment
  • Comment's Avatar Topic Author
  • Offline
  • Junior Member
  • Junior Member
More
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.
Moderators: catux
Time to create page: 0.145 seconds
Powered by Kunena Forum