关于地震波的传播方向不同时apply ff的设置

herman

1.我的第一个问题：
今天后来仔细检查了下我FLAC数据文件，原来是有些节点不对应，后来改了下，但有点问题我不知道还能不能继续算
gen zone cshell p0 0 0 0 p1 8.7 0 0 p2 0 40 0 p3 0 0 8.7 p4 8.7 40 0 p5 0 40 8.7 p8 5.2 0 0  p9 0 0 5.2 p10 5.2 40 0 p11 0 40 5.2 size 5 20 9 ;
gen zon cshell p0 0 0 0 p1 5.2 0 0 p2 0 40 0 p3 0 0 5.2 p4 5.2 40 0 p5 0 40 5.2 p8 4.85 0 0 p9 0 0 4.85 p10 4.85 40 0 p11 0 40 4.85 size 1 20 9 group rock1 ;
gen zon cshell p0 0 0 0 p1 4.85 0 0 p2 0 40 0 p3 0 0 4.85 p4 4.85 40 0 p5 0 40 4.85 p8 4.5 0 0 p9 0 0 4.5 p10 4.5 40 0 p11 0 40 4.5 size 1 20 9 group rock2 ;
gen zon cshell p0 0 0 0 p1 4.200223 0 -2.425 p2 0 40 0 p3 4.85 0 0 p4 4.200223 40 -2.425 p5 4.85 40 0 p8 3.897114 0 -2.25 p9 4.5 0 0 p10 3.897114 40 -2.25 p11 4.5 40 0 size 1 20 3 group rock2 ;      
gen zon cshell p0 0 0 0 p1 4.503332 0 -2.600 p2 0 40 0 p3 5.20 0 0 p4 4.503332 40 -2.600 p5 5.20 40 0 p8 4.200223 0 -2.425 p9 4.85 0 0 p10 4.200223 40 -2.425 p11 4.85 40 0 size 1 20 3 group rock1 ;
gen zon cshell p0 0 0 0 p1 7.534421 0 -4.350 p2 0 40 0 p3 8.70 0 0 p4 7.534421 40 -4.350 p5 8.70 40 0 p8 4.503332 0 -2.600 p9 5.20 0 0 p10 4.503332 40 -2.600 p11 5.20 40 0 size 5 20 3 ;
gen zon cshell p0 0 0 0 p1 0 0 -3.550 p2 0 40 0 p3 4.200223 0 -2.425 p4 0 40 -3.550 p5 4.200223 40 -2.425 p8 0 0 -3.200 p9 3.897114 0 -2.25 p10 0 40 -3.200 p11 3.897114 40 -2.25 size 1 20 6 group rock2 ;
gen zon cshell p0 0 0 0 p1 0 0 -3.900 p2 0 40 0 p3 4.503332 0 -2.600 p4 0 40 -3.900 p5 4.503332 40 -2.600 p8 0 0 -3.550 p9 4.200223 0 -2.425 p10 0 40 -3.550 p11 4.200223 40 -2.425 size 1 20 6 group rock1;
gen zon cshell p0 0 0 0 p1 0 0 -7.400 p2 0 40 0 p3 7.534421 0 -4.350 p4 0 40 -7.400 p5 7.534421 40 -4.350 p8 0 0 -3.900 p9 4.503332 0 -2.600 p10 0 40 -3.900 p11 4.503332 40 -2.600 size 5 20 6
gen zone cyl p0 0. 0. 0. p1 4.5 0. 0. p2 0. 40. 0. p3 0. 0. 4.5 p4 4.5 40. 0. p5 0. 40. 4.5 size 5 20 9 ratio 1 1 1 group rock3
gen zone cyl p0 0. 0. 0. p1 3.897114 0. -2.25 p2 0. 40. 0. p3 4.5 0. 0. p4 3.897114 40. -2.25 p5 4.5 40. 0 size 5 20 3 ratio 1 1 1 group rock3
gen zone cyl p0 0. 0. 0. p1 0. 0 -3.20 p2 0. 40. 0. p3 3.897114 0. -2.25 p4 0 40. -3.2 p5 3.897114 40. -2.25 size 5 20 6 ratio 1 1 1 group rock3
gen zon radcyl p0 0 0 0 p1 20.0 0 0 p2 0 40 0 p3 0 0 20 size 5 20 9 8 dim 8.7. 8.7. 8.7. 8.7 ratio 1 1 1 1.2
gen zon radcyl p0 0 0 0 p1 20.0 0 -11.547005 p2 0 40 0 p3 20 0 0 p4 20 40 -11.547005 p5 20 40 0 p6 20 0 -4 p7 20 40 -4 p8 7.534421 0 -4.350 p9 8.7 0 0 p10 7.534421 40 -4.35 p11 8.7 40 0 size 5 20 3 8 ratio 1 1 1 1.2
gen zon radcyl p0 0 0 0 p1 0 0 -20 p2 0 40 0 p3 20 0 -11.547005 p4 0 40 -20 p5 20 40 -11.547005 p6 20 0 -20 p7 20 40 -20 p8 0 0 -7.4 p9 7.534421 0 -4.35 p10 0 40 -7.40 p11 7.534421 40 -4.35 size 5 20 6 8 ratio 1 1 1 1.2
gen zon reflect normal -1 0 0
mo mohr
pro bulk 1.86e8 shear 1.75e8 fric 29 coh 5e5 ten 0.2e6 range z -20 20
ini dens=2780
set grav 0 0 -10
fix z range z -20.0001 -19.9999
fix x range x -20.0001 -19.9999
fix x range x 19.9999 20.0001
fix y range y -0.0001 0.0001
fix y range y 39.9999 40.0001
apply szz -1.0286e6 range z 19.9999 20.0001
ini szz -1.0286e6 grad 0 0 27800
ini syy -3.39438e5 grad 0 0 9174
ini sxx -3.39438e5 grad 0 0 9174
hist unbal
plot hist 1
solve
save step1.sav
算的时候收敛，而且SYY和YDISP也比较好
开挖计算
model null range group rock3;开挖核心部分
model null range group rock2;二期支护处开挖围岩
model null range group rock1;初期支护处开挖围岩
def a1
x1=0.0
x2=0.0
z1=0.0
z2=0.0
y1=0.0
loop i(1,7)
y1=5.0*i
loop n(-2,22)
  x1=5.2*cos((n-1)*0.17454167)
  z1=5.2*sin((n-1)*0.17454167)
  x2=8.7*cos((n-1)*0.17454167)
  z2=8.7*sin((n-1)*0.17454167)
  command
   sel cable id i begin x1 y1 z1 end x2 y1 z2 nseg 5
  end_command
