Finite deformation analysis of crack tip fields in plastically compressible hardening–softening–hardening solids

Crack tip fields are calculated under plane strain small scale yielding conditions. The material is characterized by a finite strain elastic-viscoplastic constitutive relation with various hardening-softening-hardening hardness functions. Both plastically compressible and plastically incompressible solids are considered. Displacements corresponding to the isotropic linear elastic mode I crack field are prescribed on a remote boundary. The initial crack is taken to be a semi-circular notch and symmetry about the crack plane is imposed. Plastic compressibility is found to give an increased crack opening displacement for a given value of the applied loading. The plastic zone size and shape are found to depend on the plastic compressibility, but not much on whether material softening occurs near the crack tip. On the other hand, the near crack tip stress and deformation fields depend sensitively on whether or not material softening occurs. The combination of plastic compressibility and softening (or softening-hardening) has a particularly strong effect on the near crack tip stress and deformation fields.

EXACT SOLUTIONS OF NEAR CRACK LINE FIELDS FOR MODE I CRACK UNDER PLANE STRESS CONDITION IN AN ELASTIC－PERFECTLY PLASTIC SOLID

The near crack line analysis method has been used in the present paper,The classical small scale yielding conditions have been completely abandoned in the analyses and one inappropriate matching condition used to be used at the elasticplastic boundary has been corrected.The reasonable solution of the plastic stresses near the crack line region has been established.By matching the plastic stresses with the exact elastic stresses at the elastic-plastic boundary,the plastic stresses the length of the plastic zone and the unit normal vector of the elastic-plastic boundary near the crock line region have been obtained for a mode I crack under uniaxial tension,as well as a mode I crack under biaxial tension,which shows that for both conditions the plastic stress componentsσy, and σsy.he length of the plastic zone and the unit normal vector of the elastic-plastic boundary are quite the same while the plastic stress σs is different.

MODE Ⅰ AND MODE Ⅱ CRACK TIP ASYMPTOTIC FIELDS WITH STRAIN GRADIENT EFFECTS

The strain gradient effect becomes significant when the size of frac- ture process zone around a crack tip is comparable to the intrinsic material length l, 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 dom- inant strain field is irrotational. For mode Ⅰ plane stress crack tip asymptotic field, the stress asymptotic field and the couple stress asymptotic field can not exist si- multaneously. 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 Ⅱ 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 Ⅱ plane stress and plane strain HRR fields, respectively. The increase in stresses is not observed in strain gradient plasticity for mode Ⅰ and mode Ⅱ, because the present theory is based only on the rotational gradi- ent of deformation and the crack tip asymptotic fields are irrotational and dominated by the stretching gradient.

THE CRACK TIP FIELDS OF ORTHOTROPIC COMPOSITE PLATE UNDER BENDING

In this paper, mechanics analysis for crack tip fields of linear elastic orthotropic composite plate under symmetric bending load was done. By using a complex variable method, the equations for bending moment, twisting moment, stress, strain and displacement fields near crack tip are derived.

Higher-order crack tip fields for functionally graded material plate with transverse shear deformation

The crack tip fields are investigated for a cracked functionally graded material （FGM） plate by Reissner＇s linear plate theory with the consideration of the transverse shear deformation generated by bending. The elastic modulus and Poisson＇s ratio of the functionally graded plates are assumed to vary continuously through the coordinate y, according to a linear law and a constant, respectively. The governing equations, i.e., the 6th-order partial differential equations with variable coefficients, are derived in the polar coordinate system based on Reissner＇s plate theory. Furthermore, the generalized displacements are treated in a separation-of-variable form, and the higher-order crack tip fields of the cracked FGM plate are obtained by the eigen-expansion method. It is found that the analytic solutions degenerate to the corresponding fields of the isotropic homogeneous plate with Reissner＇s effect when the in-homogeneity parameter approaches zero.

The crack tip fields of Mode Ⅱ stationary growth crack on bimaterial interface

A mechanical model is established for mode II interfacial crack static growing along an elastic-elastic power law creeping bimaterial interface. For frictional contact of boundary conditions on crack faces, asymptotic solutions of the stresses and strains of near tip-crack are got. It was shown that in stable creep growing phase, elastic deformation and viscous deformation are equally dominant in the near-tip field, the stress and strain have the same singularity and there is not the oscillatory singularity the field. Through numerical calculation , it is shown that the frictional coefficient η notably influence the crack-tip field.

ANISOTROPIC PLASTIC STRESS FIELDS AT A SLOWLY PROPAGATING CRACK TIP

Under the condition that any perfeetly plastic stress components at a crack tip are nothing but the Junctions of 0 only, making use of equilibriumequations,Hill ani.sutropic yield condition and unloading stress-strain relations, in this paper, we derive the general analytical expressions of anisotropic plastiestress Jields at the slowly steadyhe slowly steady propagatin tips of plane and anti-phane strain,Applying these general analytical expressions to the concrete cracks the attchvtical expressions of anisotropie plastic stress fields at the slowly steady propagating tips of Motle I and Motle III cracks are obtained. For the isolropic plastic material, the anisotropic plastic stress fields at a slowly propagating crack tip become the perfeeby plastic mress fields

CRACK TIP FIELD AND J-INTEGRAL WITH STRAIN GRADIENT EFFECT

