三维加筋土桥台

mxzpig369

以下是本人做的加筋土桥台三维模型命令流 ，只考虑面板与填土的摩擦接触，小弟才学习不久，模型有很多问题，总算不下去，望各位前辈指点纠正，希望给些三维加筋桥台建模的意见！！！  万分感激！小弟先行谢过！！！
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