24小时热门版块排行榜

返回列表

zdy_derive

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 9
帖子: 6
在线: 1小时
虫号: 4245834
注册: 2015-11-25
专业: 岩土与基础工程

[求助] 跪求大神看看这个命令流哪里有问题，好人一生平安！！！！已有1人参与

命令流如下：
finish
/clear
/PREP7
!*
ET,1,SOLID45
ET,2,SOLID65
et,3,SOLID65
!*
MP,EX,1,4800
MP,EY,1,4800
MP,EZ,1,51500
MP,PRXY,1,0.02
MP,PRYZ,1,0.02
MP,PRXZ,1,0.27
MP,GXY,1,2353
MP,GYZ,1,2353
MP,GXZ,1,20275
mp,ex,2,20000
mp,nuxy,2,0.25
tb,concr,2
tbdata,,,,2,20
mp,ex,3,1.5e4
mp,nuxy,3,0.2
tb,dp,3
tbdata,1,4,30,10
r,1$r,2$r,3
!
CYLIND,16,0,-300,1500,0,90,
CYLIND,50,16,0,1500,0,90,
CYLIND,500,50,0,1500,0,90,
!
allsel,all
vsel,s,,,1
aslv,s
lsla,s
lplot
lesize,1,,,4
lesize,2,,,4
lesize,3,,,6
lesize,4,,,6
lesize,5,,,4
lesize,6,,,4
lesize,7,25
lesize,8,25
lesize,9,25
type,1
mat,1
real,1
mshape,0,3d
mshkey,0
vsweep,all
!*
!
allsel,all
vsel,s,,,2
aslv,s
lsla,s
lplot
lesize,10,,,4
lesize,11,,,6
lesize,12,,,4
lesize,13,,,6
lesize,14,,,6
lesize,15,,,4
lesize,16,,,6
lesize,17,,,4
lesize,18,25
lesize,19,25
lesize,20,25
lesize,21,25
type,2
mat,2
real,2
mshape,0,3d
mshkey,0
vsweep,all
!*
!
allsel,all
vsel,s,,,3
aslv,s
lsla,s
lplot
lesize,23,,,6
lesize,25,,,6
lesize,26,,,6
lesize,28,,,6
lesize,22,50
lesize,24,50
lesize,27,50
lesize,29,50
lesize,30,50
lesize,31,50
lesize,32,50
lesize,33,50
type,3
mat,3
real,3
mshape,0,3d
mshkey,0
vsweep,all
!*
allsel,all
!
save
!
!
mp,mu,4,0.4 ! 修改
mat,4
r,4,4
real,4
et,4,170
et,5,174
KEYOPT,5,9,0
KEYOPT,5,10,2
nropt,unsym    ! 不对称矩阵完全牛顿—拉普森法
!*generate the target sunface
allsel,all
asel,s,,,9
cm,target1,area
type,4
nsla,s,1
esln,s,0
esurf
!*generate the contact sunface
asel,s,,,3
cm,contact1,area
type,5
nsla,s,1
esln,s,0
esurf
!
allsel,all
eplot
!
!*CONTACT PAIR CREATION - END
mp,mu,5,0.7    ! 修改
mat,5
r,5,4
real,5
et,6,170
et,7,174
KEYOPT,7,9,0
KEYOPT,7,10,2
nropt,unsym
!*generate the target sunface
allsel,all
asel,s,,,15
cm,target2,area
type,6
nsla,s,1
esln,s,0
esurf
!*generate the contact sunface
asel,s,,,8
cm,contact2,area
type,7
nsla,s,1
esln,s,0
esurf
!
allsel,all
eplot
!*CONTACT PAIR CREATION - END
save
!
!*
/solu
da,4,symm !对称约束
da,5,symm
da,10,symm
da,11,symm
da,16,symm
da,17,symm
da,13,all
da,14,all
da,7,all
!
antype,0
nlgeom,1       !大变形
autots,on  !自动时间步
neqit,50 !子步最大迭代次数
pred,on !非线性分析预测器（加速收敛）
!
solcontrol,on,on
!
nropt,full,,on    !定义H法选项（270）
nsubst,20,10000,2 !子步数
time,100          !时间
allsel,all
solve
!*
nsubst,100,10000,10
outres,all,all
autots,1
time,200
asel,s,loc,z,-300
sfa,all,,pres,-10       !修改荷载值
allsel,all
solve

错误提示为：
*** ERROR ***                         CP =    51.231 TIME= 12:26:52
The value of UZ at node 67 is 1.342765716E+09.  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 =    51.387 TIME= 12:30:06
*** 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 check contact
status in the SOLUTION module for the converged solutions using
CNCHECK.

*** ERROR ***                         CP =    51.480 TIME= 12:30:22
*** MESSAGE CONTINUATION ---- DIAGNOSTIC INFORMATION ***
Rigid body motion can also occur when net section yielding has
occurred resulting in large displacements for small increments of load
or when buckling has occurred.  You can plot the time history curve
for node 67 in the UZ direction to check for stiffness (slope of the
curve) approaching zero.

回复此楼