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

flac3d 中pile单元弯矩输出的fish程序问题 [复制链接]

上一主题 下一主题
离线blueprint27
 

发帖
11
土币
3
威望
6
原创币
0
只看楼主 倒序阅读 使用道具 楼主  发表于: 2015-07-08

def force_outp
jj=0
loop cid (1,10)
sp=s_find(cid)
   if s_type(sp)=1 then
      jj=jj+1
;;;第一个节点相关数据
      np=s_node(sp,1)
      xx1=nd_pos(np,1,1)
      yy1=nd_pos(np,1,2)
      zz1=nd_pos(np,1,3)
      ffx1=sp_force(sp,1,1)
      ffy1=sp_force(sp,1,2)
      ffz1=sp_force(sp,1,3)
      mmx1 =sp_mom(sp,1,1)
      mmy1 =sp_mom(sp,1,2)
      mmz1 =sp_mom(sp,1,3)
;;;第二个节点相关数据
      np=s_node(sp,2)
      xx1=nd_pos(np,2,1)
      yy1=nd_pos(np,2,2)
      zz1=nd_pos(np,2,3)
      ffx1=sp_force(sp,2,1)
      ffy1=sp_force(sp,2,2)
      ffz1=sp_force(sp,2,3)
      mmx1=sp_mom(sp,2,1)
      mmy1=sp_mom(sp,2,2)
      mmz1=sp_mom(sp,2,3)
;;;形心相关数据
      xx=(xx1+xx2)/2
      yy=(yy1+yy2)/2
      zz=(zz1+zz2)/2
      ffx=(ffx1+ffx2)/2
      ffy=(ffy1+ffy2)/2
      ffz=(ffz1+ffz2)/2
      mmx=(mmx1+mmx2)/2
      mmy=(mmy1+mmy2)/2
      mmz=(mmz1+mmz2)/2
     xtable(1,jj)=(1,jj)
     ytable(1,jj)=mmy
   endif
endloop
end
force_outp
set logfile pile_moment.txt
set log on
set pagelength 1000
pri table 1
请高手指导,每次都提示 table 1 not found ,意思就是表1不存在,弯矩没有输入到表中,请问为什么?命令错误在哪里呢?谢谢!
离线龙biscuit

发帖
322
土币
2507
威望
45
原创币
0
只看该作者 1楼 发表于: 2015-07-18
xtable(1,jj)=(1,jj)
右边应该是一个数吧
邮箱:blzhu@buaa.edu.cn ,欢迎交流!
离线blueprint27

发帖
11
土币
3
威望
6
原创币
0
只看该作者 2楼 发表于: 2015-07-18
回 龙biscuit 的帖子
龙biscuit:xtable(1,jj)=(1,jj)
右边应该是一个数吧 (2015-07-18 06:49) 

谢谢您的回答,这个问题我花了一下午的时间看教程解决了!问题是在pile单元的判断和后面命令有冲突上!请问你能不能看见结构单元模拟后,应力的影响啊?
离线307000063

发帖
7
土币
6
威望
0
原创币
0
只看该作者 3楼 发表于: 2016-05-31
if s_type(sp)=3
离线307000063

发帖
7
土币
6
威望
0
原创币
0
只看该作者 4楼 发表于: 2016-06-14
弯矩和剪力均可以由这个程序提取
快速回复
限100 字节
温馨提示:欢迎交流讨论,请勿纯表情、纯引用!
 
上一个 下一个

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

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