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 constraintparameters and fracture parameters near the crack tip of Mode I CT specimens in twodifferent thicknesses are carried out by employing ADINA program.The results revealthat the constraints along the thickness direction are obviously separated into twoparts:the keeping similar high constraint field(Z_1)and rapid reducing constraintsone(Z_2).The two fields are experimentally confiremed to correspond to the smoothregion and the shear lip on the fracture face respectively.So the three dimensionalstress structure of Mode I specimens can be derived through discussing the two fieldsrespectively.The distribution of the Crack Tip Opening Displacement(CTOD)alongthe thickness direction and the three dimensional distribution of the void growth ratio(V_9 )near the crack tip are also obtained.The two fracture parameters are in similartrends along the thickness direction,and both of them can reflect the effect ofthickness and that of the loading l     

Buckling of an orthotropic graded coating with an embedded crack bonded to a homogeneous substrate

To simulate buckling of nonuniform coatings, we consider the problem of an embedded crack in a graded orthotropic coating bonded to a homogeneous substrate subjected to a compressive loading. The coating is graded in the thickness direction and the material gradient is orthogonal to the crack direction which is parallel with the free surface. The elastic properties of the material are assumed to vary continuously along the thickness direction. The principal directions of orthotropy are parallel and perpendicular to the crack orientation. The loading consists of a uniform compressive strain applied away from the crack region. The graded coating is modeled as a nonhomogeneous medium with an orthotropic stress–strain law. Using a nonlinear continuum theory and a suitable perturbation technique, the plane strain problem is reduced to an eigenvalue problem describing the onset of buckling. Using integral transforms, the resulting plane elasticity equations are converted analytically into singular integral equations which are solved numerically to yield the critical buckling strain. The Finite Element Method was additionally used to model the crack problem. The main objective of the paper is to study the influence of material nonhomogeneity on the buckling resistance of the graded layer for various crack positions, coating thicknesses and different orthotropic FGMs.     

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.     

THE ANALYSES OF THE THREE-DIMENSIONAL STRESS STRUCTURE NEAR THE CRACK TIP OF MODE I CT SPECIMENS IN ELASTICPLASTIC STATE(II)—THE ANALYSES OF THE STRESS STRUCTURE

Based on [1], the stress structures of the smooth region and shear lip of the specimens have been investigated in the paper. The characteristics of the stress structure in the smooth region have been found that the variable z can separated out; the stresses in the midsection can be obtained by the plane strain FEM results or HRR structure modified by the stress triaxiality. The effects of load level and thickness on the stress structure can be reflected by the distribution of CTOD along the thickness direction. The obtained expressions of the stresses are very simple and visualized. The analyses of the stress structure in the shear lip show that the stresses can be obtained by different methods of interpolation to a certain precise degree. A new degree parameter of the plane strain state has been put forward and studied. The parameter can reflect relatively well the variation of the kind and thickness of the specimen as well as the load level. The fracture parameter has also been investigated to be sure that it can be obtained by modified CTOD with the stress triaxiality.     

THE ANALYSES OF THE THREE-DIMENSIONAL STRESS STRUCTURE NEAR THE CRACK TIP OF MODE Ⅰ CT SPECIMENS IN ELASTICPLASTIC STATE(Ⅱ)—THE ANALYSES OF THE STRESS STRUCTURE

Based on[1],the stress structures of the smooth region and shear lip of thespecimens have been investigated in the paper.The characteristics of the stress structure inthe smooth region have been found that the variable z can separated out;the stressesin the midsection can be obtained by the plane strain FEM results or HRR structuremodified by the stress triaxiality.The effects of load level and thickness on the stressstructure can be reflected by the distribution of CTOD along the thickness direction.The obtained expressions of the stresses are very simple and visualized.The analysesof the stress structure in the shear lip show that the stresses can be obtained bydifferent methods of interpolation to a certain precise degree.A new degree parameter of the plane strain state has been put forward andstudied.The parameter can reflect relatively well the variation of the kind andthickness of the specimen as well as the load level.The fracture parameter has alsobeen investigated to be sure that it can be obta     

