不定边界问题的边界元和有限元数学模型

对于不可压缩有势流动,有两类典型的不定边界或可动边界问题,即定常型不定边界如常水位闸门出流、过水坝或过水闸、水堰绕流及射流等,非定常型的不定边界如装液容器内的流体晃动问题.它们的共同特点是自由面位置与形状事先是的,需在计算过程中调整网格作自由面拟合.且由于自由面条件呈强烈非线性,给不确定数值计算带来困难.本文综述了两类不定边界问题的有限元和边界元模式,简述了笔者的一些计算经验.      

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

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

用面力边界积分方程求解断裂力学J积分

证明面力边界积分方程被积函数的散度等于零，应用Stokes公式，对平面线弹性问题，将面力边界积分的求解转化为边界点的位移势函数的点值计算。应用边界积分方程的求解结果，推导出J积分亦可表示为边界点的积分势函数的点值计算，无需进行数值积分，实例计算说明该方法的有效性。     

潜堤上波流传播的完全非线性数值模拟

利用时域高阶边界元方法建立模拟波流混合作用的完全非线性数值水槽模型, 其中自由水面满足完全非线性自由水面条件. 采用混合欧拉-拉格朗日方法追踪流体瞬时水面, 运用4阶Runge-Kutta方法更新下一时间步的波面和速度势. 为了减少计算域, 提高计算速度, 采用同时消除底面和侧面的镜像格林函数; 在每一时间步内, 对网格进行重新划分以避免由于网格运动变形而引起的数值不稳定问题. 对水流中淹没潜堤上的波浪变形在水槽中开展了物理模型试验, 并把试验结果和数值结果进行了对比, 吻合得很好. 进一步研究了水流及潜堤的存在对高阶谐波产生的影响.      

NUMERICAL SIMULATION OF 2D FIBER-REINFORCED COMPOSITES USING BOUNDARY ELEMENT METHOD

Introduction Withthedevelopmentofmodernindustry,compositesareincreasinglybeingappliedto agreatnumberofimportantstructures.Todeterminethemacroscopicaleffective characteristicsofcompositesisanessentialprobleminmanyengineeringapplications.The macroscopicalef…     

Boundary element analysis of contact problems using artificial boundary node approach

An improved version of the regular boundary element method, the artificial boundary node approach, is derived. A simple contact algorithm is designed and implemented into the direct boundary element, regular boundary element and artificial boundary node approaches. The exisiting and derived approaches are tested using some case studies. The results of the artificial boundary node approach are compared with those of the existing boundary element program, the regular element approach, ANSYS and analytical solution whenever possible. The results show the effectiveness of the artificial boundary node approach for a wider range of boundary offsets.     

快速求解叶栅流动的边界元方法

本文利用边界元方法，通过ｋｒｉｃｈｈｏｆｆ变换将描述叶栅流动的控制方程转换成线性方程。并将广义ｋ－Ｊ条件与边界积分方程联立求解，避免了非线性项和叶片出气角的迭代计算。完成了一种快速求解任意迥转面叶栅流场的计算程序，实用表明与其它数值方法及实验结果符合较好，具有快带、简明、实用的特点。     

Analytical Treatment of Boundary Integrals in Direct Boundary Element Analysis of Plan Potential and Elasticity Problems

ＩｎｔｒｏｄｕｃｔｉｏｎＡｃｏｍｐａｒｉｓｏｎ[1,2 ]ｂｅｔｗｅｅｎｔｈｅＤＢＥＭａｎｄｔｈｅｉｎｄｉｒｅｃｔｂｏｕｎｄａｒｙｅｌｅｍｅｎｔｍｅｔｈｏｄ (ＩＢＥＭ)ｈａｓｂｅｅｎｍａｄｅｂｙｅａｒｌｉｅｒｓｏｍｅａｕｔｈｏｒｓ.ＴｈｅｉｒｍａｉｎｖｉｅｗｐｏｉｎｔｓａｒｅｔｈａｔｔｈｅＤＢＥＭｈａｓａｎａｄｖａｎｔａｇｅｏｆｇｅｔｔｉｎｇｏｂｊｅｃｔｉｖｅｌｙｐｈｙｓｉｃａｌｑｕａｎｔｉｔｉｅｓｏｔｈｅｒｔｈａｎｉｍａｇｉｎａｒｙｏｎｅｓ ,ｂｕｔｉｔｒｅｑｕｉｒｅｓａｎｉｎｔｅｇｒａｔｅｄｎｕｍｅｒ…     

