压电材料三维问题的虚边界元——最小二乘配点法

从压电材料三维问题的基本方程出发,利用已有的压电材料三维问题的基本解以及弹性力学虚边界元方法的基本思想和线性叠加原理,提出了压电材料三维问题的虚边界元——最小二乘配点解法。虚边界元解法继承了传统边界元方法的优点,并且有效避免了传统边界元方法中可能遇到的边界积分奇异性问题。最后,文章给出了压电材料三维问题的几个数值算例,并且与解析解做了比较,结果表明本文的方法具有较高的精度,是解决该问题一个十分有效的数值求解方法。     

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

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

虚边界元最小二乘配点法

本方法是在虚边界上进行数值积分，在实边界上有限个点处满足给定问题边界条件的一种数值算法。在虚边界上积分是为能寻求到使原问题得到正确解的较好的分布虚体力；在实边界上配点，是依据加权残数法中超额配点，即最小二乘配点法的思想。文中给出了本文方法的基本思想，并给出了壳体的算例。由数值结果表明本文方法的计算精度是令人满意的。     

求解位势问题的虚边界元法

本文提出了求解位势问题的虚边界元法，建立了位势问题的虚边界元的离散方程式，推导了离散化求系数的积分解析式。该方法与传统边界元法相比具有不存在奇异积分和边界附近精度较高等优点，可用来计算真空静电场、稳定温度场、流体绕流、介质中的渗流等各类位势问题。大量算例均获得了满意的结果。     

求解位势问题的虚边界元法

本文提出了求解位势问题的虚边界元法，建立了位势问题的虚边界元的离散方程式，推导了离散化求系数的积分解析式。该方法与传统边界元法相比具有不存在奇异积分和边界附近精度较高等优点，可用来计算真空静电场，稳定温度场，流体绕流，介质中的渗流等各类位势问题，大量算例均获得了满意的结果。     

SPH方法的运动边界条件研究

提出一种设置运动边界条件的方法,研究边界附近流体粒子积分截断和非物理穿透边界的问题。边界外的虚粒子在每个时间步由边界附近流体粒子对称生成,赋予相应的物理量,并在虚粒子中引入排斥力,利用拉格朗日形式的N-S方程自编SPH程序,参考一维激波管的精确解验证边界方法的适用性,研究运动边界条件在计算模型中应用。激波管的模拟结果与精确解基本一致,且在运动边界模型中也计算获得合理的结果。文中提出的运动边界条件,避免了边界附近流体粒子积分截断问题,阻止流体粒子在边界处发生非物理穿透的现象。     

压电材料平面问题的虚边界元-等额配点解法

利用压电材料平面问题的基本解和弹性力学虚边界元方法的基本思想，提出了压电材料平面问题的虚边界元-等额配点解法。该解法继承了传统边界元方法的优点，而避免了传统边界元方法遇到的边界积分奇异性问题。最后给出了压电材料平面问题的一些具体算例，并与解析解作了比较。结果表明本文的方法有很高的精度，是该问题一个十分有效的数值求解方法。     

二维弹性新型快速多极虚边界元的最小二乘法

在虚边界元最小二乘法的方程求解中采用新型的快速多极展开和广义极小残值法,提出了一种二维弹性新型快速多极虚边界元最小二乘法的求解思想。基于二维弹性问题原有的快速多极虚边界元最小二乘法的展开格式,通过引入对角化的概念,以更新展开传递格式;相对于原有快速多极算法,该方法可进一步提高计算效率且仍能保证具有较高的计算精度。数值算例说明了该方法的可行性、计算效率、计算精度均较高。     

快速多极虚边界元法对含圆孔薄板有效弹性模量的模拟分析

针对虚边界元法,引入快速多极展开和广义极小残值法(GMRES)的思想,以形成快速多极虚边界元法的求解思想,并将此方法用于含圆孔薄板有效弹性模量的模拟分析.由于本文方法采用了"源点"多极展开和"场点"局部展开的组合处理方案,从而使得原问题方程组求解的计算耗时量和储存量降至与所求问题的计算自由度数成线性比例.本文工作的研究目的在于:提高虚边界元法在普通台式机上的运算能力和拓宽虚边界元法对大规模复杂问题的求解(或数值模拟).文中给出了均布圆孔的正方形薄板和之字形分布圆孔薄板二个算例,以验证该方法的可行性,计算精度和计算效率.     

SPH方法的运动边界条件研究Research on the moving boundary condition in SPH method

提出一种设置运动边界条件的方法，研究边界附近流体粒子积分截断和非物理穿透边界的问题。边界外的虚粒子在每个时间步由边界附近流体粒子对称生成，赋予相应的物理量，并在虚粒子中引入排斥力，利用拉格朗日形式的N‐S方程自编SPH程序，参考一维激波管的精确解验证边界方法的适用性，研究运动边界条件在计算模型中应用。激波管的模拟结果与精确解基本一致，且在运动边界模型中也计算获得合理的结果。文中提出的运动边界条件，避免了边界附近流体粒子积分截断问题，阻止流体粒子在边界处发生非物理穿透的现象。     

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

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