晶体取向对单晶体断裂特征影响的模拟分析

结合Ni基单晶合金制三种不同晶体取向的紧凑拉伸试样试验,本文利用考虑有限变形和品格转动效应的晶体滑移有限元程序对单晶体三维断裂特征进行了模拟计算分析,详细考察了裂纹尖端三维应力场特征和断裂特征,结果表明:晶体取向对裂纹尖端应力场有较大影响,但应力沿试样厚度方向明显分成两个部分,在试样心部,应力沿厚度方向变化不大,在试样外表面则明显变化。裂纹尖端张开位移(CTDD)沿厚度方向类似分成两个部分。垂直于滑移面的应力分量致单晶体的准解理断裂,即裂纹的起裂和扩展途径均与该应力分量有关。     

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.     

可压缩理想塑性体中逐步扩展裂纹顶端的弹塑性场

本文利用理想塑性固体平面应变问题的基本方程,分析了可压缩理想塑性体中逐步扩展裂纹顶端的弹塑性场,得到了关于应力的渐近场,分析了弹性卸载区的演变过程和修正的中心扇形区的发展过程,预示了出现二次塑性区的可能性,弹性可压缩性的影响明显表现在经典的中心扇形区必需加以修正,垂直于板面方向的应力偏量不再为零,而且随着新裂纹面的形成,裂纹前方的均匀应力场和紧连着的修正的中心扇形区的应力偏量将发生变化,这种变化是由于垂直于板面方向的应力偏量发生变化造成的。     

Mode III fracture problem of an arbitrarily oriented crack in an FGPM strip bonded to a homogeneous piezoelectric half plane

This paper studies the Mode III electric-elastic field of a cracked functionally graded piezoelectric strip bonded to a homogeneous piezoelectric half plane. The crack is oriented in arbitrary direction. The material properties of the strip vary along the strip thickness in exponential forms. By using the Fourier transform, the problem can be formulated to a system of singular integral equations and solved by applying the Gauss-Chebyshev integration formula. The effects come from the edge, crack orientations and the nonhomogeneous material parameter on intensity factors are discussed graphically.     

An Application of the Kane and Mindlin Theory to Crack Problems in Plates of Arbitrary Thickness

Classical plane solutions of the theory of elasticity, which are sometimes more than 100 years old, are still used today and provide a framework for the analysis of many practical problems. But, strictly speaking, these analytical solutions are only applicable to plates with vanishing thickness or infinite thickness, where the stress state could be classified as plane stress or plane strain, respectively. However, the through-the-thickness stresses that exist in a plate of given thickness have a significant impact in a number of practical applications; and these stresses are often inevitably ignored due to the lack of analytical tools. This paper presents new analytical results for crack tip opening displacement (CTOD) for the through-the-thickness crack in infinite plates with various thicknesses. These results are based on the solution for an edge dislocation in infinite plate of arbitrary thickness and an application of the distributed dislocation technique. The analytical predictions of the CTOD and the constraint factor are compared with the three-dimensional elasto-plastic finite element (FE) results. It is shown that both analytical and numerical results are in good agreement when the numerical calculations are not affected by the size of the FE mesh and by the boundaries of the FE model.     

Analyses of steady quasistatic crack growth in plane strain tension in elastic-plastic materials with non-isotropic hardening

Steady state crack propagation problems of elastic-plastic materials in Mode I, plane strain under small scale yielding conditions were investigated with the aid of the finite element method. The elastic-perfectly plastic solution shows that elastic unloading wedges subtended by the crack tip in the plastic wake region do exist and that the stress state around the crack tip is similar to the modified Prandtl fan solution. To demonstrate the effects of a vertex on the yield surface, the small strain version of a phenomenological J₂, corner theory of plasticity (Christoffersen, J. and Hutchinson, J. W. J. Mech. Phys. Solids,27, 465 C 1979) with a power law stress strain relation was used to govern the strain hardening of the material. The results are compared with the conventional J₂ incremental plasticity solution. To take account of Bauschinger like effects caused by the stress history near the crack tip, a simple kinematic hardening rule with a bilinear stress strain relation was also studied. The results are again compared with the smooth yield surface isotropic hardening solution for the same stress strain curve. There appears to be more potential for steady state crack growth in the conventional J₂ incremental plasticity material than in the other two plasticity laws considered here if a crack opening displacement fracture criterion is used. However, a fracture criterion dependent on both stress and strain could lead to a contrary prediction.     

