Posts by @@@duffy@@@

@@@duffy@@@ · Code_Aster usage

Hi, all.
I have changed the title of the post.
I attach a mess file in which contact iteration does not converge after 10 iterations.
I would like to get convergence for this plastic contact problem and to proceed to more deformation.
Any suggestion would be appreciated.
Thanks in advance.
Regards,
@@@duffy@@@

@@@duffy@@@ · Code_Aster usage

Hi, all,
I'm performing plastic contact analysis.
But the deformation does not proceed as I expected.
Analysis won't converge at some load level.

Here's a convergence report of my analysis at some load level. At this step, analysis has converged.
How do I interpret this report, and what can I do to prevent from non-convergence to happen.
I am not sure which residual matters in my case.
Any suggestion would be greatly appreciated.
Please find attached a *.comm file for your reference.
Regards,
@@@duffy@@@

Instant deWrite failed: Broken pipe-02 - Niveau de découpe: 1
-----------------------------------------------------------------------------------------------------------------------------------------
miyoshi@murakumo ~/workspace/TetMesh_P RESIDU | RESIDU | OPTION | CONTACT | CONTACT | CONTACT |
$ BCL. GEOM. | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | NEWTON GENE | INFOS | CRITERE |
| ITERATION | | RESI_GLOB_RELA | RESI_GLOB_MAXI | | VARI. CONT. | CYCLAGES | VALEUR |
-----------------------------------------------------------------------------------------------------------------------------------------
| 1 X | 0 X | 1.51072E-01 X | 9.85599E-05 |TANGENTE | 0 | | - SANS OBJET - |
| 1 X | 1 X | 3.19511E-04 X | 2.08246E-07 |TANGENTE | 0 | | - SANS OBJET - |
| 1 X | 2 X | 6.78194E-06 X | 4.42018E-09 |TANGENTE | 0 | | - SANS OBJET - |
| 1 X | 3 | 3.17702E-09 | 2.07065E-12 |TANGENTE | 0 | | 1.00000E-01 |
-----------------------------------------------------------------------------------------------------------------------------------------
| 2 X | 0 X | 1.51075E-01 X | 9.85622E-05 |TANGENTE | 0 | | - SANS OBJET - |
| 2 X | 1 X | 3.18884E-04 X | 2.07839E-07 |TANGENTE | 0 | | - SANS OBJET - |
| 2 X | 2 X | 6.77883E-06 X | 4.41818E-09 |TANGENTE | 0 | | - SANS OBJET - |
| 2 | 3 | 3.18725E-09 | 2.07732E-12 |TANGENTE | 0 | | 4.86842E-04 |
-----------------------------------------------------------------------------------------------------------------------------------------

@@@duffy@@@ · Code_Aster usage

Hi, AsterO'dactyle,
Thank you for your suggestion.

>1/ Why didn't you choose your field instead of "Solid Color" in field selection in paravis ?
I tried this, and I found colors of Gaussian integration points reflect the values of respective variable. The equivalent plastic strain was not included in the components of EPSI_ELGA or EPSP_ELGA.

>2/ The plastic equivalent strain is an internal variable of behavior law:
I tried this both the first option and the second option, and I successfully have shown just what I wanted to see.
Regards,
@@@duffy@@@

@@@duffy@@@ · Code_Aster usage

Hi all,
I would like to view plastic equivalent strain color contour plot either on the mesh of or on the integration points for my plastic contact analysis.
It seems contour plot is available only for variables evaluated on nodes, not for variable evaluated on Gaussian integration points.

I tried EPSP_NOEU, but there is no plastic equivalent strain in the component list. Magnitude is available. But that is not what I want to see.
I wrote EPSP_ELNO to an rmed file. But I cannot find how to view it in Paravis.
I tried EPSP_ELGA following the instructions below:

>to visualize ELGA fields:
>1 select in proprieties the result with this symbol Σ-
>2 in filters -> Integration points -> selection Gauss Points
>3 in the gauss points filter select your elga field in display "color by"
in
http://code-aster.org/forum2/viewtopic.php?id=19773

