二维桩土分析不知道对错,欢迎批评指正

hongyong78

各位高手:小弟做了份二维的桩土分析,没有考虑接触模型,仅仅考虑了用非线型弹簧去约束土的边界,不知道这样做行不行,欢迎大家给出看法,小弟感激不尽
  这是关于桩的横向受力的分析,请看命令流:
/prep7
!定义单元类型和常数
et,1,plane42
et,2,combin14
keyopt,1,3,2
keyopt,2,3,2
r,1,3.1416*9.6e6/19
!定义材料属性
!桩身
mp,ex,1,8.76E+011   
mp,prxy,1,0.3   
!第一层土
mp,ex,2,9.2e6   
mp,prxy,2,0.3
tb,dp,2
tbdata,1,11000,10,10
!第二层土
mp,ex,3,9.8e6   
mp,prxy,3,0.3   
tb,dp,3
tbdata,1,12000,6.8,6.8
!第三层土
mp,ex,4,10.8e6  
mp,prxy,4,0.3   
tb,dp,4
tbdata,1,16000,10,10
!生成模型
rectng,0,100,0,80
wpoffs,0,74.2,0
wprota,0,90,0   
asbw,1
wpoffs,0,0,6.15
asbw,2
wprota,0,-90,0  
wpoffs,0,-68.05,0
rectng,49.4,50.6,37.05,87.9
aovlap,all  
numcmp,node
numcmp,line
numcmp,area
lplot,all   
APLOT   
LPLOT   
adele,6,,keep   
adele,8,,keep  
adele,7,,keep
adele,1,,keep
ldele,25
ldele,26
ldele,12
al,9,24,18,20,10,4,11,22,15,23  
APLOT   
numcmp,node
numcmp,line
numcmp,area
lcomb,23,17   
lcomb,19,24  
lcomb,25,10
lcomb,11,21
lcomb,27,14
lcomb,28,22

numcmp,node
numcmp,line
numcmp,area
!选择并划分面积
asel,s,,,1
aesize,1,0.3
allsel
asel,s,,,2
asel,a,,,3
asel,a,,,4
asel,a,,,5
aesize,all,1
allsel
asel,s,,,6
aesize,all,2
allsel
!划分网格
mshape,0,2d
mshkey,2
mat,1
amesh,1
mat,2   
APLOT   
asel,s,,,2,3
amesh,all
allsel   
mat,3   
asel,s,,,4,5
amesh,all
allsel
mat,4   
asel,s,,,6  
amesh,all
allsel
!生成节点
n,11916,-2,80  
n,11967,-2,79.71,
ngen,18,1,11967,,,,-0.29
n,11985,-2,74.49
n,11966,-2,74.2
n,12967,-2,73.91
ngen,19,1,12967,,,,-0.29
n,12966,-2,68.05
n,12986,-2,68.34
n,15184,-2,66.11,
ngen,33,-1,15184,,,,-1.94
n,15151,-2,1.944
n,15066,-2,0
n,25066,0,-2
n,25068,2,-2
ngen,49,1,25068,,,2
n,25067,100,-2
n,35067,102,0
n,35117,102,1.944
ngen,34,1,35117,,,,1.944
n,34016,102,68.05
n,34036,102,68.34
ngen,20,-1,34036,,,,0.29
n,30867,102,74.2
n,30886,102,74.49
ngen,19,-1,30886,,,,0.29
n,30866,102,80
!由节点生成单元
type,2
real,1
e,1916,11916
e,1967,11967
*do,i,1,18,1
e,1967+i,11967+i
*enddo
e,1966,11966
e,2967,12967
*do,i,1,19,1
e,2967+i,12967+i
*enddo
e,2966,12966
e,5184,15184
*do,i,0,32,1
e,5151+i,15151+i
*enddo
e,5066,25066
*do,i,0,48,1
e,5068+i,25068+i
*enddo
e,5067,25067
e,5067,35067
e,5117,35117
*do,i,1,33,1
e,5117+i,35117+i
*enddo
e,4016,34016
*do,i,0,19,1
e,4017+i,34017+i
*enddo
e,867,30867
*do,i,0,18,1
e,868+i,30868+i
*enddo
e,866,30866
e,5066,15066
!重新设置单元属性
allsel
esel,s,,,7163,7183
esel,a,,,7366
esel,a,,,7347,7365
esel,a,,,7346
mpchg,2,all
allsel
esel,s,,,7326,7345
esel,a,,,7325
esel,a,,,7184,7204
mpchg,3,all
allsel        
!施加荷载
/SOLU
allsel
nsel,s,loc,x,-2
nsel,a,loc,x,102
nsel,a,loc,y,-2
d,all,all
allsel
nsel,s,,,2
f,2,fx,300000
allsel
time,1  
deltim,0.1
auto,on  
nlgeom,on   
pred,on
lnsrch,on   
outpr,all,all   
outres,all,all  
solve   
!进入后处理观看结果
/post1
set ,last
plnsol,u,x
plnsol,u,sum
plnsol,s,x
plnsol,s,xz

7木

哥们还在不？我也是有这方面问题，不过我的问题简单 ，就是不知道桩与弹簧之间接触建立，能不能指导下啊