反平面V形切口塑性应力奇异性分析

论文提出了用插值矩阵法计算幂硬化塑性材料反平面V形切口和裂纹尖端区域的应力奇异性.首先在切口和裂纹尖端区域采用自尖端径向度量的渐近位移场假设,将其代入塑性全量理论的基本微分方程后,推导出包含应力奇异性特征指数和特征角函数的非线性常微分方程特征值问题.然后采用插值矩阵法迭代求解导出的控制方程,得到一般的塑性材料反平面V形切口和裂纹的前若干阶应力奇异阶和相应的特征角函数,该法的重要优点是以上求解的特征角函数和它们各阶导函数具有同阶精度,并且一次性地求出前若干阶特征对.同时,插值矩阵法计算量小,易于和其他方法联合使用,这些优点在后续求解尖端区域完全应力场非常优越.论文方法的计算结果与现有结果对照,发现吻合良好,表明了论文方法的有效性.     

Finite strain analysis of crack tip fields in incompressible hyperelastic solids loaded in plane stress

A new method is developed to determine the dominant asymptotic stress and deformation fields near the tip of a Mode-I traction free plane stress crack. The analysis is based on the fully nonlinear equilibrium theory of incompressible hyperelastic solids. We show that the dominant singularity of the near tip stress field is governed by the asymptotic solution of a linear second order ordinary differential equation. Our method is applicable to any hyperelastic material with a smooth work function that depends only on the trace of the Cauchy-Green tensor and is particularly useful for materials that exhibit severe strain hardening. We apply this method to study two types of soft materials: generalized neo-Hookean solids and a solid that hardens exponentially. For the generalized neo-Hookean solids, our method is able to resolve a difficulty in the previous work by Geubelle and Knauss (1994a). Our theoretical results are compared with finite element simulations.     

Antiplane problem of circular ARC interfacial rigid line inclusions

The antiplane problem of circular arc rigid line inclusions under antiplane concentrated force and longitudinal shear loading was dealt with. By using RiemannSchwaz‘s symmetry principle integrated with the singularity analysis of complex functions, the general solution of the problem and the closed form solutions for some important practical problems were presented. The stress distribution in the immediate vicinity of circular arc rigid line end was examined in detail. The results show that the singular stress fields near the rigid inclusion tip possess a square-root singularity similar to that for the corresponding crack problem under antiplane shear loading, but no oscillatory character. Furthermore, the stresses are found to depend on geometrical dimension, loading conditions and materials parameters. Some practical results concluded are in agreement with the previous solutions.     

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.     

Hypocycloidal Inclusions in Nonuniform Out-of-Plane Elasticity: Stress Singularity vs. Stress Reduction

Stress field solutions and Stress Intensity Factors (SIFs) are found for \(n\)-cusped hypocycloidal shaped voids and rigid inclusions in an infinite linear elastic plane subject to nonuniform remote antiplane loading, using complex potential and conformal mapping. It is shown that a void with hypocycloidal shape can lead to a higher SIF than that induced by a corresponding star-shaped crack; this is counter intuitive as the latter usually produces a more severe stress field in the material. Moreover, it is observed that when the order \(m\) of the polynomial governing the remote loading grows, the stress fields generated by the hypocycloidal-shaped void and the star-shaped crack tend to coincide, so that they become equivalent from the point of view of a failure analysis. Finally, special geometries and loading conditions are discovered for which there is no stress singularity at the inclusion cusps and where the stress is even reduced with respect to the case of the absence of the inclusion. The concept of Stress Reduction Factor (SRF) in the presence of a sharp wedge is therefore introduced, contrasting with the well-known definition of Stress Concentration Factor (SCF) in the presence of inclusions with smooth boundary. The results presented in this paper provide criteria that will help in the design of ultra strong composite materials, where stress singularities always promote failure. Furthermore, they will facilitate finding the special conditions where resistance can be optimized in the presence of inclusions with non-smooth boundary.     

复合材料切口应力奇性指数计算

