SOLUTION OF MULTIPLE CURVED RIGID LINE AND ANTIPLANE CIRCULAR INCLUSION PROBLEM BY WEAKLY SINGULAR INTEGRAL EQUATION

IntroductionMany investigators have investigated rigid lineinclusion and elastic inclusion problem.Early inthe 60 s,there were investigations aboutthe ex-perimental respects[1,2 ] and theoretical re-spects[3,4] of the inclusion problem . Wang etal[5 ]discussed the elastic stress field near the rigidline tips carefully.They put forward three frac-ture models of rigid lines corresponding to thethree fracture models of cracks and they pro-posed to use stress singularity coefficient( SSC)to desc…     

含有裂纹水泥混凝土路面的理论分析

以断裂力学理论为基础,将水泥混凝土路面视为Winkler地基上的弹性板,利用傅立叶积分变换,并引入位错密度函数,建立了奇异积分方程,推导出水泥混凝土路面表面存在裂纹时裂纹尖端应力强度因子的解析表达式,然后应用Lobatto-Chebyshev数值方法求解奇异积分方程,得到了应力强度因子的数值解.以实际路面为例,计算了裂...     

含裂纹水泥混凝土路面的解析算法

水泥混凝土路面在使用过程中出现裂纹是一种常见的损坏形式。以断裂力学理论为基础,对含裂纹水泥混凝土路面进行了理论分析。将水泥混凝土路面视为Winkler地基上的弹性板,利用傅立叶积分变换,并通过引入位错密度函数建立奇异积分方程,推导出水泥混凝土路面含有垂直裂纹时裂纹尖端的应力强度因子的解析表达式。然后,再应用Lobatto—Chebyshev法求解奇异积分方程,得到应力强度因子的数值解。为了说明问题,以实际路面为例,给出水泥混凝土路面内部存在裂纹时裂纹尖端应力强度因子的计算结果,并讨论了影响应力强度因子大小的因素。     

电磁材料弧形界面裂纹SH波散射的远场特性

分析了电磁材料弧形界面裂纹SH波散射的远场特性。根据波函数展开法将所研究的混合边值问题转化为对偶级数方程组,通过引入位错密度函数,进一步将问题归结为第一类Cauchy奇异积分方程。在此基础上,导出了远场散射模式（SFFP）和散射截面（SCS）的表达式,并详细考查了入射方向、裂纹张开角及载荷组合对远场散射特性的影响。研究成果对电磁构件的设计及电磁材料内裂纹的识别具有理论与应用价值。     

远处载荷作用下自由表面附近裂纹问题

采用连续分布位错模型讨论和分析了受远处载荷下半无限平面自由边界附近任意方位的裂纹问题,问题求解最后归结为求解一组Cauchy型奇异积分方程。在数值求解这组奇异积分方程中未知位错分布密度近似的基本密度函数与Chebyshev多项式之积,计算表明,本文方法的数值结果不但收敛快而且精度高。     

断裂力学中奇异积分方程的数值求解方法

对断裂力学中奇异积分方程的数值求解技术进行了综述。重点论述了第一类和第二类Cauchy型奇异积分方程以及超奇异积分方程的数值解法。这些方法的主要思想都是通过将奇异积分方程中的未知函数表示为多项式形式连续函数与特定形式权函数的乘积,然后借助Cauchy主值积分定义与超奇异积分的有限部积分定义,将奇异积分方程的求解归结为求解一组线性代数方程。本文拟结合一些具体的数值算例,对奇异积分方程中未知函数的不同表达方式、特点进行了评述,并比较了各种算法的优缺点。最后,指出了求解奇异积分方程的数值解法研究的未来发展方向。     

Fundamental formulation for frictional contact with graded materials

