六面体单元刚度矩阵，单点积分和沙漏控制原理，欢迎来讨论

mxlzhenzhu · 发表于 2013-8-13 23:52

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

本帖最后由 mxlzhenzhu 于 2013-8-14 00:24 编辑

由于计算结果和FEM对比不上，很容易就看到原来输入的原料就不一样，所以计算结果对不上。
排除了各种可能，包括：数值积分点的选取，形函数检查等等不下5、6个，还和matlab 的FEM toolbox对比了一下，计算结果跟他一模一样。

问题是什么？

哦，老王的书里头提到过缩减积分精度可能比“精确积分”要更精确，原因是什么？零能模式怎么理解？为什么加了约束可能会变好？

搜了一下，确信商业有限元软件里面不是用普通教科书里面的Gauss积分方案，而是某种变形====我的目前的认识。

附上标准六面体的计算结果的对比：

Nastran计算结果；
FEM toolbox计算结果；
我的计算结果

Three Methods.rar (36.72 KB, 下载次数: 9) ，其中E=74000,niu=0.3

再附上影响貌似很大的一篇论文，他讨论了两种阻尼方案，其中粘性阻尼不如弹性阻尼要好，对FEM感兴趣的朋友可以看看：

Hourglass control.part1.rar (500 KB, 下载次数: 5)

Hourglass control.part2.rar (500 KB, 下载次数: 6)

Hourglass control.part3.rar (421.13 KB, 下载次数: 5)

来此请教各位大虾对此有什么认识，也欢迎拍砖啊，不怕疼。

mxlzhenzhu · 发表于 2013-8-14 00:19

本帖最后由 mxlzhenzhu 于 2013-8-14 00:21 编辑

知网上有两篇：

动力显式算法有限元分析中的砂漏控制_解春雷.pdf (151.85 KB, 下载次数: 12)

&
变薄拉深过程模拟的有限元动力显式算法

superxjw · 发表于 2013-8-14 18:23

mxlzhenzhu 研究得很深入嘛！

账号		自动登录	找回密码
密码			我要加入

[弹性力学] 六面体单元刚度矩阵，单点积分和沙漏控制原理，欢迎来讨论

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。