一种确定结合材料界面端应力强度因子的数值方法

基于弹性力平面问题的基本方程，给出了结合材料界面端的应力奇异性特征方程以及位移场和奇异应力场。提出了一种确定结合材料界面端应力强度因子的数值外插方法。对界面端区域进行了有限元网格单元划分。经过具体实例检验进一步确定了求解应力强度因子的最佳方向，该数值外插法的计算结果精度符合工程应用的要求，为工程材料强度的评价提供了有效的计算途径。     

压电材料中切口/接头端部平面电弹性场奇异性有限元分析

为了对平面载荷作用下压电材料中切口或接头端部附近电弹性场奇异性问题进行分析，首先以应力平衡方程、Maxwell方程和和边界条件为基础，得到一种求解压电材料特征问题的弱式方程；其次，假定楔形切口或接头端部附近单元内位移和电势沿径向分布为指数形式，而周向方向分布则采用泡函数插值，将其代入弱式方程，建立一种只需对楔形切口或接头端部附近周边进行离散的一维简单有限元方法．压电材料的极化轴可以是任意方向．利用该有限元模型讨论了楔形切口角度、极化轴方向和边界条件对奇性场的影响．通过和其它特定情况下的现有解相比，证实了该文有限元数值方法的有效性，而且精度很高．     

A unified definition for stress intensity factors of interface corners and cracks

Based upon linear fracture mechanics, it is well known that the singular order of stresses near the crack tip in homogeneous materials is a constant value −1/2, which is nothing to do with the material properties. For the interface cracks between two dissimilar materials, the near tip stresses are oscillatory due to the order of singularity being −1/2 ± iε and −1/2. The oscillation index ε is a constant related to the elastic properties of both materials. While for the general interface corners, their singular orders depend on the corner angle as well as the elastic properties of the materials. Owing to the difference of the singular orders of homogeneous cracks, interface cracks and interface corners, their associated stress intensity factors are usually defined separately and even not compatibly. Since homogenous cracks and interface cracks are just special cases of interface corners, in order to build a direct connection among them a unified definition for their stress intensity factors is proposed in this paper. Based upon the analytical solutions obtained previously for the multibonded anisotropic wedges, the near tip solutions for the general interface corners have been divided into five different categories depending on whether the singular order is distinct or repeated, real or complex. To provide a stable and efficient computing approach for the general mixed-mode stress intensity factors, the path-independent H-integral based on reciprocal theorem of Betti and Rayleigh is established in this paper. The complementary solutions needed for calculation of H-integral are also provided in this paper. To illustrate our results, several different kinds of examples are shown such as cracks in homogenous isotropic or anisotropic materials, central or edge notches in isotropic materials, interface cracks and interface corners between two dissimilar materials.     

Specific features of solving the problem of compression of an orthotropic plastic material between rotating plates

An instantaneous flow with compression of a wedge-shaped layer of a rigid-plastic orthotropic material between rotating plates is considered under the assumption that the principal axes of anisotropy are rays emanating from the wedge angle and lines orthogonal to them and that the maximum friction law is valid on the plate surfaces. The solution is reduced to quadratures, and its asymptotic analysis is performed. It is found that the solution is singular near the friction surface in the general case, and conditions at which the singularity disappears are given. It is demonstrated that a rigid area can arise near the friction surface. The behavior of the resultant solution near the friction surfaces is compared with the behavior of known solutions for other models of rigid-plastic materials.     

On the asymptotic crack-tip stress fields in nonlocal orthotropic elasticity

The crack-tip stress fields in orthotropic bodies are derived within the framework of Eringen’s nonlocal elasticity via the Green’s function method. The modified Bessel function of second kind and order zero is considered as the nonlocal kernel. We demonstrate that if the localisation residuals are neglected, as originally proposed by Eringen, the asymptotic stress tensor and its normal derivative are continuous across the crack. We prove that the stresses attained at the crack tip are finite in nonlocal orthotropic continua for all the three fracture modes (I, II and III). The relative magnitudes of the stress components depend on the material orthotropy. Moreover, non-zero self-balanced tractions exist on the crack edges for both isotropic and orthotropic continua. The special case of a mode I Griffith crack in a nonlocal and orthotropic material is studied, with the inclusion of the T-stress term.     

