!底板变化为二次抛物线
!三维梁单元beam188单元
/prep7
et,1,beam188
et,2,plane82
mp,ex,1,2e6
mp,prxy,1,0.167
*do,i,1,49
y1=(i**2)/200+5
k,1,5,0
k,2,5,y1
k,3,10,y1
k,4,10,y1+1
k,5,-10,y1+1
k,6,-10,y1
k,7,-5,y1
k,8,-5,0
a,1,2,3,4,5,6,7,8
blc4,-4,1,8,y1-1
asba,1,2
wpoffs,0,1
wprota,0,90
asbw,all
wpoffs,0,,-(y1-1)
asbw,all
wpcsys
wpstyl
mshkey,1
mshape,0,2d
esize,0.5
amap,5,4,3,6,5
amap,1,14,13,1,8
lsel,s,line,,10
lsel,a,line,,12
lsel,a,line,,17
lsel,a,line,,18
lesize,all,,,5
asel,s,area,,3,4
amesh,all
allsel,all
secwrite,%i%,sect,,2
gplot
aclear,all
adele,all
ldele,all
kdele,all
y2=(i+1)**2/200+5
k,1,5,0
k,2,5,y2
k,3,10,y2
k,4,10,y2+1
k,5,-10,y2+1
k,6,-10,y2
k,7,-5,y2
k,8,-5,0
a,1,2,3,4,5,6,7,8
blc4,-4,1,8,y2-1
asba,1,2
wpoffs,0,1
wprota,0,90
asbw,all
wpoffs,0,,-(y2-1)
asbw,all
wpcsys
wpstyl
mshkey,1
mshape,0,2d
esize,0.5
amap,5,4,3,6,5
amap,1,14,13,1,8
lsel,s,line,,10
lsel,a,line,,12
lsel,a,line,,17
lsel,a,line,,18
lesize,all,,,5
allsel,all
asel,s,area,,3,4
amesh,all
allsel,all
t=i+1
secwrite,%t%,sect,,2
gplot
aclear,all
adele,all
ldele,all
kdele,all
*enddo
*do,i,1,50
n,i,0,0,i
*enddo
*do,i,1,49
sectype,i,beam,mesh
secread,%i%,sect,,mesh
secoffset,user,0,(i+1)*(i-1)/200
sectype,i+1,beam,mesh
t=i+1
secread,%t%,sect,,mesh
secoffset,user,0,i*(i+2)/200
sectype,50+i,taper
secdata,i,0,0,i
secdata,i+1,0,0,i+1
type,1
secnum,50+i
mat,1
e,i,i+1
/eshape,1
eplot
*enddo |