Boundary element method to solve torsion problem of cracked cylinder

Separating the discontinuous solution by use of the single crack solution, together with the regular solution of harmonic function, the torsion problem of a cracked cylinder is reduced to solving a set of mixed-type integral equations and its numerical technique is then proposed by combining the numerical method of singular integral equation with the boundary element method. Several numerical examples are calculated which will be useful to engineering practice. The method proposed is characterized by its fine accuracy and convenience for using, which can be extended to the cases of multiple crack.The project supported by National Natural Science Foundation of China.     

Waves generated by an inclined-plate wave generator

This paper describes the characteristics of small-amplitude waves generated by a sinusoidally oscillating, inclined paddle-type wavemaker operating in a constant-depth channel. Two-dimensional, linearized potential flow is assumed. A semi-analytical method, the boundary collocation method, is used to establish the relationship between wave amplitude and paddle stroke. The numerical results are compared with the numerical results of the boundary integral equation method. It is found that the boundary collocation method is simpler and more flexible to implement and faster to compute. In addition, the numerical results are in reasonably good agreement with the laboratory experimental data. For the vertical wavemaker, which is a special case of the inclined wavemaker, an analytical series solution can be found. By using the boundary collocation method and the boundary integral equation method to solve the vertical wavemaker problem and comparing the results with the analytical series solution, it is found that the boundary collocation method yields a solution which is much more accurate than that from the boundary integral equation method. Finally, the relationships between wave amplitude and paddle stroke are established for different inclinations of the paddle-type wavemaker, based on the boundary collocation method.     

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

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

On the strong singular method for the problems of subsonic flow past obstacles

The purpose of this paper is to give a new investigation on the doublet elementary solution method of solving the problem of subsonic flow. A strong singular integral equation is produced. The definition and the calculating formula of the effective principle value of strong singular integral are given. From this integral equation, some numerical methods corresponding to the whole continuous distribution of elementary solution can be deduced. This method is applicable to calculating subsonic aerodynamic force.     

用样条边界元与奇异积分方程法联合求解多裂纹柱体的扭转

通过间解的分离，本文将径向多裂纹柱体的导曲函两个调和函数表示，使问题归为解一组混混合型积分方程。针对方程的特点，本文联合使用三次样条边界法与奇异积分方程的数值方法对所得方程建立了数值法，并对裂纹相交情形作了特殊处理。最后对工程中感兴趣的一些典型的多裂纹柱体的扭转作了例题计算，结果表明，本文方法具有收敛快，精度高的特点。     

Helmholtz边界积分方程中奇异积分间接求解方法

提出了间接求解传统Helmholtz边界积分方程CBIE的强奇异积分和自由项系数,以及Burton-Miller边界积分方程BMBIE中的超强奇异积分的特解法。对于声场的内域问题,给出了满足Helmholtz控制方程的特解,间接求出了CBIE中的强奇异积分和自由项系数。对于声场外域对应的BMBIE中的超强奇异积分,按Guiggiani方法计算其柯西主值积分需要进行泰勒级数展开的高阶近似,公式繁复,实施困难。本文给出了满足Helmholtz控制方程和Sommerfeld散射条件的特解,提出了间接求出超强奇异积分的方法。推导了轴对称结构外场问题的强奇异积分中的柯西主值积分表达式,并通过轴对称问题算例证明了本文方法的高效性。数值结果表明,对于内域问题,采用本文特解法的计算结果优于直接求解强奇异积分和自由项系数的结果,且本文的特解法可避免针对具体几何信息计算自由项系数,因而具有更好的适用性。对于外域问题,两者精度相当,但本文的特解法可避免对核函数进行高阶泰勒级数展开,更易于数值实施。     

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

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

饱和两相介质动力固结问题时域求解的精细时程积分方法

针对u-p形式的饱和两相介质波动方程,采用精细时程积分方法计算固相位移u,采用向后差分算法求解流体压力p,建立了饱和两相介质动力固结问题时域求解的精细时程积分方法。针对标准算例,对该方法的计算精度进行了校核。开展了该方法相关算法特性的研究,对采用不同数值积分方法计算非齐次波动方程特解项计算精度的差异进行了对比研究,并对采用不同积分点数目的高斯积分法计算特解项条件下计算精度的差异进行了对比研究。研究结果表明,(1)该方法具有良好的计算精度。(2)计算非齐次波动方程特解项的数值积分方法中,梯形积分法的计算精度最差,高斯积分法、辛普生积分法和科茨积分法都具有较好的计算精度。(3)增加高斯积分点数目对于提高计算精度的作用并不显著。     

