|
/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)-kb
w/sin(a2)+kbw*cos(a2)/sin(a2),ky(4)-kbh
a,1,2,4,3
CYL4,kx(5),ky(5),kr
CYL4,kx ,ky ,kr
asba,1,2
asba,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 ,ky
/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 -0.001,kx +0.001
nsel,r,loc,y,ky -0.001,ky +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
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 ,ky ,kr
CYL4,kx(7),ky(7),kr
CYL4,kx ,ky ,kr
asba,1,2
asba,6,3
asba,1,4
asba,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 ,ky
*get,nnn3,node,0,count
egen,2,nnn3+kld,1,kld,1,,,,,,kx(7),ky(7)
|