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

You are not logged in.

#1 2021-08-11 10:57:01

Perceval Zhou
Member
Registered: 2021-08-11
Posts: 12

Comment faire un analyse simple de fatigue sur tuyau?

Monsieurs,
J'ai rencontré un problème lors d'un étude dans lequel je voudrais réaliser un analyse de fatigue de type de Wohler en utilisant CALC_FATIGUE, et je me réfère sur l'exemple szlz105a pour l'édition du fichier .comm, mais il appaît toujous des erreurs dans les executions.
J'ai créé le maillage de la structure de tuyau comme un arête avec deux coudes( en fait c'est le maillage de l'exemple sdlx02d), et j'utilise la modélisation de TUYAU_6M, et les autres commandes sont tous similaires aux celles de szlz105a.
Voici mon fichier .comm et .med, j'utilise Salome-Meca2018 de version Windows.


Attachments:
ssll101d_mesh1.med, Size: 21.4 KiB, Downloads: 130

Offline

#2 2021-08-11 10:58:41

Perceval Zhou
Member
Registered: 2021-08-11
Posts: 12

Re: Comment faire un analyse simple de fatigue sur tuyau?

Merci de m'en débarrasser!


Attachments:
szlz105a_2.comm, Size: 3.75 KiB, Downloads: 149

Offline

#3 2021-08-15 09:05:03

Perceval Zhou
Member
Registered: 2021-08-11
Posts: 12

Re: Comment faire un analyse simple de fatigue sur tuyau?

Y a quelqu'un peut m'aider?

Offline

#4 2021-08-15 15:38:48

miib
Member
From: Near Goethe's Palm.
Registered: 2021-03-23
Posts: 248

Re: Comment faire un analyse simple de fatigue sur tuyau?

Perceval Zhou wrote:

Comment faire un analyse simple de fatigue sur tuyau?

1) draw a cylinder

2) define 3 geometrical groups : " fix "= a base of cylinder  ; "load" = the other base : "c" = the volume.

3) generate mesh and groups of mesh from geometry.

4) open aster module

5) add a new stage

6) select,  copy ( CTRL+C) and paste ( CTRL+V) this code in the stage :


DEBUT(LANG='EN')


tstep = 0.5

mesh = LIRE_MAILLAGE(identifier='1:1',
                     FORMAT='MED',
                     UNITE=20)

model = AFFE_MODELE(identifier='2:1',
                    AFFE=_F(MODELISATION=('3D', ),
                            PHENOMENE='MECANIQUE',
                            TOUT='OUI'),
                    MAILLAGE=mesh)

ys = DEFI_CONSTANTE(identifier='3:1',
                    VALE=210000000000.0)

ps = DEFI_CONSTANTE(identifier='4:1',
                    VALE=0.3)

ds = DEFI_CONSTANTE(identifier='5:1',
                    VALE=8000.0)

wholer = DEFI_FONCTION(identifier='6:1',
                       INTERPOL=('LIN', 'LOG'),
                       NOM_PARA='SIGM',
                       PROL_DROITE='LINEAIRE',
                       PROL_GAUCHE='CONSTANT',
                       VALE=(86200000.0, 1000000.0, 114000000.0, 200000.0, 138000000.0, 100000.0, 214000000.0, 20000.0, 262000000.0, 10000.0, 441000000.0, 2000.0, 1070000000.0, 200.0, 1410000000.0, 100.0, 1900000000.0, 50.0, 2830000000.0, 20.0, 4000000000.0, 10.0))

steel = DEFI_MATERIAU(identifier='7:1',
                      ECRO_LINE=_F(D_SIGM_EPSI=1100000000.0,
                                   SY=700000000.0),
                      ELAS_FO=_F(E=ys,
                                 NU=ps,
                                 RHO=ds),
                      FATIGUE=_F(WOHLER=wholer),
                      RCCM=_F(SU=900000000.0))

fieldmat = AFFE_MATERIAU(identifier='8:1',
                         AFFE=_F(MATER=(steel, ),
                                 TOUT='OUI'),
                         MODELE=model)

