Analysis on the cohesive stress at half infinite crack tip

The nonlinear fracture behavior of quasi-brittle materials is closely related with the cohesive force distribution of fracture process zone at crack tip. Based on fracture character of quasi-brittle materials, a mechanical analysis model of half infinite crack with cohesive stress is presented. A pair of integral equations is established according to the superposition principle of crack opening displacement in solids, and the fictitious adhesive stress is unknown function . The properties of integral equations are analyzed, and the series function expression of cohesive stress is certified. By means of the data of actual crack opening displacement, two approaches to gain the cohesive stress distribution are proposed through resolving algebra equation. They are the integral transformation method for continuous displacement of actual crack opening, and the least square method for the discrete data of crack opening displacement. The calculation examples of two approaches and associated discussions are give     

Boundary element analysis of interaction between an elastic rectangular inclusion and a crack

The interaction between an elastic rectangular inclusion and a kinked crack inan infinite elastic body was considered by using boundary element method. The new complexboundary integral equations were derived. By introducing a complex unknown function H(t)related to the interface displacement density and traction and applying integration by parts,the traction continuous condition was satisfied automatically. Only one complex boundaryintegral equation was obtained on interface and involves only singularity of order l/ r. Toverify the validity and effectiveness of the present boundary element method, some typicalexamples were calculated. The obtained results show that the crack stress intensity factorsdecrease as the shear modulus of inclusion increases. Thus, the crack propagation is easiernear a softer inclusion and the harder inclusion is helpful for crack arrest.     

An effective boundary element method for analysis of crack problems in a plane elastic plate

A simple and effective boundary element method for stress intensity factor calculation for crack problems in a plane elastic plate is presented. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements proposed by YAN Xiangqiao. In the boundary element implementation the left or the right crack-tip displacement discontinuity element was placed locally at the corresponding left or right each crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Test examples (i. e. , a center crack in an infinite plate under tension, a circular hole and a crack in an infinite plate under tension) are included to illustrate that the numerical approach is very simple and accurate for stress intensity factor calculation of plane elasticity crack problems. In addition, specifically, the stress intensity factors of branching cracks emanating from a square hole in a rectangular plate under biaxial loads were analysed. These numerical results indicate the present numerical approach is very effective for calculating stress intensity factors of complex cracks in a 2-D finite body, and are used to reveal the effect of the biaxial loads and the cracked body geometry on stress intensity factors.     

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.     

APPLICATION OF THE DIGITAL MOIRE METHOD IN FRACTURE ANALYSIS OF A CRACKED RUBBER SHEET

This paper deals with the mode I crack problem of a cracked rubber sheet under plane stress condition using the delicate digital moire technique. Through the four step phase-shifting method of automated fringe analysis, the displacement fields in the Cartesian coordinate system are given. By the coordinate-transform equation, the radial and circular displacement distributions in the polar coordinate system are obtained. The displacement isoline distributions and displacement vector distributions near the crack tip are discussed. The strain isoline distributions near the crack tip are also analyzed in this paper. Finally, the distribution rules for the mechanical fields near the crack tip are discussed with the sector division method.     

THE ELASTIC FIELD INDUCED BY A HEMISPHERICAL INCLUSION IN THE HALF—SPACE

The elastic field induced by a hemispherical inclusion with uniform eigeustralns in asemi-infinite elastic medium is solved by using the Green‘s function method and series expansion tech-nique. The exact solutions axe presented for the displacement and stress fields which can be expressedby complete elliptic integrals of the first, second, and third kinds and hypergeometric functions. Thepresent method can be used to determine the corresponding elastic fields when the shape of the inclusionis a spherical crown or a spherical segment. Finally, numerical results axe given for the displacementand stress fields along the axis of symmetry (x3-axis).     

CALCULATION OF CONFORMAL MAPPING FUNCTION OF THE TOOTH PROFILE ON EVOLVENT GEAR BY COMPUTER

The exact solution of the stress and the displacement on spur gear can be gained by the method of conformal mapping of complex variables in plane elasticity. But it is difficult to get the comformal mapping function for the tooth profile with different parameters. They used to be obtained through trial method. This is time-consuming and expensive, In this paper a computer program is drawn up for the conformal mapping function. A large amount of calculation has gi yen proof of its success and of the precision of the mapping function. thereby the main obstacle has been removed for the practical application of the conformal mapping method to the stress mid displacement of tooth on evolvent gear.     

Torsional impact response of a penny-shaped crack in a functional graded strip

