矩形薄板Kirchhoff板的单元刚度矩阵中的整体刚度矩阵

问题描述

我有一块薄的“基尔霍夫”板的单元刚度矩阵。 平板为3 [m] x 5 [m],仅在所有边缘上支撑。其厚度为0.2 [m]。在板上施加恒定的载荷p = 5 [kN /m²]。当mu = 0,2和E = 30000 [MPa]时,中点位移等于0,168 [mm]。

最有效或最快的方法是什么,以及如何实现或编码此板的0.5 [m] x 0.5 [m]的有限元装配并获得整体刚度矩阵,其结果为节点位移矢量和节点力矢量?

提示表示赞赏。谢谢。

3 DOFs per Node. 4 Nodes. 12 x 12 Element stiffness matrix.

enter image description here

enter image description here

enter image description here

enter image description here

Same question on Engineering stackexchange

Same question on Computational Science stackexchange

编辑:如果有人知道用Java,C ++甚至VBA编写的非常特定的程序的开源免费代码或实现,那么这也将是该问题的答案!

编辑:我发现here应该可以正常工作的代码。将此代码转换为另一种过程编程语言有多困难?有多少痛苦?

解决方法

暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!

如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@)