求救！加爆破地震波以后模型很快就可以平衡，算都算不出来，怎么回事呢

shanhonglan

急！求救各位斑竹和学习flac 的朋友
我做了一个flac 3d 的爆破振动影响下的地下巷道围岩稳定性的例子。原本是期望在加上爆破地震波以后，模型的应力应变、速度加速度、位移和塑性区等能有所变化。但是发现加上爆破波以后，在不使用 set dynamic multi  on 命令时，模型的不平衡很大，在10的正4次方左右，但是不平衡率又很小，在10的负6次方，所以用solve命令只运行了两个时步就不再计算了，结果当然是什么都没有；后来用step 命令运行了40000个时步，塑性区依旧没有变化。在使用了set dynamic multi on 命令后，模型的不平衡力从10的正4次方很快衰减到0，但是不平衡率又很大，为1。时间步长也很小，系统默认的为10的负9次方。不知道出了什么问题，请各位帮忙看一看！
由于学习的时间不长，爆破波的加载和边界条件都是依据例题加的，也不知道加的对不对。
附 动力分析命令流
new
restore kaiwa5.sav
set large dyn=on
；定义荷载,读入速度时程
table 1 read xvel1_hist.txt
；将位移和速度清零
ini xvel=0 yvel=0 zvel=0
ini xdisp=0 ydisp=0 zdisp=0
；set dynamic dt 1.25e-4   ；设置动力时步
；定义阻尼参数
set dynamic damp rayleigh 0.05 24.1
；定义边界并施加速度动荷载
free x y z range z -50.1 -49.9
apply nquiet squiet dquiet range z -50.1 -49.9
apply xvel=3178 hist table 1 range z=-50.1 -49.9
apply ff

；整理时程纪录
his gp xvel 0 0 5
his gp xacc 0 0 5
his dytime

solve
save  x动力.sav

时间的主人

我用了你的地震波数据换了一个模型算了一下,也是只算了20几步就结束了.我感觉是你的地震波数据是不是有问题!!

lakewater

将这一句去掉
free x y z range z -50.1 -49.9

freddie_1999

看了一下模型：
1.模型取的太大了，如果没个隧洞间距比较大的话，可以认为相互不受影响；
2.网格大小比率太大了，使用了很多的attach，估计会影响一定的精度；
3.命令
new
restore kaiwa5.sav
set large dyn=on
；定义荷载,读入速度时程
table 1 read xvel1_hist.txt
；将位移和速度清零
ini xvel=0 yvel=0 zvel=0
ini xdisp=0 ydisp=0 zdisp=0
；set dynamic dt 1.25e-4   ；设置动力时步
；定义阻尼参数
set dynamic damp rayleigh 0.05 24.1
；定义边界并施加速度动荷载
free x y z range z -50.1 -49.9
apply nquiet squiet dquiet range z -50.1 -49.9
apply xvel=3178 hist table 1 range z=-50.1 -49.9
apply ff

；整理时程纪录
his gp xvel 0 0 5
his gp xacc 0 0 5
his dytime

step 200
save  x动力.sav

[ 本帖最后由 freddie_1999 于 2007-1-14 21:03 编辑 ]

时间的主人

我有一个问题:不是说当底部采用静边界的时候应该采用应力时程输入,而程序中好象施加的不是应力时程啊.是否应该将APPLY XVEL=......修改为APPLY SXY=............?
我的意见,请斑竹能解答一下!!!
密切关注中!

时间的主人

请问lakewater斑竹,为什么要将free x y z range z -50.1 -49.9这一句去掉??????

freddie_1999

原帖由 时间的主人 于 2007-1-15 13:06 发表
我有一个问题:不是说当底部采用静边界的时候应该采用应力时程输入,而程序中好象施加的不是应力时程啊.是否应该将APPLY XVEL=......修改为APPLY SXY=............?
我的意见,请斑竹能解答一下!!!
密切关注中!

底部采用粘性边界时，动力激励采用速度，加速度以及应力形式都可以，只不过前两种形式下，运动方向的粘性边界不起作用罢了

zhuhaha

高手，又学一招，谢谢，我也在学习动力问题，不过FLAC施加的是速度等，不象LS-DYNA可以直接加爆炸单元，感觉FLAC模拟波对物体的影响还行，但是真实模拟钻爆开挖，比如模拟微差爆破可能不是很好哦

图拉

原帖由 lakewater 于 2007-1-14 19:46 发表
将这一句去掉
free x y z range z -50.1 -49.9

是不是因为施加的x方向的速度    所以不需要free x y z这个语句？？

zhuhaha

能否解释一下他的速度文件的格式啊：

wave1.vel
300 0.005
0
0.01029
0.02058
0.01029
0
-0.02058
0
0.0063
0
-0.0063
0
0.0147
0
-0.01008
0.00504
-0.00336
0.00441
0
.
.
.

baby20023402

第一行为注释行

第二行为离散点个数和时间间隔
后面为加速度。。。。

yinyunteng

加速度少几行或多几行可以吗？是不是每个加速度值要和模型中的位置相对应呢？而这里怎么看不出来啊？

chensha

建议：
1.从最简单的模型开始，不要第一步就把大变形、动力、ATTACH全部放倒一起；首先去掉这么多ATTACH。
2. 检查模型的网格尺寸是否满足波的传输要求；

一米阳光！

学习学习

一米阳光！

我运行完建模程序后，再运行动力部分时，怎么提示"not configured for dynamic caculation"啊？？？网高手指点一二，想好好学习学习。

yanqing63

请问，这个动力文件的输入还有别的方法没？
比如我有一组实测数据：格式如下   这个文件怎么改成适合加载到FLAC3D里运算？
实测数据：格式
  3.0000000
通道1
通道2
通道3
m/s
m/s
m/s
   101024.000000    10000.000000    -1024.000000
       10.000000       10.000000       10.000000
        0.000000        0.000000        0.000000
       28.000000       28.000000       28.000000
       -0.102300       -0.000011        0.000000        0.000000
       -0.102200        0.000000        0.000000        0.000022
       -0.102100       -0.000011        0.000000       -0.000011
       -0.102000        0.000011        0.000000        0.000022
       -0.101900        0.000000       -0.000011        0.000000
       -0.101800        0.000011        0.000011       -0.000033
       -0.101700        0.000000       -0.000011        0.000011
       -0.101600        0.000000       -0.000022        0.000000
       -0.101500        0.000011        0.000000        0.000000
       -0.101400       -0.000011       -0.000022        0.000011
       -0.101300        0.000011       -0.000011        0.000011
       -0.101200       -0.000022        0.000000       -0.000011
       -0.101100        0.000000        0.000000        0.000000
       -0.101000       -0.000011       -0.000011        0.000011
       -0.100900        0.000022        0.000011        0.000000
       -0.100800        0.000000        0.000011        0.000011
       -0.100700        0.000011        0.000000       -0.000011
       -0.100600        0.000022        0.000033        0.000011
                            中间省略  格式一样
                  末尾为
10.005811        0.000011        0.000000        0.000000
       10.005911        0.000000       -0.000011        0.000022
       10.006011        0.000000        0.000000        0.000011
       10.006111       -0.000011        0.000000        0.000022
       10.006211        0.000011        0.000011        0.000011
       10.006311        0.000000        0.000000        0.000022
关切中等待！

wangshuyang

收获很多！

fidelia1106

ddddddddddddddddddddddddddddddd

golgota

学了不少东西

zbhualala

应该是去掉粘性边界还是将FIX X  Y Z RANGE 这句去掉啊？