纤维与基体的轴对称界面端附近的奇异应力场

基于各向性弹性力学空间轴对称问题的基本方程,研究了纤维与基体的轴对称界面端的应力奇异性,并给出了界最佳 近的奇异应力场。研究结果表明,该轴对称界面端的应力奇异性与平面应变状态下相应模型的应力奇异性完全相同,材料对界面端附近奇异应力场的影响可用丰个双材料组合参数描述。     

Singular stress fields near contact boundaries in a composite bolted joint

Singular stresses arising in the neighborhood of contact surfaces introduced in laminated orthotropic plates by mechanical joining with clamp-up were investigated by using local asymptotic solutions and full-field numerical analysis. Three-dimensional B-spline approximation of the displacements and a penalty function-based contact solution was used in the numerical analysis. Recent work has shown that fracture in bolted composite joints may initiate near the outer edge of the bolt head or washer away from the hole edge, particularly if the joint is preloaded. Material and geometric discontinuities exist in these regions, resulting in singular stress behavior. Asymptotic stress analysis was performed to obtain the power of singularity in these regions as a function of the bolt-head (washer) stiffness. Frictionless contact conditions were assumed. It was found that the characteristics of the stress singularity for such practically important combinations as titanium bolt-head and carbon fiber composite plate are similar to a crack in terms of the power of singularity and uniqueness of the singular term. Coefficients of the singular terms of the asymptotic expansion were determined by comparison with the numerical solution in the close vicinity of the singular contour. Good agreement between the asymptotic and numerical solution in the transition regions was observed.     

Crack edge displacement and elastic constant determination for an orthotropic material

A method for determining the elastic constants of an isotropic material, based on crack edge displacement data, is extended to an orthotropic material. Complex potentials are used to obtain the stresses and displacements for plane strain. Mode I crack problems in three mutually orthogonal planes are considered and solved. In particular, the expressions of crack edge displacements are obtained in an explicit form. An iterative statistical identification method, based on a Bayesan approach, is used to identify the elastic constants of an orthotropic medium from the Mode I crack displacements measured from the mid-point of the crack. Some graphics are displayed to illustrate the convergence of the pertinent parameters and the approach of the analytical displacements to their experimental values.     

A new boundary element approach of modeling singular stress fields of plane V-notch problems

In this paper, a new boundary element (BE) approach is proposed to determine the singular stress field in plane V-notch structures. The method is based on an asymptotic expansion of the stresses in a small region around a notch tip and application of the conventional BE in the remaining region of the structure. The evaluation of stress singularities at a notch tip is transformed into an eigenvalue problem of ordinary differential equations that is solved by the interpolating matrix method in order to obtain singularity orders (degrees) and associated eigen-functions of the V-notch. The combination of the eigen-analysis for the small region and the conventional BE analysis for the remaining part of the structure results in both the singular stress field near the notch tip and the notch stress intensity factors (SIFs).Examples are given for V-notch plates made of isotropic materials. Comparisons and parametric studies on stresses and notch SIFs are carried out for various V-notch plates. The studies show that the new approach is accurate and effective in simulating singular stress fields in V-notch/crack structures.     

Higher order fields for damaged nonlinear antiplane shear notch, crack and inclusion problems

Damaged nonlinear antiplane shear problems with a variety of singularities are studied analytically. A deformation plasticity theory coupled with damage is employed in analysis. The effect of microscopic damage is considered in terms of continuum damage mechanics approach. An exact solution for the general damaged nonlinear singular antiplane shear problem is derived in the stress plane by means of a hodograph transformation, then corresponding higher order asymptotic solutions are obtained by reversing the stress plane solution to the physical plane. As example, traction free sharp notch and crack, rigid sharp wedge and flat inclusion, and mixed boundary sharp notch problems are investigated, respectively. Consequently, higher order fields are obtained, in which analytical expressions of the dominant and second order singularity exponents and angular distribution functions of the near tip fields are derived. Effects of the damage and hardening exponents of materials and the geometric angle of notch/wedge on the near tip quantities are discussed in detail. It is found that damage leads to a weaker dominant singularity of stress, but to little stronger singularities of the dominant and second order terms of strain compared to that for undamaged material. It is also seen that damage has important effect on the angular distribution functions of the near tip stress and strain fields. As special cases, higher order analytical solutions of the crack and rigid flat inclusion tip fields are obtained, respectively, by reducing the notch/wedge tip solutions. Effects of damage and hardening exponents on the dominant and second order terms in the solutions of the crack and inclusion tip fields are discussed.     

