A Theory of Crack Initiation and Growth in Viscoelastic Media !II

Total Page:16

File Type:pdf, Size:1020Kb

A Theory of Crack Initiation and Growth in Viscoelastic Media !II International Journal of Fracture, Vol. 11. No 4 August 1975 Noordhoff International Publishing Leyden 549 Printed in The Netherlands A theory of crack initiation and growth in viscoelastic media !II. Analysis of continuous growth* R. A. SCHAPERY Civil and Aerospace Engineerin9 Departments, Texas A & M University, College Station, Texas 77843, U.S.A. (Received April 19, 1973; in revised form February 21, 1975) ABSTRACT The theory of crack growth developed in Parts I and II is used to predict crack velocity and failure time for an elastomer under simple uniaxial and biaxial stress states. Included are consideration of the effect of specimen size on failure time for initially small cracks, experimental determination of fracture properties, the effect of strain level on crack propagation, and the validity of the plane strain assumption. The second half of the paper examines the effect upon crack velocity of nonlinear viscoelasticityof failing material at the crack tip in media which is otherwise linearly viscoelastic. The dependence of fracture behavior on environmental changes and aging is also considered. I. Introduction In this paper we examine the validity of the opening-mode crack growth theory developed in Parts I [1] and II [2] by comparing predictions of crack velocity and failure time with experimental data obtained by others. Also, certain practical implications of the theory are explored. In Section 2, experimental and theoretical results for the time required for complete failure of centrally cracked elastomeric sheets (Solithane) are compared and found to be in good agreement over a range of constant applied loads and temperatures. Similar agreement is shown in Section 3 for the crack velocity in long, narrow strips of the same elastomer under constant applied deformation; the same mechanical and fracture properties are used in the analysis of both geometries. At relatively large overall strains some discre- pancy between theory and experiment is noted, and a strain-modified stress intensity factor is found to remove this error. In both of these problems, stress intensity factors for plane stress are used in a crack growth equation which is based on the assumption of plane strain at points close to the crack tip; justification for this procedure is given in Section 4. For the particular elastomer investigated, fracture energy and a measure of the strength of the failing material at the crack tip are taken to be constant. However, the underlying theory developed in Parts I and II is not restricted to constant properties. In Section 5 we examine the influence of nonlinear viscoelasticity of the failing material on crack growth behavior in viscoelastic media. Effects of environmental changes and aging are briefly described in Section 6. For the prediction of crack velocity and failure time of the elastomer, we shall use a result obtained in [2, Eqn. (56)]: da/dl=(TC/2) F C2 ][/m ~I/m{]~2(l+i/m,,/~2 12 1_8F(1-K?/K?,)_] "~m ,--t ,, ~,.-, (1) This expression for crack velocity da/dt is based on the use of a creep compliance Cv(t) which is written in the form of a generalized power law : Cv(t) = Co --}- C2 t m (2) where Co=C~(0); C2 and m are positive but not necessarily constant. For a constant * Parts I and II are published, respectively, in Vol. 11, No. 1 (1975) 141-159, and Vol. 11, No. 3 (1975) 369-388. Int. Journ. of Fracture, 11 (1975) 549-562 550 R.A. Schapery Poisson's ratio v, this compliance is simply related to the creep compliance D (t) under uni- axial stress [1, Eqn. (38)], c (t) = 4(1 -,,2)o (t) (3) where D (t) is the axial strain for a unit uniaxial stress applied at t = 0. Additionally, F = frac- ture energy, am= maximum stress within the failing material near the tip, I x =integral measure of the shape of the stress distribution in the failing material (0 < 11 =< 2, with I1 = 2 when this stress is uniformily distributed) and 2~/"=parameter depending only on m (:1/,,~ 1 for O<m< 1). The opening-mode stress intensity factor is denoted by K~, as defined by the relation K, = lim~,~o (2n¢1) ~ a ° (4) where ayo= ao (~1) is the normal stress acting across the plane of crack prolongation and ~1 is the distance ahead of the crack tip; by definition, a ° is the singular stress corresponding to a mathematically sharp crack tip without a finite zone of failing material. Also, K,g =_ (st~Co) ~ (5) is identical to the critical stress intensity factor for an elastic body having a constant compli- ance Co. The compliance parameters, Co, C2 and m, and the fracture properties, F and am I1, will be assumed constant in Sections 2 and 3. 2. Crack growth and failure of a sheet under constant load Consider a centrally cracked plate, such as shown in Fig. 1. Assuming the instantaneous crack length, 2a, is large compared to sheet thickness but small compared to the in-plane dimensions, the stress intensity factor is approximately [e.g., 3], K, = (an) ~ a (6) where a is the applied stress. Substituting Eqn. (6) into Eqn. (1), assuming a is applied at t = 0 and is constant thereafter, and integrating yields t = (2/n 2) (a~ I21/a 2) (Co/2,, C2) 1/rn Ial/a°[(Oy/O0) ') -- 1]l/m(d~}/) ') (7a) where 2a o is the initial crack length, and a o is defined as a o =-- 8F/nCo a 2 . (7b) When the central crack reaches length 2ag (assuming 2a o < 2ao) complete failure of the sheet occurs; this point follows from Eqns. (1), (5), (6) and (7b) which imply fi--*~ as a--*ag. The time at which a = a 9 is called the failure time, t/, which, after changing the integration variable in Eqn. (7a), becomes t r = (2/n2)(a2 12/a2) (C0/2," C2)1/,, J'~o/o,(1 - u) 1/m U -(1 + 1/m, du. (8) The integral can be expressed in terms of incomplete Beta functions. It is of interest to notice that when ao/a o >> 1 and 0_< m < 1 the integrand in Eqn. (7a) goes quickly to zero as 7 increases. In fact, the current time is already 90 ~o of t f when Int. Journ. of Fracture, 11 (1975) 549 562 Theory of crack initiation and gru~th m ui.s~,ul, t,tiu media. III 551 a/a o -~ I(Y". (9) One can easily show that if the assumption ao/a o >> 1 is removed, the crack growth at t = 0.90 t/will be even less than given in Eqn. (9). These observations are important because they tell us if the stress intensity factor, Eqn. (6), is valid over most of the time period 0_< t _< t/ or if the effect of finite sheet width must be included when predicting failure time. For three typical values of m, Eqn. (9) yields 10 , m= 1 a/a° ~- i 3.3, m=0.5 (10) ( 1.6, m=0.2. Surprisingly, most of the time required for failure is consumed while the crack is still relatively small, especially for 0< m< 0.5. As a point of comparison, if the bracketed term in the integrand in Eqn. (7a) is replaced by (ag/ao7) TM for all 0_< t< t/, the failure time is derived by setting a/ao= oo to find tf= [(2ma2~I~)/Tz(2+ 1/"'](8r/2mCzao)t/ma 2,1 + l/m) (11) which brings out very simply the effect of stress on failure time. This is the same result as would be obtained by neglecting (KZ/KZg) in Eqn. (1) at the outset of the analysis. ///////// ~ 0.25" oJ 13 0 ° C "~l t° • 5 ° C ~eO I - 0 I0 ° C ~Z~e ./--The°ry' -- O • 15 ° C ~/ (..9 O A 25 ° C ._1 • :30 ° C 0 40" C '%. I I I \l O( 2 4 6 8 LOGio tf/aT, SEC. Figure 1. Experimental and theoretical failure time t~. as a function of the applied gross stress a for Solithane 50/50. Theory: Eqn. (13). Experimental data after Knauss [3]. Equations (8) and (11) will be compared to experimental data on Solithane 50/50, which is a crosslinked, amorphous, polyurethane rubber. The experimental data for failure time, which are shown in Fig. 1, have been normalized with respect to the stress cry, where ao~ =- (8FE~c/3~zao) ~ (12) Int. Journ. of Fracture, 11 (1975) 549 562 552 R. A. Schapery is the critical stress for the onset of crack growth in an elastic plate having Poisson's ratio ~' = ½ and Young's modulus Eoo = 4 (1 - v2)/Cv (~) = 3/C~ (~); thus, no crack growth occurs for log (a/aoo)2< O. Also, all data have been reduced to a common temperature of 0°C by recognizing that the material is essentially thermorheologically simple [3]; according to Section 6, the major effect of temperature can be introduced by replacing t/with t//aT in Eqn. (8). Since log (t/aT)---- log t-- log aT, the theory, if valid, implies experimental failure data obtained at different temperatures can be superposed by means of horizontal translations of magnitude log aT SO as to form a single curve; the failure data at the temperatures indicated in Fig. 1 were shifted to the data at 0°C. (One can interpret the reduced failure data as being the results of tests conducted entirely at 0~ C.) We should add that the uniaxial creep compliance for this material, which is shown in Fig. 2, was formed in the same way; in fact, the values of log a T used to shift the failure data are those obtained by shifting the creep compliance.
Recommended publications
  • How Similitude Can Improve Our Understanding of Fatigue Delamination Growth
    Misinterpreting the Results: How Similitude can Improve our Understanding of Fatigue Delamination Growth *Submitted for publishing in Composites Science and Technology on March 25 2010 Calvin Rans, René Alderliesten, Rinze Benedictus Structural Integrity Group, Faculty of Aerospace Engineering, Delft University of Technology, P.O. Box 5058, 2600 GB Delft, The Netherlands Abstract: Use of the strain energy release rate for characterizing delamination growth in composite and bonded structures is now commonplace. Although its use is common, a consensus on the best formulation of strain energy release rate has not been reached for describing fatigue delamination growth. Most commonly selected formulations are the maximum strain energy release rate and the strain energy release rate range. This paper examines the use strain energy release rate range and how a misconception of its definition is misleading our understanding of fatigue delamination growth. Keywords: Delamination, fatigue, strain energy release rate, similitude 1. INTRODUCTION The study and characterization of delamination growth behaviour has gained significant attention in the research community due to the growing use of composite structures in a wide range of structural applications. A widely utilized approach is to apply a linear elastic fracture mechanics (LEFM) description using the strain energy release rate (G) in a manner analogous to metal crack growth in terms of the stress intensity factor (K). The primary reason for the use of G rather than K is that the local crack-tip stresses used to determine K are difficult to obtain in inhomogeneous composite laminates [1, 2]. Use of G avoids the need to determine crack-tip stresses.
    [Show full text]
  • Cohesive Crack with Rate-Dependent Opening and Viscoelasticity: I
    International Journal of Fracture 86: 247±265, 1997. c 1997 Kluwer Academic Publishers. Printed in the Netherlands. Cohesive crack with rate-dependent opening and viscoelasticity: I. mathematical model and scaling ZDENEKÏ P. BAZANTÏ and YUAN-NENG LI Department of Civil Engineering, Northwestern University, Evanston, Illinois 60208, USA Received 14 October 1996; accepted in revised form 18 August 1997 Abstract. The time dependence of fracture has two sources: (1) the viscoelasticity of material behavior in the bulk of the structure, and (2) the rate process of the breakage of bonds in the fracture process zone which causes the softening law for the crack opening to be rate-dependent. The objective of this study is to clarify the differences between these two in¯uences and their role in the size effect on the nominal strength of stucture. Previously developed theories of time-dependent cohesive crack growth in a viscoelastic material with or without aging are extended to a general compliance formulation of the cohesive crack model applicable to structures such as concrete structures, in which the fracture process zone (cohesive zone) is large, i.e., cannot be neglected in comparison to the structure dimensions. To deal with a large process zone interacting with the structure boundaries, a boundary integral formulation of the cohesive crack model in terms of the compliance functions for loads applied anywhere on the crack surfaces is introduced. Since an unopened cohesive crack (crack of zero width) transmits stresses and is equivalent to no crack at all, it is assumed that at the outset there exists such a crack, extending along the entire future crack path (which must be known).
    [Show full text]
  • An Investigation Into the Effects of Viscoelasticity on Cavitation Bubble Dynamics with Applications to Biomedicine
    An Investigation into the Effects of Viscoelasticity on Cavitation Bubble Dynamics with Applications to Biomedicine Michael Jason Walters School of Mathematics Cardiff University A thesis submitted for the degree of Doctor of Philosophy 9th April 2015 Summary In this thesis, the dynamics of microbubbles in viscoelastic fluids are investigated nu- merically. By neglecting the bulk viscosity of the fluid, the viscoelastic effects can be introduced through a boundary condition at the bubble surface thus alleviating the need to calculate stresses within the fluid. Assuming the surrounding fluid is incompressible and irrotational, the Rayleigh-Plesset equation is solved to give the motion of a spherically symmetric bubble. For a freely oscillating spherical bubble, the fluid viscosity is shown to dampen oscillations for both a linear Jeffreys and an Oldroyd-B fluid. This model is also modified to consider a spherical encapsulated microbubble (EMB). The fluid rheology affects an EMB in a similar manner to a cavitation bubble, albeit on a smaller scale. To model a cavity near a rigid wall, a new, non-singular formulation of the boundary element method is presented. The non-singular formulation is shown to be significantly more stable than the standard formulation. It is found that the fluid rheology often inhibits the formation of a liquid jet but that the dynamics are governed by a compe- tition between viscous, elastic and inertial forces as well as surface tension. Interesting behaviour such as cusping is observed in some cases. The non-singular boundary element method is also extended to model the bubble tran- sitioning to a toroidal form.
    [Show full text]
  • Analysis of Viscoelastic Flexible Pavements
    Analysis of Viscoelastic Flexible Pavements KARL S. PISTER and CARL L. MONISMITH, Associate Professor and Assistant Professor of Civil Engineering, University of California, Berkeley To develop a better understandii^ of flexible pavement behavior it is believed that components of the pavement structure should be considered to be viscoelastic rather than elastic materials. In this paper, a step in this direction is taken by considering the asphalt concrete surface course as a viscoelastic plate on an elastic foundation. Assumptions underlying existing stress and de• formation analyses of pavements are examined. Rep• resentations of the flexible pavement structure, im• pressed wheel loads, and mechanical properties of as• phalt concrete and base materials are discussed. Re• cent results pointing toward the viscoelastic behavior of asphaltic mixtures are presented, including effects of strain rate and temperature. IJsii^ a simplified viscoelastic model for the asphalt concrete surface course, solutions for typical loading problems are given. • THE THEORY of viscoelasticity is concerned with the behavior of materials which exhibit time-dependent stress-strain characteristics. The principles of viscoelasticity have been successfully used to explain the mechanical behavior of high polymers and much basic work U) lias been developed in this area. In recent years the theory of viscoelasticity has been employed to explain the mechanical behavior of asphalts (2, 3, 4) and to a very limited degree the behavior of asphaltic mixtures (5, 6, 7) and soils (8). Because these materials have time-dependent stress-strain characteristics and because they comprise the flexible pavement section, it seems reasonable to analyze the flexible pavement structure using viscoelastic principles.
    [Show full text]
  • Engineering Viscoelasticity
    ENGINEERING VISCOELASTICITY David Roylance Department of Materials Science and Engineering Massachusetts Institute of Technology Cambridge, MA 02139 October 24, 2001 1 Introduction This document is intended to outline an important aspect of the mechanical response of polymers and polymer-matrix composites: the field of linear viscoelasticity. The topics included here are aimed at providing an instructional introduction to this large and elegant subject, and should not be taken as a thorough or comprehensive treatment. The references appearing either as footnotes to the text or listed separately at the end of the notes should be consulted for more thorough coverage. Viscoelastic response is often used as a probe in polymer science, since it is sensitive to the material’s chemistry and microstructure. The concepts and techniques presented here are important for this purpose, but the principal objective of this document is to demonstrate how linear viscoelasticity can be incorporated into the general theory of mechanics of materials, so that structures containing viscoelastic components can be designed and analyzed. While not all polymers are viscoelastic to any important practical extent, and even fewer are linearly viscoelastic1, this theory provides a usable engineering approximation for many applications in polymer and composites engineering. Even in instances requiring more elaborate treatments, the linear viscoelastic theory is a useful starting point. 2 Molecular Mechanisms When subjected to an applied stress, polymers may deform by either or both of two fundamen- tally different atomistic mechanisms. The lengths and angles of the chemical bonds connecting the atoms may distort, moving the atoms to new positions of greater internal energy.
    [Show full text]
  • Oversimplified Viscoelasticity
    Oversimplified Viscoelasticity THE MAXWELL MODEL At time t = 0, suddenly deform to constant displacement Xo. The force F is the same in the spring and the dashpot. F = KeXe = Kv(dXv/dt) (1-20) Xe is the displacement of the spring Xv is the displacement of the dashpot Ke is the linear spring constant (ratio of force and displacement, units N/m) Kv is the linear dashpot constant (ratio of force and velocity, units Ns/m) The total displacement Xo is the sum of the two displacements (Xo is independent of time) Xo = Xe + Xv (1-21) 1 Oversimplified Viscoelasticity THE MAXWELL MODEL (p. 2) Thus: Ke(Xo − Xv) = Kv(dXv/dt) with B. C. Xv = 0 at t = 0 (1-22) (Ke/Kv)dt = dXv/(Xo − Xv) Integrate: (Ke/Kv)t = − ln(Xo − Xv) + C Apply B. C.: Xv = 0 at t = 0 means C = ln(Xo) −(Ke/Kv)t = ln[(Xo − Xv)/Xo] (Xo − Xv)/Xo = exp(−Ket/Kv) Thus: F (t) = KeXo exp(−Ket/Kv) (1-23) The force from our constant stretch experiment decays exponentially with time in the Maxwell Model. The relaxation time is λ ≡ Kv/Ke (units s) The force drops to 1/e of its initial value at the relaxation time λ. Initially the force is F (0) = KeXo, the force in the spring, but eventually the force decays to zero F (∞) = 0. 2 Oversimplified Viscoelasticity THE MAXWELL MODEL (p. 3) Constant Area A means stress σ(t) = F (t)/A σ(0) ≡ σ0 = KeXo/A Maxwell Model Stress Relaxation: σ(t) = σ0 exp(−t/λ) Figure 1: Stress Relaxation of a Maxwell Element 3 Oversimplified Viscoelasticity THE MAXWELL MODEL (p.
    [Show full text]
  • Fastran an Advanced Non-Linear Crack-Closure Based Life-Prediction Code
    FASTRAN AN ADVANCED NON-LINEAR CRACK-CLOSURE BASED LIFE-PREDICTION CODE J. C. Newman, Jr. Department of Aerospace Engineering Mississippi State University AFGROW WORKSHOP Layton, Utah September 15, 2015 ffa OUTLINE OF PRESENTATION • Brief History on Fatigue-Crack Growth • Plasticity-Induced Crack-Closure Model • Crack Initiation and Small-Crack Behavior • Fatigue-Crack Growth and Fracture • Concluding Remarks fastran # 2 Stress Concentration Factor for an Elliptical Hole in an Infinite Plate Inglis (1913) c se = S KT 2c c fastran # 3 Notch Strength Analysis – Fracture Mechanics c c Paul Kuhn George Irwin Notch Strength Analysis Fracture Mechanics (Neuber ) (Griffith) fastran # 4 Father of “Modern” Fracture Mechanics Irwin, 1957 George Rankin Irwin (1907-1998) + T fastran # 5 5 25 Notch-Strength Analyses: McEvily and Illg (LaRC), NACA TN-4394, 1958 7075-T6 KNSnet against da/dN fastran # 6 Fracture Mechanics: Paris, Gomez, and Anderson, Trends in Engineering, Seattle, WA, 1961 LEFM: K against d(2a)/dN Paris (1970): KNSnet ~ Kmax fastran # 7 Plasticity-Induced Fatigue-Crack Closure: Elber, 1968 fastran # 8 DOMINANT MECHANISMS OF FATIGUE-CRACK CLOSURE Plastic wake Oxide debris Elber, 1968 Beevers, 1979 Paris et al., 1972 Newman, 1976 Suresh & Ritchie, 1982 Suresh & Ritchie, 1981 (a)(FASTRAN) Plasticity-induced (b) Roughness-induced (c) Oxide/corrosion product- closure closure fastran induced # 9 closure OUTLINE OF PRESENTATION • Brief History on Fatigue-Crack Growth • Plasticity-Induced Crack-Closure Model • Crack Initiation and Small-Crack
    [Show full text]
  • 4.2 Failure Criteria
    4.2 Failure Criteria The determination of residual strength for uncracked structures is straightforward because the ultimate strength of the material is the residual strength. A crack in a structure causes a high stress concentration resulting in a reduced residual strength. When the load on the structure exceeds a certain limit, the crack will extend. The crack extension may become immediately unstable and the crack may propagate in a fast uncontrollable manner causing complete fracture of the component. Figure 4.2.1 illustrates the results obtained from a series of tests conducted on a lug geometry containing a crack. The lug geometry shown in Figure 4.2.1a is a single-load-path structure. Figure 4.2.1b indicates that the cracks in each of the three tests extended abruptly at a critical level of load, which is noted to be a function of a crack length. The crack length-critical load level data shown in Figure 4.2.1b provide the basis for establishing the residual strength capability curve. The locus of critical load levels as a function of crack length is shown in Figure 4.2.1c, where the residual strength capability of the lug structure is shown to decrease with increasing crack length. Figure 4.2.1. Description of Crack Geometry and Residual Strength Results Considering the preceding in terms of applied stress (σ) rather than load gives the σ versus a and σc versus ac plots as shown in Figure 4.2.2 a and b. Schematically, the plots exhibit the same abrupt fracture behavior as the curves presented in Figure 4.2.1.
    [Show full text]
  • A Review on Stress Intensity Factor B.D Bhalekar1, R
    International Research Journal of Engineering and Technology (IRJET) e-ISSN: 2395 -0056 Volume: 03 Issue: 03 | March-2016 www.irjet.net p-ISSN: 2395-0072 A Review on Stress Intensity Factor B.D Bhalekar1, R. B. Patil2 1 Post Graduate student, Mechanical Engineering department, JNEC, Maharashtra, India 2 Associate Professor., Mechanical Engineering department, JNEC, Maharashtra, India ---------------------------------------------------------------------***--------------------------------------------------------------------- Abstract - In this paper effort are taken to review on predictions with real life failures fractography is widely the Stress Analysis of cracked plate for the used with fracture mechanics. The prediction of crack growth is very important in damage tolerance discipline. determination of stress intensity factor near the crack The objective of fracture mechanics is to determine what tip. Cracks in plates with different configurations often conditions will create and drive a crack. Engineers can observed in both modern and classical aerospace, expertly design against this particular mode of failure by mechanical engineering structure. To understand effect understanding the phenomena of fracture. of loading mode and crack configuration on load bearing capacity of such plate is very significant in 2. FRACTURE FAILURE designing of structure. There are number of Analytical, Consider a cracked plate on which a force is applied which Numerical, and experimental technique available for might enable the crack to propagate. Irwin projected a stress analysis to determine the stress intensity factor classification matching to the three situations represented near crack tip under different loading modes in an in Fig .1. Thus, three distinct modes are considered: mode infinite and finite cracked plate made up of different I, mode II and mode III.
    [Show full text]
  • Elasticity and Viscoelasticity
    CHAPTER 2 Elasticity and Viscoelasticity CHAPTER 2.1 Introduction to Elasticity and Viscoelasticity JEAN LEMAITRE Universite! Paris 6, LMT-Cachan, 61 avenue du President! Wilson, 94235 Cachan Cedex, France For all solid materials there is a domain in stress space in which strains are reversible due to small relative movements of atoms. For many materials like metals, ceramics, concrete, wood and polymers, in a small range of strains, the hypotheses of isotropy and linearity are good enough for many engineering purposes. Then the classical Hooke’s law of elasticity applies. It can be de- rived from a quadratic form of the state potential, depending on two parameters characteristics of each material: the Young’s modulus E and the Poisson’s ratio n. 1 c * ¼ A s s ð1Þ 2r ijklðE;nÞ ij kl @c * 1 þ n n eij ¼ r ¼ sij À skkdij ð2Þ @sij E E Eandn are identified from tensile tests either in statics or dynamics. A great deal of accuracy is needed in the measurement of the longitudinal and transverse strains (de Æ10À6 in absolute value). When structural calculations are performed under the approximation of plane stress (thin sheets) or plane strain (thick sheets), it is convenient to write these conditions in the constitutive equation. Plane stress ðs33 ¼ s13 ¼ s23 ¼ 0Þ: 2 3 1 n 6 À 0 7 2 3 6 E E 72 3 6 7 e11 6 7 s11 6 7 6 1 76 7 4 e22 5 ¼ 6 0 74 s22 5 ð3Þ 6 E 7 6 7 e12 4 5 s12 1 þ n Sym E Handbook of Materials Behavior Models Copyright # 2001 by Academic Press.
    [Show full text]
  • Viscoelastic Creep Crack Growth: a Review of Fracture Mechanical Analyses
    Mechanics of Time-Dependent Materials 1: 241–268, 1998. 241 c 1998 Kluwer Academic Publishers. Printed in the Netherlands. Viscoelastic Creep Crack Growth: A Review of Fracture Mechanical Analyses W. BRADLEY1, W.J. CANTWELL2 and H.H. KAUSCH2 1Department of Mechanical Engineering, Texas A&M University, College Station, TX, U.S.A.; 2Materials Science Department, EPFL, DMX-LP, Swiss Federal Institute of Technology, 1015 Lausanne, Switzerland (Received 8 February 1997; accepted in revised form 2 October 1997) Abstract. The study of time dependent crack growth in polymers using a fracture mechanics approach has been reviewed. The time dependence of crack growth in polymers is seen to be the result of the viscoelastic deformation in the process zone, which causes the supply of energy to drive the crack to occur over time rather than instantaneously, as it does in metals. Additional time dependence in the crack growth process can be due to process zone behavior, where both the flow stress and the critical crack tip opening displacement may be dependent on the crack growth rate. Instability leading to slip-stick crack growth has been seen to be the consequence of a decrease in the critical energy release rate with increasing crack growth rate due to adiabatic heating causing a reduction in the process zone flow stress, a decrease in the crack tip opening displacement due to a ductile to brittle transition at higher crack growth rates, or an increase in the rate of fracture work due to more rapid viscoelastic deformation. Finally, various techniques to experimentally characterize the crack growth rate as a function of stress intensity have been critiqued.
    [Show full text]
  • Viscous-To-Viscoelastic Transition in Phononic Crystal and Metamaterial Band Structures
    Viscous-to-viscoelastic transition in phononic crystal and metamaterial band structures Michael J. Frazier and Mahmoud I. Husseina) Department of Aerospace Engineering Sciences, University of Colorado Boulder, Boulder, Colorado 80309, USA (Received 3 March 2015; revised 8 September 2015; accepted 16 October 2015; published online 19 November 2015) The dispersive behavior of phononic crystals and locally resonant metamaterials is influenced by the type and degree of damping in the unit cell. Dissipation arising from viscoelastic damping is influenced by the past history of motion because the elastic component of the damping mechanism adds a storage capacity. Following a state-space framework, a Bloch eigenvalue problem incorpo- rating general viscoelastic damping based on the Zener model is constructed. In this approach, the conventional Kelvin–Voigt viscous-damping model is recovered as a special case. In a continuous fashion, the influence of the elastic component of the damping mechanism on the band structure of both a phononic crystal and a metamaterial is examined. While viscous damping generally narrows a band gap, the hereditary nature of the viscoelastic conditions reverses this behavior. In the limit of vanishing heredity, the transition between the two regimes is analyzed. The presented theory also allows increases in modal dissipation enhancement (metadamping) to be quantified as the type of damping transitions from viscoelastic to viscous. In conclusion, it is shown that engineering the dissipation allows one to control the dispersion (large versus small band gaps) and, conversely, engineering the dispersion affects the degree of dissipation (high or low metadamping). VC 2015 Acoustical Society of America.[http://dx.doi.org/10.1121/1.4934845] [ANN] Pages: 3169–3180 I.
    [Show full text]