But I could not find display "color by".
What I got is the attached screen capture. Gaussian integration points are shown, but not color contoured.
I paste the *.comm file I used for you.
I am using Salome 7(2015.1).
Any suggestion would be greatly appreciated.
Regards,
@@@duffy@@@

DEBUT(LANG='EN',);

elast_pl=DEFI_FONCTION(
NOM_PARA='EPSI',
NOM_RESU='SIGM',
VALE=(0.001072157 ,73.6 ,
0.005 ,75.1 ,),
PROL_DROITE='LINEAIRE',);

elast_ma=DEFI_FONCTION(
NOM_PARA='EPSI',
NOM_RESU='SIGM',
VALE=(0.001066839 ,205.9 ,
0.005 ,217.3 ,),
PROL_DROITE='LINEAIRE',);

MA=DEFI_MATERIAU(ELAS=_F(E=193000.0,
NU=0.3,),
TRACTION=_F(SIGM=elast_ma,),);

MA2=DEFI_MATERIAU(ELAS=_F(E=68600.0,
NU=0.33,),
TRACTION=_F(SIGM=elast_pl,),);

MAIL=LIRE_MAILLAGE(FORMAT='MED',);

MAIL=MODI_MAILLAGE(reuse =MAIL,
MAILLAGE=MAIL,
ORIE_PEAU_3D=_F(GROUP_MA='face_load',),);

MODE=AFFE_MODELE(MAILLAGE=MAIL,
AFFE=_F(TOUT='OUI',
PHENOMENE='MECANIQUE',
MODELISATION='3D',),);

MATE=AFFE_MATERIAU(MAILLAGE=MAIL,
AFFE=(_F(GROUP_MA='block_low',
MATER=MA,),
_F(GROUP_MA='block_high',
MATER=MA2,),),);

CHAR=AFFE_CHAR_MECA(MODELE=MODE,
DDL_IMPO=(_F(GROUP_MA='face_fix',
DX=0.0,
DY=0.0,
DZ=0.0,),
_F(GROUP_MA='face_fixX',
DX=0.0,),
_F(GROUP_MA='face_fixY',
DY=0.0,),),);

loadP=AFFE_CHAR_MECA(MODELE=MODE,
DDL_IMPO=_F(GROUP_MA='face_load',
DZ=-0.00035,),);

contact=DEFI_CONTACT(MODELE=MODE,
FORMULATION='CONTINUE',
ITER_GEOM_MAXI=6,
ALGO_RESO_CONT='NEWTON',
ZONE=_F(GROUP_MA_MAIT='face_contact_low',
GROUP_MA_ESCL='face_contact_high',),);

ramp=DEFI_FONCTION(
NOM_PARA='INST',
VALE=(0.0 ,0.0 ,
1.0 ,1.0 ,),);

inst=DEFI_LIST_REEL(DEBUT=0.0,
INTERVALLE=(_F(JUSQU_A=0.03,
NOMBRE=1,),
_F(JUSQU_A=0.048,
NOMBRE=12,),
_F(JUSQU_A=0.05,
NOMBRE=1,),
_F(JUSQU_A=0.1,
NOMBRE=1,),
_F(JUSQU_A=0.4,
NOMBRE=3,),
_F(JUSQU_A=0.6,
NOMBRE=1,),),);

RESU=STAT_NON_LINE(MODELE=MODE,
CHAM_MATER=MATE,
EXCIT=(_F(CHARGE=CHAR,),
_F(CHARGE=loadP,
FONC_MULT=ramp,),),
CONTACT=contact,
COMP_INCR=_F(RELATION='VMIS_ISOT_TRAC',
DEFORMATION='SIMO_MIEHE',
TOUT='OUI',),
INCREMENT=_F(LIST_INST=inst,),
NEWTON=_F(REAC_INCR=1,
PREDICTION='TANGENTE',
MATRICE='TANGENTE',
REAC_ITER=1,),
CONVERGENCE=_F(ITER_GLOB_MAXI=30,),
ARCHIVAGE=_F(PAS_ARCH=1,),);

