请教问题:我要做边坡上的预应力锚索,可是我建的模型在加预应力后就不收敛,我把预应力改的很小了,还是不收敛,请教大神帮忙看一下。很急!万分感谢!
finish
/clear
/config,nres,300000
/COM,Structural ! 指定结构分析
/TITLE,bianpo
/FILNAM,bianpo,1 ! 定义工作文件名
/prep7
!****************************各种单元参数的设置***************************************************
H=16
L=24
ET,1,solid45
ET,2,link8
ET,3,combin14
r,1
r,2,0.139*4,5e-6
r,3,15e4,6e4
!定义材料属性
MP,EX,1,1.5E9 !Block1
MP,PRXY,1,0.4
MP,DENS,1,1800
TB,DP,1
TBDATA,1,0.1e6,25
MP,EX,2,1.95E11 !maosuo
MP,PRXY,2,0.3
MP,DENS,2,7921
!****************************创建模型图形***************************************************
k,1,,
k,2,L,
k,3,L,H
k,4,H,H
k,5,L/3,L/3
k,6,,L/3
*do,i,1,5
l,i,i+1
*enddo
l,6,1
al,1,2,3,4,5,6
voffst,1,8
!****************************体划分***************************************************
wprota,,,90
wprota,,90
wpoff,,,-8
vsbw,all
wprota,,90
wpoff,,,-8
vsbw,3
wpoff,,,-6.8
vsbw,4
wpoff,,-8,-8
wprota,,45
vsbw,2
!****************************划分体1***************************************************
!赋予体属性
vsel,all
vatt,1,1,1
asel,s,,,1
lsla,s,1
lesize,all,,,20
allsel
lesize,17,,,20
vsweep,1,1,2
!****************************划分体5
lesize,8,,,20
lesize,35,,,20
lesize,37,,,17
lesize,38,,,17
vsweep,5,22,24
!****************************划分体5
lesize,20,,,17
lesize,4,,,25
lesize,22,,,25
vsweep,4,12,17
!****************************划分体5
lesize,42,,,6
lesize,32,,,6
lesize,19,,,25
lesize,23,,,20
vsweep,6,27,28
!****************************划分体3
lesize,1,,,6
lesize,24,,,20
vsweep,3,15,20
allsel
nummrg,kp
nummrg,node
numcmp,all
!****************************选取z=4m处的节点***************************************************
nsel,s,loc,z,4
TYPE,2
MAT,2
REAL,2
!第一条锚索
e,18444,20784
e,20784,21145
e,21145,21506
e,21506,21867
e,21867,22228
!第二条锚索
e,18311,22912
e,22912,23273
e,23273,23634
e,23634,23995
e,23995,24356
!第三条锚索
e,18178,25040
e,25040,25401
e,25401,25762
e,25762,26123
e,26123,10198
!****************************施加约束***************************************************
asel,s,loc,x
nsla,s,1
d,all,ux
asel,s,loc,y
nsla,s,1
d,all,all
asel,s,loc,x,24
nsla,s,1
d,all,ux
asel,s,loc,z
nsla,s,1
d,all,uz
asel,s,loc,z,8
nsla,s,1
d,all,uz
allsel
finish
!****************************重力下的求解***************************************************
/solu
antype,static !静态分析
nsubst,100 !荷载子步的设置
!pred,on !打开预测
nropt,full !使用完全牛顿拉普森算法。
nlgeom,on !打开大变形
!lnsrch,on !线性搜索
outres,all,all !控制求解的数值写入到数据库:所有
!cnvtol,f,,0.005,2,0.5 !控制非线性的收敛值
acel,,9.8
solve
finish |