我是的新手。我正在练习一个浅基础的稳定性分析。有几行命令不明白。请个为大哥帮忙解释。
==========================================> 我的问题
;---------------------------------------------------------------------
; -二维的条形基础在Tresca 材料中的塑性流动
;---------------------------------------------------------------------
gen zone brick size 20 1 10
;建立矩形的网格区域，其大小为20×1×10
model mohr
;采用摩尔-库仑模型
prop bul 2.e8 shea 1.e8 cohesion 1.e5
;模型的材料性质：体积模量2.0×108 剪切模量 1.0×108 粘聚力1.0×105
prop friction 0. dilation 0. tension 1.e10
;内摩擦角 0° 剪胀角 0° 抗拉强度1.0×1010
fix x range x -.1 .1
;在x 方向上固定边界 x =0（即只在x 方向上受约束）
fix x y z range z -.1 .1
;在x y z 方向上固定边界 z =0（即为固定约束）
fix x y z range x 19.9 20.1
;在x y z 方向上固定边界 x =20
fix y
;约束所有y 方向上运动
fix x y z range x -.1 3.1 z 9.9 10.1 ==========================================>为什么有这个fixity?
;在x y z 方向上固定平面 x =0，3 z= 10
ini zvel -0.5e-5 range x -.1 3.1 z 9.9 10.1==========================================>什么意识?
;在平面 x =0，3 z= 10 上初始化z 方向的速度为−0.5×10−5
def p_cons
;用fish 语言定义函数p_cons==========================================>什么是函数?
pdis1 = gp_near(0.,0.,10.)
;靠近坐标（0，0，10）的节点的地址赋予pdis1 ==========================================>为什么两个同样的地址?
pdis2 = gp_near(0.,1.,10.) ;
;靠近坐标（0，0，10）的节点的地址赋予pdis2 ==========================================>为什么两个同样的地址?
p_sol = (2. + pi)
; 2 +π 赋予 p_sol
end
;用于fish 语言的结尾
p_cons
;运行函数p_cons
;---------------------------------------------------------------------
; p_load : average footing pressure / c
;场地的平均压力
; c_disp : magnitude of vertical displacement at footing center / a
;场地中心的垂直沉降量
;---------------------------------------------------------------------
def actLoad
;定义函数 actLoad（活荷载）（fish 语言的起始句）==========================================>什么是活荷载?
pnt = gp_head==========================================>gp_head是什么value?
;把网格点指针的首个地址（gp_head）赋予变量 pnt==========================================>不明白
pload = 0.0
;定义变量 pload 为0
n = 0
;定义变量 n 为0
loop while pnt # null==========================================>可以解释一下这个loop的运做吗？
;loop 是循环命令，其条件是pnt ≠ 0
if gp_zpos(pnt) > 9.9 then
;当变量 pnt 所对应 的z 坐标大于9.9
if gp_xpos(pnt) < 3.1 then
;且x 坐标小于3.1 时
pload = pload + gp_zfunbal(pnt)
;变量 pload 为上次循环中的值与变量 pnt 所对应的z 方向上的不平衡力的和
n = n+ 1
;n 自我加1
endif
;结束里面的if
endif
; 结束外面的if
pnt = gp_next(pnt)
; 把网格点指针的下一个地址（gp_next）赋予变量 pnt
endloop
;结束循环
actLoad = pload / (3.5 * z_prop(zone_head,'cohesion'))
;这里函数actLoad 的值，为后面的式子所赋予，其中pload 为循环结束后的最终
;值z_prop(zone_head,'cohesion')为首个单元的粘聚力的数值
c_disp = -(gp_zdisp(pdis1) + gp_zdisp(pdis2)) / 7.0
; 同理，这里函数c_disp 的值，也为后面的式子所赋予，其中gp_zdisp(pdis1)为
;变量pdis1 所对应的网格点在z 方向的位移，gp_zdisp(pdis2) 变量pdis2 所对应
;的网格点在z 方向的位移
end
; fish 语言的结束句
def p_err
;定义函数 p_err
p_err = 100. * (actLoad - p_sol) / p_sol
;函数 p_err 的值为后面式子所定义
end
;fish 语言的结束句
hist n 50
;对n 做历史记录
hist actLoad
;对actLoad 做历史记录
hist p_sol
;对p_sol 做历史记录
hist c_disp
;对c_disp 做历史记录
hist unbal
;对unbal（不平衡力）做历史记录
plot sketch
;绘出轮廓线
plot add vel
;添加矢量
plot add his 1 2 vs 3
;绘出历史记录n，actLoad ， p_sol
cyc 15000==========================================不明白要几个cycle. 有reference 吗？
;运行15000 时步
print p_err
;输出函数 p_err
save pran.sav
;形成sav 文件，并保存为pran.sav
ret
;放在批处理文件的最后，以返回FLAC3D 的控制状态
谢谢
我的email
oyh48@hotmail.com
