Dynamic crack growth along a polymer composite-Homalite interface

Dynamic crack growth along the interface of a fiber-reinforced polymer composite-Homalite bimaterial subjected to impact shear loading is investigated experimentally and numerically. In the experiments, the polymer composite-Homalite specimens are impacted with a projectile causing shear dominated interfacial cracks to initiate and subsequently grow along the interface at speeds faster than the shear wave speed of Homalite. Crack growth is observed using dynamic photoelasticity in conjunction with high-speed photography. The calculations are carried out for a plane stress model of the experimental configuration and are based on a cohesive surface formulation that allows crack growth, when it occurs, to emerge as a natural outcome of the deformation history. The effect of impact velocity and loading rate is explored numerically. The experiments and calculations are consistent in identifying discrete crack speed regimes within which crack growth at sustained crack speeds is possible. We present the first conclusive experimental evidence of interfacial crack speeds faster than any characteristic elastic wave speed of the more compliant material. The occurrence of this crack speed was predicted numerically and the calculations were used to design the experiments. In addition, the first experimental observation of a mother-daughter crack mechanism allowing a subsonic crack to evolve into an intersonic crack is documented. The calculations exhibit all the crack growth regimes seen in the experiments and, in addition, predict a regime with a pulse-like traction distribution along the bond line.     

Sudden deceleration or acceleration of an intersonic shear crack

Sudden jumps in the crack tip velocity were revealed by numerical simulation (in both continuum/cohesive element and molecular dynamics approaches) and experiments for rapid shear cracking. The cracking velocity may accelerate from a sub-Rayleigh speed to the intersonic range, or from an intersonic speed to a higher one, when the reflected impact wave reloads the crack tip. On the other hand, the cracking velocity may decelerate from an intersonic speed to a lower one or recede to the sub-Rayleigh range when the fracture driving force declines. The velocity change encountered during intersonic cracking plays a different role from that in the acceleration or deceleration of a subsonic crack. A crack propagating at an intersonic speed would leave a shear wave trailing behind. When the crack decelerates or accelerates, the effect of the trailing wave will lead to a transition period from one steady-state solution of crack tip singularity to another. This investigation aims at quantifying these processes. The full field solution of an intersonic mode II crack whose speed changed suddenly from one velocity (intersonic or subsonic) to another (intersonic or subsonic) is given in closed form. The solution is facilitated via superposing a series of propagating crack problems that are loaded by dislocations to seal the unwanted crack-face sliding or by concentrated forces moving at various speeds to negate the crack-face traction. In contrast to the subsonic solution, the results in the intersonic case indicate that the elastic fields around the crack tip depend on the deceleration or acceleration history that is traced back over a long time. Singularity matching dictates the jump that may actually take place.     

Subsonic and intersonic mode II crack propagation with a rate-dependent cohesive zone

A recent experimental study has demonstrated the attainability of intersonic shear crack growth along weak planes in otherwise homogeneous, isotropic, linear elastic solids subjected to remote loading conditions (Rosakis et al., Science 284 (5418) (1999) 1337). The relevant experimental observations are summarized briefly here and the conditions governing the attainment of intersonic crack speeds are examined. Motivated by experimental observations, subsonic and intersonic mode II crack propagation with a rate-dependent cohesive zone is subsequently analyzed. A cohesive law is assumed, wherein the cohesive shear traction is either a constant or varies linearly with the local sliding rate. Complete decohesion is assumed to occur when the crack tip sliding displacement reaches a material-specific critical value. Closed form expressions are obtained for the near-tip fields. With a cohesive zone of finite size, it is found that the dynamic energy release rate is finite through out the intersonic regime. Crack tip stability issues are addressed and favorable speed regimes are identified. The influence of shear strength of the crack plane and of a rate parameter on crack propagation behavior is also investigated. The isochromatic fringe patterns predicted by the analytical solution are compared with the experimental observations of Rosakis et al. (1999) and comments are made on the validity of the proposed model.     

