24小时热门版块排行榜

返回列表

kakaqun

新虫 (小有名气)

应助: 4 (幼儿园)
金币: 902.9
帖子: 57
在线: 16.7小时
虫号: 3496291
注册: 2014-10-24
专业: 岩土与基础工程

[求助] 隧道开挖模拟，不收敛，请帮忙看一下原因已有1人参与

命令流如下，帮帮忙吧~~很急很关键！
finish
/clear
/prep7
k,1, 5.500000,3.5
k,2,-5.500000,3.5
k,3, 0.000000,4
k,5,-5.500000,-4.5
k,6, 5.500000,-4.5
k,8,-4.500000,-4.5
k,9,-4.500000,-8.5
k,10,-3.500000,-4.5
k,11,-3.500000,-8.5
k,12,-2.500000,-4.5
k,13,-2.500000,-8.5
k,14,-1.500000,-4.5
k,15,-1.500000,-8.5
k,16,-0.500000,-4.5
k,17,-0.500000,-8.5
k,18, 0.500000,-4.5
k,19, 0.500000,-8.5
k,20, 1.500000,-4.5
k,21, 1.500000,-8.5
k,22, 2.500000,-4.5
k,23, 2.500000,-8.5
k,24, 3.500000,-4.5
k,25, 3.500000,-8.5
k,26, 4.500000,-4.5
k,27, 4.500000,-8.5
k,100,-5.500000,-8.5
k,101,5.500000,-8.5
k,28, 15.00000, 15.00000
k,29,-15.00000, 15.00000
k,31,-15.00000,-15.00000
k,33, 15.00000,-15.00000
k,37, 45.00000, 15.00000
k,39, 45.00000,-15.00000
k,43, 15.00000, 45.00000
k,45, 45.00000, 45.00000
k,49,-15.00000, 45.00000
k,53,-45.00000, 45.00000
k,55,-45.00000, 15.00000
k,59,-45.00000,-15.00000
k,63,-45.00000,-45.00000
k,65,-15.00000,-45.00000
k,69,15.00000,-45.00000
k,73,45.00000,-45.00000
larc,1,2,3,5.5
l,2,5
l,5,8
l,8,10
l,10,12
l,12,14
l,14,16
l,16,18
l,18,20
l,20,22
l,22,24
l,24,26
l,26,6
l,1,6
l,8,9
l,10,11
l,12,13
l,14,15
l,16,17
l,18,19
l,20,21
l,22,23
l,24,25
l,26,27
l,31,33
l,33,28
l,28,29
l,29,31
l,59,31
l,31,65
l,63,65
l,59,63
l,65,69
l,33,69
l,33,39
l,39,73
l,69,73
l,28,37
l,37,39
l,43,28
l,43,45
l,45,37
l,49,43
l,49,29
l,53,49
l,53,55
l,55,29
l,55,59
l,5,100
l,6,101
l,100,9
l,9,11
l,11,13
l,13,15
l,15,17
l,17,19
l,19,21
l,21,23
l,23,25
l,25,27
l,27,101
lsel,s,line,,1,14,1
al,all
lsel,s,line,,1,14,1
lsel,a,line,,25,28,1
al,all
allsel,all,all
al,29,30,31,32
al,25,30,33,34
al,34,35,36,37
al,26,35,39,38
al,38,40,41,42
al,27,40,43,44
al,44,45,46,47
al,28,29,47,48
!围岩
E_1=1.4e15
U_1=0.4
DENS_1=2620
C_1=1.6e5
FI_1=28
MP,EX,1,E_1
MP,PRXY,1,U_1
MP,DENS,1,DENS_1
TB,DP,1
TBDATA,1,C_1,FI_1

!锚杆
E_2=17E20
U_2=0.3
DENS_2=7900
MP,EX,2,E_2
MP,PRXY,2,U_2
MP,DENS,2,DENS_2

!衬砌
E_3=21e15
U_3=0.2
DENS_3=2900
MP,EX,3,E_3
MP,PRXY,3,U_3
MP,DENS,3,DENS_3

ET,1,42
KEYOPT,1,3,2
et,2,3
keyopt,2,6,1
ET,3,1
r,1,0.3,0.3/12,0.3
r,2,3.14*(25e-3/2)**2

lsel,s,line,,1,2,1
lsel,a,line,,14
cm,beam,line
lsel,s,line,,15,24,1
cm,anchor,line
asel,s,area,,3,10,1
cm,A_ROCK1,area
asel,s,area,,1
cm,B_ROCK,area
asel,s,area,,2
cm,A_ROCK2,area

CMSEL,S,anchor,LINE
LATT,2,2,3
mshkey,0
LESIZE,ALL,0.1
mshape,1,2D
LMESH,ALL

