定义单元体能量为微单元体的应变比能,即当应力和应变满足线性关系时,微单元体在三向应力状态下的应变比能为:

(3.1)

下面代码为用户自定义云图显示变量。

Flac3d Code

 1 new
 2 gen zone brick size 10 10 10 group lay-1
 3 mo el
 4 prop bulk 3e8 shear 1e8
 5 ini den 20e3
 6 set gravity 0 0 -9.8
 7 fix z range z 0
 8 solve
 9 config zextra 2
10 def _udf_zextra
11  p_z = zone_head
12  _u=0.3
13  _E=2.7e8 ;_E=9*bulk*shear/(3*bulk+shear)
14 loop while p_z # null
15    if z_group(p_z)='lay-1';只计算在组内的单元体
16     _sig1=z_sig3(p_z)
17     _sig2=z_sig2(p_z)
18     _sig3=z_sig1(p_z)
19     _a=_sig1^2+_sig2^2+_sig3^2
20     _b=2*_u*(_sig1*_sig2+_sig2*_sig3+_sig3*_sig1)
21     z_extra(p_z,1) =(_a-_b)/(2*_E);微单元体应变比能
22     z_extra(p_z,2) = z_sxx(p_z);自定义为单元体sxx应力
23    endif
24    p_z = z_next(p_z)
25 endloop
26 end
27 _udf_zextra
28 plot add axes
29 plot set plane nor 0 1 0 ori 0 5 0
30 plot con zextra 1 plane ;range 5 7,不能加range
31 ;plot add con zextra 2 plane

第二个变量取为sxx只是为了显示在一段代码中如何同时自定义两个变量。多个变量与此类似,不再祥述。

   

a Szz应力云图

a The szz stress cloud

b 应变比能应力云图

b Strain specific stress cloud

图3.3 自定义单元应力有关云图

Fig.3.3 The custom cloud map with zone stress

下面是节点变量:

Flac3d Code

 1 new
 2 gen zone brick size 10 10 10 group lay-1
 3 mo el
 4 prop bulk 3e8 shear 1e8
 5 ini den 20e3
 6 set gravity 0 0 -9.8
 7 fix z range z 0
 8 solve
 9 config gpextra 2
10 def _udf_gpextra
11  p_gp = gp_head
12 loop while p_gp # null
13  gp_extra(p_gp,1) = sqrt((gp_xpos(p_gp)-xc)^ 2+(gp_zpos(p_gp)-zc)^ 2)
14  gp_extra(p_gp,2) = atan2((gp_xpos(p_gp)-xc),(zc-gp_zpos(p_gp)))
15  p_gp = gp_next(p_gp)
16 endloop
17 end
18 set xc=5 zc=5
19 _udf_gpextra
20 plot set background white
21 ;plot block gr
22 plot con gpextra 1
23 plot con gpextra 2

   

a 自定义节点云图1

a custom node cloud 1

b 自定义节点云图2

b custom node cloud 2

图3.4 自定义节点运动量有关云图

Fig.3.4 The custom cloud map with gridpoint movement

想看详细的.sav和fish代码请关注下面的二维码(flac3d自定义变量输出云图显示)。

转载于:https://www.cnblogs.com/zhubinglong/p/6863851.html