The torsional impact response of a penny-shaped crack in a nonhomogeneous strip is considered. The shear modulus is assumed to be functionally graded such that the mathematics is tractable. Laplace and Hankel transforms were used to reduce the problem to solving a Fredholm integral equation. The crack tip stress field is obtained by considering the asymptotic behavior of Bessel function. Explicit expressions of both the dynamic stress intensity factor and the energy density factor were derived. And it is shown that, as crack driving force, they are equivalent for the present crack problem. Investigated are the effects of material nonhomogeneity and strip‘s highness on the dynamic fracture behavior.Numerical results reveal that the peak of the dynamic stress intensity factor can be suppressed by increasing the nonhomogeneity parameter of the shear modulus, and that the dynamic behavior varies little with the adjusting of the strip‘ s highness.     

THE THREE-DIMENSIONAL STRESS INTENSITY FACTOR UNDER MOVING LOADS ON THE CRACK FACES

The dynamic stress intensity factor history for a half plane crack in anotherwise unbounded elastic body,with the crack faces subjected to a tractiondistribution consisting of two pairs of combined mode point loads that move in adirection perpendicular to the crack edge is considered.The analytic expression for thecombined mode stress intensity factors as a function of time for any point along thecrack edge is obtained.The method of solution is based on the application of integraltransform together with the Wiener-Hopf technique and the Cagniard-de Hoop method.Some features of the solution are discussed and graphical results for various point loadspeeds are presented.     

压电材料平面问题的尖端场和应力强度因子的求解

利用Lekhnitskii理论和Stroh理论的相互联系,把已知的基于Lekhnitskii理论平面应变结果转化为Stroh理论形式的结果,直接获得Stroh公式中A,B的显式表达式,此方法可扩展到平面应力情况,然后导出压电材料平面应变问题的尖端场Williams形式的展开式,采用半权函数法计算有限大压电体平面问题应力和电位移强度因子。对无穷大板含中心裂纹的情况下本文结果和已有结果进行了比较,表明本文方法得到的结果精度可靠。本文方法的最大优点是可以求解有限压电体的应力强度因子,并且需要的单元少,精度高,实用性好。     

A Crack with an Electric Displacement Saturation Zone in an Electrostrictive Material

A crack with an electric displacement saturation zone in an electrostrictive material under purely electric loading is analyzed. A strip saturation model is here employed to investigate the effect of the electrical polarization saturation on electric fields and elastic fields. A closed form solution of electric fields and elastic fields for the crack with the strip saturation zone is obtained by using the complex function theory. It is found that the K _I -dominant region is very small compared to the strip saturation zone. The generalized Dugdale zone model is also employed in order to investigate the effect of the saturation zone shape on the stress intensity factor. Using the body force analogy, the stress intensity factor for the asymptotic problem of a crack with an elliptical saturation zone is evaluated numerically.     

A MODE Ⅱ CRACK IN A TWO-DIMENSIONAL OCTAGONAL QUASICRYSTALS

The present study develops the fracture theory for a two-dimensional octagonal quasicrystals. The exact analytic solution of a Mode Ⅱ Griffith crack in the material was obtained by using the Fourier transform and dual integral equations theory, then the displacement and stress fields, stress intensity factor and strain energy release rate were determined, the physical sense of the results relative to phason and the difference between mechanical behaviors of the crack problem in crystal and quasicrystal were figured out. These provide important information for studying the deformation and fracture of the new solid phase.     

Scattering of anti-plane shear waves by a single crack in an unbounded trasversely isotropic electro-magneto-elastic medium

A theoretical treatment of the scattering of anti-plane shear (SH) waves is provided by a single crack in an unbounded transversely isotropic electro-magneto-elastic medium. Based on the differential equations of equilibrium, electric displacement and magnetic induction intensity differential equations, the governing equations for SH waves were obtained. By means of a linear transform, the governing equations were reduced to one Helmholtz and two Laplace equations. The Cauchy singular integral equations were gained by making use of Fourier transform and adopting electro-magneto impermeable boundary conditions. The closed form expression for the resulting stress intensity factor at the crack was achieved by solving the appropriate singular integral equations using Chebyshev polynomial. Typical examples are provided to show the loading frequency upon the local stress fields around the crack tips. The study reveals the importance of the electro-magneto-mechanical coupling terms upon the resulting dynamic stress intensity factor. Contributed by SHEN Ya-peng Foundation item: the National Natural Science Foundation of China (10132010, 50135030) Biographies: DU Jian-ke (1970∼)     

An analysis on three-dimensional effects near the tip of a crack in an elastic plate

