Composite materials dynamic fracture studies by generalized Shmuely difference algorithm

The generalized Shmuely Difference Algorithm (GSDA) is presented here to analyze the dynamic fracture performance of orthogonal-anisotropic composite materials, such as glass fibre reinforced phenolplast. The difference recurrence Formulae and boundary condition difference extrapolation formulae are derived and programmed. The dynamic stress intensity factors (DSIF) of the isotropic and anisotropic centrally cracked plates are computed respectively using GSDA and compared with that published previously. GSDA is proved effective and reliable. Copyright (C) 1996 Elsevier Science Ltd.

Dynamic damage and failure in viscoplastic materials

In this paper, a dynamic damage model in ductile solids under the application of a dynamic mean tensile stress is developed. The proposed model considers void nucleation and growth as parts of the damage process under intense dynamic loading (strain rates epsilon greater than or equal to 10(3) s(-1)). The evolution equation of the ductile void has the closed form, in which work-hardening behavior, rate-dependent contribution and inertial effects are taken into account. Meanwhile, a plate impact test is performed for simulating the dynamic fracture process in LY12 aluminum alloy. The damage model is incorporated in a hydrodynamic computer code, to simulate the first few stress reverberations in the target as it spalls and postimpact porosity in the specimen. Fair agreement between computed and experimental results is obtained. Numerical analysis shows that the influence of inertial resistance on the initial void growth in the case of high loading rate can not be neglected. It is also indicated that the dynamic growth of voids is highly sensitive to the strain rates.

Unloading characteristics of anti-plane shear crack in softening materials

The local characteristics of the anti-plane shear stress and strain field are determined for a material where the stress increases linearly with strain up to a limit and then softens nonlinearly. Two unloading models are considered such that the unloading path always returns to the origin while the other assumes the unloading modulus to be that of the initial shear modulus. As the applied shear increases, an unloading zone is found to prevail between a zone in which the material softens and another zone in which the material is linear-elastic even though the crack does not propagate. The divisions of these zones are displayed graphically.

Influence of inertial and thermal effects on the dynamic growth of voids in porous ductile materials

The influence of inertial, thermal and rate - sensitive effects on the void growth at high strain rate in a thermal - viscoplastic solid is investigated by means of a theoretical model presented in the present paper. Numerical analysis of the model suggests that inertial, thermal and rate - sensitive effects are three major factors which greatly influence the behavior of void growth in the high strain rate case. Comparison of the mathematical model proposed in the present work and Johnson's model shows that if the temperature - dependence is considered, material viscosity eta can take the experimentally measured values.

Dynamic failure in ductile porous materials

In this paper, a mathematical model of dynamic fracture in porous ductile materials under intense dynamic general loading is developed. The mathematical model includes the influence of inertial effects and material rate sensitivity, as well as the contribution of surface energy of a void and material work-hardening. In addition, the condition of the void compaction is considered as well. The threshold stresses for the void growth and compaction are obtained. A simple criterion for ductile fracture which is associated with material distention and plastic deformation is adopted. As an application of the theoretical model, the processes of two-dimensional spallation in LY12 aluminum alloy are successfully simulated by means of two-dimensional finite-difference Lagrangian code.

Higher-order analysis of crack-tip fields in power-law hardening materials

In this paper, we present an exact higher-order asymptotic analysis on the near-crack-tip fields in elastic-plastic materials under plane strain, Mode I. A four- or five-term asymptotic series of the solutions is derived. It is found that when 1.6 < n less-than-or-equal-to 2.8 (here, n is the hardening exponent), the elastic effect enters the third-order stress field; but when 2.8< n less-than-or-equal-to 3.7 this effect turns to enter the fourth-order field, with the fifth-order field independent. Moreover, if n>3.7, the elasticity only affects the fields whose order is higher than 4. In this case, the fourth-order field remains independent. Our investigation also shows that as long as n is larger than 1.6, the third-order field is always not independent, whose amplitude coefficient K3 depends either on K1 or on both K1 and K2 (K1 and K2 arc the amplitude coefficients of the first- and second-order fields, respectively). Firmly, good agreement is found between our results and O'Dowd and Shih's numerical ones[8] by comparison.