The mode I plane strain crack tip field with strain gradient effects is presented in this paper based on a simplified strain gradient theory within the framework proposed by Acharya and Bassani.The theory retains the essential structure of the incremental version of the conventional J_2 deformation theory.No higher-order stress is introduced and no extra boundary value conditions beyond the conventional ones are required.The strain gradient effects are considered in the constitutive relation only through the instantaneous tangent modulus.The strain gradient measures are included into the tangent modulus as internal parameters.Therefore the boundary value problem is the same as that in the conventional theory.Two typical crack problems are studied:(a)the crack tip field under the small scale yielding condition induced by a linear elastic mode-I K-field and(b)the complete field for a compact tension specimen.The calculated results clearly show that the stress level near the crack tip with strain gradient effects is considerable higher than that in the classical theory.The singularity of the strain field near the crack tip is nearly equal to the square-root singularity and the singularity of the stress field is slightly greater than it.Consequently,the J-integral is no longer path independent and increases monotonically as the radius of the calculated circular contour decreases.

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 lip 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 analvsis based on the calculation both for the crack tip field and for the whole crack body field is provided.

Field structure at mode Ⅲ dynamically propagating crack tip in elastic-viscoplastic materials

An elastic-viscoplastic mechanics model is used to investigate asymptotically the mode Ⅲ dynamically propagating crack tip field in elastic-viscoplastic materials. The stress and strain fields at the crack tip possess the same power-law singularity under a linear-hardening condition. The singularity exponent is uniquely determined by the viscosity coefficient of the material. Numerical results indicate that the motion parameter of the crack propagating speed has little effect on the zone structure at the crack tip. The hardening coefficient dominates the structure of the crack-tip field. However, the secondary plastic zone has little influence on the field. The viscosity of the material dominates the strength of stress and strain fields at the crack tip while it does have certain influence on the crack-tip field structure. The dynamic crack-tip field degenerates into the relevant quasi-static solution when the crack moving speed is zero. The corresponding perfectly-plastic solution is recovered from the linear-hardening solution when the hardening coefficient becomes zero.

The stress field and rupture feature in compound model containing typical cracks system and hard-inclusion

Aiming at the research on mechanical mechanism of hard-inclusion earthquake preparation model, in this paper,experimental and contrast research on stress field and rupture feature of hard-inclusion model has been made respectively, which contained en echelon and composite cracks systems in models, and was loaded under uniaxial compressive stress. The result shows that reverse en echelon and T-shape cracks systems in hard-inclusion are the favorable geological structures to trigger earthquakes.

Characteristics of the crack tip field in high-speed railway tunnel linings under train-induced aerodynamic shockwaves

High-speed railway tunnels in various countries have continuously reported accidents of vault falling concrete blocks.Once the concrete block falling occurs,serious consequences follow,and traffic safety may be endangered.The aerodynamic shockwave evolves from the initial compression wave may be an important inducement causing the tunnel lining cracks to grow and form falling concrete blocks.A joint calculation framework is established based on ANSYS Fluent,ABAQUS,and FRANC3D for calculating the crack tip field under the aerodynamic shockwave.The intensification effect of aerodynamic shockwaves in the crack is revealed,and the evolution characteristics of the crack tip field and the influence factors of stress intensity factor(SIF)are analyzed.Results show that(1)the aerodynamic shockwave intensifies after entering the crack,resulting in more significant pressure in the crack than the input pressure.The maximum pressure of the inclined and longitudinal cracks is higher than the corresponding values of the circumferential crack,respectively.(2)The maximum SIF of the circumferential,inclined,and longitudinal crack appears at 0.5,0.68,and 0.78 times the crack front length.The maximum SIF of the circumferential crack is higher than that of the inclined and longitudinal crack.The possibility of crack growth of the circumferential crack is the highest under aerodynamic shockwaves.(3)The influence of train speed on the SIF of the circumferential crack is more than 40%.When the train speed,crack depth,and crack length change,the change of pressure in the crack is the direct cause of the change of SIF.

Transient near tip fields in crack dynamics

Transient effects of stress-strain fields in the vicinity of a stationary crack tip under high rate loads are discussed.Exact analytical solutions to near tip stresses are compared to fields prescribed by leading terms(one or several) of Williams asymptotic expansion.Influence of load application mode,time(or,which is the same,distance from a crack tip) and Poisson's ratio on this discrepancy is extensively examined.Some effects connected with crack tip propagation speed are also discussed.Significant inconsistencies between real(or received in numerical solutions of state equations-e.g.finite element computations) crack tip fields and stress intensity factor(SIF) singular field observed by numerous researchers are explained.The scope of problems where SIF field can be used for correct prediction of dynamic stress-strain fields in the crack tip region is established.Possibility to correctly approximate fields that are not SIF dominated,accounting additional terms of Williams expansion,is studied.

Incompatible numerical manifold method for fracture problems

The incompatible numerical manifold method （INMM） is based on the finite cover approximation theory, which provides a unified framework for problems dealing with continuum and discontinuities. The incompatible numerical manifold method employs two cover systems as follows. The mathematical cover system provides the nodes for forming finite covers of the solution domain and the weighted functions, and the physical cover system describes geometry of the domain and the discontinuous surfaces therein. In INMM, the mathematical finite cover approximation theory is used to model cracks that lead to interior discontinuities in the process of displacement. Therefore, the discontinuity is treated mathematically instead of empirically by the existing methods. However, one cover of a node is divided into two irregular sub-covers when the INMM is used to model the discontinuity. As a result, the method sometimes causes numerical errors at the tip of a crack. To improve the precision of the INMM, the analytical solution is used at the tip of a crack, and thus the cover displacement functions are extended with higher precision and computational efficiency. Some numerical examples are given.

Stability of Quark Star Models

In this paper,we investigate the stability of quark stars with four different types of inner matter configurations;isotropic,charged isotropic,anisotropic and charged anisotropic by using the concept of cracking.For this purpose,we have applied local density perturbations technique to the hydrostatic equilibrium equation as well as on physical parameters involved in the model.We conclude that quark stars become potentially unstable when inner matter configuration is changed and electromagnetic field is applied.