lsel,s,line,,3,13,1
LESIZE,ALL,,,10
lsel,s,line,,49,50,1
LESIZE,ALL,,,40
lsel,s,line,,51,61,1
LESIZE,ALL,,,10
allsel,all,all
al,3,15,51,49
al,4,15,52,16
al,5,16,53,17
al,6,17,54,18
al,7,18,55,19
al,8,19,56,20
al,9,20,57,21
al,10,21,58,22
al,11,22,59,23
al,12,23,60,24
al,13,24,61,50
ASEL,S,AREA,,11,21,1
AATT,1,,1
mshkey,1
mshape,0
AMESH,ALL

CMSEL,S,beam,LINE
LATT,3,1,2
mshkey,0
LESIZE,ALL,0.1
mshape,1,2D
LMESH,ALL
lsel,s,line,,49,61,1
lsel,a,line,,1,2,1
lsel,a,line,,25,28,1
lsel,a,line,,14
al,all

ASEL,S,AREA,,22
AATT,1,,1
mshkey,0
ESIZE,0.8
mshape,0
AMESH,ALL

CMSEL,S,A_ROCK1,AREA
AATT,1,,1
mshkey,0
ESIZE,2
mshape,0,2D
AMESH,ALL

allsel,all,all

LSEL,S,LINE,,36,42,3
dl,all,,ux
lsel,s,line,,46,48,2
lsel,a,line,,32
dl,all,,ux
lsel,s,line,,31,33,2
lsel,a,line,,37
dl,all,,uy
/solu
ACEL,,9.8
antype,static
pred,on
lnsrch,on
nlgeom,on
nropt,full
OUTRES,ALL,ALL
alls
nsub,6,10
SOLVE

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

一个隧道底板加锚杆的算例，加上锚杆不收敛，请帮我看下命令流已经有2人回复
求助ANSYS 塑性分析时求解不收敛的问题已经有9人回复

1楼 2015-11-12 17:43:23

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

liusizhe66

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 3758.6
红花: 3
帖子: 623
在线: 226.1小时
虫号: 3591631
注册: 2014-12-12
专业: 岩土力学与岩土工程

从新划分一下网格，尽量使网格接近正方形，

发自小木虫Android客户端

赞一下

回复此楼

2楼2015-11-13 00:29:45

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

专家经验: +809
仿真EPI: 1
应助: 1311 (讲师)
金币: 13381.4
红花: 267
帖子: 3531
在线: 473.4小时
虫号: 3644138
注册: 2015-01-14
专业: 动力学与控制
管辖: 仿真模拟

【答案】应助回帖

感谢参与，应助指数 +1

看你的命令流，感觉材料数据有问题，比如，你的三种材料的弹性模量分别是： 1.4E15、17E20 和 21E15，如果你的长度单位是 m，那么钢材的弹性模量也只有 2E11 左右，你的三种材料不知是什么材料，弹性模量如此之大？
另外，材料 1 的非线性部分的数据可能也有问题。

赞一下

回复此楼

3楼2015-11-13 08:34:12

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

kakaqun

新虫 (小有名气)

应助: 4 (幼儿园)
金币: 902.9
帖子: 57
在线: 16.7小时
虫号: 3496291
注册: 2014-10-24
专业: 岩土与基础工程

引用回帖:

3楼: Originally posted by shgao20 at 2015-11-13 08:34:12
看你的命令流，感觉材料数据有问题，比如，你的三种材料的弹性模量分别是： 1.4E15、17E20 和 21E15，如果你的长度单位是 m，那么钢材的弹性模量也只有 2E11 左右，你的三种材料不知是什么材料，弹性模量如此之大？ ...

