[转帖]APDL-------幕墙企业玻璃简化计算

NASA · 发表于 2005-8-13 10:00

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

/encrypt,yuhui4,Magic Glass,mac,d:\ /nopr finish /CLEAR,NOSTART multipro,'start',10 *cset,1,3,H,'The height of glass(unit:m)',2 *cset,4,6,W,'The width of glass(unit:m)',2 *cset,7,9,a1,'Angle left side to down(unit:C)',90 *cset,10,12,a2,'Angle right side to down(unit:C)',90 *cset,13,15,hole,'the number of hole',4 *cset,16,18,type,'Layered Glass Type(!SeePeak!)',1 *cset,19,21,thi,'Thickness of inner Glass(unit:m)',0.01 *cset,22,24,tho,'Thickness of outer Glass(unit:m)',0.008 *cset,25,27,Wk,'Designed wind pressure(unit:N/m2)',1000 *cset,61,62,'Please Fill in related blanks to',' build FEA model And apply loads' *cset,28,30,Seismic,'Designed seismic accelerate(m/s2)',0.08 *cset,63,64,'FOR LAYERED GLASS TYPE:1--[dan','pian],2--[jiajiao],3--[ganghua]' multipro,'end' *if,hole,eq,0,then multipro,'start',8 *cset,1,3,doftop,'constrain type on the top side',2 *cset,4,6,dofdown,'constrain type on the down side',1 *cset,7,9,dofleft,'constrain type on the left side',2 *cset,10,12,dofright,'constrain type on the right side',1 *cset,13,15,Wdiv,'divide number on the width side',40 *cset,16,18,Hdiv,'divide number on the heigth side',40 *cset,19,21,gravity,'The gravity Acceleration(unit:m/s2)',0 *cset,22,24,switch,'Choose solution type',1 *cset,61,62,'Please Fill in related blanks to',' Get FEA solution' *cset,63,64,'SWITCH:0--[By Youself]1--[Small',',Deformed] 2--[LargeDeformed]' multipro,'end' *elseif,hole,eq,2,then multipro,'start',9 *cset,1,3,kr,'Radium of glass hole(unit:m)',0.0175 *cset,4,6,kbw,'Distance to width side(unit:m)',0.106 *cset,7,9,kbh,'Distance to height side(unit:m)',0.106 *cset,10,12,dofdown,'constrain type on the down side',1 *cset,13,15,Wdiv,'divide number on the width side',40 *cset,16,18,Hdiv,'divide number on the heigth side',40 *cset,19,21,holediv,'divide number on the hole side',32 *cset,22,24,gravity,'The gravity Acceleration(unit:m/s2)',0 *cset,25,27,switch,'Choose solution type',1 *cset,61,62,'Please Fill in related blanks to',' Get FEA solution' *cset,63,64,'SWITCH:0--[By Youself]1--[Small',',Deformed] 2--[LargeDeformed]' multipro,'end' *elseif,hole,eq,4,then multipro,'start',8 *cset,1,3,kr,'Radium of glass hole(unit:m)',0.0175 *cset,4,6,kbw,'Distance to width side(unit:m)',0.106 *cset,7,9,kbh,'Distance to height side(unit:m)',0.106 *cset,10,12,Wdiv,'divide number on the width side',40 *cset,13,15,Hdiv,'divide number on the heigth side',40 *cset,16,18,holediv,'divide number on the hole side',32 *cset,19,21,gravity,'The gravity Acceleration(unit:m/s2)',0 *cset,22,24,switch,'Choose solution type',1 *cset,61,62,'Please Fill in related blanks to',' Get FEA solution' *cset,63,64,'SWITCH:0--[By Youself]1--[Small',',Deformed] 2--[LargeDeformed]' multipro,'end' *elseif,hole,eq,6,then multipro,'start',9 *cset,1,3,kr,'Radium of glass hole(unit:m)',0.0175 *cset,4,6,kbw,'Distance to width side(unit:m)',0.106 *cset,7,9,kbh,'Distance to height side(unit:m)',0.106 *cset,10,12,Dbottom,'hole central distance to bottom',1 *cset,13,15,Wdiv,'divide number on the width side',40 *cset,16,18,Hdiv,'divide number on the heigth side',40 *cset,19,21,holediv,'divide number on the hole side',32 *cset,22,24,gravity,'The gravity Acceleration(unit:m/s2)',0 *cset,25,27,switch,'Choose solution type',1 *cset,61,62,'Please Fill in related blanks to',' Get FEA solution' *cset,63,64,'SWITCH:0--[By Youself]1--[Small',',Deformed] 2--[LargeDeformed]' multipro,'end' *endif *if,hole,eq,0,then!!!!!!!!!!平板建模 /prep7 csys,0 k,1,1,1 k,2,1+w,1 *afun,deg k,3,1+H*cos(a1)/sin(a1),1+H k,4,1+w-H*cos(a2)/sin(a2),1+H a,1,2,4,3 et,1,shell63 *if,type,eq,1,then thforcal=1.1*thi *elseif,type,eq,2,then thforcal=1.1*1.25*thi *elseif,type,eq,3,then thforcal=1.1*1.20*thi *endif R,1,thforcal, MP,EX,1,7.2e10 MP,PRXY,1,0.21 mp,dens,1,2560 lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01 lesize,all,,,wdiv lsel,all lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01 lesize,all,,,wdiv lsel,all lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01 lesize,all,,,hdiv lsel,all lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01 lesize,all,,,hdiv lsel,all amesh,1 *if,doftop,eq,0,then *elseif,doftop,eq,1,then lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01 dl,all,,uz dl,all,,uX allsel *elseif,doftop,eq,2,then lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01 dl,all,,uz dl,all,,uy allsel *endif *if,dofdown,eq,0,then *elseif,dofdown,eq,1,then lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01 dl,all,,uz dl,all,,uX allsel *elseif,dofdown,eq,2,then lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01 dl,all,,uz dl,all,,uy allsel *endif *if,dofleft,eq,0,then *elseif,dofleft,eq,1,then lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01 dl,all,,uz dl,all,,uY allsel *elseif,dofleft,eq,2,then lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01 dl,all,,uz dl,all,,ux allsel *endif *if,dofright,eq,0,then *elseif,dofright,eq,1,then lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01 dl,all,,uz dl,all,,uY allsel *elseif,dofright,eq,2,then lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01 dl,all,,uz dl,all,,ux allsel *endif finish *elseif,hole,eq,2,then!!!!!!!!!!!!两点建模 /prep7 /UIS,MSGPOP,3 et,1,shell63 *if,type,eq,1,then thforcal=1.1*thi *elseif,type,eq,2,then thforcal=1.1*1.25*thi *elseif,type,eq,3,then thforcal=1.1*1.20*thi *endif r,2,thforcal mp,ex,2,2.06e11 mp,prxy,2,0.3 csys,1 n,1, n,2,kr, kld=holediv ngen,kld,1,2,,,,360/kld *creat,make,mac et,1,shell63 real,2 mat,2 ! :do loop1 *do,i,2,kld,1 e,1,i,i+1 *enddo finish /prep7 e,1,kld+1,2 *end make finish /prep7 csys,0 k,1,1,1 k,2,1+w,1 *afun,deg k,3,1+H*cos(a1)/sin(a1),1+H k,4,1+w-H*cos(a2)/sin(a2),1+H k,5,kx(3)+kbw/sin(a1)-kbh*cos(a1)/sin(a1),ky(3)-kbh k,6,kx(4)-kbw/sin(a2)+kbw*cos(a2)/sin(a2),ky(4)-kbh a,1,2,4,3 CYL4,kx(5),ky(5),kr CYL4,kx(6),ky(6),kr a***a,1,2 a***a,4,3 aplot /UIS,MSGPOP,3 R,1,thforcal, MP,EX,1,7.2e10 MP,PRXY,1,0.21 mp,dens,1,2560 SMRT,6 TYPE, 1 MAT, 1 REAL, 1 ESYS, 0 lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01 lesize,all,,,wdiv lsel,all lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01 lesize,all,,,wdiv lsel,all lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01 lesize,all,,,hdiv lsel,all lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01 lesize,all,,,hdiv lsel,all LSEL,S,LENGTH,,2*kr,2*(h+w) LSEL,INVE lesize,all,,,kld/4 smrtsize,6 mshape,1,2d mshkey,0 amesh,1 eplot *get,nnn1,node,0,count egen,2,nnn1+kld,1,kld,1,,,,,,kx(5),ky(5) *get,nnn2,node,0,count egen,2,nnn2+kld,1,kld,1,,,,,,kx(6),ky(6) /UIS,MSGPOP,3 edele,1,kld,1 ndele,1,kld+1 /auto,1 /rep nsel,s,loc,x,kx(5)-0.001,kx(5)+0.001 nsel,r,loc,y,ky(5)-0.001,ky(5)+0.001 d,all,ux,0,,,,uy,uz allsel nsel,s,loc,x,kx(6)-0.001,kx(6)+0.001 nsel,r,loc,y,ky(6)-0.001,ky(6)+0.001 d,all,uy,0,,,,uz allsel CPINTF,ALL,0.0001, *if,dofdown,eq,0,then *elseif,dofdown,eq,1,then lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01 dl,all,,uz allsel *elseif,dofdown,eq,2,then lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01 dl,all,,uz dl,all,,uy allsel *endif finish *elseif,hole,eq,4,then!!!!!!!!!!!!!!!!!四点建模 /prep7 /UIS,MSGPOP,3 et,1,shell63 *if,type,eq,1,then thforcal=1.1*thi *elseif,type,eq,2,then thforcal=1.1*1.25*thi *elseif,type,eq,3,then thforcal=1.1*1.20*thi *endif

