浅谈沙漏问题及其解决方法
沙漏是中文对 Hourglass 的翻译，有限元里的 Hourglass 就是非物理的零能变形模式，或简而言之有变形没有应力或应变。
1. 什么是沙漏
在有限元分析中，一般以节点的位移作为基本变量，单元内节点的位移以及应变均采用形函数对各点位移进行插值计算得到。应力根据本构方程由应变计算得到，然后通过积分计算单元的内能。如果采用单点积分（积分点在等参元中心），在某些情况下节点位移不为零，即单元有形变，但插值得到的应变却为零。比如，一个正方体单元变形为等腰梯形，节点位移相等却方向相反，各点的形函数为零，所以插值结果为零，这样内能计算结果也为零，即计算认识单元没有变形，与事实相左。LS-DYNA 里的 Hourglass 现象基本上只发生在单点积分的单元里。比如，咱们随手画一个四个边的壳单元，只有一个积分点，这四个边只要围绕着这个积分点即可，具体这四个边是什么样子，是正方形，菱形，还是什么其他有创意的图形，积分点就管不了了，如果开始形状是长方形，在计算过程中变成了其他莫明奇妙的形状，这就很有可能是产生了 Hourglass。
2. 如何判断是否出现沙漏模式
方法一：查看单元变形过程：如果有单元变形明显异常，或有单元变成交替出现的梯形形状，一般是出现沙漏模式。
方法二：查看沙漏能在总内能中所占比例：当沙漏能约占总内能的 1% 时，表明沙漏模式对计算结果的影响不大；当其超过总内能的 10% 时，分析就是无效的，必须采取措施加以解决。
3. 控制沙漏的方法
Hourglass 可通过引入内部节点力来控制，相关的关键字有 *CONTROL_Hourglass, *Hourglass 等，但是这就有了另一个问题，既然是人为的引用了节点力，那么这个人为的力所产生的能量（Hourglass energy) 就得越小越好，一般要求沙漏能要小于总能量的 5% 才认为结果是可靠的。
- 分享于 · 2019.09.16 17:09 · 阅读 · 13131
[版权声明] ：本文系网友分享，仅以非商业性的交流和科研为目的，并不代表本网赞同其观点和对其真实性负责。如涉及作品内容、版权和其它问题，请在30日内与本网（ friends@stuch.cn ）联系！我们将协调给予处理。转载授权码：确权中，请原作者尽快与我们取得联系，阅读原文(请登录)..
相关推荐:
- 文章 LS DYNA 共享式计算（SMP）和分布式计算（MPP）计算效率对比 11889 浏览
- 文章 Consistent system of units 单位制度（LS DYNA) 2909 浏览
- 文章 比较有限元中复杂结构的网格划分方法 4349 浏览
- 文章 有限元软件发展趋势分析 3376 浏览
- 文章 LS DYNA 二进制计算结果文件 d3plot 通用文件格式转换器 6354 浏览
- 问答 ANSYS/LS-DYNA 中怎么给模型施加位移荷载呢？ 10219 浏览
- 问答 无法找到 LS971 R5.1.1 32 位安装包，在做二次开发，目前只有 64 位安装包 2890 浏览
- 问答 流固耦合方式速度和加速度耦合（2）和罚函数（4、5）怎么选取 4884 浏览
- 问答 您好，请问一下我再用 LS-DYNA R7.0.0 版本二次开发的时候，利用 Fortran 编译器 nmake 后一直出现无法读取 libguide.lib 的情况，请问是不是版本的问题（VS2008+IVF2011），多谢多谢 3622 浏览
- 问答 使用 vendor 关键字加密 LS DYNA 的 K 文件时，提示错误 input conversion error 5954 浏览
- 问答 数值模拟对单个单元进行双轴压缩模拟，加载方式应该用力加载还是用位移加载 3908 浏览