Transition of mode II cracks from sub-Rayleigh to intersonic speeds in the presence of favorable heterogeneity

Understanding sub-Rayleigh-to-intersonic transition of mode II cracks is a fundamental problem in fracture mechanics with important practical implications for earthquake dynamics and seismic radiation. In the Burridge-Andrews mechanism, an intersonic daughter crack nucleates, for sufficiently high prestress, at the shear stress peak traveling with the shear wave speed in front of the main crack. We find that sub-Rayleigh-to-intersonic transition and sustained intersonic propagation occurs in a number of other models that subject developing cracks to intersonic loading fields. We consider a spontaneously expanding sub-Rayleigh crack (or main crack) which advances, along a planar interface with linear slip-weakening friction, towards a place of favorable heterogeneity, such as a preexisting subcritical crack or a small patch of higher prestress (similar behavior is expected for a small patch of lower static strength). For a range of model parameters, a secondary dynamic crack nucleates at the heterogeneity and acquires intersonic speeds due to the intersonic stress field propagating in front of the main crack. Transition to intersonic speeds occurs directly at the tip of the secondary crack, with the tip accelerating rapidly to values numerically equal to the Rayleigh wave speed and then abruptly jumping to an intersonic speed. Models with favorable heterogeneity achieve intersonic transition and propagation for much lower prestress levels than the ones implied by the Burridge-Andrews mechanism and have transition distances that depend on the position of heterogeneity. We investigate the dependence of intersonic transition and subsequent crack propagation on model parameters and discuss implications for earthquake dynamics.     

Theoretical development and experimental validation of a thermally dissipative cohesive zone model for dynamic fracture of amorphous polymers

A thermally dissipative cohesive zone model is developed for predicting the temperature increase at the tip of a crack propagating dynamically in a nominally brittle material exhibiting a cohesive-type failure such as crazing. The model assumes that fracture energy supplied to the crack tip region that is in excess of that needed for the creation of new free surfaces during crack advance is converted to heat within the cohesive zone. Bulk dissipation mechanisms, such as plasticity, are not accounted for. Several cohesive traction laws are examined, and the model is then used to make predictions of crack tip heating at various crack propagation speeds in the nominally brittle amorphous polymer PMMA, observed to fail by a crazing-type mechanism. The heating predictions are compared to experimental data where the temperature field surrounding a high speed crack in PMMA was measured. Measurements are made in real time using a multi-point high speed HgCdTe infrared radiation detector array. At the same time as temperature, simultaneous measurement of fracture energy is made by a strain gauge technique, and crack tip speed is monitored through a resistance ladder method. Material strength can be estimated through uniaxial tension tests, thus minimizing the need for parameter fitting in the stress-opening traction law. Excellent agreement between experiments and theory is found for two of the cohesive traction law temperature predictions, but only for the case where a single craze is active during the dynamic fracture of PMMA, i.e. crack tip speed up to approximately 0.2c_R. For higher speed fracture where subsurface damage becomes prominent, the line dissipation model of a cohesive zone is inadequate, and a distributed damage model is needed.     

Dynamic behaviors of mode III interfacial crack under a constant loading rate

In an earlier study on intersonic crack propagation, Gao et al. (J. Mech. Phys. Solids 49: 2113–2132, 2001) described molecular dynamics simulations and continuum analysis of the dynamic behaviors of a mode II dominated crack moving along a weak plane under a constant loading rate. The crack was observed to initiate its motion at a critical time after the onset of loading, at which it is rapidly accelerated to the Rayleigh wave speed and propagates at this speed for a finite time interval until an intersonic daughter crack is nucleated at a peak stress at a finite distance ahead of the original crack tip. The present article aims to analyze this behavior for a mode III crack moving along a bi-material interface subject to a constant loading rate. We begin with a crack in an initially stress-free bi-material subject to a steadily increasing stress. The crack initiates its motion at a critical time governed by the Griffith criterion. After crack initiation, two scenarios of crack propagation are investigated: the first one is that the crack moves at a constant subsonic velocity; the second one is that the crack moves at the lower shear wave speed of the two materials. In the first scenario, the shear stress ahead of the crack tip is singular with exponent ?1/2, as expected; in the second scenario, the stress singularity vanishes but a peak stress is found to emerge at a distance ahead of the moving crack tip. In the latter case, a daughter crack supersonic with respect to the softer medium can be expected to emerge ahead of the initial crack once the peak stress reaches the cohesive strength of the interface.     

