马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
急求帮助!在模型计算的过程中总是出现节点刚度位移过大而导致计算终止,希望各位高手多多指点!问题如下: *** 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
|