ymcheng123 发表于 2012-4-15 23:44

求助:管桩复合地基桩土接触问题

急求帮助!在模型计算的过程中总是出现节点刚度位移过大而导致计算终止,希望各位高手多多指点!问题如下: *** ERROR ***                           CP =   430.922   TIME= 23:20:57
The value of UZ at node 22202 is 56142870.6.It is greater than the   
current limit of 1000000.This generally indicates rigid body motion   
as a result of an unconstrained model.Verify that your model is      
properly constrained.                                                   

*** ERROR ***                           CP =   431.140   TIME= 23:21:11
*** MESSAGE CONTINUATION ---- DIAGNOSTIC INFORMATION ***               
If one or more parts of the model are held together only by contact   
verify that the contact surfaces are closed.You can use the CNCHECK   
command to check the initial contact status in the SOLUTION module.
我的命令流:

fini
/clear
/PLOPTS,date,0
/prep7
/UNITS,SI

!设置单元类型

ET,1,SOLID45

!设置参数

mp,ex,1,1.5e7 $mp,nuxy,1,0.31 $mp,dens,1,1950    !粉质粘土1参数

mp,ex,2,3e7   $mp,nuxy,2,0.28 $mp,dens,2,1900    !细砂参数

mp,ex,3,2.1e7 $mp,nuxy,3,0.35 $mp,dens,3,1950   !粉质粘土2参数

mp,ex,4,8e7   $mp,nuxy,4,0.3 $mp,dens,4,2180      !卵石参数
mp,ex,5,6e7    $mp,nuxy,5,0.15 $mp,dens,5,2040    !路基土参数

mp,mu,6,0.2   $mp,ex,6,3e10$mp,nuxy,6,0.2 $mp,dens,6,3140   !PHC桩参数

mp,ex,7,3e10   $mp,nuxy,7,0.2 $mp,dens,7,2500   !钢筋混凝土板参数

mp,ex,8,8e7    $mp,nuxy,8,0.25 $mp,dens,8,2200   !碎石垫层参数



tb,dp,1
tbdata,1,45e3,18.5,0    !粉质粘土1粘聚力,内摩擦角
tb,dp,2            
tbdata,1,15e3,32,0       !细砂粘聚力,内摩擦角
tb,dp,3            
tbdata,1,40e3,21,0       !粉质粘土2粘聚力,内摩擦角
tb,dp,4            
tbdata,1,10e3,45,0       !卵石粘聚力,内摩擦角
tb,dp,5         
tbdata,1,14.7e3,33,0      !路基土粘聚力,内摩擦角

K,1,0,0,0
K,2,2,0,0
K,3,2,0,-1.5
K,4,0,0,-1.5

A,1,2,3,4

VOFFST,1,-0.5,,      !混凝土板
VOFFST,2,-0.15,,   !碎石垫层
VOFFST,7,-12,,       !土层1
VOFFST,12,-3,,       !土层2
VOFFST,17,-20,,      !土层3
VOFFST,22,-9,,       !土层3
VOFFST,27,-26,,      !土层4


wprota,,-90                     !创建与管桩相同直径的实心桩
wpoffs,,,-0.65
CYL4,1,0.75,0,0,0.2,360,-12
wpoffs,,,-12
CYL4,1,0.75,0,0,0.2,360,-3
wpoffs,,,-3
CYL4,1,0.75,0,0,0.2,360,-20


VSBV,3,8                     ! 删除圆柱体
VSBV,4,9                     ! 删除圆柱体
VSBV,5,10                     ! 删除圆柱体

wpoffs,,,15
CYL4,1,0.75,0.12,0,0.2,360,-35   !创建管桩体





VGLUE,2,11,3,4,6,7    !布尔运算粘结体所有的除桩之外的部分,即土层
NUMCMP,ALL    !压缩编号消除空号
/PNUM,VOLU,1!显示体号码
VPLOT         !显示所选体
SAVE   !

!采用扫掠方法划分网格

VSEL,S,,,4
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,13,0
LSEL,,,,53,56,1
LESIZE,all,,,4
VSEL,,,,4
TYPE,1
MAT,4
VSWEEP,4,29,24    !划分第5层土层
allsel
vplot



VSEL,S,,,6
AADD,36,46      
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,6,0
LSEL,,,,45,48,1
LESIZE,all,,,4
VSEL,,,,6
TYPE,1
MAT,3
VSWEEP,6,24,51    !划分第4层土层
allsel
vplot


VSEL,S,,,8   
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,20,0
LSEL,,,,37,40,1
LESIZE,all,,,2
LSEL,,,,75,78,1
LESIZE,all,,,2

VSEL,,,,8
TYPE,1
MAT,3

mshape,0,3d

mshkey,1