flac3d自定义变量输出云图相关推荐

  1. awk 内置变量与自定义变量

    点击上方"生信科技爱好者",马上关注真爱,请置顶或星标 作者:ghostwu 原文:https://www.cnblogs.com/ghostwu/p/9085653.html A ...

  2. 自定义变量 配置文件_「系统架构」Nginx调优之变量的使用(3)

    在上一篇文章「系统架构」Nginx调优之变量的使用(2)中我们介绍了自定义变量和内置变量,下面我们继续接着介绍Nginx中变量的可见性和动态内置变量. 变量的可见性 nginx中的变量虽然不全是全局变 ...

  3. e0312 不存在用户定义的_VistaPro创建自定义变量

    IESVE软件中VistaPro模块中的结果变量众多,很多用户反应变量数据不容易获得,其实在此模块中用户可以根据现有变量(天气数据.房间设置和负荷.建筑负荷.能源.碳排放)进行组合,利用VE软件自带的 ...

  4. 打开高效文本编辑之门_Linux awk之自定义变量与操作符

    Linux awk之自定义变量与操作符 声明与简介 AWK 是一种处理文本文件的语言,是一个强大的文本分析工具.awk通常用来处理结构化(固定格式)的文本文件, awk每接收文件的一行,然后执行相应的 ...

  5. 打开高效文本编辑之门_Linux Awk自定义变量与操作符

    Linux awk之变量与操作符应用案例 声明与简介 AWK 是一种处理文本文件的语言,是一个强大的文本分析工具.awk通常用来处理结构化(固定格式)的文本文件, awk每接收文件的一行,然后执行相应 ...

  6. 文件自定义变量_awk 内置变量与自定义变量

    点击上方"生信科技爱好者",马上关注真爱,请置顶或星标 作者:ghostwu 原文:https://www.cnblogs.com/ghostwu/p/9085653.html A ...

  7. MATLAB代码编写中自定义变量的问题

    1. 吃过亏,上过当 自己在进行数据处理的过程中,遇到过很多次自己编写的变量存在问题的情况,到目前为止,遇到过的情况如下: 情况 问题 详细描述 解决办法 1 自定义变量和程序中的函数名称相同 比如: ...

  8. [CSS] 自定义变量带你随心所欲,一键换肤

    认识CSS自定义变量 CSS自定义属性?听着怎么那么神奇呢,属性还可以自定义,那不是可以放肆地玩耍?我自己定义的属性浏览器都能认识? 一.基础 (一)名字和用途 其实CSS自定义属性还有很多小名,比如 ...

  9. Shell脚本,环境变量,位置变量,预定义变量,自定义变量,单引号双引号反引号用法,一般变量和交互式变量read定义,密文输入变量信息

    Shell脚本,变量,系统预设的环境变量,位置变量,预定义变量,用户的自定义变量 一.定义变量 1.以固定的名称存放可以变化的值 2.格式为 变量名=变量值 3.取消变量的格式 unset 变量名 4 ...

最新文章

  1. 美多商城之验证码(异步方案)
  2. 微软笔试题 2013暑期实习笔试题目
  3. scrapy-splash抓取动态数据例子六
  4. WebStorm 2019.3.1安装教程
  5. 配置ORACLE 客户端连接到数据库
  6. CCNA笔记之第十九节:RIP协议(大实验3)
  7. C#将一个excel工作表根据指定范围拆分为多个excel文件
  8. 基于QGIS初探PostgreSQL的PostGIS插件,包括YUM和编译安装PostGIS
  9. 交错排列(Alternating Permutation)问题详解
  10. 生产线改造升级 蔚来部分10月份车辆订单将推迟交付
  11. 【Kafka】Kafka 修改某个消费组的偏移量
  12. java8 foreach 伟参_【java8】为java8的foreach正名
  13. 阿里大佬浅谈大型项目前端架构设计
  14. 通信原理及系统系列19—— 锁相环(鉴相器分析_2)
  15. sql语句的增删改操作附加代码
  16. 阿里云推高校YY计划助推大学生创业潮
  17. 前端入门-HTML篇
  18. 【Unity】TimeLine常见问题:如何实现人形角色动画的平滑切换
  19. svg + image/svg+xml 文件上传实现xss
  20. Python matplotlib 批量绘图内存不够问题(Out of memory)

热门文章

  1. w10本地安装RocketMq(插件)以及整合springboot
  2. SAP UI5 应用开发教程之一百零一 - SAP UI5 应用的 Locale 决定机制试读版
  3. MIT 6.S081 实验5 笔记与心得
  4. 2018可信区块链峰会丨《区块链溯源白皮书(版本1.0)》重磅发布
  5. 树莓派编译jssc-2.9.2 实现自定义波特率
  6. UUCTF2022-Writeup
  7. rstp 小米网络摄像头_常见网络摄像机的端口及RTSP地址
  8. 标准工业控制编程语言IEC61131-3简述
  9. 一次关于缓存行的切身体验
  10. Win10系统下怎么开HDR?