新手上路

- 积分
- 39
- 极客币
- 73
- 主题
- 22
- 帖子
- 35
- 注册时间
- 2018-10-12
- 在线时间
- 2 小时
- 性别
- 保密
|
应该是用Fish语言编写的额外的调用函数,我也在做这方面的工作,堤坝的地震液化分析,这是fish语言编号的代码:
config zextra 20
def get_old_stress
p_z=zone_head
loop while p_z#null
z_extra(p_z,5)=z_pp(p_z)
p_z=z_next(p_z)
endloop
end
get_old_stress
def get_ppr
whilestepping
p_z=zone_head
loop while p_z#null
z_extra(p_z,7)= z_szz(p_z)+z_pp(p_z)
z_extra(p_z,4)=z_szz(p_z)
z_extra(p_z,6)=z_pp(p_z)-z_extra(p_z,5)
z_extra(p_z,9)=z_szz(p_z)
z_extra(p_z,2)=(-z_extra(p_z,6)/(z_extra(p_z,4)+z_extra(p_z,5))
z_extra(p_z,3)=1.0-z_extra(p_z,7)/(z_extra(p_z,4)+z_extra(p_z,5)
p_z=z_next(p_z)
endloop
end
get_ppr
这个也是陈育民老师书中的,但是我在调用的时候总出错,我在找原因,看到你问就顺便告诉你,如果你知道原因还请回复我,谢谢 |
|