RESU=CALC_CHAMP(reuse =RESU,
MODELE=MODE,
CHAM_MATER=MATE,
RESULTAT=RESU,
CONTRAINTE='SIGM_NOEU',
DEFORMATION=('EPSI_ELGA','EPSP_ELGA','EPSI_ELNO','EPSP_NOEU','EPSP_ELNO','EPSI_NOEU',),
CRITERES='SIEQ_NOEU',
FORCE='REAC_NODA',);

IMPR_RESU(FORMAT='MED',
UNITE=80,
RESU=_F(MAILLAGE=MAIL,
RESULTAT=RESU,
NOM_CHAM=
('SIGM_NOEU','SIEQ_NOEU','DEPL','EPSI_ELGA','EPSP_ELGA','EPSI_ELNO','EPSI_NOEU',
'EPSP_ELNO','EPSP_NOEU',),),);

sumReacF=POST_RELEVE_T(ACTION=_F(OPERATION='EXTRACTION',
INTITULE='ReacForc',
RESULTAT=RESU,
NOM_CHAM='REAC_NODA',
GROUP_NO='face_load',
RESULTANTE=('DX','DY','DZ',),),);

FIN();
#CHECKSUM:b9afb6ae7322ecc2d6ccc19c36945646 -:FIN CHECKSUM

@@@duffy@@@ · Code_Aster usage

Hi, reiteres:
Thank you again for your valuable comment.
I made a mesh group named 'fix_c' on the bottom of the part 'Base' and changed the .comm file as follows.
Now the calculation results show what I wanted.