Intersonic crack propagation in homogeneous media under shear-dominated loading: theoretical analysis

The mechanics of cohesive failure under mixed-mode loading is investigated for the case of a steadily propagating subsonic and intersonic dynamic crack subjected to a follower tensile and shear distributed load. The cohesive failure model chosen in this study is rate independent but accounts for the coupling between normal and tangential damage. Special emphasis is placed here on mixed-mode cases with predominantly shear loading. The analysis shows that the size of the mixed-mode cohesive zone is smaller than that obtained in the pure shear case. The relative extent of the shear and tensile cohesive damage zones depends on the crack speed and the mode mixity. In the intersonic regime, the failure process takes place exclusively in shear, even under remote mixed-mode loading conditions.     

Molecular-dynamics simulation-based cohesive zone representation of intergranular fracture processes in aluminum

A traction-displacement relationship that may be embedded into a cohesive zone model for microscale problems of intergranular fracture is extracted from atomistic molecular-dynamics (MD) simulations. An MD model for crack propagation under steady-state conditions is developed to analyze intergranular fracture along a flat Σ99 [1 1 0] symmetric tilt grain boundary in aluminum. Under hydrostatic tensile load, the simulation reveals asymmetric crack propagation in the two opposite directions along the grain boundary. In one direction, the crack propagates in a brittle manner by cleavage with very little or no dislocation emission, and in the other direction, the propagation is ductile through the mechanism of deformation twinning. This behavior is consistent with the Rice criterion for cleavage vs. dislocation blunting transition at the crack tip. The preference for twinning to dislocation slip is in agreement with the predictions of the Tadmor and Hai criterion. A comparison with finite element calculations shows that while the stress field around the brittle crack tip follows the expected elastic solution for the given boundary conditions of the model, the stress field around the twinning crack tip has a strong plastic contribution. Through the definition of a Cohesive-Zone-Volume-Element—an atomistic analog to a continuum cohesive zone model element—the results from the MD simulation are recast to obtain an average continuum traction-displacement relationship to represent cohesive zone interaction along a characteristic length of the grain boundary interface for the cases of ductile and brittle decohesion.     

Dynamic crack deflection and penetration at interfaces in homogeneous materials: experimental studies and model predictions

We examine the deflection/penetration behavior of dynamic mode-I cracks propagating at various speeds towards inclined weak planes/interfaces of various strengths in otherwise homogeneous isotropic plates. A dynamic wedge-loading mechanism is used to control the incoming crack speeds, and high-speed photography and dynamic photoelasticity are used to observe, in real-time, the failure mode transition mechanism at the interfaces. Simple dynamic fracture mechanics concepts used in conjunction with a postulated energy criterion are applied to examine the crack deflection/penetration behavior and, for the case of interfacial deflection, to predict the crack tip speed of the deflected crack. It is found that if the interfacial angle and strength are such as to trap an incident dynamic mode-I crack within the interface, a failure mode transition occurs. This transition is characterized by a distinct, observable and predicted speed jump as well as a dramatic crack speed increase as the crack transitions from a purely mode-I crack to an unstable mixed-mode interfacial crack.     

Analysis of dynamic growth of a tensile crack in an elastic-plastic material

