finish
/clear
!ز
*DIM,MK,ARRAY,1,4,, , ,     !ˮƽϵmľ
*SET,MK(1,1,1) , 6000       !Ĳmֵ
*SET,MK(1,2,1) , 5000       !mֵ
*SET,MK(1,3,1) , 4000       !ڶmֵ
*SET,MK(1,4,1) , 2000       !һmֵ
*DIM,FQ,ARRAY,1,4,, , ,     !׮ܵλĦ
*SET,FQ(1,1,1) , 40.53      !ԵµĵλĦ
*SET,FQ(1,2,1) , -33.78     !ԵϵĵλĦ
*SET,FQ(1,3,1) , -16.14  
*SET,FQ(1,4,1) , -4.8
*SET,QX, 100                !ˮƽ
*SET,QZ, 100
*SET,b0, 1
*SET,h0, 0.5


!ǰ
/prep7
et,1,solid45                !׮Ԫ
et,2,combin14               !ɵԪ
et,3,surf154                !ЧӦԪ
keyopt,1,1,0
keyopt,2,3,0                !Ϊά
keyopt,3,4,1                !мڵ
r,1
mp,ex,1,3.25e10             !׮Ĳ
mp,nuxy,1,0.3
mp,dens,1,2700

!׮ģ͵Ľڵ
n,1                         
n,6,1
fill,1,6
ngen,6,6,1,6,1,,,0.2
ngen,21,36,1,36,1,,0.5
/view,1,1,1,1               !Ȳʾ
/ang,1
/REP,FAST
GPLOT

!׮ļʵģ
k,1
k,2,1
k,3,1,,1
k,4,,,1
k,5,,10,
k,6,1,10,
k,7,1,10,1
k,8,,10,1
/PNUM,kp,1
kplot
v,1,2,3,4,5,6,7,8
lsel,s,loc,y,0
lsel,a,loc,y,10
lesize,all,,,5
allsel
lsel,s,loc,y,5
lesize,all,,,20
vmesh,1
nummrg,all                  !ѹظԪ
numcmp,all                  !ϲظԪ

!һ鵯ģ,7*16=102
nsel,s,loc,x,1
nsel,r,loc,y,-0.1,8.1
nplot
ngen,2,1000,all,,,2
type,2
*do,i,1,17,1
 *if,i,le,4,then            !жһ
   kk=1
   *elseif,i,le,8,then
   kk=2
   *elseif,i,le,12,then
   kk=3
   *elseif,i,le,17,then
   kk=4
 *endif
 *do,j,6,36,6
 *get,vy,node,j+(i-1)*36,loc,y       !ȡڵYֵ
    r,kk+1,mk(1,kk,1)*abs(vy-8)*b0*h0/6      !Yֵˮƽϵɵĸն
    real,kk+1
    e,j+(i-1)*36,1000+j+(i-1)*36
 *enddo
*enddo
allsel
!Ԫ602

!ڶ鵯ɵԪ,7*16=102
allsel
nsel,s,loc,x,0
nsel,r,loc,y,-0.1,8.1
nsel,r,loc,z,-0.1,1.1
nplot
ngen,2,2000,all,,,-2
type,2
*do,i,1,17,1
 *if,i,le,4,then
   kk=1
   *elseif,i,le,8,then
   kk=2
   *elseif,i,le,12,then
   kk=3
   *elseif,i,le,17,then
   kk=4
 *endif
 *do,j,1,31,6
 *get,vy,node,j+(i-1)*36,loc,y 
    r,kk+10,mk(1,kk,1)*abs(vy-8)*b0*h0/6
    real,kk+10
    e,j+(i-1)*36,2000+j+(i-1)*36
 *enddo
*enddo
allsel
!Ԫ704

!鵯ɵԪ,7*16=102
local,11,0,,,,,90
csys,0
allsel
nsel,s,loc,z,1
nsel,r,loc,y,-0.1,8.1
nsel,r,loc,x,-0.1,1.1
nplot
ngen,2,3000,all,,,,,2
type,2
*do,i,1,17,1
 *if,i,le,4,then
   kk=1
   *elseif,i,le,8,then
   kk=2
   *elseif,i,le,12,then
   kk=3
   *elseif,i,le,17,then
   kk=4
 *endif
 *do,j,31,36,1
 *get,vy,node,j+(i-1)*36,loc,y 
    r,kk+20,mk(1,kk,1)*abs(vy-8)*b0*h0/6
    real,kk+20
    esys,0
    e,j+(i-1)*36,3000+j+(i-1)*36
 *enddo
*enddo
allsel
!Ԫ806

!鵯ɵԪ,7*16=102
allsel
nsel,s,loc,z,0
nsel,r,loc,y,-0.1,8.1
nsel,r,loc,x,-0.1,1.1
nplot
ngen,2,4000,all,,,,,-2
type,2
*do,i,1,17,1
 *if,i,le,4,then
   kk=1
   *elseif,i,le,8,then
   kk=2
   *elseif,i,le,12,then
   kk=3
   *elseif,i,le,17,then
   kk=4
 *endif
 *do,j,1,6,1
 *get,vy,node,j+(i-1)*36,loc,y 
    r,kk+20,mk(1,kk,1)*abs(vy-8)*b0*h0/6
    real,kk+20
    esys,0
    e,j+(i-1)*36,4000+j+(i-1)*36
 *enddo
