Crack growth under plane stress conditions in an elastic perfectly-plastic material

Mode-I crack growth under conditions of generalized plane stress has been investigated. It has been assumed that near the plane of the crack in the loading zone, the simple stress components corresponding to a central fan field maintain validity up to the elastic-plastic boundary. By the use of expansions of the particle velocities in the coordinate y, and by matching of the relevant stress components and particle velocities to the dominant terms of appropriate elastic fields at the elastic-plastic boundary, a complete solution has been obtained for ε_y in the plane of the crack. The solution applies from the propagating crack tip up to the moving elastic-plastic boundary. The strain fields for a self-similar crack nucleating at a point and for steady-state propagation of a crack have been considered as special cases.     

Compressive response and failure of braided textile composites: Part 1—experiments

Experimental results obtained by examining the planar biaxial compression/tension response of carbon 2D triaxial braided composites (2DTBC) are reported in this paper. These experiments were motivated by a need to examine the failure of 2DTBC in a state of stress that would be similar to what is experienced by the walls of a tubular member under compressive crush loads. Results obtained from a series of biaxial tests that were conducted with different proportional displacement loading ratio combinations of compression and tension are reported. In all cases, the dominant failure mechanism under such a stress state is the buckling of the bias and axial tows within the composite. Full field surface displacement data is acquired concurrently during all biaxial and some uniaxial tests using the technique of digital speckle photography. Digital images of the specimen surface that is illuminated with a He-Ne laser are acquired at discrete time intervals during the loading history using a high-resolution digital camera. These images are stored and analyzed to obtain the incremental inplane surface displacement field, Δu(x,y) and Δv(x,y). From these, the incremental inplane surface strains Δε_x, Δε_y and Δγ_xy are obtained by numerical differentiation. The present paper, which is the first in a two part series, is devoted to the biaxial experimental results pertaining to 2DTBC failure.     

Stress singularity due to traction discontinuity on an anisotropic elastic half-plane

In a half-plane problem with x₁ paralleling with the straight boundary and x₂ pointing into the medium, the stress components on the boundary whose acting plane is perpendicular to x₁ direction may be denoted by t₁ = [σ₁₁, σ₁₂, σ₁₃]^T. Stress components σ₁₁ and σ₁₃ are of more interests since σ₁₂ is completely determined by the boundary conditions. For isotropic materials, it is known that under uniform normal loading σ₁₁ is constant in the loaded region and vanishes in the unloaded part. Under uniform shear loading, σ₁₁ will have a logarithmic singularity at the end points of shear loading. In this paper, the behavior of the stress components σ₁₁ and σ₁₃ induced by traction-discontinuity on general anisotropic elastic surfaces is studied. By analyzing the problem of uniform tractions applied on the half-plane boundary over a finite loaded region, exact expressions of the stress components σ₁₁ and σ₁₃ are obtained which reveal that these components consist of in general a constant term and a logarithmic term in the loaded region, while only a logarithmic term exists in unloaded region. Whether the constant term or the logarithmic term will appear or not completely depends on what values of the elements of matrices Ω and Γ will take for a material under consideration. Elements for both matrices are expressed explicitly in terms of elastic stiffness. Results for monoclinic and orthotropic materials are all deduced. The isotropic material is a special case of the present results.     

A penny-shaped crack in transversely isotropic magneto-electro-thermo-elastic medium subjected to uniform symmetric heat flux

The problem of a penny-shaped crack subjected to symmetric uniform heat flux in an infinite transversely isotropic magneto-electro-thermo-elastic medium is investigated. The exact solution in the full space is in terms of line integrals and the exact solution in the crack plane also is obtained. Although we start our derivations with magneto-electro-thermo-elastic, the solution presented in this paper is also applicable for linear transversely isotropic thermopiezoelectric, thermomagnetoelastic,thermoelastic materials (see Appendix E). The solution in the crack plane, which shows a great agreement with the solution for a transversely isotropic medium obtained by Tsai (1983), indicates that _σx,_σy,_Dx,_Dy,_Bx${\sigma }_x\text{,}{\sigma }_y\text{,}{D}_x\text{,}{D}_y\text{,}{B}_x$, and _By$B_y$ along the crack rim are of the same singularity of the normal stress or its equivalent quantities. To illustrate how the applied symmetric heat fluxes affect the whole fields, a numerical example is also given.     

