动力分析的时候calculating geometric constants 后出现 unknown fatal ...
本帖最后由 曾锦秀 于 2015-2-11 01:15 编辑RT
动力分析的时候calculating geometric constants 后出现 unknown fatal exception
我尝试了删除接触面后计算，可以计算
我还是了其他简单模型，建立接触面动力计算，完全没有问题。
我想问的是我的边坡模型为何独独不能进行动力计算呢？
为了大家交流方便，我把命令流贴出来，还望高手不吝指教
本帖最后由 曾锦秀 于 2015-3-3 12:05 编辑new
impgrid 36万.Flac3d
attach face
gen merge 1e-6
ini y add -10
config dynamic
set dynamic off
pl bl gr
plot set background white
plot set plane normal (0.000,0.000,1.000)
plot set plane origin (0.0000e+000,0.0000e+000,0.0000e+000)
plot set mode model
plot set center (2.5371e+002,1.2053e+002,1.7480e+002)
plot set rotation (20.00, 0.00,50.00)
plot set distance 1.5855e+003
plot set angle 22.50
plot set magnification1.00e+000
group "含砾黏土" range group 50
group "粘土层" range group 46
group "漂卵石层" range group 45
group "强风化闪长岩" range group 43
group "中风化闪长岩" range group 42
group "强风化蚀变二长花岗岩" range group 44
group "中风化蚀变二长花岗岩" range group 41
group "承台1+含砾黏土" range group 48
group "承台2+强风化闪长岩" range group 47
group "承台3" range group 49
group "左隧1+含砾黏土" range group 13
group "左隧2+强风化蚀变二长花岗岩" range group 15
group "左隧3" range group 17
group "左隧4" range group 19
group "左隧5" range group 21
group "左隧6" range group 23
group "左隧7" range group 25
group "左隧8" range group 27
group "左隧9" range group 29
group "左隧10" range group 31
group "左隧11" range group 33
group "左隧12" range group 35
group "左隧13" range group 37
group "左隧14" range group 39
group "右隧1+含砾黏土" range group 14
group "右隧2+强风化蚀变二长花岗岩" range group 16
group "右隧3" range group 18
group "右隧4" range group 20
group "右隧5" range group 22
group "右隧6" range group 24
group "右隧7" range group 26
group "右隧8" range group 28
group "右隧9" range group 30
group "右隧10" range group 32
group "右隧11" range group 34
group "右隧12" range group 36
group "右隧13" range group 38
group "右隧14" range group 40
group "左锚1+含砾黏土" range group 11
group "左锚2+强风化蚀变二长花岗岩" range group 9
group "左锚3" range group 7
group "左锚4" range group 5
group "左锚5+锚塞体" range group 3
group "左锚6" range group 1
group "右锚1+含砾黏土" range group 12
group "右锚2+强风化蚀变二长花岗岩" range group 10
group "右锚3" range group 8
group "右锚4" range group 6
group "右锚5+锚塞体" range group 4
group "右锚6" range group 2
model elas
;>>>>>>>含砾黏土
prop shear 7.25e8 bulk 27.78e8range group 含砾黏土 any &
gr 承台1+含砾黏土 any &
gr 左隧1+含砾黏土 any &
gr 右隧1+含砾黏土 any &
gr 左锚1+含砾黏土 any &
gr 右锚1+含砾黏土 any
ini dens 1900 range group 含砾黏土 any &
gr 承台1+含砾黏土 any &
gr 左隧1+含砾黏土 any &
gr 右隧1+含砾黏土 any &
gr 左锚1+含砾黏土 any &
gr 右锚1+含砾黏土 any
;>>>>>>>漂卵石层
prop shear 125e8 bulk 166.7e8 range group 漂卵石层
ini dens 2200 range group 漂卵石层
;>>>>>>>粘土层
prop shear 5.6e8 bulk 15.63e8range group 粘土层
ini dens 1900 range group 粘土层
;>>>>>>>强风化闪长岩
prop shear 188e8 bulk 490e8 range group 强风化闪长岩 any gr 承台2+强风化闪长岩 any
ini dens 2500 range group 强风化闪长岩 any gr 承台2+强风化闪长岩 any
;>>>>>>>中风化闪长岩
prop shear 7857e8 bulk 13750e8 range group 中风化闪长岩
ini dens 2610 range group 中风化闪长岩
;>>>>>>>强风化蚀变二长花岗岩
prop shear 225.6e8 bulk 588.2e8 range group 强风化蚀变二长花岗岩 any &
gr 左隧2+强风化蚀变二长花岗岩 any &
gr 右隧2+强风化蚀变二长花岗岩 any &
gr 左锚2+强风化蚀变二长花岗岩 any &
gr 右锚2+强风化蚀变二长花岗岩 any
ini dens 2550 range group 强风化蚀变二长花岗岩 any &
gr 左隧2+强风化蚀变二长花岗岩 any &
gr 右隧2+强风化蚀变二长花岗岩 any &
gr 左锚2+强风化蚀变二长花岗岩 any &
gr 右锚2+强风化蚀变二长花岗岩 any
;>>>>>>>中风化蚀变二长花岗岩
prop shear 6496e8 bulk 11956e8 range group 中风化蚀变二长花岗岩 any &
gr 左隧3 any gr 左隧4 any gr 左隧5 any gr 左隧6 any group 左隧7 any &
gr 左隧8 any gr 左隧9 any gr 左隧10 any gr 左隧11 any gr 左隧12 any &
gr 左隧13 any gr 左隧14 any &
gr 右隧3 any gr 右隧4 any gr 右隧5 any gr 右隧6 any group 右隧7 any &
gr 右隧8 any gr 右隧9 any gr 右隧10 any gr 右隧11 any gr 右隧12 any &
gr 右隧13 any gr 右隧14 any &
gr 左锚3 any gr 左锚4 any gr 左锚5+锚塞体 any gr 左锚6 any &
gr 右锚3 any gr 右锚4 any gr 右锚5+锚塞体 any gr 右锚6 any &
ini dens 2610 range group 中风化蚀变二长花岗岩 any &
gr 左隧3 any gr 左隧4 any gr 左隧5 any gr 左隧6 any group 左隧7 any &
gr 左隧8 any gr 左隧9 any gr 左隧10 any gr 左隧11 any gr 左隧12 any &
gr 左隧13 any gr 左隧14 any &
gr 右隧3 any gr 右隧4 any gr 右隧5 any gr 右隧6 any group 右隧7 any &
gr 右隧8 any gr 右隧9 any gr 右隧10 any gr 右隧11 any gr 右隧12 any &
gr 右隧13 any gr 右隧14 any &
gr 左锚3 any gr 左锚4 any gr 左锚5+锚塞体 any gr 左锚6 any &
gr 右锚3 any gr 右锚4 any gr 右锚5+锚塞体 any gr 右锚6 any &
;将来的填土
mo nu ran gr 承台3
;pause
;设置接触面1、2、3=========================================================
gen separate 粘土层
interface 3 wrap 粘土层 漂卵石层
; int 3 maxedge 4
;gen merge 1e-6 ran gr 含砾黏土 any gr 漂卵石层 any
;gen merge 1e-6 ran gr 粘土层 any gr 漂卵石层 any
;===========接触面参数，可以调整=========kn一般比ks的值大一个数量级========
interface 3 prop kn 40e10 ks 20e10 coh 220e3 fric 29 ten 0e3 bslip on
;>>>>>>>设置重力加速度
set gravity 0 0 -9.8
;>>>>>>>设置大变形
;set large
;>>>>>>>>fix boundary边界
fix x y z ran z-0.1 0.1
fix x range x -0.1 0.1
fix x range x 529.9 530.1
fix y range y -0.1 0.1
fix y range y 239.9 240.1
;>>>>>>>>监测
hist n 5
hist unbal
plot his 1
pl add con xdis ou on shade on
solve
;save 1-0地震弹性计算结果.sav
;=======================================
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
model mohr
mo nu ran gr 承台3
;>>>>>>>含砾黏土 提高coh 25e3到40e3，提高fric 33.5到35.5
prop shear 7.25e6 bulk 27.78e6 coh 30e3fric 25 ten 0 range group 含砾黏土 any &
gr 承台1+含砾黏土 any &
gr 左隧1+含砾黏土 any &
gr 右隧1+含砾黏土 any &
gr 左锚1+含砾黏土 any &
gr 右锚1+含砾黏土 any
ini dens 1900 range group 含砾黏土 any &
gr 承台1+含砾黏土 any &
gr 左隧1+含砾黏土 any &
gr 右隧1+含砾黏土 any &
gr 左锚1+含砾黏土 any &
gr 右锚1+含砾黏土 any
;>>>>>>>块石土
prop shear 18.54e6 bulk 50.35e6 coh 10e3fric 39.5 ten 0 range group 含砾黏土 x 0 190
ini dens 2020 range group 含砾黏土 x 0 190
prop shear 18.54e6 bulk 50.35e6 coh 10e3fric 39.5 ten 0 range group 承台1+含砾黏土
ini dens 2020 range group 承台1+含砾黏土
;>>>>>>>块石土，增加范围
prop shear 18.54e6 bulk 50.35e6 coh 10e3fric 39.5 ten 0 range &
group 含砾黏土 x 0 240 y 209 240.5
ini dens 2020 range group 含砾黏土 x 0 240 y 209 240.5
prop shear 18.54e6 bulk 50.35e6 coh 10e3fric 39.5 ten 0 range &
group 含砾黏土 sphere center 192.9 207.3 181.6 radius 46
ini dens 2020 range group 含砾黏土 sphere center 192.9 207.3 181.6 radius 50
prop shear 18.54e6 bulk 50.35e6 coh 10e3fric 39.5 ten 0 range &
group 含砾黏土 sphere center 200.7 40.03 200.9 radius 50
ini dens 2020 range group 含砾黏土 sphere center 200.7 40.03 200.9 radius 50
;>>>>>>>漂卵石层 原先的coh 5e3 fric 45，根据试验资料修改后
prop shear 125e6 bulk 166.7e6 coh 150e3 fric 38 ten 0 range group 漂卵石层
ini dens 2200 range group 漂卵石层
;>>>>>>>粘土层(半成岩)根据试验抗压强度挺大，ten的值提升 coh 30e3 提示到40e3
prop shear 21.8e6 bulk 60.9e6 coh 60e3 fric 26 ten 0.8e3 range group 粘土层
ini dens 1900 range group 粘土层
;>>>>>>>强风化闪长岩ten 由0变成25e3
prop shear 188e6 bulk 490e6 coh 250e3 fric 37 ten 25e3 range &
group 强风化闪长岩 any gr 承台2+强风化闪长岩 any
ini dens 2500 range group 强风化闪长岩 any gr 承台2+强风化闪长岩 any
;>>>>>>>中风化闪长岩
prop shear 7857e6 bulk 13750e6 coh 1500e3 fric 45 ten 6.2e6 range group 中风化闪长岩
ini dens 2610 range group 中风化闪长岩
;>>>>>>>强风化蚀变二长花岗岩ten 由0变成25e3
prop shear 225.6e6 bulk 588.2e6 coh 250e3 fric 30 ten 25e3 range &
gr 强风化蚀变二长花岗岩 any &
gr 左隧2+强风化蚀变二长花岗岩 any &
gr 右隧2+强风化蚀变二长花岗岩 any &
gr 左锚2+强风化蚀变二长花岗岩 any &
gr 右锚2+强风化蚀变二长花岗岩 any
ini dens 2550 range &
gr 强风化蚀变二长花岗岩 any &
gr 左隧2+强风化蚀变二长花岗岩 any &
gr 右隧2+强风化蚀变二长花岗岩 any &
gr 左锚2+强风化蚀变二长花岗岩 any &
gr 右锚2+强风化蚀变二长花岗岩 any
;>>>>>>>中风化蚀变二长花岗岩
prop shear 6496e6 bulk 11956e6 coh 800e3 fric 38 ten 4.3e6 &
range group 中风化蚀变二长花岗岩 any &
gr 左隧3 any gr 左隧4 any &
gr 左隧5 any gr 左隧6 any &
gr 左隧7 any gr 左隧8 any &
gr 左隧9 any gr 左隧10 any &
gr 左隧11 any group 左隧12 any &
gr 左隧13 any gr 左隧14 any &
gr 右隧3 any gr 右隧4 any &
gr 右隧5 any gr 右隧6 any &
gr 右隧7 any gr 右隧8 any &
gr 右隧9 any gr 右隧10 any &
gr 右隧11 any group 右隧12 any &
gr 右隧13 any gr 右隧14 any &
gr 左锚3 any gr 左锚4 any gr 左锚5+锚塞体 any gr 左锚6 any &
gr 右锚3 any gr 右锚4 any gr 右锚5+锚塞体 any gr 右锚6 any
ini dens 2610 &
range group 中风化蚀变二长花岗岩 any &
gr 左隧3 any gr 左隧4 any &
gr 左隧5 any gr 左隧6 any &
gr 左隧7 any gr 左隧8 any &
gr 左隧9 any gr 左隧10 any &
gr 左隧11 any group 左隧12 any &
gr 左隧13 any gr 左隧14 any &
gr 右隧3 any gr 右隧4 any &
gr 右隧5 any gr 右隧6 any &
gr 右隧7 any gr 右隧8 any &
gr 右隧9 any gr 右隧10 any &
gr 右隧11 any group 右隧12 any &
gr 右隧13 any gr 右隧14 any &
gr 左锚3 any gr 左锚4 any gr 左锚5+锚塞体 any gr 左锚6 any &
gr 右锚3 any gr 右锚4 any gr 右锚5+锚塞体 any gr 右锚6 any
plot his 1
pl add con xdis ou on shade on
;>>>>>>>>fix boundary边界
fix x y z ran z-0.5 0.5
fix x range x -0.5 5.1
fix x range x 525.9 530.5
fix y range y -0.5 5.1
fix y range y 235.9 240.5
;>>>>>>>设置大变形=========================================================
set large
plot his 1
pl add con xdis ou on shade on
set log off
solve
;save 1-0地震初始应力场.sav
set dyn on
;set large
pl con xdi ou on
;初始化条件
ini state=0
ini xdis=0 ydis=0 zdis=0
ini xvel=0 yvel=0 zvel=0
free x y z ran z -0.1 0.1 ;去掉模型底部原有的静力条件
apply nquiet squiet dquiet ran z -0.1 0.1;静态边界条件
table 101 read zali_xia_50_10处理后.dat
;table 102 read accy.txt
;table 103 read accz.txt
apply xa 1.0 hist table 101 ran z -0.1 0.1 ;加动力荷载
;apply ya 1.0 hist table 102 ran z -0.1 0.1 ;加动力荷载
;apply za 1.0 hist table 103 ran z -0.1 0.1 ;加动力荷载
set dyn damp local 0.157 ;施加局部阻尼
apply ff ;施加自由场边界条件
group ff_corner
group ff_side ran x 0 530;
group ff_side ran y 0 240;
group main_grid ran x 0 530 y 0 240 ;
set dyn time = 0 ;设置动力计算从0s开始
set dyn multi on;加快运算速度
;;监测位移，速度，加速度
hist reset ;清空已有的历史信息
;>>>>>>>>监测
hist n 5
hist id=1 unb
hist id=2 dytime
plot his 1
;pause
pl add con xdis ou on shade on
solve age 0.2
去掉这段命令流（即不设置接触面）后可以正常计算
;设置接触面1、2、3=========================================================gen separate 粘土层
interface 3 wrap 粘土层 漂卵石层
; int 3 maxedge 4
;gen merge 1e-6 ran gr 含砾黏土 any gr 漂卵石层 any
;gen merge 1e-6 ran gr 粘土层 any gr 漂卵石层 any
;===========接触面参数，可以调整=========kn一般比ks的值大一个数量级========
interface 3 prop kn 40e10 ks 20e10 coh 220e3 fric 29 ten 0e3 bslip on 还是没人啊 我的是结构面的问题，不知道为什么，加了结构面就是不能进行动力运算，不加结构面完全可以。我也不知道问题处到哪了。请问你的问题解决了吗 smilehzg 发表于 2015-8-30 20:55
我的是结构面的问题，不知道为什么，加了结构面就是不能进行动力运算，不加结构面完全可以。我也不知道问 ...
你用什么模拟结构面，如果是接触面模拟，那问题一样 曾锦秀 发表于 2015-9-21 20:55
你用什么模拟结构面，如果是接触面模拟，那问题一样
是的，接触面 学习学习
页:
[1]