基坑开挖

dzsky

想做一个深基坑开挖的实例，先用简单的算了下，不知道桩与土的接触对建立的对不对，
算到后面有warning.....
不知道为什么是？

1/4模型，基坑开挖面14m*9m，开挖深度20m；地连墙厚1m，桩深50m，上20米为钢桩(开挖面以上也有桩，为钢桩)，开挖面以下为钢筋混凝土桩，是不是上面的就不需要建立接触对了？

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
附命令流：
/prep7
et,1,plane42   !作为拉伸单元的辅助单元使用
et,2,solid45   !其它单元（土体、地下连续墙、结构等）
et,3,SHELL63   !地下连续墙的单元
et,4,beam188
mp,ex,1,2.5e6   !土层1材料属性         
mp,nuxy,1,0.2            
mp,dens,1,2500
tb,dp,1             !土层1的D-P参数
tbdata,1,19000，31，30
  !定义C30混凝土的材料
mp,ex,2,3e10
mp,nuxy,2,0.2
mp,dens,2,2500
  !定义钢材Q345
mp,ex,3,2.1e11
mp,nuxy,3,0.3
mp,dens,3,7850
!内部网格围轮廓线
  !(1-4内轮廓，5-7外轮廓)
  !点
k,1,0,0,0
k,2,15,0,0
k,3,15,10,0
k,4,0,10,0
　!线
l,1,2
l,2,3
l,3,4
l,4,1
!形成内部线
lgen,2,4,,,3
lgen,2,5,,,1
lgen,2,6,,,6
lgen,2,7,,,1
lgen,2,8,,,3

lgen,2,1,,,,1
lgen,2,10,,,,1
lgen,2,11,,,,1
lgen,2,12,,,,1
lgen,2,13,,,,2
lgen,2,14,,,,1
lgen,2,15,,,,2
alls
lsel,all
lovlap,all

lsel,all
!形成面
!第一排
al,23,24,26,67
al,25,26,30,123
al,29,30,34,124
al,33,34,38,125
al,37,38,41,126
al,21,22,41,68
!第二排
al,63,64,67,74
al,74,80,119,123
al,80,86,120,124
al,86,92,121,125
al,92,98,122,126
al,98,65,66,68
!第三排
al,59,60,64,73
al,73,79,115,119
al,79,85,116,120
al,85,91,117,121
al,91,97,118,122
al,97,61,62,66
!第四排
al,55,72,56,60
al,72,78,111,115
al,78,84,112,116
al,84,90,113,117
al,90,96,114,118
al,96,57,62,58
!第五排
al,51,71,52,56
al,71,77,107,111
al,77,83,108,112
al,83,89,109,113
al,89,95,110,114
al,95,53,54,58
!第六排
al,47,70,48,52
al,70,76,103,107
al,76,82,104,108
al,82,88,105,109
al,88,94,106,110
al,94,49,50,54
!第七排
al,43,69,44,48
al,69,75,99,103
al,75,81,100,104
al,81,87,101,105
al,87,93,102,106
al,93,45,46,50
!第八排
al,18,28,17,44
al,28,32,27,99
al,32,36,31,100
al,36,40,35,101
al,40,42,39,102
al,42,20,19,46

!外围土体定义点、线（并分网）、面
k,,60,0,0
kgen,5,64,,,,1
kgen,2,68,,,,2
kgen,2,69,,,,1
kgen,2,70,,,,2
kgen,2,71,,,,1
kgen,2,72,,,,40

kgen,2,73,,,-45
kgen,2,74,,,-1
kgen,2,75,,,-3
kgen,2,76,,,-1
kgen,2,77,,,-6
kgen,2,78,,,-1
kgen,2,79,,,-3

*do,i,1,16
l,63+i,64+i
*enddo
l,2,64
l,16,65
l,18,66
l,20,67
l,22,68
l,24,69
l,26,70
l,28,71
l,3,72
l,3,74
l,13,75
l,11,76
l,9,77
l,7,78
l,5,79
l,4,80
l,72,73
l,73,74

