我在做一个热--应力耦合的分析,现在得出的结果(见附件云图)很奇怪,大家帮我分析看看,为什么进行热应力耦合计算后,结构的变形方向和没考虑温度是完全相反,而且应力大小相差很大!!!
下面是我进行热应力耦合分析时的命令流文件,大家帮我看看是不是那里有问题:
*SET,a1,0.55
*SET,b1,0.37
*SET,r1,2.4
*SET,r2,1.1
*SET,dens,7.8e3
/PREP7
ET,1,SOLID87
MP,KXX,1,49.3
CYLIND,3.0,0.5,0,a1,0,360,
CYL4,0,2.627,0.13, , , ,a1
CSYS,1
FLST,3,1,6,ORDE,1
FITEM,3,2
VGEN,12,P51X, , , ,30, , ,0
FLST,3,12,6,ORDE,2
FITEM,3,2
FITEM,3,-13
VSBV, 1,P51X
CYLIND,r1,r2,0,b1,0,20,
FLST,3,1,6,ORDE,1
FITEM,3,1
VGEN, ,P51X, , , ,5, , , ,1
FLST,3,1,6,ORDE,1
FITEM,3,1
VGEN,12,P51X, , , ,30, , ,0
FLST,3,12,6,ORDE,2
FITEM,3,1
FITEM,3,-12
VSBV, 14,P51X
CYLIND,0.15,0,0,b1,0,360,
CSYS,0
FLST,3,1,6,ORDE,1
FITEM,3,1
VGEN, ,P51X, , , ,0.85, , , ,1
CSYS,1
FLST,3,1,6,ORDE,1
FITEM,3,1
VGEN, ,P51X, , , ,15, , , ,1
FLST,3,1,6,ORDE,1
FITEM,3,1
VGEN,12,P51X, , , ,30, , ,0
FLST,3,12,6,ORDE,2
FITEM,3,1
FITEM,3,-12
VSBV, 13,P51X
CYLIND,1.85,0,a1,a1+0.37,0,360,
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,14
VGLUE,P51X
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL
/VIEW,1,,1
/ANG,1
/REP,FAST
TYPE, 1
MAT, 1
REAL,
ESYS, 0
SECNUM,
ESIZE,0.2,0,
MSHAPE,1,3D
MSHKEY,0
FLST,5,2,6,ORDE,2
FITEM,5,1
FITEM,5,-2
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
FINISH
/SOL
FLST,2,1,5,ORDE,1
FITEM,2,131
/GO
SFA,P51X,1,CONV,7,27
ASEL,S,LOC,Z,0.910009,0.920001
APLOT
FLST,2,1,5,ORDE,1
FITEM,2,50
/GO
SFA,P51X,1,CONV,300,160
ALLSEL,ALL
/REPLOT
SAVE
ANTYPE,0
/STATUS,SOLU
SOLVE
FINISH
/POST1
/NOPR
/PMETH,OFF,0
KEYW,PR_SET,1
KEYW,PR_STRUC,1
KEYW,PR_THERM,0
KEYW,PR_FLUID,0
KEYW,PR_ELMAG,0
KEYW,MAGNOD,0
KEYW,MAGEDG,0
KEYW,MAGHFE,0
KEYW,MAGELC,0
KEYW,PR_MULTI,0
KEYW,PR_CFD,0
/GO
/COM,
/COM,Preferences for GUI filtering have been set to display:
/COM, Structural
FINISH
/PREP7
LSCLEAR,SOLID
ETCHG,TTS
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2.06e11
MPDATA,PRXY,1,,0.3
MPTEMP,,,,,,,,
MPTEMP,1,0
UIMP,1,REFT,,,27
MPDATA,ALPX,1,,1.62e-5
FINISH
/SOL
FLST,2,24,5,ORDE,24
FITEM,2,8
FITEM,2,-9
FITEM,2,12
FITEM,2,-13
FITEM,2,16
FITEM,2,-17
FITEM,2,20
FITEM,2,-21
FITEM,2,24
FITEM,2,-25
FITEM,2,28
FITEM,2,-29
FITEM,2,32
FITEM,2,-33
FITEM,2,36
FITEM,2,-37
FITEM,2,40
FITEM,2,-41
FITEM,2,44
FITEM,2,-45
FITEM,2,48
FITEM,2,-49
FITEM,2,52
FITEM,2,-53
/GO
DA,P51X,ALL,
LDREAD,TEMP,,, , ,'file','rth',' '
ASEL,S,LOC,Z,0.919999,0.920001
APLOT
FLST,2,1,5,ORDE,1
FITEM,2,50
/GO
SFA,P51X,1,PRES,-3.56e6
ALLSEL,ALL
/REPLOT |