二维各向异性位势薄体问题的虚边界元法

传统边界元法分析各向异性薄体问题时涉及奇异边界积分和拟奇异边界积分的处理,估计这些积分具有相当的难度而且耗时.提出了求解二维各向异性位势薄体问题的虚边界元方法,给出了求解此类问题的新途径,同时拓展了虚边界元法的应用范围.数值算例表明,虚边界元法可有效求解二维各向异性位势薄体问题,且方法简单、精度高、易于程序设计.     

二维弹性薄体问题的虚边界元法

拓展了虚边界元方法的应用范围,将其应用于二维弹性薄体问题,避免了奇异边界积分和几乎奇异边界积分的计算.通过数值算例验证了虚、实边界的距离公式,公式的特点是距离与边界离散单元数有关,表明公式对于二维薄体结构同样适用.按照张耀明等虚边界元法的理论分析公式选择虚、实边界间的距离,即使结构狭窄到纳米级(10-9 m),依然可获得高精度的数值解.     

非定常扩散方程的虚边界元解法

对于二维齐次和非齐次非定常扩散方程问题,利用与时间有关的基本解,基于单层位势的延拓,建立虚边界积分方程,然后用虚边界元法求解.最后,给出了数值算例验证了虚边界元法求解非定常扩散方程问题的可行性和有效性.     

多域组合非耦合热弹性问题的虚边界元法

采用边界元法求解热弹性力学问题通常涉及到关于温度作用的域内体积分,使其在求解此类问题时失去了可降维的优点.为此,应用虚边界元法思想分别考虑热传导问题和与之对应的弹性力学问题的数值格式,并将两者的求解思路结合起来,从而形成解多域组合非耦合热弹性问题时无需计算域内体积分的虚边界元法思想.该方法具有一般性,既适用于二维问题又适用于三维问题,而且可将多域求解思想蜕化到单域问题.按单域定义的方板、厚壁圆筒热应力的计算和按多域定义的含圆形夹杂方板有效热膨胀系数的数值模拟结果已充分表明该方法具有较好的计算效率和较高的计算精度.     

二维位势问题快速多极虚边界元解

快速多极虚边界元法是近期发展起来的一种数值算法;其对大规模复杂问题的计算,能在保证求解精度的前提下,使计算量和存储量均比常规虚边界元法具有在数量级上的减少.本文给出了二维位势问题快速多极虚边界元法的求解思想,并进行了数值论证;由文中的数值结果可知,本文方法具有可行性,且有较好的计算精度.     

虚边界元法在二维涂层结构温度场中的应用

将虚边界元法应用于平面涂层结构温度场问题,并发展了多域虚边界元法。给出求解涂层问题的新的途径,同时也拓展了虚边界元法的应用范围。对狭长比为10-1～10-10的涂层结构进行了研究,所取得的数值结果与精确解高度吻合,表明虚边界元法是求解二维涂层结构温度场问题的强有力工具,且方法简单、易于程序设计。     

正交各向异性不同材料组合体的虚边界元解法

依据多域组合问题虚边界元法思想,采用适用于正交各向异性介质和各向同性介质平面弹性问题基本解的统一模式,提出了虚边界元法求解正交各向异性弹性体与不同材料性质弹性体的组合问题的数值算法思想.文中给出了带孔的正交各向异性板和正交各向异性材料与各向同性材料结合体的数值算例.数值结果表明,该方法具有较高的计算精度和较好的计算效率.提出的数值思想具有较好的通用性,其不但能求解正交各向异性材料的多域组合问题及正交各向异性材料与各向同性材料的结合体问题,而且也能蜕化求解各向同性材料的多域组合问题.     

样条虚边界元法的数值稳定性与误差估计

样条虚边界元法是针对传统间接奇异边界元法存在的问题而提出的一种半解析半数值方法。它既保留了边界元法的优点，也避开了求解奇异积分方程的问题，在试函数和权函数的选取方面也作出了改进，具有精度好、效率高等优点。本文主要针对弹性力学平面问题样条虚边界元法在数值稳定性与误差估计方面的问题展开讨论，获得了虚边界的布设规律及方法误差的直观度量，为该法的实际应用打下了更好的基础。     

场点与载荷点分离时的边界元法误差分析

本文推导了场点与载荷点分离时的边界积分方程。误差分析的结果证实边界元系统方程中,各组元的误差和模的大小对数值求解精度起决定作用。以分析场点与载荷点的距离对方程中组元误差和的大小的影响,可得出非奇异间接边界元法的求解精度比奇异边界元法的求解精度要高得多,间接边界元法的数值稳定性好于直接边界元法的数值稳定性等有益的结论。这种分析方法能够定性的解释文中的计算结果。     

各向异性板应力集中问题最小二乘虚边界元解法

针对各向异性板的应力集中问题,依据虚边界元法的求解思路,以复变函数表达的基本解作为权函数,建立了相应最小二乘虚边界元的数学模式;其可求解正交各向异性或一般各向异性材料的平面问题.文中给出了含圆孔的各向异性板应力集中问题的数值算例;通过与边界元直接法、有限元法的数值比较可知,本文方法的数值结果具有较高的计算精度.此外,相对其它数值方法本文方法对于各向异性板应力集中问题的求解,具有较好的适用性和数值计算的稳定性.     

随机分布圆孔板有效弹性模量快速多极虚边界元法模拟

将快速多极算法和广义极小残值法(GMRES)结合于虚边界元法的方程求解,形成了快速多极虚边界元法的求解思想.本方法采用了"源点"多极展开和"场点"局部展开的组合处理方案,使得原问题方程组求解的计算耗时量和储存量均降至与所求问题的计算自由度数成线性比例.文中分析了含随机分布多圆孔板的有效弹性模量,并与其它数值方法的结果进行了比较,同时数值验证了本方法的可行性、计算精度及计算效率.     

三维快速多极虚边界元配点法

以三维弹性力学问题为研究背景,提出了一种三维快速多极虚边界元配点法的求解思想,即将三维快速多极展开的基本思想和广义极小残值法运用于求解传统虚边界元配点法方程.文中将三维弹性问题的基本解推导为适合于虚边界元快速多极算法的展开格式,经数值计算格式的演变,使求解方程的计算量和储存量与所求问题的计算自由度数成线性比例,以达到数值模拟大规模自由度问题的目的.算例说明了该方法的可行性、计算效率和计算精度.此外,该方法的思想具有一般性,应用上具有扩展性.     

三维Helmholtz方程外边值问题的虚边界元法

基于位势的延拓,推导出三维虚边界积分方程.通过选择不同的虚边界,避免相应内问题的特征值与波数重合,从而保证解的唯一性.数值算例验证了该方法求解任意波数三维Helmholtz方程外边值问题的有效性.     

厚壳三维分析的虚边界元最小二乘法

给出虚边界元最小二乘法的基本思想，并计算了厚壳问题，与边界元直接法相比，避免了奇异积分的数值处理，且系数阵是对称的，程序实现较容易，节省近一半机子内存。由数值算例表明，计算精度是令人满意的。     

无单元法在平面热应力分析中的应用

采用无单元法对平面热应力问题进行了研究。首先从热应力边值问题出发,给出了与其等价的能量泛函的弱变分形式,并导出了无单元方法的计算公式。在此基础上,求解了平面稳态温度场的第一类边值条件的温度分布。并对其进行了热应力分析,取得了较好的数值结果。