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

You are not logged in.

#1 2018-11-06 05:45:08

Ceduardolino
Member
Registered: 2017-03-19
Posts: 17

Calc_Fatigue number of total cycles!!!!!!!!!!!!!

Help me, I can not calculate the total number of cycles on this body of evidence.
How to generate an image in the post-processor in Paravis, from this process of fatigue, I can not.
Help me please. The mesh is attached.
Thank you!!

DEBUT();

EM  = 73000.0
NU = 0.33

a = 6.5
ce = 34.0

Pressao = 66.43

F_WOHLER=DEFI_FONCTION(NOM_PARA='SIGM',
                        PROL_DROITE='LINEAIRE',
                        PROL_GAUCHE='LINEAIRE',
                        VALE=(500.,  10.E6,
                             1000.,     0.,
			      ),                                      
                        TITRE='FONCTION DE WOHLER');
								   
								   
ALU=DEFI_MATERIAU(ELAS=_F(E=EM,
                         NU=NU,
                         ALPHA=0.,),
                         ECRO_LINE=_F(D_SIGM_EPSI=17.5,
                         SY=503.,),
			 FATIGUE=_F(WOHLER=F_WOHLER,),);


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


MAIL=MODI_MAILLAGE(reuse =MAIL,
                   MAILLAGE=MAIL,
                   ORIE_PEAU_2D=_F(GROUP_MA=('S_UP','S_DW',),),);

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

CRACK=DEFI_FISS_XFEM(MAILLAGE=MAIL,
                     DEFI_FISS=_F(FORM_FISS='DEMI_DROITE',
                                 PFON=(ce+a,0.0,0.0,),
                                 DTAN=(1.,0.0,0.0,),),);


MODEX=MODI_MODELE_XFEM(MODELE_IN=MODE,
                       FISSURE=CRACK,);

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


CONST=AFFE_CHAR_MECA(MODELE=MODEX,
                     DDL_IMPO=(
                        _F(GROUP_NO='DX_0',
                           DX=0.0,),
                        _F(GROUP_NO='DY_0',
                          DY=0.0,),
                        ),
                     PRES_REP=(
                        _F(GROUP_MA='S_UP',
                           PRES=Pressao,),
                        _F(GROUP_MA='S_DW',
                           PRES=Pressao,),
                        ),);


L_INST=DEFI_LIST_REEL(DEBUT=0.,
                      INTERVALLE=(_F(JUSQU_A=1.0,NOMBRE=1),
                                  ),
                      );

RESU=STAT_NON_LINE(MODELE=MODEX,
                    CHAM_MATER=MATE,
                    EXCIT=(_F(CHARGE=CONST,),),
                    COMPORTEMENT=_F(RELATION='ELAS_VMIS_LINE'),
                    INCREMENT=_F(LIST_INST=L_INST),
                    NEWTON=_F(REAC_ITER=1,),
                    );

G_1=CALC_G(THETA=_F(FISSURE=CRACK,
                    R_INF = 0.1,
                    R_SUP = 5.0),
           RESULTAT=RESU,
           );
                  
IMPR_TABLE(TABLE=G_1,);

MA_XFEM=POST_MAIL_XFEM(MODELE=MODEX,);

MOD_VISU=AFFE_MODELE(MAILLAGE=MA_XFEM,
                     AFFE=_F(TOUT='OUI',
                             PHENOMENE='MECANIQUE',
                             MODELISATION='C_PLAN',),);

RES_XFEM=POST_CHAM_XFEM(MODELE_VISU=MOD_VISU,
                        RESULTAT=RESU,);

   
RES_XFEM=CALC_CHAMP(reuse=RES_XFEM,
                    RESULTAT=RES_XFEM,
                    CONTRAINTE=('SIGM_NOEU'),
                    FORCE='FORC_NODA',
                    CRITERES=('SIEQ_ELNO','SIEQ_NOEU','SIEQ_ELGA'),
                );

IMPR_RESU(FORMAT='MED', UNITE=80,
          RESU=(
              _F(MAILLAGE=MAIL,
                   RESULTAT=RES_XFEM,
                   NOM_CHAM='SIEQ_NOEU',
                   NOM_CMP=('VMIS','VMIS_SG'),
                ),
             _F(MAILLAGE=MAIL,
                   RESULTAT=RES_XFEM,
                   NOM_CHAM='DEPL',
                ),
          ),
)

#IMPR_RESU(FORMAT='MED',RESU=_F(RESULTAT=RES_XFEM))

 
FADIGA=CALC_FATIGUE(TYPE_CALCUL='CUMUL_DOMMAGE',
		      OPTION='DOMA_ELNO_SIGM',
                      HISTOIRE=_F(RESULTAT=RES_XFEM,
                                  EQUI_GD='VMIS_SG',),
                      DOMMAGE='WOHLER',
                      MATER=ALU,
                      INFO=2,);	     
				
IMPR_RESU(FORMAT='MED', UNITE=81,
          RESU=_F(CHAM_GD=FADIGA),);

FIN();

Thank you!!

Carlos Eduardo


Attachments:
cdp_9.mmed, Size: 431.82 KiB, Downloads: 13

Offline

#2 2018-11-08 18:55:13

xavierdestev
Member
Registered: 2017-01-12
Posts: 9

Re: Calc_Fatigue number of total cycles!!!!!!!!!!!!!

Filters->Alphabetical->ELNO_MESH to see the result field of CALC_FATIGUE.

FACY_R NBRUP1 for the number of cycles before dislocation of the structure.

Offline

Board footer