vsweep,all,,                   !划分第3土层
allsel
vplot

VSEL,S,,,2
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,3,0
LSEL,,,,29,32,1
LESIZE,all,,,2
LSEL,,,,65,68,1
LESIZE,all,,,2
VSEL,,,,2
TYPE,1
MAT,2

VSWEEP,2,35,34   !划分第2土层
allsel
vplot


VSEL,S,,,5         !划分垫层
AADD,42,45
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,1,0
LSEL,,,,13,16,1
LESIZE,all,,,4
VSEL,,,,5
TYPE,1
MAT,8
VSWEEP,5,46,2   
allsel
vplot



VSEL,S,,,7      !划分第1土层
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,12,0
LSEL,,,,21,24,1
LESIZE,all,,,2
LSEL,,,,65,68,1
LESIZE,all,,,2
VSEL,,,,7
TYPE,1
MAT,1
VSWEEP,7,34,42   
allsel
vplot


VSEL,S,,,3          !划分桩
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,35,0
LSEL,,,,79,82,1
LESIZE,all,,,2
LSEL,,,,83,86,1
LESIZE,all,,,2
VSEL,,,,3
TYPE,1
MAT,6
VSWEEP,3,15,37   
allsel
vplot


VSEL,S,,,1         !钢筋混凝土板
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,1,0
LSEL,,,,5,8,1
LESIZE,all,,,4
VSEL,,,,1
TYPE,1
MAT,7

VSWEEP,1,2,1   
allsel
vplot

!接触命令流

!allsel
ET,2,170! 定义目标单元类型
ET,3,173! 定义接触单元类型

keyopt,3,4,1
keyopt,3,5,1
keyopt,3,9,1    !消除初始渗透

keyopt,3,12,2    !不分开接触,允许滑动,法向不分开

mp,mu,9,0.2         !接触单元的摩擦系数

allsel,all
r,1,,,30,2,0.00001,,   !定义单元实常数
rmore,,,,-0.00001,,5      !切向接触刚度

! 设置管桩桩侧目标面   

VSEL,S,,,3
ASLV,S
ASEL,S,,,40,41,1
!cm,target,area
aplot
NSLA,S,1
TYPE,2
mat,9
!NSLA,S,1
cm,target,node
ESURF,ALL   
esel,s,type,,2
eplot
/psymb,esys,1   
allsel

! 创建土层2,建立接触面

ALLSEL,ALL
VSEL,S,,,2
ASLV,S
ASEL,S,,,38,39,1
!cm,contact1,area
aplot
NSLA,S,1
TYPE,3
mat,9
!NSLA,S,1
cm,contact1,node
ESURF,ALL
esel,s,type,,3
eplot
/psymb,esys,1
allsel

! 创建土层1接触面

ALLSEL,ALL
VSEL,s,,,7
ASLV,S
ASEL,S,,,47,48,1
!cm,contact2,area
aplot
NSLA,S,1
TYPE,3
mat,9
!NSLA,S,1
ESURF,ALL
cm,contact2,node
esel,s,type,,3
eplot
/psymb,esys,1
allsel

! 创建土层3接触面
allsel,all
VSEL,s,,,8
ASLV,S
ASEL,S,,,49,50,1
!cm,contact3,area
aplot
NSLA,S,1
TYPE,3
mat,9
!NSLA,S,1
cm,contact3,node
ESURF,ALL
esel,s,type,,3
eplot
/psymb,esys,1
allsel



!(桩端与土的接触不知道用不用设置)

!ET,4,170! 定义目标单元类型
!ET,5,173! 定义接触单元类型
!keyopt,5,4,1
!keyopt,5,5,1
!keyopt,5,9,1    !消除初始渗透
!keyopt,5,12,2    !不分开接触,允许滑动,法向不分开

!mp,mu,10,0.2         !接触单元的摩擦系数


!r,2,,,30,1,,,       !定义单元实常数
!rmore,,,,,,5      !切向接触刚度


!桩端,建立目标面


!allsel,all


!VSEL,S,,,3
!ASLV,S
!ASEL,S,,,15,37,22
!cm,target2,area
!TYPE,4
!NSLA,S,1
!ESURF,ALL   
!esel,s,type,,2
!eplot
!/psymb,esys,1   
!allsel               

!桩端(上部),建立接触面

!ALLSEL,ALL
!VSEL,S,,,5
!ASLV,S
!ASEL,S,,,46
!cm,contact4,area
!TYPE,5
!NSLA,S,1
!ESURF,ALL
!esel,s,type,,3
!eplot
!/psymb,esys,1
!allsel

!桩端(下部),建立接触面