正交异性材料平面裂纹尖端应力场

根据正交各向异性材料力学性能确定出了用应力函数表示的弹性力学基本方程,利用坐标变换和复变函数方法求解了正交异性材料平面裂纹体的应力边值问题。借鉴一般断裂力学解法构造了I型和II型裂纹问题的应力函数,推导出了正交各向异性板裂纹尖端区的奇异应力场。通过数值计算说明了裂纹尖端应力表达式的正确性,验证了裂尖前沿应力变化规律,即σx与材料特征参数h2成正比,而σy和τxy不随材料特性变化。     

The collinear crack problem for an orthotropic functionally graded coating–substrate structure

A theoretical treatment of antiplane crack problem of two collinear cracks on the two sides of and perpendicular to the interface between a functionally graded orthotropic strip bonded to an orthotropic homogeneous substrate is put forward. Various internal cracks and crack terminating at the interface and crack crossing the interface configurations are investigated, respectively. The problem is formulated in terms of a singular integral equation with the crack face displacement as the unknown variable. The asymptotic stress field near the tip of a crack crossing the interface is examined, and it is shown that, unlike the corresponding stress field in piecewise homogeneous materials, in this case, the “kink” in material property at the interface does not introduce any singularity. Numerical calculations are carried out, and the influences of the orthotropy and nonhomogeneous parameters and crack interactions on the mode III stress intensity factors are investigated.     

基于非协调元特征法的奇异性问题分析

提出了一个基于位移的、分析平面尖劈尖端奇性应力场和位移场问题的非协调FE特征分析法．该方法与过去原有求解裂纹尖端近似场的有限元特征分析方法导出公式的出发点不同，并且采用的单元形式为非协调元，尖劈尖端邻域内的位移场假定没有采用奇异变换技术，运用该方法处理了若干尖劈和接头的算例，所有的计算结果表明，该方法较原有方法使用的单元少而且精度高，具有应用广泛性。     

Cracking of an orthotropic substrate reinforced by an orthotropic plate

A solution method is derived to determine the stress intensity factors for both an internal crack and an edge crack in an orthotropic substrate that is reinforced on its boundary by a finite-length orthotropic plate. The method utilizes the Green’s functions for a pair of dislocations and a concentrated force on the boundary while invoking the concept of superposition. Enforcing the traction-free boundary condition along the crack surfaces and the continuity of displacement gradients along the plate/substrate interface results in a coupled system of singular integral equations. An asymptotic analysis of the kernels in these equations for the region of the junction point between the plate corner and the substrate boundary reveals the strength of the singularity in the case of an edge crack. The numerical solution of the integral equations provides results for the stress intensity factors for both an internal crack and an edge crack perpendicular to the substrate boundary and aligned with one of the corners of the plate. The present results have been validated against previously published stress intensity factors for an internal crack and an edge crack in an isotropic substrate.     

QFP结构中引脚焊缝界面端的奇异性热应力问题研究

由于引脚、印制电路板和焊接剂的热-机材料属性不同,在受到热载荷或机械载荷时,引脚焊接界面端会产生奇异性应力,有可能产生界面开裂.为了基于界面端奇异场来评价QFP结构引脚界面端力学行为,本文拟采用数值方法求解引脚焊缝任意角度尖劈界面端的应力强度系数.具体步骤为:首先,基于高次内插有限元特征分析法确定两相任意角度尖劈界面端的奇异性指数和应力角分布函数,并引入常数热应力项,获得热-机耦合奇异性应力场表达式;采用有限元分析技术和最小二乘拟合法来获得应力强度系数的数值解.文中考察了热-机材料属性对热载荷下焊接剂/印制电路板界面端应力强度系数的影响,并给出改善界面端热应力状态的建议.     

