Influence of temperature on a carbon–fibre epoxy composite subjected to static and fatigue loading under mode-I delamination

In this paper, interlaminar crack initiation and propagation under mode-I with static and fatigue loading of a composite material are experimentally assessed for different test temperatures. The material under study is made of a 3501-6 epoxy matrix reinforced with AS4 unidirectional carbon fibres, with a symmetric laminate configuration [0°]_16/S. In the experimental programme, DCB specimens were tested under static and fatigue loading. Based on the results obtained from static tests, fatigue tests were programmed to analyse the mode-I fatigue behaviour, so the necessary number of cycles was calculated for initiation and propagation of the crack at the different temperatures. G–N curves were determined under fatigue loading, N being the number of cycles at which delamination begins for a given energy release rate. G_ICmax–a, a–N and da/dN–a curves were also determined for different G_cr rates (90%, 85%, 75%, etc.) and different test temperatures: 90 °C, 50 °C, 20 °C, 0 °C, ?30 °C and ?60 °C.     

Experimental quantification of the plastic blunting process for stage II fatigue crack growth in one-phase metallic materials

The plastic blunting process during stage II fatigue crack growth was studied in pure polycrystalline Ni to investigate effects of strain localization and inelastic behavior on the kinematics of crack advance. Correlations were obtained between strain fields ahead of a fatigue crack, crack advance per cycle and crack growth kinetics. Strain fields were quantified using a combination of in situ loading experiments, scanning electron microscopy and digital image correlation for 8 < ΔK < 20 MPa m^1/2 and a fixed load ratio of 0.1. Results indicate that strain localized along a dominant deformation band, which was usually crystallographic and carried mostly pure shear for large loads and was of mixed character for lower loads. Instances of double deformation bands were observed, with bands acting either in a simultaneous or alternating fashion. It was found that the area integral of the opening strain for values larger than a given threshold, an “integrated” strain, had a power-law relationship with ΔK, with the exponent approximately equal to the Paris exponent (m). Therefore, the crack growth rate was proportional to the integrated strain. An analysis based on this correlation and the presence of dominant shear bands indicated that the integrated strain is related to the accumulated displacement in the band. This, in turn, is proportional to the product of the cyclic plastic zone radius and the average shear strain ahead of the tip, which represents a basic length scale for plastic blunting. Assumptions on the load dependence of these quantities, based on their observed spatial variation, allowed estimating $m=2\left(1+\frac{1}{1+n^{\prime }}\right)$, where n′ is the cyclic hardening exponent (0 < n′ < 1). This gives 3 < m < 4, which accounts for about 50% of the observed values of m between 1.5 and 6 for a wide variety of metallic materials.     

Numerical study on coalescence of two pre-existing coplanar flaws in rock

Crack propagation and coalescence processes are the fundamental mechanisms leading to progressive failure processes in rock masses, in which parallel non-persistent rock joints are commonly involved. The coalescence behavior of the latter, which are represented as pre-existing coplanar flaws (cracks), is numerically investigated in the present study. By using AUTODYN as the numerical tool, the present study systematically simulates the coalescence of two pre-existing coplanar flaws in rock under compression. The cumulative damage failure criterion is adopted in the numerical models to simulate the cumulative damage process in the crack initiation and propagation. The crack types (shear or tensile) are identified by analyzing the mechanics information associated with the crack initiation and propagation processes. The simulation results, which are generally in a good accordance with physical experimental results, indicate that the ligament length and the flaw inclination angle have a great influence on the coalescence pattern. The coalescence pattern is relatively simple for the flaw arrangements with a short ligament length, which becomes more complicated for those with a long ligament length. The coalescence trajectory is composed of shear cracks only when the flaw inclination angle is small (such as β ? 30°). When the pre-existing flaws are steep (such as β ? 75°), the coalescence trajectory is composed of tensile cracks as well as shear cracks. When the inclination angle is close to the failure angle of the corresponding intact rock material, and the ligament length is not long (such as L ? 2a), the direct shear coalescence is the more favorable coalescence pattern. In the special case that the two pre-existing flaws are vertical, the model will have a direct tensile coalescence pattern when the ligament length is short (L ? a), while the coalescence between the two inner flaw tips is not easy to achieve if the ligament length is long (L ? 2a).     

