钢框架受热分析

ccppss

此例子是选自《ansys高级工程有限元分析范例精选》，感觉很好，贴出来
!以下是防火间受热分析
/prep7
/title,part 1:heat transfer analysis
et,1,solid70
w_col=0.2039          !定义建模参数
h_col=0.2062
tf_col=0.0125
tw_col=0.0083
b_col=(w_col-tw_col)/2
d_col=h_col-2*tf_col
w_beam=0.165
h_beam=0.3038
tf_beam=0.0102
tw_beam=0.0061
b_beam=(w_beam-tw_beam)/2
d_beam=h_beam-2*tf_beam
dis_hor=5.5
dis_ver=3.0
pp=(w_col-w_beam)/2
mptemp,,20,800,900,1000
mpdata,kxx,1,,53.334,27.36,27.36,27.36
mptemp
mptemp,,20,100,180,260,380
mpdata,c,1,,439.8,487.62,522.33,550.75,596.52
mptemp,,500,600,640,720,735
mpdata,c,1,,666.5,759.92,798.67,1388,5000
mptemp,,750,830,900,1000
mpdata,c,1,,1483,725,650,650
mp,dens,1,7850
n,1,-h_col/2,,-w_col/2
n,2,-h_col/2,,-w_col/2+pp
ngen,4,1,2,,,,,(b_col-pp)/3
n,6,-h_col/2,,tw_col/2
ngen,4,1,6,,,,,(b_col-pp)/3
n,10,-h_col/2,,w_col/2
ngen,2,10,1,10,1,tf_col
ngen,10,6,15,16,1,d_col/10
ngen,2,60,11,20,1,d_col
ngen,2,10,71,80,1,tf_col
ngen,2,100,all,,,,(dis_ver-h_beam/2)/60
e,1,2,12,11,101,102,112,111
egen,9,1,1
e,15,16,22,21,115,116,122,121
egen,10,6,10
e,71,72,82,81,171,172,182,181
egen,9,1,20
egen,60,100,1,28,1,,,,,,0,(dis_ver-h_beam/2)/60
nsel,s,node,,6001,6090,1
ngen,2,100,all,,,,tf_beam
egen,2,100,1653,1680,1
nsel,all
nsel,s,node,,6101,6190
ngen,11,100,all,,,,d_beam/10
egen,11,100,1681,1708,1
nsel,all
nsel,s,node,,7101,7190
ngen,2,100,all,,,,tf_beam
egen,2,100,1961,1988,1
nsel,all
nsel,s,node,,7201,7290
ngen,7,100,all,,,,h_beam/6
egen,7,100,1989,2016,1
nsel,all
ngen,2,10000,all,,,dis_hor
egen,2,10000,1,2184,1
n,20001,h_col/2,dis_ver-h_beam/2,-w_beam/2
ngen,4,1,20001,,,,,b_beam/3
n,20005,h_col/2,dis_ver-h_beam/2,tw_beam/2
ngen,4,1,20005,,,,,b_beam/3
ngen,2,10,20001,20008,1,,tf_beam
ngen,10,6,20014,20015,1,,d_beam/10
ngen,2,60,20011,20018,1,,d_beam
ngen,2,10,20071,20078,1,,tf_beam
ngen,2,100,20001,20090,,(dis_hor-h_col)/100
e,20001,20002,20012,20011,20101,20102,20112,20111
egen,7,1,4369
e,20014,20015,20021,20020,20114,20115,20121,20120
egen,10,6,4376
e,20071,20072,20082,20081,20171,20172,20182,20181
egen,7,1,4386
egen,100,100,4369,4392,1,,,,,,(dis_hor-h_col)/100
cpintf,all,0.002
finish
/solu
antype,trans
tunif,20
nsel,s,node,,71,6071,100
*do,i,80,90,1
nsel,a,node,,i,6000+i,100
*enddo
cm,htbound1,node
nsel,all
nsel,s,node,,10020,16020,100
*do,i,10001,10011,1
nsel,a,node,,i,6000+i,100
*enddo
cm,htbound2,node
nsel,all
nsel,s,node,,20011,30011,100
nsel,a,node,,20018,30018,100
*do,i,20001,20008,1
nsel,a,node,,i,10000+i,100
*enddo
cm,htbound3,node
nsel,all
nsel,s,node,,20011,30011,100
nsel,a,node,,20012,30012,100
nsel,a,node,,20013,30013,100
nsel,a,node,,20071,30071,100
nsel,a,node,,30072,30072,100
nsel,a,node,,20073,30073,100
*do,i,20014,20074,6
nsel,a,node,,i,10000+i,100
*enddo
cm,htbound4,node
nsel,all
nsel,s,node,,20016,30016,100
nsel,a,node,,20017,30017,100
nsel,a,node,,20018,30018,100
nsel,a,node,,20076,30076,100
nsel,a,node,,20077,30077,100
nsel,a,node,,20078,30078,100
*do,i,20015,20075,6
nsel,a,node,,i,10000+i,100
*enddo
cm,htbound5,node
nsel,all
*do,tm,60,180,60
time,tm
deltim,20
autots,on
temp=20+345*log10(8*tm/60+1)
sf,htbound1,conv,25,temp    !施加对流作用
sf,htbound2,conv,25,temp
sf,htbound3,conv,25,temp
sf,htbound4,conv,25,temp
sf,htbound5,conv,25,temp
sf,htbound1,rdsf,0.9,1     !定义辐射场
sf,htbound2,rdsf,0.9,2
sf,htbound3,rdsf,0.9,3
sf,htbound4,rdsf,0.9,4
sf,htbound5,rdsf,0.9,5
stef,5.6696e-8
toffst,273
spctemp,1,temp
spctemp,2,temp
spctemp,3,temp
spctemp,4,temp
spctemp,5,temp
solve
*enddo
finish
/post1
plnsol,temp,,0

