弹性板边界元中计算边界高阶导数的渐近收敛积分方程

本文针对板弯曲边界元方法中计算边界曲率等高阶导数项时边界积分方程中出现的高阶奇异积分项,通过对未知挠曲函数作渐近展开并加以适当摄动,获得了渐近收敛的边界积分方程。采用这一方法计算板边界上的曲率分布,获得了满意的数值结果。     

平头物体三维带空泡入水的数值模拟

本文用时间步进法和边界积分方程方法数值求解平头物体的垂直及斜向入水过程,这是一个在非线性自由面条件下,物体与流体有耦合作用,三维、非定常、理想不可压流体的运动问题,自由面用Lagrangian参数描述,物面用固结在物体上的Eulerian坐标描述,数值计算上提出了物面动力学条件(流-固耦合运动方程)和自由面动力学条件的二阶精度的时间差分隐格式,最后给出若干入水情况的详细计算结果。     

薄板弯曲问题的虚边界元-最小二乘法

本文提出求解任意形状的薄板弯曲问题的虚边界元-最小二乘法。本法首先利用薄板弯曲平衡方程的格林函数和离开实际边界上分布的未知的横向荷载和法向弯矩函数建立满足实际边界条件的积分方程;然后采用最小二乘法和沿虚边界分段离散化的待定的分布横向荷载和法向弯矩函数得到求上述积分方程离散化数值解的线性代数方程组。导出了一系列的数值积分的公式,并求解了许多例题,数值结果说明本法完全避免了奇异积分及其复杂的处理方法和耗时的运算,而且在边界及其附近区域解的精度比普通边界元(以后简称边界元)法大大地提高了。     

线性船舶兴波阻力边界元新方法

采用常数边界元对船舶与流体界面进行离散，求解船舶兴波势及船舶兴波阻力。这种方法可避免在船舶与流体自由面交线上安置节点，因而避免了这些节点建立补充方程。因为满足自由面条件的Ｈａｖｅｌｏｃｋ源函数的源点和场点不能同时在自由面上，使得自由面上的节点无法用Ｈａｖｅｌｏｃｋ源函数的建立方程。如对自由面交线上的节点建立补充方程，则要对线性自由面条件中包含的未知势函数的二阶导数用差分形式表达，引入较大误差。     

直接边界元法三维自由面兴波时域模拟

应用直接边界元法在时域中求解稳定航速运动的三维自由面兴波问题.基于格林定理,在所有边界面上划分网格,对边界积分方程进行数值离散,采用线性自由面边界条件,随时间步进更新自由面势.由于物体空间位置移动辐射条件不需要单独表述,迭代过程中自由面计算域保持不变.以割划水面NACA0024为例,计算模拟了自由面兴波稳定波形;提出了求解矩阵方程组奇异性的处理方法和解决割划问题的动网格技术.本文计算结果和有限体积法及有关试验结果对比表明,该方法是可靠的.     

反平面弹性中边缘内分叉裂纹的奇异积分方程解法

利用复变函数和奇异积分方程方法,求解反平面弹性中半平面边缘内分叉裂纹问题。提出了满足半平面边界自由的由分布位错密度表示的半平面中单裂纹的基本解,此基本解由主要部分和辅助部分组成。将半平面边缘内分叉裂纹问题看作是许多单裂纹问题的叠加,建立了以分布位错密度为未知函数的Cauchy型奇异积分方程组。然后,利用半开型积分法则求解奇异积分方程,得到了裂纹端处的应力强度因子。文中给出两个数值算例的计算结果。     

常航速船波势及船波阻力的边界元法

利用满足Laplace方程,线性化自由面条件及无穷远处条件的Havelock兴波源涵数,建立了关于常航速稳态船波势函数的边界积分方程.针对这个积分方程,建立了相应的数值计算方法,编制了一般三维问题的边界元法计算机程序,可用来计算全潜和半潜物体的稳态绕流场及船舶兴波阻力.     

横观各向同性压电材料中裂纹问题的边界元法

采用Somigiliana公式给出了三维横观各向同性压电材料中的非渗漏裂纹问题的一般解和超奇异积分方程，其中未知函数为裂纹面上的位移间断和电势间断．在此基础上，使用有限部积分和边界元结合的方法，建立了超奇异积分方程的数值求解方法，并给出了一些典型数值算例的应力强度因子和电位移强度因子的数值结果，结果令人满意．     

弹性力学平面问题的等价边界积分方程的边界轮廓法

基于边界积分方程中被积函数散度为零的特性,提出了弹性力学平面问题的等价边界积分方程的边界轮廓法,该方法无需进行数值积分,只需要计算单元两结点势函数值之差。实例计算说明,基于传统的边界积分方程的边界轮廓法所得到的面力结果是错误,而本文建立的边界轮廓法则可给出精确的结果。     

弹性力学的一种边界无单元法

首先对移动最小二乘副近法进行了研究，针对其容易形成病态方程的缺点，提出了以带权的正交函数作为基函数的方法－改进的移动最小二乘副近法，改进的移动最小二乘逼近法比原方法计算量小，精度高，且不会形成病态方程组，然后，将弹性力学的边界积分方程方法与改进的移动最小二乘逼近法结合，提出了弹性力学的一种边界无单元法，这种边界无单元法法是边界积分方程的无网格方法，与原有的边界积分方程的无网格方法相比，该方法直接采用节点变量的真实解为基本未知量，是边界积分方程无网格方法的直接解法，更容易引入界条件，且具有更高的精度，最后给出了弹性力学的边界无单元法的数值算例，并与原有的边界积分方程的无网格方法进行了较为详细的比较和讨论。     

