- 积分
- 0
- 注册时间
- 2007-10-17
- 仿真币
-
- 最后登录
- 1970-1-1
|
对一固支梁进行温度加载,然后进行热应力分析,载荷为三步,命令流如下,希望对你有些帮助。
finish
/clear,start
/FILENAME,heating
/prep7
et,1,plane55,,,1 !轴对称
mp,kxx,1,170
mp,c,1,691
mp,dens,1,2330
mat,1
*dim,wendu,arrar,13
wendu(1)=30
wendu(2)=150
wendu(3)=150
wendu(4)=150
wendu(5)=150
wendu(6)=150
wendu(7)=150
wendu(8)=150
wendu(9)=150
wendu(10)=150
wendu(11)=150
wendu(12)=150
wendu(13)=30
beamlegh=3e-4
beamdept=2e-6
rectng,0,beamlegh,0,beamdept
/replot
asel,s,area,,1
aatt,1
asel,all
meshsize=8e-7
lesize,all,meshsize
mshkey,1
amesh,all
/solu
Antype,trans
timint,on
!开始瞬态分析
trnopt,full !设置计算方法
lumpm,0 !使用默认单元质量矩阵
esel,s,mat,,1
nsle,s
tunif, 30 !初始温度为30度
*do,minute,2,24,2
*if,minute,le,2,then
nsel,s,loc,y,0
nsel,a,loc,y,beamdept
d,all,temp,wendu(1+minute/2)
alls
time,minute
autots,0 !不用自动载荷步功能
deltim,0.1 !让温度很接近于30度,一直没能做到初始温度为30度,或许哪位能改进一下
kbc,0 !子步载荷线性插值,0为ramped
outres,all,all
solve
*elseif,minute,gt,2,then
*if,minute,le,22,then
nsel,s,loc,y,0
nsel,a,loc,y,beamdept
d,all,temp,wendu(1+minute/2)
alls
time,minute
autots,0
deltim,1
kbc,1
outres,all,all
solve
*elseif,minute,gt,22,then
nsel,s,loc,y,0
nsel,a,loc,y,beamdept
d,all,temp,wendu(1+minute/2)
alls
time,minute
autots,0 !不用自动载荷步功能
deltim,1 !分成了两个子步
kbc,0 !子步载荷线性插值,0为ramped
outres,all,all
solve
*endif
*endif
*enddo
esel,s,mat,,1
nsle,s
!ddelt,all,temp
Tintpr,,,,1.0
Finish
/PREP7
etchg,tts
! plane55变plane42
mp,ex,1,1.71e11
mp,nuxy,1,0.28
mp,alpx,1,2.6e-6
mp,dens,1,2330
/SOLU
nropt,full !使用全Newton-Raphson选项来进行静态或全瞬态分析
tref,20 !零应变温度
nsel,s,loc,x,0 !选择y轴
d,all,all,0 !对上面两轴进行固支
alls
*do,minute,2, 24,2
asel,s,mat,,1
ldread,temp,,,minute,,'heating','rth',''
time,minute
autots,0
kbc,0
outres,all,all
solve
*enddo
Finish
[ 本帖最后由 badx2004 于 2008-3-15 18:12 编辑 ] |
|