三维动态裂纹问题的超奇异积分方程法

基于弹性材料的动态基本方程，结合广义Betti-Rayleigh互易等式与时域下的边界积分方程，推导得到时域下的超奇异积分方程组。引入Laplace域下的动态基本解，将经过主部分析的积分核函数分解为静态和动态部分，其中动态积分核不具有奇异性。在裂纹前沿附近单元，采用与理论分析一致的平方根位移模型。结合Lubich时间卷积实现拉氏变换，采用配置点法计算超奇异积分，获得问题的数值解。并针对椭圆裂纹算例编写Fortran程序，得到冲击荷载作用下张开型裂纹的动态应力强度因子变化规律，数值结果稳定且收敛速度快。     

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

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

含圆形嵌体弹性平面中径向裂纹问题的超奇异积分方程方法

根据含圆形嵌体平面问题在极坐标下的弹性力学基本解,使用Betti互换定理,在有限部积分意义下将问题归结为两个以裂纹岸位移间断为基本未知量、对于Ⅰ型和Ⅱ型问题相互独立的超奇异积分方程,对含圆形嵌体弹性平面中的径向裂纹问题进行了研究.根据有限部积分原理,建立了问题的数值算法.计算结果表明,嵌体半径、裂纹位置及材料剪切弹性模量等都对裂纹应力强度因子具有较为明显的影响.     

势问题的重构核粒子边界无单元法

将重构核粒子法和势问题的边界积分方程方法结合,提出了势问题的重构核粒子边界无单元法. 推导了势问题的重构核粒子边界无单元法的公式,研究其数值积分方案,建立了重构核粒子边界无单元法的离散化边界积分方程,并推导了重构核粒子边界无单元法的内点位势的积分公式. 重构核粒子法形成的形函数具有重构核函数的光滑性,且能再现多项式在插值点的精确值,所以该方法具有更高的精度. 最后给出了数值算例,验证了所提方法的有效性和正确性. }      

Solution of the laminar duct flow problem by a boundary integral equation method

Summary A boundary integral equation method is proposed for approximate numerical and exact analytical solutions to fully developed incompressible laminar flow in straight ducts of multiply or simply connected cross-section. It is based on a direct reduction of the problem to the solution of a singular integral equation for the vorticity field in the cross section of the duct. For the numerical solution of the singular integral equation, a simple discretization of it along the cross-section boundary is used. It leads to satisfactory rapid convergency and to accurate results. The concept of hydrodynamic moment of inertia is introduced in order to easily calculate the flow rate, the main velocity, and the fRe-factor. As an example, the exact analytical and, comparatively, the approximate numerical solution of the problem of a circular pipe with two circular rods are presented. In the literature, this is the first non-trivial exact analytical solution of the problem for triply connected cross section domains. The solution to the Saint-Venant torsion problem, as a special case of the laminar duct-flow problem, is herein entirely incorporated.     

Finite-part integral and boundary element method to solve three-dimensional crack problems in piezoelectric materials

Using the hypersingular integral equation method based on body force method, a planar crack in a three-dimensional transversely isotropic piezoelectric solid under mechanical and electrical loads is analyzed. This crack problem is reduced to solve a set of hypersingular integral equations. Compare with the crack problems in elastic isotropic materials, it is shown that for the impermeable crack, the intensity factors for piezoelectric materials can be obtained from those for elastic isotropic materials. Based on the exact analytical solution of the singular stresses and electrical displacements near the crack front, the numerical method of the hypersingular integral equation is proposed by the finite-part integral method and boundary element method, which the square root models of the displacement and electric potential discontinuities in elements near the crack front are applied. Finally, the numerical solutions of the stress and electric field intensity factors of some examples are given.     

垂直磁压电材料界面三维裂纹的超奇异积分法