Isotropic/orthotropic two-layered strips containing cracks terminating at and going through an interface

Crack problems for isotropic/orthotropic two-layered strips have been investigated. A system of two singular integral equations can be derived by using Fourier integral transformation and boundary conditions of crack problems. After stress singularities at crack tips or other special points are determined for internal and edge cracks, and for cracks terminating at and going through the interface, the system of singular integral equations is solved numerically by Gauss-Jacobi or Gauss-Chebyshev integration formulas for stress intensity factors at the tips and other singular points of cracks. Finally, possible crack growth behavior for cracks approaching and going through the interface is discussed.     

Singular plastic fields in wedge indentation of pressure sensitive solids

The paper examines singular plastic fields induced near the tip of a wedge indentating a pressure sensitive solid. Plane strain conditions are assumed and material response is modelled by the small strain Drucker–Prager rigid/plastic constitutive law. A standard separation of variables solution is numerically generated for pure power-law hardening. Three possible measures of wall roughness are studied with an attempt to expose the coupling between wall friction and material pressure sensitivity. Sample calculations illustrate that stress singularity decreases with increasing friction, wedge angle and hardening exponent, but increases with pressure sensitivity. At large values of the hardening exponent, when the material is nearly perfectly plastic, effective stress contours approach the slip line limit. The concept of indentation index is introduced as a possible estimate for average indentation pressure.     

Influence of a thin isotropic coating on the edge effect zone in isotropic and transversely isotropic materials

The effect of a thin isotropic coating on the edge effect zone in a representative element of a coated material is examined. Isotropic and transversely isotropic materials are considered. The transversely isotropic material has the elastic properties of unidirectional glass-fiber-reinforced plastic. The decay of the edge effect in the directions perpendicular to the coating plane and to the plane of isotropy is studied. A boundary-value problem of elasticity for piecewise-homogeneouse orthotropic bodies and a quantitative edge effect decay criterion for normal stresses are used as a design model. The problem is solved using the finite-difference method and base schemes. The results of evaluation of the edge effect zone in homogeneous and inhomogeneous materials are presented. It is shown that the presence of a thin isotropic coating blocks the edge effect, that is, decreases the edge effect zone in both isotropic and transversely isotropic materials __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 12, pp. 61–67, December 2007.     

垂直于界面V型缺口尖端的应力奇异性及其应用

基于双材料垂直于界面V型缺口理论,给出了单一材料和双材料裂纹问题、V型缺口问题应力强度因子的统一定义,得到了应力外推法计算双材料K_I的公式,数值算例验证了本文方法的有效性.以双材料单向拉伸和三点弯曲模型为对象,深入研究了双材料中弹性模量、泊松比、缺口深度、缺口张角对缺口尖端奇异应力场的影响,获得了一定范围内各种参数变化对缺口尖端奇异应力场的影响规律,为异体材料形成的V型缺口在应力断料中的应用提供了必要的参考依据.     

High-order asymptotics and perturbation problems for 3D interfacial cracks

We present an asymptotic algorithm for analysis of a singularly perturbed problem in a domain containing an interfacial crack. The crack is assumed to be flat and its front, initially straight, is perturbed in the plane containing the crack. The aim of the work is to determine the asymptotic representation of the stress-intensity factors near the edge of the crack. Mathematically, the limit problem is reduced to the analysis of a matrix, 3×3, Wiener-Hopf problem, and its solution generates the “weight matrix-function” characterised by a special singular solution near the crack edge. The two-term asymptotic representation for the weight function components is required by the asymptotic algorithm, together with two-term asymptotics for stress components associated with the physical fields near the edge of the crack. The particular feature of the solution is the coupling between the normal opening mode (Mode-I), and the shear modes (Mode-II and Mode-III), and the oscillatory behaviour of certain stress components near the crack edge. Explicit asymptotic formulae for the stress-intensity factors are obtained at the edge of a “wavy crack” at an interface.