关于flac3D书上“15-1冰碛土三轴数值模拟试验”代码的疑问。
前提:按此代码运行确实可以得到预期的结果
但我按书上的代码自己写了一段,却得不到想要的结果
其中有几句关于边界条件的代码我深有疑问,先把书上的代码附上:
n
;导入分组网格
impgrid 1
group soil range group rock not
;定义材料性质
model mohr
pro den 2200 bulk 4e7 she 1.5e7 co 8e4 &
fric 28 ten 9e3 dil 10 range group soil
pro den 2900 bulk 2.5e10 she 1.2e10 co 4e6 &
fric 45 ten 2e6 dil 15 range group rock
;设置边界条件
fix x y z range z -0.1 0.1
fix x y z range z 7.9 8.1; =========Q1=========
fix y
;加载
app nstress -5e5 range x -0.1 0.1
app nstress -5e5 range x 4.9 5.1;=========Q2======
;设置初始条件
ini zvel -2.5e-5 range z 7.9 8.1
;定义轴向应力求解函数
def axi_stress
。。。。。
。。。。。
关于Q1:此处将实验模型用于加载的上端固定住,貌似不能继续加载了啊?为啥要有这一句呢?
关于Q2:这才是我最纳闷儿的。明明所建模型x方向的尺度为0.1m, 他怎么在固定x值为5的地方呢?
谢谢各位达人