!ALLSEL,ALL
!VSEL,S,,,6
!ASLV,S
!ASEL,S,,,51
!cm,contact5,area
!TYPE,5
!NSLA,S,1
!ESURF,ALL
!esel,s,type,,3
!eplot
!/psymb,esys,1
!allsel




vgen,14,1,8,1,2         
vgen,3,1,120,1,,,-1.5

ALLSEL,ALL    !选择所有实体
NUMCMP,ALL    !压缩编号消除空号


k,3001,58,,-4.5       !建立桩周土
k,3002,58,,,
a,2627,3001,3002,834


vext,2101,,,,-0.5
vext,2102,,,,-0.15
vext,2107,,,,-12
vext,2112,,,,-3
vext,2117,,,,-20
vext,2122,,,,-9
vext,2127,,,,-26





!划分体337-343

vsel,s,,,337
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,1,0
LSEL,,,,4285,4287,2
LESIZE,all,,,30
LSEL,,,,4286,4288,2
LESIZE,all,,,12
VSEL,,,,337
TYPE,1
MAT,8
VSWEEP,337,2101,2102   
allsel
vplot


vsel,s,,,338
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,1,0
LSEL,,,,4289,4291,2
LESIZE,all,,,30
LSEL,,,,4290,4292,2
LESIZE,all,,,12
VSEL,,,,338
TYPE,1
MAT,8
VSWEEP,338,2102,2107   
allsel
vplot

vsel,s,,,339
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,12,0
LSEL,,,,4297,4299,2
LESIZE,all,,,30
LSEL,,,,4298,4300,2
LESIZE,all,,,12
VSEL,,,,339
TYPE,1
MAT,1
VSWEEP,339,2107,2112   
allsel
vplot


vsel,s,,,340
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,3,0
LSEL,,,,4305,4307,2
LESIZE,all,,,30
LSEL,,,,4306,4308,2
LESIZE,all,,,12
VSEL,,,,340
TYPE,1
MAT,2
VSWEEP,340,2112,2117   
allsel
vplot


vsel,s,,,341
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,20,0
LSEL,,,,4313,4315,2
LESIZE,all,,,30
LSEL,,,,4314,4316,2
LESIZE,all,,,12
VSEL,,,,341
TYPE,1
MAT,3
VSWEEP,341,2117,2122   
allsel
vplot


vsel,s,,,342
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,6,0
LSEL,,,,4321,4323,2
LESIZE,all,,,30
LSEL,,,,4322,4324,2
LESIZE,all,,,12
VSEL,,,,342
TYPE,1
MAT,3
VSWEEP,342,2122,2127   
allsel
vplot


vsel,s,,,343
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,13,0
LSEL,,,,4329,4331,2
LESIZE,all,,,30
LSEL,,,,4330,4332,2
LESIZE,all,,,12
VSEL,,,,343
TYPE,1
MAT,4
VSWEEP,343,2127,2132   
allsel
vplot

vsel,s,loc,x,26,58
vplot
vgen,2,337,343,1,-58

NUMCMP,ALL               !压缩编号消除空号
allsel


!!!!!生成路基

asel,s,loc,x,0,28
asel,r,loc,y,0
local,11,0,14
csys,11
esize,,4
TYPE,1
mat,5
vext,all,,,,4.5,,0.51785757   
allsel
vplot

!csys,0    !笛卡尔坐标


/solu
asel,s,loc,x,-44      !模型左边界
da,all,uz
da,all,ux

asel,s,loc,x,44         !模型右边界
da,all,uz
da,all,ux

asel,s,loc,z,0            !模型前边界
da,all,uz
asel,s,loc,z,-4.5         !模型后边界
da,all,uz
asel,s,loc,y,-70.65       !模型下边界
da,all,all
allsel

!!!施加重力和预加力

/solu
antype,static
nlgeom,on
nropt,full
pred,on                !打开时间步长预测器
lnsrch,on
outres,all,all
cncheck,detail      !检查接触状态
cncheck,summary      

acel,0,9.8,0,
time,1
solve


ymcheng123 发表于 2012-4-15 23:47

希望多多指点啊,万分感激

ymcheng123 发表于 2012-4-16 12:12

怎么没有人回复啊,求回复!!

ymcheng123 发表于 2012-4-16 13:41

求回复啊,求回复

ymcheng123 发表于 2012-4-16 23:06

怎么没人气啊,求回复

ymcheng123 发表于 2012-4-17 15:13

还是没回复,哎

ymcheng123 发表于 2012-4-17 20:47

回复,回复,回复吧

qusongzhao 发表于 2012-4-18 11:03

你可以尝试在接触面加弹簧单元~~

ymcheng123 发表于 2012-4-18 22:50

回复 8 # qusongzhao 的帖子

好的,我试一下,太感谢了
页: [1]
查看完整版本: 求助:管桩复合地基桩土接触问题