Experimental GH singularity field around growing crack-tip

The displacement fieldsu _x ,u _y at growing crack tip of LY12-M specimens with double edge cracks are measured using moire method. The experimental singularity fields are compared with GH theoretical field [12–14]. The size and shape of the experimental GH singularity fields are obtained. The error in both the experimental and theoretical evaluations is controlled within ±10%. The experiments show that there is singularity dominant around a growing crack tip. The shape of this dominant region ranges from butterfly wing to oblate and circular. Inside GH-field, there is a 3-D deformed damage zone where no GH singularity exists. The project suppoted by National Natural Science Foundation of China     

Finite deformation analysis of COD,J-integral and crack tip intense strain region in plane strain large-scale yielding

A finite element analysis was performed to simulate crack tip blunting and the development of the intense strain region in a small compact tension specimen (0.4 T CT) of SA533B-1 under plane strain large-scale yielding, with the condition of large-geometry change around the crack tip taken into consideration. The region where the equivalent plastic strain $\text{}\text{?}$g3_p is greater than 0.15 was defined as the intense strain region, which corresponded to the recrystallized-etched zone delineated experimentally around the blunting crack tip. The development of the intense strain region was discussed as a function of the J-integral and the crack opening displacement. A linear relationship was obtained between the plastic work W_p dissipated within the intense strain region and (J/σ_y)² or b², where b is the crack opening displacement, defined as the separation of the two points at which the boundary of the intense strain region surrounding the crack tip intersects with the free surfaces of the crack.     

Crack separation energy-rates for the DBCS model under biaxial modes of loading

A Modified version of the Dugdale-Bilby-Cottrell-Swinden (DBCS) model simulating the effect of plasticity at the tip of a crack in an infinite region was used by kfouri and rice (1978) to calculate the crack separation energy-rate G^Δ corresponding to a finite crack growth step Δa during plane strain mode I crack extension. The loading consisted of a remote uniaxial tension σp applied normally to the plane of the crack. Using Rice's path-independent integral J to characterize the applied load in the crack tip region, and assuming the length R of the crack tip plastic zone to be small compared with the length of the crack, an analytical expression was derived relating the ratios (J/G^Δ) and (2a/R) for small values of (2a/R). The analytical solution was incomplete in itself in that the value assumed in the plastic region of the DBCS model for the normal stress Y acting on the extending crack surfaces was the product of the yield stress in uniaxial tension σ_Y and an unknown parameter C, the value of which depends on the effect of the local hydrostatic stresses in the case of plane strain conditions. The analytical solution was compared with a numerical solution obtained from a plane strain elastic-plastic finite element analysis on a centre-cracked plate (CCP) of material obeying the von Mises yield criterion. The value used for the yield stress was 310 MN/m² and moderate isotropic linear hardening was applied with a tangent modulus of 4830 MN/m². A uniaxial tension σp was applied on the two appropriate sides of the plate. The comparisons showed that the analytical and finite element solutions were mutually consistent and they enabled the value of C to be established at 1.91. In the present paper similar comparisons are made between the analytical solution and the finite element solutions for the CCP of the same material under different biaxial modes of loading. By assuming further that the form of the analytical solution does not depend on the details of the geometry and of the loading at remote boundaries, a comparison has also been made with the results of a finite element analysis on a compact tension specimen (CTS) made of the same material as the CCP. The different values of C obtained in each case are consistent with investigations by other authors on the effect of load biaxiality on crack tip plasticity.     

A stress analytical solution for Mode III crack within modified gradient elasticity

The strain gradient exists near a crack tip may significantly influence the near-tip stress field. In this paper, the strain gradient and the internal length scales are introduced into the basic equations of mode III crack by the modified gradient elasticity (MGE). By using a complex function approach, the analytical solution of stress fields for mode III crack problem is derived within MGE. When the internal length scales vanish, the stress fields can be simplified to the stress fields of classical linear elastic fracture mechanics. The results show that the singularity of the shear stress is made up of two parts, r^−1/2 part and r^−3/2 part, and the sign of the stress σ_yz changes. With the increase of l_x, the peak value of σ_yz decrease and its location moves farther from the fracture vertex. The influence of strain gradient for mode III crack problem cannot be ignored.     