In the paper, we develop the fundamental solutions for a graded half-plane subjected to concentrated forces acting perpendicularly and parallel to the surface. In the solutions, Young's modulus is assumed to vary in the form of E(y) = E0e(ay) and Poisson's ratio is assumed to be constant. On the basis of the fundamental solutions, the singular integral equations are formulated for the unknown traction distributions with Green's function method. From the fundamental integral equations, a series of integral equations for special cases may be deduced corresponding to practical contact situations. The validity of the fundamental solutions and the integral equations is demonstrated with the degenerate solutions and two typical numerical examples.     

整数一次不定方程组的矩阵解法与程序设计

用欧几里德算法和整数环上的可逆线性变换，从理论上对整数一次不定方程组的解进行了深入研究，提出了用矩阵的初等变换求解整数一次不定方程组的矩阵解法，并利用MATLAB数学软件开发了相应的计算机程序．     

Steady state analysis of towed marine cables

Efficient numerical schemes were presented for the steady state solutions of towed marine cables. For most of towed systems, the steady state problem can be resolved into two-point boundary-value problem, or initial value problem in some special cases where the initial values are available directly. A new technique was proposed and attempted to solve the two-point boundary-value problem rather than the conventional shooting method due to its algorithm complexity and low efficiency. First, the boundary conditions are transformed into a set of nonlinear governing equations about the initial values, then bisection method is employed to solve these nonlinear equations with the aid of 4th order Runge-Kutta method. In common sense, non-uniform (sheared) current is assumed, which varies in magnitude and direction with depth. The schemes are validated through the DE Zoysa's example, then several numerical examples are also presented to illustrate the numerical schemes.     

Numerical Investigation of the Dynamics for Low Tension Marine Cables

This paper presents a numerical investigation into the dynamics of marine cables which are extensively used in offshore industry. In this numerical study, the Euler-Bernoulli beam model is adopted to develop the governing equations of the cable. Bending stiffness is considered to cope with the low tension problem in local area of towing cable, and thus a more accurate solution with the consideration of the axial elongation can be given.The derived strongly-coupled and nonlinear governing equations are solved by a second-order accurate, implicit,and large time step stable central finite difference method. The quadratically convergent Newton-Raphson iteration method is applied to solving the discrete nonlinear algebraic equations. Then a towed array sonar system(TASS)problem is studied. The numerical solutions agree reasonably well with the experimental data and the simulated results of the references. The specified program of the present paper shows great robustness with high efficiency.     

Fracture Mechanics for a Mode Ⅲ Crack in a Piezoelectromagnetic Material

Analytical solutions for a Griffith crack inside an infinite piezoelectromagnetic medium under combined mechanical-electrical-magnetic loadings are formulated using integral transform method. The singular stress, electric and magnetic fields in the piezoelectromagnetic material are obtained by the theory of linear piezoelectromagneticity. Fourier transforms are used to reduce the mixed boundary value problems of the crack, which is assumed to be permeable, to dual integral equations. The solution of the dual integral equations is then expressed analytically. Expressions for strains, stresses, electric fields, electric displacements, magnetic fields and magnetic inductions in the vicinity of the crack tip are derived. Field intensity factors and energy release rate for piezoelectromagnetic material are obtained. The stresses, electric displacements and magnetic inductions at the crack tip show the traditional square root singularities; and the electric field intensity factor (EFIF) and the magnetic field intensity factor (MFIF) are always zero.     

子矩阵约束下三类矩阵方程的对称正交对称迭代解法

讨论了子矩阵约束下三类矩阵方程的对称正交对称迭代解,利用广义共轭梯度法构造了迭代算法,并证明了算法的有限终止性.该算法能自动判定解的情况:当矩阵方程(组)相容时,得到矩阵方程(组)的解;当矩阵方程(组)不相容时,得到矩阵方程(组)的最小二乘解.     

Fracture Mechanics for a Mode Ⅲ Crack in a Piezoelectromagnetic Material