NASA · 发表于 2005-8-13 10:00

r,2,thforcal mp,ex,2,2.06e11 mp,prxy,2,0.3 csys,1 n,1, n,2,kr, kld=holediv ngen,kld,1,2,,,,360/kld *creat,make,mac et,2,shell63 real,2 mat,2 ! :do loop1 *do,i,2,kld,1 e,1,i,i+1 *enddo finish /prep7 e,1,kld+1,2 *end make finish /prep7 csys,0 k,1,1,1 k,2,1+w,1 *afun,deg k,3,1+H*cos(a1)/sin(a1),1+H k,4,1+w-H*cos(a2)/sin(a2),1+H k,5,kbw/sin(a1)+kbh*cos(a1)/sin(a1)+kx(1),kbh+ky(1) k,6,kx(2)-kbw/sin(a2)+kbh*cos(a2)/sin(a2),kbh+ky(2) k,7,kx(3)+kbw/sin(a1)-kbh*cos(a1)/sin(a1),ky(3)-kbh k,8,kx(4)-kbw/sin(a2)+kbw*cos(a2)/sin(a2),ky(4)-kbh a,1,2,4,3 CYL4,kx(5),ky(5),kr CYL4,kx(6),ky(6),kr CYL4,kx(7),ky(7),kr CYL4,kx(8),ky(8),kr a***a,1,2 a***a,6,3 a***a,1,4 a***a,2,5 aplot /UIS,MSGPOP,3 R,1,thforcal, MP,EX,1,7.2e10 MP,PRXY,1,0.21 mp,dens,1,2560 SMRT,6 TYPE, 1 MAT, 1 REAL, 1 ESYS, 0 lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01 lesize,all,,,wdiv lsel,all lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01 lesize,all,,,wdiv lsel,all lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01 lesize,all,,,hdiv lsel,all lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01 lesize,all,,,hdiv lsel,all LSEL,S,LENGTH,,2*kr,l+w LSEL,INVE lesize,all,,,kld/4 smrtsize,6 mshape,1,2d mshkey,0 amesh,1 eplot *get,nnn1,node,0,count egen,2,nnn1+kld,1,kld,1,,,,,,kx(5),ky(5) *get,nnn2,node,0,count egen,2,nnn2+kld,1,kld,1,,,,,,kx(6),ky(6) *get,nnn3,node,0,count egen,2,nnn3+kld,1,kld,1,,,,,,kx(7),ky(7) *get,nnn4,node,0,count egen,2,nnn4+kld,1,kld,1,,,,,,kx(8),ky(8) /UIS,MSGPOP,3 edele,1,kld,1 ndele,1,kld+1 /auto,1 /rep finish /solu allsel nsel,s,loc,x,kx(5)-0.001,kx(5)+0.001 nsel,r,loc,y,ky(5)-0.001,ky(5)+0.001 d,all,ux,0,,,,uy,uz allsel nsel,s,loc,x,kx(6)-0.001,kx(6)+0.001 nsel,r,loc,y,ky(6)-0.001,ky(6)+0.001 d,all,uy,0,,,,uz allsel nsel,s,loc,x,kx(7)-0.001,kx(7)+0.001 nsel,r,loc,y,ky(7)-0.001,ky(7)+0.001 d,all,ux,0,,,,uz allsel nsel,s,loc,x,kx(8)-0.001,kx(8)+0.001 nsel,r,loc,y,ky(8)-0.001,ky(8)+0.001 d,all,uz,0,,,, allsel CPINTF,ALL,0.0001, finish gplot *elseif,hole,eq,6,then!!!!!!!!!!!!!!!!!六点建模 /prep7 /UIS,MSGPOP,3 et,1,shell63 *if,type,eq,1,then thforcal=1.1*thi *elseif,type,eq,2,then thforcal=1.1*1.25*thi *elseif,type,eq,3,then thforcal=1.1*1.20*thi *endif r,2,thforcal mp,ex,2,2.06e11 mp,prxy,2,0.3 csys,1 n,1, n,2,kr, kld=holediv ngen,kld,1,2,,,,360/kld *creat,make,mac et,2,shell63 real,2 mat,2 ! :do loop1 *do,i,2,kld,1 e,1,i,i+1 *enddo finish /prep7 e,1,kld+1,2 *end make finish /prep7 csys,0 k,1,1,1 k,2,1+w,1 *afun,deg k,3,1+H*cos(a1)/sin(a1),1+H k,4,1+w-H*cos(a2)/sin(a2),1+H k,5,kbw/sin(a1)+kbh*cos(a1)/sin(a1)+kx(1),kbh+ky(1) k,6,kx(2)-kbw/sin(a2)+kbh*cos(a2)/sin(a2),kbh+ky(2) k,7,kx(3)+kbw/sin(a1)-kbh*cos(a1)/sin(a1),ky(3)-kbh k,8,kx(4)-kbw/sin(a2)+kbw*cos(a2)/sin(a2),ky(4)-kbh hp=ky(7)-ky(5) xx9=(hp-dbottom)*(kx(7)-kx(5))/hp yy9=ky(5)+dbottom xx10=(hp-dbottom)*(kx(6)-kx(8))/hp yy10=ky(6)+dbottom k,9,kx(7)-xx9,yy9 k,10,kx(8)+xx10,yy10 a,1,2,4,3 CYL4,kx(5),ky(5),kr CYL4,kx(6),ky(6),kr CYL4,kx(7),ky(7),kr CYL4,kx(8),ky(8),kr CYL4,kx(9),ky(9),kr CYL4,kx(10),ky(10),kr a***a,1,2 a***a,8,3 a***a,1,4 a***a,2,5 a***a,1,6 a***a,2,7 aplot /UIS,MSGPOP,3 R,1,thforcal, MP,EX,1,7.2e10 MP,PRXY,1,0.21 mp,dens,1,2560 SMRT,6 TYPE, 1 MAT, 1 REAL, 1 ESYS, 0 lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01 lesize,all,,,wdiv lsel,all lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01 lesize,all,,,wdiv lsel,all lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01 lesize,all,,,hdiv lsel,all lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01 lesize,all,,,hdiv lsel,all LSEL,S,LENGTH,,2*kr,l+w LSEL,INVE lesize,all,,,kld/4 smrtsize,6 mshape,1,2d mshkey,0 amesh,1 eplot *get,nnn1,node,0,count egen,2,nnn1+kld,1,kld,1,,,,,,kx(5),ky(5) *get,nnn2,node,0,count egen,2,nnn2+kld,1,kld,1,,,,,,kx(6),ky(6) *get,nnn3,node,0,count egen,2,nnn3+kld,1,kld,1,,,,,,kx(7),ky(7) *get,nnn4,node,0,count egen,2,nnn4+kld,1,kld,1,,,,,,kx(8),ky(8) *get,nnn5,node,0,count egen,2,nnn5+kld,1,kld,1,,,,,,kx(9),ky(9) *get,nnn6,node,0,count egen,2,nnn6+kld,1,kld,1,,,,,,kx(10),ky(10) /UIS,MSGPOP,3 edele,1,kld,1 ndele,1,kld+1 /auto,1 /rep finish /solu allsel nsel,s,loc,x,kx(5)-0.001,kx(5)+0.001 nsel,r,loc,y,ky(5)-0.001,ky(5)+0.001 d,all,ux,0,,,,uy,uz allsel nsel,s,loc,x,kx(6)-0.001,kx(6)+0.001 nsel,r,loc,y,ky(6)-0.001,ky(6)+0.001 d,all,uy,0,,,,uz allsel nsel,s,loc,x,kx(7)-0.001,kx(7)+0.001 nsel,r,loc,y,ky(7)-0.001,ky(7)+0.001 d,all,ux,0,,,,uz allsel nsel,s,loc,x,kx(8)-0.001,kx(8)+0.001 nsel,r,loc,y,ky(8)-0.001,ky(8)+0.001 d,all,uz,0,,,, allsel nsel,s,loc,x,kx(9)-0.001,kx(9)+0.001 nsel,r,loc,y,ky(9)-0.001,ky(9)+0.001 d,all,ux,0,,,,uz allsel nsel,s,loc,x,kx(10)-0.001,kx(10)+0.001 nsel,r,loc,y,ky(10)-0.001,ky(10)+0.001 d,all,uz,0,,,, allsel CPINTF,ALL,0.0001, finish gplot *endif finish