allsel,all
al,1,20,127,128
al,2,45,128,129
al,3,49,129,130
al,4,53,130,131
al,5,57,131,132
al,6,61,132,133
al,7,65,133,134
al,8,22,134,135
al,9,10,135,136
al,11,21,136,137
al,12,37,137,138
al,13,33,138,139
al,14,29,139,140
al,15,25,140,141
al,16,23,141,142

!形成桩位图
allsel,all
asel,s,area,,14
asel,a,area,,16
asel,a,area,,26
asel,a,area,,28
asel,a,area,,38
asel,a,area,,40
cm,Pile_first,area

　!坑内土体（板面,无地连墙）
asel,s,area,,7,11,1
asel,a,area,,13,17,1
asel,a,area,,19,23,1
asel,a,area,,25,29,1
asel,a,area,,31,35,1
asel,a,area,,37,41,1
asel,a,area,,43,47,1
cm,InSoil_1,area

　!坑外土体面
asel,all
asel,u,area,,InSoil_1
cm,Outsoil_1,area

!平面网格剖分
lsel,all
lsel,s,line,,127,142,1
lsel,a,line,,9,10,1
cm,ten,line
lesize,all,,,10
allsel
lsel,s,line,,14
lsel,a,line,,29
lsel,a,line,,124
lsel,a,line,,120
lsel,a,line,,116
lsel,a,line,,112
lsel,a,line,,108
lsel,a,line,,104
lsel,a,line,,100
lsel,a,line,,31
cm,six,line
lesize,all,,,6
allsel
lsel,s,line,,16
lsel,a,line,,23
lsel,a,line,,67
lsel,a,line,,64
lsel,a,line,,60
lsel,a,line,,56
lsel,a,line,,52
lsel,a,line,,48
lsel,a,line,,44
lsel,a,line,,17
lsel,a,line,,12
lsel,a,line,,37
lsel,a,line,,126
lsel,a,line,,122
lsel,a,line,,118
lsel,a,line,,114
lsel,a,line,,110
lsel,a,line,,106
lsel,a,line,,102
lsel,a,line,,39
cm,three,line
lesize,all,,,3

allsel
lsel,s,line,,5
lsel,a,line,,63
lsel,a,line,,74
lsel,a,line,,80
lsel,a,line,,86
lsel,a,line,,92
lsel,a,line,,98
lsel,a,line,,65
lsel,a,line,,7
lsel,a,line,,55
lsel,a,line,,72
lsel,a,line,,78
lsel,a,line,,84
lsel,a,line,,90
lsel,a,line,,96
lsel,a,line,,57
cm,two,line
lesize,all,,,2
lsel,all
lsel,u,line,,ten
lsel,u,line,,six
lsel,u,line,,three
lsel,u,line,,two
lesize,all,,,1

asel,ALL
aatt,1,,1,0
amesh,ALL




type,2
extopt,esize,1
extopt,aclear,1
extopt,ATTR,0,0,0
mat,1                 !指定单元的材料属性指针
esys,0                !单元坐标系号
asel,s,area,,InSoil_1
vext,InSoil_1,,,0,0,-1
vsel,s,loc,z,0,-1
cm,IN_E1,volu
asel,s,loc,z,0
asel,r,area,,Outsoil_1
vext,Outsoil_1,,,0,0,-1
vsel,s,loc,z,0,-1
vsel,u,volu,,IN_E1
cm,OUT_E1,volu

