24小时热门版块排行榜    

查看: 1769  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 290求调剂 +6 材料专硕调剂; 2026-02-28 7/350 2026-03-01 09:21 by L135790
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
[考研] 284求调剂 +5 天下熯 2026-02-28 5/250 2026-03-01 08:49 by ms629
[考研] 材料类求调剂 +7 wana_kiko 2026-02-28 7/350 2026-03-01 07:55 by ms629
[考研] 298求调剂 +5 axyz3 2026-02-28 5/250 2026-03-01 06:45 by 刘兵
[考研] 285求调剂 +6 满头大汗的学生 2026-02-28 6/300 2026-03-01 06:29 by Trying]
[考研] 317一志愿华南理工电气工程求调剂 +5 Soliloquy_Q 2026-02-28 8/400 2026-02-28 23:36 by xyx2012xyx
[考研] 化工专硕348,一志愿985求调剂 +4 弗格个 2026-02-28 6/300 2026-02-28 22:00 by wang_dand
[考研] 292求调剂 +3 yhk_819 2026-02-28 3/150 2026-02-28 21:57 by gaoxiaoniuma
[考研] 295求调剂 +5 19171856320 2026-02-28 5/250 2026-02-28 21:39 by gaoxiaoniuma
[考研] 311求调剂 +8 南迦720 2026-02-28 8/400 2026-02-28 21:30 by gaoxiaoniuma
[考研] 求调剂 +4 repeatt?t 2026-02-28 4/200 2026-02-28 21:16 by gaoxiaoniuma
[考研] 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
[考研] 0856材料求调剂 +10 hyf hyf hyf 2026-02-28 11/550 2026-02-28 18:50 by 无际的草原
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
[考研] 寻找调剂 +3 LYidhsjabdj 2026-02-28 3/150 2026-02-28 12:59 by miniwendy
[考研] 304求调剂 +5 曼殊2266 2026-02-28 6/300 2026-02-28 12:44 by 迷糊CCPs
[基金申请] 面上可以超过30页吧? +12 阿拉贡aragon 2026-02-22 13/650 2026-02-26 22:09 by Hahaxia
信息提示
请填处理意见