manbuzhe 发表于 2009-3-14 17:09

?热力耦合分析竟然不破坏

做钢管混凝土柱的热力耦合分析时,温度场分析正常,加上力耦合时总是烧不坏怎么回事啊?一直算到3小时还不破坏(没有保护层)。

请各位大侠赐教!谢过!



/solu
nsel,s,loc,y,0   !柱脚固端约束
d,all,uz,0
d,all,ux,0
d,all,uy,0
allsel,all

asel,s,loc,x,0      !定义对称边
da,all,symm
allsel,all

nsel,s,loc,x,0         !对称边约束,一半模型
d,all,ux,0
allsel,all

nsel,s,loc,y,L      !柱顶约束,并加力
d,all,uz,0
d,all,ux,0
f,all,fy,-5021.7
allsel,all
nsel,s,loc,y,L
cp,next,uy,all
allsel,all   

!施加荷载与边界条件并求解------------------------------------------------------------

ANTYPE,0            !静力分析
eqslv,sparse
tref,20
outres,all,50         !50步保存一次数据



!准备求解--------------------------------------------------------------------------
!常温下的反应对应于时间为0,在此,设置一个非常小的时间
NLGEOM,ON                !设置大变形效应
TIME,1                   !第一步常温下的反应分析,时间为1

sstif,off                !关闭应力刚化
!nsubst,50            !荷载步设置为50

cnvtol,u,,0.05             !定义收敛条件
neqit,50                   !每一子步中方程的迭代次数限值
cutcontrol,noiterpred,1
CUTCONTROL,PLSLIMIT,0.9!一次迭代内的最大允许塑性应变
ncnv,2                               !如果不收敛时结束而不退出
SOLVE                           !求解常温下的反应

autots,on
*Do,tm,30,10800,30   !设置时间从60到7200,步长60,即每1分钟求解一次
TIME,tm                  !当前时间为tm
sstif,off                         !关闭应力刚化
LDREAD,TEMP,,,tm,1,final,RTH    !读入时间tm时的温度分布
nsubst,500                        !荷载步设置为500
cnvtol,u,0.024,0.01         !位移控制收敛,柱高为2.4m
neqit,50                              !每一子步中方程的迭代次数限值
cutcontrol,noiterpred,1
CUTCONTROL,PLSLIMIT,0.9    !一次迭代内的最大允许塑性应变
ncnv,2                                    !如果不收敛时结束而不退出
SOLVE                  

*ENDDO
页: [1]
查看完整版本: ?热力耦合分析竟然不破坏