24小时热门版块排行榜    

查看: 1770  |  回复: 4

kakaqun

新虫 (小有名气)

[求助] 隧道开挖模拟,不收敛,请帮忙看一下原因 已有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
回复此楼

» 猜你喜欢

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

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

liusizhe66

木虫 (正式写手)

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

发自小木虫Android客户端
2楼2015-11-13 00:29:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
看你的命令流,感觉材料数据有问题,比如,你的三种材料的弹性模量分别是:  1.4E15、17E20 和 21E15,如果你的长度单位是 m,那么钢材的弹性模量也只有 2E11 左右,你的三种材料不知是什么材料,弹性模量如此之大?
  另外,材料 1 的非线性部分的数据可能也有问题。
3楼2015-11-13 08:34:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kakaqun

新虫 (小有名气)

引用回帖:
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的回帖

kakaqun

新虫 (小有名气)

引用回帖:
2楼: Originally posted by liusizhe66 at 2015-11-13 00:29:45
从新划分一下网格,尽量使网格接近正方形,

问题可能出在锚杆,不加锚杆只开挖就能收敛~~
5楼2015-11-13 09:29:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kakaqun 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 279求调剂 +3 dua1 2026-03-01 4/200 2026-03-02 00:23 by 大脸蛋子
[考研] 292求调剂 +6 yhk_819 2026-02-28 6/300 2026-03-01 23:23 by 向上的胖东
[考研] 0856材料与化工,270求调剂 +6 YXCT 2026-03-01 6/300 2026-03-01 23:21 by 向上的胖东
[考研] 0854复试调剂 276 +3 wmm9 2026-03-01 3/150 2026-03-01 23:13 by 热情沙漠
[考研] 265分求调剂不调专业和学校有行学上就 +6 礼堂丁真258 2026-02-28 8/400 2026-03-01 22:50 by jian_
[基金申请] 成果系统访问量大,请一小时后再尝试。---NSFC啥时候好哦,已经两天这样了 +4 NSFC2026我来了 2026-02-28 4/200 2026-03-01 22:37 by 铁门栓
[考研] 0856调剂 +5 刘梦微 2026-02-28 5/250 2026-03-01 22:30 by wang_dand
[考研] 299求调剂 +3 Y墨明棋妙Y 2026-02-28 5/250 2026-03-01 21:01 by tangxiaotian
[考研] 化工299分求调剂 一志愿985落榜 +5 嘻嘻(*^ω^*) 2026-03-01 5/250 2026-03-01 19:47 by 无际的草原
[考研] 298求调剂 +6 axyz3 2026-02-28 6/300 2026-03-01 19:00 by 18137688336
[考研] 0856材料求调剂 +11 hyf hyf hyf 2026-02-28 12/600 2026-03-01 18:57 by 18137688336
[考研] 281求调剂 +4 2026计算机_诚心 2026-03-01 7/350 2026-03-01 17:20 by 2026计算机_诚心
[考研] 311求调剂 +6 亭亭亭01 2026-03-01 6/300 2026-03-01 15:41 by 324616
[考研] 303求调剂 +4 今夏不夏 2026-03-01 4/200 2026-03-01 14:46 by 嘟嘟小浣熊
[考研] 298求调剂 +9 人间唯你是清欢 2026-02-28 12/600 2026-03-01 14:23 by Ducount.Y
[考研] 材料284求调剂,一志愿郑州大学英一数二专硕 +10 想上岸的土拨鼠 2026-02-28 10/500 2026-03-01 14:12 by yc258
[考研] 311求调剂 +9 南迦720 2026-02-28 10/500 2026-03-01 10:55 by sunny81
[考研] 085600材料工程一志愿中科大总分312求调剂 +8 吃宵夜1 2026-02-28 10/500 2026-02-28 20:27 by L135790
[考研] 276求调剂 +3 路lyh123 2026-02-28 4/200 2026-02-28 19:45 by 路lyh123
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
信息提示
请填处理意见