Analytical solutions for a Griffith crack inside an infinite piezoelectromagnetic medium under combined mechanical-electrical-magnetic loadings are formulated using integral transform method. The singular stress, elec tric and magnetic fields in the piezoelectromagnetic material are obtained by the theory of linear piezoelectromagneticity. Fourier transforms are used to reduce the mixed boundary value problems of the crack, which is as sumed to be permeable, to dual integral equations. The solution of the dual integral equations is then expressed analytically. Expressions for strains, stresses, electric fields, electric displacements, magnetic fields and magnetic inductions in the vicinity of the crack tip are derived. Field intensity factors and energy release rate for piezoelectromagnetic material are obtained. The stresses, electric displacements and magnetic inductions at the crack tip show the traditional square root singularities; and the electric field intensity factor （EFIF） and the magnetic field intensity factor （MFIF） are always zero.     

An Iterative Stabilized Scheme for Unsteady Incompressible Viscous Flow

An efficient iterative algorithm is presented for the numerical solution of viscous incompressible NavierStokes equations based on Taylor-Galerkin like split and pressure correction method in this paper. Taylor-Hood element is introduced to overcome the numerical difficulties arising from the fluid incompressibility. In order to confirm the properties of the algorithm, the numerical simulation on plane Poisseuille flow problem and liddriven cavity flow problem with different Reynolds numbers is presented. The numerical results indicate that the proposed iterative version can be effectively applied to the simulation of viscous incompressible flows. Moreover, the proposed iterative version has a better overall performance in maximum time step size allowed, under comparable convergence rate, stability and accuracy, than other tested versions in numerical solutions of the plane PoisseuiUe flow with different Reynolds numbers ranging from low to high viscosities.     

GPS-边值问题的自然边界元数值解法

对GPS-边值问题进行了讨论.在球近似下将GPS-边值问题转化为Neumann外问题,用自然边界元法进行自然边界归化,得到自然边界积分方程,再求自然边界积分方程的数值解.与其他边界元方法相比,自然边界元法的计算量大大减小,并且具有很好的逼近性质,能有效地处理奇异积分,是求解GPS-边值问题有效方法.     

一类双曲方程解的存在唯一性

对在等离子体和一维波动方程中出现的具体有初始分布数据的一类双曲微分方程正问题，将其化成第二类Volterra积分方程后，用压缩映像原理给出了解的存在唯一性。     

一类微分方程反问题解的存在唯一性

对在等离子体和一维波动方程中出现的具有初始分布数据的一类双曲微分方程反问题，用压缩映像原理给出了解的存在唯一性。     

物理大地测量边值问题的自然边界元解法

物理大地测量边值问题是物理大地测量学的核心问题和理论基础.在数学上表示为Laplace方程的以大地水准面为边界面的第三外部边值问题,即Stokes边值问题.文中在物理大地测量Stokes边值问题的球近似下,将Stokes边值问题转化为Neumann外问题,用自然边界元法,对Neumann外问题进行自然边界归化,得到自然边界积分方程,再求自然边界积分方程的数值解,并给出了一个算例.     

用面元法预报船舶的兴波阻力及波型

采用基于开尔文源的面元法求解航行于均匀流场中的船舶兴波问题．在船体湿表面上布置开尔文源，使其严格满足船体表面条件．数值计算中，首先采用分部积分对被积函数进行简化处理，然后采用高斯积分实现面元和线元上的积分，避免了被积函数为高频振荡函数所带来的数值计算的复杂性和不确定性．计算结果表明，其波型及兴波阻力系数与试验结果和其他数值方法相比均吻合较好．     

半无限平面压电双材料中导电界面边裂纹与螺位错的相互作用

研究了半无限大压电双材料中导电界面边裂纹与螺位错之间的相互作用问题。首先利用无限大压电双材料的位错解和扰动的概念,推导了具有导电边界半无限大压电双材料的位错解,然后使用保角变换方法,得到了含有限长导电界面边裂纹半无限大压电双材料的位错解。基于所得到的解,给出了应力强度因子和作用在位错上的像力的显函表达式,并通过数值算例分析了相互作用机理。