A Zener–Stroh crack interacting with a coated inclusion with generalized Irwin plastic zone correction

The elastic–plastic fracture behavior of a Zener–Stroh crack interacting with a coated inclusion in composite materials has been investigated with crack tip plastic zone corrections. With the distributed dislocation method, the crack problem is formulated into a set of singular integral equations which are solved numerically. The plastic zone sizes at the both crack tips are determined by a generalized Irwin model where Von Mises stress yielding criterion is used. The stress intensity factor (SIF), the plastic zone size (PZS), the crack tip opening displacement (CTOD) and the effective stress intensity factor have been evaluated. In the numerical examples, the influence of the inclusion shear modulus, the coating-layer thickness and shear modulus, as well as the distance between the crack and inclusion, on the SIF, the PZS and the CTOD are discussed in detail. Numerical examples show that increasing the shear modulus or the thickness of the coating phase, the influence of the inclusion on the normalized SIF and the normalized PZS will be shielded.     

Influence of non-singular stress terms and specimen geometry on small-scale yielding at crack tips in elastic-plastic materials

The plane strain elastic-plastic state at a crack tip is determined for compact tension, bend, double edge-cracked and centre-cracked specimens using a finite element method with triangular constant-strain elements. The solutions are found to differ by 10 to 30 per cent at the ASTM-limit as regards fracture surface displacement, normal stress and plastic zone size. In order to bring the boundary layer solution for the crack problem into agreement with the solution for a specific specimen one has to modify this solution. The modification consists of an addition to the boundary tractions for the boundary layer problem of tractions corresponding to the non-singular, constant second term in a series expansion of the normal stress parallel to the crack plane.     

STRESS-STRAIN STATE OF ELASTOPLASTIC BODIES WITH CRACK

A physical cut model is used to describe the changes in the stress-strain state(SSS)in elastoplastic bodies weakened by cracks. The distance between the crack edges is considered to be finite in contrast to the mathematical cut. The interactive layer with a thickness limited by the possibility of using the hypothesis of continuity is distinguished on the physical cut extension.Distribution of stresses and strains over the layer thickness is constant and does not depend on the geometry of the boundary between the cut and the interactive layer. The relationship between stresses and strains is determined by the deformation plasticity theory. The problem of plane strain or plane stress state of an arbitrary finite body weakened by a physical cut is reduced to solving a system of two variational equations for displacement fields in the body parts adjacent to the interactive layer. The proposed approach eliminates the singularity in stress distribution in contrast to the mathematical cut model. Use of local strength criteria allows us to determine the time, point and direction of the fracture initiation. Possibilities of the proposed model are illustrated by solving the problems of determining the SSS of a rectangular body weakened by a physical cut under symmetric and antisymmetric loadings.     

An embedded crack in a functionally graded orthotropic coating bonded to a homogeneous substrate under a frictional Hertzian contact

In this paper, we consider the elasto-static problem of an embedded crack in a graded orthotropic coating bonded to a homogeneous substrate subject to statically applied normal and tangential surface loading. The crack direction is parallel to the free surface. The coating is graded in the thickness direction and is orthogonal to the crack direction. This coating is modelled as a non-homogeneous medium with an orthotropic stress–strain law. The equivalent crack surface stresses are first obtained and substituted in the plane elasticity equations. Using integral transforms, the governing equations are converted into singular integral equations which are solved numerically to yield the displacement field as well as the crack-tip stress intensity factors. This study presents a complete theoretical formulation for the problem in the static case. A numerical predictive capability for solving the singular integral equations and computing the crack-tip stress intensity factors is proposed. Since the loading is compressive, a previously developed crack-closure algorithm is applied to avoid interpenetration of the crack faces. The main objective of the paper is to investigate the effects of the material orthotropy and non-homogeneity of the graded coating on the crack-tip stress intensity factors, with and without using the crack-closure algorithm, for the purpose of gaining better understanding on the behavior and design of graded coatings.     