23c23
< DDL_IMPO=(_F(GROUP_MA='Base_c',
---
> DDL_IMPO=(_F(GROUP_MA='fix_c',

I attach the whole .comm file contents for reference.
Best regards,
@@@duffy@@@

DEBUT();

MA=DEFI_MATERIAU(ELAS=_F(E=193000.0,
NU=0.3,),);

MA2=DEFI_MATERIAU(ELAS=_F(E=68600.0,
NU=0.33,),);

MAIL=LIRE_MAILLAGE(FORMAT='MED',);

MAIL=MODI_MAILLAGE(reuse =MAIL,
MAILLAGE=MAIL,
ORIE_PEAU_3D=_F(GROUP_MA='press_c',),);

MODE=AFFE_MODELE(MAILLAGE=MAIL,
AFFE=_F(TOUT='OUI',
PHENOMENE='MECANIQUE',
MODELISATION='3D',),);

MATE=AFFE_MATERIAU(MAILLAGE=MAIL,
AFFE=(_F(GROUP_MA='Base_c',
MATER=MA,),
_F(GROUP_MA='Bar_v',
MATER=MA2,),),);

CHAR=AFFE_CHAR_MECA(MODELE=MODE,
DDL_IMPO=(_F(GROUP_MA='fix_c',
DX=0.0,
DY=0.0,
DZ=0.0,),
_F(GROUP_MA='press_c',
DX=0.0,
DY=0.0,),),);

loadP=AFFE_CHAR_MECA(MODELE=MODE,
DDL_IMPO=_F(GROUP_MA='press_c',
DZ=-0.0002,),);

contact=DEFI_CONTACT(MODELE=MODE,
FORMULATION='CONTINUE',
ITER_GEOM_MAXI=100,
ALGO_RESO_CONT='NEWTON',
ZONE=_F(GROUP_MA_MAIT='contBase_c',
GROUP_MA_ESCL='contBar_c',),);

ramp=DEFI_FONCTION(
NOM_PARA='INST',
VALE=(0.0 ,0.0 ,
1.0 ,1.0 ,),);

inst=DEFI_LIST_REEL(DEBUT=0.0,
INTERVALLE=(_F(JUSQU_A=0.01,
NOMBRE=5,),
_F(JUSQU_A=0.02,
NOMBRE=2,),
_F(JUSQU_A=0.06,
NOMBRE=4,),
_F(JUSQU_A=0.1,
NOMBRE=2,),
_F(JUSQU_A=0.2,
NOMBRE=2,),
_F(JUSQU_A=0.6,
NOMBRE=4,),
_F(JUSQU_A=1.0,
NOMBRE=2,),),);

RESU=STAT_NON_LINE(MODELE=MODE,
CHAM_MATER=MATE,
EXCIT=(_F(CHARGE=CHAR,),
_F(CHARGE=loadP,
FONC_MULT=ramp,),),
CONTACT=contact,
COMP_ELAS=_F(RELATION='ELAS',
DEFORMATION='PETIT',
TOUT='OUI',),
INCREMENT=_F(LIST_INST=inst,),
NEWTON=_F(REAC_INCR=1,
MATRICE='TANGENTE',
REAC_ITER=1,),
CONVERGENCE=_F(ITER_GLOB_MAXI=30,),
ARCHIVAGE=_F(PAS_ARCH=1,),);

RESU=CALC_CHAMP(reuse =RESU,
MODELE=MODE,
CHAM_MATER=MATE,
RESULTAT=RESU,
CONTRAINTE=('SIGM_ELNO','SIGM_NOEU',),
DEFORMATION=('EPSI_NOEU','EPSI_ELNO','EPSP_NOEU',),
CRITERES=('SIEQ_NOEU','SIEQ_ELNO',),
FORCE='REAC_NODA',);

IMPR_RESU(FORMAT='MED',
UNITE=80,
RESU=_F(MAILLAGE=MAIL,
RESULTAT=RESU,
NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','DEPL',),),);

FIN();
#CHECKSUM:4e8fe363d334a9fa76acd54f2467ae9a -:FIN CHECKSUM

@@@duffy@@@ · Code_Aster usage

Hi,
After I switched type of groups from volume to face for groups involved in contact, Code_Aster calculation finished successfully.
But the Base, the part compressed indirectly through the part Bar does not seem to deform nor has meaningful stresses.
The part Bar seems to deformed and has non-zero stresses.
What is wrong with my comm file?
Any suggestion would be appreciated.
Thank you in advance.
Best regards,
@@@duffy@@@

@@@duffy@@@ · Code_Aster usage

Hi, reiteres:
Thank you for your quick reply.
I will modify the group from volume to contact faces.
Best regards,
@@@Duffy@@@

@@@duffy@@@ · Code_Aster usage

Hi,
I'm novice Salome-meca/code_Aster user.
I try to solve a simple contact problem in which a part 'bar' is on top of another part 'base' and pressed downward.

The essential part of the error message is as follows:
!-------------------------------------------------------!
! <EXCEPTION> <DVP_1> !
! !
! Erreur de programmation : condition non respectée. !

I do not know how to fixed the problem definition although I searched the forum.
I will be glad if somebody will give me suggestion to handle this error.
Please find attache the *.mess file, *.comm file and *.mmes file archived in a file.

I paste the .comm file in the post.
Thank you in advance.
Best regards,
@@@duffy@@@

DEBUT();

MA=DEFI_MATERIAU(ELAS=_F(E=130300.0,
NU=0.343,),);

MAIL=LIRE_MAILLAGE(FORMAT='MED',);

MAIL=MODI_MAILLAGE(reuse =MAIL,
MAILLAGE=MAIL,
ORIE_PEAU_3D=_F(GROUP_MA='press_c',),);

MODE=AFFE_MODELE(MAILLAGE=MAIL,
AFFE=_F(TOUT='OUI',
PHENOMENE='MECANIQUE',
MODELISATION='3D',),);

MATE=AFFE_MATERIAU(MAILLAGE=MAIL,
AFFE=_F(TOUT='OUI',
MATER=MA,),);

CHAR=AFFE_CHAR_MECA(MODELE=MODE,
DDL_IMPO=(_F(GROUP_MA='Base_c',
DX=0.0,
DY=0.0,
DZ=0.0,),
_F(GROUP_MA='press_c',
DX=0.0,
DY=0.0,),),);

loadP=AFFE_CHAR_MECA(MODELE=MODE,
DDL_IMPO=_F(GROUP_MA='press_c',
DZ=-0.0002,),);

contact=DEFI_CONTACT(MODELE=MODE,
FORMULATION='CONTINUE',
ITER_GEOM_MAXI=100,
ALGO_RESO_CONT='NEWTON',
ZONE=_F(GROUP_MA_MAIT='Base_c',
GROUP_MA_ESCL='Bar_c',),);

ramp=DEFI_FONCTION(
NOM_PARA='INST',
VALE=(0.0 ,0.0 ,
1.0 ,1.0 ,),);

inst=DEFI_LIST_REEL(DEBUT=0.0,
INTERVALLE=(_F(JUSQU_A=0.01,
NOMBRE=5,),
_F(JUSQU_A=0.02,
NOMBRE=2,),
_F(JUSQU_A=0.06,
NOMBRE=4,),
_F(JUSQU_A=0.1,
NOMBRE=2,),
_F(JUSQU_A=0.2,
NOMBRE=2,),
_F(JUSQU_A=0.6,
NOMBRE=4,),
_F(JUSQU_A=1.0,
NOMBRE=2,),),);

RESU=STAT_NON_LINE(MODELE=MODE,
CHAM_MATER=MATE,
EXCIT=(_F(CHARGE=CHAR,),
_F(CHARGE=loadP,
FONC_MULT=ramp,),),
CONTACT=contact,
COMP_ELAS=_F(RELATION='ELAS',
DEFORMATION='PETIT',
TOUT='OUI',),
INCREMENT=_F(LIST_INST=inst,),
NEWTON=_F(REAC_INCR=1,
MATRICE='TANGENTE',
REAC_ITER=1,),
CONVERGENCE=_F(ITER_GLOB_MAXI=30,),
ARCHIVAGE=_F(PAS_ARCH=1,),);

RESU=CALC_CHAMP(reuse =RESU,
MODELE=MODE,
CHAM_MATER=MATE,
RESULTAT=RESU,
CONTRAINTE=('SIGM_ELNO','SIGM_NOEU',),
DEFORMATION=('EPSI_NOEU','EPSI_ELNO','EPSP_NOEU',),
CRITERES=('SIEQ_NOEU','SIEQ_ELNO',),
FORCE='REAC_NODA',);

IMPR_RESU(FORMAT='MED',
UNITE=80,
RESU=_F(MAILLAGE=MAIL,
RESULTAT=RESU,
NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','DEPL',),),);

FIN();
#CHECKSUM:4309bf144e0487bf0ab2f823b62b9a87 -:FIN CHECKSUM

#1 Re: Code_Aster usage » How to achieve convergence for plastic contact » 2016-07-29 03:04:54

#2 Code_Aster usage » How to achieve convergence for plastic contact » 2016-07-28 16:14:27

#3 Re: Code_Aster usage » (Closed)How to view color contor of plas equiv strain in paravis? » 2016-07-28 12:15:16

#4 Code_Aster usage » (Closed)How to view color contor of plas equiv strain in paravis? » 2016-07-27 06:02:08

#5 Re: Code_Aster usage » Closed:condition non respectée in elastic contact problem » 2016-04-27 01:35:38

#6 Re: Code_Aster usage » Closed:condition non respectée in elastic contact problem » 2016-04-25 08:28:51

#7 Re: Code_Aster usage » Closed:condition non respectée in elastic contact problem » 2016-04-22 11:28:00

#8 Code_Aster usage » Closed:condition non respectée in elastic contact problem » 2016-04-21 14:12:55

Board footer