论坛首页    职业区    学术与写作    工程技术区    软件区    资料区    商务合作区    社区办公室
 您好! 欢迎 登录注册 最新帖子 邀请注册 活动聚焦 统计排行 社区服务 帮助
 
  • 帖子
  • 日志
  • 用户
  • 版块
  • 群组
帖子
  • 8010阅读
  • 24回复

请各位高手解答一下FLAC-3D的一个小问题,谢谢 [复制链接]

上一主题 下一主题
离线fengzm
 

发帖
99
土币
1167
威望
845
原创币
0
只看楼主 倒序阅读 使用道具 楼主  发表于: 2010-08-13
在FLAC-3D建模中,老出现这样的错误,是什么意思?
Mach error:zero divide invalid
File:.\Solve.cpp
Line:1276

;建立隧道衬砌结构及内部土体部分模型
gen zone cshell p0 0 0 0 p1 3.1 0 0 p2 0 180 0 &
p3 0 0 3.1 size 1 20 8 4 dim 2.75 2.75 2.75 2.75 &
rat 1 1 1 1 group concretliner fill group insiderock
;形成隧道衬砌结构及内部土体部分模型
gen zone reflect dip 90 dd 90 orig 0 0 0 range &
group concretliner
gen zone reflect dip 90 dd 90 orig 0 0 0 range &
group insiderock
gen zone reflect dip 0 dd 0 orig 0 0 0 range &
group concretliner
gen zone reflect dip 0 dd 0 orig 0 0 0 range &
group insiderock
;隧道周围土体模型
;上方覆
gen zone radcyl p0 0 45 0 p1 21.3 45 0 p2 0 135 0 &
p3 0 45 23.1 P4 21.3 135 0 P5 0 135 3.1 &
P6 21.3 45 23.1 P7 21.3 135 3.1 &
size 4 10 8 4 dim 3.1 3.1 3.1 3.1 &
rat 1 1 1 1.2 group outsiderock1
;下方覆土
gen zone radcyl p0 0 0 0 p1 0 0 -12.4 p2 0 180 0 &
p3 21.3 0 0 size 4 20 8 4 dim 3.1 3.1 3.1 3.1 &
rat 1 1 1 1.2 group outsiderock2
;镜像右边到左边
gen zone reflect dip 90 dd 90 orig 0 0 0 range &
group outsiderock1
gen zone reflect dip 90 dd 90 orig 0 0 0 range &
group outsiderock2
plot set center auto
plot set rotation 0 0 90
plot block group
plot add axes black
;边界条件
set gra 0 0 -10
fix z range z -12.41 -12.39
fix x range x -21.31 -21.29
fix x range x 21.29 21.31
fix y range y -0.01 0.01
fix y range y 179.99 180.01
;隧道周围土体的模型及材料参数
model mo
ini density 1800
prop bulk=1.3e8 shear=1.2e7 fric=30 coh=5.0e4 tension=1.0e4
;自重应力场求解
set mech ratio 1e-4
solve
plot cont zdisp outl on

[ 此帖被fengzm在2010-08-14 19:58重新编辑 ]
离线haiqyantu

发帖
186
土币
394
威望
226
原创币
0
只看该作者 1楼 发表于: 2010-08-14
没有命令流,很难判断。
莫非你的是导入的网格?
可能是存在零厚度单元,或节点不耦合。
可以pl sketch或pl bl gr 看下所显示的网格形状是否你所想要的效果。
半生闲吟今终止,
一步江湖无尽期。
离线wind328

发帖
332
土币
528
威望
1143
原创币
0
只看该作者 2楼 发表于: 2010-08-14
基本上是材料参数没有赋全!
我是念水利的,我的专长是能依据水沟内杂草的生长状况判断这条水沟到底有多久没疏浚,却不能一眼看出女孩子到底有多长时间没有谈恋爱了。
离线fengzm

