用Tchebychev多项式逼近边界元法中的基本解

正交各向异性板动力问题的边界元法之难点在于找不到其封闭形式的基本解。因此。寻找一种适用性强,精度高、使用方便的求近似基本解的方法,在边界元法的应用中具有重要意义。本文用Tchebychev多项式组成广义双Fourier级数逼近问题的基本解,由算例表明该方法是令人满意的。     

边界元理论在复杂外边界油藏水平井渗流中的应用

受构造作用的影响,实际油藏的外边界往往是复杂多样的.本文从渗流理论出发建立了复杂外边界油藏水平井渗流数学模型,并采用Lord Kelvin点源解、贝塞尔函数积分和泊松叠加公式等方法求解了复杂外边界油藏水平井的边界元基本解,利用边界元的理论建立了复杂外边界油藏水平井井底压力响应数学模型.通过计算得到了无因次压力和压力导数双对数理论图版,并在其基础上分析了复杂外边界油藏水平井渗流特征及其影响因素.     

正交各向异性弹性问题的规则化边界元法

论文致力于平面正交各向异性弹性问题的规则化边界元法研究,提出了新的规则化边界元法的理论和方法.对问题的基本解的特性进行了研究,确立基本解的积分恒等式,提出一种基本解的分解技术,在此基础上,结合转化域积分方程为边界积分方程的极限定理,建立了新颖的规则化边界积分方程.和现有方法比,论文不必将问题变换为各向同性的去处理,从而不含反演运算,也有别于Galerkin方法,无需计算重积分,因此所提方法不仅效率高,而且程序设计简单.特别是,所建方程可计算任何边界位移梯度,进而可计算任意边界应力,而不仅限于面力.数值实施时,采用二次单元和椭圆弧精确单元来描述边界几何,使用不连续插值逼近边界函数.数值算例表明,论文算法稳定、效率高,所取得的边界量数值结果与精确解相当接近.     

薄壳问题的三维虚边界元解法

直接从三维弹性力学微分方程出发,依据三维的Kelvin解,应用最小二乘法建立了三维虚边界元法解薄壳问题的一般方法。本方法的显著优点是：不论求解何种壳体问题,方法的思想是不变的,均以三维的Kelvin解来建立方程,而勿需对不同几何形状的壳体采用不同的基本解。文中给出了数值算例,以作为本方法的应用。本文方法与边界元直接法相比,优点在于无需处理奇异积分,且系数阵是对称的;再者,本文方法思想简单,程序实现容易。     

杆件结构强迫振动分析时域基本解

基本解是边界元法、格林函数法、基本解法等数值方法必备的理论基础之一.大多问题的基本解已经被推导出来,但杆件结构强迫振动分析的时域基本解未见相关报道.现有方法在求解强迫振动问题时通常采用一定技巧将问题等效为静力问题,并利用静力基本解进行求解,这就增加了列式过程的复杂程度和方法的计算工作量.论文通过振型分解、杜哈梅积分等步骤推导出杆件结构强迫振动分析的时域基本解,该基本解可以直接用于结构的动力响应分析,提高运算效率.论文工作可以为边界元法等数值方法在进行杆系结构强迫振动分析时提供新的思路.     

薄壳问题的一种数值解法

本文直接从三维弹性力学微分方程出发,依据三维的Kelvin解,应用最小二乘法建立了三维虚边界元法解薄壳问题的一般方法。本方法的显著优点是：不论求解何种壳体问题,思想是不变的,均以三维的Kelvin解来建立方程,而勿需对不同几何形状的壳本采用不同的基本解。文中给出了数值算例,以作为本方法的应用。本文方法与边界元直接法相比,优点在于无需处理奇异积分,且系数阵是对称的;再者,本文方法思想简单,程序实现容易。     

弹性半空间动力问题的基本解

近年来,以Lamb解为基本解用边界元法求解了大量的弹性半空间表面基础的动力响应问题,但是,Lamb解却无助于用边界元法处理弹性半空间的内部问题。本文根据文献[1]所提出的弹性动力学的完备解,应用叠加原理导出了弹性半空间在内部集中动载作用下的基本解。从而为用边界元方法求解基础埋入和地下工程等弹性半空间内部问题开拓了一条新的途径。1.弹性半空间在非轴对称表面动载作用下的一般解当弹性半空间仅受到零初始条件的非轴对称表面动载作用时,总可以将表面载荷     

