请教遍布节理模型的参数怎么定义
在遍布节理模型中,有几个参数不是很明白,jnx jny jnz 这三个参数是指什么?是不是软弱面的法向和这三个坐标轴的余弦?还是就是只是角度?还有就是手册上有一个例子
例子如下
def test0
tim0=clock
end
def tim
tim=(clock-tim0)/100.0
end
gen zone cyl p0 0 0 0 p1 1 0 0 p2 0 4 0p3 0 0 1 size 2 2 3
gen zone cyl p0 0 0 0 p1 0 0 1 p2 0 4 0p3 -1 0 0size 2 2 3
gen zone cyl p0 0 0 0 p1 -1 0 0p2 0 4 0p3 0 0 -1size 2 2 3
gen zone cyl p0 0 0 0 p1 0 0 -1p2 0 4 0p3 1 0 0 size 2 2 3
;
def p_cons
eps = 0.001
rad = 1.
c_area = 3.
betaj = 1
if beta = 0. then
betaj = 0
end_if
if beta = 90. then
betaj = 0
end_if
end
;
def trisol
loop k (0,18)
kk = k
beta = k * 5.
p_cons
command
model null
model ubiquitous
prop bul 1.e8 shea 7.e7 cohesion 2.e3
prop friction 40. dilation 0. tension 2400.
prop jdip beta jdd 0. jcoh 1.e3 jfric 30. jdil 0. jten 2000.
fix y range y -.1 .1
fix y range y 3.9 4.1
ini yvel1.e-7 range y -.1 .1
ini yvel -1.e-7 range y 3.9 4.1
cyc 4500
print beta
print sigmav
print anasol
print p_err
print tim
end_command
end_loop
end
在上面的这个例子中,我不是很明白那个软弱面的位置是怎么确定的,是不是若面不需要划分网格?
请各位高手给予指教 非常感谢,我也在弄关于遍布节理的模型问题,向各位多多向各位学习了..
页:
[1]