三维扰动波的非平行边界层稳定性研究

导出了三维扰动波的原始变量形式的抛物化稳定性方程（PSE）,研究了三维空间模态TS波的非平行边界层稳定性问题.采用了法向四阶紧致格式,以提高计算精度.通过给出不会导致奇性的坐标变换、修改外边界条件以及克服平行流初始值的瞬态影响和推进步长的限制,保证了计算的数值稳定.用补全元素带状矩阵法求解块三对角矩阵,大大提高了速度.计算结果清楚地显示了三维扰动波的演化过程和非平行性对边界层稳定性的影响,特别是,观察到非平行性对三维扰动波的影响,有时会使其稳定性出现逆转的现象.还研究了逆压梯度的作用.算例的结果与其他结果符合良好.     

边界积分方程中近奇异积分计算的一种变量替换法

准确估计近奇异边界积分是边界元分析中一项很重要的课题,其重要性仅次于对奇异积分的处理. 近年来已发展了许多方法,都取得了一定程度的成功,但这个问题至今仍未得到彻 底的解决. 基于一种新的变量变换的思想和观点,提交了一种通用的积分变换法, 它非常有效地改善了被积函数的震荡特性,从而消除了积分的近奇异性,在不增加计算量的情况 下, 极大地改进了近奇异积分计算的精度. 数值算例表明,其算法稳定,效率高, 并可达到很高的计算精度,即使区域内点非常地靠近边界,仍可取得很理想的结果.     

Non-homogeneous radiation and outflow boundary conditions simulating incoming acoustic and vorticity waves for exterior computational aeroacoustics problems

A set of non-homogeneous radiation and outflow boundary conditions which automatically generate prescribed incoming acoustic or vorticity waves and, at the same time, are almost transparent to outgoing sound waves produced internally in a finite computation domain is proposed. This type of boundary condition is needed for the numerical solution of many exterior aeroacoustics problems. In computational aeroacoustics, the computation scheme must be as non-dispersive and non-dissipative as possible. It must also support waves with wave speeds which are nearly the same as those of the original linearized Euler equations. To meet these requirements, a high-order/large-stencil scheme is often necessary. The proposed non-homogeneous radiation and outflow boundary conditions are designed primarily for use in conjunction with such high-order/large-stencil finite difference schemes. © 1998 John Wiley & Sons, Ltd.     

确定裂纹体等效弹性模量的边界元方法

采用边界元方法计算含有序分布裂纹的裂纹体在压缩载荷作用下的等效弹性模量，利用一种能适当考虑裂纹有间相作用的自洽理论，建立了相应的迭代格式，通过算例研究了裂纹方向，裂纹面间摩擦系数对裂纹体等效弹性模量的影响。     

水坝瞬态热传导分析中的拟初始条件边界元法

使用一种时域边界元方法对混凝土水坝进行瞬态热传导分析。在对时间积分进行离散计算时,采用一种拟初始条件法,即在时间步迭代计算的过程中,将之前计算结果对当前时间步的影响都视作当前时间步的初始条件。在所取时间步长较小的情况下,这种处理方法容易导致数值结果不稳定,即每一步的计算误差会累计放大,最终导致计算崩溃。本文提出一种虚拟时刻方法以缓解这类数值不稳定现象,在该方法中,时间步长首先放大至合适尺度,计算某个虚拟时刻(往往在真实计算时刻之后)的温度和流量分布,再通过插值方法换算出真实时刻的温度和流量分布。在虚拟时刻点上的温度和流量计算过程中,边界已知温度或流量由真实时刻的温度或流量进行外插得到。本文简单证明了该方法在温度和流量随时间呈线性变化情况下的正确性,最后给出了两个分析实例,验证了该方法的准确性和稳定性。     

Time‐accurate stabilized finite‐element model for weakly nonlinear and weakly dispersive water waves