回转体弹塑性扭转问题的一种边界积分方程解法

本文利用调和函数性质和格林公式,建立以应力函数和边界合剪应力为场变量的边界积分方程.对奇异积分,按离开奇异点远近不同分别采用不同阶次的高斯积分.处理了角点处五种可能的情况和零轴上积分分母为零的现象.对于弹性和弹塑性问题都给出了部分实例.并与理论解,由Kelvin解建立的边界元法和松弛法作了比较.均得出十分满意的结果.     

电磁弹性固体三维问题的虚边界元-等额配点法

依据弹性力学虚边界元法的基本思想和电磁弹性固体的基本解,提出了电磁弹性固体三维问题的虚边界元一等额配点法.该方法继承传统边界元法优点的同时,有效地避免了传统边界元法的边界积分奇异性的问题.算例表明该方法有很高的精度,是求解电磁弹性固体三维问题的一个有效的数值方法.     

基于Erdogan基本解边界元法计算应力强度因子

引入含裂纹问题基本解(Erdogan基本解),提出了基于Erdogan基本解的样条虚边界 元法,并阐述了该法在实施过程中的特点与具体做法. 采用该方法详细分析了若干 典型裂纹问题,全面考察了方法的计算精度和收敛情况,以及在求解复杂裂纹问题方面 的能力. 结果显示,该方法具有精度高、收敛快、计算能力强等优点,是裂纹问题分析中 一种具有竞争力的通用计算方法.     

弹塑性杂交/混合有限元法与实施

本文提出两种结构弹塑性分析的有限元方法——杂交/混合非协调元法。该法采用场变量分解,文[1]对增量形式的杂交应力元能量泛函进行修正,简化了高阶矩阵的求逆运算,从而提高了多变量非线性有限元分析的计算效率和收敛精度。文中按文[2]给出的第二种能量泛函的修正形式建立了弹塑性平面问题中的杂交/混合非协调四边形单元。最后给出算例,并与实验解及文[3]解进行比较。表明该方法分解弹塑性问题是十分有效的。     

非线性大变形问题边界元分析的自校正方法

针对用增量法求解非线性方程解的漂移问题，在非线性问题边界元法计算中建立了自我校正方法，对在拖带坐标上建立的增量形式的基本方程，引入Ｌａｎｇｒａｎｇｅ校正因子，以全量形式的基本方程作为其辅助方程，在此基础上导出含校正项的边界积分方程，边界元自我校正方法的建立有效地保证了在非线性问题的计算中最终收敛在其解附近，提高了计算精度和运算效率。     

三维弹性问题Taylor展开多极边界元法的误差分析

Taylor展开多极边界元法有效的提高了边界元法的求解效率,使之可用于大规模问题的计算。然而,由于计算中对基本解进行了Taylor级数展开,与传统边界元方法相比计算精度有所下降。本文主要针对三维弹性问题Taylor展开多极边界元法的计算精度和误差进行研究。文中对两种方法的计算精度进行了比较;研究了核函数的Taylor展开性质;推导了三维弹性问题基本解的误差估计公式;给出了Taylor展开多极边界元法中远近场的划分原则。通过具体的算例,证明了该方法的正确性和误差估计公式的有效性,说明了影响Taylor展开多极边界元法求解精度的因素。     

基于奇异分解和自适应BEM积分算法的水下航行体机械噪声预报

提出一种改进的声学边界元法（M-BEM）用于准确计算水下航行体发动机振动引起的近场辐射噪声。分别采用奇异分解技术和自适应边界元积分算法解决了Helmholtz积分方程在求解近场声压时出现的超奇异积分和奇异积分问题。采用一脉动球源的声辐射算例对方法进行验证,数值解与精确解误差小于1．5dB。结合有限元方法并考虑流固耦合作...     

一种处理弹性动力边界元法中奇异积分的方法