An infinite plate containing a finite through crack under tensile loading is analysed by Fourier transform based on the Kane-Mindlin kinematic assumptions for the quasi-three-dimensional deformation of plates in extension. The asymptotic expressions of stress and displacement fields near the crack tip, the variation of the stress intensity factor with the plate-thickness and the three-dimensional deformation zone near the crack tip are investigated. The results of the analysis show that, (a) the crack-tip stress and displacement fields accounting for the plate-thickness effects are different from the plane stress solutions and this is true even for extremely small parameter (=1–vh/6 a). In a very small region near the crack tip, plane strain solutions prevail; (b) the ratio of the stress intensity factor K_I to the corresponding plane stress one K_I, K_I/K _I ^o , approaches 1/(1–v²) as tends to zero; (c) plane stress solutions can give satisfactory results for points a distance from the crack tip greater than about three-fourths of the plate-thickness; (d) the linear elastic result for the zone of three-dimensional effects is approximately valid for an elasto-plastic material with linear strain-hardening when the plastic tangential mudulus E_t is not very small.The Project Supported by National Natural Science Foundation of China.     

Crack problems of piezoelectric materials

This paper presents an analysis of crack problems in homogeneous piezoelectrics or on the interfaces between two dissimilar piezoelectric materials based on the continuity of normal electric displacement and electric potential across the crack faces. The explicit analytic solutions are obtained for a single crack in piezoelectrics or on the interfaces of piezoelectric bimaterials. A class of boundary problems involving many cracks is also solved. For homogeneous materials it is found that the normal electric displacementD ₂ induced by the crack is constant along the crack faces which depends only on the applied remote stress field. Within the crack slit, the electric fields induced by the crack are also constant and not affected by the applied electric field. For the bimaterials with realH, the normal electric displacementD ₂ is constant along the crack faces and electric fieldE ₂ has the singularity ahead of the crack tip and a jump across the interface. The project is supported by the National Natural Science Foundation of China(No. 19704100) and the Natural Science Foundation of Chinese Academy of Sciences(No. KJ951-1-201).     

Uniformly valid asymptotic solutions of the nonlinear unsymmetrical bending for orthotropic rectangular thin plate of four clamped edges with variable thickness

By using “the method of modified two-variable”,“the method of mixing perturbation ”and introducing four small parameters, the problem of the nonlinear unsymmetrical bending for orthotropic rectangular thin plate with linear variable thickness is studied. And the uniformly valid asymptotic solution of Nth- order for ε1 and Mth- order for ε2 of the deflection functions and stress function are obtained.     

A study ofJ-integral of the orthotropic composite material

The relation between J-integral near model I crack tip in the orthotropic plate and displacement derivative is derived in this paper. Meanwhile, the relation between stress intensity factor K _I and displacement is also given in this paper. With sticking film moire interferometry method, the three-point bending beam is tested, thus the values of J-integral and K _I can be obtained from the displacement field, and then the truth of relation formula between J-integral and K _I in the orthotropic composite materials is experimentally verified.     

Mixed Mode Dynamic Fracture in Particulate Reinforced Functionally Graded Materials

A detailed analytical and experimental investigation is presented to understand the dynamic fracture behavior of functionally graded materials (FGMs) under mode I and mixed mode loading conditions. Crack-tip stress, strain and displacement fields for a mixed mode crack propagating at an angle from the direction of property gradation were obtained through an asymptotic analysis coupled with a displacement potential approach. This was followed by a comprehensive series of experiments to gain further insight into the behavior of propagating cracks in FGMs. Dynamic photoelasticity coupled with high-speed photography was used to obtain crack tip velocities and dynamic stress fields around the propagating cracks. Birefringent coatings were used to conduct the photoelastic study due to the opaqueness of the FGMs. Dynamic fracture experiments were performed using different specimen geometries to develop a dynamic constitutive fracture relationship between the mode I dynamic stress intensity factor (K _ID ) and crack-tip velocity ( ) for FGMs with the crack moving in the direction of increasing fracture toughness. A similar -K _ID relation was also obtained for matrix material (polyester) for comparison purposes. The results obtained show that crack propagation velocities in FGMs were about 80% higher than the polyester matrix. Crack arrest toughness was found to be about 10% lower than the value of local fracture toughness in FGMs.     

Penny-shaped crack in transversely isotropic piezoelectric materials

Using a method of potential functions introduced successively to integrate the field equations of three-dimensional problems for transversely isotropic piezoelectric materials, we obtain the so-called general solution in which the displacement components and electric potential functions are represented by a singular function satisfying some special partial differential equations of 6th order. In order to analyse the mechanical-electric coupling behaviour of penny-shaped crack for above materials, another form of the general solution is obtained under cylindrical coordinate system by introducing three quasi-harmonic functions into the general equations obtained above. It is shown that both the two forms of the general solutions are complete. Furthermore, the mechanical-electric coupling behaviour of penny-shaped crack in transversely isotropic piezoelectric media is analysed under axisymmetric tensile loading case, and the crack-tip stress field and electric displacement field are obtained. The results show that the stress and the electric displacement components near the crack tip have (r ^−1/2) singularity. The project supported by the Natural Science Foundation of Shaanxi Province, China