Asymptotic solution of mode-iii crack in damaged softening materials

In this paper, a constitutive model of elasticity coupled with damage suggested by Lemaitre et al, [1] is used. The macroscopic stress-strain response of the model includes two stages: strain hardening and strain softening. The basic equation is derived for the anti-plane shear problem. Several lowest order asymptotic solutions are obtained, and assembled for the crack-tip fields.

Higher-order analysis of crack tip fields in elastic power-law hardening materials

A HIGHER-ORDER asymptotic analysis of a stationary crack in an elastic power-law hardening material has been carried out for plane strain, Mode 1. The extent to which elasticity affects the near-tip fields is determined by the strain hardening exponent n. Five terms in the asymptotic series for the stresses have been derived for n = 3. However, only three amplitudes can be independently prescribed. These are K1, K2 and K5 corresponding to amplitudes of the first-, second- and fifth-order terms. Four terms in the asymptotic series have been obtained for n = 5, 7 and 10; in these cases, the independent amplitudes are K1, K2 and K4. It is found that appropriate choices of K2 and K4 can reproduce near-tip fields representative of a broad range of crack tip constraints in moderate and low hardening materials. Indeed, fields characterized by distinctly different stress triaxiality levels (established by finite element analysis) have been matched by the asymptotic series. The zone of dominance of the asymptotic series extends over distances of about 10 crack openings ahead of the crack tip encompassing length scales that are microstructurally significant. Furthermore, the higher-order terms collectively describe a spatially uniform hydrostatic stress field (of adjustable magnitude) ahead of the crack. Our results lend support to a suggestion that J and a measure of near-tip stress triaxiality can describe the full range of near-tip states.

The strength distribution of brittle materials with a high-concentration of cavities

The maximum stress concentration factor in brittle materials with a high concentration of cavities is obtained. The interaction between the nearest cavities, in addition to the far field interactions, is taken into account to evaluate the strength distribution based on the statistical analysis of the nearest distance distribution. Through this investigation, it is found that the interaction between the nearest neighbors is much more important than the far field interactions, and one has to consider it in calculating the strength of brittle materials even if the volume fraction of cavities it contains is small. The other important conclusion is that the maximum stress concentration factor has a wide scattered distribution.

Crack tip behavior and crack-propagation in ductile materials

Dilatational plastic equations, which can include the effects of ductile damage, are derived based on the equivalency in expressions for dissipated plastic work. Void damage developed internally at the large-strain stage is represented by an effective continuum being strain-softened and plastically dilated. Accumulation of this local damage leads to progressive failure in materials. With regard to this microstructural background, the constitutive parameters included for characterizing material behaviour have the sense of internal variables. They are not able to be determined explicitly by macroscopic testing but rather through computer simulation of experimental curves and data. Application of this constitutive model to mode-I cracking examples demonstrates that a huge strain concentration accompanied by a substantial drop of stress does occur near the crack tip. Eventually, crack propagation is simulated by using finite elements in computations. Two numerical examples show good accordance with experimental data. The whole procedure of study serves as a justification of the constitutive formulation proposed in the text.

On dilatational plastic constitutive equation of ductile materials and plastic loading paths at bifurcation

The dilatational plastic constitutive equation presented in this paper is proved to be in a form of generality. Based on this equation, the constitutive behaviour of materials at the moment of bifurcation is demonstrated to follow a loading path with the response as "soft" as possible.

Crack deflection at an interface between dissimilar elastic-materials

A crack intersecting an interface between two dissimilar materials may advance by either penetrating through the interface or deflecting into the interface. The competition between deflection and penetration can be assessed by comparison of two ratios: (i) the ratio of the energy release rates for interface cracking and crack penetration; and (ii) the ratio of interface to material fracture energies. Residual stresses caused by thermal expansion misfit can influence the energy release rates of both the deflected and penetrating crack. This paper analyses the role of residual stresses. The results reveal that expansion misfit can be profoundly important in systems with planar interfaces (such as layered materials, thin film structures, etc.), but generally can be expected to be of little significance in fiber composites. This paper corrects an earlier result for the ratio of the energy release rate for the doubly deflected crack to that for the penetrating crack in the absence of residual stress.