Solution of multiple crack problem in a finite plate using coupled integral equations

This paper studies a numerical solution of multiple crack problem in a finite plate using coupled integral equations. After using the principle of superposition, the multiple crack problem in a finite plate can be converted into two problems: (a) the multiple crack problem in an infinite plate and (b) a usual boundary value problem for the finite plate. For the former problem, the Fredholm integral equation is used. For the latter problem, a BIE based on complex variable is suggested in which a Cauchy singular kernel exists. For the proposed BIE, after using the inverse matrix technique, the dependence of the traction at a domain point from the boundary tractions is formulated indirectly. This is a particular advantage of the present study. Several numerical examples are provided and the computed results for stress intensity factor and T-stress at crack tips are given.     

A boundary element solution to scattering of elastic SH wave by arbitrarily shaped inclusions embedded in an infinite anisotropic medium

Scattering problems for inhomogeneous bodies are investigated by the integral equation method. The boundary integral equation (BIE) for the scattered displacement field associated with finite inhomogeneities in an anisotropic medium are derived with the help of the generalized Green's identity. The discretization of BIE is based upon the constant element, linear element and quadratic element. Several numerical examples for calculating the scattering displacement, stress and scattering cross section from a cylinder, an interface crack, and two elliptic cylinders are given. Results show that the present method can be advantageously applied to a wide range of scattering problems of elastic waves.     

含有非圆形双孔的无限平板中应力的解析解研究

无限平板中含有任意形状单个孔的问题可以使用复变函数方法获得其应力解析解.对于无限平板中含有两个圆孔或两个椭圆孔的双连通域问题,也可以利用多种方法进行求解,比如双极坐标法、应力函数法、复变函数法以及施瓦茨交替法等.其中复变函数中的保角变换方法是获得应力解析解的一个重要方法.但目前尚未见到用此方法求解无限板中含有一个正方形孔和一个椭圆孔的问题.当板在无穷远处受有均布载荷和孔边作用垂直均布压力时,利用保角变换方法可以求解板中含有两个特定形状孔的问题.该方法将所讨论的区域映射成象平面里的一个圆环,其中最关键的一步是找出相应的映射函数.基于黎曼映射定理,提出了该映射函数一般形式,并利用最优化方法,找到了该问题的具体映射函数,然后通过孔边应力边界条件建立了求解两个解析函数的基本方程,获得了该问题的应力解析解.运用ANSYS有限单元法与结果进行了对比.研究了孔距、椭圆形孔大小和两孔布置方位对边界切向应力的影响,以及不同载荷下两孔中心线上应力分布规律.      

Numerical analysis of the stress concentration near holes originating in previously loaded viscoelastic bodies at finite strains

The state of stress and strain of previously loaded viscoelastic bodies with holes originating in them, successively or simultaneously, is analyzed under finite plane deformations. The problem statement and solution are based on the theory of repeatedly superimposed large deformations. The material mechanical properties are described using integral relations of the convolution type over time with a weakly singular kernel. The problem solving is based on the finite-element method. To calculate the integral of the convolution type, a recurrence formula is used that can be obtained by approximating the initial kernel with a linear combination of exponential functions (the truncated Prony’s series). The nonlinear effects and the effect of the interaction between holes on the stress concentration are analyzed. For the dynamic problems, the results for incompressible and weakly compressible materials are compared.     

用有限积分计算曲率复杂分布下的结构挠度

有限积分法是Brown和Trahair在求解微分方程时采用的数值解法, 其核心环节是已知函数z= z(x)的导函数z′ =z′(x)的某些值的情况下数值分析z的方法. 由曲率φ计算挠度, 实质意义上是由z″计算z的数学问题. 基于有限积分法给出的z与z″之间及z′与z″之间的数值关系, 通过矩阵运算推导得到了挠曲矩阵,通过引入转换式φ= -z″得到了曲率挠度关系式, 讨论了几种常见边界条件下的曲率挠度关系, 提出了曲率复杂分布情况下结构挠度计算的有限积分方法.     

周期界面裂纹反平面问题的动态应力强度因子

在研究动载荷作用下复合材料层板结构的安全与可靠性问题以及在抗震设计中关于地层裂缝的运动等问题中,都与界面裂纹有关。本文研究了分布于两个半空间之间的周期界面裂纹在反平面剪切波作用下裂纹尖端应力强度因子的动态特性。文中利用有限 Pourier变换,将在一个周期带内的边值问题转化成求解一个带周期性奇异核的积分方程,再借助于Chebyshev 多项式求得问题的级数解,最后分析了应力场在裂纹尖端的奇异性,得到了裂纹尖端动态应力强度因子的计算公式,并通过数值计算给出了应力强度因子随入射波频率变化的特性曲线。     

The scattering of harmonic elastic anti-plane shear waves by a finite crack in infinitely long strip using the non-local theory

