Analysis to the stress and strain structures of mode I 3-D crack in elastic-plastic state

In this paper, the stress and strain structures of Mode I 3-D crack in power hardening material are studied by analyzing the fundamental equations of elastic-plastic mechanics. It is shown that three regions, Z1,Z2 and Z3 can be divided in the thickness direction according to the stress characteristic. In region Z1, the stress components in the plane Perpendicular to z axis (thickness direction) can be solved first using the fundamental equations of plane strain state; in region Z3, they can be solved first by the equations of plane stress state. The region Z2 is defined as a transition layer. It is shown that the transition layer is the characteristic of Mode I 3-D crack in elastic-plastic state, and it is significant to the research on 3-D fracture. The crack tip opening displacement CTOD is chosen to describe the amplitude coefficient of the local stress field, and the distribution of CTOD in 3-D state is investigated.The project supported by National Natural Science Foundation of China.     

Mode I and mode II crack tip asymptotic fields with strain gradient effects

The strain gradient effect becomes significant when the size of fracture process zone around a crack tip is comparable to the intrinsic material lengthl, typically of the order of microns. Using the new strain gradient deformation theory given by Chen and Wang, the asymptotic fields near a crack tip in an elastic-plastic material with strain gradient effects are investigated. It is established that the dominant strain field is irrotational. For mode I plane stress crack tip asymptotic field, the stress asymptotic field and the couple stress asymptotic field can not exist simultaneously. In the stress dominated asymptotic field, the angular distributions of stresses are consistent with the classical plane stress HRR field; In the couple stress dominated asymptotic field, the angular distributions of couple stresses are consistent with that obtained by Huang et al. For mode II plane stress and plane strain crack tip asymptotic fields, only the stress-dominated asymptotic fields exist. The couple stress asymptotic field is less singular than the stress asymptotic fields. The stress asymptotic fields are the same as mode II plane stress and plane strain HRR fields, respectively. The increase in stresses is not observed in strain gradient plasticity for mode I and mode II, because the present theory is based only on the rotational gradient of deformation and the crack tip asymptotic fields are irrotational and dominated by the stretching gradient. The project supported by the National Natural Science Foundation of China (19704100), National Natural Science Foundation of Chinese Academy of Sciences (KJ951-1-20), CAS K.C. Wong Post-doctoral Research Award Fund and Post-doctoral Science Fund of China     

Anisotropic plastic stress field at a mixed-mode crack tip under plane and anti-plane strain

On condition that any perfectly plastic stress component at a crack tip is nothingbut the function ofθ.by making use of equilibrium equations,anisotropic plastic stress-strain-rate relations,compatibility equations and Hill anisotropic plastic yieldcondition,in the present paper,we derive the generally analytical expressions of theanisotropic plastic stress field at a mixed-mode crack tip under plane and anti-planestrain.Applying these generally analytical expressions to the mixed-mode cracks,wecan obtain the analytical expressions of anisotropic plastic stress fields at the tips ofmixed-modeⅠ-Ⅲ,Ⅱ-ⅢandⅠ-Ⅱ-Ⅲcracks.     

Determination of elastic-plastic stresses and strains from measured surface strain data

A rigorous approach founded in the fundamental principles of plasticity is used to develop an accurate numerical algorithm for the determination of stresses and elastic and plastic strains from total strain data measured on a structure surface. The approach used to develop the algorithm and its relationship to both the flow theory of plasticity and recent advances in tangent stiffness-based numerical solution procedures for elastic-plastic boundary value problems are presented. Verification of the method for plane stress problems is demonstrated. A discussion of how the method can be used with measured surface displacement data is proved.     

航空结构铝合金Ⅰ/Ⅱ复合型断裂的三维弹塑性数值分析

在LC4CS铝合金三维Ⅰ/Ⅱ复合型断裂实验的基础上,对Ⅰ/Ⅱ复合型加载模型进行了细致的三维弹塑性有限元计算,系统分析了不同复合度φ和厚度B对裂纹端部应力场和起裂角的影响.结果表明,对于不同φ和B的试样,周向应力极大值和三维离面约束因子极小值出现在相同的角度;根据最大周向应力准则,在Ⅰ型载荷占主导(βeq＞65.45°)...     

A weighted residual method for elastic-plastic analysis near a crack tip and the calculation of the plastic stress intensity factors

