基于自适应交叉近似边界元法的快速声学灵敏度分析

基于自适应交叉近似边界元法构造一组快速声学灵敏度分析方法,其中,声学灵敏度分析分别采用直接微分法和伴随变量法;而自适应交叉近似算法被用以克服常规边界元法的高计算量和高存储量的固有缺点。自适应交叉近似算法在迭代求解之前对边界元系数矩阵进行压缩存储,可以在降低存储量的同时提高求解效率。在声学灵敏度分析中,通过直接使用求解未知边界状态值时保存的压缩系数矩阵,可以进一步提高求解效率。数值算例验证了所构造的方法的计算精度和求解效率,以及在大规模声场问题的最优化分析中的应用潜力。     

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

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

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

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

弹性力学问题的虚边界元-配点法

本文提出虚边界方法,建立了离散化虚边界元-配点法,给出了离散化求系数的积分解析式。本文方法完全避免了边界奇异积分及其复杂耗时的运算,成功地提高了普通边界元法(以下简称边界元法)中边界附近区域内包括边界上解的精度,保留了边界元法的优点并扬弃了其弱点。边界元间接法是本文方法中的一个特例。数值算例表明,程序可靠,节省机时,计算精度较高。     

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

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

动态断裂力学的无限相似边界元法

对弹性动力学的相似边界元法进行了进一步研究,推导了相应的计算公式,并在此基础上提出了动态断裂力学的无限相似边界元法.与传统的边界元法相比,相似边界元法由于只需在少数单元上进行数值积分,大大减少了计算量.对动态断裂力学问题,无限相似边界元法由于在裂纹尖端的边界上设置了逼近于裂纹尖端的无限个相似边界单元,可直接得到裂纹尖端具有奇异性的应力,而不需要设置奇异单元,从而突破了奇异单元对应力奇异性阶次的局限.另外,还讨论了无限相似边界元法得到的无限阶的线性代数方程组的求解方法.     

断裂力学的相似边界元法及其应用

首先对弹性力学的相似边界元法进行了研究,推导了相应的计算公式。与传统的边界元法相比,相似边界元法由于只需在少数单元上进行数值积分,当边界单元数目较多时大大减少了计算量。在此基础上,将相似边界元法应用于断裂力学,对路面断裂力学问题进行了计算,与有限元法的结果比较,说明了本文方法在减少计算量的情况下仍能较好地保证精度。     

动力分析中的边界元法(Ⅱ)

<正> 3 其他边界方法 第1,2两节较详细介绍了频域中和时域中的直接边界元法。这种边界元法称为直接的,是因为它讨论了具有重要物理意义的量,例如位移和拉力等,从而成为不仅在弹性动力学中而且在这种方法的所有其他力学应用中最广泛采用的边界元法。反之,以没有物理意义的某些虚拟量或源密度表述问题的间接边界元法是不太流行的,尽管事实上它具有比直接法更长的历史,而且与积分方程法方面的经典工作密切相关。     

无单元Galerkin法和边界元耦合法

无网格方法与有限元法或边界元法耦合是无网格方法处理边界条件的方法之一,在无网格方法中研究无网格方法与有限元法或边界元法耦合的研究显得非常重要.本文在无单元Galerkin法和边界元法的基础上,基于无单元Galerkin法子域和边界元法子域的界面上位移连续和面力平衡条件,提出了一种新的无单元Galerkin法和边界元法的直接耦合方法,对弹性力学问题详细推导了在整个求解域上的耦合公式.与以往的耦合法相比,这种方法简单直观,不需要增加新的耦合区域,也不需要建立新的逼近函数来保证界面位移的连续性.算例结果表明,该方法具有较好的计算精度.     

热弹性平面问题的规则化边界积分方程

对于热弹性平面问题,过去广泛集中在直接变量边界元法研究,本文研究间接变量规则化边界元法,建立了间接变量规则化边界积分方程。和直接边界元法相比,间接法具有降低密度函数的连续性要求、位移梯度方程中的热载荷体积分具有较弱奇异性等优点。数值实施中,用精确单元描述边界几何,不连续插值函数逼近边界量。算例表明,本文方法效率高,所得数值结果与精确解相当吻合。     

瞬态热传导问题的精细积分-双重互易边界元法

采用双重互易边界元法结合精细积分法求解二维含热源的瞬态热传导问题。针对边界积分方程中热源项和温度关于时间导数项引起的域积分,采用双重互易法处理,将域积分转换为边界积分。采用边界元法将边界积分方程离散后,得到关于时间的微分方程组,并利用精细积分法处理其中的指数型矩阵;对于微分方程组中由边界条件和热源项引起的非齐次项,采用解析的方法计算。为了比较精细积分-双重互易边界元法的计算效果,同时使用有限差分法计算温度对时间的导数项。通过数值算例验证了本文方法的有效性和精确性。计算结果表明:时间步长对于精细积分-双重互易边界元法的结果影响较小,而有限差分法对时间步长比较敏感且只在时间步长选取较小时有效;当选取较大时间步长时,精细积分-双重互易边界元法依然具有良好的计算精度。     

边界元法分析薄形层合结构

