FLAC3D500 只能提取全部单元体应力命令_猪冰龙

http://blog.sina.com.cn/u/2379298071

首页博文目录关于我

个人资料

微博

加好友发纸条

写留言加关注

博客等级：
博客积分：

博客访问：
关注人气：
获赠金笔：0支
赠出金笔：0支
荣誉徽章：

正文字体大小：大中小

FLAC3D500 只能提取全部单元体应力命令

(2015-04-16 23:46:54)

标签：

flac3d500

stress

zone

分类： flac3d

 C++ Code 

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21

def GET_SZZ_
array buf4(1)
buf4(1) = 'x  y  z\n'
file='1fish-resul.txt'
status=open(file,1,1)
p_z = zone_head
loop while p_z # null
buf4(1) =' '
     if z_zcen(p_z)>=4  then
      if z_zcen(p_z)<=5  then
buf4(1) = buf4(1) + string(z_xcen(p_z)) + ' '
buf4(1) = buf4(1) + string(z_ycen(p_z)) + ' '
buf4(1) = buf4(1) + string(z_szz(p_z)*(-0.000001))
      endif
     endif
p_z = z_next(p_z)
  status=write(buf4,1)
endloop
status=close
end
@GET_SZZ_

阅读┊ 收藏 ┊ 喜欢 ▼ ┊打印┊举报/Report

前一篇：flac3d500 仅显示塑性区

后一篇：郑海棠老师Flac3D2Tecplot_new.dat代码及个人理解

新浪BLOG意见反馈留言板　欢迎批评指正