*enddo
allsel
!ԪΪ908

!ײ㵯
allsel
csys,0
ngen,2,5000,1,36,1,,-1
type,2
r,1000,mk(1,1,1)*8*b0*h0/6*20
real,1000
*do,i,1,36,1
  e,i,5000+i
*enddo
allsel
gplot

!һЧӦԪ
allsel
r,1001
real,1001
local,13,0,,,,90
type,3
*do,i,1,16,1
 *do,j,6,30,6
  nsel,s,node,,j+(i-1)*36
  nsel,a,node,,j+6+(i-1)*36
  nsel,a,node,,j+42+(i-1)*36
  nsel,a,node,,j+36+(i-1)*36
  esys,13
  esurf
 *enddo
*enddo

!ڶЧӦԪ
allsel
local,13,0,,,,90
type,3
*do,i,1,16,1
 *do,j,1,29,6
  nsel,s,node,,j+(i-1)*36
  nsel,a,node,,j+6+(i-1)*36
  nsel,a,node,,j+42+(i-1)*36
  nsel,a,node,,j+36+(i-1)*36
  esys,13
  esurf
 *enddo
*enddo

!ЧӦԪ
allsel
local,14,0,,,,90
type,3
*do,i,1,16,1
 *do,j,1,5,1
  nsel,s,node,,j+(i-1)*36
  nsel,a,node,,j+1+(i-1)*36
  nsel,a,node,,j+37+(i-1)*36
  nsel,a,node,,j+36+(i-1)*36
  esys,14
  esurf
 *enddo
*enddo
csys,0

!ЧӦԪ
allsel
local,15,0,,,,90
type,3
*do,i,1,16,1
 *do,j,31,35,1
  nsel,s,node,,j+(i-1)*36
  nsel,a,node,,j+1+(i-1)*36
  nsel,a,node,,j+37+(i-1)*36
  nsel,a,node,,j+36+(i-1)*36
  esys,15
  esurf
 *enddo
*enddo
save
finish

!
/solu
csys,0
allsel
antype,static
alls
fcum,add,,
nsel,s,loc,x,3                
nsel,a,loc,x,-2
nsel,a,loc,z,3
nsel,a,loc,z,-2
d,all,all                   !Լ߽
nsel,a,loc,y,-1
d,all,all
nsel,a,loc,y,0
d,all,ux
d,all,uz
ACEL,0,9.81,0,              !ʩٶ

allsel
TIME,1
deltim,0.2,0.05,0.4  
autots,on !ʹԶʱ䲽 
pred,on !ʱ䲽Ԥ 
lnsrch,on ! 
nlgeom,on !򿪴λЧ 
nropt,full !趨ţ-ɭѡ 
cutcontrol,plslimit,0.20
KBC,0               !Ӻ
cnvtol,f,,0.02,2,1 
!cnvtol,f,,0.7,2
!NSUBST,30
OUTRES,basic,1   
solve 
rescontrol,file_summany


/solu
antype,,rest,1,4
!׮ƽֲغ
allsel
nsel,s,loc,y,10
nsel,r,loc,x,-0.1,1.1
nsel,r,loc,z,-0.1,1.1
nplot
esln,s
sfe,all,6,pres,,5000        !׮ֲ
nsel,r,loc,x,1
f,all,fx,QX                  !ˮƽغ
nsel,s,loc,y,10
nsel,r,loc,z,1
f,all,fz,QZ                  !ˮƽغ
allsel
gplot

!ĲıĦ
csys,0
allsel
nsel,s,loc,y,0,1.5
esln,s
esel,r,type,,3
sfe,all,2,pres,,FQ(1,1,1)
/PSF,PRES,TANx,2,0,1
gplot

!ı渺Ħ
csys,0
allsel
nsel,s,loc,y,2.5,3.5
esln,s
esel,r,type,,3
sfe,all,2,pres,,FQ(1,2,1)
/PSF,PRES,TANx,2,0,1
gplot

!ڶı渺Ħ
csys,0
allsel
nsel,s,loc,y,4.5,5.5
esln,s
esel,r,type,,3
sfe,all,2,pres,,FQ(1,3,1)
/PSF,PRES,TANx,2,0,1
gplot

!һı渺Ħ
csys,0
allsel
nsel,s,loc,y,6.5,7.5
esln,s
esel,r,type,,3
sfe,all,2,pres,,FQ(1,4,1)
/PSF,PRES,TANx,2,0,1
gplot

allsel
gplot
allsel,all
TIME,2
deltim,0.2,0.05,0.4  
autots,on !ʹԶʱ䲽 
pred,on !ʱ䲽Ԥ 
lnsrch,on ! 
nlgeom,on !򿪴λЧ 
nropt,full !趨ţ-ɭѡ 
cutcontrol,plslimit,0.20
KBC,0               !Ӻ
cnvtol,f,,0.02,2,1 
!cnvtol,f,,0.7,2
!NSUBST,30
OUTRES,basic,1   
solve 
save
rescontrol,file_summany

/POST1
esel ,s,type,,1
PLNSOL,S,y,1,1

/POST1
LCDEF,1,1, ,          !1㲽Ϊ1
LCDEF,2,2, ,          !6㲽(ͶΪ100%)Ϊع2
LCASE,2,             !ع2
LCOPER,SUB,1, , ,     !ع2ȥع1
     plnsol,u,z