Constant moving crack in a magnetoelectroelastic material under anti-plane shear loading

Analytical solutions for an anti-plane Griffith moving crack inside an infinite magnetoelectroelastic medium under the conditions of permeable crack faces are formulated using integral transform method. The far-field anti-plane mechanical shear and in-plane electrical and magnetic loadings are applied to the magnetoelectroelastic material. Expressions for stresses, electric displacements and magnetic inductions in the vicinity of the crack tip are derived. Field intensity factors for magnetoelectroelastic material are obtained. The stresses, electric displacements and magnetic inductions at the crack tip show inverse square root singularities. The moving speed of the crack have influence on the dynamic electric displacement intensity factor (DEDIF) and the dynamic magnetic induction intensity factor (DMIIF), while the dynamic stress intensity factor (DSIF) does not depend on the velocity of the moving crack. When the crack is moving at very lower or very higher speeds, the crack will propagate along its original plane; while in the range of M_c1 < M < M_c2, the propagation of the crack possibly brings about the branch phenomena in magnetoelectroelastic media.     

Periodic group crack problems in an infinite plate

This paper investigates periodic group crack problems in an infinite plate. The periodic group crack is composed of infinite groups with numbering from j = −∞, …, −2, −1, 0, 1, 2, …, to j = ∞, and the groups are placed periodically. The same loading condition and the same geometry are assumed for cracks in all groups. A singular integral equation is used to solve the problems. The singular integral equation is formulated on cracks of the 0th group (or the central group) with the collection of influences from the infinite groups. The influences of many neighboring groups to the central group are evaluated exactly. Meantime, the influences of many remote groups to the central group can be summed up into one term approximately. The stress intensity factors at crack tips can be evaluated from the solution of the singular integral equation. It is found from some sample problems that the obtained results are very accurate. Finally, several numerical examples are presented and interaction among the group cracks is addressed.     

A cohesive zone model for fatigue crack growth allowing for crack retardation

A damage-based cohesive model is developed for simulating crack growth due to fatigue loading. The cohesive model follows a linear damage-dependent traction–separation relation coupled with a damage evolution equation. The rate of damage evolution is characterized by three material parameters corresponding to common features of fatigue behavior captured by the model, namely, damage accumulation, crack retardation and stress threshold. Good agreement is obtained between finite element solutions using the model and fatigue test results for an aluminum alloy under different load ratios and for the overload effect on ductile 316 L steel.     

A reexamination of plasticity-induced crack closure in fatigue crack propagation

The crack closure concept is often used to consider the R-ratio and overload effects on fatigue crack growth. The presumption is that when the crack is closed, the external load produces negligible fatigue damage in the cracked component. The current investigation provides a reassessment of the frequently used concept with an emphasis on the plasticity-induced crack closure. A center cracked specimen made of 1070 steel was investigated. The specimen was subjected to plane-stress mode I loading. An elastic–plastic stress analysis was conducted for the cracked specimens using the finite element method. By applying the commonly used one-node-per-cycle debonding scheme for the crack closure simulations, it was shown that the predicted crack opening load did not stabilize when the extended crack was less than four times of the plastic zone size. The predicted opening load was strongly influenced by the plasticity model used. When the elastic–perfectly plastic (EPP) stress–strain relationship was used together with the kinematic hardening plasticity theory, the predicted crack opening load was found to be critically dependent on the element size of the finite element mesh model. For R = 0, the predicted crack opening load was greatly reduced when the finite element size became very fine. The kinematic hardening rule with the bilinear (BL) stress–strain relationship predicted crack closure with less dependence on the element size. When a recently developed cyclic plasticity model was used, the element size effect on the predicted crack opening level was insignificant. While crack closure may occur, it was demonstrated that cyclic plasticity persisted in the material near the crack tip. The cyclic plasticity was reduced but not negligible when the crack was closed. The traditional approaches may have overestimated the effect of crack closure in fatigue crack growth predictions.     

