Ansys讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断

AaronSpark 发表于 2006-5-12 02:08

[转帖]钢框架受热分析

此例子是选自《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

页: [1]

声振论坛's Archiver

[转帖]钢框架受热分析