In-plane and out-of-plane constraint effects under pressurized thermal shocks

Transferability of fracture toughness data obtained on small scale specimens to a full-scale cracked structure is one of the key issues in integrity assessment of engineering structures. In order to transfer fracture toughness under different constraints, both in-plane and out-of-plane constraint effect should be considered for the specimens and structures. In this paper both in-plane and out-of-plane constraint effects of a crack in a reference reactor pressure vessel (RPV) subjected to pressurized thermal shocks (PTSs) are analyzed by two-parameter and three-parameter methods. The comparison between elastic and elastic–plastic analysis shows that the constraint effect varies with the material property. T₁₁ (the second term of William’s extension acting parallel to the crack plane) generally displays a reversed relation to the stress intensity factor (SIF) with the transient time, which indicates that the loading (SIF) plays an important role on the in-plane constraint effect. The thickness at the crack tip contributes more than the loading to the out-of-plane constraint, such that T₃₃ (the second term of William’s extension acting along the thickness) displays a similar relation to ε₃₃ (strain along the thickness direction) and a different relation to T₁₁ during the transient. The results demonstrate that both in-plane and out-of-plane constraint effect should be analyzed separately in order to describe precisely the stress distribution ahead of the crack tip.     

An experimental investigation of crack-path directional stability

Using the criterion that a crack will extend perpendicular to the maximum circumferential stress,σ _θ, we show that the directional stability of crack growth is governed by the location of microcrack initiation ahead of the crack tip. At distances greater than a geometrical radiusr _o, the maximum value ofσ _θ deviates from the position of symmetry. Thus, if we assume that the physical processes involved in fracture lead to crack initiation at a distancer _c ahead of the crack tip, the criterion for directional stability isr _o>r _c. Experimental and theoretical values ofr _o verify that, asr _o becomes small, the crack's directional stability deteriorates. Observing that a lengthwise compressive stress increasesr _o, a center-cracked specimen was developed which allows the application of controlled lengthwise compression independently of the opening-mode load. A detailed photoelastic analysis of the specimen has provided the value ofr _o as a function of the crack length. The value ofr _o is then compared with the expected microcrack initiation distances in ductile fracture. By applying sufficient lengthwise compression, we are able to make the crack grow straight and obtain numerous data points from this specimen which would otherwise be directionally unstable. The results indicate that, as the total lengthwise tensile stress at the crack tip increases, the fracture toughness also increases. Using this information we can then adjustK _Ic for zero lengthwise loading and obtain a geometry independent fracture toughness.     

Effects of pre-stress on crack-tip fields in elastic,incompressible solids

A closed-form asymptotic solution is provided for velocity fields and the nominal stress rates near the tip of a stationary crack in a homogeneously pre-stressed configuration of a nonlinear elastic, incompressible material. In particular, a biaxial pre-stress is assumed with stress axes parallel and orthogonal to the crack faces. Two boundary conditions are considered on the crack faces, namely a constant pressure or a constant dead loading, both preserving an homogeneous ground state. Starting from this configuration, small superimposed Mode I or Mode II deformations are solved, in the framework of Biot's incremental theory of elasticity. In this way a definition of an incremental stress intensity factor is introduced, slightly different for pressure or dead loading conditions on crack faces. Specific examples are finally developed for various hyperelastic materials, including the J₂-deformation theory of plasticity. The presence of pre-stress is shown to strongly influence the angular variation of the asymptotic crack-tip fields, even if the nominal stress rate displays a square root singularity as in the infinitesimal theory. Relationships between the solution with shear band formation at the crack tip and instability of the crack surfaces are given in evidence.     

In Plane Displacement Formulation for Finite Cracked Plates Under Mode I Using Grid Method and Finite Element Analysis

