7.3六角星型

六角星型穹顶的几何非线性分析 finish /clear /prep7 elestif=3 r1=500 r2=250 b=62.16 a=20

*if,elestif,eq,1,then et,1,link8 *else

et,2,beam4 *endif

mp,ex,1,3030.0

mp,prxy,1,3030/1096/2-1 r,1,317

r,2,317,23770,2950,1,1,, rmore,,9180

r,3,317,8370,8370,1,1,, rmore,,14110

csys,1 k,1,r1,30

kgen,6,1,1,,,60 k,7,r2,,b

kgen,6,7,7,,,60 k,13,,,a+b csys,0

*do,i,1,5,1 l,i,i+6 l,i,i+7 *enddo l,6,7 l,6,12

*do,i,7,11,1 l,i,i+1 *enddo l,12,7

*do,i,7,12,1 l,i,13

*enddo

*if,elestif,eq,1,then latt,1,1,1

lesize,all,,,1

*elseif,elestif,eq,2 latt,1,2,2

lesize,all,,,10

*elseif,elestif,eq,3 latt,1,3,2

lesize,all,,,10 *endif lmesh,all

ksel,s,loc,z,0 dk,all,all ksel,all finish /solu antype,0 nlgeom,1

outres,all,all nsubst,200

fk,13,fz,-1500 arclen,1 solve finish

/post26

node1=node(0,0,kz(13)) nsol,2,node1,u,z,dduz prod,3,2,,,,,,-1,1,1 prod,4,1,,,,,,1500,1,1 xvar,3 plvar,4

/post1 pldisp,1

etable,mi,smisc,6 etable,mj,smisc,12 plls,mi,mj,-1

etable,fi,smisc,1 etable,fj,smisc,7 plls,fi,fj,1

pldisp,1

PRETAB,MI,FI