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

求塑性区体积,大神来帮帮忙 [复制链接]

上一主题 下一主题
离线panpantao
关闭
个人中心可以申请新版勋章哦
立即申请知道了
 

发帖
8
土币
42
威望
1
原创币
0
只看楼主 倒序阅读 使用道具 楼主  发表于: 2014-03-12


def get_plast   ;定义函数get_plast
  shearnow = 1                                     ;shearnow、tensionnow
  tensionnow = 2                                   ;shearpast、tensionpast
  shearpast = 4                                    ;的state value分别为1、2、4、8
  tensionpast = 8                                  ;具体内容参见fish reference 2.5.3.3

  v_shear_now = 0                           ;需求的shearnow、tensionnow
  v_tension_now = 0                         ;shearpast、tensionpast的体积,初始值为0
  v_shear_past = 0
  v_tension_past = 0   ;flac默认的值

;xx=z_state(find_zone(59840),0)   ;通过本步可验证

  p_z = zone_head                  ;单元指针
  loop while p_z # null            ;开始循环,如果单元指针不为空,执行以下操作

    if and(z_state(p_z,0),shearnow) = shearnow  then   ;判断语句,判断单元如果是shearnow状态
        v_shear_now = v_shear_now + z_volume(p_z)      ;则把该单元的体积加进v_shear_now
    endif
    if and(z_state(p_z,0),tensionnow) = tensionnow then  ;同上判断语句,判断单元是否为tensionnow状态
       v_tension_now = v_tension_now + z_volume(p_z)
    endif
    if and(z_state(p_z,0),shearpast) = shearpast then       ;同上,判断单元是否shearpast
       v_shear_past = v_shear_past + z_volume(p_z)
    endif
    if and(z_state(p_z,0),tensionpast) = tensionpast  then
       v_tension_past = v_tension_past + z_volume(p_z)
    endif
    p_z = z_next(p_z)          ;单元指针指向下一个单元
   endloop

   ii = out('shear_now : ' + string(v_shear_now))         ;输出shear_now的体积值
   ii = out('tension_now : ' + string(v_tension_now))     ;输出tension_now的体积值
   ii = out('shear_past : ' + string(v_shear_past))       ;输出shear_past的体积值
   ii = out('tension_past : ' + string(v_tension_past))   ;输出tension_past的体积值

end

get_plast                   ;执行函数


在论坛里找的帖子,但是运行以后 出现:cannot interpret function or array argument  这个错误,求帮忙解释下
离线chenpf09

发帖
334
土币
1142
威望
29
原创币
0
只看该作者 1楼 发表于: 2014-03-12
路过,抢个沙发!!!
离线wanglunan

发帖
328
土币
8783
威望
29
原创币
0
只看该作者 2楼 发表于: 2014-03-12
def get_plast   ;定义函数get_plast
  shearnow = 1                                     ;shearnow、tensionnow
  tensionnow = 2                                   ;shearpast、tensionpast
  shearpast = 4                                    ;的state value分别为1、2、4、8
  tensionpast = 8                                  ;具体内容参见fish reference 2.5.3.3

  v_shear_now = 0                           ;需求的shearnow、tensionnow
  v_tension_now = 0                         ;shearpast、tensionpast的体积,初始值为0
  v_shear_past = 0
  v_tension_past = 0   ;flac默认的值

  ;xx=z_state(find_zone(11779),0)   ;通过本步可验证

  p_z = zone_head                  ;单元指针
  loop while p_z # null            ;开始循环,如果单元指针不为空,执行以下操作

    if and(z_state(p_z,0),shearnow) = shearnow then ;判断语句,判断单元如果是shearnow状态
        v_shear_now = v_shear_now + z_volume(p_z)      ;则把该单元的体积加进v_shear_now
    endif
    if and(z_state(p_z,0),tensionnow) = tensionnow then  ;同上判断语句,判断单元是否为tensionnow状态
       v_tension_now = v_tension_now + z_volume(p_z)
    endif
    if and(z_state(p_z,0),shearpast) = shearpast then       ;同上,判断单元是否shearpast
       v_shear_past = v_shear_past + z_volume(p_z)
    endif
    if and(z_state(p_z,0),tensionpast) = tensionpast  then
       v_tension_past = v_tension_past + z_volume(p_z)
    endif
    p_z = z_next(p_z)          ;单元指针指向下一个单元
   endloop

   ii = out('shear_now : ' + string(v_shear_now))         ;输出shear_now的体积值
   ii = out('tension_now : ' + string(v_tension_now))     ;输出tension_now的体积值
   ii = out('shear_past : ' + string(v_shear_past))       ;输出shear_past的体积值
   ii = out('tension_past : ' + string(v_tension_past))   ;输出tension_past的体积值

end
get_plast                   ;执行函数
;;;运行fish后,执行pri fish就可以显示shearnow(当前剪切破坏),tensionnow(当前拉帐破坏)
;;;shearpast(过去剪切;破坏),tensionpast(过去拉帐破坏)的体积值
离线chenshunman

发帖
62
土币
6
威望
0
原创币
0
只看该作者 3楼 发表于: 2014-03-16
很不错资料
快速回复
限100 字节
温馨提示:欢迎交流讨论,请勿纯表情、纯引用!
 
上一个 下一个

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

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