以下是本人做的加筋土橋臺三維模型命令流 ，只考慮面板與填土的摩擦接觸，小弟才學習不久，模型有很多問題，總算不下去，望各位前輩指點糾正，希望給些三維加筋橋臺建模的意見！！！ 萬分感激！小弟先行謝過！！！
finish
/clear
/prep7
!!------------定義單元類型與材料性質--------------
et,1,solid65 !面板基礎的單元土體的單元
et,2,solid45
et,3,link10 !筋帶的單元
!et,4,170 !定義面板接觸面
!et,5,173 !定義土體目標面
mp,ex,1,2.4e10 !定義混凝土面板，樁的性質
mp,prxy,1,0.2
mp,dens,1,2400
mp,ex,2,4.0e7 !定義土體性質
mp,prxy,2,0.35
mp,dens,2,1800
tb,dp,2 !定義 土體本構模型
tbdata,,18000,30,0 !定義土體內摩擦角
mp,ex,3,2.0e10 !定義筋帶性質
mp,prxy,3,0.3
mp,dens,3,1500
tbtemp,0 !第二種材料即筋帶用biso模型
tbdata,,200e6,2.1e9
mp,mu,4,0.4 !摩擦系數
r,1
r,2
r,3,0.0006,1.0e-6
!r,4,,,,, ,,
RMORE,,,1000000, ,,
RMORE,,,, ,,
RMORE,,,, ,,
!!r,4,,,100,0.2,0,rmore,,,1e20,0.0,1.0,RMORE,0.0,0,1.0,,1.0,0.5 RMORE,0,0.5,1.0,0.0,
!keyopt,5,5,1
!KEYOPT,5,7,0
!KEYOPT,5,9,0
!KEYOPT,5,12,2
!!第一層筋帶模擬
*do,i,1,13
n,i,i/2-0.25,1,i/2-0.25
*enddo
type,3
mat,3
real,3
*do,i,0,11,1
e,1+i,2+i
*enddo
*do,i,14,26
n,i,i/2-6.75,1,i/2-6.25
*enddo
type,3
mat,3
real,3
*do,i,0,11,1
e,14+i,15+i
*enddo
*do,i,27,39
n,i,i/2-13.25,1,i/2-12.25
*enddo
type,3
mat,3
real,3
*do,i,0,11,1
e,27+i,28+i
*enddo
*do,i,40,52
n,i,i/2-19.75,1,i/2-18.25
*enddo
type,3
mat,3
real,3
*do,i,0,11,1
e,40+i,41+i
*enddo
*do,i,53,65
n,i,i/2-26.25,1,i/2-24.25
*enddo
type,3
mat,3
real,3
*do,i,0,11,1
e,53+i,54+i
*enddo
*do,i,66,77
n,i,i/2-32.25,1,i/2-32.75
*enddo
type,3
mat,3
real,3
*do,i,0,10,1
e,66+i,67+i
*enddo
*do,i,78,88
n,i,i/2-37.75,1,i/2-38.75
*enddo
type,3
mat,3
real,3
*do,i,0,9,1
e,78+i,79+i
*enddo
*do,i,89,98
n,i,i/2-42.75,1,i/2-44.25
*enddo
type,3
mat,3
real,3
*do,i,0,8,1
e,89+i,90+i
*enddo
*do,i,99,107
n,i,i/2-47.25,1,i/2-49.25
*enddo
type,3
mat,3
real,3
*do,i,0,7,1
e,99+i,100+i
*enddo
wpoff,0,0,5.5
CSYS,4
FLST,4,107,1,ORDE,2
FITEM,4,1
FITEM,4,-107
NSYM,Z,107,P51X
CSYS,0
type,3
mat,3
real,3
*do,i,0,11,1
e,160+i,161+i
*enddo
*do,i,0,11,1
e,147+i,148+i
*enddo
*do,i,0,11,1
e,134+i,135+i
*enddo
*do,i,0,11,1
e,121+i,122+i
*enddo
*do,i,0,11,1
e,108+i,109+i
*enddo
*do,i,0,10,1
e,175+i,174+i
*enddo
*do,i,0,9,1
e,185+i,186+i
*enddo
*do,i,0,8,1
e,196+i,197+i
*enddo
*do,i,0,7,1
e,206+i,207+i
*enddo
numcmp,all
EGEN,10,300,1,196,1,,,,,,,0.5,,
*do,i,1,20
k,i,0.25,1,i/2+0.25
*enddo
*do,i,21,40
k,i,8.25,1,i/2-9.75
*enddo
*do,i,1,20
l,i,i+20
*enddo
*do,i,41,61
k,i,i/2-19.75,1,10.75
*enddo
*do,i,62,82
k,i,i/2-30.25,1,3.75
*enddo
*do,i,41,61
l,i,i+21
*enddo
*do,i,83,103
k,i,i/2-40.75,1,0.25
*enddo
*do,i,104,124
k,i,i/2-51.25,1,7.25
*enddo
*do,i,83,103
l,i,i+21
*enddo
cm,lsoil1,line
!!第11層筋帶建立
*do,i,125,139
k,i,i/2-58.75,6,10.75
*enddo
*do,i,140,154
k,i,i/2-66.25,6,3.75
*enddo
*do,i,125,139
l,i,i+15
*enddo
*do,i,155,169
k,i,i/2-73.75,6,0.25
*enddo
*do,i,170,184
k,i,i/2-81.25,6,7.25
*enddo
*do,i,155,169
l,i,i+15
*enddo
cm,lsoil2,line
!!第2_10層筋帶建立
Lgen,10,lsoil1,,0.5,,0.5
!!第12_13層筋帶建立
Lgen,3,lsoil2,,0.5,,0.5
allsel,all
!!劃分筋帶
elemsiz=0.5
lsel,s,loc,y,,1
lsel,a,loc,y,,1.5
latt,3,3,3
lsel,s,loc,y,,1
lsel,a,loc,y,,1.5
mshkey,1
esize,elemsiz
lmesh,all
elemsiz=0.5
lsel,s,loc,y,,2
lsel,a,loc,y,,2.5
latt,3,3,3
lsel,s,loc,y,,2
lsel,a,loc,y,,2.5
mshkey,1
esize,elemsiz
lmesh,all
elemsiz=0.5
lsel,s,loc,y,,3
lsel,a,loc,y,,3.5
latt,3,3,3
lsel,s,loc,y,,3
lsel,a,loc,y,,3.5
mshkey,1
esize,elemsiz
lmesh,all
elemsiz=0.5
lsel,s,loc,y,,4
lsel,a,loc,y,,4.5
latt,3,3,3
lsel,s,loc,y,,4
lsel,a,loc,y,,4.5
mshkey,1
esize,elemsiz
lmesh,all
elemsiz=0.5
lsel,s,loc,y,,5
lsel,a,loc,y,,5.5
latt,3,3,3
lsel,s,loc,y,,5
lsel,a,loc,y,,5.5
mshkey,1
esize,elemsiz
lmesh,all
allsel,all
!!第12_13層筋帶劃分
elemsiz=0.5
lsel,s,loc,y,,6
lsel,a,loc,y,,6.5
latt,3,3,3
lsel,s,loc,y,,6
lsel,a,loc,y,,6.5
mshkey,1
esize,elemsiz
lmesh,all
elemsiz=0.5
lsel,s,loc,y,,7
latt,3,3,3
lsel,s,loc,y,,7
mshkey,1
esize,elemsiz
lmesh,all
allsel,all
!幾何建模
k,600000,,,
k,600001,,5.5
k,600002,3.25,5.5
k,600003,3.25,7.5,0
k,600004,11.25,7.5,0
k,600005,11.25,0,0
a,600000,600001,600002,600003,600004,600005
lsel,none !!劃分混凝土
vsel,none
voffst,1,-11
wpoff,,,-5.25
vsbw,all
wpoff,,,10.5
vsbw,all
wprota,,,90
wpoff,,,0.25
vsbw,all
wpcsys,-1
wprota,,-90,
wpoff,,,5.5
vsbw,all
wpcsys,-1
wprota,,-90,
wpoff,,,0.5
vsbw,all
wpcsys,-1
FLST,5,68,4,ORDE,9
FITEM,5,836
FITEM,5,840
FITEM,5,-844
FITEM,5,846
FITEM,5,-861
FITEM,5,863
FITEM,5,-874
FITEM,5,876
FITEM,5,-909
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1,0.5, , , , , , ,1
vsel,s,,,1,2,1
vsel,a,,,8
vsel,a,,,13,15,2
vsel,a,,,19,20,1
elemsiz=0.5
vatt,1,1,1
mshkey,1
esize,elemsiz
vmesh,all
vsel,s,,,7,9,2
vsel,a,,,4,14,10
vsel,a,,,16,18,1
vsel,a,,,21
elemsiz=0.5
vatt,2,2,2
mshkey,1
esize,elemsiz
vmesh,all
nummrg,all
numcmp,all
lsel,none !!劃分土
vsel,none
voffst,1,-0.5
voffst,3,5
voffst,6,2
vatt,2,2,2
elemsiz=0.5
mshkey,1
esize,elemsiz
vmesh,all
nummrg,all
numcmp,all
allsel,all !接觸對
asel,s,,,10,12,1
asel,a,,,2 !!定義面板目標面
asel,a,,,35
asel,a,,,41
asel,a,,,26,32,6
asel,a,,,15
asel,a,,,22
nsla,s,1
real,4
type,4
mat,4
esurf,all
allsel,all
asel,s,,,48,50,1 !!定義土體接觸面
asel,a,,,53,55,2
asel,a,,,1
nsla,s,1
real,4
type,5
mat,4
esurf,all
allsel,all
nummrg,node
allsel,all
!!---------------定義位移約束-----------------
!!橋臺后側約束
nsel,s,loc,y,0,7.5
nsel,r,loc,x,11.25
d,all,ux
!!橋臺兩側約束
nsel,s,loc,y,0,1
nsel,r,loc,z,0
d,all,uz
nsel,s,loc,y,0,1
nsel,r,loc,z,11
d,all,uz
!!橋臺前側約束
nsel,s,loc,y,0,1
nsel,r,loc,x,0
d,all,ux
!!橋臺蓋梁處約束
nsel,s,loc,y,5.5,7.5
nsel,r,loc,x,3.25
d,all,ux
!!地面約束
nsel,s,loc,x,0,11.25
nsel,r,loc,y,0
d,all,all
allsel,all
!!施加荷載-----------------
ASEL,S,,,49
NSLA,S,1
SF,all,PRES,-36000
!求解控制設置
allsel,all
/solu
antype,static
nlgeom,on
autots,on
lnsrch,1 !打開線性搜索
pred,on !打來預測器
acel,,9.8
time,1
nsubst,50
outres,all,all
neqit,100 !每一子步中方程的迭代次數限值
cnvtol,u,,0.02!定義位移收斂條件
ncnv,2 !如果不收斂時結束而不退出
allsel,all
solve