[转帖]APDL-------幕墙企业玻璃简化计算
/encrypt,yuhui4,Magic Glass,mac,d:\ <BR>/nopr <BR>finish <BR>/CLEAR,NOSTART <BR>multipro,'start',10 <BR>*cset,1,3,H,'The height of glass(unit:m)',2 <BR>*cset,4,6,W,'The width of glass(unit:m)',2 <BR>*cset,7,9,a1,'Angle left side to down(unit:C)',90 <BR>*cset,10,12,a2,'Angle right side to down(unit:C)',90 <BR>*cset,13,15,hole,'the number of hole',4 <BR>*cset,16,18,type,'Layered Glass Type(!SeePeak!)',1 <BR>*cset,19,21,thi,'Thickness of inner Glass(unit:m)',0.01 <BR>*cset,22,24,tho,'Thickness of outer Glass(unit:m)',0.008 <BR>*cset,25,27,Wk,'Designed wind pressure(unit:N/m2)',1000 <BR>*cset,61,62,'Please Fill in related blanks to',' build FEA model And apply loads' <BR>*cset,28,30,Seismic,'Designed seismic accelerate(m/s2)',0.08 <BR>*cset,63,64,'FOR LAYERED GLASS TYPE:1--,2--,3--' <BR>multipro,'end' <BR>*if,hole,eq,0,then <BR> multipro,'start',8 <BR> *cset,1,3,doftop,'constrain type on the top side',2 <BR> *cset,4,6,dofdown,'constrain type on the down side',1 <BR> *cset,7,9,dofleft,'constrain type on the left side',2 <BR> *cset,10,12,dofright,'constrain type on the right side',1 <BR> *cset,13,15,Wdiv,'divide number on the width side',40 <BR> *cset,16,18,Hdiv,'divide number on the heigth side',40 <BR> *cset,19,21,gravity,'The gravity Acceleration(unit:m/s2)',0 <BR> *cset,22,24,switch,'Choose solution type',1 <BR> *cset,61,62,'Please Fill in related blanks to',' Get FEA solution' <BR> *cset,63,64,'SWITCH:0--1-- 2--' <BR> multipro,'end' <BR>*elseif,hole,eq,2,then <BR> multipro,'start',9 <BR> *cset,1,3,kr,'Radium of glass hole(unit:m)',0.0175 <BR> *cset,4,6,kbw,'Distance to width side(unit:m)',0.106 <BR> *cset,7,9,kbh,'Distance to height side(unit:m)',0.106 <BR> *cset,10,12,dofdown,'constrain type on the down side',1 <BR> *cset,13,15,Wdiv,'divide number on the width side',40 <BR> *cset,16,18,Hdiv,'divide number on the heigth side',40 <BR> *cset,19,21,holediv,'divide number on the hole side',32 <BR> *cset,22,24,gravity,'The gravity Acceleration(unit:m/s2)',0 <BR> *cset,25,27,switch,'Choose solution type',1 <BR> *cset,61,62,'Please Fill in related blanks to',' Get FEA solution' <BR> *cset,63,64,'SWITCH:0--1-- 2--' <BR> multipro,'end' <BR>*elseif,hole,eq,4,then <BR> multipro,'start',8 <BR> *cset,1,3,kr,'Radium of glass hole(unit:m)',0.0175 <BR> *cset,4,6,kbw,'Distance to width side(unit:m)',0.106 <BR> *cset,7,9,kbh,'Distance to height side(unit:m)',0.106 <BR> *cset,10,12,Wdiv,'divide number on the width side',40 <BR> *cset,13,15,Hdiv,'divide number on the heigth side',40 <BR> *cset,16,18,holediv,'divide number on the hole side',32 <BR> *cset,19,21,gravity,'The gravity Acceleration(unit:m/s2)',0 <BR> *cset,22,24,switch,'Choose solution type',1 <BR> *cset,61,62,'Please Fill in related blanks to',' Get FEA solution' <BR> *cset,63,64,'SWITCH:0--1-- 2--' <BR> multipro,'end' <BR>*elseif,hole,eq,6,then <BR> multipro,'start',9 <BR> *cset,1,3,kr,'Radium of glass hole(unit:m)',0.0175 <BR> *cset,4,6,kbw,'Distance to width side(unit:m)',0.106 <BR> *cset,7,9,kbh,'Distance to height side(unit:m)',0.106 <BR> *cset,10,12,Dbottom,'hole central distance to bottom',1 <BR> *cset,13,15,Wdiv,'divide number on the width side',40 <BR> *cset,16,18,Hdiv,'divide number on the heigth side',40 <BR> *cset,19,21,holediv,'divide number on the hole side',32 <BR> *cset,22,24,gravity,'The gravity Acceleration(unit:m/s2)',0 <BR> *cset,25,27,switch,'Choose solution type',1 <BR> *cset,61,62,'Please Fill in related blanks to',' Get FEA solution' <BR> *cset,63,64,'SWITCH:0--1-- 2--' <BR> multipro,'end' <BR>*endif <BR>*if,hole,eq,0,then!!!!!!!!!!平板建模 <BR> /prep7 <BR> csys,0 <BR> k,1,1,1 <BR> k,2,1+w,1 <BR> *afun,deg <BR> k,3,1+H*cos(a1)/sin(a1),1+H <BR> k,4,1+w-H*cos(a2)/sin(a2),1+H <BR> a,1,2,4,3 <BR> et,1,shell63 <BR>*if,type,eq,1,then <BR>thforcal=1.1*thi <BR>*elseif,type,eq,2,then <BR>thforcal=1.1*1.25*thi <BR>*elseif,type,eq,3,then <BR>thforcal=1.1*1.20*thi <BR>*endif <BR>R,1,thforcal, <BR>MP,EX,1,7.2e10 <BR>MP,PRXY,1,0.21 <BR>mp,dens,1,2560 <BR>lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01 <BR> lesize,all,,,wdiv <BR> lsel,all <BR>lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01 <BR> lesize,all,,,wdiv <BR> lsel,all <BR>lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01 <BR> lesize,all,,,hdiv <BR> lsel,all <BR>lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01 <BR> lesize,all,,,hdiv <BR> lsel,all <BR>amesh,1 <BR>*if,doftop,eq,0,then <BR> *elseif,doftop,eq,1,then <BR> lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01 <BR> dl,all,,uz <BR> dl,all,,uX <BR> allsel <BR> *elseif,doftop,eq,2,then <BR> lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01 <BR> dl,all,,uz <BR> dl,all,,uy <BR> allsel <BR>*endif <BR>*if,dofdown,eq,0,then <BR> *elseif,dofdown,eq,1,then <BR> lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01 <BR> dl,all,,uz <BR> dl,all,,uX <BR> allsel <BR> *elseif,dofdown,eq,2,then <BR> lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01 <BR> dl,all,,uz <BR> dl,all,,uy <BR> allsel <BR>*endif <BR>*if,dofleft,eq,0,then <BR> *elseif,dofleft,eq,1,then <BR> lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01 <BR> dl,all,,uz <BR> dl,all,,uY <BR> allsel <BR> *elseif,dofleft,eq,2,then <BR> lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01 <BR> dl,all,,uz <BR> dl,all,,ux <BR> allsel <BR>*endif <BR>*if,dofright,eq,0,then <BR> *elseif,dofright,eq,1,then <BR> lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01 <BR> dl,all,,uz <BR> dl,all,,uY <BR> allsel <BR> *elseif,dofright,eq,2,then <BR> lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01 <BR> dl,all,,uz <BR> dl,all,,ux <BR> allsel <BR>*endif <BR>finish <BR>*elseif,hole,eq,2,then!!!!!!!!!!!!两点建模 <BR>/prep7 <BR>/UIS,MSGPOP,3 <BR>et,1,shell63 <BR>*if,type,eq,1,then <BR>thforcal=1.1*thi <BR>*elseif,type,eq,2,then <BR>thforcal=1.1*1.25*thi <BR>*elseif,type,eq,3,then <BR>thforcal=1.1*1.20*thi <BR>*endif <BR>r,2,thforcal <BR>mp,ex,2,2.06e11 <BR>mp,prxy,2,0.3 <BR>csys,1 <BR>n,1, <BR>n,2,kr, <BR>kld=holediv <BR>ngen,kld,1,2,,,,360/kld <BR>*creat,make,mac <BR>et,1,shell63 <BR>real,2 <BR>mat,2 <BR>! :do loop1 <BR>*do,i,2,kld,1 <BR>e,1,i,i+1 <BR>*enddo <BR>finish <BR>/prep7 <BR>e,1,kld+1,2 <BR>*end <BR>make <BR>finish <BR>/prep7 <BR>csys,0 <BR>k,1,1,1 <BR>k,2,1+w,1 <BR>*afun,deg <BR>k,3,1+H*cos(a1)/sin(a1),1+H <BR>k,4,1+w-H*cos(a2)/sin(a2),1+H <BR>k,5,kx(3)+kbw/sin(a1)-kbh*cos(a1)/sin(a1),ky(3)-kbh <BR>k,6,kx(4)-kbw/sin(a2)+kbw*cos(a2)/sin(a2),ky(4)-kbh <BR>a,1,2,4,3 <BR>CYL4,kx(5),ky(5),kr <BR>CYL4,kx(6),ky(6),kr <BR>a***a,1,2 <BR>a***a,4,3 <BR>aplot <BR>/UIS,MSGPOP,3 <BR>R,1,thforcal, <BR>MP,EX,1,7.2e10 <BR>MP,PRXY,1,0.21 <BR>mp,dens,1,2560 <BR>SMRT,6 <BR>TYPE, 1 <BR>MAT, 1 <BR>REAL, 1 <BR>ESYS, 0 <BR>lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01 <BR> lesize,all,,,wdiv <BR> lsel,all <BR>lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01 <BR> lesize,all,,,wdiv <BR> lsel,all <BR>lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01 <BR> lesize,all,,,hdiv <BR> lsel,all <BR>lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01 <BR> lesize,all,,,hdiv <BR> lsel,all <BR>LSEL,S,LENGTH,,2*kr,2*(h+w) <BR>LSEL,INVE <BR>lesize,all,,,kld/4 <BR>smrtsize,6 <BR>mshape,1,2d <BR>mshkey,0 <BR>amesh,1 <BR>eplot <BR>*get,nnn1,node,0,count <BR>egen,2,nnn1+kld,1,kld,1,,,,,,kx(5),ky(5) <BR>*get,nnn2,node,0,count <BR>egen,2,nnn2+kld,1,kld,1,,,,,,kx(6),ky(6) <BR>/UIS,MSGPOP,3 <BR>edele,1,kld,1 <BR>ndele,1,kld+1 <BR>/auto,1 <BR>/rep <BR>nsel,s,loc,x,kx(5)-0.001,kx(5)+0.001 <BR>nsel,r,loc,y,ky(5)-0.001,ky(5)+0.001 <BR>d,all,ux,0,,,,uy,uz <BR>allsel <BR>nsel,s,loc,x,kx(6)-0.001,kx(6)+0.001 <BR>nsel,r,loc,y,ky(6)-0.001,ky(6)+0.001 <BR>d,all,uy,0,,,,uz <BR>allsel <BR>CPINTF,ALL,0.0001, <BR>*if,dofdown,eq,0,then <BR> *elseif,dofdown,eq,1,then <BR> lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01 <BR> dl,all,,uz <BR> allsel <BR> *elseif,dofdown,eq,2,then <BR> lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01 <BR> dl,all,,uz <BR> dl,all,,uy <BR> allsel <BR>*endif <BR>finish <BR>*elseif,hole,eq,4,then!!!!!!!!!!!!!!!!!四点建模 <BR>/prep7 <BR>/UIS,MSGPOP,3 <BR>et,1,shell63 <BR>*if,type,eq,1,then <BR>thforcal=1.1*thi <BR>*elseif,type,eq,2,then <BR>thforcal=1.1*1.25*thi <BR>*elseif,type,eq,3,then <BR>thforcal=1.1*1.20*thi <BR>*endif <BR>回复：（oscar32）[转帖]APDL-------幕墙企业玻璃简化...
r,2,thforcal <BR>mp,ex,2,2.06e11 <BR>mp,prxy,2,0.3 <BR>csys,1 <BR>n,1, <BR>n,2,kr, <BR>kld=holediv <BR>ngen,kld,1,2,,,,360/kld <BR>*creat,make,mac <BR>et,2,shell63 <BR>real,2 <BR>mat,2 <BR>! :do loop1 <BR>*do,i,2,kld,1 <BR>e,1,i,i+1 <BR>*enddo <BR>finish <BR>/prep7 <BR>e,1,kld+1,2 <BR>*end <BR>make <BR>finish <BR>/prep7 <BR>csys,0 <BR>k,1,1,1 <BR>k,2,1+w,1 <BR>*afun,deg <BR>k,3,1+H*cos(a1)/sin(a1),1+H <BR>k,4,1+w-H*cos(a2)/sin(a2),1+H <BR>k,5,kbw/sin(a1)+kbh*cos(a1)/sin(a1)+kx(1),kbh+ky(1) <BR>k,6,kx(2)-kbw/sin(a2)+kbh*cos(a2)/sin(a2),kbh+ky(2) <BR>k,7,kx(3)+kbw/sin(a1)-kbh*cos(a1)/sin(a1),ky(3)-kbh <BR>k,8,kx(4)-kbw/sin(a2)+kbw*cos(a2)/sin(a2),ky(4)-kbh <BR>a,1,2,4,3 <BR>CYL4,kx(5),ky(5),kr <BR>CYL4,kx(6),ky(6),kr <BR>CYL4,kx(7),ky(7),kr <BR>CYL4,kx(8),ky(8),kr <BR>a***a,1,2 <BR>a***a,6,3 <BR>a***a,1,4 <BR>a***a,2,5 <BR>aplot <BR>/UIS,MSGPOP,3 <BR>R,1,thforcal, <BR>MP,EX,1,7.2e10 <BR>MP,PRXY,1,0.21 <BR>mp,dens,1,2560 <BR>SMRT,6 <BR>TYPE, 1 <BR>MAT, 1 <BR>REAL, 1 <BR>ESYS, 0 <BR>lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01 <BR> lesize,all,,,wdiv <BR> lsel,all <BR>lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01 <BR> lesize,all,,,wdiv <BR> lsel,all <BR>lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01 <BR> lesize,all,,,hdiv <BR> lsel,all <BR>lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01 <BR> lesize,all,,,hdiv <BR> lsel,all <BR>LSEL,S,LENGTH,,2*kr,l+w <BR>LSEL,INVE <BR>lesize,all,,,kld/4 <BR>smrtsize,6 <BR>mshape,1,2d <BR>mshkey,0 <BR>amesh,1 <BR>eplot <BR>*get,nnn1,node,0,count <BR>egen,2,nnn1+kld,1,kld,1,,,,,,kx(5),ky(5) <BR>*get,nnn2,node,0,count <BR>egen,2,nnn2+kld,1,kld,1,,,,,,kx(6),ky(6) <BR>*get,nnn3,node,0,count <BR>egen,2,nnn3+kld,1,kld,1,,,,,,kx(7),ky(7) <BR>*get,nnn4,node,0,count <BR>egen,2,nnn4+kld,1,kld,1,,,,,,kx(8),ky(8) <BR>/UIS,MSGPOP,3 <BR>edele,1,kld,1 <BR>ndele,1,kld+1 <BR>/auto,1 <BR>/rep <BR>finish <BR>/solu <BR>allsel <BR>nsel,s,loc,x,kx(5)-0.001,kx(5)+0.001 <BR>nsel,r,loc,y,ky(5)-0.001,ky(5)+0.001 <BR>d,all,ux,0,,,,uy,uz <BR>allsel <BR>nsel,s,loc,x,kx(6)-0.001,kx(6)+0.001 <BR>nsel,r,loc,y,ky(6)-0.001,ky(6)+0.001 <BR>d,all,uy,0,,,,uz <BR>allsel <BR>nsel,s,loc,x,kx(7)-0.001,kx(7)+0.001 <BR>nsel,r,loc,y,ky(7)-0.001,ky(7)+0.001 <BR>d,all,ux,0,,,,uz <BR>allsel <BR>nsel,s,loc,x,kx(8)-0.001,kx(8)+0.001 <BR>nsel,r,loc,y,ky(8)-0.001,ky(8)+0.001 <BR>d,all,uz,0,,,, <BR>allsel <BR>CPINTF,ALL,0.0001, <BR>finish <BR>gplot <BR>*elseif,hole,eq,6,then!!!!!!!!!!!!!!!!!六点建模 <BR>/prep7 <BR>/UIS,MSGPOP,3 <BR>et,1,shell63 <BR>*if,type,eq,1,then <BR>thforcal=1.1*thi <BR>*elseif,type,eq,2,then <BR>thforcal=1.1*1.25*thi <BR>*elseif,type,eq,3,then <BR>thforcal=1.1*1.20*thi <BR>*endif <BR>r,2,thforcal <BR>mp,ex,2,2.06e11 <BR>mp,prxy,2,0.3 <BR>csys,1 <BR>n,1, <BR>n,2,kr, <BR>kld=holediv <BR>ngen,kld,1,2,,,,360/kld <BR>*creat,make,mac <BR>et,2,shell63 <BR>real,2 <BR>mat,2 <BR>! :do loop1 <BR>*do,i,2,kld,1 <BR>e,1,i,i+1 <BR>*enddo <BR>finish <BR>/prep7 <BR>e,1,kld+1,2 <BR>*end <BR>make <BR>finish <BR>/prep7 <BR>csys,0 <BR>k,1,1,1 <BR>k,2,1+w,1 <BR>*afun,deg <BR>k,3,1+H*cos(a1)/sin(a1),1+H <BR>k,4,1+w-H*cos(a2)/sin(a2),1+H <BR>k,5,kbw/sin(a1)+kbh*cos(a1)/sin(a1)+kx(1),kbh+ky(1) <BR>k,6,kx(2)-kbw/sin(a2)+kbh*cos(a2)/sin(a2),kbh+ky(2) <BR>k,7,kx(3)+kbw/sin(a1)-kbh*cos(a1)/sin(a1),ky(3)-kbh <BR>k,8,kx(4)-kbw/sin(a2)+kbw*cos(a2)/sin(a2),ky(4)-kbh <BR>hp=ky(7)-ky(5) <BR>xx9=(hp-dbottom)*(kx(7)-kx(5))/hp <BR>yy9=ky(5)+dbottom <BR>xx10=(hp-dbottom)*(kx(6)-kx(8))/hp <BR>yy10=ky(6)+dbottom <BR>k,9,kx(7)-xx9,yy9 <BR>k,10,kx(8)+xx10,yy10 <BR>a,1,2,4,3 <BR>CYL4,kx(5),ky(5),kr <BR>CYL4,kx(6),ky(6),kr <BR>CYL4,kx(7),ky(7),kr <BR>CYL4,kx(8),ky(8),kr <BR>CYL4,kx(9),ky(9),kr <BR>CYL4,kx(10),ky(10),kr <BR>a***a,1,2 <BR>a***a,8,3 <BR>a***a,1,4 <BR>a***a,2,5 <BR>a***a,1,6 <BR>a***a,2,7 <BR>aplot <BR>/UIS,MSGPOP,3 <BR>R,1,thforcal, <BR>MP,EX,1,7.2e10 <BR>MP,PRXY,1,0.21 <BR>mp,dens,1,2560 <BR>SMRT,6 <BR>TYPE, 1 <BR>MAT, 1 <BR>REAL, 1 <BR>ESYS, 0 <BR>lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01 <BR> lesize,all,,,wdiv <BR> lsel,all <BR>lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01 <BR> lesize,all,,,wdiv <BR> lsel,all <BR>lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01 <BR> lesize,all,,,hdiv <BR> lsel,all <BR>lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01 <BR> lesize,all,,,hdiv <BR> lsel,all <BR>LSEL,S,LENGTH,,2*kr,l+w <BR>LSEL,INVE <BR>lesize,all,,,kld/4 <BR>smrtsize,6 <BR>mshape,1,2d <BR>mshkey,0 <BR>amesh,1 <BR>eplot <BR>*get,nnn1,node,0,count <BR>egen,2,nnn1+kld,1,kld,1,,,,,,kx(5),ky(5) <BR>*get,nnn2,node,0,count <BR>egen,2,nnn2+kld,1,kld,1,,,,,,kx(6),ky(6) <BR>*get,nnn3,node,0,count <BR>egen,2,nnn3+kld,1,kld,1,,,,,,kx(7),ky(7) <BR>*get,nnn4,node,0,count <BR>egen,2,nnn4+kld,1,kld,1,,,,,,kx(8),ky(8) <BR>*get,nnn5,node,0,count <BR>egen,2,nnn5+kld,1,kld,1,,,,,,kx(9),ky(9) <BR>*get,nnn6,node,0,count <BR>egen,2,nnn6+kld,1,kld,1,,,,,,kx(10),ky(10) <BR>/UIS,MSGPOP,3 <BR>edele,1,kld,1 <BR>ndele,1,kld+1 <BR>/auto,1 <BR>/rep <BR>finish <BR>/solu <BR>allsel <BR>nsel,s,loc,x,kx(5)-0.001,kx(5)+0.001 <BR>nsel,r,loc,y,ky(5)-0.001,ky(5)+0.001 <BR>d,all,ux,0,,,,uy,uz <BR>allsel <BR>nsel,s,loc,x,kx(6)-0.001,kx(6)+0.001 <BR>nsel,r,loc,y,ky(6)-0.001,ky(6)+0.001 <BR>d,all,uy,0,,,,uz <BR>allsel <BR>nsel,s,loc,x,kx(7)-0.001,kx(7)+0.001 <BR>nsel,r,loc,y,ky(7)-0.001,ky(7)+0.001 <BR>d,all,ux,0,,,,uz <BR>allsel <BR>nsel,s,loc,x,kx(8)-0.001,kx(8)+0.001 <BR>nsel,r,loc,y,ky(8)-0.001,ky(8)+0.001 <BR>d,all,uz,0,,,, <BR>allsel <BR>nsel,s,loc,x,kx(9)-0.001,kx(9)+0.001 <BR>nsel,r,loc,y,ky(9)-0.001,ky(9)+0.001 <BR>d,all,ux,0,,,,uz <BR>allsel <BR>nsel,s,loc,x,kx(10)-0.001,kx(10)+0.001 <BR>nsel,r,loc,y,ky(10)-0.001,ky(10)+0.001 <BR>d,all,uz,0,,,, <BR>allsel <BR>CPINTF,ALL,0.0001, <BR>finish <BR>gplot <BR>*endif <BR>finish <BR>回复：（oscar32）[转帖]APDL-------幕墙企业玻璃简化...
/solu!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!求解部分 <BR>Ek=3*seismic*(thi+tho)*26.0*1000*9.8 <BR>w=1.4*wk <BR>e=1.3*ek <BR>q1=0.6*ek+wk <BR>q2=w+0.6*e <BR>finish <BR>/solu <BR>/UIS,MSGPOP,3 <BR>esel,s,real,,1,1,1 <BR>eplot <BR>sfe,all,1,pres,,q1 <BR>acel,,gravity, <BR>allsel <BR>gplot <BR>*if,switch,eq,1,then <BR>nlgeom,off <BR>kbc,1 <BR>*elseif,switch,eq,2,then <BR>nlgeom,on <BR>sstif,on <BR>nsubst,10 <BR>kbc,1 <BR>*endif <BR>solv <BR>finish <BR>/post1 <BR>esel,s,real,,1,1,1 <BR>AVPRIN <BR>plnsol,u,sum,0,1 <BR>*get,numall,node,0,count <BR>*get,numstart,node,0,num,min <BR> umax=0 <BR> inuse=numstart <BR>*do,i,1,numall,1 <BR> *get,usum,node,inuse,u,sum <BR> *if,umax,ge,usum,then <BR> umax=umax <BR> *elseif,umax,lt,usum,then <BR> umax=usum <BR> *endif <BR>inuse=ndnext(inuse) <BR>*enddo <BR>finish <BR>/solu <BR>/UIS,MSGPOP,3 <BR>esel,s,real,,1,1,1 <BR>eplot <BR>sfe,all,1,pres,,q2 <BR>allsel <BR>gplot <BR>*if,switch,eq,1,then <BR>nlgeom,off <BR>kbc,1 <BR>solv <BR>finish <BR>/post1 <BR>esel,s,real,,1,1,1 <BR>AVPRIN <BR>plnsol,s,eqv,0,1 <BR>allsel <BR>esel,s,real,,1,1,1 <BR>NSLE,S <BR>*get,numall,node,0,count <BR>*get,numstart,node,0,num,min <BR> smax=0 <BR> inuse=numstart <BR>*do,i,1,numall,1 <BR> *get,ssum,node,inuse,s,eqv <BR> *if,smax,ge,ssum,then <BR> smax=smax <BR> *elseif,smax,lt,ssum,then <BR> smax=ssum <BR> *endif <BR>inuse=ndnext(inuse) <BR>*enddo <BR>allsel <BR>*if,thforcal,le,0.012,then <BR> dmq=84 <BR> kbq=58.8 <BR>*elseif,thforcal,gt,0.012,then <BR> dmq=59 <BR> kbq=41.3 <BR>*endif <BR>*if,h,lt,w,then <BR>duanb=h <BR>*elseif,h,ge,w,then <BR>duanb=w <BR>*endif <BR>str=smax/1e6 <BR>esel,s,real,,1,1,1 <BR>/UIS,MSGPOP,1 <BR>*msg,note,str,umax*1000,dmq,10*duanb <BR>Ultimate Limit States Glass STRESS=%gMa,%/& <BR>Serviceabilitylimits DISPLACEMENT=%gmm,%/& <BR>=%gMa,%/& <BR>=%gmm, <BR>/UIS,MSGPOP,3 <BR>*elseif,switch,eq,2,then <BR>nlgeom,on <BR>sstif,on <BR>nsubst,10 <BR>kbc,1 <BR>solv <BR>finish <BR>/post1 <BR>esel,s,real,,1,1,1 <BR>AVPRIN <BR>plnsol,s,eqv,0,1 <BR>allsel <BR>esel,s,real,,1,1,1 <BR>NSLE,S <BR>*get,numall,node,0,count <BR>*get,numstart,node,0,num,min <BR> smax=0 <BR> inuse=numstart <BR>*do,i,1,numall,1 <BR> *get,ssum,node,inuse,s,eqv <BR> *if,smax,ge,ssum,then <BR> smax=smax <BR> *elseif,smax,lt,ssum,then <BR> smax=ssum <BR> *endif <BR>inuse=ndnext(inuse) <BR>*enddo <BR>allsel <BR>*if,thforcal,le,0.012,then <BR> dmq=84 <BR> kbq=58.8 <BR>*elseif,thforcal,gt,0.012,then <BR> dmq=59 <BR> kbq=41.3 <BR>*endif <BR>*if,h,lt,w,then <BR>duanb=h <BR>*elseif,h,ge,w,then <BR>duanb=w <BR>*endif <BR>str=smax/1e6 <BR>esel,s,real,,1,1,1 <BR>/UIS,MSGPOP,1 <BR>*msg,note,str,umax*1000,dmq,10*duanb <BR>Ultimate Limit States Glass STRESS=%gMa,%/& <BR>Serviceabilitylimits DISPLACEMENT=%gmm,%/& <BR>=%gMa,%/& <BR>=%gmm, <BR>/UIS,MSGPOP,3 <BR>*elseif,swtich,eq,0 <BR>allsel <BR>gplot <BR>*endif <BR>/gopr <BR>finish <BR>/encrypt <BR>
页:
[1]