跟大家请教并探讨输出某个单元应变的问题
new;建立模型
gen zone brick p0=(0,0,0) p1=(5,0,0) p2=(0,5,0) p3=(0,0,5) size 5 5 5
;模型参数
model mohr
def look_d
p_z=zone_head
loop while p_z # null
z_prop(p_z,'bulk')=urand*1e8
z_prop(p_z,'shear')=urand*0.3e8
z_prop(p_z,'friction')=urand*35
z_prop(p_z,'cohesion')=urand*1e10
z_prop(p_z,'tension')=urand*1e10
p_z=z_next(p_z)
end_loop
end
look_d
;边界条件
fix z range z -0.1 0.1
;初始条件
set gravity 0,0,-9.81
initial density=1000
apply szz=-2e7 range z 4.9 5.1
;解
solve
;寻找应力应变最大块
def find_i
array ar(6)
p_z=zone_head
loop while p_z # null
dum=z_fsi(p_z,ar)
max=ar(3)
li=z_szz(p_z)
smax=max*li
if tt<smax then
tt=smax
zid=z_id(p_z）
end if
p_z=z_next(p_z)
end loop
ii=out('单元体zid号=’string(zid))
kk=out('应力乘应变最大='tt)
end
find_i
这事源代码，一下附文本文档，但当程序运行至
;寻找应力应变最大块
def find_i
array ar(6)
p_z=zone_head
loop while p_z # null
dum=z_fsi(p_z,ar)
max=ar(3)
li=z_szz(p_z)
smax=max*li
处时会出现如图问题
当我输入pr fish后命令栏显示
我想请教各位牛人，这是什么情况
该如何解决
我程序的目的是赋给每个单元格随即的属性
然后寻找到z轴方向应力乘应变最大的那个单元格
很感谢
页:
[1]