A grid method is used to experimentally determine the in-plane displacement fields around a crack tip in a Single-Edge-Notch (SEN) tensile polyurethane specimen. Horizontal displacement u _x-exp and vertical displacement u _y-exp are expressed as functions of circular coordinates centred on the crack tip. These are compared with the approximate solutions of linear elastic fracture mechanics with a view to studying the applicability to polymers. The results show that this solution is not in agreement with the experiments at the focused on the vicinity of a crack tip. Taking this into account, an FEA program is developed with CAST3M for the purpose of comparing the experimental displacements and the numerical data. New formulations of displacements u _x and u _y are then developed. These formulations are derived from the principle of superposition and based on Arakawa’s formulation. With the displacement gradients obtained from the FEA and the new formulations, the determination of J-integrals is found to be in very good agreement with those derived from numerical calculation. Consequently, the proposed formulations can give displacement fields compatible with the J-integral calculation for the region near the crack tip. An application based on an experimental test is proposed to evaluate the performances of the proposed formulations.     

Asymptotic crack-tip fields for dynamic crack growth in compressive power-law hardening materials

The effect of material compressibility on the stress and strain fields for a mode-I crack propagating steadily in a power-law hardening material is investigated under plane strain conditions. The plastic deformation of materials is characterized by the J₂ flow theory within the framework of isotropic hardening and infinitesimal displacement gradient. The asymptotic solutions developed by the present authors [Zhu, X.K., Hwang K.C., 2002. Dynamic crack-tip field for tensile cracks propagating in power-law hardening materials. International Journal of Fracture 115, 323–342] for incompressible hardening materials are extended in this work to the compressible hardening materials. The results show that all stresses, strains, and particle velocities in the asymptotic fields are fully continuous and bounded without elastic unloading near the dynamic crack tip. The stress field contains two free parameters σ_eq0 and s₃₃₀ that cannot be determined in the asymptotic analysis, and can be determined from the full-field solutions. For the given values of σ_eq0 and s₃₃₀, all field quantities around the crack tip are determined through numerical integration, and then the effects of the hardening exponent n, the Poisson ratio ν, and the crack growth speed M on the asymptotic fields are studied. The approximate behaviors of the proposed solutions are discussed in the limit of ν → 0.5 or n → ∞.     

Two-parameter failure criterion for elastoplastic bodies with mode I cracks

The generalized Dugdale crack model is used to formulate two-parameter failure criteria for the cases of quasibrittle state and developed plastic zones at a mode I crack tip. The failure criteria relate the fracture strength characteristics and the stress mode at the crack tip through the plastic constraint factor. The critical state of bodies with cracks under uni-and biaxial loading is analyzed in the cases of plane stress and plane strain using the Tresca and von Mises yield criteria. A small-scale yield criterion, which is an analytic relation between the critical stress intensity factor and T-stresses, is established __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 7, pp. 47–57, July 2007.     

A crack in a piezoelectric layer bonded to a dissimilar elastic layer under transient load

Summary  A piezoelectric layer bonded to the surface of an elastic structure is considered. The piezoelectric and the elastic layers are infinite along the x-axis and have finite thickness in the y-direction. The polarization direction of the piezoelectric material is along the y-axis. By means of the method of singular integral equations, the solution in a Laplace transform plane is demonstrated. Laplace inversion yields the results in the time domain. Numerical values of the crack tip fields under in-plane transient electromechanical loading are obtained. The influence of layers thickness on stress and electric displacement intensity factors is investigated. Received 16 March 2000; accepted for publication 16 August 2000     

Crack-tip fields in steady crack-growth with linear strain-hardening

Singular stress and strain fields are found at the tip of a crack growing steadily and quasi-statically into an elastic-plastic strain-hardening material. The material is characterized byJ₂ flow theory together with a bilinear effective stress-strain curve. The cases of anti-plane shear, plane stress and plane strain are each considered. Numerical results are given for the order of the singularity, details of the stress and strain-rate fields, and the near-tip regions of plastic loading and elastic unloading.     

Stress and finite-strain analysis of a circular ring under diametral compression

