/title,part 1:heat transfer ananlysis ET,1,SOLID70 !定义单元类型 !-----------------------------------------
!设置参数!单位制为：N；M；摄氏度；s 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
!将已生成的节点沿y偏移(dis_ver-h_beam/2)/60拷贝一层 !节点号加100
ngen,2,100,all,,,,(dis_ver-h_beam/2)/60 !从已生成的两层节点生成第一层单元
e,1,2,12,11,101,102,112,111 !定义翼缘第一个单元
egen,9,1,1 !沿翼缘长度拷贝9个单元
e,15,16,22,21,115,116,122,121 !定义腹板第一个单元，编号为10 egen,10,6,10 !拷贝单元10形成腹板
e,71,72,82,81,171,172,182,181 !定义另一翼缘一个单元，编号为20 egen,9,1,20 !拷贝单元20形成另一腹板 !将第一层单元沿y拷贝60层，自动生成所有节点，编号每层加100 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 !并生成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
!将实体模型的柱向上延伸h_beam的高度，避免梁单元和实体单元在梁柱节点处切换 nsel,s,node,,7201,7290 !生成6层单元，每层高h_beam/6 ngen,7,100,all,,,,h_beam/6 egen,7,100,1989,2016,1 nsel,all
!第一根实体模型柱完成
!共计：节点79层，每层编号1-90，逐层加100，顶层编号7801-7890 !单元78层，自动编号。每层28个，共28*78=2184个 !---------------------------------------- !拷贝第一根柱生成第二根柱
ngen,2,10000,all,,,dis_hor !拷贝所有节点，节点号加10000 egen,2,10000,1,2184,1 !拷贝所有单元
!----------------------------------------
!生成梁
!梁被加在两根柱之间，实际长度dis_hor-h_col
!生成左侧形成一个梁截面的所有节点，梁的节点编号从20001开始 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
!沿x偏移(dis_hor-h_col)/100拷贝一层节点，节点编号加100 ngen,2,100,20001,20090,,(dis_hor-h_col)/100
!生成梁的第一层截面单元
!两根柱单元总数为4368，故梁单元编号从4369开始
e,20001,20002,20012,20011,20101,20102,20112,20111 !定义单元4369 egen,7,1,4369 !拷贝单元4369形成一个翼缘
e,20014,20015,20021,20020,20114,20115,20121,20120 !定义单元4376 egen,10,6,4376 !拷贝单元4376形成腹板
e,20071,20072,20082,20081,20171,20172,20182,20181 !定义单元4386 egen,7,1,4386 !拷贝单元4386形成另一翼缘 !沿x拷贝100层生成整根梁，每层高度(dis_hor-h_col)/100 egen,100,100,4369,4392,1,,,,,,(dis_hor-h_col)/100
!梁的实体模型完成。
!总计：梁的节点为101层，每层编号1-88，从20001开始，逐层加100 !左端截面的节点为20001-20088，右端截面的节点为30001-30088 !每层单元数为24个，总计24*100=2400个，单元编号为4369-6768。 !----------------------------------------
!建立梁和柱连接处的耦合关系。用cpintf自动耦合所有节点坐标重合的节点 !梁翼缘的节点和柱的侧面完全重合，可以自动耦合
!梁腹板节点距离柱侧面相应节点距离为(tw_col-tw_beam)/2=0.0011 !因此，设置耦合误差为0.002时，也能自动耦合 !---------------------------------------- cpintf,all,0.002
!---------------------------------------- finish
!---------------------------------------- !定义边界条件，并求解
!---------------------------------------- /solu
antype,trans !定义分析类型 tunif,20 !定义初始温度 !定义受火边界
!选择第一根柱右侧翼缘的节点，定义为htbound1
nsel,s,node,,71,6071,100 !选择节点71到6071，间距100
*do,i,80,90,1 !定义参数i从80到90间距1 nsel,a,node,,i,6000+i,100 !另选择节点i到6000+i，间距100 *enddo !结束do循环
cm,htbound1,node !定义以上节点为htbound1 nsel,all
!选择第2根柱左侧翼缘的节点，定义为htbound2 nsel,s,node,,10020,16020,100 *do,i,10001,10011,1 nsel,a,node,,i,6000+i,100 *enddo
cm,htbound2,node nsel,all
!选择梁下翼缘表面，定义为htbound3 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
!选择梁的一侧表面，定义为htbound4 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,,20072,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
!选择梁的另一侧面，定义为htbound5 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 !定义时间参数tm从60到600秒 time,tm !当前时间为tm deltim,20 !定义初始时间步长 autots,on !打开自动步长控制 !----------------------------------------
temp=20+345*log10(8*tm/60+1) !计算环境空气温度 sf,htbound1,conv,25,temp !对边界1施加对流作用 sf,htbound2,conv,25,temp !对边界2施加对流作用 sf,htbound3,conv,25,temp !对边界3施加对流作用 sf,htbound4,conv,25,temp !对边界4施加对流作用 sf,htbound5,conv,25,temp !对边界5施加对流作用
sf,htbound1,rdsf,0.9,1 !定义边界1为第1个热辐射场 sf,htbound1,rdsf,0.9,2 !定义边界2为第2个热辐射场 sf,htbound1,rdsf,0.9,3 !定义边界3为第3个热辐射场 sf,htbound1,rdsf,0.9,4 !定义边界4为第4个热辐射场 sf,htbound1,rdsf,0.9,5 !定义边界5为第5个热辐射场 stef,5.6696e-8 !定义stefan-boltzmann常数 toffst,273 !定义绝对温度偏差
spctemp,1,temp !定义第1个热辐射场的环境温度 spctemp,2,temp !定义第2个热辐射场的环境温度 spctemp,3,temp !定义第3个热辐射场的环境温度 spctemp,4,temp !定义第4个热辐射场的环境温度 spctemp,5,temp !定义第5个热辐射场的环境温度 !----------------------------------------
solve !求解 *enddo finish
!---------------------------------------- /post1
plnsol,temp,,0, !画出温度分布图 finish
!---------------------------------------- !结构分析 /prep7
/title,Part 2:structural analysis et,1,solid45,1,1 ET,2,beam188
!--------------------------------------- !定义结构分析材料特性 ！-------------------------------
fy=275e6 ！常温下屈服应力 exx=2.1e11 !常温下杨氏模量 mptemp
mptemp,20,,100,250,300,400
mpdata,ex,1,,exx,exx,0.9*exx,0.8*exx,0.7*exx mptemp,,500,500,700,800,900
mpdata,ex,1,,0.6*exx,0.31*exx,0.31*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 temp,,0.02,fy temp,,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.11*fy tbpt,,0.15,0.11*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,2*dis_ver k,6,dis_hor,3*dis_ver
k,7,dis_hor+r_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
!shengchengxian 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,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,linemm11,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
*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,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
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,-4,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,n4,rotz,-dx num=num+1
dx=nx(k+1)-dis_hor
ce,num,0,k+1,uy,1,n4,uy,-1,n4,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,rotxz,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,m16285,17085,100 nsel,a,node,m16286,17086,100 cp,nextmuy,all finish /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,py,-75500 fk,6,py,-151000 lsel,s,line,,12,18,1 sfl,all,pres,-25.00 lsel,all
nsel,s,node,,20084,30084,100 nsel,a,node,,20085,30085,100 sp,all,pres,25400/tw_beam nsel,all time,1
deltim,0,2,1.0e-3,0.5 solve
*do,tm,60,180,60 time,tm
ldread,temp,,,tm,,,rth deltim,20,1,20 solve *enddo finish /post1 plinsol,u,y finish /post26
nsol,2,25005,u,y nsol,3,20004,u,x plnar,2,3
因篇幅问题不能全部显示，请点此查看更多更全内容