Welcome to the forums. Please post in English or French.

You are not logged in.

#1 2021-10-20 11:26:23

TomekBeginner
Member
Registered: 2021-10-20
Posts: 7

Shells - DKT - Plastsicity - (In)correct Mises

Hello,

I'm very new to code-aster and i can't find solution to my problem.

I don't know how to obtain correct VMIS stresses from my model (or model is incorrect).
Results show max 493 MPa while material is set to max about 355 MPa.

Deformation and cumulated plastic deformation seems to be ok.

Please help.

Regards,
Tomek

P.S.: .comm file:

#####Units: mm, N, N/mm2=MPa

#####Parameters:
###Shell:
T=12.5		#Shell thickness,

###Load:
N=200e3		#Total force
b=75		#Width of load area
l=200		#Length of load area
A=b*l		#Area of load
q=N/A		#Surface load

###Iterations:
T1=10		#Number of iterations

###Materials:
#Steel:
Es=210e+3
Ep=2.1
fy=355

######Start: 
DEBUT(LANG='EN')


#####Mesh:
###Loading:
mesh = LIRE_MAILLAGE(FORMAT='MED',
                     UNITE=20)
###Groups:
mesh = DEFI_GROUP(reuse=mesh,
				  MAILLAGE=mesh,
				  CREA_GROUP_MA=(_F(TOUT="OUI",
									NOM='lgma'),),
				  CREA_GROUP_NO=(_F(GROUP_MA='bottom'),),
									)
									
#####Model:
model = AFFE_MODELE(AFFE=(_F(MODELISATION=('DKT', ),
							PHENOMENE='MECANIQUE',
                            TOUT="OUI"),),
							MAILLAGE=mesh)
					
######Elements:
elemprop = AFFE_CARA_ELEM(COQUE=_F(EPAIS=T,
                                   GROUP_MA='lgma',
								   COQUE_NCOU=3,
								   ),		
                          MODELE=model)
						  
#####Materials:
###Properties:						
mater = DEFI_MATERIAU(ELAS=_F(E=Es,
                              NU=0.3),
					  ECRO_LINE=_F(D_SIGM_EPSI=Ep,
								   SY=fy))
###Assign:
fieldmat = AFFE_MATERIAU(AFFE=_F(MATER=(mater, ),
                                 TOUT="OUI"),
                         MAILLAGE=mesh)

#####Loads and BC:
###Loads:
load = AFFE_CHAR_MECA(FORCE_COQUE=_F(FZ=-q,
                                     GROUP_MA=('top', )),
                      MODELE=model)
###BC's:					  
support = AFFE_CHAR_MECA(DDL_IMPO=_F(DX=0.0,
                                     DY=0.0,
									 DZ=0.0,
                                     GROUP_MA=('bottom', )),
                         MODELE=model)

#####Functions:
###Loads:
loadf = DEFI_FONCTION(NOM_PARA='INST',
					  VALE=(0.0,  0.0,
							T1,   1.0),
					  PROL_DROITE='EXCLU',
					  PROL_GAUCHE='EXCLU')

#####Lists
###Iterations:
listiter = DEFI_LIST_REEL(DEBUT=0.0,
						  INTERVALLE=_F(JUSQU_A=T1,
										NOMBRE=T1))
						 
#####Calculation:
###Nonlinear:
results = STAT_NON_LINE(CARA_ELEM=elemprop,
                        CHAM_MATER=fieldmat,
                        EXCIT=(_F(CHARGE=load,
								  FONC_MULT=loadf),
                               _F(CHARGE=support)),
						COMPORTEMENT=_F(RELATION='VMIS_ISOT_LINE',
									 DEFORMATION='PETIT',
									 TOUT='OUI'),
                        INCREMENT=_F(LIST_INST=listiter),
						CONVERGENCE=_F(ITER_GLOB_MAXI=20,
                                       ITER_GLOB_ELAS=25),
						NEWTON=_F(REAC_INCR=1,
                                  MATRICE='TANGENTE',
                                  REAC_ITER=1,
                                  REAC_ITER_ELAS=1,
                                  MATR_RIGI_SYME='NON'),
                        MODELE=model)

#####Prepare of results:
results = CALC_CHAMP(reuse=results,
                     CONTRAINTE=('SIGM_ELGA',
								 ),
				     DEFORMATION=('EPSI_ELGA',
								  ),
					 CRITERES=('SIEQ_ELGA',
							   ),
                     RESULTAT=results)
	
