【讨论】请各位帮我看看这个热应力分析的问题在哪？
命令流：/FILNAME,test
/title,an
/prep7
!指定单元类型
et,1,87,0
et,2,90,0
!指定单元属性
! 1号材料是bean Q235BH
! 2号材料是cloumn Q345BH-16Mn
! 3号材料是weldQ235BH
mptemp,1,20,200,500,800,1500
mpdata,ex,1,1,2.064e11,1.997e11,1.348e11,8.106e10,1.0e8
mpdata,ex,2,1,2.1e11,1.95e11,1.55e11,1.0e11,5e7
mp,alpx,1,11e-6
mp,alpx,2,12e-6
mp,kxx,1,68
mp,kxx,2,46
mp,c,1,480
mp,c,2,450
mp,dens,1,7840
mp,dens,2,7815
mp,prxy,1,0.285
mp,prxy,2,0.29
TB,BKIN,1,5
TBTEMP,20,1
TBDATA,1,2.45E8,2.064e10
TBTEMP,200,2
TBDATA,1,1.6E8,1.997e10
TBTEMP,600,3
TBDATA,1,0.7E8,1.348e10
TBTEMP,800,4
TBDATA,1,5E7,8.106e9
TBTEMP,1500,5
TBDATA,1,1E6,1.0e7
TB,BKIN,2,5
TBTEMP,20,1
TBDATA,1,3.9E8,2.1e10
TBTEMP,200,2
TBDATA,1,2.58E8,1.95e10
TBTEMP,600,3
TBDATA,1,1.37E8,1.55e10
TBTEMP,800,4
TBDATA,1,7.5E7,1.0e10
TBTEMP,1500,5
TBDATA,1,1E6,5e6
BLOCK,0,130e-3,0,4e-3,0,10e-3
wprota,,,90
wpoffs,,,70e-3
vsbw,1
wpoffs,,,5e-3
vsbw,3
wpstyl,defa
WPSTYLE,,,,,,,,0
mat,2
type,2
esize,1e-3
MSHAPE,1,3D
MSHKEY,0
vmesh,4
mat,1
type,1
esize,3e-3
MSHAPE,1,3D
MSHKEY,0
vmesh,1
vmesh,2
/SOLU
ANTYPE,TRANS
TREF,24
TIMINT,on
tintp,0.005,,,1,0.5,0.2
VSEL,S,,,4
NSLV,S,1
D,ALL,TEMP,1450
NSEL,INVE
D,ALL,TEMP,30
TIME,5
KBC,1
NSUBST,2,2,2
OUTRES,all,all
ALLSEL,all
SOLVE
OUTRES,all,all
DDELE,ALL,TEMP
TIME,40 !phase change
DELTIME,0.2,0.1,2
AUTOTS,ON
SOLVE
finish
/prep7
etchg,tts
finish
/SOLu
ANTYPE,0
nsel,s,loc,x,0 ,0
nsel,a,loc,x,130e-3 ,130e-3
d,all,ux,0
LDREAD,TEMP,,,, ,'test','rth',' '
time,1
nlgeom,on
nsubst,15
autots,on
pred,on
nropt,full,,on
TREF,24
allsel,all
solve
热应力的最大值应该在中间，可这个刚好相反，居然在两边。
请各位帮我看看出了什么问题？？
页:
[1]