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

充填采场flac3d初始应力平衡不收敛 [复制链接]

上一主题 下一主题
离线csuhwf
 

发帖
6
土币
19
威望
1
原创币
0
只看楼主 倒序阅读 使用道具 楼主  发表于: 2014-05-23
在用flac3d模拟充填采场的时候,初始应力平衡一直不收敛,算了几万次最大不平衡应力还有几个兆帕。
希望大家能帮我看下命令流哪里出错了
以下是命令流

new;建立矿体模型
gen zone brick p0=(0,0,0) size 56 40 8 ;底柱
gen zone brick p0=(0,0,48) size 56 40 8 ;顶柱
group dingdizhu ;组合成顶底柱对象
gen zone brick p0=(0,0,8) size 8 40 40 ;矿柱
gen zone brick p0=(16,0,8) size 8 40 40 ;矿柱
gen zone brick p0=(32,0,8) size 8 40 40 ;矿柱
gen zone brick p0=(48,0,8) size 8 40 40 ;矿柱
group kuangzhu range group dingdizhu not ;组合成矿柱对象

;建立充填矿房

;建立人工假底
gen zone brick p0=(8,0,8) size 8 40 2 ;人工假底
gen zone brick p0=(40,0,8) size 8 40 2 ;人工假底
group rengongjiadi range group dingdizhu not group kuangzhu not;组合人工假底对象

;建立胶结充填体
gen zone brick p0=(8,0,10) size 8 40 3 ;胶结充填
gen zone brick p0=(8,0,14) size 8 40 3 ;胶结充填
gen zone brick p0=(8,0,18) size 8 40 3 ;胶结充填
gen zone brick p0=(8,0,22) size 8 40 3 ;胶结充填
gen zone brick p0=(8,0,26) size 8 40 3 ;胶结充填
gen zone brick p0=(8,0,30) size 8 40 3 ;胶结充填
gen zone brick p0=(8,0,34) size 8 40 3 ;胶结充填
gen zone brick p0=(8,0,38) size 8 40 3 ;胶结充填
gen zone brick p0=(8,0,42) size 8 40 3 ;胶结充填

gen zone brick p0=(40,0,10) size 8 40 3 ;胶结充填
gen zone brick p0=(40,0,14) size 8 40 3 ;胶结充填
gen zone brick p0=(40,0,18) size 8 40 3 ;胶结充填
gen zone brick p0=(40,0,22) size 8 40 3 ;胶结充填
gen zone brick p0=(40,0,26) size 8 40 3 ;胶结充填
gen zone brick p0=(40,0,30) size 8 40 3 ;胶结充填
gen zone brick p0=(40,0,34) size 8 40 3 ;胶结充填
gen zone brick p0=(40,0,38) size 8 40 3 ;胶结充填
gen zone brick p0=(40,0,42) size 8 40 6 ;胶结充填
group jiaojie range group dingdizhu not group kuangzhu not group rengongjiadi not ;组合胶结充填对象

;建立浇面充填
gen zone brick p0=(8,0,13) size 8 40 1 ;浇面充填
gen zone brick p0=(8,0,17) size 8 40 1 ;浇面充填
gen zone brick p0=(8,0,21) size 8 40 1 ;浇面充填
gen zone brick p0=(8,0,25) size 8 40 1 ;浇面充填
gen zone brick p0=(8,0,29) size 8 40 1 ;浇面充填
gen zone brick p0=(8,0,33) size 8 40 1 ;浇面充填
gen zone brick p0=(8,0,37) size 8 40 1 ;浇面充填
gen zone brick p0=(8,0,41) size 8 40 1 ;浇面充填

gen zone brick p0=(40,0,13) size 8 40 1 ;浇面充填
gen zone brick p0=(40,0,17) size 8 40 1 ;浇面充填
gen zone brick p0=(40,0,21) size 8 40 1 ;浇面充填
gen zone brick p0=(40,0,25) size 8 40 1 ;浇面充填
gen zone brick p0=(40,0,29) size 8 40 1 ;浇面充填
gen zone brick p0=(40,0,33) size 8 40 1 ;浇面充填
gen zone brick p0=(40,0,37) size 8 40 1 ;浇面充填
gen zone brick p0=(40,0,41) size 8 40 1 ;浇面充填
group jiaomian range group dingdizhu not group kuangzhu not group rengongjiadi not group jiaojie not;组合浇面充填对象

;建立初始应力平衡矿房
gen zone brick p0=(24,0,8) size 8 40 40 ;建立初始模拟矿房




;初始应力平衡
;模型参数
mo mohr
pro bulk 14.58e9 she 8.75e9 fric 43.9 coh 3.56e6 ten 2.18e6 dil 3 range group dingdizhu
pro bulk 14.58e9 she 8.75e9 fric 43.9 coh 3.56e6 ten 2.18e6 dil 3 range group kuangzhu
pro bulk 14.58e9 she 8.75e9 fric 43.9 coh 3.56e6 ten 2.18e6 dil 3 range x 24 32 z 8 48
pro bulk 1.83e9 she 0.9e9 fric 32.5 coh 0.55e6 ten 1e6 dil 3 range group rengongjiadi
pro bulk 1.23e9 she 0.57e9 fric 30 coh 0.45e6 ten 0.5e6 dil 3 range group jiaojie
pro bulk 1.67e9 she 0.77e9 fric 30 coh 0.31e6 ten 1.9e6 dil 3 range group jiaomian

;设置边界条件
fix x range x -.1 .1
fix x range x 59.9 56.1
fix y range y -.1 .1
fix y range y -40.1 -39.9
fix z range z -.1 .1
;设置密度
ini density 3800 range group dingdizhu
ini density 3800 range group kuangzhu
ini density 3800 range x 24 32 z 8 48
ini density 2400 range group rengongjiadi
ini density 2200 range group jiaojie
ini density 2200 range group jiaomian
hist unbal
;设置初始应力条件
apply szz=-22e6 range z=(-.1,.1)
apply szz=-21e6 range z=(55.9,56.1)
apply sxx=-18e6 range x=(-.1, .1)
apply sxx=-18e6 range x=(55.9,56.1)
apply syy=-24e6 range y=(-.1,.1)
apply syy=-24e6 range y=(-39.9,-40.1)
set mec fo 50
set grav 0 0 -10
solve

;绘制初始应力图
plot
add bcontour szz
add axes black
show
离线林海雪狼

发帖
114
土币
2079
威望
3
原创币
0
只看该作者 1楼 发表于: 2016-10-22
只能表示路过,flac真心没用过
三十年来寻刀剑,几回落叶又抽枝。自从一见桃花后,直至如今更不移.
离线x891208

发帖
11
土币
41
威望
0
原创币
0
只看该作者 2楼 发表于: 2016-10-26
我也是想搞采矿充填的模拟 可惜自己还不懂
快速回复
限100 字节
温馨提示:欢迎交流讨论,请勿纯表情、纯引用!
 
上一个 下一个

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

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