Numerical solutions of Navier-Stokes equations with an integrated compartment method (ICM)

The most common numerical methods that are used by physical scientists to approximate partial differential equations employ finite differences and/or finite elements. In addition, compartment analyses have been adopted by ecological system analysts to simulate the evolution of processes governed by differential equations without spatial derivatives. An integrated compartment method (ICM) is proposed to combine the merits of these three numerical techniques. The basic procedures of the ICM are first to discretize the region of interest into compartments, then to apply three integral theorems of vectors to transform the volume integral to the surface integral, and finally to use interpolation to relate the interfacial values in terms of compartment values to close the system. These procedures are applied to the Navier-Stokes equations to yield the computational algorithm from which computer programs can be coded. The computer code is designed to solve one-, two-, or three-dimensional problems as desired. The program is applied to two simple cases: wake formation behind an obstacle in a channel and circulatory motion of a body of fluid in the square cavity. These preliminary applications have shown promising results.     

层状介质中多个非共面Griffith裂纹的弹性波散射问题研究

本文利用积分变换的方法,研究了层状介质中多个非共面Griffith裂纹的弹性波散射问题,导出了当入射波分别为P(SV)波及SH波时的散射对偶积分方程,并以双覆盖层半空间中的双裂纹为例,对反平面剪切波的散射场进行了求解和较为详尽的分析。最后通过数值计算,得到了许多有关的动力学特性曲线,从而揭示了层状介质中多裂纹弹性波散射问题中的某些内在规律。     

Rolling of a Rigid Cylinder over a Half-Plane with Initial Stresses

The Mehler–Fock transformation method is used to study the rolling of a rigid cylinder over a half-plane with initial stresses. The problem is reduced to a system of two dual integral equations, which are then reduced to a system of two integral Fredholm equations of the second kind. The system of integral equations is solved by the method of degenerate kernels. The dependences of the normal and tangential stresses on the elongation are plotted.     

层状介质中的双Griffith交界裂纹的SH波散射(反平面运动)——近场解

本文研究层状介质中的双Griffith交界裂纹的SH波散射。文中采用积分变换法将双裂纹的弹性波散射化为对偶积分方程,进而将其转化为第一类奇异积分方程组,借助于第一类Chebyshev多项式,给出了奇异积分方程组的解答,并得到了双裂纹的动态应力强度因子的计算公式     

Transient response of a rectangular piezoelectric medium with a center crack

The dynamic field intensity factors and energy release rates in a rectangular piezoelectric ceramic medium containing a center crack are obtained for boundary conditions of a permeable and an impermeable crack under electro-mechanical impact loading. An integral transform method is used to reduce the problem to two pairs of dual integral equations, which are then expressed as Fredholm integral equations of the second kind. Numerical values on the dynamic energy release rate are obtained to show the dependences upon the geometry and electric field.     

On the elastic strip with an internal crack

The paper presents a method to deal with an inclined crack in an elastic strip. No assumptions of symmetry are made. The method involves the solutions for a cracked plane and an uncracked strip and results in two coupled singular integral equations with finite interval of integration. A crack in a half-plane arises as a limiting case. For internal cracks the integral equations are of a standard type and do not present any numerical difficulties. Results are presented for loads according to the technical beam theory.     

关于加权残值法的充要性问题及其权函数的选择

本文论证了偏微分方程初-边值混合问题的加权残值方程的充要性条件和有关定理。关于充要性条件的讨论也为正确选择解析函数作为权函数,建立充要的加权残值方程提供了理论指导。重点讨论了边界积分方程的充要性。指出了两种满足充要条件的权函数和相应的边界积分方程。     

三维无限体中两平行平片裂纹相互干扰的超奇异积分方程解法

本文利用三维断裂力学的超奇异积分方程求解理论，对三维无限体中两平行平片裂纹在任意载荷作用下的相互干扰问题作了研究。首先导出了以裂纹面移间断（位借）为未知函数的超奇异积分方程组，然后为其建立了有限积分边界元法；在此基础上，讨论用了裂纹面位移间断计算应力强度因子的方法，最后用此计算了两平行平片裂纹的相对位置对裂前沿应力强度因子的影响，其数值结果令人满意。     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

