ansys桩的沉降命令流

*dim,mk,array,1,4,,,,
*set,mk(1,1,1),6000
*set,mk(1,2,1),5000
*set,mk(1,3,1),4000
*set,mk(1,4,1),2000
*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,qx,100
*set,bo,1
*set,ho,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

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
r,kk+1,mk(1,kk,1)*abs(vy-8)*bo*ho/6
real,kk+1
e,j+(i-1)*36,1000+j+(i-1)*36
*enddo
*enddo
allsel
!第二组弹簧单元
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)*bo*ho/6
real,kk+10
e,j+(i-1)*36,2000+j+(i-1)*36
*enddo
*enddo
allsel
!第三组弹簧单元
local,11,0,,,,,90
csys,0
allsel
nsel,s,loc,z,1
nsel,r,loc,y,-0.1,8.1
nsel,r,loc,z,-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)*bo*ho/6
real,kk+20
esys,0
e,j+(i-1)*36,3000+j+(i-1)*36
*enddo
*enddo
allsel
!第四组弹簧单元
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)*bo*ho/6
real,kk+20
esys,0
e,j+(i-1)*36,4000+j+(i-1)*36
*enddo
*enddo
allsel
!底层弹簧
allsel
csys,0
ngen,2,5000,1,36,1,,-1
type,2
r,1000,mk(1,1,1)*8*bo*ho/6
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,25,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
es

urf
*enddo
*enddo
!第3组
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
!第4组
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
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
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,press,,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,2
esln,s
esel,r,type,,3
sfe,all,2,press,,FQ(1,1,1)
/psf,pres,TANx,2,0,1
gplot
!第3层
csys,0
allsel
nsel,s,loc,y,2,4
esln,s
esel,r,type,,3
sfe,all,2,press,,FQ(1,2,1)
/psf,pres,TANx,2,0,1
gplot
!第2层
csys,0
allsel
nsel,s,loc,y,4,6
esln,s
esel,r,type,,3
sfe,all,2,press,,FQ(1,3,1)
/psf,pres,TANx,2,0,1
gplot
!第1层
csys,0
allsel
nsel,s,loc,y,6,8
esln,s
esel,r,type,,3
sfe,all,2,press,,FQ(1,4,1)
/psf,pres,tanx,2,0,1
gplot
allsel
gplot
time,10
autot,on
nlgeom,on
nropt,full
csys,0
solve
finish