发帖
99
土币
1167
威望
845
原创币
0
只看该作者 3楼 发表于: 2010-08-14
回 1楼(haiqyantu) 的帖子
你好,谢谢你的回答,我附上命令流了,您再帮我看看吧,看我是不是建模就根本出现了问题。谢谢
离线fengzm

发帖
99
土币
1167
威望
845
原创币
0
只看该作者 4楼 发表于: 2010-08-14
回 2楼(wind328) 的帖子
您好,我附上命令流了,您再帮我看看是不是建模出了问题啊。
离线haiqyantu

发帖
186
土币
394
威望
226
原创币
0
只看该作者 5楼 发表于: 2010-08-14
  运行了一下LZ的命令流,说实话,我是没大明白LZ到底是要模拟什么样的工况,直接出来就是下图所示的模型。




  
半生闲吟今终止,
一步江湖无尽期。
离线haiqyantu

发帖
186
土币
394
威望
226
原创币
0
只看该作者 6楼 发表于: 2010-08-14
我觉得问题可能出在下图中红色线圈中的部分,应该出现了零厚度单元,因为此处上覆土的高度与隧道衬砌的高度一样啊,都是z=3.1m吧。
不了解实际工况,故此难以理解LZ如何建立了这样的模型。
半生闲吟今终止,
一步江湖无尽期。
离线haiqyantu

发帖
186
土币
394
威望
226
原创币
0
只看该作者 7楼 发表于: 2010-08-14
假如将上图红色圈中部分,即上覆土y=135m处,改为与y=45m处同样高度的网格,则不再存在此问题。


;上方覆
gen zone radcyl p0 0 45 0 p1 21.3 45 0 p2 0 135 0 &
p3 0 45 23.1 P4 21.3 135 0 P5 0 135 23.1 &
P6 21.3 45 23.1 P7 21.3 135 23.1 &
size 4 10 8 4 dim 3.1 3.1 3.1 3.1 &
rat 1 1 1 1.2 group outsiderock1

得到如下模型,则运行没有问题。可见,就是由于上图中y=135m处的零厚度单元导致错误。

半生闲吟今终止,
一步江湖无尽期。
离线fengzm

发帖
99
土币
1167
威望
845
原创币
0
只看该作者 8楼 发表于: 2010-08-15
回 7楼(haiqyantu) 的帖子
谢谢你的回答。
我就是把隧道看做无限长梁,上边覆土为三角形荷载,所以就建成这样的了。你所说的红色的部分,就是与隧道一样高,应该是出现你所说的零厚度单元,所以计算自重应力场下隧道的的竖向变形时,程序就进行不下去了。那这个问题我应该怎么解决呢?
离线fengzm

发帖
99
土币
1167
威望
845
原创币
0
只看该作者 9楼 发表于: 2010-08-15
回 7楼(haiqyantu) 的帖子
如果没有解决的办法,我只能改变我的论文实例:隧道埋深10m,上方加有三角形荷载,应该不会出现零厚度单元了。谢谢你啊。FLAC高手,想加你为好友。QQ420235837
离线haiqyantu

发帖
186
土币
394
威望
226
原创币
0
只看该作者 10楼 发表于: 2010-08-15
如果你非要建立这个样子的模型的话,恐怕需要分成几块来建立上覆土的网格。
我中午再试下,你也可以多尝试,多考虑下其他简化方法。
半生闲吟今终止,
一步江湖无尽期。
离线haiqyantu

发帖
186
土币
394
威望
226
原创币
0
只看该作者 11楼 发表于: 2010-08-15
试了一下,LZ可考虑将上覆土分成两部分进行建模。
一部分与隧道等高,采用radcyl单元;
另一部分采用wedge单元(楔形体单元),模拟三角形荷载。

其中,与隧道等高部分网格,参考图如下

半生闲吟今终止,
一步江湖无尽期。
快速回复
限100 字节
温馨提示:欢迎交流讨论,请勿纯表情、纯引用!
 
上一个 下一个

      浙公网安备 33010602003799号 浙ICP备14021682号-1

      工业和信息化部备案管理系统网站