找回密码
 注册
Simdroid-非首页
查看: 632|回复: 9

(原创)自己做了一个热应力分析的例子,有温度加载过程,发出来吧

[复制链接]
发表于 2008-3-15 16:49:07 | 显示全部楼层 |阅读模式 来自 江苏南京
对一固支梁进行温度加载,然后进行热应力分析,载荷为三步,命令流如下,希望对你有些帮助。


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             !子步载荷线性插值,0ramped
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            !子步载荷线性插值,0ramped
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 编辑 ]
发表于 2008-10-19 19:09:10 | 显示全部楼层 来自 河南郑州
Simdroid开发平台
谢谢楼主
回复 不支持

使用道具 举报

发表于 2008-10-22 17:04:00 | 显示全部楼层 来自 重庆南岸区
不错,谢谢楼主。。。
回复 不支持

使用道具 举报

发表于 2008-10-22 20:32:45 | 显示全部楼层 来自 辽宁沈阳

(原创)自己做了一个热应力分析的例子,有温度加载过程,发出来吧

多谢热心的楼主了
回复 不支持

使用道具 举报

发表于 2009-11-29 15:34:22 | 显示全部楼层 来自 陕西西安
很不错的例子,楼主费心了.
回复 不支持

使用道具 举报

发表于 2009-12-4 10:18:21 | 显示全部楼层 来自 北京航空航天大学
正在学习中,对我有帮助,谢谢
回复 不支持

使用道具 举报

发表于 2009-12-9 14:20:58 | 显示全部楼层 来自 吉林长春
zhichiyixia very
回复 不支持

使用道具 举报

发表于 2009-12-10 11:35:41 | 显示全部楼层 来自 吉林长春
感谢楼主分享,正好搜寻有关热应力分析的命令流呢!!
回复 不支持

使用道具 举报

发表于 2009-12-15 16:16:58 | 显示全部楼层 来自 湖南长沙
多谢楼主,
回复 不支持

使用道具 举报

发表于 2009-12-18 10:14:28 | 显示全部楼层 来自 江苏淮安
非常好 谢谢了
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Simapps系列直播

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-9-20 21:25 , Processed in 0.059481 second(s), 12 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表