源代码如下finish$/clear$/prep7
k,1,0,0$k,2,0,-1$k,3,-0.2,-1.3$k,4,0,-1.5$k,5,0,-2$k,6,-0.2,-2.3$k,7,0,-2.5$k,8,,-3.5$k,9,0.25,-3.5$k,10,0.25,0
*do,i,1,9$l,i,i+1$*enddo
l,10,1
al,all
vrotat,1,,,,,,10,9,90
et,1,solid45$et,2,170$et,3,173
keyopt,3,9,0$keyopt,3,12,2
mp,ex,1,2.6e10$mp,prxy,1,0.167$mp,dens,1,2500
mp,ex,2,5e6$mp,prxy,2,0.42$mp,dens,2,1900
tb,dp,2$tbdata,,19000,31,29
mp,mu,3,0.2$r,1$r,2$r,3
blc4,0.25,,-2,-7,-2$vsbv,2,1
wpoff,-0.3$wprota,,,90$vsbw,all$wpcsys,-1
wpoff,,,-0.55$vsbw,all$wpcsys,-1
wpoff,,-3.6$wprota,,90$vsbw,all$wpcsys,-1$numcmp,all
wpoff,,-1$wprota,,90$vsbw,all
wpoff,,,0.3$vsbw,all
wpoff,,,0.2$vsbw,all
wpoff,,,0.5$vsbw,all
wpoff,,,0.3$vsbw,all
wpoff,,,0.2$vsbw,all
wpoff,,,1$vsbw,all
aadd,8,166
vsel,all$vsel,u,,,15,39,4$vsel,u,,,40
vatt,2,2,1$esize,0.2$mshape,0,3d$mshkey,1$vmesh,all
vsel,all
*do,t,1,7
vatt,2,2,1$vsweep,11+t*4,52+16*t,53+16*t,18+t
*enddo
vatt,2,2,1$vsweep,40,167,155,77
k,100,,0.1$k,101,,-1$k,102,-0.2,-1.3$k,103,,-1.5$k,104,,-2$k,105,-0.2,-2.3$k,106,,-2.5$k,107,,-3.5$k,108,0.25,-3.5$k,109,0.25,0.1
*do,i,100,108$l,i,i+1$*enddo
l,109,100
lsel,s,,,32,33$lsel,a,,,187,188$lsel,a,,,227,228$lsel,a,,,192,194,2$lsel,a,,,109,149,40
al,all
vrotat,8,,,,,,108,109,-90
wpoff,,,-1$vsbw,1
wpoff,,,-0.2$vsbw,5
wpoff,,,-0.3$vsbw,36
wpoff,,,-0.5$vsbw,41
wpoff,,,-0.2$vsbw,42
wpoff,,,-0.3$vsbw,43
vatt,1,1,1$esize,0.2$mshape,0,3d$mshkey,1$vsweep,3,169,170,244
vatt,1,1,1$esize,0.2$mshape,0,3d$mshkey,1$vsweep,1,151,173,243
vatt,1,1,1$esize,0.2$mshape,0,3d$mshkey,1$vsweep,5,172,176,242
vatt,1,1,1$esize,0.2$mshape,0,3d$mshkey,1$vsweep,36,175,179,241
vatt,1,1,1$esize,0.2$mshape,0,3d$mshkey,1$vsweep,41,178,182,240
vatt,1,1,1$esize,0.2$mshape,0,3d$mshkey,1$vsweep,44,186,187,239
vatt,1,1,1$esize,0.2$mshape,0,3d$mshkey,0$vsweep,42,181,185,238
asel,s,,,11,13,2$asel,a,,,33,36,3$asel,a,,,73,88,15$asel,a,,,14$nsla,s,1
real,3$type,2$mat,3$esurf,all$allsel,all
asel,s,,,1,7$nsla,s,1$real,3,type,3$esurf,all
asel,s,loc,z,0$asel,a,loc,x,0.25$da,all,symm
asel,s,loc,x,-1.75$asel,a,loc,z,-2$da,all,all
asel,s,loc,y,-7$da,all,uy
asel,s,,,120$sfa,all,1,pres,50e2$allse,all
/solu$antype,static$nlgeom,on$autots,on
acel,,,9.8$time,1$nsubst,100$outres,all,all
solve$finish
|