end_loop
end_loop
end
a1
sel cable prop emod 2.0e10 ytension 8.0e4 xcarea 3.8e-4 &   
gr_coh 2e5 gr_k 2.0e7 gr_per 0.3                              
range name ex1 group rock1 any
model elastic range ex1
prop s=1.15e10 b=9.5e8 range ex1  
range name ex2 group rock2 any
model elastic range ex2
prop s=1.64e10 b=1.36e10 range ex2  
hist reset;
hist unbal;   
solve
save excav.sav

我算动力时，是在计算了开挖之后，算从模型底部往上传播的地震波，也就是沿着Z轴正方向
conf dyn
;进行动力分析
set dyn on   
set log on   
ini xvel 0 yvel 0 zvel 0
ini xdis 0 ydis 0 zdis 0
ini state 0
free x y z ran z -20.0001 -19.9999                                                      
;动力边界条件
tab 1 read zv.txt
ini dens=2780
sel cable prop dens=3000
apply ff
set dyn multi on ;加快运算速度
apply zacc 1.0 hist table 1 ran z -20.0001 -19.9999
hist reset
hist dytime       ;                        
set dyn damp local 0.2
set dyn dt 9.999e-6
sol age 5
save dynam1.sav
但是lakewater你不是说加了apply ff后其他情况的约束自动消除，只剩下波的传播方向才有约束吗，但是