Efficient methods free of irregular frequencies in wave and solid/porous structure interactions

The method of boundary integral equation is widely applied to compute and analyze wave–structure interactions in marine and offshore engineering, and the application is also seen in marine aquaculture to deal with waves and porous structure interactions. The application of the Fredholm integral equation of the second kind together with the free-surface Green function for a surface-piercing body suffers from irregular frequencies which may be confused with resonance peaks. A simple and efficient method to remove irregular frequencies in the wave–structure interactions is developed via enforcing null potential (and horizontal derivatives) on discrete points on the interior water-plane area and is referred to as overdetermined integral equations (and enhanced overdetermined integral equations), respectively. Structures with solid surface, porous surface and their blending are considered, and numerical results demonstrate the effectiveness of this method. In contrast to extended integral equations, the overdetermined integral equations are easy to implement and more time-efficient.     

Quasi-Green＇s function method for free vibration of simply-supported trapezoidal shallow spherical shell

The idea of quasi-Green＇s function method is clarified by considering a free vibration problem of the simply-supported trapezoidal shallow spherical shell. A quasi- Green＇s function is established by using the fundamental solution and boundary equation of the problem. This function satisfies the homogeneous boundary condition of the prob- lem. The mode shape differential equations of the free vibration problem of a simply- supported trapezoidal shallow spherical shell are reduced to two simultaneous Fredholm integral equations of the second kind by the Green formula. There are multiple choices for the normalized boundary equation. Based on a chosen normalized boundary equa- tion, a new normalized boundary equation can be established such that the irregularity of the kernel of integral equations is overcome. Finally, natural frequency is obtained by the condition that there exists a nontrivial solution to the numerically discrete algebraic equations derived from the integral equations. Numerical results show high accuracy of the quasi-Green＇s function method.     

圆形域多圆孔多裂纹反平面问题研究

本文运用复变函数及积分方程方法，求解了圆形域多圆孔多裂纹反平面问题，建立了两种类型的基本解。复叠加原理和所得的基本解并沿国圆孔和裂纹表面取待定的基本解密度函数，可得到一组以基本解密度函数为未知函数的Ｆｒｅｄｈｏｌｍ积分方程。通过该积分方程组的数值可以得到密度函数的离散值，进而得到了裂纹尖端的应力强度因子。     

The first integral method for solving some important nonlinear partial differential equations

Exact solutions of some important nonlinear partial differential equations are obtained by using the first integral method. The efficiency of the method is demonstrated by applying it for two selected equations.     

反平面圆形夹杂和多圆孔多裂纹相互作用问题

动用复变函数及积分方法方法求解了反平面圆形夹杂和多圆孔多裂纹相互作用问题。为解决该问题,建立了两种类型的基本解。利用叠加原理和所得的基本解没圆孔和裂纹表面取待定的基本解密度函数,可得一组Ｆｒｅｄｈｏｌｍ积分方程,通过积分方程组的数值求解,可以得到密度函数的离散值,进而得到应力强度因子。     

The Problem of an External Circular Crack Under Asymmetric Loadings

ＩｎｔｒｏｄｕｃｔｉｏｎＷｈｅｎｔｈｅｂｏｕｎｄａｒｙｉｎｔｅｇｒａｌｅｑｕａｔｉｏｎｍｅｔｈｏｄｉｓａｐｐｌｉｅｄｔｏｃｒａｃｋａｎａｌｙｓｉｓ,ｏｎｌｙｎｕｍｅｒｉｃａｌｓｏｌｕｔｉｏｎｓｃａｎｂｅｏｂｔａｉｎｅｄ ,ｓｕｃｈａｓ:ｔｈｅｔｙｐｉｃａｌｗｏｒｋｓｏｆＳｎｙｄｅｒａｎｄＣｒｕｓｅ[1],Ｃｒｏｕｃｈ[2 ],Ｂｌａｎｄｆｏｒｄｅｔａｌ.[3],Ｐｏｒｔｅｌａｅｔａｌ.[4 ],Ｂｕｉ[5 ],Ｗｅａｖｅｒ[6 ]ａｎｄＷＡＮＧｅｔａｌ.[7- 9].Ｉｔｉｓｖｅｒｙｄｉｆｆｉｃｕｌｔｔｏａｐｐｌｙｔｈｅｂｏｕｎｄａｒ…     

