谁有用fish建不规则几何模型的例子啊？彭文斌141页运行之后和书上给的模型不一样啊。命令流如下：
;不规则表面初始应力状态
new
gen zone brick size 15 15 10 p0 0 0 0 edge = 100.0
model elastic
prop shear 3e8 bulk 5e8
def mountain
gp =gp_head
loop while gp # null
zz = sqrt(gp_xpos (gp)^2 + gp_ypos (gp)^2)
;傅里叶条件求和
dz = 0.06 * sin(0.2 * zz +100.0)
;伪随机表面
dz = dz + 0.06 * sin(0.22 * zz -20.3)
dz = dz - 0.04 * sin(0.33 * zz +33.3)
gp = gp_next (gp)
end_loop
end
mountain
fix x range x -0.1 0.1
fix x range x 99.9 100.1
fix y range y -0.1 0.1
fix y range y 99.9 100.1
fix z range z -0.1 0.1
set grav 0 0 -10
ini density = 2000
ini szz = -2.0e6(grad 0 0 2.0e4)
ini sxx = -4.0e6 (grad 0 0 4.0e4) syy = -4.0e6(grad 0 0 4.0e4)
solve
