Problem in computing thermal stresses
- Student
- Topic Author
- Offline
- New Member
-
Less
More
- Posts: 2
- Thank you received: 0
9 years 11 months ago - 9 years 11 months ago #8258
by Student
Problem in computing thermal stresses was created by Student
Hello Friends,
I am trying to do thermal analysis. For that I have divided my comm file in two sections. I do get thermal stresses in resu file but I am not able to get it on rmed file.
Can any of you please share your knowledge to solve this problem?
Thank you very much in advance.
I am trying to do thermal analysis. For that I have divided my comm file in two sections. I do get thermal stresses in resu file but I am not able to get it on rmed file.
Can any of you please share your knowledge to solve this problem?
Thank you very much in advance.
DEBUT();
mesh1=LIRE_MAILLAGE(UNITE= 20, FORMAT='MED',)
TH_MOD=AFFE_MODELE(MAILLAGE=mesh1,
AFFE=(_F(GROUP_MA=('Block1_vol',),
PHENOMENE='THERMIQUE',
MODELISATION='3D',),
_F(GROUP_MA=('Block1_contact','Block1_top','Block_bottom',),
PHENOMENE='THERMIQUE',
MODELISATION='3D',),),);
MAT1=DEFI_MATERIAU(ELAS=_F(E=900e3,
NU=0.2,
RHO=2100,
ALPHA=3.0E-6,),
THER=_F(LAMBDA=26,
RHO_CP=720000,
),);
MATAF=AFFE_MATERIAU(MAILLAGE=mesh1,
MODELE=TH_MOD,
AFFE=_F(GROUP_MA='Block1_vol',
MATER=MAT1,),
INFO=1,);
T950=DEFI_FONCTION(NOM_PARA='INST',
VALE=(
1728, 15,
3456, 30,
5184, 45,
6912, 60,
8640, 75,
10368, 90,
12096, 105,
13824, 120,
15552, 135,
17280, 150,
19008, 165,
20736, 180,
22464, 195,
24192, 210,
25920, 225,
27648, 240,
29376, 255,
31104, 270,
32832, 285,
34560, 300,
36288, 315,
38016, 320,
39744, 328,
41472, 336,
43200, 344,
44928, 352,
46656, 360,
48384, 368,
50112, 376,
51840, 384,
53568, 392,
55296, 400,
57024, 408,
58752, 416,
60480, 424,
62208, 432,
63936, 440,
65664, 448,
67392, 456,
69120, 464,
70848, 472,
72576, 480,
74304, 488,
76032, 496,
77760, 504,
79488, 512,
81216, 520,
82944, 528,
84672, 536,
86400, 544,
88128, 552,
89856, 560,
91584, 568,
93312, 576,
95040, 584,
96768, 592,
98496, 600,
100224, 608,
101952, 616,
103680, 624,
105408, 632,
107136, 640,
108864, 648,
110592, 656,
112320, 664,
114048, 672,
115776, 680,
117504, 688,
119232, 696,
120960, 704,
122688, 712,
124416, 720,
126144, 728,
127872, 736,
129600, 744,
131328, 752,
133056, 760,
134784, 768,
136512, 776,
138240, 784,
139968, 792,
141696, 800,
143424, 808,
145152, 816,
146880, 824,
148608, 832,
150336, 840,
152064, 848,
153792, 856,
155520, 864,
157248, 872,
158976, 880,
160704, 888,
162432, 896,
164160, 904,
165888, 912,
167616, 920,
169344, 928,
171072, 936,
172800, 944,
259200, 920,
432120, 950,
864000, 900,
8640000,950,
8812800,920,
12960000,950,
43200000,900,
60480000,950,
86400000,920,
172800000,950,
172972800,20,
),
PROL_DROITE='EXCLU',
PROL_GAUCHE='EXCLU', );
load1=AFFE_CHAR_THER_F( MODELE= TH_MOD,
TEMP_IMPO=_F(GROUP_MA='Block1_top',TEMP=T950,),
);
list= DEFI_LIST_REEL( DEBUT=0., INTERVALLE=(_F( JUSQU_A=172800, NOMBRE=100,),
_F( JUSQU_A=864000, NOMBRE = 10,),
_F( JUSQU_A=8640000, NOMBRE = 10,),
_F( JUSQU_A=8812800, NOMBRE = 10,),
_F( JUSQU_A=12960000, NOMBRE = 10,),
_F( JUSQU_A=43200000, NOMBRE = 10,),
_F( JUSQU_A=60480000, NOMBRE = 10,),
_F( JUSQU_A=86400000, NOMBRE = 10,),
_F( JUSQU_A=172800000, NOMBRE = 10,),
_F( JUSQU_A=172972800, NOMBRE = 10,),
),)
Result1=THER_LINEAIRE(MODELE=TH_MOD,
CHAM_MATER=MATAF,
EXCIT=_F(CHARGE=load1,),
INCREMENT=_F(LIST_INST=list,),
ETAT_INIT =_F(VALE=20.0,),
);
IMPR_RESU(
FORMAT='MED',
UNITE=80,
RESU=_F(MAILLAGE=mesh1,
RESULTAT=Result1,
NOM_CHAM='TEMP',),);
#-----------------------------------------------------------------------------------------------------------------------
mesh2=LIRE_MAILLAGE(UNITE= 21, FORMAT='MED',)
ME_MOD=AFFE_MODELE(MAILLAGE=mesh2,
AFFE=(_F(GROUP_MA=('Block1_vol',),
PHENOMENE='MECANIQUE',
MODELISATION='3D',),
_F(GROUP_MA=('Block1_contact','Block1_top','Block_bottom',),
PHENOMENE='MECANIQUE',
MODELISATION='3D',),),);
Proj=PROJ_CHAMP(#METHODE='ELEM',
RESULTAT=Result1,
MODELE_1=TH_MOD,
MODELE_2=ME_MOD,
VIS_A_VIS=_F(TOUT_1='OUI',
TOUT_2='OUI',
CAS_FIGURE='3D',),);
MATAF1=AFFE_MATERIAU(MAILLAGE=mesh2,
MODELE=ME_MOD,
AFFE=_F(GROUP_MA='Block1_vol',
MATER=MAT1,),
AFFE_VARC =_F( TOUT='OUI',
NOM_VARC = 'TEMP',
EVOL = Proj,
VALE_REF = 20
),
INFO=1,);
BC1=AFFE_CHAR_MECA(MODELE =ME_MOD,
DDL_IMPO =_F(GROUP_MA ='Block_bottom',
DX=0,
DY=0,
DZ=0,),);
BC2=AFFE_CHAR_MECA(MODELE = ME_MOD,
FORCE_FACE =_F(GROUP_MA ='Block1_top',
FZ = 12943.72,),);
Result2=MECA_STATIQUE(MODELE=ME_MOD,
CHAM_MATER=MATAF1,
EXCIT=(_F(CHARGE=BC1,),
_F(CHARGE=BC2,),),);
Result2=CALC_CHAMP(reuse =Result2,
RESULTAT=Result2,
CONTRAINTE=('SIGM_ELNO',),
#FORCE='REAC_NODA',
);
Result2=CALC_CHAMP(reuse =Result2,
RESULTAT=Result2,
DEFORMATION=('EPSI_ELGA',),);
Result2=CALC_CHAMP(reuse =Result2,
RESULTAT=Result2,
CRITERES=('SIEQ_ELGA',),);
table01=CREA_TABLE(RESU=_F( RESULTAT= Result2,
GROUP_MA='Block1_vol',
NOM_CHAM='EPSI_ELGA',
NOM_CMP=('EPXX','EPYY','EPZZ',),
#INST=1.0,
),);
IMPR_TABLE (TABLE=table01);
table02=CREA_TABLE(RESU=_F( RESULTAT= Result2,
GROUP_MA='Block1_vol',
NOM_CHAM='SIGM_ELGA',
NOM_CMP=('SIXX','SIYY','SIZZ',),
#INST=1.0,
),);
IMPR_TABLE (TABLE=table02);
IMPR_RESU(FORMAT='MED',
UNITE=81,
RESU=_F(
MAILLAGE=mesh2,
RESULTAT=Result2,
NOM_CHAM='SIGM_ELNO',
#NOM_CMP=('VMIS',),
),);
IMPR_RESU(FORMAT='MED',
UNITE=82,
RESU=_F(
MAILLAGE=mesh2,
RESULTAT=Result2,
#TOUT_CHAM='OUI',
NOM_CHAM='EPSI_ELGA',
),);
FIN();
Last edit: 9 years 11 months ago by Student.
- Student
- Topic Author
- Offline
- New Member
-
Less
More
- Posts: 2
- Thank you received: 0
9 years 11 months ago - 9 years 11 months ago #8259
by Student
Replied by Student on topic Re: Problem in computing thermal stresses
This problem has been solved by following procedure,
Filters -> Integration points -> (e.g) ELNO Mesh
Filters -> Integration points -> (e.g) ELNO Mesh
Last edit: 9 years 11 months ago by Student.
Moderators: catux
Time to create page: 0.130 seconds