提出一种计算广义平面应交状态下复合材料切口应力奇性指数的新方法.在切口尖端的位移幂级数渐近展开式被引入正交各向异性材料的物理方程后,将用位移表示的应力分量代入切口端部柱状邻域的线弹性理论控制方程,切口应力奇性指数的计算被转化为常微分方程组特征值的求解.采用插值矩阵法求解该常微分方程组,可一次性地获取切口尖端多阶应力奇性指数.本法适合平面和反平面应力场耦合或解耦的情形,并可退化计算裂纹或各向同性材料切口的应力奇性指数.算例表明,所提方法对分析复合材料切口应力奇性指数是一种准确有效的手段.     

The singular field near a sharp notch with mixed boundary conditions for hardening materials under longitudinal shear

The field behavior near a sharp notch tip with mixed homogeneous stress and displacement boundary conditions is examined for a power law hardening material. Using the hodograph transformation, the singularity and the angular distribution of the fields are determined. Special cases as those for linear elastic and perfect plastic materials are discussed.     

Wedge indentation into elastic-plastic single crystals, 1: Asymptotic fields for nearly-flat wedge

Asymptotic stress and deformation fields under the contact point singularities of a nearly-flat wedge indenter and of a flat punch are derived for elastic ideally-plastic single crystals with three effective in-plane slip systems that admit a plane strain deformation state. Face-centered cubic (FCC), body-centered cubic (BCC), and hexagonal-close packed (HCP) crystals are considered. The asymptotic fields for the flat punch are analogous to those at the tip of a stationary crack, so a potential solution is that the deformation field consists entirely of angular constant stress plastic sectors separated by rays of plastic deformation across which stresses change discontinuously. The asymptotic fields for a nearly-flat wedge indenter are analogous to those of a quasistatically growing crack tip fields in that stress discontinuities can not exist across sector boundaries. Hence, the asymptotic fields under the contact point singularities of a nearly-flat wedge indenter are significantly different than those under a flat punch. A family of solutions is derived that consists entirely of elastically deforming angular sectors separated by rays of plastic deformation across which the stress state is continuous. Such a solution can be found for FCC and BCC crystals, but it is shown that the asymptotic fields for HCP crystals must include at least one angular constant stress plastic sector. The structure of such fields is important because they play a significant role in the establishment of the overall fields under a wedge indenter in a single crystal. Numerical simulations—discussed in detail in a companion paper—of the stress and deformation fields under the contact point singularity of a wedge indenter for a FCC crystal possess the salient features of the analytical solution.     

脆性损伤材料裂纹尖端的局部解

本文利用Mazars和Lemaitre提出的混凝土脆性损伤模型,求得了裂纹尖端应力、应变及损伤的局部解.对手Ⅲ型及不可压缩平面应变Ⅰ型裂纹,其尖端场的构造和理想塑性材料相类似.指出由于丧失了应力全连续条件,从而损伤边界不能由局部解定出.     

Crack-tip field on mode Ⅱ interface crack of double dissimilar orthotropic composite materials

Two systems of non-homogeneous linear equations with 8 unknowns are obtained.This is done by introducing two stress functions containing 16 undetermined coefficients and two real stress singularity exponents with the help of boundary conditions.By solving the above systems of non-homogeneous linear equations,the two real stress singularity exponents can be determined when the double material parameters meet certain conditions.The expression of the stress function and all coefficients are obtained based on the uniqueness theorem of limit.By substituting these parameters into the corresponding mechanics equations,theoretical solutions to the stress intensity factor,the stress field and the displacement field near the crack tip of each material can be obtained when both discriminants of the characteristic equations are less than zero.Stress and displacement near the crack tip show mixed crack characteristics without stress oscillation and crack surface overlapping.As an example,when the two orthotropic materials are the same,the stress singularity exponent,the stress intensity factor,and expressions for the stress and the displacement fields of the orthotropic single materials can be derived.     

Analytical representation of the non-square-root singular stress field at a finite angle sharp notch

