六面体单元刚度矩阵，单点积分和沙漏控制原理，欢迎来讨论
本帖最后由 mxlzhenzhu 于 2013-8-14 00:24 编辑由于计算结果和FEM对比不上，很容易就看到原来输入的原料就不一样，所以计算结果对不上。
排除了各种可能，包括：数值积分点的选取，形函数检查等等不下5、6个，还和matlab 的FEM toolbox对比了一下，计算结果跟他一模一样。
问题是什么？
哦，老王的书里头提到过缩减积分精度可能比“精确积分”要更精确，原因是什么？零能模式怎么理解？为什么加了约束可能会变好？
搜了一下，确信商业有限元软件里面不是用普通教科书里面的Gauss积分方案，而是某种变形====我的目前的认识。
附上标准六面体的计算结果的对比：
Nastran计算结果；
FEM toolbox计算结果；
我的计算结果
，其中E=74000,niu=0.3
再附上影响貌似很大的一篇论文，他讨论了两种阻尼方案，其中粘性阻尼不如弹性阻尼要好，对FEM感兴趣的朋友可以看看：
来此请教各位大虾对此有什么认识，也欢迎拍砖啊，不怕疼。
本帖最后由 mxlzhenzhu 于 2013-8-14 00:21 编辑
知网上有两篇：
&
变薄拉深过程模拟的有限元动力显式算法
mxlzhenzhu 研究得很深入嘛！
页:
[1]