type,2
extopt,esize,34
extopt,aclear,1
extopt,ATTR,0,0,0
mat,1                 !指定单元的材料属性指针
esys,0                !单元坐标系号
asel,s,LOC,Z,-1
asel,r,LOC,X,0,14
asel,r,LOC,Y,0,9
cm,IN2_A,AREA
vext,IN2_A,,,0,0,-34
vsel,s,loc,z,0,-34
cm,IN_E2,volu
asel,S,LOC,Z,-1
asel,U,AREA,,IN2_A
cm,OUT2_A,AREA
vext,OUT2_A,,,0,0,-34
vsel,s,loc,z,-1,-35
vsel,u,volu,,IN_E2
cm,OUT_E2,volu
type,2
extopt,esize,15
extopt,aclear,1
extopt,ATTR,0,0,0
mat,1                 !指定单元的材料属性指针
esys,0                !单元坐标系号
asel,s,LOC,Z,-35
asel,r,LOC,X,0,14
asel,r,LOC,Y,0,9
cm,IN3_A,AREA
vext,IN3_A,,,0,0,-15
vsel,s,loc,z,-35,-50
cm,IN_E3,volu
asel,S,LOC,Z,-35
asel,U,AREA,,IN3_A
cm,OUT3_A,AREA
vext,OUT3_A,,,0,0,-15
vsel,s,loc,z,-35,-50
vsel,u,volu,,IN_E3
cm,OUT_E3,volu     

type,2
extopt,esize,30
extopt,aclear,1
extopt,ATTR,0,0,0

mat,1                  
esys,0               
asel,S,LOC,Z,-50
vext,ALL,,,0,0,-30
vsel,s,loc,z,-50,-80
cm,E4,volu     

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!
!!删除桩体位置处的土体单元，再建立桩体单元
!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!先选择桩体体位置处的土体单元
allsel,ALL
vsel,S,LOC,Z,-1,-50
vsel,r,LOC,x,3,4
vsel,r,LOC,y,1,2
cm,PILE11_V,VOLU
vsel,S,LOC,Z,-1,-50
vsel,r,LOC,x,3,4
vsel,r,LOC,y,3,4
cm,PILE12_V,VOLU
vsel,S,LOC,Z,-1,-50
vsel,r,LOC,x,3,4
vsel,r,LOC,y,6,7
cm,PILE13_V,VOLU
vsel,S,VOLU,,PILE11_V
vsel,a,VOLU,,PILE12_V
vsel,a,VOLU,,PILE13_V
cm,PILE1_V,VOLU

vsel,ALL
vsel,S,LOC,Z,-1,-50
vsel,r,LOC,x,10,11
vsel,r,LOC,y,1,2
cm,PILE21_V,VOLU
vsel,S,LOC,Z,-1,-50
vsel,r,LOC,x,10,11
vsel,r,LOC,y,3,4
cm,PILE22_V,VOLU
vsel,S,LOC,Z,-1,-50
vsel,r,LOC,x,10,11
vsel,r,LOC,y,6,7
cm,PILE23_V,VOLU
vsel,S,VOLU,,PILE21_V
vsel,a,VOLU,,PILE22_V
vsel,a,VOLU,,PILE23_V
cm,PILE2_V,VOLU
vsel,S,VOLU,,PILE1_V
vsel,a,VOLU,,PILE2_V
cm,Soilpile_V,VOLU

!!!!!!!!!!!!!!!!!!!!!!!!!已选出
vsel,S,VOLU,,Soilpile_V
aslv
asel,U,LOC,Z,-1
asel,U,LOC,Z,-35
cm,Target_Soilp,AREA
vsel,S,VOLU,,Soilpile_V
aslv
asel,r,LOC,Z,-50
cm,Soilptip_A,AREA
vsel,S,VOLU,,Soilpile_V
aslv
asel,r,LOC,Z,-1
cm,Soilptop_A,AREA
!!删除桩体所在处的土体单元
vsel,S,VOLU,,Soilpile_V
vclear,ALL
vdele,ALL
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
!  建立桩体单元
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!建立桩体单元(拉伸PILETOP_A单元)
!分3段拉伸
!第1段 （钢材，拉伸至坑底）
asel,S,AREA,,Soilptop_A
type,4
extopt,esize,20
extopt,aclear,1
extopt,ATTR,0,0,0
mat,3
esys,0
vext,Soilptop_A,,,0,0,-19
cm,Pile_Gang,volu

vsel,s,volu,,Pile_Gang
aslv
asel,U,LOC,Z,-1
asel,U,LOC,Z,-20
cm,Contact_Pile1,AREA
vsel,s,volu,,Pile_Gang
aslv
asel,r,LOC,Z,-20
cm,Pilesecond_A,AREA