薄形层合结构由于几何形状的特殊性,其力学分析是数值计算的难点.边界元法分析层合结构具有较大的优势,但对于薄形层合结构,边界源点和对边上的积分单元距离很近,边界积分方程中存在几乎奇异积分,常规的数值积分方法已经失效.文章引入一种半解析化方法,计算薄形层合结构边界元法中的几乎奇异积分,使边界元法能成功分析三维薄形层合结构的层间界面应力和各层内点力学参量.     

模拟相变问题的精细积分边界元法

将精细积分边界元法和界面追踪法相结合求解相变问题。因为边界元法只需要将待求解空间域的边界离散,方便连续追踪移动界面位置和重构网格,所以边界元法适合应用于移动边界问题的模拟。首先,利用精细积分边界元法在固相区域和液相区域分别求解相应的瞬态热传导控制方程,从而求得温度场和边界热流密度。然后,根据固-液相变界面上的能量平衡方程,利用热流密度求得相变界面的移动速度,再采用界面追踪法预测移动相变界面的位置变化。最后,给出了几个数值算例,并通过与参考解的对比验证本文方法的准确性。     

一种求解瑞利波散射问题的修正边界元方法

边界元法的一大优势是用于求解半空间等无限域问题,然而对于弹性波的传播问题,传统边界元法在采用全平面或全空间格林函数时,在截断边界处仍会产生虚假的反射回波,直接影响到散射场的求解准确性。因此,本文在传统边界元法基础上提出一种修正边界元法,用于计算无限大半平面中的弹性波场问题。该方法以瑞利波形式的远端散射场代替原本因截断而舍去的部分,通过互易定理建立单位瑞利波和全平面格林函数的积分方程,求得修正系数,并代入修正边界元矩阵,计算出瑞利波的散射场。为验证本文所提方法,文中将多个算例的结果与解析解对比,并用该方法计算了不同缺陷的散射场。这些对比结果表明,本文所提修正边界元法可准确求解瑞利波散射场,为基于表面波的缺陷反演问题研究提供了有效的正演途径。     

固体力学边界元分析程序包BESMAP——基本结构功能

要使边界元法象有限元法那样得到广泛应用，必须开发与有限元法相媲美的大型边界元法程序包。国外虽出现少数几个边界元法程序包，但由于结构、功能等方面的原因，还不能成为边界元法广泛应用的工具。针对这种情况，本文借鉴现有的先进的有限元法程序包，研制了大型固体力学边界元法程序包ＢＥＳＭＡＰ，ＢＥＳＭＡＰ在总体设计、计算能力、使用方便和功能齐全等方面作了研究和改进。     

扩展边界元法分析切口和裂纹结构应力场的准确性

在线弹性理论中,切口/裂纹结构尖端区域存在奇异应力场,数值方法不易求解。本文建立的扩展边界元法(XBEM)对围绕尖端区域位移函数采用自尖端径向距离 的渐近级数展开式表达,其级数项的幅值系数作为基本未知量,而外部区域采用常规边界元法离散方程。两者方程联立求解可获得切口和裂纹结构完整的位移和应力场。扩展边界元法具有半解析法特征,适用于一般的切口和裂纹结构应力场分析,其解可精细描述从尖端区域到整体结构区域的应力场。作者研制了扩展边界元法程序,文中给出了两个算例,通过计算结果分析,表明扩展边界元法求解切口和裂纹结构应力场的准确性和有效性。     

快速多极子展开技术在高阶边界元方法中的实现

高阶边界元法以较常数元方法计算精度高存储低而在工程计算中得到了广泛的应用,但由于其平方存储和计算量的本质,无法应用于大型工程问题中。本文将快速多极子方法（FMM）应用于高阶边界元中从而使其计算量和存储量分别降为O（Nlog N）和O（N）。通过无限区域中水流绕射算例的数值计算,对FMM高阶边界元法与传统高阶边界元法的运算速度和内存消耗进行了分析对比,结果表明对于大型计算问题FMM高阶边界元算法更加有效。     

三维边界元法求解浮体水动力系数及与试验结果的比较

简要介绍了三维边界元法求解浮体水动力系数的方法；着重介绍了用三维边界元法求解程序的编制方程，给出了比较具体的程序框图，同时给出了理论计算曲线和试验结果的比较。     

书刊评介

本书系作者根据1981年以来讲授边界元法的讲义编写成,主要讲述直接边界元法在求解非弹性问题方面的应用,可以认为本书是边界元法发展现状的一个简明介绍。本书分三个主要部分,第一部分是弹性连续介质的基础知识概述,并介绍了无粘性和粘性材料的性质,还试图阐明非弹性介质的初始应变/初始应力表达式所涉及的符号,第二部分是弹性问题边界元法的引论,它既涉及边界元法的理论背景,也涉及它的数值实现,若干      

断裂力学的半解析相似边界元法

本文首先对弹性力学的相似边界元法进行了研究，推导了相应的计算公式。与传统的边界元法相比，相似边界元法由于只需在少数单元上进行数值积分，大大减少了计算量。在此基础上，对断裂力学问题，利用裂纹尖端位移场的解析表达式将裂纹尖端节点未知量转化为几个待定常数，提出了半解析相似边界元法，可大大减少最终形成的线性代数方程组的系数矩阵的阶数，进一步减小计算量。最后给出了算例，说明了本文方法的有效性。