Stress and energy density fields near a blunted crack front in power hardening material

Using Jaumann and Dienes rates of Euler stress in elastic-plastic constitutive equations of finite deformation, plane strain finite element analysis for a compact tension specimen with a blunted crack front is made. The Euler stress, Kirchhoff stress and volume strain energy density near a blunted crack tip are computed. Constitutive relations with different deformation rates affect the the near crack tip solution in a region within an order of magnitude of the crack opening displacement. The results differed from the corresponding solution of deformation plasticity (or nonlinear elasticity) with increasing deformation. They are smaller in a local region of about 2 to 10 times of the crack opening distance.The volume energy density near the crack tip is computed, the stationary values of which determine the locations of extensive yielding and possible sites of crack initiation. It remained nearly constant with increasing deformation. Such a character tends to support the volume energy density criterion as a means for quantifying the ductile fracture behavior of metals.     

Bonded half planes containing an arbitrarily oriented crack

The plane elastostatic problem for two bonded half planes containing an arbitrarily oriented crack in the neighborhood of the interface is considered. Using Mellin Transforms the problem is formulated as a system of singular integral equations. The equations are solved for various crack orientations, material combinations, and external loads. The numerical results given in the paper include the stress intensity factors, the strain energy release rates, and the probable cleavage angles giving the direction of crack propagation.     

Logarithm strain singularity at tip of mode I growing crack in elasticand perfectly-plastic material

Reanalyzed in detail is the stress and strain distribution near the tip of a Mode I steadily growing crack in an elastic and perfectly-plastic material. The crack tip region is divided into five angular sectors, one of which is singular in character and represents a rapid transition zone that becomes a line of strain discontinuity in the limit as crack tip is approached. It is shown for an incompressible material (ν=0.5) under plane strain that the local strain in all the angular sectors possesses the same logarithm singularity, i.e., In r where r is the radial distance measured from the crack tip. This result also prevails for the compressible material ( v < 0.5) and resolves a long standing controversy concerning the strain singularity in the sector just ahead of the crack tip.     

Reduction of dimensions in random,elastic soil medium

The paper deals with an application of the plane strain analysis in a stochastic three-dimensional soil medium. In a framework of random elasticity theory, the geostatical state of stresses and the problem of a unit force acting in a statistically homogeneous half-space are considered. Only the modulus of elasticity is considered to be random and is modelled as a three-dimensional (3-D) homogeneous random field. As the result of imposed constrains due to the plane strain assumption the additional body and surface forces are induced. In order to determine them, additional equations must be introduced. The equations in a form of constrain relations are proposed in this paper. These equations are also valid for a case of uniformly distributed external loading.First, the two-dimensional (2-D) problem and its reduction to the uni-axial strain state, for the gravity forces and uniform, unlimited surface loading is considered. Then, it is generalised into a 2-D schematization of the 3-D state. Next, the problem of a unit force acting in a statistically homogeneous half-space is considered. For a 3-D state of stress and strain the resulting stresses are compared with those for a 2-D state. These stresses for the multidimensional state of strain and stress are presented as a sum of two components. The first one reflects plane strain state stresses and is given in a form of a 3-D random field. This term allows for incorporating a spatial, 3-D soil variability into a two-dimensional analysis. The second component can be treated as a correction term and it represents the longitudinal influence of a 3-D analysis.Some numerical results are presented in this paper. The proposed method can be regarded as a framework for further research aiming at application to a variety of geotechnical problems, for which the plane strain state is assumed.