h-层状自适应边界元方法的算法

摘要

给出了h-层状自适应边界元方法的基本算法,并介绍了应用这一算法编写的求解平面弹性静力学问题的计算机程序Elquahhbe.这个程序中包含一些新的内容:它的自适应细分过程通过一个新型的误差指示因子加以驱动;整个细分过程的信息都压缩到两组"指针"数组中,对"指针"数组进行简单运算即可得到生成离散网格和装配系统矩阵所需的信息,因此较好地解决了数据管理中存储与导出的矛盾.数值算例证明这些方法是行之有效的.