In this paper, the scattering of harmonic anti-plane shear waves by a finite crack in infinitely long strip is studied using the non-local theory. The Fourier transform is applied and a mixed boundary value problem is formulated. Then a set of dual integral equations is solved using the Schmidt method instead of the first or the second integral equation method. A one-dimensional non-local kernel is used instead of a two-dimensional one for the anti-plane dynamic problem to obtain the stress occurring at the crack tips. Contraty to the classical elasticity solution, it is found that no stress singularity is present at the crack tip. The non-local dynamic elastic solutions yield a finite hoop stress at the crack tip, thus allowing for a fracture criterion based on the maximum dynamic stress hypothesis. The finite hoop stress at the crack tip depends on the crack length, the width of the strip and the lattice parameter. Supported by the Post Doctoral Science Foundation of Heilongjiang Province, the Natural Science Foundation of Heilongjiang Province and the National Foundation for Excellent Young Investigators.     

Degenerate scale for multiply connected Laplace problems

The degenerate scale in the boundary integral equation (BIE) or boundary element method (BEM) solution of multiply connected problem is studied in this paper. For the mathematical analysis, we use the null-field integral equation, degenerate kernels and Fourier series to examine the solvability of BIE for multiply connected problem in the discrete system. Two treatments, the method of adding a rigid body term and CHEEF concept (Combined Helmholtz Exterior integral Equation Formulation), are applied to remedy the non-unique solution due to the critical scale. The efficiency and accuracy of the two regularizations are also addressed. For simplicity without loss of generality, the eccentric case is considered to demonstrate the occurring mechanism of degenerate scale.     

A boundary element method for calculating the SH waves scattering from arbitrarily shaped holes in an anisotropic medium

The scattering problem of elastic wave by arbitrarily shaped cavities in an infinite anisotropic medium is investigated by the boundary integral equation (BIE) method. The formulations of BIE are derived with the help of generalized Green's formula. The discretization of BIE is based upon constant elements. After confirmation of the accuracy of the present method, some numerical examples are given for various cavities in a full space, in which an isotropic body with a circular cylinder hole is used for comparison and good agreement is observed. It has been proved that the method developed in this paper is effective.     

Scattering of flexural wave in a thin plate with multiple circular holes by using the multipole Trefftz method

The scattering of flexural wave by multiple circular holes in an infinite thin plate is analytically solved by using the multipole Trefftz method. The dynamic moment concentration factor (DMCF) along the edge of circular holes is determined. Based on the addition theorem, the solution of the field represented by multiple coordinate systems centered at each circle can be transformed into one coordinate system centered at one circle, where the boundary conditions are given. In this way, a coupled infinite system of simultaneous linear algebraic equations is derived as an analytical model for the scattering of flexural wave by multiple holes in an infinite plate subject to the incident flexural wave. The formulation is general and is easily applicable to dealing with the problem containing multiple circular holes. Although the number of hole is not limited in our proposed method, the numerical results of an infinite plate with three circular holes are presented in the truncated finite system. The effects of both incident wave number and the central distance among circular holes on the DMCF are investigated. Numerical results show that the DMCF of three holes is larger than that of one, when the space among holes is small and meanwhile the specified direction of incident wave is subjected to the plate.     

Quasi-Green’s function method for free vibration of clamped thin plates on Winkler foundation

The quasi-Green’s function method is used to solve the free vibration problem of clamped thin plates on the Winkler foundation. Quasi-Green’s function is established by the fundamental solution and the boundary equation of the problem. The function satisfies the homogeneous boundary condition of the problem. The mode-shape differential equation of the free vibration problem of clamped thin plates on the Winkler foundation is reduced to the Fredholm integral equation of the second kind by Green’s formula. The irregularity of the kernel of the integral equation is overcome by choosing a suitable form of the normalized boundary equation. The numerical results show the high accuracy of the proposed method.     

有限长轴承非稳态油膜力自由-移动边界问题的有限条解法

讨论了不可压缩流体润滑的动载径向滑动轴承油膜压力分布的自由移动边界问题的有限条解法．将自由边界问题转化为全域(矩形域)的具有不等式约束的微分方程的边值问题，进一步化为具有不等式约束的泛函极值问题。借助有限条法在矩形域上离散这个泛函，得到了一个特殊的二次泛函的规化问题。通过变量平移变换，使其化为标准的二次规划问题。然后借助于牛顿非光滑算法，迭代求解非线性的互补方程。给出了有限长轴承真实的油膜应力分布。对于所求解方程的系数矩阵的高度稀疏性。给出了紧缩存储算法。节省了存储空间和减少了计算量。算例表明该方法是有效的。     

各向异性体内含任意孔洞对反平面波散射的边界元方法

本文借助于广义格林公式导出了用位移表示的各向异性介质中SH波入射时的边界积分方程.根据本文作者在文献[8]给出的基本解,求解了各向异性介质中孔洞对SH波的散射问题.边界积分方程的离散基于常数元模式.文中给出了一个圆柱、一个椭圆柱和两个椭圆柱形式的孔洞周围的位移场和应力场的数值结果.最后,对入射波频率较高时的情形作了说明.     

