可能你没理解重启动分析中 TIME 的含义,下面是略微修改过的命令流:
/PREP7
ET,1,PLANE13 !定义单元类型
MP,DENS,1,1400 !定义型材的密度
MPTEMP,,-3,27,47,67 !定义随温度变化的型材的导热性
MPDATA,KXX,1,,0.167,0.170,0.171,0.172
MPTEMP !清除温度场
MPTEMP,,2,23,87,107,147 !定义随温度变化的型材的比热
MPDATA,C,1,,865,960,1475,1569,1741
MPTEMP !清除原来的温度场
MPTEMP,,20,60,80,100,140,180
MPDATA,EX,1,,2.18E+9,1.30E+9,0.86E+8,7.2E6 ,4.22E6 ,8.31E5 !定义随温度变化的弹性模量
MP,PRXY,1,0.3 !定义泊松比
MPTEMP !清除原来的温度场
MP,ALPX,1,1e-4 !定义热膨胀系数
!这里定义了型材的各项参数,采用耦合单元PLANE13
K,1
K,2,0.002
K,3,0.058
K,4,0.06
K,5,,0.002
K,6,0.002,0.002
K,7,0.058,0.002
K,8,0.06,0.002
K,9,,0.058
K,10,0.002,0.058
K,11,0.058,0.058
K,12,0.06,0.058
K,13,,0.06
K,14,0.002,0.06
K,15,0.058,0.06
K,16,0.06,0.06
a,1,2,6,5
a,2,3,7,6
a,3,4,8,7
a,5,6,10,9
a,7,8,12,11
a,9,10,14,13
a,10,11,15,14
a,11,12,16,15
esize,0.001
MSHKEY,1
amesh,all
finish
/solu
antype,trans,new
trnopt,full !full方法
cnvtol,f,,0.1,,0.1 !设定以力作为收敛标准,不设定时在错误文件中会出现相关的错误信息。
tunif,190
tref,190
lsel,s,,,1
lsel,a,,,5
lsel,a,,,8
lsel,a,,,9
lsel,a,,,14
lsel,a,,,23
lsel,a,,,24
lsel,a,,,22
lsel,a,,,18
lsel,a,,,19
lsel,a,,,13
lsel,a,,,4
sfl,all,conv,2000,,15
allsel,all
lsel,s,,,1
lsel,a,,,5
lsel,a,,,8
lsel,a,,,24
lsel,a,,,22
lsel,a,,,18
dl,all,,uy,0
allsel,all
lsel,s,,,9
lsel,a,,,14
lsel,a,,,23
lsel,a,,,19
lsel,a,,,13
lsel,a,,,4
dl,all,,ux,0
allsel,all
timint,on
time,12
deltim,0.2,0.02,0.2,1
kbc,0
solve
finish
!为接着计算的程序
/solu
antype,trans,rest
trnopt,full
cnvtol,f,,0.1,,0.1
tref,190
lsel,s,,,1
lsel,a,,,5
lsel,a,,,8
lsel,a,,,9
lsel,a,,,14
lsel,a,,,23
lsel,a,,,24
lsel,a,,,22
lsel,a,,,18
lsel,a,,,19
lsel,a,,,13
lsel,a,,,4
sfldele,all,conv
allsel,all
lsel,s,,,1
lsel,a,,,5
lsel,a,,,8
lsel,a,,,24
lsel,a,,,22
lsel,a,,,18
dldele,all,uy
allsel,all
lsel,s,,,9
lsel,a,,,14
lsel,a,,,23
lsel,a,,,19
lsel,a,,,13
lsel,a,,,4
dldele,all,ux
allsel,all
d,1,ux
d,1,uy
d,178,uy
lsel,s,,,1
lsel,a,,,5
lsel,a,,,8
lsel,a,,,9
lsel,a,,,14
lsel,a,,,23
lsel,a,,,24
lsel,a,,,22
lsel,a,,,18
lsel,a,,,19
lsel,a,,,13
lsel,a,,,4
sfl,all,conv,3000,,15
allsel,all
time,12.2
deltim,0.05,0.005,0.05,1
kbc,0
solve
finish
这里面对模型的修改只是为了分网的需要,时间的定义是主要修改之处。 |