我有薄“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] 的有限元组件并获得全局刚度矩阵，结果作为节点位移矢量和节点力向量？

提示表示赞赏。谢谢。

关于工程堆栈交换的相同问题

关于stackoverflow的同样问题

编辑：如果有人知道这个用 Java、C++ 甚至 VBA 编写的非常具体的程序的开源免费代码或实现，那么这也将是问题的答案！

编辑：我在这里找到了一个应该可以正常工作的代码。将这段代码翻译成另一种过程编程语言有多难？涉及多少痛苦？