论坛首页    职业区    学术与写作    工程技术区    软件区    资料区    商务合作区    社区办公室
 您好! 欢迎 登录 或 注册 最新帖子 邀请注册 活动聚焦 统计排行 社区服务 帮助
 
  • 帖子
  • 日志
  • 用户
  • 版块
  • 群组
帖子
  • 3241阅读
  • 0回复

fish程序中关于zextra的问题 [复制链接]

上一主题 下一主题
离线lxz6510151
 

发帖
30
土币
159
威望
11
原创币
0
只看楼主 倒序阅读 使用道具 楼主  发表于: 2012-12-15
各位高手,下面是我的一段程序,计算地震液化。里面有一个额外变量zextra_1。另外三个变量yl_1, yl_2, c_ppr 计算结束时,我用pri fish命令显示fish程序内变量值,得到yl_2=1.5e6 yl_1=1.3e6 c_ppr=0.13.但是我在flac命令行中输入pl contour zextra_1 时,显示的云图数据竟然是100多~这让我很困惑,我程序想法是让z_extra 1的值等于c_ppr的值想(小于1),但是结果却很大。请高手帮我看看,有什么不对的地方,不胜感激!!!
config zextra 1
def get_old_stress
  p_z = zone_head
  loop while p_z # null
    yl_1 = (z_sig1(p_z) + z_sig2(p_z) + z_sig3(p_z)) / 3.0 + 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
    yl_2 = (z_sig1(p_z) + z_sig2(p_z) + z_sig3(p_z)) / 3.0 + z_pp(p_z)  
    c_ppr =1 - (yl_1)/(yl_2)
    z_extra(p_z,1) = c_ppr
    p_z = z_next(p_z)
  endloop
end
get_ppr

快速回复
限100 字节
温馨提示:欢迎交流讨论,请勿纯表情、纯引用!
 
上一个 下一个

      https://beian.mps.gov.cn/ 粤公网安备 44010602012919号 广州半山岩土网络科技有限公司 粤ICP备2024274469号

      工业和信息化部备案管理系统网站