A class of two-dimensional dual integral equations and its application

Because exact analytic solution is not available,we use double expansion and boundary collocation to construct an approximate solution for a class of two-dimensional dual integral equations in mathematical physics.The integral equations by this procedure are reduced to infinite algebraic equations.The accuracy of the solution lies in the boundary collocation technique.The application of which for some complicated initial- boundary value problems in solid mechanics indicates the method is powerful.     

On the dynamics of rotating and rolling structures

The dynamics of rotating flexible bodies is influenced by gyroscopic effects, which yields complex-valued eigenvectors and a split of eigenfrequencies in comparison with a resting body depending on rotational speed. These phenomena, which have been well investigated for simple structures like elastic rings, also affect more complicated systems such as rolling drums in printing machines. The aim of this contribution is to give a clear introduction into the physics of rotating bodies, starting with a rather simple analytical model, followed by some basic experimental investigations and ending with a finite element approach of more complicated three-dimensional structures. Special care is taken for the numerical solution of complex eigenvalue problems. The results obtained in each step are discussed in detail regarding their physical meaning.     

Stresses in heterogeneous plates

Conclusion The numerical solution described in the preceding section shows that the heterogeneity of a plate affects seriously the strain distribution but not so much the stress distribution (although this need not be true in general). This conclusion which cannot be immediately introduced from the known knowledge of elasticity, tends to support the assumption made by Müller and Macherauch.Electronic computers allow to find numerical solutions of complicated differential equations. But direct application of an electronic computer to the evaluation of a problem such as treated here is not possible, since the boundary conditions are not prescribed concretely.The solution obtained in this paper also applies to more complicated problems, such as the problems of heterogeneous plates containing holes or notches, because it represents the solution at points distant from the holes or notches. Using this solution, one can prescribe the boundary conditions for those problems concretely, hence numerical solutions can be obtained by the use of an electronic computer.A grant-in-aid for scientific research was given for this study by the Educational Department of Japan.     

Amended influence matrix method for removal of rigid motion in the interior BVP for plane elasticity

A conventional complex variable boundary integral equation(CVBIE) in plane elasticity is provided. After using the Somigliana identity between a particular fundamental stress field and a physical stress field, an additional integral equality is obtained. By adding both sides of this integral equality to both sides of the conventional CVBIE, the amended boundary integral equation(BIE) is obtained. The method based on the discretization of the amended BIE is called the amended influence matrix method.With this method, for the Neumann boundary value problem(BVP) of an interior region,a unique solution for the displacement can be obtained. Several numerical examples are provided to prove the efficiency of the suggested method.     

带源参数的二维热传导反问题的无网格方法

利用无网格有限点法求解带源参数的二维热传导反问题，推导了相应的离散方程. 与 其它基于网格的方法相比，有限点法采用移动最小二乘法构造形函数，只需要节点信息，不 需要划分网格，用配点法离散控制方程，可以直接施加边界条件，不需要在区域内部求积分. 用有限点法求解二维热传导反问题具有数值实现简单、计算量小、可以任意布置节点等优点. 最后通过算例验证了该方法的有效性.     

BIEM analysis of dynamically loaded anti-plane cracks in graded piezoelectric finite solids

Anti-plane cracks in finite functionally graded piezoelectric solids under time-harmonic loading are studied via a non-hypersingular traction based boundary integral equation method (BIEM). The formulation allows for a quadratic variation of the material properties in two directions. The boundary integral equation (BIE) system is treated by using the frequency dependent fundamental solution based on Radon transforms. Its numerical solution provides the displacements and tractions on the external boundary as well as the crack opening displacements from which the mechanical stress intensity factor (SIF) and the electrical displacement intensity factor (EDIF) are determined. Several examples for single and multiple straight and curved cracks demonstrate the applicability of the method and show the influence of the different system parameters.     

Certain effects occurring during boundary-layer control on a section of a surface

We examine the boundary layer on a plate upon which a circular cylinder is mounted. In investigations of this problem, the solution has successfully been found by numerical methods [1–3] only outside a certain domain whose location depends on the selection of the coordinate system, the spacings of the difference mesh [1], and, possibly, the computation method. Hence, a question arises as to whether it is impossible to solve the problem in a broader domain by modifying the algorithm; in other words, what is the maximum domain in which the solution of the boundary-layer equation can be found if the flow in the separation zone ahead of the obstacle is unknown? Application of the influence principle permitted the work set forth below to provide an answer to this question. The influence of cooling a section of surface on the boundary-layer characteristic was studied in [4] for parabolic external flow streamlines. The influence of cooling strips near the plane of symmetry on separation of a three-dimensional boundary layer is investigated numerically in the present paper. Experimental results that agree qualitatively with the flow singularities in the boundary layer exposed in the computation during boundary-layer control on a section of the surface are presented.