!以下是受力分析，热分析结果作为载荷读入
/prep7
/title,part2:structrural analysis
et,1,solid45
et,2,beam188
fy=275e+6
exx=2.1e+11
mptemp
mptemp,,20,100,200,300,400
mpdata,ex,1,,exx,exx,0.9*exx,0.8*exx,0.7*exx
mptemp,,500,600,700,800,900
mpdata,ex,1,,0.6*exx,0.31*exx,0.13*exx,0.09*exx,0.0675*exx
mp,nuxy,1,0.3
mp,alpx,1,1.4e-5
tb,miso,1,10,3
tbtemp,20
tbpt,,fy/exx,fy
tbpt,,0.02,fy
tbpt,,0.15,fy
tbtemp,100
tbpt,,fy/exx,fy
tbpt,,0.02,fy
tbpt,,0.15,fy
tbtemp,200
tbpt,,0.807*fy/(0.9*exx),0.807*fy
tbpt,,0.02,fy
tbpt,,0.15,fy
tbtemp,300
tbpt,,0.613*fy/(0.8*exx),0.613*fy
tbpt,,0.02,fy
tbpt,,0.15,fy
tbtemp,400
tbpt,,0.420*fy/(0.7*exx),0.420*fy
tbpt,,0.02,fy
tbpt,,0.15,fy
tbtemp,500
tbpt,,0.360*fy/(0.6*exx),0.360*fy
tbpt,,0.02,0.780*fy
tbpt,,0.15,0.780*fy
tbtemp,600
tbpt,,0.180*fy/(0.310*exx),0.180*fy
tbpt,,0.02,0.470*fy
tbpt,,0.15,0.470*fy
tbtemp,700
tbpt,,0.075*fy/(0.130*exx),0.075*fy
tbpt,,0.02,0.230*fy
tbpt,,0.15,0.230*fy
tbtemp,800
tbpt,,0.05*fy/(0.09*exx),0.05*fy
tbpt,,0.02,0.110*fy
tbpt,,0.15,0.110*fy
tbtemp,900
tbpt,,0.0375*fy/(0.0675*exx),0.0375*fy
tbpt,,0.02,0.060*fy
tbpt,,0.15,0.060*fy
sectype,1,beam,i,column
secdata,w_col,w_col,h_col,tf_col,tf_col,tw_col
sectype,2,beam,i,beam
secdata,w_beam,w_beam,h_beam,tf_beam,tf_beam,tw_beam
k,1,,dis_ver+h_beam*1.5
k,2,,2*dis_ver
k,3,,3*dis_ver
k,4,dis_hor,dis_ver+h_beam*1.5
k,5,dis_hor,dis_ver*2
k,6,dis_hor,dis_ver*3
k,7,dis_hor+h_col/2,dis_ver
k,8,2*dis_hor
k,9,2*dis_hor,dis_ver
k,10,2*dis_hor,2*dis_ver
k,11,2*dis_hor,3*dis_ver
k,12,3*dis_hor
k,13,3*dis_hor,dis_ver
k,14,3*dis_hor,2*dis_ver
k,15,3*dis_hor,3*dis_ver
k,100,-3,3
k,200,5,20
l,1,2
l,2,3
l,4,5
l,5,6
l,8,9
l,9,10
l,10,11
l,12,13
l,13,14
l,14,15
l,2,5
l,3,6
l,7,9
l,5,10
l,6,11
l,9,13
l,10,14
l,11,15
lsel,s,line,,1,10,1
latt,1,,2,,100,,1
lsel,all
lsel,s,line,,11,18,1
latt,1,,2,,200,,2
lsel,all
lesize,all,0.3
lmesh,all
cpintf,all,0.002
n1=node(0,dis_ver+h_beam*1.5,0)
num=0
*do,k,7801,7820,1
num=num+1
dx=nx(k)
ce,num,0,k,uy,1,n1,uy,-1,n1,rotz,-dx
*enddo
*do,k,7871,7890,1
num=num+1
dx=nx(k)
ce,num,0,k,uy,1,n1,uy,-1,n1,rotz,-dx
*enddo
*do,k,7821,7869,6
num=num+1
dx=nx(k)
ce,num,0,k,uy,1,n1,uy,-1,n1,rotz,-dx
num=num+1
dx=nx(k+1)
ce,num,0,k+1,uy,1,n1,uy,-1,n1,rotz,-dx
*enddo
nsel,s,node,,n1
nsel,a,node,,7821,7869,6
nsel,a,node,,7822,7870,6
cp,next,ux,all
nsel,all
n4=node(dis_hor,dis_ver+h_beam*1.5,0)
*do,k,17801,17820,1
num=num+1
dx=nx(k)-dis_hor
ce,num,0,k,uy,1,n4,uy,-1,n4,rotz,-dx
*enddo
*do,k,17871,17890,1
num=num+1
dx=nx(k)-dis_hor
ce,num,0,k,uy,1,n4,uy,-1,n4,rotz,-dx
*enddo
*do,k,17821,17869,6
num=num+1
dx=nx(k)-dis_hor
ce,num,0,k,uy,1,n4,uy,-1,nr,rotz,-dx
num=num+1
dx=nx(k+1)-dis_hor
ce,num,0,k+1,uy,1,n4,uy,-1,nr,rotz,-dx
*enddo
nsel,s,node,,n4
nsel,a,node,,17821,17869,6
nsel,a,node,,17822,17870,6
cp,next,ux,all
nsel,all