利用边界元法求解瞬态弹性动力学问题时,时域基本解函数的分段连续性和奇异性为该问题的求解带来很大的困难。为了解决时域基本解中的奇异性问题,本文依据柯西主值的定义,对经过时间解析积分之后的时域基本解进行奇异值分解,将其分成奇异和正则积分两部分;其中正则部分可通过采用常规高斯积分方法来计算,而奇异部分具有简单的形式,可以利用解析积分计算。经过上述操作之后,就可以达到直接消除时域基本解中奇异积分的目的。和传统方法相比,本文方法并不依赖静力学基本解来消除奇异性,是一种直接求解方法。最后给定两个数值算例来验证本文提出方法的正确性和可行性,结果表明使用本文算法可以解决弹性动力学边界积分方程中的奇异性问题。     

A HIGH ORDER KERNEL INDEPENDENT FAST MULTIPOLE BOUNDARY ELEMENT METHOD FOR ELASTODYNAMICS

基于核无关的快速多极方法, 发展了一种弹性动力学问题的快速、高精度边界元分析方法. 采用基于二次曲面单元的Nyström 离散, 将边界积分方程转化为求和形式, 可以方便地进行加速计算;由于采用二次元, 边界元分析精度很高. 将一种新型快速多极方法用于Nyström 边界元法的加速计算, 该方法的数值实现简便、不依赖于积分方程基本解的表达式, 因此通用性很好;该方法还具有最优的计算量和存储量、精度高且可以控制. 结合Nyström 边界元系数矩阵和快速多极方法转换矩阵的特点, 提出一种大幅度降低边界元内存消耗的策略. 数值结果表明, 该方法无论在分析精度, 还是计算速度和内存消耗上, 都大大优于同类方法, 是一种快速、通用的工程弹性动力学问题大规模数值分析方法.     

弹性动力学高阶核无关快速多极边界元法

基于核无关的快速多极方法, 发展了一种弹性动力学问题的快速、高精度边界元分析方法. 采用基于二次曲面单元的Nystr?m 离散, 将边界积分方程转化为求和形式, 可以方便地进行加速计算;由于采用二次元, 边界元分析精度很高. 将一种新型快速多极方法用于Nystr?m 边界元法的加速计算, 该方法的数值实现简便、不依赖于积分方程基本解的表达式, 因此通用性很好;该方法还具有最优的计算量和存储量、精度高且可以控制. 结合Nystr?m 边界元系数矩阵和快速多极方法转换矩阵的特点, 提出一种大幅度降低边界元内存消耗的策略. 数值结果表明, 该方法无论在分析精度, 还是计算速度和内存消耗上, 都大大优于同类方法, 是一种快速、通用的工程弹性动力学问题大规模数值分析方法.      

一种改进的边界元法在弹性扭转问题中的应用

本文用一种改进的边界元法分析与计算了椭圆截面等直杆的扭转问题,并与正规的边界元法的解进行比较,其结果完全一致.然而,改进边界元法较正规边界元法需要准备的数据大大减少,计算时间更加缩短.因此,本文方法对求解 Poisson 方程问题是一种经济而行之有效的数值计算方法.     

弹性力学柯西问题奇异值选取的傅里叶系数法

截断奇异值分解和边界元法用于重建二维线弹性力学问题的边界条件.提出了傅里叶系数法得到奇异值截断数.将傅里叶系数法与L-曲线法和偏差原理进行了对比,发现傅里叶系数法是得到奇异值截断数的一种有效的方法.对弹性力学Cauchy问题,数值结果表明,随着输入已知边界条件数据随机偏差的减小,基于傅里叶系数法的截断奇异值分解方法可以得到稳定精确的解.     

动力边界元法的正交多项式函数的近似基本解研究

为了克服在用边界元法求解正交各向异性板的动力问题时寻求相应的基本解的困难,本文探讨了二种类型的近似基本解:级数形式及迭加修正形式的近似基本解。前者由正交函数系组成,本文中采用了Hermite正交多项式。后者由相近问题(广义各向同性板静力弯曲问题)的解析基本解迭加上用以修正的级数项组成。一些算例说明了近似基本解方法的可行性和有效性。迭加修正形式的解具有较高精度和较快收敛性,比单纯的级数形式解来得优越。