handsomewww 发表于 2010-7-8 16:50:53

三维三角形常单元的 C系数 到底是多少？

最近我用matlab编写了一个求解三维弹性力学问题的边界元程序，用的是三角形常单元，其中Hu=Gp中的 H矩阵的主对角线系数的确定给我带来疑惑，书中都说系数是1/2*I加上节点对自身单元的解析积分之和, I为3x3的单位矩阵，但是我这这样做的时候 发现结果有一定误差，而我用刚体位移法，却得到了正确的结果，H对角线上的系数，也没有都等于1/2*I，只要三角形单元有一条边在几何模型的边上（曲率不连续的边，如立方体棱边）系数就有一个不是1/2，但是常单元，取节点为三角形形心，应该都算平滑的面，应该就是1/2，不明白 ，其中缘由 ，不知道有没有大牛指点一下。

handsomewww 发表于 2010-7-10 10:29:25

自己顶自己一下 这板块果然冷清

绿色的农民 发表于 2011-1-5 18:47:17

我帮顶啊 ，希望有高手帮你

handsomewww 发表于 2011-4-21 16:01:31

问题早已解决 程序编写有误
页: [1]
查看完整版本: 三维三角形常单元的 C系数 到底是多少？