The stress field near the tip of a finite angle sharp notch is singular. However, unlike a crack, the order of the singularity at the notch tip is less than one-half. Under tensile loading, such a singularity is characterized by a generalized stress intensity factor which is analogous to the mode I stress intensity factor used in fracture mechanics, but which has order less than one-half. By using a cohesive zone model for a notional crack emanating from the notch tip, we relate the critical value of the generalized stress intensity factor to the fracture toughness. The results show that this relation depends not only on the notch angle, but also on the maximum stress of the cohesive zone model. As expected the dependence on that maximum stress vanishes as the notch angle approaches zero. The results of this analysis compare very well with a numerical (finite element) analysis in the literature. For mixed-mode loading the limits of applicability of using a mode I failure criterion are explored.     

Antiplane interaction of line crack with an arbitrarily located elliptical inclusion

The antiplane problem of the interaction between a main crack and an arbitrarily located elastic elliptical inclusion near its tip is addressed in the current study. The analysis is based on the use of the complex potentials for the antiplane problem, Laurent series expansion method and an appropriate superposition scheme. The stress intensity factor at the main crack is obtained in a general series form. Explicit asymptotic solutions are also derived by using a perturbation technique and retaining the leading order terms in series expansion. The present solutions are shown to coincide with the Taylor expansion of exact solutions for special cases available in the literature. Discussed are changes in the crack tip stress intensity which can be enhanced or suppressed depending on the location of the elliptical inclusion. The explicit solutions provided herein are well suited for the further quantitative analysis of toughening mechanisms in ceramic composite materials.     

The near crack tip fields of stress and damage for concrete

The crack tip fields of stress, strain and damage for concrete under both antiplane shear and plane strain conditions are investigated based on the damage model proposed by Mazars and Lemaitre [2]. The structures of near tip fields obtained are similar to those for an elastic-perfectly-plastic material. It has been found that damage boundaries can not be determined by the near-tip analysis due to the discontinuities of stresses on the damage boundaries induced by the damage model used in the present paper.The Project is Supported by National Natural Science Fundation of China.     

STRESS SINGULARITY ANALYSIS OF MULTI-MATERIAL WEDGES UNDER ANTIPLANE DEFORMATION

With the help of the coordinate transformation technique, the symplectic dual solving system is developed for multi-material wedges under antiplane deformation. A virtue of present method is that the compatibility conditions at interfaces of a multi-material wedge are expressed directly by the dual variables, therefore the governing equation of eigenvalue can be derived easily even with the increase of the material number. Then, stress singularity on multi-material wedges under antiplane deformation is investigated, and some solutions can be presented to show the validity of the method. Simultaneously, an interesting phenomenon is found and proved strictly that one of the singularities of a special five-material wedge is independent of the crack direction.     

Antiplane electro-mechanical field of a piezoelectric finite wedge under shear loading and at fixed-grounded boundary conditions

This paper presents the general solutions of antiplane electro-mechanical field solutions for a piezoelectric finite wedge subjected to a pair of concentrated forces and free charges. The boundary conditions on the circular segment are considered as fixed and grounded. Employing the finite Mellin transform method, the stress and electrical displacement at all fields of the piezoelectric finite wedge are derived analytically. In addition, the singularity orders and intensity factors of stress and electrical displacement can also be obtained. These parameters can be applied to examine the fracture behavior of the wedge structure. After being reduced to the problem of an antiplane edge crack or an infinite wedge in a piezoelectric medium, the results compare well with those of previous studies.     

Interface crack problems for mode II of double dissimilar orthotropic composite materials

The fracture problems near the similar orthotropic composite materials are interface crack tip for mode Ⅱ of double disstudied. The mechanical models of interface crack for mode Ⅱ are given. By translating the governing equations into the generalized hi-harmonic equations, the stress functions containing two stress singularity exponents are derived with the help of a complex function method. Based on the boundary conditions, a system of non-homogeneous linear equations is found. Two real stress singularity exponents are determined be solving this system under appropriate conditions about bimaterial engineering parameters. According to the uniqueness theorem of limit, both the formulae of stress intensity factors and theoretical solutions of stress field near the interface crack tip are derived. When the two orthotropic materials are the same, the stress singularity exponents, stress intensity factors and stresses for mode II crack of the orthotropic single material are obtained.     

