什么matlab和ansys计算的单元刚度矩阵不一致？

lorenxjtu

为什么matlab和ansys计算的单元刚度矩阵不一致？

ansys的solid187十节点四面体单元,单元刚度矩阵matlab和ansys的计算结果每个单元的差距在10%左右，基本如此。是什么原因造成的？难道是matlab程序的精度有问题？
我用的是体积坐标的积分，在积分这部分应该没有问题？

ansys help文件中的solid187单元形状函数是不是有问题

u=u1(2L1-1)L1+U2(2L2-1)L2+u3(2L3-1)L3+U4(2L4-1)L4

+4u5L1L2+U6L2L3+u7L1L3+u8L1L4+u9L2L4+u10L3L4

我怎么觉得这个式子有问题，ui各项系数之和不为1啊

是不是这样才对?：

u=u1(2L1-1)L1+U2(2L2-1)L2+u3(2L3-1)L3+U4(2L4-1)L4

+4(u5L1L2+U6L2L3+u7L1L3+u8L1L4+u9L2L4+u10L3L4)

pasuka

为什么matlab和ansys计算的单元刚度矩阵不一致？

ansys的solid187十节点四面体单元,单元刚度矩阵matlab和ansys的计算结果每个单元的差距在10%左右，基本如此。是什么原因造成的？难道是matlab程序的精度有问题？
...
lorenxjtu 发表于 2010-11-21 20:02

lz用的是几点高斯积分呢？照理2点高斯积分和3点高斯积分的结果肯定是不一样的

lorenxjtu

积分过程我是直接按照cook书上说明的，体积坐标按单元体积积分的公式计算的，书上介绍说积分变量的不同多项式按体积积分直接等于体积乘上不同系数即可。
这么做有什么不妥么？

pasuka

积分过程我是直接按照cook书上说明的，体积坐标按单元体积积分的公式计算的，书上介绍说积分变量的不同多项式按体积积分直接等于体积乘上不同系数即可。
这么做有什么不妥么？
lorenxjtu 发表于 2010-11-21 20:59

形函数第5-10项顺序问题
库克和王勖成的书上，顺序应该是
4L1L2
4L1L3
4L1L4
4L2L3
4L3L4
4L2L4

lorenxjtu

形函数第5-10项顺序问题
库克和王勖成的书上，顺序应该是
4L1L2
4L1L3
4L1L4
4L2L3
4L3L4
4L2L4
pasuka 发表于 2010-11-22 09:21

好，谢谢
我去查书，刚开始摸索着弄，好多东西不清楚
这个顺序是直接从ansys帮助文件上抄下来的

lorenxjtu

4# pasuka
多谢指教！

lorenxjtu

4# pasuka
查清楚了，ansys help与教材不符的原因是因为编号顺序不一致，排除这个因素，ansys help和教材中的顺序是一样的

pasuka

4# pasuka  
多谢指教！
lorenxjtu 发表于 2010-11-22 15:01

不客气！