Introduction of a time‐accurate stabilized finite‐element approximation for the numerical investigation of weakly nonlinear and weakly dispersive water waves is presented in this paper. To make the time approximation match the order of accuracy of the spatial representation of the linear triangular elements by the Galerkin finite‐element method, the fourth‐order time integration of implicit multistage Padé method is used for the development of the numerical scheme. The streamline‐upwind Petrov–Galerkin (SUPG) method with crosswind diffusion is employed to stabilize the scheme and suppress the spurious oscillations, usually common in the numerical computation of convection‐dominated flow problems. The performance of numerical stabilization and accuracy is addressed. Treatments of various boundary conditions, including the open boundary conditions, the perfect reflecting boundary conditions along boundaries with irregular geometry, are also described. Numerical results showing the comparisons with analytical solutions, experimental measurements, and other published numerical results are presented and discussed. Copyright © 2007 John Wiley & Sons, Ltd.     

Interaction of an explosion bubble with a fixed rigid structure

In the recent times, the boundary integral method has been utilised extensively in the study of bubble dynamics. This paper presents a modified form of the boundary integral method to model the motion of a bubble close to a fixed rigid structure. The resulting integral equations are solved using the boundary element method, and the system is integrated through time using a simple Euler scheme. Finally, the results of the model are presented to predict the motion of a bubble in a number of typical axisymmetric situations. Copyright © 1999 John Wiley & Sons, Ltd.     

The determination of an equivalent elastic modulus for bodies with cracks by boundary element method

The equivalent elastic modulus of cracked bodies with orderly distributed cracks was computed with the boundary element method. A practical self-consistent scheme has been proposed in consideration of the mutual interaction effects of the cracks. The influence of friction coefficients and orientation of cracks has been investigated. Some computational examples have been given, and the results show that the proposed method is adequate and the scheme is efficient.This project is supported by the National Natural Science Foundation of China.     

Equivalence between the Local Boundary Integral Equation and the Mean Value Theorem in the Theory of Elasticity

The boundary element method based on a boundary integral equation has been very successful in computational mechanics. Atluri et al. [4] recently developed a new meshless method using the local boundary integral equations. It eliminates the tedious step of mesh generation and thus greatly simplifies the numerical computation process. This paper shows the equivalence between the local boundary integral equation and the mean value theorem in the theory of elasticity. In addition, it gives new proofs for the mean value theorem of elasticity and its converse based on the concept of a companion solution. This revised version was published online in June 2006 with corrections to the Cover Date.     

基于多重多级子结构声子能带与传输特性分析

基于多重多级子结构方法提出一种快速的声子晶体能带与传输特性的计算策略. 主要思想是将声子晶体划分成多层级子结构有限元模型,在能带计算中采用静凝聚和子结构周游树技术将子结构的内部刚度阵凝聚到Bloch 边界上. 由于内部刚度阵并不随着简约波矢变化,所以这种计算策略可以大大降低求解规模并提高计算效率,并不对整体有限元模型引入近似. 在传输特性计算中同样采用该策略,由于声子晶体单胞具有周期性,所以各个单胞的系数矩阵是相同的,从而减少计算量,并且可以灵活地选择是否回代求解单胞内部自由度. 数值算例以三维局域共振型声子晶体和二维Bragg 散射型声子晶体为例,计算结果验证了这种求解策略的正确性和高效性,并适用于复杂声子晶体分析.      

BAND STRUCTURE AND TRANSMISSION CHARACTERISTICS CALCULATION OF PHONONIC CRYSTALS BASED ON MULTI-LEVEL SUBSTRUCTURE TECHNIQUE

基于多重多级子结构方法提出一种快速的声子晶体能带与传输特性的计算策略. 主要思想是将声子晶体划分成多层级子结构有限元模型,在能带计算中采用静凝聚和子结构周游树技术将子结构的内部刚度阵凝聚到Bloch 边界上. 由于内部刚度阵并不随着简约波矢变化,所以这种计算策略可以大大降低求解规模并提高计算效率,并不对整体有限元模型引入近似. 在传输特性计算中同样采用该策略,由于声子晶体单胞具有周期性,所以各个单胞的系数矩阵是相同的,从而减少计算量,并且可以灵活地选择是否回代求解单胞内部自由度. 数值算例以三维局域共振型声子晶体和二维Bragg 散射型声子晶体为例,计算结果验证了这种求解策略的正确性和高效性,并适用于复杂声子晶体分析.     

非连续边界元——有限元耦合方法分析

对边界元－有限元耦合方法进行了分析，采用非连续元离散边界积分方程，解决了耦合分析中自由度约束问题，给出了非连续边界元同有限元耦合的具体实施步骤，通过对二维弹性力学和流＝固耦合问题分析，表明了该文方法的有效性。