Using a moiré, large-strain analysis method, a complete solution is shown in this paper of the fields of strain and stress for a circular ring subjected to diametral compression between two flat platens. The isotheticsu andv, obtained using 1000-lines-per-inch gratings, were differentiated photographically by the shifting technique (moiré-of-moiré) to determine ?u/?x, ?v/?y, ?u/?y and ?v/?x. Using the exact finite strain-displacement relationship, the Eulerian strains ? _x ^E , ? _y ^E and γ _xy ^E were computed. From these, the principal Eulerian strains were obtained. These results were verified with the isochromatics obtained from a large-deformation photoelasticity analysis. The ring was made of a polyurethane rubber which exhibits a linear relationship between natural strain and a newly introduced concept of “natural stress”. The Eulerian strains were converted to natural strains, and from these natural stresses were computed using the newly developed concept. Results are presented graphically for the whole field of the ring.     

State vector approach to analysis of multilayered magneto-electro-elastic plates

The state vector equations for three dimensional, orthotropic and linearly magneto-electro-elastic media are derived from the governing equations by eliminating σ_x, σ_y·τ_xy, B_x, B_y, D_x and D_y. An efficient method is presented for analysis of multilayered magneto-electro-elastic plates. The methodology is based on the mixed formulation, in which basic unknowns are formed by collecting not only displacements, electrical potential and magnetic potential but also some of stresses, electrical displacements, and magnetic induction. As special case, simply supported and multilayered rectangular plate is analyzed under the surface loading. Numerical results are presented graphically. The procedure of numerical calculation shows that the formulation presented here is simple and direct.     

On the relationship between critical tensile stress and fracture toughness in mild steel

An analysis is presented which relates the critical value of tensile stress (σ_f) for unstable cleavage fracture to the fracture toughness (K_Ic) for a high-nitrogen mild steel under plane strain conditions. The correlation is based on (i) the model for cleavage cracking developed by E. Smith and (ii) accurate plastic-elastic solutions for the stress distributions ahead of a sharp crack derived by J. R. Rice and co-workers. Unstable fracture is found to be consistent with the attainment of a stress intensification close to the tip such that the maximum principal stress σ_yy exceeds σ_f over a characteristic distance, determined as twice the grain size. The model is seen to predict the experimentally determined variation of K_Ic with temperature over the range -150 to -75°C from a knowledge of the yield stress and hardening properties. It is further shown that the onset of fibrous fracture ahead of the tip can be deduced from the position of the maximum achievable stress intensification. The relationship between the model for fracture ahead of a sharp crack, and that ahead of a rounded notch, is discussed in detail.     

Magneto elastic stress subjected to uniform magnetic field over a thin infinite plate containing an elliptical hole with an edge crack

Two dimensional solutions of the magnetic field and magneto elastic stress are presented for a magnetic material of a thin infinite plate containing an elliptical hole with an edge crack subjected to uniform magnetic field. Using a rational mapping function, each solution is obtained as a closed form. The linear constitutive equation is used for these analyses. According to the electro-magneto theory, only Maxwell stress is caused as a body force in a plate. In the present paper, it raises a plane stress state for a thin plate, the deformation of the plate thickness and the shear deflection. Therefore the magneto elastic stress is analyzed using Maxwell stress. No further assumption of the plane stress state that the plate is thin is made for the stress analysis, though Maxwell stress components are expressed by nonlinear terms. The rigorous boundary condition expressed by Maxwell stress components is completely satisfied without any linear assumptions on the boundary. First, magnetic field and stress analyses for soft ferromagnetic material are carried out and then those analyses for paramagnetic and diamagnetic materials are carried out. It is stated that those plane stress components are expressed by the same expressions for those materials and the difference is only the magnitude of the permeability, though the magnetic fields H_x, H_y are different each other in the plates. If the analysis of magnetic field of paramagnetic material is easier than that of soft ferromagnetic material, the stress analysis may be carried out using the magnetic field for paramagnetic material to analyze the stress field, and the results may be applied for a soft ferromagnetic material. It is stated that the stress state for the magnetic field H_x, H_y is the same as the pure shear stress state. Solutions of the magneto elastic stress are nonlinear for the direction of uniform magnetic field. Stresses in the direction of the plate thickness and shear deflection are caused and the solutions are also obtained. Figures of the magnetic field and stress distribution are shown. Stress intensity factors are also derived and investigated for the crack length.