The influence of inertia on the stress and deformation fields near the tip of a crack growing in an elastic-plastic material is studied. The material is characterized by the von Mises yield criterion and J₂ flow theory of plasticity. The crack grows steadily under plane strain conditions in the tensile opening mode. Features of the stress and deformation state at points near the moving crack tip are described for elastic-perfectly plastic response and for several crack propagation speeds. It is found that inertia has a significant effect on the elastic-plastic response of material particles near the crack tip, and that elastic unloading may occur behind the crack tip for higher speeds. The relationship between the applied crack driving force, represented by a remote stress intensity factor, and the crack tip speed is examined on the basis of a critical crack tip opening angle growth criterion. The calculated result is compared with dynamic fracture toughness versus crack speed data for a 4340 steel.     

Stability of an intersonic shear crack to a perturbation of its edge

A weight function matrix is developed for obtaining the stress singularity coefficients at the edge of a plane crack, moving uniformly at an intersonic speed while subjected to arbitrary shear loading. This is then utilised for deriving, to first order, the perturbations of these coefficients associated with a small spatially and temporally varying perturbation of its edge. The perturbation solution is employed, in conjunction with a simple fracture criterion, to investigate the stability of a uniformly moving intersonic crack, subjected to following loads.     

A refined cohesive zone model that accounts for inertia of cohesive zone of a moving crack

Existing cohesive zone models assume that actual fracture zone of non-zero mass can be modeled by a line segment (cohesive zone) with no mass and inertia. In the present work, a simplified mass-spring model is presented to study inertia effect of cohesive zone on a mode-I steady-state moving crack. It is showed that fracture energy predicted by the present model increases dramatically when a finite limiting crack speed is approached. Reasonable agreement with known experiments indicates that the present model has the potential to catch the inertia effect of cohesive zone which has been ignored in existing cohesive zone models and better simulate dynamic fracture at high crack speed.     

界面裂纹萌生与扩展的分子动力学模拟

运用分子动力学模拟方法研究了裂纹在界面端处萌生与沿界面扩展的临界条件. 模拟考虑了一双相材料的3种模型，即构成90°/90°和 90°/180°夹角的两个界面端和一个界面裂纹. 模拟采用了包含原子区域与连续区域的并发型多尺度模型，即在界面端尖端和裂纹尖端附近 采用分子动力学(MD)方法，MD区域之外则按照线弹性有限元方法分析. 结果表明，在断裂启动时刻，3个模型沿界面的最大应力均达到界面理想强度；而且，其界 面能恰好足以克服界面材料的本征内聚能. 因此，界面端裂纹萌生与沿界面扩展的断裂条件可以通过界面理想强度和内聚能联系起来. 并基于模拟计算结果提出了界面断裂启动的统一准则.     

Effect of T-stress on branch angle of moving cracks

The effects of the T-stress on Yoffe crack propagation are analyzed. Using a maximum k_I fracture criterion near the kink of a moving crack tip, a branch angle is determined via asymptotic crack-tip field containing two fracture parameters related to singular and constant terms. Results indicate that crack speeds decrease the T-stress. The crack-tip field and the branch angle depend on the T-stress, especially for higher crack velocities. The critical speed for crack bifurcation is independent of remote transverse loading if neglecting the T-stress. Otherwise, the crack branch speed is reduced or raised, depending on positive or negative transverse loading, respectively.     

高铁锻钢制动盘热疲劳裂纹扩展仿真及寿命评价

巨大的制动热能导致高铁制动盘面出现热斑或热点,形成沿径向扩展的热疲劳裂纹,严重危及行车安全.通过仿真获得热疲劳裂纹扩展规律,据此制定出合适、经济的检测周期,具有重要意义.基于断裂力学理论,采用扩展有限单元法和自主研发的裂尖网格加密技术,得到热点形成规律,植入半椭圆形裂纹进行热疲劳断裂仿真,得到制动速度为400 km/h...     

On the path of a crack near a graded interface under large scale yielding