loadf = DEFI_FONCTION(identifier='9:1',
                      NOM_PARA='INST',
                      PROL_DROITE='CONSTANT',
                      PROL_GAUCHE='EXCLU',
                      VALE=(1.0, 70000000.0, 2.0, -70000000.0, 3.0, 70000000.0, 4.0, -70000000.0, 5.0, 70000000.0, 6.0, -70000000.0, 7.0, 70000000.0, 8.0, -70000000.0, 9.0, 70000000.0, 10.0, -70000000.0))

listr = DEFI_LIST_REEL(identifier='10:1',
                       DEBUT=1.0,
                       INTERVALLE=_F(JUSQU_A=10.0,
                                     PAS=tstep))

times = DEFI_LIST_INST(identifier='11:1',
                       DEFI_LIST=_F(LIST_INST=listr),
                       METHODE='AUTO')

load = AFFE_CHAR_MECA(identifier='12:1',
                      FACE_IMPO=_F(DX=0.0,
                                   DY=0.0,
                                   DZ=0.0,
                                   GROUP_MA=('fix', )),
                      MODELE=model)

loadd = AFFE_CHAR_MECA_F(identifier='13:1',
                         FORCE_FACE=_F(FY=loadf,
                                       GROUP_MA=('load', )),
                         MODELE=model)

resnonl = DYNA_NON_LINE(identifier='14:1',
                        CHAM_MATER=fieldmat,
                        COMPORTEMENT=_F(DEFORMATION='SIMO_MIEHE',
                                        RELATION='VMIS_ISOT_LINE'),
                        CONVERGENCE=_F(ITER_GLOB_MAXI=50),
                        EXCIT=(_F(CHARGE=load),
                               _F(CHARGE=loadd)),
                        INCREMENT=_F(LIST_INST=times),
                        MODELE=model,
                        NEWTON=_F(MATRICE='TANGENTE',
                                  PREDICTION='TANGENTE'),
                        SCHEMA_TEMPS=_F(FORMULATION='DEPLACEMENT',
                                        SCHEMA='HHT'),
                        SOLVEUR=_F(METHODE='MUMPS'))

resnonl = CALC_CHAMP(identifier='15:1',
                     reuse=resnonl,
                     CRITERES=('SIEQ_ELGA', ),
                     RESULTAT=resnonl)

fati = CALC_FATIGUE(identifier='16:1',
                    DOMMAGE='WOHLER',
                    HISTOIRE=_F(EQUI_GD='VMIS_SG',
                                RESULTAT=resnonl),
                    MATER=steel,
                    OPTION='DOMA_ELGA_SIGM',
                    TYPE_CALCUL='CUMUL_DOMMAGE')

IMPR_RESU(identifier='17:1',
          FORMAT='MED',
          RESU=(_F(RESULTAT=resnonl),
                _F(CHAM_GD=fati,
                   NOM_CHAM_MED=('DOMA', ),
                   NOM_CMP=('DOMA', ))),
          UNITE=2)

FIN()

7) in the tree pasted delete deprecated

8) reassigning the mesh

9) reassigning the directory of output results.

10 ) run.


SE,Dio è il creatore di tutto, e tu credi in Dio, crederai anche che le formule della matematica e della fisica sono "linguaggi" che consentono di spiegare ciò che è stato creato da Dio (la natura) e che tali "linguaggi" sono essi stessi creazioni di Dio. Sarei stupito del contrario.

Offline

#5 2021-08-16 08:27:27

miib
Member
From: Near Goethe's Palm.
Registered: 2021-03-23
Posts: 248

Re: Comment faire un analyse simple de fatigue sur tuyau?

miib wrote:

1) draw a cylinder

the procedure is the same drawing a tube (or a plate with hole)


Attachments:
domage.jpg, Size: 122.04 KiB, Downloads: 128

SE,Dio è il creatore di tutto, e tu credi in Dio, crederai anche che le formule della matematica e della fisica sono "linguaggi" che consentono di spiegare ciò che è stato creato da Dio (la natura) e che tali "linguaggi" sono essi stessi creazioni di Dio. Sarei stupito del contrario.

Offline

#6 2021-08-16 09:34:42

miib
Member
From: Near Goethe's Palm.
Registered: 2021-03-23
Posts: 248

Re: Comment faire un analyse simple de fatigue sur tuyau?

the procedure is the same if you use STAT_NON_LINE