Asymptotic fields around an interfacial crack with a cohesive zone ahead of the crack tip

This paper considers an interfacial crack with a cohesive zone ahead of the crack tip in a linearly elastic isotropic bi-material and derives the mixed-mode asymptotic stress and displacement fields around the crack and cohesive zone under plane deformation conditions (plane stress or plane strain). The field solution is obtained using elliptic coordinates and complex functions and can be represented in terms of a complete set of complex eigenfunction terms. The imaginary portion of the eigenvalues is characterized by a bi-material mismatch parameter ε = arctanh(β)/π, where β is a Dundurs parameter, and the resulting fields do not contain stress singularity. The behaviors of “Mode I” type and “Mode II” type fields based on dominant eigenfunction terms are discussed in detail. For completeness, the counterpart for the Mode III solution is included in an appendix.     

A comparison of cohesive zone modeling and classical fracture mechanics based on near tip stress field

A mode III crack with a cohesive zone in a power-law hardening material is studied under small scale yielding conditions. The cohesive law follows a softening path with the peak traction at the start of separation process. The stress and strain fields in the plastic zone, and the cohesive traction and separation displacement in the cohesive zone are obtained. The results show that for a modest hardening material (with a hardening exponent N = 0.3), the stress distribution in a large portion of the plastic zone is significantly altered with the introduction of the cohesive zone if the peak cohesive traction is less than two times yield stress, which implies the disparity in terms of the fracture prediction between the classical approach of elastic–plastic fracture mechanics and the cohesive zone approach. The stress distributions with and without the cohesive zone converge when the peak cohesive traction becomes infinitely large. A qualitative study on the equivalency between the cohesive zone approach and the classical linear elastic fracture mechanics indicates that smaller cracks require a higher peak cohesive traction than that for longer cracks if similar fracture initiations are to be predicted by the two approaches.     

Effect of gravity on the mechanical properties of lunar regolith tested using a low gravity simulation device

Simulations of the bearing capacity and shear strength of regolith under Earth’s gravity produce different results from those under low gravity. A low-gravity simulation device was developed in this study, and an internal stress model of regolith simulant was established to correct the errors. The model revealed additional force on both shear plane in the shear test and the press plate area in the pressure–sinkage test. The sinkage and shear test results showed that low gravity decreased the deformable index n, frictional modulus k_φ and cohesion c, whereas there were no obvious changes to the cohesive modulus k_c and internal friction angle φ. The sinkage generally increased as the gravity decreased under a consistent normal load larger than 50 N, but when the wheel load was lower than 50 N, the sinkage of the TYII-1 simulant was larger under 1 G than 1/6 G. Gravity had little effect on the shear strength of the regolith. However, the tractive thrust of the TYII-1 simulant was lower under 1/6 G than 1 G. The smaller difference was due to differences in the way the soils responded to changes in the gravity level for the TYII-2 simulant.     

Fatigue properties of age-hardened Al alloy 2017-T4 under ultrasonic loading

Fatigue properties of age-hardened Al alloy 2017-T4 under ultrasonic loading frequency (20 kHz) were investigated and compared with the results under conventional loading of rotating bending (50 Hz). The growth of a crack retarded at about 500 µm in surface length under ultrasonic loading, while at about 20 µm under rotating bending. Although striations being a typical fracture mechanism were observed under conventional loading, most of fracture surface was covered with many facets under ultrasonic loading. These facets were also observed under rotating bending in nitrogen gas. The difference in growth mechanism depending on the loading frequency and the retardation of a crack growth under ultrasonic loading may be caused by the environment at the crack tip due to high crack growth rate under ultrasonic loading.     

Use of equivalent body forces for acoustic emission from a crack in a plate