NASA · 发表于 2005-8-13 10:00

/solu!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!求解部分 Ek=3*seismic*(thi+tho)*26.0*1000*9.8 w=1.4*wk e=1.3*ek q1=0.6*ek+wk q2=w+0.6*e finish /solu /UIS,MSGPOP,3 esel,s,real,,1,1,1 eplot sfe,all,1,pres,,q1 acel,,gravity, allsel gplot *if,switch,eq,1,then nlgeom,off kbc,1 *elseif,switch,eq,2,then nlgeom,on sstif,on nsubst,10 kbc,1 *endif solv finish /post1 esel,s,real,,1,1,1 AVPRIN plnsol,u,sum,0,1 *get,numall,node,0,count *get,numstart,node,0,num,min umax=0 inuse=numstart *do,i,1,numall,1 *get,usum,node,inuse,u,sum *if,umax,ge,usum,then umax=umax *elseif,umax,lt,usum,then umax=usum *endif inuse=ndnext(inuse) *enddo finish /solu /UIS,MSGPOP,3 esel,s,real,,1,1,1 eplot sfe,all,1,pres,,q2 allsel gplot *if,switch,eq,1,then nlgeom,off kbc,1 solv finish /post1 esel,s,real,,1,1,1 AVPRIN plnsol,s,eqv,0,1 allsel esel,s,real,,1,1,1 NSLE,S *get,numall,node,0,count *get,numstart,node,0,num,min smax=0 inuse=numstart *do,i,1,numall,1 *get,ssum,node,inuse,s,eqv *if,smax,ge,ssum,then smax=smax *elseif,smax,lt,ssum,then smax=ssum *endif inuse=ndnext(inuse) *enddo allsel *if,thforcal,le,0.012,then dmq=84 kbq=58.8 *elseif,thforcal,gt,0.012,then dmq=59 kbq=41.3 *endif *if,h,lt,w,then duanb=h *elseif,h,ge,w,then duanb=w *endif str=smax/1e6 esel,s,real,,1,1,1 /UIS,MSGPOP,1 *msg,note,str,umax*1000,dmq,10*duanb Ultimate Limit States Glass STRESS=%gMa,%/& Serviceabilitylimits DISPLACEMENT=%gmm,%/& [Toughened Glass Allowble Stress]=%gMa,%/& [Control Displacement]=%gmm, /UIS,MSGPOP,3 *elseif,switch,eq,2,then nlgeom,on sstif,on nsubst,10 kbc,1 solv finish /post1 esel,s,real,,1,1,1 AVPRIN plnsol,s,eqv,0,1 allsel esel,s,real,,1,1,1 NSLE,S *get,numall,node,0,count *get,numstart,node,0,num,min smax=0 inuse=numstart *do,i,1,numall,1 *get,ssum,node,inuse,s,eqv *if,smax,ge,ssum,then smax=smax *elseif,smax,lt,ssum,then smax=ssum *endif inuse=ndnext(inuse) *enddo allsel *if,thforcal,le,0.012,then dmq=84 kbq=58.8 *elseif,thforcal,gt,0.012,then dmq=59 kbq=41.3 *endif *if,h,lt,w,then duanb=h *elseif,h,ge,w,then duanb=w *endif str=smax/1e6 esel,s,real,,1,1,1 /UIS,MSGPOP,1 *msg,note,str,umax*1000,dmq,10*duanb Ultimate Limit States Glass STRESS=%gMa,%/& Serviceabilitylimits DISPLACEMENT=%gmm,%/& [Toughened Glass Allowble Stress]=%gMa,%/& [Control Displacement]=%gmm, /UIS,MSGPOP,3 *elseif,swtich,eq,0 allsel gplot *endif /gopr finish /encrypt

账号		自动登录	找回密码
密码			我要加入

[综合讨论] [转帖]APDL-------幕墙企业玻璃简化计算

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（oscar32）[转帖]APDL-------幕墙企业玻璃简化...

回复：（oscar32）[转帖]APDL-------幕墙企业玻璃简化...