您好，材料是我试的时候改的，忘了改回原样，我把正常材料的发给您了
finish
/clear
/prep7
k,1, 5.500000,3.5
k,2,-5.500000,3.5
k,3, 0.000000,4
k,5,-5.500000,-4.5
k,6, 5.500000,-4.5
k,8,-4.500000,-4.5
k,9,-4.500000,-8.5
k,10,-3.500000,-4.5
k,11,-3.500000,-8.5
k,12,-2.500000,-4.5
k,13,-2.500000,-8.5
k,14,-1.500000,-4.5
k,15,-1.500000,-8.5
k,16,-0.500000,-4.5
k,17,-0.500000,-8.5
k,18, 0.500000,-4.5
k,19, 0.500000,-8.5
k,20, 1.500000,-4.5
k,21, 1.500000,-8.5
k,22, 2.500000,-4.5
k,23, 2.500000,-8.5
k,24, 3.500000,-4.5
k,25, 3.500000,-8.5
k,26, 4.500000,-4.5
k,27, 4.500000,-8.5
k,100,-5.500000,-8.5
k,101,5.500000,-8.5
k,28, 15.00000, 15.00000
k,29,-15.00000, 15.00000
k,31,-15.00000,-15.00000
k,33, 15.00000,-15.00000
k,37, 45.00000, 15.00000
k,39, 45.00000,-15.00000
k,43, 15.00000, 45.00000
k,45, 45.00000, 45.00000
k,49,-15.00000, 45.00000
k,53,-45.00000, 45.00000
k,55,-45.00000, 15.00000
k,59,-45.00000,-15.00000
k,63,-45.00000,-45.00000
k,65,-15.00000,-45.00000
k,69,15.00000,-45.00000
k,73,45.00000,-45.00000
larc,1,2,3,5.5
l,2,5
l,5,8
l,8,10
l,10,12
l,12,14
l,14,16
l,16,18
l,18,20
l,20,22
l,22,24
l,24,26
l,26,6
l,1,6
l,8,9
l,10,11
l,12,13
l,14,15
l,16,17
l,18,19
l,20,21
l,22,23
l,24,25
l,26,27
l,31,33
l,33,28
l,28,29
l,29,31
l,59,31
l,31,65
l,63,65
l,59,63
l,65,69
l,33,69
l,33,39
l,39,73
l,69,73
l,28,37
l,37,39
l,43,28
l,43,45
l,45,37
l,49,43
l,49,29
l,53,49
l,53,55
l,55,29
l,55,59
l,5,100
l,6,101
l,100,9
l,9,11
l,11,13
l,13,15
l,15,17
l,17,19
l,19,21
l,21,23
l,23,25
l,25,27
l,27,101
lsel,s,line,,1,14,1
al,all
lsel,s,line,,1,14,1
lsel,a,line,,25,28,1
al,all
allsel,all,all
al,29,30,31,32
al,25,30,33,34
al,34,35,36,37
al,26,35,39,38
al,38,40,41,42
al,27,40,43,44
al,44,45,46,47
al,28,29,47,48
!围岩
E_1=2.8e7
U_1=0.4
DENS_1=2620
C_1=1e7
FI_1=38
MP,EX,1,E_1
MP,PRXY,1,U_1
MP,DENS,1,DENS_1
TB,DP,1
TBDATA,1,C_1,FI_1

!锚杆
E_2=17E10
U_2=0.3
DENS_2=7900
MP,EX,2,E_2
MP,PRXY,2,U_2
MP,DENS,2,DENS_2

!衬砌
E_3=21e9
U_3=0.2
DENS_3=2900
MP,EX,3,E_3
MP,PRXY,3,U_3
MP,DENS,3,DENS_3

ET,1,42
KEYOPT,1,3,2
et,2,3
keyopt,2,6,1
ET,3,1
r,1,0.3,0.3/12,0.3
r,2,3.14*(25e-3/2)**2

lsel,s,line,,1,2,1
lsel,a,line,,14
cm,beam,line
lsel,s,line,,15,24,1
cm,anchor,line
asel,s,area,,3,10,1
cm,A_ROCK1,area
asel,s,area,,1
cm,B_ROCK,area
asel,s,area,,2
cm,A_ROCK2,area

CMSEL,S,anchor,LINE
LATT,2,2,3
mshkey,0
LESIZE,ALL,0.1
mshape,1,2D
LMESH,ALL

lsel,s,line,,3,13,1
LESIZE,ALL,,,10
lsel,s,line,,49,50,1
LESIZE,ALL,,,40
lsel,s,line,,51,61,1
LESIZE,ALL,,,10
allsel,all,all
al,3,15,51,49
al,4,15,52,16
al,5,16,53,17
al,6,17,54,18
al,7,18,55,19
al,8,19,56,20
al,9,20,57,21
al,10,21,58,22
al,11,22,59,23
al,12,23,60,24
al,13,24,61,50
ASEL,S,AREA,,11,21,1
AATT,1,,1
mshkey,1
mshape,0
AMESH,ALL

CMSEL,S,beam,LINE
LATT,3,1,2
mshkey,0
LESIZE,ALL,0.1
mshape,1,2D
LMESH,ALL
lsel,s,line,,49,61,1
lsel,a,line,,1,2,1
lsel,a,line,,25,28,1
lsel,a,line,,14
al,all

ASEL,S,AREA,,22
AATT,1,,1
mshkey,0
ESIZE,0.8
mshape,0
AMESH,ALL

CMSEL,S,A_ROCK1,AREA
AATT,1,,1
mshkey,0
ESIZE,2
mshape,0,2D
AMESH,ALL

allsel,all,all
nummrg,all
numcmp,all

LSEL,S,LINE,,36,42,3
dl,all,,ux
lsel,s,line,,46,48,2
lsel,a,line,,32
dl,all,,ux
lsel,s,line,,31,33,2
lsel,a,line,,37
dl,all,,uy
/solu
ACEL,,9.8
antype,static
pred,on
lnsrch,on
nlgeom,on
nropt,full
OUTRES,ALL,ALL
alls
nsub,6,10
SOLVE

赞一下

回复此楼

4楼2015-11-13 09:28:25

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