!第2段 （混凝土，拉伸至地下墙底面处）
asel,S,AREA,,Pilesecond_A
type,4
extopt,esize,14
extopt,aclear,1
extopt,ATTR,0,0,0
mat,4
esys,0
vext,Pilesecond_A,,,0,0,-15
cm,Pile_Hun1,volu
vsel,s,volu,,Pile_Hun1
aslv
asel,r,LOC,Z,-20,-35
asel,U,LOC,Z,-20
asel,U,LOC,Z,-35
cm,Contact_Pile2,AREA
vsel,s,volu,,Pile_Hun1
aslv
asel,r,LOC,Z,-35
cm,Pilethird_A,AREA
!第三段，拉伸至桩底
asel,S,AREA,,Pilethird_A
type,4
extopt,esize,10
extopt,aclear,1
extopt,ATTR,0,0,0
mat,4
esys,0
vext,Pilethird_A,,,0,0,-15  
cm,Pile_Hun2,volu
vsel,s,volu,,Pile_Hun2
aslv
asel,r,LOC,Z,-35,-50
asel,U,LOC,Z,-35
asel,U,LOC,Z,-50
cm,Contact_Pile3,AREA
vsel,s,volu,,Pile_Hun2
aslv
asel,r,LOC,Z,-50
cm,Pilebottom_A,AREA
allsel,ALL
vsel,s,volu,,Pile_Gang
vsel,a,volu,,Pile_Hun1
vsel,a,volu,,Pile_Hun2
cm,Pile,volu
!!!综合以上三段桩体单元的接触面，得到contact_pile的接触面
asel,S,AREA,,Contact_Pile1
asel,a,AREA,,Contact_Pile2
asel,a,AREA,,Contact_Pile3
cm,Contact_Pile,AREA

et,5,170
et,6,174
mp,MU,8,0.6

r,12,,,10,0.1,,                  !!!!!桩体与土体的接触用R参数号为12
mat,8
real,13
keyopt,6,12,2

asel,s,area,,Target_Soilp
type,5
real,1
mat,8
nsla,s,1
esurf,all

asel,s,area,,Contact_Pile
type,6
real,1
mat,8
nsla,s,1
esurf,all
!桩尖
asel,s,area,,Soilptip_A
type,5
real,1
mat,8
nsla,s,1
esurf,all

asel,s,area,,Pilebottom_A
type,6
real,1
mat,8
nsla,s,1
esurf,all

!#######################################
!调整部分土体模量（调整实参数号）
!##########################################
!!!!!!!墙体
vsel,S,LOC,Z,0,-35
vsel,r,loc,x,0,14
vsel,r,LOC,Y,0,9
cm,VV,VOLU
vsel,S,LOC,Z,0,-35
vsel,r,LOC,X,0,15
vsel,r,LOC,Y,0,10
vsel,U,VOLU,,VV
cm,soilwall,volu
eslv
emodif,ALL,mat,3


!!!!!!!!!!!!!!!!!!!!!!!
!!                   !!
!!    加载和求解     !!
!!                   !!
!!!!!!!!!!!!!!!!!!!!!!!
!施加边界条件及重力加速度
/solu
!!!!对土体进行约束
allsel
asel,S,LOC,Z,-80
da,all,all,0
asel,s,LOC,X,0
da,all,symm
asel,s,LOC,X,60
da,all,ux,0
asel,s,LOC,Y,0
da,all,symm
asel,s,LOC,Y,50
da,all,uy,0

deltim,0.1,0.05,0.2
autots,on        !使用自动时间步
pred,on
lnsrch,on
nlgeom,on
nropt,full
cnvtol,f,,0.05,,0.5
!进行初始地应力计算
vsel,s,volu,,soilwall
vsel,s,volu,,Pile
eslv
ekill,all
esel,all
esel,s,live
nsel,invert
d,all,all
nsel,all
esel,all
/pbc,all,,1
gplot
acel,,,-9.8
solve

[ 本帖最后由 dzsky 于 2007-11-10 19:35 编辑 ]

ljz0702

桩多的话，用接触很难收敛，就是收敛，所得到的摩阻力也不太正确