Interference-fit analysis
- theCapitain
- Topic Author
- Offline
- Junior Member
-
Less
More
- Posts: 21
- Thank you received: 0
14 years 9 months ago #4653
by theCapitain
Replied by theCapitain on topic Re:Interference-fit analysis
You are right, the case comply with small displacements theory, I started from an exiting command file and I didn't realize that the displacement setting was on GREEN, anyway I had no convergence problem at that time.
Yesterday trying to reproduce the results using STA9.4 I noticed the inaccuracy and changed the parameter to PETIT without improvement of any kind.
Yesterday trying to reproduce the results using STA9.4 I noticed the inaccuracy and changed the parameter to PETIT without improvement of any kind.
- Joël Cugnoni
-
- Offline
- Moderator
-
14 years 9 months ago #4655
by Joël Cugnoni
Joël Cugnoni - a.k.a admin
www.caelinux.com
Replied by Joël Cugnoni on topic Re:Interference-fit analysis
After trying to run your case, I see that actually the first increment fails to converges because the relative residual (=the force "mis" equilibrium divided by the total applied force) is too high; but the absolute value of the residue is very low which means that the system is in equilibirum. To converge, by default, Aster requires that both absolute and relative residuals are below a certain value. This convergence problem is typical of cases where you don't have external load applied to the system and thus, relative residual is always something like 1/0=inf. and thuis does not converge.
So in your case, you should ask Aster to consider only the absolute residual (in unit force, so Newton) in the convergence check. So try to use the following settings which will enable ONLY global absolute residual magnitude verification.
CONVERGENCE=_F(RESI_GLOB_MAXI=1e-4,
ITER_GLOB_MAXI=60,),
Then your problem should converge.
In the following file, I also have increased the number of iterations for contact and the number of increments, but it should not be necessary... This proble has been running successfully on the Aster 10.1 OpenMP included in Salome-Meca 2010.
[file]<br /><br />Post edited by: Administrator, at: 2010/09/01 19:14
So in your case, you should ask Aster to consider only the absolute residual (in unit force, so Newton) in the convergence check. So try to use the following settings which will enable ONLY global absolute residual magnitude verification.
CONVERGENCE=_F(RESI_GLOB_MAXI=1e-4,
ITER_GLOB_MAXI=60,),
Then your problem should converge.
In the following file, I also have increased the number of iterations for contact and the number of increments, but it should not be necessary... This proble has been running successfully on the Aster 10.1 OpenMP included in Salome-Meca 2010.
[file]<br /><br />Post edited by: Administrator, at: 2010/09/01 19:14
Joël Cugnoni - a.k.a admin
www.caelinux.com
- Joël Cugnoni
-
- Offline
- Moderator
-
14 years 9 months ago #4656
by Joël Cugnoni
Joël Cugnoni - a.k.a admin
www.caelinux.com
Replied by Joël Cugnoni on topic Re:Interference-fit analysis
Attachment interferencefit-a419439f57afb3f2e17356202fe710e3.gz not found
Joël Cugnoni - a.k.a admin
www.caelinux.com
- JMB
- Offline
- Elite Member
-
Less
More
- Posts: 166
- Thank you received: 0
14 years 9 months ago #4658
by JMB
Replied by JMB on topic Re:Interference-fit analysis
Hello all,
I still believe there is a regression, since it will not run in Ver 10.1.27, 10.2.0 or 10.2.09. I had to make a minor change to move the CONTACT option outside AFFE_CHAR_MECA to DEFI_CONTACT, since it is no longer permitted under AFFE_CHAR_MECA. Other than that I have made no other changes. My *.mess file(s) are attached for review.
What does work on V10.2 is:
[code:1]
Fit=DEFI_CONTACT(MODELE=Mod_Lin,
FORMULATION='CONTINUE',
REAC_GEOM='CONTROLE',
NB_ITER_GEOM=2,
ITER_CONT_MAXI=50,
ZONE=_F(APPARIEMENT='MAIT_ESCL',
GROUP_MA_MAIT='cEX',
GROUP_MA_ESCL='cIN',
DIST_MAIT=Rad_int,
CONTACT_INIT='OUI',),);
[/code:1]
i.e. adding the REAC_GEOM='CONTROLE' and NB_ITER_GEOM=3. I have not examined the results carefully under Salome 5.1.4 yet.
Regards,
JMB
I still believe there is a regression, since it will not run in Ver 10.1.27, 10.2.0 or 10.2.09. I had to make a minor change to move the CONTACT option outside AFFE_CHAR_MECA to DEFI_CONTACT, since it is no longer permitted under AFFE_CHAR_MECA. Other than that I have made no other changes. My *.mess file(s) are attached for review.
What does work on V10.2 is:
[code:1]
Fit=DEFI_CONTACT(MODELE=Mod_Lin,
FORMULATION='CONTINUE',
REAC_GEOM='CONTROLE',
NB_ITER_GEOM=2,
ITER_CONT_MAXI=50,
ZONE=_F(APPARIEMENT='MAIT_ESCL',
GROUP_MA_MAIT='cEX',
GROUP_MA_ESCL='cIN',
DIST_MAIT=Rad_int,
CONTACT_INIT='OUI',),);
[/code:1]
i.e. adding the REAC_GEOM='CONTROLE' and NB_ITER_GEOM=3. I have not examined the results carefully under Salome 5.1.4 yet.
Regards,
JMB
Attachment Interffit.gz not found
- theCapitain
- Topic Author
- Offline
- Junior Member
-
Less
More
- Posts: 21
- Thank you received: 0
14 years 9 months ago #4664
by theCapitain
Replied by theCapitain on topic Re:Interference-fit analysis
Administrator wrote:
By the way, have you any idea of why, even with this measure, there are problems using the linear-quad conversion in Aster (this is also the subject of my first post)? In this case we get a solution from which it's clear that some slave nodes penetrate the master surface.
Thank you very much
@JMB
I'm sorry, I'm unable to try CA10.2 right now and I'm not experienced enough to judge from your command file.
Perfect, it works for me (in the meantime I downloaded CAElinux 2010). After the first load step it also converges on the relative residual. Remains a mystery the convergence on my home PC.Then your problem should converge.
By the way, have you any idea of why, even with this measure, there are problems using the linear-quad conversion in Aster (this is also the subject of my first post)? In this case we get a solution from which it's clear that some slave nodes penetrate the master surface.
Thank you very much
@JMB
I'm sorry, I'm unable to try CA10.2 right now and I'm not experienced enough to judge from your command file.
Moderators: catux
Time to create page: 0.160 seconds