Atom topic feed | site map | contact | login | Protection des données personnelles | Powered by FluxBB | réalisation artaban
You are not logged in.
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.
Offline
Merci de m'en débarrasser!
Offline
Y a quelqu'un peut m'aider?
Offline
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.
I ever used code_aster for study.
Code_Aster ??? is the Novacula Occami of the FEM world.
“Non è bene cercare di fermare il progresso della conoscenza. L'ignoranza non è mai meglio della conoscenza.” Enrico Fermi.
Offline
1) draw a cylinder
the procedure is the same drawing a tube (or a plate with hole)
I ever used code_aster for study.
Code_Aster ??? is the Novacula Occami of the FEM world.
“Non è bene cercare di fermare il progresso della conoscenza. L'ignoranza non è mai meglio della conoscenza.” Enrico Fermi.
Offline
the procedure is the same if you use STAT_NON_LINE
I ever used code_aster for study.
Code_Aster ??? is the Novacula Occami of the FEM world.
“Non è bene cercare di fermare il progresso della conoscenza. L'ignoranza non è mai meglio della conoscenza.” Enrico Fermi.
Offline
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
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 ??? :-)
I ever used code_aster for study.
Code_Aster ??? is the Novacula Occami of the FEM world.
“Non è bene cercare di fermare il progresso della conoscenza. L'ignoranza non è mai meglio della conoscenza.” Enrico Fermi.
Offline
attachment
I ever used code_aster for study.
Code_Aster ??? is the Novacula Occami of the FEM world.
“Non è bene cercare di fermare il progresso della conoscenza. L'ignoranza non è mai meglio della conoscenza.” Enrico Fermi.
Offline
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"
Offline
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.
the procedure is the same drawing a tube (or a plate with hole)
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
I ever used code_aster for study.
Code_Aster ??? is the Novacula Occami of the FEM world.
“Non è bene cercare di fermare il progresso della conoscenza. L'ignoranza non è mai meglio della conoscenza.” Enrico Fermi.
Offline
attached here an example of capability
I ever used code_aster for study.
Code_Aster ??? is the Novacula Occami of the FEM world.
“Non è bene cercare di fermare il progresso della conoscenza. L'ignoranza non è mai meglio della conoscenza.” Enrico Fermi.
Offline