伸缩虚拟边界元法解二维Helmholtz外问题

以位势理论为基础，提出了求解Helmholtz外问题的伸缩虚拟边界元法．给出了该方法在全波数域内获得唯一解的严格数学证明，其核心是通过伸缩虚拟边界使对偶内问题的特征频率(本征值)避开与波数重合，从而保证了解的唯一性，同以往前人提出的几种解法途径相比，该法简单得多；通过诸多边界曲线形状和不同边界量的声辐射算例，从计算精度、稳定性以及克服解的非唯一性等方面，对该方法进行了检验．计算结果表明：对远场或近场辐射声压，该方法都具有非常高的效率和精度．     

复合材料周期结构数学均匀化方法的一种新型单胞边界条件

数学均匀化方法是计算周期复合材料结构的有效方法之一,单胞边界条件施加的合理性直接决定了影响函数控制方程的计算效率和精度,进而影响均匀化弹性参数和摄动位移的计算精度.本文首先将单胞影响函数作为虚拟位移处理,给出了单胞在结构中真实的边界条件,结果表明,四边固支适合作为二维结构单胞边界条件;其次,针对二维结构提出了超单胞周期边界条件,有效提高了影响函数的计算精度,并使用与虚拟位移相对应的虚拟势能泛函验证超单胞周期边界条件的有效性;最后,利用数值分析验证多尺度渐进展开方法的计算精度,强调了二阶摄动的必要性.     

The boundary element method applied to elastic contact problems with friction

In the current paper the boundary integral equations (BIE) for elastic contact problems with friction are derived from the incremental virtual work principle. After introducing contact conditions of adhesion and slip into BIE all variants on boundary are made to discretize by quadratic isoparametric boundary element. In the current paper not only an auto-increment loading law is presented but also the iterative calculation laws for open, slip and adhesion condition are given. The results of numerical examples are satisfactory.     

Imposing virtual origins on the velocity components in direct numerical simulations

The relative wall-normal displacement of the origin perceived by different components of near-wall turbulence is known to produce a change in drag. This effect is produced for instance by drag-reducing surfaces of small texture-size like riblets and superhydrophobic surfaces. To facilitate the research on how these displacements alter near-wall turbulence, this paper studies different strategies to model such displacement effect through manipulated boundary conditions. Previous research has considered the effect of offsetting the virtual origins perceived by the tangential components of the velocity from the reference, boundary plane, where the wall-normal velocity was set to zero. These virtual origins are typically characterised by slip-length coefficients in Robin, slip-like boundary conditions. In this paper, we extend this idea and explore several techniques to define and implement virtual origins for all three velocity components on direct numerical simulations (DNSs) of channel flows, with special emphasis on the wall-normal velocity. The aim of this work is to provide a suitable foundation to extend the existing understanding on how these virtual origins affect the near-wall turbulence, and ultimately aid in the formulation of simplified models that capture the effect of complex surfaces on the overlying flow and on drag, without the need to resolve fully the turbulence and the surface texture. From the techniques tested, Robin boundary conditions for all three velocities are found to be the most satisfactory method to impose virtual origins, relating the velocity components to their respective wall-normal gradients linearly. Our results suggest that the effect of virtual origins on the flow, and hence the change in drag that they produce, can be reduced to an offset between the virtual origin perceived by the mean flow and that perceived by the overlying turbulence, and that turbulence remains otherwise smooth-wall-like, as proposed by Luchini (1996). The origin for turbulence, however, would not be set by the spanwise virtual origin alone, but by a combination of the spanwise and wall-normal origins. These observations suggest the need for an extension of Luchini’s virtual-origin theory to predict the change in drag, accounting for the wall-normal transpiration when its effect is not negligible.     

THE MESHLESS VIRTUAL BOUNDARY METHOD AND ITS APPLICATIONS TO 2D ELASTICITY PROBLEMS

A novel numerical method for eliminating the singular integral and boundary effect is processed. In the proposed method, the virtual boundaries corresponding to the numbers of the true boundary arguments are chosen to be as simple as possible. An indirect radial basis function network (IRBFN) constructed by functions resulting from the indeterminate integral is used to construct the approaching virtual source functions distributed along the virtual boundaries. By using the linear superposition method, the governing equations presented in the boundaries integral equations (BIE) can be established while the fundamental solutions to the problems are introduced. The singular value decomposition (SVD) method is used to solve the governing equations since an optimal solution in the least squares sense to the system equations is available. In addition, no elements are required, and the boundary conditions can be imposed easily because of the Kronecker delta function properties of the approaching functions. Three classical 2D elasticity problems have been examined to verify the performance of the method proposed. The results show that this method has faster convergence and higher accuracy than the conventional boundary type numerical methods.