#####Results:
###Files units:
#MED=80
#RESU=19
#TABLE=8
###
IMPR_RESU(FORMAT='RESULTAT',
          RESU=(_F(RESULTAT=results,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='SIEQ_ELGA',
                   NOM_CMP='VMIS',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),
				_F(RESULTAT=results,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='DEPL',
                   NOM_CMP=('DX','DY','DZ'), 
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),				   	
				_F(RESULTAT=results,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='VARI_ELGA',
                   #NOM_CMP='VMIS',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI')					
				),
          UNITE=19)

#####End:
FIN()			

Last edited by TomekBeginner (2021-10-20 12:56:30)


Attachments:
rhs_plasticity.7z, Size: 16.36 KiB, Downloads: 28

Offline

#2 2021-10-22 07:19:57

tpa
Member
From: Denmark
Registered: 2009-04-14
Posts: 61

Re: Shells - DKT - Plastsicity - (In)correct Mises

I have not tried your case but I have one observation.

When postprocessing shell/plate elements you need to extract the result at a certain position across the thickness with post_champ. Something like the below code which is taken from the book (highly recommended) by Jean-Pierre Aubry, "Beginning with Code Aster", p90. You can also look in u2.02.01 section 2.4.2.

stat2=POST_CHAMP(
RESULTAT=stat ,
GROUP_MA= (’panel’ , ) ,
EXTR_COQUE=_F(
NUME_COUCHE=1 ,
NIVE_COUCHE=’SUP’ ,
NOM_CHAM= (’SIGM_ELNO’ , ) ,
) ,);

Last edited by tpa (2021-10-22 07:56:11)

Offline

#3 2021-10-22 11:46:38

TomekBeginner
Member
Registered: 2021-10-20
Posts: 7

Re: Shells - DKT - Plastsicity - (In)correct Mises

Thank you for advise tpa.

Well...
- it worked for DKT,

 --------------------------------------------------------------------------------
 ASTER 14.04.00 CONCEPT results CALCULE LE 22/10/2021 A 12:30:09 DE TYPE         
 EVOL_NOLI                                                                       


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  SIEQ_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
LA VALEUR MAXIMALE DE VMIS    EST:  0.49367E+03
OBTENUE DANS LA MAILLE M2759   AU SOUS_POINT   5 DU POINT   1
  
LA VALEUR MINIMALE DE VMIS    EST:  0.32177E-01
OBTENUE DANS LA MAILLE M2838   AU SOUS_POINT   1 DU POINT   1

 --------------------------------------------------------------------------------
 CONCEPT   vmbot     DE TYPE  EVOL_NOLI                                         


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  SIEQ_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
 LA VALEUR MAXIMALE DE VMIS             EST     3.550E+02 EN    3 MAILLE(S) : M5504   
 LA VALEUR MINIMALE DE VMIS             EST     3.218E-02 EN    1 MAILLE(S) : M2838   

 --------------------------------------------------------------------------------
 CONCEPT   vmtop     DE TYPE  EVOL_NOLI                                         


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  SIEQ_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
 LA VALEUR MAXIMALE DE VMIS             EST     3.550E+02 EN    3 MAILLE(S) : M5504   
 LA VALEUR MINIMALE DE VMIS             EST     3.218E-02 EN    1 MAILLE(S) : M2838   

 --------------------------------------------------------------------------------
 ASTER 14.04.00 CONCEPT results CALCULE LE 22/10/2021 A 12:30:09 DE TYPE         
 EVOL_NOLI                                                                       


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  VARI_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
LA VALEUR MAXIMALE DE V1      EST:  0.94364E-03
OBTENUE DANS LA MAILLE M5505   AU SOUS_POINT   1 DU POINT   2
  
LA VALEUR MINIMALE DE V1      EST:  0.00000E+00
OBTENUE DANS LA MAILLE M969    AU SOUS_POINT   1 DU POINT   1

 --------------------------------------------------------------------------------
 CONCEPT   vmbot     DE TYPE  EVOL_NOLI                                         


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  VARI_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
 LA VALEUR MAXIMALE DE V1               EST     9.436E-04 EN    1 MAILLE(S) : M5505   
 LA VALEUR MINIMALE DE V1               EST     0.000E+00 EN **** MAILLE(S) : M969    

 --------------------------------------------------------------------------------
 CONCEPT   vmtop     DE TYPE  EVOL_NOLI                                         


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  VARI_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
 LA VALEUR MAXIMALE DE V1               EST     7.981E-04 EN    1 MAILLE(S) : M5505   
 LA VALEUR MINIMALE DE V1               EST     0.000E+00 EN **** MAILLE(S) : M969    

- but not for COQUE_3D,

 --------------------------------------------------------------------------------
 ASTER 14.04.00 CONCEPT results CALCULE LE 22/10/2021 A 12:18:54 DE TYPE         
 EVOL_NOLI                                                                       


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  SIEQ_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
LA VALEUR MAXIMALE DE VMIS    EST:  0.44450E+03
OBTENUE DANS LA MAILLE M3039   AU SOUS_POINT   9 DU POINT   3
  
LA VALEUR MINIMALE DE VMIS    EST:  0.12453E+00
OBTENUE DANS LA MAILLE M2843   AU SOUS_POINT   5 DU POINT   1

 --------------------------------------------------------------------------------
 CONCEPT   vmbot     DE TYPE  EVOL_NOLI                                         


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  SIEQ_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
 LA VALEUR MAXIMALE DE VMIS             EST     4.223E+02 EN    1 MAILLE(S) : M2918   
 LA VALEUR MINIMALE DE VMIS             EST     2.328E+00 EN    1 MAILLE(S) : M2843   

 --------------------------------------------------------------------------------
 CONCEPT   vmtop     DE TYPE  EVOL_NOLI                                         


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  SIEQ_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
 LA VALEUR MAXIMALE DE VMIS             EST     4.445E+02 EN    1 MAILLE(S) : M3039   
 LA VALEUR MINIMALE DE VMIS             EST     3.376E-01 EN    1 MAILLE(S) : M1920   

 --------------------------------------------------------------------------------
 ASTER 14.04.00 CONCEPT results CALCULE LE 22/10/2021 A 12:18:54 DE TYPE         
 EVOL_NOLI                                                                       


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  VARI_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
LA VALEUR MAXIMALE DE V1      EST:  0.10154E-02
OBTENUE DANS LA MAILLE M5504   AU SOUS_POINT   1 DU POINT   5
  
LA VALEUR MINIMALE DE V1      EST:  0.00000E+00
OBTENUE DANS LA MAILLE M969    AU SOUS_POINT   1 DU POINT   1

 --------------------------------------------------------------------------------
 CONCEPT   vmbot     DE TYPE  EVOL_NOLI                                         


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  VARI_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
 LA VALEUR MAXIMALE DE V1               EST     1.015E-03 EN    1 MAILLE(S) : M5504   
 LA VALEUR MINIMALE DE V1               EST     0.000E+00 EN **** MAILLE(S) : M969    

 --------------------------------------------------------------------------------
 CONCEPT   vmtop     DE TYPE  EVOL_NOLI                                         


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  VARI_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
 LA VALEUR MAXIMALE DE V1               EST     8.535E-04 EN    1 MAILLE(S) : M5504   
 LA VALEUR MINIMALE DE V1               EST     0.000E+00 EN **** MAILLE(S) : M969    

I have also problem with export VARI_ELGA V1 from results from shell layers to med file
(MISES export was successfully):

JDC.py : ERREUR WITH THE EXECUTION - INTERRUPTION
>> JDC.py: DEBUT CR of execution of JDC in MIXTE
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   ! <S> Exception user raised but not interceptee.                                  !
   ! The bases are fermees.                                                          !
   ! Type of the exception: error                                                    !
   !                                                                                 !
   !  Erreur utilisateur :                                                           !
   !    Le programme a besoin d'accéder au champ COMPORTEMENT de la structure de     !
   ! données résultat vmbot                                                          !
   !    pour le                                                                      !
   ! NUME_ORDRE: 0                                                                   !
   !    Mais ce champ n'existe pas dans la structure de données fournie.             !
   !    On ne peut pas continuer.                                                    !
   !                                                                                 !
   !  Risques & conseils                                                             !
   ! :                                                                               !
   !  Vérifiez que la structure de données vmbot est bien celle qu'il faut utiliser. !
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
fine CR of execution of JDC in MIXTE

The med.file is ok and i can open it but there are no results of VARI_ELGA.

Regards,
Tomek

COMM files:
- dkt:

#####Units: mm, N, N/mm2=MPa

#####Parameters:
###Shell:
T=12.5		#Shell thickness,
LN=3		#Number of layers,

###Load:
N=200e3		#Total force
b=75		#Width of load area
l=200		#Length of load area
A=b*l		#Area of load
q=N/A		#Surface load

###Iterations:
T1=5		#Number of iterations

###Materials:
#Steel:
Es=210e+3
Ep=2.1
fy=355

######Start: 
DEBUT(LANG='EN')

#####Mesh:
###Loading:
mesh = LIRE_MAILLAGE(FORMAT='MED',
                     UNITE=20)
###Groups:
mesh = DEFI_GROUP(reuse=mesh,
				  MAILLAGE=mesh,
				  CREA_GROUP_MA=(_F(TOUT="OUI",
									NOM='lgma')),
				  CREA_GROUP_NO=(_F(GROUP_MA='bottom'),),
				  )

#####Model:
model = AFFE_MODELE(AFFE=(_F(MODELISATION='DKT',
							PHENOMENE='MECANIQUE',
                            TOUT="OUI")),
					MAILLAGE=mesh)
					
#####Elements:
elemprop = AFFE_CARA_ELEM(COQUE=_F(EPAIS=T,
                                   GROUP_MA='lgma',
								   COQUE_NCOU=LN),		
                          MODELE=model)
						  
#####Materials:
###Properties:						
mater = DEFI_MATERIAU(ELAS=_F(E=Es,
                              NU=0.3),
					  ECRO_LINE=_F(D_SIGM_EPSI=Ep,
								   SY=fy))
###Assign:
fieldmat = AFFE_MATERIAU(AFFE=_F(MATER=mater,
                                 TOUT="OUI"),
                         MAILLAGE=mesh)

#####Loads and BC:
###Loads:
load = AFFE_CHAR_MECA(FORCE_COQUE=_F(FZ=-q,
                                     GROUP_MA='top'),
                      MODELE=model)
###BC's:					  
support = AFFE_CHAR_MECA(DDL_IMPO=_F(DX=0.0, DY=0.0, DZ=0.0,
                                     GROUP_MA='bottom'),
                         MODELE=model)

#####Functions:
###Loads:
loadf = DEFI_FONCTION(NOM_PARA='INST',
					  VALE=(0.0,  0.0,
							T1,   1.0),
					  PROL_DROITE='EXCLU',
					  PROL_GAUCHE='EXCLU')

#####Lists
###Iterations:
listiter = DEFI_LIST_REEL(DEBUT=0.0,
						  INTERVALLE=_F(JUSQU_A=T1,
										NOMBRE=T1))

#####Calculation:
###Nonlinear:
results = STAT_NON_LINE(CARA_ELEM=elemprop,
                        CHAM_MATER=fieldmat,
                        EXCIT=(_F(CHARGE=load,
								  FONC_MULT=loadf),
                               _F(CHARGE=support)),
						COMPORTEMENT=_F(RELATION='VMIS_ISOT_LINE',
										DEFORMATION='PETIT',
										TOUT='OUI'),
                        INCREMENT=_F(LIST_INST=listiter),
						CONVERGENCE=_F(ITER_GLOB_MAXI=20,
                                       ITER_GLOB_ELAS=25),
						NEWTON=_F(REAC_INCR=1,
                                  MATRICE='TANGENTE',
                                  REAC_ITER=1,
                                  REAC_ITER_ELAS=1,
                                  MATR_RIGI_SYME='NON'),
                        MODELE=model)

#####Prepare of results:
###Global:
results = CALC_CHAMP(reuse=results,
					 CRITERES=('SIEQ_ELGA',),
                     RESULTAT=results)
###Bottom layer:
vmbot = POST_CHAMP(EXTR_COQUE=_F(NIVE_COUCHE='INF',
                                 NOM_CHAM=('SIEQ_ELGA','VARI_ELGA'),
                                 NUME_COUCHE=1),
                   RESULTAT=results)
###Top layer:
vmtop = POST_CHAMP(EXTR_COQUE=_F(NIVE_COUCHE='SUP',
                                 NOM_CHAM=('SIEQ_ELGA','VARI_ELGA'),
                                 NUME_COUCHE=LN),
                   RESULTAT=results)
				   
#####Results:
###Time points for results:
listresu = DEFI_LIST_ENTI(VALE=(T1))

###Files units:
#MED=80
#RESU=19
#TABLE=8

###Results:
IMPR_RESU(FORMAT='RESULTAT',
          RESU=(
				#MISES:
				_F(RESULTAT=results,		#Global,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='SIEQ_ELGA',
                   NOM_CMP='VMIS',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),
				_F(RESULTAT=vmbot,			#Bottom,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='SIEQ_ELGA',
                   NOM_CMP='VMIS',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),
				_F(RESULTAT=vmtop,			#Top,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='SIEQ_ELGA',
                   NOM_CMP='VMIS',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),
				#Cumulated plastic deformation:   
				_F(RESULTAT=results,		#Global,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='VARI_ELGA',
                   NOM_CMP='V1',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),
				_F(RESULTAT=vmbot,			#Bottom,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='VARI_ELGA',
                   NOM_CMP='V1',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),
				_F(RESULTAT=vmtop,			#Top,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='VARI_ELGA',
                   NOM_CMP='V1',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),
					   
				),
          UNITE=19)
###MED:
IMPR_RESU(FORMAT='MED',
          RESU=(_F(RESULTAT=results,
                   NOM_CHAM='DEPL',
				   NOM_CMP=('DX','DY','DZ'),
                   ),
				_F(RESULTAT=results,
				   NOM_CHAM='SIEQ_ELGA',
                   NOM_CMP='VMIS',
                   ),				
				_F(RESULTAT=vmbot,
				   NOM_CHAM='SIEQ_ELGA',
                   NOM_CMP='VMIS',
                   ),
                _F(RESULTAT=vmtop,
				   NOM_CHAM='SIEQ_ELGA',
                   NOM_CMP='VMIS',
                   ),
                _F(RESULTAT=results,
				   NOM_CHAM='VARI_ELGA',
                   NOM_CMP='V1',
                   ),
                _F(RESULTAT=vmbot,
				   NOM_CHAM='VARI_ELGA',
                   NOM_CMP='V1',
                   ),
                _F(RESULTAT=vmtop,
				   NOM_CHAM='VARI_ELGA',
                   NOM_CMP='V1',
                   ),				   
				),
          UNITE=80)

#End:
FIN()			

- coque_3d:

#####Units: mm, N, N/mm2=MPa

#####Parameters:
###Shell:
T=12.5		#Shell thickness,
LN=3		#Number of layers,

###Load:
N=200e3		#Total force
b=75		#Width of load area
l=200		#Length of load area
A=b*l		#Area of load
q=N/A		#Surface load

###Iterations:
T1=5		#Number of iterations

###Materials:
#Steel:
Es=210e+3
Ep=2.1
fy=355

######Start: 
DEBUT(LANG='EN')

#####Mesh:
###Loading:
mesh = LIRE_MAILLAGE(FORMAT='MED',
                     UNITE=20)
###Groups:
mesh = DEFI_GROUP(reuse=mesh,
				  MAILLAGE=mesh,
				  CREA_GROUP_MA=(_F(TOUT="OUI",
									NOM='lgma')),
				  CREA_GROUP_NO=(_F(GROUP_MA='bottom'),),
				  )

#####Model:
model = AFFE_MODELE(AFFE=(_F(MODELISATION='COQUE_3D',
							PHENOMENE='MECANIQUE',
                            TOUT="OUI")),
					MAILLAGE=mesh)
					
#####Elements:
elemprop = AFFE_CARA_ELEM(COQUE=_F(EPAIS=T,
                                   GROUP_MA='lgma',
								   COQUE_NCOU=LN),		
                          MODELE=model)
						  
#####Materials:
###Properties:						
mater = DEFI_MATERIAU(ELAS=_F(E=Es,
                              NU=0.3),
					  ECRO_LINE=_F(D_SIGM_EPSI=Ep,
								   SY=fy))
###Assign:
fieldmat = AFFE_MATERIAU(AFFE=_F(MATER=mater,
                                 TOUT="OUI"),
                         MAILLAGE=mesh)

#####Loads and BC:
###Loads:
load = AFFE_CHAR_MECA(FORCE_COQUE=_F(FZ=-q,
                                     GROUP_MA='top'),
                      MODELE=model)
###BC's:					  
support = AFFE_CHAR_MECA(DDL_IMPO=_F(DX=0.0, DY=0.0, DZ=0.0,
                                     GROUP_MA='bottom'),
                         MODELE=model)

#####Functions:
###Loads:
loadf = DEFI_FONCTION(NOM_PARA='INST',
					  VALE=(0.0,  0.0,
							T1,   1.0),
					  PROL_DROITE='EXCLU',
					  PROL_GAUCHE='EXCLU')

#####Lists
###Iterations:
listiter = DEFI_LIST_REEL(DEBUT=0.0,
						  INTERVALLE=_F(JUSQU_A=T1,
										NOMBRE=T1))

#####Calculation:
###Nonlinear:
results = STAT_NON_LINE(CARA_ELEM=elemprop,
                        CHAM_MATER=fieldmat,
                        EXCIT=(_F(CHARGE=load,
								  FONC_MULT=loadf),
                               _F(CHARGE=support)),
						COMPORTEMENT=_F(RELATION='VMIS_ISOT_LINE',
										DEFORMATION='PETIT',
										TOUT='OUI'),
                        INCREMENT=_F(LIST_INST=listiter),
						CONVERGENCE=_F(ITER_GLOB_MAXI=20,
                                       ITER_GLOB_ELAS=25),
						NEWTON=_F(REAC_INCR=1,
                                  MATRICE='TANGENTE',
                                  REAC_ITER=1,
                                  REAC_ITER_ELAS=1,
                                  MATR_RIGI_SYME='NON'),
                        MODELE=model)

#####Prepare of results:
###Global:
results = CALC_CHAMP(reuse=results,
					 CRITERES=('SIEQ_ELGA',),
                     RESULTAT=results)
###Bottom layer:
vmbot = POST_CHAMP(EXTR_COQUE=_F(NIVE_COUCHE='INF',
                                 NOM_CHAM=('SIEQ_ELGA','VARI_ELGA'),
                                 NUME_COUCHE=1),
                   RESULTAT=results)
###Top layer:
vmtop = POST_CHAMP(EXTR_COQUE=_F(NIVE_COUCHE='SUP',
                                 NOM_CHAM=('SIEQ_ELGA','VARI_ELGA'),
                                 NUME_COUCHE=LN),
                   RESULTAT=results)
				   
#####Results:
###Time points for results:
listresu = DEFI_LIST_ENTI(VALE=(T1))

###Files units:
#MED=80
#RESU=19
#TABLE=8

###Results:
IMPR_RESU(FORMAT='RESULTAT',
          RESU=(
				#MISES:
				_F(RESULTAT=results,		#Global,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='SIEQ_ELGA',
                   NOM_CMP='VMIS',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),
				_F(RESULTAT=vmbot,			#Bottom,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='SIEQ_ELGA',
                   NOM_CMP='VMIS',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),
				_F(RESULTAT=vmtop,			#Top,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='SIEQ_ELGA',
                   NOM_CMP='VMIS',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),
				#Cumulated plastic deformation:   
				_F(RESULTAT=results,		#Global,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='VARI_ELGA',
                   NOM_CMP='V1',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),
				_F(RESULTAT=vmbot,			#Bottom,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='VARI_ELGA',
                   NOM_CMP='V1',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),
				_F(RESULTAT=vmtop,			#Top,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='VARI_ELGA',
                   NOM_CMP='V1',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),
					   
				),
          UNITE=19)
#End:
FIN()			

Attachments:
rhs.7z, Size: 56.57 KiB, Downloads: 26

Offline

#4 2021-10-22 15:08:00

AsterO'dactyle
Administrator
Registered: 2007-11-29
Posts: 390

Re: Shells - DKT - Plastsicity - (In)correct Mises

Hello,

To print VARI_ELGA after a POST_CHAMP, use IMPR_NOM_VARI='NON' in IMPR_RESU command.
It's a known bug will been corrected in next version


Code_Asterの開発者

Offline

#5 2021-10-22 15:46:42

tpa
Member
From: Denmark
Registered: 2009-04-14
Posts: 61

Re: Shells - DKT - Plastsicity - (In)correct Mises

Good that it works for DKT modelisation. I do not use Coque3D a lot, as the results require some processing above my level of ability to be handled by PARAVIS.

Offline

#6 2021-10-22 16:14:53

TomekBeginner
Member
Registered: 2021-10-20
Posts: 7

Re: Shells - DKT - Plastsicity - (In)correct Mises

Hello,

AsterO'dactyle wrote:

To print VARI_ELGA after a POST_CHAMP, use IMPR_NOM_VARI='NON' in IMPR_RESU command.

Unfortunately it doesn't work either:

ERREUR WITH THE SYNTACTIC VERIFICATION - INTERRUPTION
>> JDC.py: DEBUT DEBUT CR validation : fort.1
   Etape : IMPR_RESU    ligne : 233    fichier : 'fort.1'
      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
      ! Unknown key words: IMPR_NOM_VARI !
      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   Fin Etape : IMPR_RESU
FIN CR validation :fort.1

In v14 documentation of IMPR_RESU there is no keyword IMPR_NOM_VARI.

Regards,
Tomek

Offline

#7 2021-10-22 17:00:39

AsterO'dactyle
Administrator
Registered: 2007-11-29
Posts: 390

Re: Shells - DKT - Plastsicity - (In)correct Mises

Look again at the documentation !


Code_Asterの開発者

Offline

#8 2021-10-22 19:30:59

TomekBeginner
Member
Registered: 2021-10-20
Posts: 7

Re: Shells - DKT - Plastsicity - (In)correct Mises

I just find out there are several documents for IMPR_RESU.
Aside I still don't know how to obtain correct MISES value from COQUE_3D everything is ok (for) now.

@AsterO'dactyle
Thank you very much and sorry for trouble.

Offline

#9 2021-10-22 21:19:25

tpa
Member
From: Denmark
Registered: 2009-04-14
Posts: 61

Re: Shells - DKT - Plastsicity - (In)correct Mises

Do you use the same mesh for DKT and Coque3D? Element type for Coque3D should be TRIA7 or QUAD9 ie. 2nd order with a center node. One way to do that in case your mesher cant handle these type (e.g. Salome) is to issue a MODI_MAILLAGE in Code Aster. This should be just after you read the mesh.

Last edited by tpa (2021-10-22 21:23:22)

Offline

#10 2021-10-22 21:38:14

TomekBeginner
Member
Registered: 2021-10-20
Posts: 7

Re: Shells - DKT - Plastsicity - (In)correct Mises

No.
For COQUE_3D i've used QUAD9 mesh.
You can download cases files from post #3.

Regards,
Tomek

Offline

#11 2021-10-27 18:19:34

TomekBeginner
Member
Registered: 2021-10-20
Posts: 7

Re: Shells - DKT - Plastsicity - (In)correct Mises

Hello,

I still wonder why there is wrong MISES value in at least one Gauss point of middle layer.
Mesh in attachment (MED3.3).

Regard,
Tomek

COMM:

#####Units: mm, N, N/mm2=MPa

#####Parameters:
###Shell:
T=12.5		#Shell thickness,
LN=3		#Number of layers,

###Load:
N=200e3		#Total force
b=75		#Width of load area
l=200		#Length of load area
A=b*l		#Area of load
q=N/A		#Surface load

###Iterations:
T1=5		#Number of iterations

###Materials:
#Steel:
Es=210e+3
Ep=2.1
fy=355

######Start: 
DEBUT(LANG='EN')

#####Mesh:
###Loading:
mesh = LIRE_MAILLAGE(FORMAT='MED',
                     UNITE=20)
###Groups:
mesh = DEFI_GROUP(reuse=mesh,
				  MAILLAGE=mesh,
				  CREA_GROUP_MA=(_F(TOUT="OUI",
									NOM='lgma')),
				  CREA_GROUP_NO=(_F(GROUP_MA='bottom'),),
				  )

#####Model:
model = AFFE_MODELE(AFFE=(_F(MODELISATION='DKT',
							PHENOMENE='MECANIQUE',
                            TOUT="OUI")),
					MAILLAGE=mesh)
					
#####Elements:
elemprop = AFFE_CARA_ELEM(COQUE=_F(EPAIS=T,
                                   GROUP_MA='lgma',
								   COQUE_NCOU=LN),		
                          MODELE=model)
						  
#####Materials:
###Properties:						
mater = DEFI_MATERIAU(ELAS=_F(E=Es, NU=0.3),
					  ECRO_LINE=_F(D_SIGM_EPSI=Ep,
								   SY=fy))
###Assign:
fieldmat = AFFE_MATERIAU(AFFE=_F(MATER=mater,
                                 TOUT="OUI"),
                         MAILLAGE=mesh)

#####Loads and BC:
###Loads:
load = AFFE_CHAR_MECA(FORCE_COQUE=_F(FZ=-q,
                                     GROUP_MA='top'),
                      MODELE=model)
###BC's:					  
support = AFFE_CHAR_MECA(DDL_IMPO=_F(DX=0.0, DY=0.0, DZ=0.0,
                                     GROUP_MA='bottom'),
                         MODELE=model)

#####Functions:
###Loads:
loadf = DEFI_FONCTION(NOM_PARA='INST',
					  VALE=(0.0,  0.0,
							T1,   1.0),
					  PROL_DROITE='EXCLU',
					  PROL_GAUCHE='EXCLU')

#####Lists
###Iterations:
listiter = DEFI_LIST_REEL(DEBUT=0.0,
						  INTERVALLE=_F(JUSQU_A=T1,
										NOMBRE=T1))

#####Calculation:
###Nonlinear:
results = STAT_NON_LINE(CARA_ELEM=elemprop,
                        CHAM_MATER=fieldmat,
                        EXCIT=(_F(CHARGE=load,
								  FONC_MULT=loadf),
                               _F(CHARGE=support)),
						COMPORTEMENT=_F(RELATION='VMIS_ISOT_LINE',
										DEFORMATION='PETIT',
										TOUT='OUI'),
                        INCREMENT=_F(LIST_INST=listiter),
						CONVERGENCE=_F(ITER_GLOB_MAXI=20,
                                       ITER_GLOB_ELAS=25),
						NEWTON=_F(REAC_INCR=1,
                                  MATRICE='TANGENTE',
                                  REAC_ITER=1,
                                  REAC_ITER_ELAS=1,
                                  MATR_RIGI_SYME='NON'),
                        MODELE=model)

#####Prepare of results:
###Global:
results = CALC_CHAMP(reuse=results,
					 CRITERES=('SIEQ_ELGA',),
                     RESULTAT=results)
###Bottom layer:
vmbot = POST_CHAMP(EXTR_COQUE=_F(NIVE_COUCHE='INF',
                                 NOM_CHAM=('SIEQ_ELGA','VARI_ELGA'),
                                 NUME_COUCHE=1),
                   RESULTAT=results)
###Top layer:
vmtop = POST_CHAMP(EXTR_COQUE=_F(NIVE_COUCHE='SUP',
                                 NOM_CHAM=('SIEQ_ELGA','VARI_ELGA'),
                                 NUME_COUCHE=LN),
                   RESULTAT=results)
				   
###Middle layer:
vmmid = POST_CHAMP(EXTR_COQUE=_F(NIVE_COUCHE='MOY',
                                 NOM_CHAM=('SIEQ_ELGA','VARI_ELGA'),
                                 NUME_COUCHE=2),
                   RESULTAT=results)
				   
#####Results:
###Time points for results:
listresu = DEFI_LIST_ENTI(VALE=(T1))

###Files units:
#MED=80
#RESU=19
#TABLE=8

###Results:
IMPR_RESU(FORMAT='RESULTAT',
          RESU=(
				###MISES###:
				_F(RESULTAT=results,		#Global,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='SIEQ_ELGA',
                   NOM_CMP='VMIS',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),
				_F(RESULTAT=vmbot,			#Bottom,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='SIEQ_ELGA',
                   NOM_CMP='VMIS',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),
				_F(RESULTAT=vmtop,			#Top,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='SIEQ_ELGA',
                   NOM_CMP='VMIS',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),
				_F(RESULTAT=vmmid,			#Mid,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='SIEQ_ELGA',
                   NOM_CMP='VMIS',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),				   
				###Cumulated plastic deformation###:   
				_F(RESULTAT=results,		#Global,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='VARI_ELGA',
                   NOM_CMP='V1',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),
				_F(RESULTAT=vmbot,			#Bottom,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='VARI_ELGA',
                   NOM_CMP='V1',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),
				_F(RESULTAT=vmtop,			#Top,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='VARI_ELGA',
                   NOM_CMP='V1',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),
				_F(RESULTAT=vmmid,			#Mid,
				   LIST_ORDRE=listresu,
				   FORMAT_R='1PE12.3',
                   NOM_CHAM='VARI_ELGA',
                   NOM_CMP='V1',
                   VALE_MAX='OUI',
                   VALE_MIN='OUI'),					   
				),
          UNITE=19)
		  
#End:
FIN()			

RESU:

 --------------------------------------------------------------------------------
 ASTER 14.04.00 CONCEPT results CALCULE LE 27/10/2021 A 19:08:19 DE TYPE         
 EVOL_NOLI                                                                       


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  SIEQ_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
LA VALEUR MAXIMALE DE VMIS    EST:  0.49367E+03
OBTENUE DANS LA MAILLE M2759   AU SOUS_POINT   5 DU POINT   1
  
LA VALEUR MINIMALE DE VMIS    EST:  0.32177E-01
OBTENUE DANS LA MAILLE M2838   AU SOUS_POINT   1 DU POINT   1

 --------------------------------------------------------------------------------
 CONCEPT   vmbot     DE TYPE  EVOL_NOLI                                         


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  SIEQ_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
 LA VALEUR MAXIMALE DE VMIS             EST     3.550E+02 EN    3 MAILLE(S) : M5504   
 LA VALEUR MINIMALE DE VMIS             EST     3.218E-02 EN    1 MAILLE(S) : M2838   

 --------------------------------------------------------------------------------
 CONCEPT   vmtop     DE TYPE  EVOL_NOLI                                         


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  SIEQ_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
 LA VALEUR MAXIMALE DE VMIS             EST     3.550E+02 EN    3 MAILLE(S) : M5504   
 LA VALEUR MINIMALE DE VMIS             EST     3.218E-02 EN    1 MAILLE(S) : M2838   

 --------------------------------------------------------------------------------
 CONCEPT   vmmid     DE TYPE  EVOL_NOLI                                         


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  SIEQ_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
 LA VALEUR MAXIMALE DE VMIS             EST     4.937E+02 EN    1 MAILLE(S) : M2759   
 LA VALEUR MINIMALE DE VMIS             EST     1.171E+00 EN    1 MAILLE(S) : M2978   

 --------------------------------------------------------------------------------
 ASTER 14.04.00 CONCEPT results CALCULE LE 27/10/2021 A 19:08:19 DE TYPE         
 EVOL_NOLI                                                                       


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  VARI_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
LA VALEUR MAXIMALE DE V1      EST:  0.94364E-03
OBTENUE DANS LA MAILLE M5505   AU SOUS_POINT   1 DU POINT   2
  
LA VALEUR MINIMALE DE V1      EST:  0.00000E+00
OBTENUE DANS LA MAILLE M969    AU SOUS_POINT   1 DU POINT   1

 --------------------------------------------------------------------------------
 CONCEPT   vmbot     DE TYPE  EVOL_NOLI                                         


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  VARI_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
 LA VALEUR MAXIMALE DE V1               EST     9.436E-04 EN    1 MAILLE(S) : M5505   
 LA VALEUR MINIMALE DE V1               EST     0.000E+00 EN **** MAILLE(S) : M969    

 --------------------------------------------------------------------------------
 CONCEPT   vmtop     DE TYPE  EVOL_NOLI                                         


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  VARI_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
 LA VALEUR MAXIMALE DE V1               EST     7.981E-04 EN    1 MAILLE(S) : M5505   
 LA VALEUR MINIMALE DE V1               EST     0.000E+00 EN **** MAILLE(S) : M969    

 --------------------------------------------------------------------------------
 CONCEPT   vmmid     DE TYPE  EVOL_NOLI                                         


 ======>

 ------>
 CHAMP PAR ELEMENT AUX POINTS DE GAUSS DE NOM SYMBOLIQUE  VARI_ELGA              
 NUMERO D'ORDRE: 5 INST:     5.000E+00                                           
  
 LA VALEUR MAXIMALE DE V1               EST     0.000E+00 EN **** MAILLE(S) : M969    
 LA VALEUR MINIMALE DE V1               EST     0.000E+00 EN **** MAILLE(S) : M969    

Attachments:
mesh_linear.mmed, Size: 611.52 KiB, Downloads: 26

Offline