三维非稳定渗流自由面边界积分项的精确数值计算

利用固定网格法分析三维非稳定渗流问题时,将要面对两项积分难题:以自由面及单元表面为边界的空间积分及以自由面为边界的曲面积分。针对常用的任意8结点6平面三维普通单元,提出采用坐标变换及等参变换技术求取空间积分项的精确数值解;至于曲面积分项,建议改用单元非饱和区部分表面作为积分边界,经过坐标变换及等参变换处理积分边界后,利用高斯数值积分可求出曲面积分项的精确数值解。通过一个普通单元及一项均质半无限边界堤坝的实例分析,表明此方法的精确性和稳定性良好。     

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

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

缓变主流中三维气泡的非线性振动

空化现象和水下噪声机制与液体中气泡的动力学行为密切相关．在无粘势流的假定下，采用多参数摄动分析，研究了缓变主流中三维气泡的非线性体积模态振动．推导了关于缓变泡形展开的各阶扰动方程，获得了一阶振动的演化方程和一些特殊情况下的解析解；并采用高阶有限元离散的边界积分方程方法，对平面固壁和自由面附近三维气泡的固有频率进行了数值计算     

双材料中平片裂纹问题的超奇异积分方程解法

利用三维断裂力学的超奇异积分方程方法,对双材料空间中重直于界面的平片裂纹Ⅰ型问题进行了研究。首先根据双材料空间的弹性力学基本解,使用边界积分方程方法,在有限部积分的意义下导出了以裂纹面位罗间断为未知函数的超奇异积分方程,并为其建立了数值法。在此基础上,讨论了用裂纹面位移问题计算应力强度因子的方法。最后用此计算了几个典型的Ⅰ型下片裂纹问题的应力强度因子,其数值结果令人满意。     

利用格林函数方法数值求解不可压粘性流非线性问题的研究

1.概述 本文对格林函数方法用以计算不可压粘性流非线性问题的能力进行了研究.该方法将定常运动的边值问题化为求解速度和边界应力的非线性积分方程,由于积分方程系完全精确推得,且在方程中可利用格林公式将速度、应力等物理量的微商转移为对基本解的微商,因而在数值计算中处理比较容易,且精度较高.     

弹性力学中一种新的边界轮廓法

利用基本解的特性,将面力积分方程化成仅含有Ｃａｕｃｈｙ主值积分的形式,基于这种边界积分方程,提出了一种新的边界轮廓法,对于三维问题,该方法只须计算沿边界单元界线的线积分,对二维问题,则只需计算边界单元两点的热函数之差,无须进行数值积分计算,实例计算说明该方法是有效的。     

受有集中力的样条边界元法

本文应用Ｂｅｔｔｉ定理，建立弹性体边界和体内受有有限个集中力时的样条边界积分方程，将集中力的影响表征为边界积分方程中的自由项，从客观实际出发，对具有两套奇性交会的积分方程给出一种方便有效的处理方法，使得集中力下的边界元法得以实施，在样插值基础上，即使稀疏剖分也能给出很高精度的位移场。应力场和未知集中反力。     

微重力环境下Marangoni对流的有限元数值模拟

本文以悬浮区为背景研究液桥中气液交界面上由表面张力所驱动的流体对流。我们采用有限元方法对轴对称、定常运动方程,能量方程,自由面切向、法向应力平衡条件迭代求解,首次考虑了边界形状的影响,得到了自洽的流区位形和流场、温度场、表面压强分布。一般说来,流区自由面呈弯月形。结果表明,只要表面张力数 C_a<1,静态平衡的流区位形就是动态情况的良好近似。本文还分析了C_a数及G_r 数对流区位形的影响,得到了不同 M_a 数及散热条件下的温度场和流函数分布。     

用Level Set方法求解具有自由面的流动问题

为采用Ｌｅｖｅｌ Ｓｅｔ方法来计算有自由的流动问题提出了一种方案,把自由水面视为水和空气的交界面,两种介质用统一的Ｎ－Ｓ方法求解,在自由面两侧分别采用各自的密度和粘性,并在自由面上给以适当的光滑;采用边界元法求解双调和方程来确定距离函数;Ｎ－Ｓ方程用投影法求解,文中给出了二维水池水面振荡和瞬时溃坝问题的算例,可以看出用ＬｅｖｅｌＳｅｔ方法求解有自由面流动问题是有效的。     

关于自由面压力脉动与摇板造波的等价性

考虑了一个逆问题，即在浮体运动和自由面压力分布之间是否存在某种等价性？作为这一方向努力的一部分，从二维半无限长水槽摇板造波的特例出发考察了这个问题．基于线性势流理论，导出了自由面压力分布应满足的带复数核的第一类Ｆｒｅｄｈｏｌｍ积分方程，并用正则化方法得到了这类不适定的积分方程的近似正则解．结果表明，自由面压力脉动与摇板造波之间的等价性可以通过解相应的逆问题来解决．