SUPPLEMENTARY INFORMATION
|
|
- Michael Heath
- 5 years ago
- Views:
Transcription
1 Supplementary material How things break under shear and tension? How does a crack propagate in three dimensions when a material is both under tension and sheared parallel to the crack front? Numerous experimental observations for over half a century have shown that mixing those two fracture modes (antiplane shear and tension) produces stepped fracture surfaces with characteristic lance-shaped markings. Those markings are ubiquitous in both engineering and geological materials including glasses, polymers, ceramics, metals, and rocks. Although it is known that stepped surfaces result from a complex segmentation of the crack front into partial fronts, the mechanism of this segmentation has remained elusive. Besides these crack-front evolutions are hard to track during the material failing, especially in three dimensions. Movies of three-dimensional crack-front instability The accompanying movies of phase-field simulations reveal for the first time the complex path followed by a crack-front in mixed mode fracture. Those movies were obtained by numerical simulations of crack propagation using the phase-field model described in this supplementary material. The first two movies show frontal (frontal.mov) and lateral (lateral.mov) views of the propagating crack-front that produced the fracture surfaces that are displayed at different times in Fig. 1(c-d) of the article. The third movie (helical.mov) is an animation of Fig. 2(a) that highlights the helical nature of the crack-front instability. It shows that the front deforms into a helix that subsequently evolves nonlinearly into a more complex sawtooth wave shape as the crack-front segments into partial fronts to produce the observed stepped surfaces. 1
2 Phase-field method Crack propagation was studied using the phase-field method 1-4. This method has the advantage that it tracks automatically the crack-front evolution by the introduction of a scalar phase field, which distinguishes between broken and unbroken states of the material. This field varies smoothly in space on the fracture process zone scale, thereby providing a smooth cut-off for the divergence of the singular stress fields. This method has recently been tested numerically in two dimensions through benchmark comparisons with predictions of continuum fracture mechanics for kink cracks under mixed mode I+II plane strain loading 1 and for wavy cracks under biaxial loading 2. In addition, a crack propagation law consistent with the standard principle of local symmetry 5 was derived analytically from the phase-field model for twodimensional propagation in isotropic media 2. The total energy of the material is represented by the functional E = d 3 x ρ 2 tu φ 2 + g(φ) ( e strain e c ) (s1) ( ) 2 + κ 2 where is the strain energy density and is the usual strain tensor of linear elasticity with corresponding to, respectively. The broken state of the material becomes energetically favoured when the strain energy exceeds the threshold and is a monotonously increasing function that describes the softening of the elastic energy at large strain. The evolution equations for and the three components of the displacement are derived variationally from equation (s1) and are given by (s2) (s3) 2
3 The resulting set of four coupled partial differential equations provides a self-consistent description of both macroscopic linear elasticity and material failure without the need to track explicitly the crack front location, which can be defined conveniently as the leading edge of the fracture surface. Moreover, this gradient dynamics guarantees that the total energy is a monotonously decreasing function of time with all the energy dissipated inside the fracture process zone on a characteristic time scale. After writing the phase-field equations in dimensionless form by measuring length in units of the fracture process zone scale and time in units of, the remaining parameters were chosen as, ( ), and, where is the shear wave speed. Slab geometry and loading conditions. Crack propagation was simulated in a large strip of width, thickness, and length along the propagation axis. Mixed mode I+III loading was imposed by choosing initially (mode I) and (mode III) corresponding to stress intensity factors and displacements fixed at all time on the, respectively, and by keeping the boundaries. Periodic boundary conditions were imposed in the direction. The strip length and width ( and ) were chosen large enough to eliminate boundary effects as much as possible and was varied over a wide range (see Figure captions) to study the development of instability for different extensions of the crack front. Aside from the strip size, the control parameters of the simulations include the ratio and the ratio where is the energy release rate at planar crack propagation and the Griffith threshold for propagation of a semi-infinite crack. The phase-field model describes an ideally brittle limit where this threshold is simply twice the surface energy,, where the surface energy is given by the expression 1,4 (s4) 3
4 For the given choice of, the above expression yields. The ratio was chosen in the range 1.25 to In this range, the speed of the parent and daughter cracks do not exceed about one third of the shear wave speed and inertial effects do not suffice to produce a crack branching instability. Therefore the crack segmentation instability studied here only occurs with the superposition of mode III for this range of. This was checked explicitly by repeating a few test simulations with the only modification that the quasistatic equations without inertia were solved with a successive over-relaxation (SOR) method instead of the wave equations. The main characteristics of the instability were unchanged. Numerical implementation. The equations were discretized with a second-order accurate finite difference scheme on a uniform mesh with grid spacing (which yields a slightly larger surface energy ) and integrated with an explicit scheme that handles accurately the energy equation with a time step. The FORTRAN code was parallelized with MPI to handle the computing and memory requirements of large strip sizes with up to grid points and the simulations were carried out on Linux Clusters at Northeastern University. Relationship of phase-field and materials parameters. The parameters that enter the phase-field model include both known material parameters and parameters specific to the phase-field model that can be related to physical quantities. The known material parameters include the density and the Lamé coefficients in the elastic energy, from which we can define Poisson s ratio and the shear wave speed. The three parameters specific to the phase-field model include the coefficient of the phase-field gradient-square term in the energy functional (equation (s1)), the threshold elastic energy density for bond breaking, and the kinetic coefficient that governs the rate of evolution of the phase-field in equation (s2). As discussed above, those last three parameters can be combined to define (i) the physical 4
5 length that measures the size of the process zone, which is the region around the crack edge where the phase-field increases from a value close to zero in the broken material to one in the unbroken material, (ii) the characteristic time scale of energy dissipation in this zone, and (iii) the surface energy via equation (s4). With length and time measured in units of and, respectively, the crack dynamics in the phase-field model is controlled solely by Poisson s ratio, the dimensionless combination, and the external loading conditions through the ratios and, where is the energy release rate. In the quasi-static limit where inertial effects are small, the geometrical evolution of the crack surface becomes insensitive to the wave speed and thus to the ratio. Consequently, this evolution is controlled predominantly by Poisson s ratio and external loading conditions with the process zone scale acting only as a scaling length. Namely, if a material A has a process zone scale a hundred times larger than some other material B, the initial scale of instability would be predicted by equation (2) to be a hundred times larger in A than B, but the crack evolutions would be identical in the two materials up to a change of scale. Thus estimates of, which can vary by several orders of magnitude for different materials, can be used to compare phase-field simulations to experiments. This was done in the main text when comparing phase-field predictions for the initial instability wavelength (equation 2) to experimental observations in glass and PMMA. Instability wavelength An analytical expression for the instability wavelength can be derived heuristically by considering the effects of both configurational and cohesive forces acting on a helical crack front. Configurational forces acting perpendicular to the local crack plane have been treated by Eshelby 6 in the traditional theoretical framework of continuum fracture mechanics. The configuration force originates from the directional dependence of the 5
6 energy release rate that is present when is non-vanishing at some point along the Figure s1. Schematic representation of the in-plane (green line) and the outof-plane (red line) projections of a helical crack-front. The arrows indicate the direction of forces acting perpendicularly to the plane of the parent crack at points (black filled circles) separating leading A and lagging B zones of the front. The forces include the destabilizing configurational forces originating from the directional dependence of the energy release rate (orange arrows) and the stabilizing cohesive forces (gray arrows). front. In this case, the stress distribution is non-symmetrical about the instantaneous crack propagation axis and the crack can release energy at a faster rate by changing its growth direction. This force has magnitude 1,6, where the prime symbols indicate that the stress intensity factors are evaluated locally along the front and are distinct from the unprimed ones that are fixed by external loading conditions. More recently, Hakim and Karma derived laws of crack motion from the phase-field model of fracture in two dimensions. They obtained a condition for at the crack tip that reflects the balance between Eshelby s configurational force 6 and cohesive forces acting perpendicularly to the crack front. The latter are non-vanishing when the surface energy is anisotropic. When the surface energy is isotropic, which is the case considered here, 6
7 this force balance reduces to the standard principal of local symmetry 5 ( ). For propagation in three dimensions, however, the crack front is a curved line. In this case, Eshelby s configurational force 6 perpendicular to the local crack plane can be balanced by a local cohesive force generated by the out-of-plane curvature of the front. This cohesive force has a magnitude where is the surface energy and is the local out-of-plane radius of curvature of the front as illustrated schematically in Figure s1 of this supplemental material. The balance of these forces yields the condition of equation (1) in the main text, which extends the principle of local symmetry ( three dimensions. A scaling relation for the marginal instability wavelength ) to, i.e. the wavelength for which a helical perturbation neither grows nor decays, can now be derived by requiring that those forces balance on the sides of protruding A facets, which correspond to the points separating A and B zones marked by filled circles in Figure s1. At those points, we have and, where is the amplitude of a helical perturbation of wavenumber. Furthermore, we use the fact that for a small amplitude perturbation close to the onset of propagation and that in the limit cohesive forces yields the equality. Balancing the configurational and. Substituting in this equality the above relations,,, and, we obtain the prediction for the marginally unstable wavelength. Finally, we use the fact that the numerically computed linear stability spectrum (Fig. 2(c) of the main text) is well fitted by a quadratic polynomial of the form. This form implies that the fastest growing wavelength is approximately twice larger than the marginally stable wavelength, yielding. 7
8 Facet coarsening The coarsening of the wavelength by the elimination of daughter cracks is directly analogous to the coarsening of the wavelength of finger fronts in other well-known interfacial pattern forming instabilities, including Saffman-Taylor viscous fingers 7, directional solidification fingers 8 and Laplacian needle growth 9. In those examples, the fastest unstable mode of an array of fingers typically has a wavelength equal to twice the finger spacing 7,8. The symmetry is spontaneously broken by small perturbations and the spatial period doubling instability has generically no threshold in situations where the system is invariant under translation along the growth direction. The absence of threshold stems from this translational invariance (i.e. the standard goldstone mode argument). With an imposed gradient as in directional solidification 8, translation symmetry is broken and coarsening stops when the finger spacing is large enough. The rule is that the number of fingers decreases roughly inversely proportionally to the total growth length of the fingers up to small logarithmic corrections 9. This inverse relation follows from the simple picture that fingers interact through a long range field (here the stress field) without any intrinsic length scale. Therefore, there are no other scales for the coarsening than the spacing between actively growing fingers. This rule appears to describe well the coarsening of daughter cracks in Figure 3. Facet rotation angle Analytical expressions for the facet rotation angle can be obtained using the expressions for the stress intensity factors at the edge of daughter cracks
9 The assumptions that facets orient to be free of shear stress (shear-free), is mathematically equivalent to the assumption that they orient to maximize the mode I stress intensity factor,. Both principles yield at once the prediction of equation (3) in the main text. The assumption that the facets orient to maximize the energy release rate yields the condition where. This condition gives the same prediction as the first two principles up to a threshold value of beyond which it bifurcates into low and high angle branches with equal maximum energy release rate as shown in Fig. 4(a) in the main text. Above this threshold, which depends on Poisson s ratio, the shear- free branch has minimum energy release rate. As shown in Figure 4(a) of the manuscript, the shear-free (SF) and maximum energyrelease-rate (maximum G') conditions give identical predictions up to a critical value of /, beyond which the latter predicts a bifurcation into separate low and high angle branches (dashed lines). Experiments to date have explored ratios of / that fall below the threshold value of this bifurcation. Hence, they have been so far unable to test whether such a bifurcation occurs. The present phase-field simulations explore a much larger range of / and show that this bifurcation is absent. The phase field facet angle increases smoothly with /, in qualitative agreement with the shearfree prediction. After reaching a maximum, it then decreases at large. This decrease, which is not predicted by previous theories, can be understood by noting that 9
10 the magnitude of the destabilizing configurational force is proportional to K I K III because along a helical crack front is induced by mode III. Since vanishes at large for fixed energy release rate, this destabilizing force also vanishes for large. Thus cohesive forces that stabilize planar crack growth ultimately dominate for large, thereby causing the facet angle to decrease. We note that Sommer s experiments in glass 11 only explored a range of / much smaller than unity. Consequently, those experiments produced very small rotation angles (less than 3 degrees) that are not shown in Fig. 4(b) of the manuscript since they do not provide an adequate basis of comparison with our phase-field results. Only experiments in steel and PMMA that explore a larger range of /, and hence a larger range of angles, are shown. Also, to make the comparison of phase-field results and experiments quantitative, we have plotted the ratio versus / to scale out the dependence on Poisson s ratio that is material-dependent. REFERENCES 1. Hakim, V. & Karma, A., Laws of crack motion and phase-field models of fracture. J. Mech. Phys. Solids 57, (2009). 2. Henry, H. & Levine, H., Dynamic instabilities of fracture under biaxial strain using a phase field model. Phys. Rev. Lett. 93, (2004). 3. Karma, A., in Handbook of Materials Modeling (ed. Yip, S.) (Springer, Netherlands, 2005). 4. Karma, A., Kessler, D. A., & Levine, H., Phase-field model of mode III dynamic fracture. Phys. Rev. Lett. 87, (2001). 5. Goldstein, R. V. & Salganik, R. L., Brittle fracture of solids with arbitrary cracks. Int. J. Fract. 10, (1974). 6. Eshelby, J.D., The elastic energy-momentum tensor. J. Elast. 5, (1975). 7. Kessler, D. A. & Levine, H.,Coalescence of Saffman-Taylor fingers: A new global instability. Phys. Rev. A 33, (1986). 10
11 8. Losert, W., Shi, B. Q. & Cummins, H. Z., Spatial period-doubling instability of dendritic arrays in directional solidification. Phys. Rev. Lett. 77, (1996). 9. Krug, J., Kassner, K., Meakin, P. & Family, F., Laplacian needle growth. Europhys. Lett. 24, (1993). 10. Cooke, M. L. & Pollard, D. D., Fracture propagation paths under mixed mode loading within rectangular blocks of polymethyl methacrylate. J. Geophys. Res. 101, (1996). 11. Sommer, E. Formation of fracture `lances' in glass. Eng. Fract. Mech. 1, (1969). 11
Linear Elastic Fracture Mechanics
Measure what is measurable, and make measurable what is not so. - Galileo GALILEI Linear Elastic Fracture Mechanics Krishnaswamy Ravi-Chandar Lecture presented at the University of Pierre and Marie Curie
More informationConfigurational Forces as Basic Concepts of Continuum Physics
Morton E. Gurtin Configurational Forces as Basic Concepts of Continuum Physics Springer Contents 1. Introduction 1 a. Background 1 b. Variational definition of configurational forces 2 с Interfacial energy.
More informationFracture mechanics fundamentals. Stress at a notch Stress at a crack Stress intensity factors Fracture mechanics based design
Fracture mechanics fundamentals Stress at a notch Stress at a crack Stress intensity factors Fracture mechanics based design Failure modes Failure can occur in a number of modes: - plastic deformation
More informationThe effect of plasticity in crumpling of thin sheets: Supplementary Information
The effect of plasticity in crumpling of thin sheets: Supplementary Information T. Tallinen, J. A. Åström and J. Timonen Video S1. The video shows crumpling of an elastic sheet with a width to thickness
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Surface and body forces Tensors, Mohr circles. Theoretical strength of materials Defects Stress concentrations Griffith failure
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lectures & 3, 9/31 Aug 017 www.geosc.psu.edu/courses/geosc508 Discussion of Handin, JGR, 1969 and Chapter 1 Scholz, 00. Stress analysis and Mohr Circles Coulomb Failure
More informationBifurcation Analysis in Geomechanics
Bifurcation Analysis in Geomechanics I. VARDOULAKIS Department of Engineering Science National Technical University of Athens Greece and J. SULEM Centre d'enseignement et de Recherche en Mecanique des
More informationEffect of Anisotropic Permeability on Band Growth
Effect of Anisotropic Permeability on Band Growth Jesse Taylor-West October 7, 24 Introduction The mechanics of partially molten regions of the mantle are not well understood- in part due to the inaccessibility
More informationStudies of Bimaterial Interface Fracture with Peridynamics Fang Wang 1, Lisheng Liu 2, *, Qiwen Liu 1, Zhenyu Zhang 1, Lin Su 1 & Dan Xue 1
International Power, Electronics and Materials Engineering Conference (IPEMEC 2015) Studies of Bimaterial Interface Fracture with Peridynamics Fang Wang 1, Lisheng Liu 2, *, Qiwen Liu 1, Zhenyu Zhang 1,
More informationFig. 1. Different locus of failure and crack trajectories observed in mode I testing of adhesively bonded double cantilever beam (DCB) specimens.
a). Cohesive Failure b). Interfacial Failure c). Oscillatory Failure d). Alternating Failure Fig. 1. Different locus of failure and crack trajectories observed in mode I testing of adhesively bonded double
More informationModeling the Dynamic Propagation of Shear Bands in Bulk Metallic Glasses
Modeling the Dynamic Propagation of Shear Bands in Bulk Metallic Glasses B.J. Edwards, K. Feigl, M.L. Morrison*, B. Yang*, P.K. Liaw*, and R.A. Buchanan* Dept. of Chemical Engineering, The University of
More informationBrittle fracture of rock
1 Brittle fracture of rock Under the low temperature and pressure conditions of Earth s upper lithosphere, silicate rock responds to large strains by brittle fracture. The mechanism of brittle behavior
More informationarxiv:cond-mat/ v1 25 Feb 1994
A Model for Fracture in Fibrous Materials A. T. Bernardes Departamento de Física - ICEB arxiv:cond-mat/9402110v1 25 Feb 1994 Universidade Federal de Ouro Preto Campus do Morro do Cruzeiro 35410-000 Ouro
More informationFinite Element simulations of a phase-field model for mode-iii fracture. by V. Kaushik
Finite Element simulations of a phase-field model for mode-iii fracture. by V. Kaushik 1 Motivation The motivation behind this model is to understand the underlying physics behind branching of cracks during
More informationVariational phase field model for dynamic brittle fracture
Variational phase field model for dynamic brittle fracture Bleyer J., Roux-Langlois C., Molinari J-F. EMMC 15, September 8th, 2016 1 / 18 Outline Mechanisms of dynamic fracture Variational phase-field
More informationCracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics
Cracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics Tobias Gasch, PhD Student Co-author: Prof. Anders Ansell Comsol Conference 2016 Munich 2016-10-12 Contents Introduction Isotropic damage
More informationCrack dynamics in elastic media
PHILOSOPHICAL MAGAZINE B, 1998, VOL. 78, NO. 2, 97± 102 Crack dynamics in elastic media By Mokhtar Adda-Bedia and Martine Ben Amar Laboratoire de Physique Statistique de l Ecole Normale Supe  rieure,
More informationMechanical properties of polymers: an overview. Suryasarathi Bose Dept. of Materials Engineering, IISc, Bangalore
Mechanical properties of polymers: an overview Suryasarathi Bose Dept. of Materials Engineering, IISc, Bangalore UGC-NRCM Summer School on Mechanical Property Characterization- June 2012 Overview of polymer
More informationSlow crack growth in polycarbonate films
EUROPHYSICS LETTERS 5 July 5 Europhys. Lett., 7 (), pp. 4 48 (5) DOI:.9/epl/i5-77-3 Slow crack growth in polycarbonate films P. P. Cortet, S. Santucci, L. Vanel and S. Ciliberto Laboratoire de Physique,
More informationMechanics of wafer bonding: Effect of clamping
JOURNAL OF APPLIED PHYSICS VOLUME 95, NUMBER 1 1 JANUARY 2004 Mechanics of wafer bonding: Effect of clamping K. T. Turner a) Massachusetts Institute of Technology, Cambridge, Massachusetts 0219 M. D. Thouless
More informationSUPPLEMENTARY INFORMATION 1
1 Supplementary information Effect of the viscoelasticity of substrate: In the main text, we indicated the role of the viscoelasticity of substrate. In all problems involving a coupling of a viscous medium
More informationA Finite Difference Implementation of Phase Field Theory
University of Connecticut DigitalCommons@UConn Master's Theses University of Connecticut Graduate School 12-16-2012 A Finite Difference Implementation of Phase Field Theory Nicholas W. Oren University
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationCritical applied stresses for a crack initiation from a sharp V-notch
Focussed on: Fracture and Structural Integrity related Issues Critical applied stresses for a crack initiation from a sharp V-notch L. Náhlík, P. Hutař Institute of Physics of Materials, Academy of Sciences
More informationarxiv: v2 [cond-mat.mtrl-sci] 25 Jul 2015
Recent developments in dynamic fracture: Some perspectives Jay Fineberg 1 and Eran Bouchbinder 2 1 Racah Institute of Physics, Hebrew University of Jerusalem, Jerusalem 91904, Israel 2 Chemical Physics
More informationBrittle fracture dynamics with arbitrary paths III. The branching instability under general loading
Journal of the Mechanics and Physics of Solids 53 (2005) 227 248 www.elsevier.com/locate/jmps Brittle fracture dynamics with arbitrary paths III. The branching instability under general loading M. Adda-Bedia
More informationANSYS Mechanical Basic Structural Nonlinearities
Lecture 4 Rate Independent Plasticity ANSYS Mechanical Basic Structural Nonlinearities 1 Chapter Overview The following will be covered in this Chapter: A. Background Elasticity/Plasticity B. Yield Criteria
More informationCHAPTER 1 MODELING DYNAMIC FRACTURE USING LARGE-SCALE ATOMISTIC SIMULATIONS
CHAPTER 1 MODELING DYNAMIC FRACTURE USING LARGE-SCALE ATOMISTIC SIMULATIONS Markus J. Buehler Massachusetts Institute of Technology, Department of Civil and Environmental Engineering 77 Massachusetts Avenue
More informationConservation of mass. Continuum Mechanics. Conservation of Momentum. Cauchy s Fundamental Postulate. # f body
Continuum Mechanics We ll stick with the Lagrangian viewpoint for now Let s look at a deformable object World space: points x in the object as we see it Object space (or rest pose): points p in some reference
More informationLinearized Theory: Sound Waves
Linearized Theory: Sound Waves In the linearized limit, Λ iα becomes δ iα, and the distinction between the reference and target spaces effectively vanishes. K ij (q): Rigidity matrix Note c L = c T in
More informationLecture #7: Basic Notions of Fracture Mechanics Ductile Fracture
Lecture #7: Basic Notions of Fracture Mechanics Ductile Fracture by Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling of Materials in Manufacturing
More informationFor an imposed stress history consisting of a rapidly applied step-function jump in
Problem 2 (20 points) MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0239 2.002 MECHANICS AND MATERIALS II SOLUTION for QUIZ NO. October 5, 2003 For
More informationHomework Problems. ( σ 11 + σ 22 ) 2. cos (θ /2), ( σ θθ σ rr ) 2. ( σ 22 σ 11 ) 2
Engineering Sciences 47: Fracture Mechanics J. R. Rice, 1991 Homework Problems 1) Assuming that the stress field near a crack tip in a linear elastic solid is singular in the form σ ij = rλ Σ ij (θ), it
More informationCrack Front Waves in Dynamic Fracture
International Journal of Fracture 119: 247 261, 2003. 2003 Kluwer Academic Publishers. Printed in the Netherlands. Crack Front Waves in Dynamic Fracture J. FINEBERG 1,, E. SHARON 2 and G. COHEN 1 1 The
More informationChapter 7. Highlights:
Chapter 7 Highlights: 1. Understand the basic concepts of engineering stress and strain, yield strength, tensile strength, Young's(elastic) modulus, ductility, toughness, resilience, true stress and true
More informationClassical fracture and failure hypotheses
: Chapter 2 Classical fracture and failure hypotheses In this chapter, a brief outline on classical fracture and failure hypotheses for materials under static loading will be given. The word classical
More informationarxiv: v1 [cond-mat.mtrl-sci] 1 Nov 2009
arxiv:0911.0173v1 [cond-mat.mtrl-sci] 1 Nov 2009 Submitted to Annual Review of Condensed Matter Physics Dynamics of Simple Cracks Eran Bouchbinder 1, Jay Fineberg 2, and M. Marder 3 1 Department of Chemical
More informationPhysics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top
Physics 106a, Caltech 4 December, 2018 Lecture 18: Examples on Rigid Body Dynamics I go through a number of examples illustrating the methods of solving rigid body dynamics. In most cases, the problem
More informationA TIME-DEPENDENT DAMAGE LAW IN DEFORMABLE SOLID: A HOMOGENIZATION APPROACH
9th HSTAM International Congress on Mechanics Limassol, Cyprus, - July, A TIME-DEPENDENT DAMAGE LAW IN DEFORMABLE SOLID: A HOMOGENIZATION APPROACH Cristian Dascalu, Bertrand François, Laboratoire Sols
More informationstrain appears only after the stress has reached a certain critical level, usually specied by a Rankine-type criterion in terms of the maximum princip
Nonlocal damage models: Practical aspects and open issues Milan Jirasek LSC-DGC, Swiss Federal Institute of Technology at Lausanne (EPFL), Switzerland Milan.Jirasek@ep.ch Abstract: The purpose of this
More informationNumerical modeling of standard rock mechanics laboratory tests using a finite/discrete element approach
Numerical modeling of standard rock mechanics laboratory tests using a finite/discrete element approach S. Stefanizzi GEODATA SpA, Turin, Italy G. Barla Department of Structural and Geotechnical Engineering,
More informationMICRO-BRANCHING AS AN INSTABILITY IN DYNAMIC FRACTURE. To appear in the proceedings of the IUTAM conference (Cambridge, 1995)
To appear in the proceedings of the IUTAM conference (Cambridge, 1995) MICRO-BRANCHING AS AN INSTABILITY IN DYNAMIC FRACTURE J. FINEBERG, S. P. GROSS AND E. SHARON The Racah Institute of Physics The Hebrew
More informationBrittle Deformation. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm
Lecture 6 Brittle Deformation Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm WW Norton, unless noted otherwise Brittle deformation EarthStructure (2 nd
More informationThe effective slip length and vortex formation in laminar flow over a rough surface
The effective slip length and vortex formation in laminar flow over a rough surface Anoosheh Niavarani and Nikolai V. Priezjev Movies and preprints @ http://www.egr.msu.edu/~niavaran A. Niavarani and N.V.
More information1. Comparison of stability analysis to previous work
. Comparison of stability analysis to previous work The stability problem (6.4) can be understood in the context of previous work. Benjamin (957) and Yih (963) have studied the stability of fluid flowing
More informationThe waves of damage in elastic plastic lattices with waiting links: Design and simulation
Mechanics of Materials 3 () 7 75 www.elsevier.com/locate/mechmat The waves of damage in elastic plastic lattices with waiting links: Design and simulation A. Cherkaev *, V. Vinogradov, S. Leelavanichkul
More informationPeridynamic model for dynamic fracture in unidirectional fiber-reinforced composites
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Mechanical & Materials Engineering Faculty Publications Mechanical & Materials Engineering, Department of 4-2012 Peridynamic
More informationCohesive Band Model: a triaxiality-dependent cohesive model inside an implicit non-local damage to crack transition framework
University of Liège Aerospace & Mechanical Engineering MS3: Abstract 131573 - CFRAC2017 Cohesive Band Model: a triaxiality-dependent cohesive model inside an implicit non-local damage to crack transition
More informationULTRASONIC REFLECTION BY A PLANAR DISTRIBUTION OF SURFACE BREAKING CRACKS
ULTRASONIC REFLECTION BY A PLANAR DISTRIBUTION OF SURFACE BREAKING CRACKS A. S. Cheng Center for QEFP, Northwestern University Evanston, IL 60208-3020 INTRODUCTION A number of researchers have demonstrated
More informationInstabilities and Dynamic Rupture in a Frictional Interface
Instabilities and Dynamic Rupture in a Frictional Interface Laurent BAILLET LGIT (Laboratoire de Géophysique Interne et Tectonophysique) Grenoble France laurent.baillet@ujf-grenoble.fr http://www-lgit.obs.ujf-grenoble.fr/users/lbaillet/
More informationarxiv:cond-mat/ v1 [cond-mat.soft] 29 May 2002
Stretching Instability of Helical Springs David A. Kessler and Yitzhak Rabin Dept. of Physics, Bar-Ilan University, Ramat-Gan, Israel (Dated: October 31, 18) arxiv:cond-mat/05612v1 [cond-mat.soft] 29 May
More informationMMJ1133 FATIGUE AND FRACTURE MECHANICS A - INTRODUCTION INTRODUCTION
A - INTRODUCTION INTRODUCTION M.N.Tamin, CSMLab, UTM Course Content: A - INTRODUCTION Mechanical failure modes; Review of load and stress analysis equilibrium equations, complex stresses, stress transformation,
More informationChapter 6: Mechanical Properties of Metals. Dr. Feras Fraige
Chapter 6: Mechanical Properties of Metals Dr. Feras Fraige Stress and Strain Tension Compression Shear Torsion Elastic deformation Plastic Deformation Yield Strength Tensile Strength Ductility Toughness
More informationNumerical modelling of induced tensile stresses in rock in response to impact loading
Numerical modelling of induced tensile stresses in rock in response to impact loading M.T. Mnisi, D.P. Roberts and J.S. Kuijpers Council for Scientific and Industrial Research (CSIR): Natural Resources
More informationNumerical simulation of delamination onset and growth in laminated composites
Numerical simulation of delamination onset and growth in laminated composites G. Wimmer, C. Schuecker, H.E. Pettermann Austrian Aeronautics Research (AAR) / Network for Materials and Engineering at the
More informationAn Experimental Characterization of the Non-linear Rheology of Rock
An Experimental Characterization of the Non-linear Rheology of Rock G. N. BorrNoTr New England Research Inc. Contract: F49620-95-C-0019 Sponsor: AFOSR ABSTRACT A laboratory experimental program is underway
More informationHyperelasticity governs dynamic fracture at a critical length scale
Hyperelasticity governs dynamic fracture at a critical length scale articles Markus J. Buehler 1 *, Farid F. Abraham 2 * & Huajian Gao 1 * 1 Max Planck Institute for Metals Research, Heisenbergstrasse
More informationModelling Localisation and Spatial Scaling of Constitutive Behaviour: a Kinematically Enriched Continuum Approach
Modelling Localisation and Spatial Scaling of Constitutive Behaviour: a Kinematically Enriched Continuum Approach Giang Dinh Nguyen, Chi Thanh Nguyen, Vinh Phu Nguyen School of Civil, Environmental and
More informationDirect Comparison of Anisotropic Damage Mechanics to Fracture Mechanics of Explicit Cracks
Direct Comparison of Anisotropic Damage Mechanics to Fracture Mechanics of Explicit Cracks John A. Nairn Wood Science and Engineering, Oregon State University, Corvallis, OR 97330, USA Tel: +1-541-737-4265
More information36. TURBULENCE. Patriotism is the last refuge of a scoundrel. - Samuel Johnson
36. TURBULENCE Patriotism is the last refuge of a scoundrel. - Samuel Johnson Suppose you set up an experiment in which you can control all the mean parameters. An example might be steady flow through
More informationMechanical Properties of Polymer Rubber Materials Based on a New Constitutive Model
Mechanical Properties of Polymer Rubber Materials Based on a New Constitutive Model Mechanical Properties of Polymer Rubber Materials Based on a New Constitutive Model J.B. Sang*, L.F. Sun, S.F. Xing,
More informationFinite Element Solution of Nonlinear Transient Rock Damage with Application in Geomechanics of Oil and Gas Reservoirs
Finite Element Solution of Nonlinear Transient Rock Damage with Application in Geomechanics of Oil and Gas Reservoirs S. Enayatpour*1, T. Patzek2 1,2 The University of Texas at Austin *Corresponding author:
More informationA METHOD OF LOAD INCREMENTS FOR THE DETERMINATION OF SECOND-ORDER LIMIT LOAD AND COLLAPSE SAFETY OF REINFORCED CONCRETE FRAMED STRUCTURES
A METHOD OF LOAD INCREMENTS FOR THE DETERMINATION OF SECOND-ORDER LIMIT LOAD AND COLLAPSE SAFETY OF REINFORCED CONCRETE FRAMED STRUCTURES Konuralp Girgin (Ph.D. Thesis, Institute of Science and Technology,
More informationChapter 4. Gravity Waves in Shear. 4.1 Non-rotating shear flow
Chapter 4 Gravity Waves in Shear 4.1 Non-rotating shear flow We now study the special case of gravity waves in a non-rotating, sheared environment. Rotation introduces additional complexities in the already
More informationElements of Rock Mechanics
Elements of Rock Mechanics Stress and strain Creep Constitutive equation Hooke's law Empirical relations Effects of porosity and fluids Anelasticity and viscoelasticity Reading: Shearer, 3 Stress Consider
More informationSize effect in the strength of concrete structures
Sādhanā Vol. 27 Part 4 August 2002 pp. 449 459. Printed in India Size effect in the strength of concrete structures B L KARIHALOO and Q Z XIAO Division of Civil Engineering School of Engineering Cardiff
More informationINFLUENCE OF THE LOCATION AND CRACK ANGLE ON THE MAGNITUDE OF STRESS INTENSITY FACTORS MODE I AND II UNDER UNIAXIAL TENSION STRESSES
INFLUENCE OF THE LOCATION AND CRACK ANGLE ON THE MAGNITUDE OF STRESS INTENSITY FACTORS MODE I AND II UNDER UNIAXIAL TENSION STRESSES Najah Rustum Mohsin Southern Technical University, Technical Institute-Nasiriya,
More informationC. R. McKee and M. E. Hanson Lawrence Livermore Laboratory University of California Livermore, California 94550
PREDICTING EXPLOSION-GENERATED PERMEABILITY AROUND GEOTHERMAL WELLS C. R. McKee and M. E. Hanson Lawrence Livermore Laboratory University of California Livermore, California 94550 The problem of stimulating
More informationViscosity of magmas containing highly deformable bubbles
Journal of Volcanology and Geothermal Research 105 (2001) 19±24 www.elsevier.nl/locate/jvolgeores Viscosity of magmas containing highly deformable bubbles M. Manga a, *, M. Loewenberg b a Department of
More informationIAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1.
IAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1.978 PDF) http://web.mit.edu/mbuehler/www/teaching/iap2006/intro.htm
More information2.1 Strain energy functions for incompressible materials
Chapter 2 Strain energy functions The aims of constitutive theories are to develop mathematical models for representing the real behavior of matter, to determine the material response and in general, to
More informationtheory, which can be quite useful in more complex systems.
Physics 7653: Statistical Physics http://www.physics.cornell.edu/sethna/teaching/653/ In Class Exercises Last correction at August 30, 2018, 11:55 am c 2017, James Sethna, all rights reserved 9.5 Landau
More informationChapter 12. Static Equilibrium and Elasticity
Chapter 12 Static Equilibrium and Elasticity Static Equilibrium Equilibrium implies that the object moves with both constant velocity and constant angular velocity relative to an observer in an inertial
More informationA modified quarter point element for fracture analysis of cracks
ndian Journal of Engineering & Materials Sciences Vol. 14, February 007, pp. 31-38 A modified quarter point element for fracture analysis of cracks Sayantan Paul & B N Rao* Structural Engineering Division,
More informationEDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d
Institute of Industrial Science, University of Tokyo Bulletin of ERS, No. 48 (5) A TWO-PHASE SIMPLIFIED COLLAPSE ANALYSIS OF RC BUILDINGS PHASE : SPRING NETWORK PHASE Shanthanu RAJASEKHARAN, Muneyoshi
More informationRUPTURE OF FRICTIONALLY HELD INCOHERENT INTERFACES UNDER DYNAMIC SHEAR LOADING
RUPTURE OF FRICTIONALLY HELD INCOHERENT INTERFACES UNDER DYNAMIC SHEAR LOADING G. Lykotrafitis and A.J. Rosakis Graduate Aeronautical Laboratories, Mail Stop 105-50, California Institute of Technology,
More informationFRICTIONAL HEATING DURING AN EARTHQUAKE. Kyle Withers Qian Yao
FRICTIONAL HEATING DURING AN EARTHQUAKE Kyle Withers Qian Yao Temperature Change Along Fault Mode II (plain strain) crack rupturing bilaterally at a constant speed v r Idealize earthquake ruptures as shear
More informationSpectral Element simulation of rupture dynamics
Spectral Element simulation of rupture dynamics J.-P. Vilotte & G. Festa Department of Seismology, Institut de Physique du Globe de Paris, 75252 France ABSTRACT Numerical modeling is an important tool,
More informationELASTICITY AND FRACTURE MECHANICS. Vijay G. Ukadgaonker
THEORY OF ELASTICITY AND FRACTURE MECHANICS y x Vijay G. Ukadgaonker Theory of Elasticity and Fracture Mechanics VIJAY G. UKADGAONKER Former Professor Indian Institute of Technology Bombay Delhi-110092
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Overview Milestones in continuum mechanics Concepts of modulus and stiffness. Stress-strain relations Elasticity Surface and body
More informationFracture Mechanics, Damage and Fatigue Linear Elastic Fracture Mechanics - Energetic Approach
University of Liège Aerospace & Mechanical Engineering Fracture Mechanics, Damage and Fatigue Linear Elastic Fracture Mechanics - Energetic Approach Ludovic Noels Computational & Multiscale Mechanics of
More informationMoving screw dislocations in piezoelectric bimaterials
phys stat sol (b) 38 No 1 10 16 (003) / DOI 10100/pssb00301805 Moving screw dislocations in piezoelectric bimaterials Xiang-Fa Wu *1 Yuris A Dzenis 1 and Wen-Sheng Zou 1 Department of Engineering Mechanics
More informationarxiv: v1 [cond-mat.mtrl-sci] 3 Jun 2008
Laws of crack motion and phase-field models of fracture arxiv:0806.0593v1 [cond-mat.mtrl-sci] 3 Jun 2008 Vincent Hakim a a Laboratoire de Physique Statistique, CNRS-UMR8550 associé aux universités Paris
More informationENG2000 Chapter 7 Beams. ENG2000: R.I. Hornsey Beam: 1
ENG2000 Chapter 7 Beams ENG2000: R.I. Hornsey Beam: 1 Overview In this chapter, we consider the stresses and moments present in loaded beams shear stress and bending moment diagrams We will also look at
More information5 ADVANCED FRACTURE MODELS
Essentially, all models are wrong, but some are useful George E.P. Box, (Box and Draper, 1987) 5 ADVANCED FRACTURE MODELS In the previous chapter it was shown that the MOR parameter cannot be relied upon
More information202 Index. failure, 26 field equation, 122 force, 1
Index acceleration, 12, 161 admissible function, 155 admissible stress, 32 Airy's stress function, 122, 124 d'alembert's principle, 165, 167, 177 amplitude, 171 analogy, 76 anisotropic material, 20 aperiodic
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 11 Dec 2002
arxiv:cond-mat/0212257v1 [cond-mat.stat-mech] 11 Dec 2002 International Journal of Modern Physics B c World Scientific Publishing Company VISCOUS FINGERING IN MISCIBLE, IMMISCIBLE AND REACTIVE FLUIDS PATRICK
More informationTristan Cambonie (Postdoc ANR Mephystar and GeoSMEC), Harold Auradou (DR CNRS), Véronique Lazarus (MdC UPMC)
Tristan Cambonie (Postdoc ANR Mephystar and GeoSMEC), Harold Auradou (DR CNRS), Véronique Lazarus (MdC UPMC) 2 Do upper-crust fault segments result from the bottom-up propagation of mode 3 loaded cracks?
More informationThe Ultimate Load-Carrying Capacity of a Thin-Walled Shuttle Cylinder Structure with Cracks under Eccentric Compressive Force
The Ultimate Load-Carrying Capacity of a Thin-Walled Shuttle Cylinder Structure with Cracks under Eccentric Compressive Force Cai-qin Cao *, Kan Liu, Jun-zhe Dong School of Science, Xi an University of
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationMeasure what is measurable, and make measurable what is not so. - Galileo GALILEI. Dynamic fracture. Krishnaswamy Ravi-Chandar
Measure what is measurable, and make measurable what is not so. - Galileo GALILEI Dynamic fracture Krishnaswamy Ravi-Chandar Lecture 4 presented at the University of Pierre and Marie Curie May, 04 Center
More informationDevelopment and application of time-lapse ultrasonic tomography for laboratory characterisation of localised deformation in hard soils / soft rocks
Development and application of time-lapse ultrasonic tomography for laboratory characterisation of localised deformation in hard soils / soft rocks Erika Tudisco Research Group: Stephen A. Hall Philippe
More informationChapter 3 LAMINATED MODEL DERIVATION
17 Chapter 3 LAMINATED MODEL DERIVATION 3.1 Fundamental Poisson Equation The simplest version of the frictionless laminated model was originally introduced in 1961 by Salamon, and more recently explored
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationChapter Two: Mechanical Properties of materials
Chapter Two: Mechanical Properties of materials Time : 16 Hours An important consideration in the choice of a material is the way it behave when subjected to force. The mechanical properties of a material
More information3D dynamic crack propagation analysis with PDS-FEM
JAMSTEC-R IFREE Special Issue, November 2009 3D dynamic crack propagation analysis with PDS-FEM 1*, Hide Sakaguchi 1, Kenji Oguni 2, Muneo Hori 2 discretizations. Keywords +81-45-778-5972 lalith@jamstec.go.jp
More informationAdvanced Strength of Materials Prof S. K. Maiti Mechanical Engineering Indian Institute of Technology, Bombay. Lecture 27
Advanced Strength of Materials Prof S. K. Maiti Mechanical Engineering Indian Institute of Technology, Bombay Lecture 27 Last time we considered Griffith theory of brittle fracture, where in it was considered
More informationAlternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering
University of Liège Aerospace & Mechanical Engineering Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE Van Dung NGUYEN Innocent NIYONZIMA Aerospace & Mechanical engineering
More informationNumerical modeling of sliding contact
Numerical modeling of sliding contact J.F. Molinari 1) Atomistic modeling of sliding contact; P. Spijker, G. Anciaux 2) Continuum modeling; D. Kammer, V. Yastrebov, P. Spijker pj ICTP/FANAS Conference
More informationAbstract. 1 Introduction
Contact analysis for the modelling of anchors in concrete structures H. Walter*, L. Baillet** & M. Brunet* *Laboratoire de Mecanique des Solides **Laboratoire de Mecanique des Contacts-CNRS UMR 5514 Institut
More information