求助：瞬态热应力问题？
高手帮忙：瞬态热应力问题？/PREP7
ET,1,SOLID70
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,1,,43.27
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,2,,0.81
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,3,,1.687
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,7850
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,1,,450
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,1830
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,2,,860
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,3,,2390
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,3,,1186
*SET,r2,0.0889
*SET,t1,0.00919
*SET,r1,r2-t1
*SET,t2,0.03
*SET,r3,r2+t2
CYLIND,r2,r1,-5,5,0,90,
CYLIND,r3,r2,-5,5,0,90,
BLOCK,,10,,10,-5,5,
FLST,2,3,6,ORDE,2
FITEM,2,1
FITEM,2,-3
VPTN,P51X
VDELE, 4, , ,1
FLST,2,3,6,ORDE,3
FITEM,2,1
FITEM,2,5
FITEM,2,-6
VGLUE,P51X
NUMCMP,VOLU
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT, 1, , 1, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
CM,_Y,VOLU
VSEL, , , , 2
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT, 2, , 1, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
CM,_Y,VOLU
VSEL, , , , 3
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT, 3, , 1, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
SMRT,6
FLST,5,3,6,ORDE,2
FITEM,5,1
FITEM,5,-3
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VSWEEP,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
/UI,MESH,OFF
FINISH
/SOL
ALLSEL,ALL
ANTYPE,4
TRNOPT,FULL
LUMPM,0
TUNIF,50,
CSYS,1
ASEL,S,LOC,X,r1
FLST,2,1,5,ORDE,1
FITEM,2,4
/GO
SFA,P51X,1,CONV,10e3,350
ALLSEL,ALL
CSYS,0
OUTRES,BASIC,ALL,
TIME,3600*24*20
AUTOTS,1
NSUBST,80, , ,1
KBC,0
TSRES,ERASE
SOLVE
FINISH
/POST1
FINISH
/PREP7
ETCHG,TTS
MPTEMP,,,,,,,,
MPTEMP,1,50
MPTEMP,2,150
MPTEMP,3,250
MPTEMP,4,350
MPDATA,EX,1,,195.7e9
MPDATA,EX,1,,169.7e9
MPDATA,EX,1,,143.7e9
MPDATA,EX,1,,117.7e9
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPTEMP,,,,,,,,
MPTEMP,1,0
UIMP,1,REFT,,,
MPDATA,ALPX,1,,12.1e-6
MPTEMP,,,,,,,,
MPTEMP,1,0
UIMP,2,REFT,,,
MPDATA,ALPX,2,,10.3e-6
MPTEMP,,,,,,,,
MPTEMP,1,0
UIMP,3,REFT,,,
MPDATA,ALPX,3,,10.3e-6
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,15e9
MPDATA,PRXY,2,,0.25
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,3,,10e9
MPDATA,PRXY,3,,0.25
FINISH
/SOL
ALLSEL,ALL
CSYS,1
ASEL,S,LOC,X,r1
FLST,2,1,5,ORDE,1
FITEM,2,4
/GO
SFA,P51X,1,PRES,12e6
ALLSEL,ALL
CSYS,0
ASEL,S,LOC,Z,5
FLST,2,3,5,ORDE,3
FITEM,2,2
FITEM,2,24
FITEM,2,30
/GO
DA,P51X,UZ,
ALLSEL,ALL
ASEL,S,LOC,Z,-5
FLST,2,3,5,ORDE,3
FITEM,2,1
FITEM,2,23
FITEM,2,29
/GO
DA,P51X,UZ,
ALLSEL,ALL
ASEL,S,LOC,X,10
FLST,2,1,5,ORDE,1
FITEM,2,18
/GO
SFA,P51X,1,PRES,24.26e6
ALLSEL,ALL
ASEL,S,LOC,Y,10
FLST,2,1,5,ORDE,1
FITEM,2,16
/GO
SFA,P51X,1,PRES,17.423e6
ALLSEL,ALL
ASEL,S,LOC,Y,0
FLST,2,3,5,ORDE,3
FITEM,2,5
FITEM,2,26
FITEM,2,28
DA,P51X,SYMM
ALLSEL,ALL
ASEL,S,LOC,X,0
FLST,2,3,5,ORDE,3
FITEM,2,6
FITEM,2,25
FITEM,2,27
DA,P51X,SYMM
ALLSEL,ALL
TREF,50,
LDREAD,TEMP,,, , ,'ZHS','rth',' '
SOLVE
FINISH
/POST1
各位高手：模型中间为一钢管，管内流体温度350度，流动时间20天，现在要求每一天所对应的钢管温度、应力。（这里建的是一个四分之一模型）
用上面的命令流进行模拟时，存在下面的问题：
1、流体流动20天时（3600*24*20）的中间时刻（10天时）对应的管子的温度值，与直接流体通过10天时（3600*24*10）算出来的温度结果不相等？而实际情况下应该是相等的。
2、用上面的命令流计算出来的管子内壁的应力结果感觉也不合理。
哪位高手对瞬态热应力求解比较熟悉的请不吝赐教啊！！！
也请各位简要介绍一下模型瞬态温度场（包括不同的载荷步时）、应力场的计算步骤。
页:
[1]