In this paper,a weighted residual method for the elastic-plastic analysis near a crack tip is systematically given by taking the model of power-law hardening under plane strain condition as a sample.The elastic-plastic solutions of the crack tip field and an approach based on the superposition of the nonlinear finite element method on the complete solution in the whole crack body field,to calculate the plastic stress intensity factors,are also developed.Therefore,a complete analysis based on the calculation both for the crack tip field and for the whole crack body field is provided.     

The brittle fracture criterion based on the maximum tensile stress on the surface of blunt crack tip

The brittle fracture criterion is developed for a blunted crack. The curvature radius of the blunt crack tip is suggested as a characteristic length for brittle materials, and then the fracture toughness of the brittle materials can be determined from the cohesion strength and the characteristic length of the materials.     

Modification of the fracture criterion for V-shaped notches (plane problem). Relationship between toughness and strength and structural parameters

A fracture criterion of the type of the Neuber-Novozhilov criterion is proposed to describe the fracture in the vicinity of the tip of a V-shaped notch under tensile and shear loading. In the proposed criterion, the limits of averaging of the stresses along the notch axis depend on the presence, location, and size of the initial defects in the material. The crystal lattice parameter of the initial material is chosen for the characteristic linear size. For a V-shaped notch subjected to tension and shear, simple equations are obtained that relate the stress intensity factors for the modified singularity coefficients, the singularity coefficients themselves, and the theoretical tensile and shear strengths of a single crystal of the material taking into account the damage to the material in the vicinity of the notch tip. The equations obtained allow a passage to the limit from a notch to a crack. It is shown that the classical critical stress intensity factor used in the strength analysis of cracked solids is not a material constant.Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 1, pp. 106–115, January–February, 2005.     

OnJ-integrals in the plane fracture of composite materials

By using a complex function method in this paper, the complex form of J-integral of mixed mode crack tip for unidirectional plate of linear-elastic orthotropic composites is obtained first by substituting crack tip stresses and displacements into general formula of J-integral. And then, the path-independence of this J-integral is proved. Finally, the computing formula of this J-integral is derived. As special examples, the complex forms, path-independence and computing formulae of J-integrals of mode Ⅰ and mode Ⅱ crack tips for unidirectional plate of linear-elastic orthotropic composites are given.     

Asymptotics of stresses and strain rates near the tip of a transverse shear crack in a material whose behavior is described by a fractional-linear law

An approximate solution of the problem of determining the fields of stresses and strain rates due to creep near the tip of a transverse shear crack in a material whose behavior is described by a fractional-linear law of the theory of steady-state creep is given. It is shown that the strain rates have a singularity of the type ∼ r^−α near the crack tip; the order of singularity α changes discretely, depending on the polar angle, and takes the values 1, 2/3, and 1/2. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 1, pp. 165–176, January–February, 2009.     

Stress field of a damaged material near the crack tip

In this paper, a deformation theory of plasticity for damaged materials is proposed. An asymptotic expression forH near a crack tip is obtained. Finally, the stress and strain fields near the crack tip are presented.     

Modelling inelastic behaviour of orthotropic metals in a unique alignment of deviatoric plane within the stress space

A finite strain constitutive model to predict the deformation behaviour of orthotropic metals is developed in this paper. The important features of this constitutive model are the multiplicative decomposition of the deformation gradient and a new Mandel stress tensor combined with the new stress tensor decomposition generalized into deviatoric and spherical parts. The elastic free energy function and the yield function are defined within an invariant theory by means of the structural tensors. The Hill’s yield criterion is adopted to characterize plastic orthotropy, and the thermally micromechanical-based model, Mechanical Threshold Model (MTS) is used as a referential curve to control the yield surface expansion using an isotropic plastic hardening assumption. The model complexity is further extended by coupling the formulation with the shock equation of state (EOS). The proposed formulation is integrated in the isoclinic configuration and allows for a unique treatment for elastic and plastic anisotropy. The effects of elastic anisotropy are taken into account through the stress tensor decomposition and plastic anisotropy through yield surface defined in the generalized deviatoric plane perpendicular to the generalized pressure. The proposed formulation of this work is implemented into the Lawrence Livermore National Laboratory-DYNA3D code by the modification of several subroutines in the code. The capability of the new constitutive model to capture strain rate and temperature sensitivity is then validated. The final part of this process is a comparison of the results generated by the proposed constitutive model against the available experimental data from both the Plate Impact test and Taylor Cylinder Impact test. A good agreement between experimental and simulation is obtained in each test.     