The trajectory of a crack lying parallel to a thin graded layer between two plastically dissimilar materials is studied using the exclusion region (ER) theory of fracture. The ER theory is a theoretical framework for surface separation within which a broad range of fracture phenomenologies can be represented. In the present study, the direction of crack advance is determined by maximizing the resolved normal-opening force on the near-tip region, whereas separation itself is governed by the intensity of plastic deformation near the tip. A computational study was undertaken using the ER theory. The special-purpose finite element analysis platform accommodates arbitrary––and a priori unknown––crack trajectories. The model problem considered herein involves two plastically dissimilar, but elastically identical, materials joined by a thin, graded interface layer. The initial crack lies parallel to the interface layer, and crack advance occurs under conditions of extensive plastic flow. It is found that the position of the initial crack relative to the interface layer has a strong influence on the fracture behavior. In general, the crack trajectories tend to curve toward the less-ductile material. Also, the presence of the interface layer leads to fracture toughnesses that significantly exceed those of either material individually for the configurations studied.     

塑性动态断裂实验研究

本文利用自制的实验装置,对韧性材料在爆炸冲击载荷作用下的塑性动态断裂特性,进行了实验研究,其中包括高塑性应变速率下,塑性区裂纹扩展过程和扩展速度的测试;塑性动态断裂韧性CTOD,及其在不同裂纹扩展速度下变化规律的测试。同时,对不同裂纹扩展速度的试件断口进行微观分析。     

内聚力模型裂纹问题分析的解析奇异单元

内聚力模型已经被广泛应用于需要考虑断裂过程区的裂纹问题当中,然而常用的数值方法应用于分析内聚力模型裂纹问题时还存在着一些不足,比如不能准确的给出断裂过程区的长度、需要网格加密等。为了克服这些缺点,论文构造了一个新型的解析奇异单元,并将之应用于基于内聚力模型的裂纹分析当中。首先将虚拟裂纹表面处的内聚力用拉格拉日插值的方法近似表示为多项式的形式,而多项式表示的内聚力所对应的特解可以被解析地给出。然后利用一个简单的迭代分析,基于内聚力模型的裂纹问题就可以被模拟出来了。最后,给出二个数值算例来证明本文方法的有效性。     

混凝土黏聚开裂模型若干进展

黏聚模型是用来描述混凝土断裂行为的基本模型, 首先介绍了混凝土的黏聚开裂模型的基本概念,总结了确定黏聚区的本构方程的各种方法,即直接单轴拉伸测试、J积分方法、R曲线法、柔度法和逆推法.然后介绍了黏聚模型在I型和复合型裂纹问题、疲劳断裂问题中的应用以及黏聚模型与混凝土尺寸效应的关系.最后对黏聚开裂模型与桥联模型、带状裂缝模型进行了比较和总结, 指出了该模型存在的问题, 并对其以后的发展方向提出了建议.      

On the importance of thermo-elastic cooling in the fracture of glassy polymers at high rates

In a previous thermo-mechanical analysis [Estevez, R., Basu, S., van der Giessen, E., 2005. Analysis of temperature effects near mode I cracks in glassy polymers. Int. J. Fract. 132, 249–273] in which shear yielding of the bulk and failure by crazing were accounted for, we examined which of these two viscoplastic processes contributed to heat in mode I fracture. The present study completes this work by investigating the conditions for thermo-elastic cooling prior to crack propagation as reported experimentally by Rittel [Rittel, D., 1998. Experimental investigation of transient thermo-elastic effects in dynamic fracture. Int. J. Solids Struct. 35, 2959–2973] and Bougaut and Rittel [Bougaut, O., Rittel, D., 2001. On crack tip cooling during dynamic crack propagation. Int. J. Solids Struct. 38, 2517–2532] on high strain rate loading of PMMA. To this end, coupled thermo-mechanical finite element simulations are carried out by accounting for the thermo-elastic source, in addition to the heat sources related to shear yielding and crazing. The bulk as well as cohesive zone parameters for crazing realistically describe PMMA as they are obtained from detailed calibration experiments. Our results show that if significant thermo-elastic cooling has to be observed in the vicinity of the crack tip of a polymeric material, suppression of shear yielding as well as suppression of crazing is necessary. It seems that at these high strain rates a brittle fracture mechanism activated at very high stresses takes over from crazing, or at least that craze initiation occurs for stress levels very different to those for quasi-static conditions.