Weakly-singular,weak-form integral equations for cracks in three-dimensional anisotropic media

Singularity-reduced integral relations are developed for displacement discontinuities in three-dimensional, anisotropic linearly elastic media. An isolated displacement discontinuity is considered first, and a systematic procedure is followed to develop relations for the displacement and stress fields induced by the discontinuity. The singularity-reduced relation for the stress is particularly important since it is in a form which allows a weakly-singular, weak-form traction integral equation to be readily established. The integral relations obtained for a general displacement discontinuity are then specialized to an isolated crack and to dislocations; the relations for dislocations are introduced to emphasize their direct connection to corresponding results for cracks and to allow earlier independent findings for these two types of discontinuities to be put into proper context. Next, the singularity-reduced integral equations obtained for an isolated crack are extended to allow treatment of cracks in a finite domain, and a pair of weakly-singular, weak-form displacement and traction integral equations is established. These integral equations can be combined to obtain a final formulation which is in a symmetric form, and in this way they serve as the basis for a weakly-singular, symmetric Galerkin boundary element method suitable for analysis of cracks in anisotropic media.     

APPLICATION OF THE HYPERSINGULAR INTEGRAL EQUATIONS METHOD TO THE INTERACTION BETWEEN TWO PARALLEL PLANAR CRACKS IN 3-D ELASTICITY

By using the analytic theory of hypersingular integral equations in three-dimensional fracture mechanics, the interactions between two parallel planar cracksunder arbitrary loads are investigated. According to the concepts and method of finite-part integrals, a set of hypersingular integral equations is derived, in which theunknown functions are the displacement discontinuities of the crack surfaces. Then itsnumerical method is proposed by combining the finite-part integral method with theboundary element method. Based on the above results, the method for calculating thestress intensity factors with the displacement discontinuities of the crack surfaces ispresented. Finally, several typical examples are calculated and the numerical resultsare satisfactory.     

A boundary integral equation formulation for the Reissner's plates resting on two-parameter foundation

In this paper a boundary integral equation formulation for the Reissner's plates resting on a two-parameter foundation is established. With the aid of the Hormander Operator method, the equations of the corresponding fundamental solutions are converted into a sixth order partial differential equation with a scale function as an unknown. In order to reduce the equation further, two auxiliary functions are introduced. They satisfy a second and a fourth order equation respectively. The expressions of the auxiliary functions can be derived easily. The fundamental solutions of the Reissnei's plates on the two-parameter foundation arc expressed by a linear combination of the auxiliary functions and their derivatives. The boundary integral equations are formulated by the use of the weighted residual procedure. The fundamental solutions obtained are taken as the kernel functions of the boundary integral equations. A few examples are studied. The numerical results show high accuacy and efficiency of the present formulation.This work was supported by the National Natural Science Foundation of China.     

Penny-shaped interfacial crack between dissimilar magnetoelectroelastic layers

A penny-shaped interfacial crack between dissimilar magnetoelectroelastic layers subjected to magnetoelectromechanical loads is investigated,where the magnetoelectrically impermeable crack surface condition is adopted. By using Hankel transform technique,the mixed boundary value problem is firstly reduced to a system of singular integral equations,which are further reduced to a system of algebraic equations. The field intensity factors and energy release rate are finally derived. Numerical results elucidate the eects of crack configuration,electric and/or magnetic loads,and material parameters of the magnetoelectroelastic layers on crack propagation and growth. This work should be useful for the design of magnetoelectroelastic composite structures.     

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

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