herman

这是约束情况图，但隧道的位置有些平行隧道轴向的约束

herman

我怀疑是不是和加了锚杆有关，所以就把算完静力后直接算动力，没算开挖的，这时也算从模型底部往上传播的波，这时就像lakewater说的。既然这个和加锚杆有关，那么我们一般也是算完了开挖才算动力，还要不要加锚杆

herman

以上是第一个问题，下面是我的第二个问题
这是算波从模型底部往上传播的，但现在如果考虑波平行洞轴线从y=0这个面往模型隧道轴向方向转播呢？
我试了下
conf dyn
;进行动力分析
set dyn on   
set log on   
ini xvel 0 yvel 0 zvel 0
ini xdis 0 ydis 0 zdis 0
ini state 0
free x y z ran y -0.0001 0.0001                                                      
;动力边界条件
tab 1 read yv.txt
ini dens=2780
sel cable prop dens=3000
apply ff
set dyn multi on ;加快运算速度
apply yacc 1.0 hist table 1 ran y -0.0001 0.0001
但这时多出来的东西和第一种情况（从模型底部传上来的地震波一样），我怀疑这是不是正确的？

herman

地震波从y=0这个面往隧道轴线传播，但也只是隧道前后左右多加了那个边界，是不是对了哦？这是此时的约束情况，在隧道后边界y=40处有一个地方有y方向约束，但模型底部怎么还有约束

herman

这是用静力计算后得到的直接算水平地震波的情况，y=40处的约束消掉了
不知道你们算动力有没有考虑到这些情况

herman

希望lakewater等人能好好看看这个问题

herman

对于平行洞轴线传播的波，我突然想到用是不是一定要释放模型的底部，然后在底部ran z -20.0001 -19.9999这个范围施加平行洞轴线的地震波free x y z ran z -20.0001 -19.9999
                           
tab 1 read yv.txt
ini dens=2780
sel cable prop dens=3000
apply ff
set dyn multi on ;加快运算速度
apply yacc 1.0 hist table 1 ran z -20.0001 -19.9999
结果是像lakewater 说的那样，我不知道考虑其他方向的地震波时，是不是应该这样做

wyb8254

关注中

lakewater

herman wrote:
但是lakewater你不是说加了apply ff后其他情况的约束自动消除，只剩下波的传播方向才有约束吗

加了ff边界以后，静力平衡的条件在动力分析之前自动转换到ff边界中（因为ff边界上增加了单元），所以称为没有其他约束（比如速度约束条件肯定是不成立的）。
但，我没提过“只剩下波的传播方向”这样的描述，这个边界条件跟波的方向是没关系的。另外，在模型底部施加动力边界（时程）时，要在ff边界设置之前，你的命令中这个顺序有错误！

lakewater

第2个问题：模型建立有错误，见下图，两个cshell命令的单元数不匹配，细看有空隙。
通过gen merge 0.01可以发现存在很多节点空隙，请检查。

lakewater

对于gpfix情况，从下图上看，只发生在开挖掉的单元上，这个对计算应该没有影响。

lakewater

herman wrote:
tab 1 read yv.txt
ini dens=2780
sel cable prop dens=3000
apply ff
set dyn multi on ;加快运算速度
apply yacc 1.0 hist table 1 ran z -20.0001 -19.9999
结果是像lakewater 说的那样，我不知道考虑其他方向的地震波时，是不是应该这样做

首先，这个app ff与app yacc的顺序有错误！
其次，施加动力边界条件与方向无关，关键是施加的位置，比如在模型底部施加加速度时程，任何方向的时程方法都是一样的。

herman

谢谢，我已经明白了许多，越问得多，越显得自己很渺小，看来光靠勤奋还不够，还要掌握方法

riken

是啊，学海无涯，人还是要谦虚啊

dubing222

谢谢楼主吧

kobe333

hao hahao

qbzhang666

学了不少东西

maocw

好东西啊

baby20023402

学习了。。。。。。。。。