finish
/clear
/prep7

!ģ
k,1,0,2.7388794
k,2,-2.0711036,1.7033277
k,3,-2.6388794,0
k,4,-2.3359835,-1.9455903
k,5,-1.6736199,-2.6205221
k,6,0,-2.8972106
k,1000,,,
larc,1,2,1000,2.589
larc,2,3,1000,2.839
larc,3,4,1000,6.4
larc,4,5,1000,1.05
larc,5,6,1000,5.2
lsymm,x,all
nummrg,all
numcmp,all
lcomb,1,2
lcomb,6,7
nummrg,all
numcmp,all
lcomb,3,4
lcomb,7,8
nummrg,all
numcmp,all
lcomb,2,3
lcomb,5,6
nummrg,all
numcmp,all

al,all
blc4,-36,-69,72,84
asba,2,1,,,keep
numcmp,all

!з
asel,s,,,2 
wprota,,,90
asbw,all
wpoffs,,,6           !Zƶ6
asbw,all
wpoffs,,,-12         !ZḺƶ6
asbw,all
wpcsys,-1            !ָʼ״̬
wprota,,90
asbw,all
wpoffs,,,6
asbw,all
wpoffs,,,-12
asbw,all
wpcsys,-1
alls
nummrg,all
numcmp,all

!

*set,th1,0.3
et,1,beam3
et,2,plane42
keyopt,2,3,2

r,1,th1,th1*th1*th1/12,th1,,

mp,ex,1,31e9       !!Ĳ
mp,prxy,1,0.3
mp,dens,1,2500

mp,ex,2,1.40e8     !!ز
mp,prxy,2,0.3
mp,dens,2,2049
tb,dp,2
tbdata,1,3.4e4,30

!
lsel,s,,,1,4
latt,1,1,1
asel,all
aatt,2,,2

!Χ
lsel,s,,,1,3,2
lesize,all,,,8
lsel,s,,,2,4,2
lesize,all,,,8
lsel,s,,,1,4
lmesh,all

lsel,s,,,24,25
lsel,a,,,35
lsel,a,,,44
lesize,all,,,8,3
lsel,s,,,29,30
lesize,all,,,4
lsel,a,,,32,34,2
lesize,all,,,4
lsel,s,,,38,39
lsel,a,,,41,43,2
lesize,all,,,4

asel,s,,,1
amesh,all
alls
amap,16,1,29,19,2
amap,17,4,17,29,1
amap,12,3,24,19,2
amap,13,4,17,24,3

!
lsel,s,,,15
lsel,a,,,18
lesize,all,,,4
lsel,a,,,16
lsel,a,,,19
lesize,all,,,4
lsel,s,,,12
lsel,a,,,13
lsel,a,,,23
lsel,a,,,27
lsel,a,,,36
lesize,all,,,20
lsel,s,,,5
lsel,a,,,7
lsel,a,,,11
lsel,a,,,14
lsel,a,,,20
lesize,all,,,20

!
lsel,s,,,9,10
lsel,a,,,37
lsel,a,,,21,22
lesize,all,,,6

lsel,s,,,42
lsel,a,,,40
lesize,all,,,4

lsel,s,,,33
lsel,a,,,31
lesize,all,,,4

lsel,s,,,17
lsel,a,,,26
lsel,a,,,28
lsel,a,,,8
lsel,a,,,6
lesize,all,,,42

alls
asel,u,,,1
asel,u,,,12,13
asel,u,,,16,17
amesh,all
alls

acel,,9.8

asel,s,,,1       !
esla,s
ekill,all
alls

nsel,s,loc,x,-36
d,all,uy,0
nsel,s,loc,x,36
d,all,uy,0
nsel,s,loc,y,-69
d,all,all
alls

!ٶʱ
alls
*DIM,shuiping,ARRAY,2,1000,0, , ,   
!*  
*CREATE,ansuitmp
*VREAD,shuiping(1,1),'shuiping','txt','E:\',1000, , , ,
(e9.3,e11.3)
*END
/INPUT,ansuitmp

!!!!!
/solu
antype,4
trnopt,full 
lumpm,0
!autots,on     !ʹԶʱ䲽 
!pred,on       !ʱ䲽Ԥ 
!lnsrch,on     ! 
!nlgeom,on     !򿪴λЧ 
!nropt,full    !趨ţ-ɭѡ 
!cnvtol,f,,0.02,2,0.5  !
!CUTCONTROL,DSPLIMIT,1000000000
!NCNV,2,        !㲻ʱ˳ 
!NEQIT,40             !ÿһӲз̵ĵ(Ĭ25)
!SOLCONTROL,OFF, ,    !  λ

*do,t,1,75,1
time,0.02*t
kbc,1
nsub,1
acel,shuiping(2,t)
alphad,0.069  !ԽӰ첻
betad,0.035
solve
*enddo

*do,t,76,150,1
time,0.02*t
kbc,1
nsub,50
acel,shuiping(2,t)
alphad,0.069  !ԽӰ첻
betad,0.035
solve
*enddo
*do,t,151,750,1
time,0.02*t
kbc,1
nsub,600
acel,shuiping(2,t)
alphad,0.069  !ԽӰ첻
betad,0.035
solve
*enddo
save
!finish




!
/post26
/axlab,y,uy          !!yΪuy
nsol,2,2811,u,y,uy_2   !ʾʱ̱仯plvar,2 
prvar,2               ı
plvar,2              ͼ

/axlab,y,vy
deriv,3,2,1,,,,,1
plvar,3
prvar,3
/axlab,y,ay
deriv,4,3,1,,,,,1
plvar,4
prvar,4
plvar,2,3,4
eep,all
            


!ˮƽ
/post26
/axlab,y,ux
nsol,2,2811,u,x,ux_2
plvar,2                 !ͼ
prvar,2                 !
/axlab,y,vx
deriv,3,2,1,,,,,1
plvar,3
prvar,3
/axlab,y,ax
deriv,4,3,1,,,,,1
plvar,4
prvar,4
plvar,2,3,4
!eep,all
