fish中对模量按照应力大小赋值的问题

zengyouyuan

这是《FLAC/FLAC3D基础与工程实例》第8.6部分应用实例的的一个修改，原来例子是土体模量随小主应力变化，我自己建的模型：建模----初始应力---开挖隧道----计算----根据小主应力大小改变模量的值，但是会报错：property ‘young’ not found in zone 1，zone 1就是开挖后的null单元，如果直接建有隧道模型就不会报错，请问这个问题出错的原在哪里，在此先谢谢了。命令：

new

; 圆柱体外环绕放射网格1
ge zo radc size 5 1 15 10 &
p0=0 0 0 p1=15 0 0 p2=0 1 0 p3=0 0 15 &
di 2 2 2 2 fill gro tunnel

;圆柱体外环绕放射网格2
ge zo radc size 5 1 15 10 &
p0=0 0 0 p1=0 0 -15 p2=0 1 0 p3=15 0 0 &
di 2 2 2 2 fill gro tunnel

;镜像
ge zo ref nor -1 0 0 or 0 0 0

mo el
prop bu 3e8 sh 1e8

ini dens 2000
set grav 0 0 -10

fix z ran z -15.01 -14.99
fix x ran x -15.01 -14.99
fix x ran x 14.99 15.01
fix y

solve

ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0
mo nu ra gro tunnel

solve

save modify.sav


rest modify.sav

def E_modify
        p_z=zone_head
        d_k=704
        d_n=0.38
        d_pa=101325.0
        loop while p_z # null
                sigma_3=-1.0*z_sig1(p_z)
                E_new=d_k*d_pa*(sigma_3/d_pa)^d_n
                z_prop(p_z,'young')=E_new
                p_z=z_next(p_z)
        endloop
end
E_modify

huangyu1452

初学，有个疑点，还有个建议，不知能否帮到你：
疑点：为什么修改材料参数在开挖之后呢。
建议：印象中用null操作单元，他的各项数据会有残余，具体可参见《FLAC3D实用教程》第10章初始条件，最后一页。不知道是不是这个带来的影响。

zengyouyuan

修改材料参数针对的是对具体的情况，就如同我提到的那本书上的例子，材料的模量能够反映主应力的影响，只是原例子中没有开挖这一步。这个问题出错的原因是null单元没有young式模量的属性，现在我想到的办法是用group来解决，只是计算相对较慢。谢谢你的回复。

zhuliang1987

楼主辛苦了。