Application of a stress triaxiality dependent fracture criterion in the finite element analysis of unnotched Charpy specimens

Nonlinear dynamic finite element analysis (FEA) is conducted to simulate the fracture of unnotched Charpy specimens of steel under pendulum impact loading by a dedicated, oversized and nonstandard Bulk Fracture Charpy Machine (BFCM). The impact energy needed to fracture an unnotched Charpy specimen in a BFCM test can be two orders of magnitude higher than the typical impact energy of a Charpy V-notch specimen. To predict material failure, a phenomenological, stress triaxiality dependent fracture initiation criterion and a fracture evolution law in the form of strain softening are incorporated in the constitutive relations. The BFCM impact energy results obtained from the FEA simulations compare favorably with the corresponding experimental data. In particular, the FEA predicts accurately the correlations of the BFCM impact energy with such factors as specimen geometry, impactor tup width and material type. The analyses show that a specimen’s progressive deterioration through the thickness dimension displays a range of shear to ductile fracture modes, demonstrating the necessity of applying a stress state dependent fracture initiation criterion. Modeling the strain softening behavior helps to capture the residual load carrying capability of a ductile metal or alloy beyond the onset of damage. The total impact energy can be significantly under predicted if a softening branch is not included in the stress-strain curve. This research supports a study of the puncture failure of railroad tank cars under dynamic impact loading. Applications of the presented fracture model in failure analyses of other structures are further discussed.     

The singularity analysis of the stress and strain fields for Mode I fracture in elastic-plastic state

The stress and strain singularities of power hardening material for Mode I fracrure are analysed according to the fundamental equations of elastic-plastic mechanics. It is found that the singularities of all stress and strain components do not change in the thick direction, and neither the six stress components nor the six strain components have the same singularity.     

Estimating the error of the beam approximation in the plane stability problem for a rectangular plate with a central crack

The plane stability problem for a rectangular, linearly elastic, isotropic plate with a central crack is solved. The dependence of the critical load on the crack length is studied using exact (the three-dimensional linearized theory of stability of elastic bodies) and approximate (beam approximation) approaches. The results produced by the beam approach are evaluated.Translated from Prikladnaya Mekhanika, Vol. 40, No. 11, pp. 117–126, November 2004.This revised version was published online in April 2005 with a corrected cover date.     

Refinement of the plastic-zone boundary in the vicinity of a crack tip for the quasiviscous and viscous types of fracture

A refined solution of the elastoplastic problem of an insulated mode I crack in a thin plate of reasonably large dimensions is obtained. Estimates of the plastic zone in the vicinity of the crack tip are given for quasiviscous and viscous types of fracture.Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 1, pp. 126–132, January–February, 2005     

On the types and number of plane waves in hypoelastic materials

General principles are formulated for modeling the elastic deformation of materials and analyzing plane waves in nonlinearly elastic materials such as hyperelastic, hypoelastic, and those governed by the general law of elasticity. The results of studying the propagation of plane waves in hypoelastic materials are further outlined. The influence of initial stresses and initial velocities on the types and number of plane waves is studied. Wave effects characteristic of hypoelastic materials are predicted theoretically. One of such effects is blocking of certain types of plane waves by initial stresses __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 11, pp. 96–107, November 2005.     

THE ANALYSES OF THE THREE-DIMENSIONAL STRESS STRUCTURE NEAR THE CRACK TIP OF MODE I CT SPECIMENS IN ELASTICPLASTIC STATE(Ⅰ)—THE ANALYSES OF CONSTRAINT PARAMETERS AND FRACTURE PARAMETERS

In the present paper, three dimensional analyses of some general constraint parameters and fracture parameters near the crack tip. of Mode I CT specimens in two different thicknesses are carried out by employing ADINA program. The results reveal that the constraints along the thickness direction are obviously separated into two parts: the keeping similar high constraint field (Z₁) and rapid reducing constraints one (Z₂). The two fields are experimentally confiremed to correspond to the smooth region and the shear lip on the fracture face respectively. So the three dimensional stress structure of Mode I specimens can be derived through discussing the two fields respectively. The distribution of the Crack Tip Opening Displacement (CTOD) along the thickness direction and the three dimensional distribution of the void growth ratio (V_g) near the crack tip are also obtained. The two fracture parameters are in similar trends along the thickness direction, and both of them can reflect the effect of thickness and that of the loading level to a certain degree.