A method to determine acoustic emission of surface waves from a crack near the free edge of a plate, is presented, in terms of the function f(t), which defines the time dependence of the crack opening process, the crack opening volume per unit thickness of the plate, and the elastic constants of the plate. The determination of the time-varying displacement is based on the use of equivalent body forces, which are shown to be two double forces. The acoustic emission of the crack, or the equivalent radiation from the double forces, has been obtained by a novel use of the elastodynamic reciprocity theorem. It is of interest that the normal surface-wave displacement at a position x₀ of the free edge comes out as depending on df/dt evaluated at x₀ for t > x₀/c_R, where c_R is the velocity of surface waves on the free edge.     

Lattice orientation effects on void growth and coalescence in fcc single crystals

Void growth and coalescence in fcc single crystals were studied using crystal plasticity under uniaxial and biaxial loading conditions and various orientations of the crystalline lattice. A 2D plane strain unit cell with one and two cylindrical voids was employed using three-dimensional 12 potentially active slip systems. The results were compared to five representative orientations of the tensile axis on the stereographic triangle. For uniaxial tension conditions, the void volume fraction increase under the applied load is strongly dependent on the crystallographic orientation with respect to the tensile axis. For some orientations of the tensile axis, such as [1 0 0] or [1 1 0], the voids exhibited a growth rate twice as fast compared with other orientations ([1 0 0], [2 1 1]). Void growth and coalescence simulations under uniaxial loading indicated that during deformation along some orientations with asymmetry of the slip systems, the voids experienced rotation and shape distortion, due mainly to lattice reorientation. Coalescence effects are shown to diminish the influence of lattice orientation on the void volume fraction increase, but noteworthy differences are still present. Under biaxial loading conditions, practically all differences in the void volume fraction for different orientations of the tensile axes during void growth vanish. These results lead to the conclusion that at microstructural length scales in regions under intense biaxiality/triaxiality conditions, such as crack tip or notched regions, the plastic anisotropy due to the initial lattice orientation has only a minor role in influencing the void growth rate. In such situations, void growth and coalescence are mainly determined by the stress triaxiality, the magnitude of accumulated strain, and the spatial localization of such plastic strains.     

Coplanar propagation paths of 3D cracks in infinite bodies loaded in shear

Bower and Ortiz, recently followed by Lazarus, developed a powerful method, based on a theoretical work of Rice, for numerical simulation of planar propagation paths of mode 1 cracks in infinite isotropic elastic bodies. The efficiency of this method arose from the need for the sole 1D meshing of the crack front. This paper presents an extension of Rice’s theoretical work and the associated numerical scheme to mixed-mode (2 + 3) shear loadings. Propagation is supposed to be channeled along some weak planar layer and to remain therefore coplanar, as in the case of a geological fault for instance. The capabilities of the method are illustrated by computing the propagation paths of cracks with various initial contours (circular, elliptic, rectangular, heart-shaped) in both fatigue and brittle fracture. The crack quickly reaches a stable, almost elliptic shape in all cases. An approximate but accurate analytic formula for the ratio of the axes of this stable shape is derived.     

Cracked infinite cylinder with two rigid inclusions under axisymmetric tension

This paper considers the problem of an axisymmetric infinite cylinder with a ring shaped crack at z = 0 and two ring-shaped rigid inclusions with negligible thickness at z = ±L. The cylinder is under the action of uniformly distributed axial tension applied at infinity and its lateral surface is free of traction. It is assumed that the material of the cylinder is linearly elastic and isotropic. Crack surfaces are free and the constant displacements are continuous along the rigid inclusions while the stresses have jumps. Formulation of the mixed boundary value problem under consideration is reduced to three singular integral equations in terms of the derivative of the crack surface displacement and the stress jumps on the rigid inclusions. These equations, together with the single-valuedness condition for the displacements around the crack and the equilibrium equations along the inclusions, are converted to a system of linear algebraic equations, which is solved numerically. Stress intensity factors are calculated and presented in graphical form.     

An experimental and numerical study of the localization behavior of tantalum and stainless steel