应用有限部积分概念和广义位移基本解，垂直于磁压电双材料界面三维复合型裂纹问题被转 化为求解一组以裂纹表面广义位移间断为未知函数的超奇异积分方程问题. 进而，通过主部 分析法精确地求得裂纹尖端光滑点附近的奇性应力场解析表达式. 然后，通过将裂纹表面 位移间断未知函数表达为位移间断基本密度函数与多项式之积，使用有限部积分法对超奇异 积分方程组建立了数值方法. 最后，通过典型算例计算，讨论了广义应力强度因子的变化规 律.     

A solution of steady-state fluid flow in multiply fractured isotropic porous media

Herein a plane, steady-state fluid flow solution for fractured porous media is first presented. The solution is based on the theory of complex potentials, the theory of Cauchy integrals, and of singular integral equations. Subsequently, a numerical method is illustrated that may be used for the accurate estimation of the pore pressure and pore pressure gradient fields due to specified hydraulic pressure or pore pressure gradient acting on the lips of one or multiple non-intersecting curvilinear cracks in a homogeneous and isotropic porous medium. It is shown that the numerical integration algorithm of the singular integral equations is fast and converges rapidly. After the successful validation of the numerical scheme several cases of multiple curvilinear cracks are illustrated.     

Regular wave integral approach to numerical simulation of radiation and diffraction of surface waves

A regular wave integral method is developed in the discretisation of a linear hydrodynamic problem on radiation and diffraction of surface waves by a floating or submerged body. The velocity potential of the problem is expressed as a solution of a body boundary integral equation involving the pulsating free surface Green function or pulsating free surface sources distributed on the body surface. With the use of a discretisation on the regular wave integral rather than discretisations on the singular wave integral of the Green function as in earlier investigations, the singular wave integral is approximated as an expansion of regular (or nonirregular) wave potentials. Influence coefficients between pulsating free surface source points are computed by the approximate expansion together with Hess–Smith panel integral formulas. Thus the velocity potential solution is evaluated by a boundary element algorithm. The numerical results produced from the proposed method agree well with semi-analytic solution results.     

三维断裂力学的超奇异积分方程方法

本文利用有限部积分的概念和方法,严格地证明了三维弹性体中受任意载荷作用的平片裂纹问题的超奇异积分方程组,并对未知解的性态作了理论分析,得到了性态指数,在此基础上通过主部分析,精确地求得了裂纹前沿光滑点附近的奇性应力场,从而找到了以裂纹面位移间断(位错)表示的应力强度因子表达式,最后对所得的超奇异积分方程组建立了数值法,并用此计算了若干典型的平片裂纹问题,数值结果令人满意。     

表面堆载作用下群桩负摩擦研究

利用Biot固结理论和积分方程方法研究了表面有堆载的群桩负摩擦问题。根据基本解得出了群桩在圆形均布载荷作用下在时间域内的第二类Fredholm积分方程组。运用Laplace变换对上述积分方程组进行简化，求解上述积分方程组并进行相应的数值逆变换就可得出群桩在表面圆形均布载荷作用下的变形、轴力、孔压和桩侧摩阻力随时间的变化情况。     

混合式CRP面元法计算对比

目前,对于混合式CRP(contra-rotating propeller)的数值研究大部分是基于黏流方法,这主要是由于混合式CRP组成构件较多,运用势流方法求解存在一定的困难.但势流方法求解效率较高,在混合式CRP的设计方面,具有独特的优势.为建立混合式CRP的高效数值分析方法,从而为混合式CRP的设计工作奠定基础,研究首先将混合式CRP分成前桨和吊舱推进器两部分,以单个螺旋桨面元法和吊舱推进器整体面元法为基础,建立了混合式CRP的迭代面元法.随后对混合式CRP的几何结构以及桨叶面元奇点强度的变化特点进行分析,提出了将混合式CRP各构件进行整体计算的面元法,推导了整体面元法的计算公式.编译完成迭代面元法以及整体面元法的计算程序,并对混合式CRP的敞水性能进行了计算分析,两种面元法的计算结果与试验值的对比表明;迭代面元法与整体面元法计算得到的推力、扭矩系数相对误差均能控制在5%以内,但整体面元法的计算时间比迭代面元法节约1/3左右,同时整体面元法省去了迭代步骤,适宜于混合式CRP的设计.最后对研究中所建立的整体面元法的理论误差进行了分析.