n7=node(dis_hor+h_col/2,dis_ver,0)
*do,i,16000,16100,100
*do,j,81,90,1
num=num+1
dy=ny(i+j)-dis_ver
ce,num,0,i+j,ux,1,n7,ux,-1,n7,rotz,dy
*enddo
*enddo

*do,i,17100,17200,100
*do,j,81,90,1
num=num+1
dy=ny(i+j)-dis_ver
ce,num,0,i+j,ux,1,n7,ux,-1,n7,rotz,dy
*enddo
*enddo

nsel,s,node,,n7
nsel,a,node,,16285,17085,100
nsel,a,node,,16285,17086,100
cp,next,uy,all
nsel,all
fini

/solu
antype,0
tref,20
nsel,s,loc,y,0
d,all,all
nsel,all

dk,13,ux
dk,14,ux
dk,15,ux
ksel,u,kp,,100,200,100
dk,all,uz
dk,all,rotx
dk,all,roty
ksel,all

fk,3,fy,-75500
fk,6,fy,-151000
fk,11,fy,-151000
fk,15,fy,-75500

lsel,s,line,,11,18,1
sfl,all,pres,-25400
lsel,all

nsel,s,node,,20084,30084,100
nsel,a,node,,20085,30085,100
sf,all,pres,25400/tw_beam
nsel,all

time,1
deltim,0.2,1e-3,0.5
solve

*do,tm,60,180,60
time,tm
ldread,temp,,,tm,,,rth
deltim,20,1,20
solve
*enddo
finish

/post1
plnsol,u,y

yywodelanqiu

谢谢分享，下来看看