Interface crack problems for mode Ⅱ of double dissimilar orthotropic composite materials

The fracture problems near the interface crack tip for mode Ⅱ of double dissimilar orthotropic composite materials are studied. The mechanical models of interface crack for mode Ⅱ are given. By translating the governing equations into the generalized bi-harmonic equations,the stress functions containing two stress singularity exponents are derived with the help of a complex function method. Based on the boundary conditions,a system of non-homogeneous linear equations is found. Two real stress singularity exponents are determined be solving this system under appropriate conditions about himaterial engineering parameters. According to the uniqueness theorem of limit,both the formulae of stress intensity factors and theoretical solutions of stress field near the interface crack tip are derived. When the two orthotropic materials are the same,the stress singularity exponents,stress intensity factors and stresses for mode Ⅱ crack of the orthotropic single material are obtained.     

Some problems in the antiplane shear deformation of bi-material wedges

The antiplane shear deformation of a bi-material wedge with finite radius is studied in this paper. Depending upon the boundary condition prescribed on the circular segment of the wedge, traction or displacement, two problems are analyzed. In each problem two different cases of boundary conditions on the radial edges of the composite wedge are considered. The radial boundary data are: traction–displacement and traction–traction. The solution of governing differential equations is accomplished by means of finite Mellin transforms. The closed form solutions are obtained for displacement and stress fields in the entire domain. The geometric singularities of stress fields are observed to be dependent on material property, in general. However, in the special case of equal apex angles in the traction–traction problem, this dependency ceases to exist and the geometric singularity shows dependency only upon the apex angle. A result which is in agreement with that cited in the literature for bi-material wedges with infinite radii. In part II of the paper, Antiplane shear deformation of bi-material circular media containing an interfacial edge crack is considered. As a special case of bi-material wedges studied in part I of the paper, explicit expressions are derived for the stress intensity factor at the tip of an edge crack lying at the interface of the bi-material media. It is seen that in general, the stress intensity factor is a function of material property. However, in special cases of traction–traction problem, i.e., similar materials and also equal apex angles, the stress intensity factor becomes independent of material property and the result coincides with the results in the literature.     

Asymptotic mechanical fields at the tip of a mode I crack in rubber-like solids

Asymptotic analyses of the mechanical fields in front of stationary and propagating cracks facilitate the understanding of the mechanical and physical state in front of crack tips, and they enable prediction of crack growth and failure. Furthermore, efficient modelling of arbitrary crack growth by use of XFEM (extended finite element method) requires accurate knowledge of the asymptotic crack tip fields. In the present work, we perform an asymptotic analysis of the mechanical fields in the vicinity of a propagating mode I crack in rubber. Plane deformation is assumed, and the material model is based on the Langevin function, which accounts for the finite extensibility of polymer chains. The Langevin function is approximated by a polynomial, and only the term of the highest order contributes to the asymptotic solution. The crack is predicted to adopt a wedge-like shape, i.e. the crack faces will be straight lines. The angle of the wedge and the order of the stress singularity depend on the hardening of the strain energy function. The present analysis shows that in materials with a significant hardening, the inertia term in the equations of motion becomes negligible in the asymptotic analysis. Hence, there is no upper theoretical limit to the crack speed.     

Permeable cracks between two dissimilar piezoelectric materials

An interfacial crack with electrically permeable surfaces between two dissimilar piezoelectric ceramics under electromechanical loading is investigated. An exact expression for singular stress and electric fields near the tip of a permeable crack between two dissimilar anisotropic piezoelectric media are obtained. The interfacial crack-tip fields are shown to consist of both an inverse square root singularity and a pair of oscillatory singularities. It is found that the singular fields near the permeable interfacial crack tip are uniquely characterized by the real valued stress intensity factors proposed in this paper. The energy release rate is obtained in terms of the stress intensity factors. The exact solution of stress and electric fields for a finite interfacial crack problem is also derived.