有人编写过三维三角形常单元的程序吗？
我写了一个程序，但是有很多问题。我建了一个边长为1的立方体小块，在Z方向上的两个面分别施加Uz=0和Uz=1的约束。相当于拉伸。其他面上都为铰支。在计算过程中遇到很多问题。
（1）Pij*的奇异积分怎么算？我采用《工程科学中的边界元法》国防工业出版社1988中的方法计算。结果施加刚体位移时，居然能算出力来。
（2）于是我采用刚体位移原理计算奇异单元的H小矩阵（3x3）。解决了刚体位移的问题，但是结果还是不理想，随着单元密度增大，有的单元的应力也不断增大。
（3）我怀疑是普通单元积分的问题。但是检验了半天觉得采用的数值积分方法很简单，不会出错。
在程序的检验过程中，发现许多书上的三角形积分把积分点坐标和权系数都写错了。正确的位置和权系数应该是王勖成最新版的《有限元方法与数值计算》。
最后希望大家交流交流。 《matlab有限元分析与应用》的最后一章也许是你需要的！ 你到底是用有限元还是边界元? 楼主加我啊 qq10746268 我也在写三维三角形常单元的程序 希望共同交流 求写过三维三角形常单元的程序的高手十万火急 加我qq10746268
页:
[1]