SE,Dio è il creatore di tutto, e tu credi in Dio, crederai anche che le formule della matematica e della fisica sono "linguaggi" che consentono di spiegare ciò che è stato creato da Dio (la natura) e che tali "linguaggi" sono essi stessi creazioni di Dio. Sarei stupito del contrario.

Offline

#7 2021-08-17 04:58:49

Perceval Zhou
Member
Registered: 2021-08-11
Posts: 12

Re: Comment faire un analyse simple de fatigue sur tuyau?

miib,
thank you for your explication, but this command file still won't work with my pipe model, maybe I should learn more from files of non-linear analysis..

Offline

#8 2021-08-17 05:35:08

miib
Member
From: Near Goethe's Palm.
Registered: 2021-03-23
Posts: 248

Re: Comment faire un analyse simple de fatigue sur tuyau?

Perceval Zhou wrote:

miib,
thank you for your explication, but this command file still won't work with my pipe model, maybe I should learn more from files of non-linear analysis..

can you show the message error of the tipe like attached file ??? :-)


SE,Dio è il creatore di tutto, e tu credi in Dio, crederai anche che le formule della matematica e della fisica sono "linguaggi" che consentono di spiegare ciò che è stato creato da Dio (la natura) e che tali "linguaggi" sono essi stessi creazioni di Dio. Sarei stupito del contrario.

Offline

#9 2021-08-17 05:36:09

miib
Member
From: Near Goethe's Palm.
Registered: 2021-03-23
Posts: 248

Re: Comment faire un analyse simple de fatigue sur tuyau?

attachment


Attachments:
message, Size: 41.99 KiB, Downloads: 137

SE,Dio è il creatore di tutto, e tu credi in Dio, crederai anche che le formule della matematica e della fisica sono "linguaggi" che consentono di spiegare ciò che è stato creato da Dio (la natura) e che tali "linguaggi" sono essi stessi creazioni di Dio. Sarei stupito del contrario.

Offline

#10 2021-08-17 08:34:37

Perceval Zhou
Member
Registered: 2021-08-11
Posts: 12

Re: Comment faire un analyse simple de fatigue sur tuyau?

miib wrote:
Perceval Zhou wrote:

miib,
thank you for your explication, but this command file still won't work with my pipe model, maybe I should learn more from files of non-linear analysis..

can you show the message error of the tipe like attached file ??? :-)

Yes, it shows that the modelisation is not compatible with the choice "SIMO_MIEHE"


Attachments:
message, Size: 44.76 KiB, Downloads: 136

Offline

#11 2021-08-17 09:31:12

miib
Member
From: Near Goethe's Palm.
Registered: 2021-03-23
Posts: 248

Re: Comment faire un analyse simple de fatigue sur tuyau?

miib wrote:

1) draw a cylinder

2) define 3 geometrical groups : " fix "= a base of cylinder  ; "load" = the other base : "c" = the volume.

3) generate mesh and groups of mesh from geometry.

miib wrote:

the procedure is the same drawing a tube (or a plate with hole)

miib wrote:

mesh = LIRE_MAILLAGE(identifier='1:1',
                     FORMAT='MED',
                     UNITE=20)

model = AFFE_MODELE(identifier='2:1',
                    AFFE=_F(MODELISATION=('3D', ),
                            PHENOMENE='MECANIQUE',
                            TOUT='OUI'),
                    MAILLAGE=mesh)

p.s. naturally you must adapt your loads to your case study


SE,Dio è il creatore di tutto, e tu credi in Dio, crederai anche che le formule della matematica e della fisica sono "linguaggi" che consentono di spiegare ciò che è stato creato da Dio (la natura) e che tali "linguaggi" sono essi stessi creazioni di Dio. Sarei stupito del contrario.

Offline

#12 2021-08-17 11:26:18

miib
Member
From: Near Goethe's Palm.
Registered: 2021-03-23
Posts: 248

Re: Comment faire un analyse simple de fatigue sur tuyau?

attached here an example of capability


Attachments:
post_process_screen_shot.zip, Size: 250.9 KiB, Downloads: 151

SE,Dio è il creatore di tutto, e tu credi in Dio, crederai anche che le formule della matematica e della fisica sono "linguaggi" che consentono di spiegare ciò che è stato creato da Dio (la natura) e che tali "linguaggi" sono essi stessi creazioni di Dio. Sarei stupito del contrario.

Offline