In general, the shear localization process involves initiation and growth. Initiation is expected to be a stochastic process in material space where anisotropy in the elastic–plastic behavior of single crystals and inter-crystalline interactions serve to form natural perturbations to the material’s local stability. A hat-shaped sample geometry was used to study shear localization growth. It is an axi-symmetric sample with an upper “hat” portion and a lower “brim” portion with the shear zone located between the hat and brim. The shear zone length is 870–890 μm with deformation imposed through a split-Hopkinson pressure bar system at maximum top-to-bottom velocity in the range of 8–25 m/s. We present experimental results of the deformation response of tantalum and 316L stainless steel samples. The tantalum samples did not form shear bands but the stainless steel sample formed a late stage shear band. We have also modeled these experiments using both conductive and adiabatic continuum models. An anisotropic elasto-viscoplastic constitutive model with damage evolution was used within the finite element code EPIC. A Mie-Gruneisen equation of state and the rate and temperature sensitive MTS flow stress model together with a Gurson flow surface were employed. The models performed well in predicting the experimental data. The numerical results for tantalum suggested a maximum equivalent strain rate on the order of 7 × 10⁴ s⁻¹ in the gage section for an imposed top surface displacement rate of 17.5 m/s. The models also suggested that for an initial temperature of 298 K a temperature in the neighborhood of 900 K was reached within the shear section. The numerical results for stainless steel suggest that melting temperature was reached throughout the shear band shortly after peak load. Due to sample geometry, the stress state in the shear zone was not pure shear; a significant normal stress relative to the shear zone basis line was developed.     

Thermo-mechanical behaviour of TRIP 1000 steel sheets subjected to low velocity perforation by conical projectiles at different temperatures

This paper presents and analyzes the behaviour of TRIP 1000 steel sheets subjected to low velocity perforation by conical projectiles. The relevance of this material resides in the potential transformation of retained austenite to martensite during impact loading. This process leads to an increase in strength and ductility of the material. However, this transformation takes place only under certain loading conditions strongly dependent on the initial temperature and deformation rate. In order to study the material behaviour under impact loading, perforation tests have been performed using a drop weight tower. Experiments were carried out at two different initial temperatures T₀ = 213 K and T₀ = 288 K, and within the range of impact velocities 2.5 m/s ? V₀ ? 4.5 m/s. The experimental setup enabled the measuring of impact velocity, residual velocity, load-time history and failure mode. In addition, dry and lubricated contacts between the striker and the plate have been investigated. Finally, by using X-ray diffraction it has been shown that no martensitic transformation takes place during the perforation process. The causes involving the none-appearance of martensite are examined.     

Dynamic fracture of concrete – compact tension specimen

The behavior of concrete structures is strongly influenced by the loading rate. Compared to quasi-static loading concrete loaded by impact loading acts in a different way. First, there is a strain-rate influence on strength, stiffness, and ductility, and, second, there are inertia forces activated. Both influences are clearly demonstrated in experiments. Moreover, for concrete structures, which exhibit damage and fracture phenomena, the failure mode and cracking pattern depend on loading rate. In general, there is a tendency that with the increase of loading rate the failure mode changes from mode-I to mixed mode. Furthermore, theoretical and experimental investigations indicate that after the crack reaches critical speed of propagation there is crack branching. The present paper focuses on 3D finite-element study of the crack propagation of the concrete compact tension specimen. The rate sensitive microplane model is used as a constitutive law for concrete. The strain-rate influence is captured by the activation energy theory. Inertia forces are implicitly accounted for through dynamic finite element analysis. The results of the study show that the fracture of the specimen strongly depends on the loading rate. For relatively low loading rates there is a single crack due to the mode-I fracture. However, with the increase of loading rate crack branching is observed. Up to certain threshold (critical) loading rate the maximal crack velocity increases with increase of loading rate, however, for higher loading rates maximal velocity of the crack propagation becomes independent of the loading rate. The critical crack velocity at the onset of crack branching is found to be approximately 500 m/s.