Convergence, Adaptive Refinement, and Scaling in 1D Peridynamics

Total Page:16

File Type:pdf, Size:1020Kb

Convergence, Adaptive Refinement, and Scaling in 1D Peridynamics University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Faculty Publications from the Department of Mechanical & Materials Engineering, Engineering Mechanics Department of 2009 Convergence, adaptive refinement, and scaling in 1D peridynamics Florin Bobaru Ph.D. University of Nebraska at Lincoln, [email protected] Mijia Yabg Ph.D. University of Nebraska at Lincoln Leonardo F. Alves M.S. University of Nebraska at Lincoln Stewart A. Silling Ph.D. Sandia National Laboratories Ebrahim Askari Ph.D. Boeing Co., Bellevue, WA See next page for additional authors Follow this and additional works at: https://digitalcommons.unl.edu/engineeringmechanicsfacpub Part of the Applied Mechanics Commons, Computational Engineering Commons, Engineering Mechanics Commons, and the Mechanics of Materials Commons Bobaru, Florin Ph.D.; Yabg, Mijia Ph.D.; Alves, Leonardo F. M.S.; Silling, Stewart A. Ph.D.; Askari, Ebrahim Ph.D.; and Xu, Jifeng Ph.D., "Convergence, adaptive refinement, and scaling in 1D peridynamics" (2009). Faculty Publications from the Department of Engineering Mechanics. 68. https://digitalcommons.unl.edu/engineeringmechanicsfacpub/68 This Article is brought to you for free and open access by the Mechanical & Materials Engineering, Department of at DigitalCommons@University of Nebraska - Lincoln. It has been accepted for inclusion in Faculty Publications from the Department of Engineering Mechanics by an authorized administrator of DigitalCommons@University of Nebraska - Lincoln. Authors Florin Bobaru Ph.D., Mijia Yabg Ph.D., Leonardo F. Alves M.S., Stewart A. Silling Ph.D., Ebrahim Askari Ph.D., and Jifeng Xu Ph.D. This article is available at DigitalCommons@University of Nebraska - Lincoln: https://digitalcommons.unl.edu/ engineeringmechanicsfacpub/68 INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2009; 77:852–877 Published online 26 August 2008 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/nme.2439 Convergence, adaptive refinement, and scaling in 1D peridynamics Florin Bobaru1,∗,†,‡, Mijia Yang1,4, Leonardo Frota Alves1, Stewart A. Silling2, Ebrahim Askari3 and Jifeng Xu3 1Department of Engineering Mechanics, University of Nebraska-Lincoln, Lincoln, NE, U.S.A. 2Multiscale Dynamic Material Modeling Department, Sandia National Laboratories, Albuquerque, NM, U.S.A. 3Boeing Co., Bellevue, WA, U.S.A. 4Department of Civil and Environmental Engineering, University of Texas - San Antonio, San Antonio, TX, U.S.A. SUMMARY We introduce here adaptive refinement algorithms for the non-local method peridynamics, which was proposed in (J. Mech. Phys. Solids 2000; 48:175–209) as a reformulation of classical elasticity for discontinuities and long-range forces. We use scaling of the micromodulus and horizon and discuss the particular features of adaptivity in peridynamics for which multiscale modeling and grid refinement are closely connected. We discuss three types of numerical convergence for peridynamics and obtain uniform convergence to the classical solutions of static and dynamic elasticity problems in 1D in the limit of the horizon going to zero. Continuous micromoduli lead to optimal rates of convergence independent of the grid used, while discontinuous micromoduli produce optimal rates of convergence only for uniform grids. Examples for static and dynamic elasticity problems in 1D are shown. The relative error for the static and dynamic solutions obtained using adaptive refinement are significantly lower than those obtained using uniform refinement, for the same number of nodes. Copyright q 2008 John Wiley & Sons, Ltd. Received 21 December 2007; Revised 3 June 2008; Accepted 8 July 2008 KEY WORDS: peridynamics; non-local methods; adaptive refinement; convergence; multiscale modeling ∗Correspondence to: Florin Bobaru, Department of Engineering Mechanics, University of Nebraska-Lincoln, Lincoln, NE, U.S.A. †E-mail: [email protected] ‡Associate Professor. Contract/grant sponsor: Boeing Co.; contract/grant number: DE-AC04-94AL85000 Contract/grant sponsor: Computer Science Research Foundation and Computer Science Research Institute at Sandia National Laboratories Copyright q 2008 John Wiley & Sons, Ltd. CONVERGENCE, ADAPTIVE REFINEMENT, AND SCALING IN 1D PERIDYNAMICS 853 1. INTRODUCTION 1.1. Literature review Adaptive grid refinement techniques have been introduced in numerical simulations in order to gain both accuracy and efficiency. Capturing stress and strain concentrations near cracks, for example, requires dense computational grids. Adaptive methods permit reduction of grid density away from steep gradients thus rendering the computational problems tractable. Error estimators are required to determine where grid refinement is necessary. Adaptive methods for finite element methods (FEM) have been studied for the past three decades (see, e.g. [1, 2]). Automatic and goal-oriented h, p,andhp-adaptive strategies for FEM have been introduced and discussed in [3–6]. While many applications of adaptive refinement methods in FEM have been developed for static problems, adaptivity for dynamic problems has also been studied (e.g. [7–10]). In 1D time-dependent problems, several additional issues have to be dealt with, such as (see [7, 8]) spurious reflections of waves over non-uniform grids. A local error estimator and an adaptive time-stepping procedure for dynamic analysis are proposed in [9]. Dispersion effects due to FE spatial discretization are also discussed in [11]. Adaptive finite element approximations for the acoustic wave equation using a global duality argument and Galerkin orthogonality are presented in [12], while in [10] automatic energy conserving space–time refinement techniques for dynamic FEM are shown. The method presented in [10] performs automatic refining and coarsening for the wave propagation problem in 1D. Adaptive discontinuous Galerkin FEM for dynamic problems are developed in [13]. In multiple dimensions, the grid generation for adaptive refinement becomes more complex and special algorithms have been proposed. One of the simplest is using quad tree structures (in 2D) and has been used in the context of finite volume method in [14–16], for FEM in [17], and with the meshfree element-free Galerkin method [18], where level-1 quad-tree structures are employed even though this is not explicitly mentioned in the paper. The level-1 (or 1-irregular) grid refinement has been originally introduced in [1]. In the present study we introduce algorithms and show convergence results regarding adap- tive grid refinement algorithms for static and dynamic elasticity problems in 1D for a non-local method, the peridynamic formulation. Non-local methods are necessary in modeling of localization phenomena, such as fracture, damage, shear bands, for reasons described in detail in the review paper [19]. Adaptive refinement for other non-local methods applied to damage problems and discretized using the finite element method has appeared in, for example, [20–23]. These methods cannot be directly applied to peridynamics because the peridynamic formulation, first introduced in [24], is different from other non-local methods. We discuss this topic briefly below. 1.2. Overview of peridynamics The peridynamic method, introduced in [24], is a non-local formulation of continuum mechanics in which each material point is connected (via peridynamic bonds) with nodes in a certain region around it and not only with its nearest neighbors. The method is particularly well suited for dealing with cracks and damage in solid mechanics (see [25, 26]) especially in situations where the crack path, for example, is not known in advance. The method also allows for easy introduction of long- range forces, such as van der Waals interactions and examples of deformation and damage at the nanoscale are treated in [26, 27]. In fact, all forces in the peridynamic formulation are long-range. Other non-local methods have been proposed over the years: see, e.g. [28–33]. For a recent review Copyright q 2008 John Wiley & Sons, Ltd. Int. J. Numer. Meth. Engng 2009; 77:852–877 DOI: 10.1002/nme 854 F. BOBARU ET AL. of other non-local methods, we refer the reader to [19]. Peridynamics differs from other integral non-local methods in at least two crucial aspects: (1) spatial derivatives are eliminated while other methods average infinitesimal stresses over a certain area and (2) the way damage and fracture are introduced in the method (see [26]). A description of a meshfree numerical implementation for the peridynamic formulation is given in [25], where bond failure is related to the classical energy release rate in brittle fracture. The method is used in [26] to simulate the tearing of non-linear membranes and failure of nanofiber networks (see also [27]). Spontaneous fracture initiation and dynamic propagation are captured and membrane wrinkling and the influence of Van der Waals forces in nanofiber networks are discussed. Classical continuum mechanics has been applied to explain the behavior of structural bodies which remain continuous as they deform under external loads. However, the effectiveness of compu- tational methods based on classical continuum mechanics, such as finite elements, in modeling material failure has lagged far behind their capabilities in traditional stress analysis. This difficulty arises because these methods solve the classical equations of continuum mechanics, which are a set of partial differential equations. These equations cannot be applied directly across disconti- nuities resulted from material
Recommended publications
  • 14Th U.S. National Congress on Computational Mechanics
    14th U.S. National Congress on Computational Mechanics Montréal • July 17-20, 2017 Congress Program at a Glance Sunday, July 16 Monday, July 17 Tuesday, July 18 Wednesday, July 19 Thursday, July 20 Registration Registration Registration Registration Short Course 7:30 am - 5:30 pm 7:30 am - 5:30 pm 7:30 am - 5:30 pm 7:30 am - 11:30 am Registration 8:00 am - 9:30 am 8:30 am - 9:00 am OPENING PL: Tarek Zohdi PL: Andrew Stuart PL: Mark Ainsworth PL: Anthony Patera 9:00 am - 9.45 am Chair: J.T. Oden Chair: T. Hughes Chair: L. Demkowicz Chair: M. Paraschivoiu Short Courses 9:45 am - 10:15 am Coffee Break Coffee Break Coffee Break Coffee Break 9:00 am - 12:00 pm 10:15 am - 11:55 am Technical Session TS1 Technical Session TS4 Technical Session TS7 Technical Session TS10 Lunch Break 11:55 am - 1:30 pm Lunch Break Lunch Break Lunch Break CLOSING aSPL: Raúl Tempone aSPL: Ron Miller aSPL: Eldad Haber 1:30 pm - 2:15 pm bSPL: Marino Arroyo bSPL: Beth Wingate bSPL: Margot Gerritsen Short Courses 2:15 pm - 2:30 pm Break-out Break-out Break-out 1:00 pm - 4:00 pm 2:30 pm - 4:10 pm Technical Session TS2 Technical Session TS5 Technical Session TS8 4:10 pm - 4:40 pm Coffee Break Coffee Break Coffee Break Congress Registration 2:00 pm - 8:00 pm 4:40 pm - 6:20 pm Technical Session TS3 Poster Session TS6 Technical Session TS9 Reception Opening in 517BC Cocktail Coffee Breaks in 517A 7th floor Terrace Plenary Lectures (PL) in 517BC 7:00 pm - 7:30 pm Cocktail and Banquet 6:00 pm - 8:00 pm Semi-Plenary Lectures (SPL): Banquet in 517BC aSPL in 517D 7:30 pm - 9:30 pm Viewing of Fireworks bSPL in 516BC Fireworks and Closing Reception Poster Session in 517A 10:00 pm - 10:30 pm on 7th floor Terrace On behalf of Polytechnique Montréal, it is my pleasure to welcome, to Montreal, the 14th U.S.
    [Show full text]
  • Peridynamics Review
    Peridynamics Review Ali Javilia,∗, Rico Morasataa, Erkan Oterkusb, Selda Oterkusb aDepartment of Mechanical Engineering, Bilkent University, 06800 Ankara, Turkey bDepartment of Naval Architecture, Ocean and Marine Engineering, University of Strathclyde, Glasgow, United Kingdom Abstract Peridynamics (PD) is a novel continuum mechanics theory established by Stewart Silling in 2000 [1]. The roots of PD can be traced back to the early works of Gabrio Piola according to dell'Isola et al. [2]. PD has been attractive to researchers as it is a nonlocal formulation in an integral form; unlike the local differential form of classical continuum mechanics. Although the method is still in its infancy, the literature on PD is fairly rich and extensive. The prolific growth in PD applications has led to a tremendous number of contributions in various disciplines. This manuscript aims to provide a concise description of the peridynamic theory together with a review of its major applications and related studies in different fields to date. Moreover, we succinctly highlight some lines of research that are yet to be investigated. Keywords: Peridynamics, Continuum mechanics, Fracture, Damage, Nonlocal elasticity 1. Introduction The analysis of material behavior due to various types of loading and boundary conditions may be carried out at different length scales. At the nanoscale, often Molecular Dynamics (MD) is used to analyze the behavior of atoms and molecules. However, simulation time, integration algorithms and the required computational power are limiting factors to the method’s capabilities. On the other hand, at the macroscale, the behavior of a body is usually analyzed using the Classical Continuum Mechanics (CCM) formalism.
    [Show full text]
  • Discretized Bond-Based Peridynamics for Solid Mechanics
    DISCRETIZED BOND-BASED PERIDYNAMICS FOR SOLID MECHANICS By Wenyang Liu A DISSERTATION Submitted to Michigan State University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Civil Engineering 2012 ABSTRACT DISCRETIZED BOND-BASED PERIDYNAMICS FOR SOLID MECHANICS By Wenyang Liu The numerical analysis of spontaneously formed discontinuities such as cracks is a long- standing challenge in the engineering field. Approaches based on the mathematical frame- work of classical continuum mechanics fail to be directly applicable to describe discontinuities since the theory is formulated in partial differential equations, and a unique spatial deriva- tive, however, does not exist on the singularities. Peridynamics is a reformulated theory of continuum mechanics. The partial differential equations that appear in the classical contin- uum mechanics are replaced with integral equations. A spatial range, which is called the horizon δ, is associated with material points, and the interaction between two material points within a horizon is formed in terms of the bond force. Since material points separated by a finite distance in the reference configuration can interact with each other, the peridynamic theory is categorized as a nonlocal method. The primary focus in this research is the development of the discretized bond-based peridynamics for solid mechanics. A connection between the classical elasticity and the dis- cretized peridynamics is established in terms of peridynamic stress. Numerical micromoduli for one- and three-dimensional models are derived. The elastic responses of one- and three- dimensional peridynamic models are examined, and the boundary effect associated with the size of the horizon is discussed. A pairwise compensation scheme is introduced in this re- search for simulations of an elastic body of Poisson ratio not equal to 1/4.
    [Show full text]
  • Mathematical and Numerical Analysis for Linear Peridynamic Boundary Value Problems
    Mathematical and Numerical Analysis for Linear Peridynamic Boundary Value Problems by Hao Wu A dissertation submitted to the Graduate Faculty of Auburn University in partial fulfillment of the requirements for the Degree of Doctor of Philosophy Auburn, Alabama May 7, 2017 Keywords: peridynamics, nonlocal boundary value problems, finite-dimensional approximations, error estimates, exponential convergence. Copyright 2017 by Hao Wu Approved by Yanzhao Cao, Chair, Professor of Department of Mathematics and Statistics Junshan Lin, Asistant Professor of Department of Mathematics and Statistics Ulrich Albrecht, Professor of Department of Mathematics and Statistics Wenxian Shen, Professor of Department of Mathematics and Statistics Abstract Peridynamics is motivated in aid of modeling the problems from continuum mechanics which involve the spontaneous discontinuity forms in the motion of a material system. By replacing differentiation with integration, peridynamic equations remain equally valid both on and off the points where a discontinuity in either displacement or its spatial derivatives is located. A functional analytical framework was established in literature to study the linear bond- based peridynamic equations associated with a particular kind of nonlocal boundary condi- tion. Investigated were the finite-dimensional approximations to the solutions of the equa- tions obtained by spectral method and finite element method; as a result, two corresponding general formulas of error estimates were derived. However, according to these formulas, one can only conclude that the optimal convergence is algebraic. Based on this theoretical framework, first we show that analytic data functions produce analytic solutions. Afterwards, we prove these finite-dimensional approximations will achieve exponential convergence under the analyticity assumption of data. At the end, we validate our results by conducting a few numerical experiments.
    [Show full text]
  • Peridynamic Theory of Solid Mechanics
    Peridynamic Theory of Solid Mechanics S. A. Silling∗ R. B. Lehoucqy Sandia National Laboratories Albuquerque, New Mexico 87185-1322 USA April 28, 2010 Dedicated to the memory of James K. Knowles Contents 1 Introduction 4 1.1 Purpose of the peridynamic theory . 4 1.2 Summary of the literature . 5 1.3 Organization of this article . 11 2 Balance laws 13 2.1 Balance of linear momentum . 13 2.2 Principle of virtual work . 18 2.3 Balance of angular momentum . 19 2.4 Balance of energy . 21 2.5 Master balance law . 24 3 Peridynamic states: notation and properties 28 4 Constitutive modeling 32 4.1 Simple materials . 32 4.2 Kinematics of deformation states . 34 4.3 Directional decomposition of a force state . 34 4.4 Examples . 35 ∗Multiscale Dynamic Material Modeling Department, [email protected] yApplied Mathematics and Applications Department, [email protected] 1 4.5 Thermodynamic restrictions on constitutive models . 36 4.6 Elastic materials . 38 4.7 Bond-based materials . 39 4.8 Objectivity . 41 4.9 Isotropy . 42 4.10 Isotropic elastic solid . 43 4.11 Peridynamic material derived from a classical material . 44 4.12 Bond-pair materials . 44 4.13 Example: a bond-pair material in bending . 48 5 Linear theory 50 5.1 Small displacements . 50 5.2 Double states . 50 5.3 Linearization of an elastic constitutive model . 52 5.4 Equations of motion and equilibrium . 53 5.5 Linear bond-based materials . 54 5.6 Equilibrium in a one dimensional model . 56 5.7 Plane waves and dispersion in one dimension .
    [Show full text]
  • AN16-LS16 Abstracts
    AN16-LS16 Abstracts Abstracts are printed as submitted by the authors. Society for Industrial and Applied Mathematics 3600 Market Street, 6th Floor Philadelphia, PA 19104-2688 USA Telephone: +1-215-382-9800 Fax: +1-215-386-7999 Conference E-mail: [email protected] Conference web: www.siam.org/meetings/ Membership and Customer Service: 800-447-7426 (US & Canada) or +1-215-382-9800 (worldwide) 2 2016 SIAM Annual Meeting Table of Contents AN16 Annual Meeting Abstracts ................................................4 LS16 Life Sciences Abstracts ...................................................157 SIAM Presents Since 2008, SIAM has recorded many Invited Lectures, Prize Lectures, and selected Minisymposia from various conferences. These are available by visiting SIAM Presents (http://www.siam.org/meetings/presents.php). 2016 SIAM Annual Meeting 3 AN16 Abstracts 4 AN16 Abstracts SP1 University of Cambridge AWM-SIAM Sonia Kovalevsky Lecture - Bioflu- [email protected] ids of Reproduction: Oscillators, Viscoelastic Net- works and Sticky Situations JP1 From fertilization to birth, successful mammalian repro- Spatio-temporal Dynamics of Childhood Infectious duction relies on interactions of elastic structures with a Disease: Predictability and the Impact of Vaccina- fluid environment. Sperm flagella must move through cer- tion vical mucus to the uterus and into the oviduct, where fertil- ization occurs. In fact, some sperm may adhere to oviduc- Violent epidemics of childhood infections such as measles tal epithelia, and must change their pattern
    [Show full text]
  • Peridynamics Review
    Article Mathematics and Mechanics of Solids 1–26 Ó The Author(s) 2018 Peridynamics review Article reuse guidelines: sagepub.com/journals-permissions DOI: 10.1177/1081286518803411 journals.sagepub.com/home/mms Ali Javili Department of Mechanical Engineering, Bilkent University, Ankara, Turkey Rico Morasata Department of Mechanical Engineering, Bilkent University, Ankara, Turkey Erkan Oterkus Department of Naval Architecture, Ocean and Marine Engineering, University of Strathclyde, Glasgow, UK Selda Oterkus Department of Naval Architecture, Ocean and Marine Engineering, University of Strathclyde, Glasgow, UK Received 29 November 2017; accepted 5 September 2018 Abstract Peridynamics (PD) is a novel continuum mechanics theory established by Stewart Silling in 2000. The roots of PD can be traced back to the early works of Gabrio Piola according to dell’Isola et al. PD has been attractive to researchers as it is a non-local formulation in an integral form, unlike the local differential form of classical continuum mechanics. Although the method is still in its infancy, the literature on PD is fairly rich and extensive. The prolific growth in PD applications has led to a tremendous number of contributions in various disciplines. This manuscript aims to provide a concise description of the PD theory together with a review of its major applications and related studies in different fields to date. Moreover, we succinctly highlight some lines of research that are yet to be investigated. Keywords Peridynamics, continuum mechanics, fracture, damage, non-local elasticity 1. Introduction The analysis of material behavior owing to various types of loading and boundary conditions may be carried out at different length scales.
    [Show full text]
  • Peridynamic Models for Fatigue and Fracture in Isotropic and in Polycrystalline Materials Guanfeng Zhang [email protected]
    University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Mechanical (and Materials) Engineering -- Mechanical & Materials Engineering, Department Dissertations, Theses, and Student Research of 5-2017 Peridynamic models for fatigue and fracture in isotropic and in polycrystalline materials Guanfeng Zhang [email protected] Follow this and additional works at: http://digitalcommons.unl.edu/mechengdiss Part of the Applied Mechanics Commons, Engineering Mechanics Commons, and the Other Materials Science and Engineering Commons Zhang, Guanfeng, "Peridynamic models for fatigue and fracture in isotropic and in polycrystalline materials" (2017). Mechanical (and Materials) Engineering -- Dissertations, Theses, and Student Research. 114. http://digitalcommons.unl.edu/mechengdiss/114 This Article is brought to you for free and open access by the Mechanical & Materials Engineering, Department of at DigitalCommons@University of Nebraska - Lincoln. It has been accepted for inclusion in Mechanical (and Materials) Engineering -- Dissertations, Theses, and Student Research by an authorized administrator of DigitalCommons@University of Nebraska - Lincoln. PERIDYNAMIC MODELS FOR FATIGUE AND FRACTURE IN ISOTROPIC AND IN POLYCRYSTALLINE MATERIALS by Guanfeng Zhang A DISSERTATION Presented to the Faculty of The Graduate College at the University of Nebraska In Partial Fulfilment of Requirements For the Degree of Doctor of Philosophy Major: Engineering (Materials Engineering) Under the Supervision of Professor Florin Bobaru Lincoln, Nebraska May, 2017 PERIDYNAMIC MODELS FOR FATIGUE AND FRACTURE IN ISOTROPIC AND IN POLYCRYSTALLINE MATERIALS Guanfeng Zhang, Ph.D. University of Nebraska, 2017 Adviser: Florin Bobaru To improve design and reliability, extensive efforts has been devoted to understanding damage and failure of materials and structures using numerical simulation, as a complement of theory and experiment.
    [Show full text]
  • Modeling Tumor Growth with Peridynamics
    Biomech Model Mechanobiol (2017) 16:1141–1157 DOI 10.1007/s10237-017-0876-8 ORIGINAL PAPER Modeling tumor growth with peridynamics Emma Lejeune1 · Christian Linder1 Received: 23 November 2016 / Accepted: 11 January 2017 / Published online: 26 January 2017 © Springer-Verlag Berlin Heidelberg 2017 Abstract Computational models of tumors have the poten- 1 Introduction tial to connect observations made on the cellular and the tissue scales. With cellular scale models, each cell can be Neoplasms, commonly referred to as tumors, are defined as treated as a discrete entity, while tissue scale models typi- abnormal growth of tissue (Deisboeck et al. 2011). Tumor cally represent tumors as a continuum. Though the discrete growth is problematic because tumors can disrupt organ approach often enables a more mechanistic and biologically function and infiltrate and invade healthy tissue (Gatenby driven description of cellular behavior, it is often compu- and Gawlinski 1996). Tumor growth is also the precursor tationally intractable on the tissue scale. Here, we adapt for metastasis, the spread of cancer from the primary tumor peridynamics, a theoretical and computational approach location to secondary locations, where patient outcomes sig- designed to unify the mechanics of discrete and continuous nificantly worsen (Poste and Fidler 1980). On the cellular media, for the growth of biological materials. The result is scale, tumor growth is driven by mitosis, the process by which a computational model for tumor growth that can represent a single-parent cell divides into two daughter cells (Bel- either individual cells or the tissue as a whole. We take advan- lomo et al. 2008). This differs from other forms of biological tage of the flexibility provided by the peridynamic framework growth, such as the progression of cardiovascular disease, to implement a cell division mechanism, motivated by the fact which are driven by changes in cell size and shape (Kuhl that cell division is the mechanism driving tumor growth.
    [Show full text]
  • Peridynamics: Convergence & Influence of Probabilistic Material
    6th ECCOMAS Thematic Conference on the Mechanical Response of Composites COMPOSITES 2017 J.J.C. Remmers and A. Turon (Editors) PERIDYNAMICS: CONVERGENCE & INFLUENCE OF PROBABILISTIC MATERIAL DISTRIBUTION ON CRACK INITIATION Martin Rädel?, Anna-Janina Bednarek?, Jochen Schmidty, Christian Willberg? ?Structural Mechanics Department yComposite Technology Department Institute for Composite Structures and Adaptive Systems, German Aerospace Center (DLR) Lilienthalplatz 7, D-38108 Braunschweig, Germany [email protected] Keywords: Peridynamics, Fracture, Mesh dependence. Summary: Simulations based on the peridynamic theory are a promising approach to un- derstand the processes involved in matrix failure inside fibre reinforced plastics. Before such complex simulations are carried out, the material behavior of bulk resin material as well as the influence of numerical parameters have to be investigated. In the present text, the linear elas- tic part of the material response is used to examine the convergence behavior of peridynamic simulations. Possibilities to minimize the effect of different discretization schemes are explored by means of a stochastic material distribution in correlation with scatter found in the material tests regarding the elastic material response and failure patterns. This procedure may also be used to investigate the nature of failure initiation and the robustness of the solution. δ horizon ρ mass density θ dilatation strain tensor σ Cauchy stress tensor G shear modulus neighborhood K bulk modulus V volume H Vw weighted volume ! influence function t pairwise force density x initial bond length y deformed bond length e bond extension b external forces u displacements ¨u acceleration x initial position of refer- x0 initial position of family y deformed position of ref- ence point point erence point y0 deformed position of T force vector state Y peridynamic state family point 1.
    [Show full text]
  • A Peridynamics-Based Micromechanical Modeling Approach for Random Heterogeneous Structural Materials
    materials Article A Peridynamics-Based Micromechanical Modeling Approach for Random Heterogeneous Structural Materials Sumeru Nayak 1 , R Ravinder 2, N M Anoop Krishnan 2,* and Sumanta Das 1,* 1 Civil and Environmental Engineering, University of Rhode Island, Kingston, RI 02881, USA; [email protected] 2 Department of Civil Engineering, Indian Institute of Technology Delhi, New Delhi, India; [email protected] * Correspondence: [email protected] (N.M.A.K.); [email protected] (S.D.) Received: 30 January 2020; Accepted: 4 March 2020; Published: 13 March 2020 Abstract: This paper presents a peridynamics-based micromechanical analysis framework that can efficiently handle material failure for random heterogeneous structural materials. In contrast to conventional continuum-based approaches, this method can handle discontinuities such as fracture without requiring supplemental mathematical relations. The framework presented here generates representative unit cells based on microstructural information on the material and assigns distinct material behavior to the constituent phases in the random heterogenous microstructures. The framework incorporates spontaneous failure initiation/propagation based on the critical stretch criterion in peridynamics and predicts effective constitutive response of the material. The current framework is applied to a metallic particulate-reinforced cementitious composite. The simulated mechanical responses show excellent match with experimental observations signifying efficacy of the peridynamics-based micromechanical
    [Show full text]
  • A Comparative Review of Peridynamics and Phase-Field Models for Engineering Fracture Mechanics
    A COMPARATIVE REVIEW OF PERIDYNAMICS AND PHASE-FIELD MODELS FOR ENGINEERING FRACTURE MECHANICS APREPRINT Patrick Diehl Center of Computation & Technology Louisiana State University Baton Rouge, Louisiana 70803 Email: [email protected] Robert Lipton Department of Mathematics Center of Computation & Technology Louisiana State University Baton Rouge, Louisiana 70808 Email: [email protected] Thomas Wick Leibniz Universitat¨ Hannover Institut fur¨ Angewandte Mathematik Welfengarten 1, 30167 Hannover, Germany Email: [email protected] Mayank Tyagi Craft & Hawkins Department of Petroleum Engineering Center for Computation & Technology Louisiana State University Baton Rouge, Louisiana, 70803 Email: [email protected] July 13, 2021 ABSTRACT Computational modeling of the initiation and propagation of complex fracture is central to the disci- pline of engineering fracture mechanics. This review focuses on two promising approaches: phase- field (PF) and peridynamic (PD) models applied to this class of problems. The basic concepts consisting of constitutive models, failure criteria, discretization schemes, and numerical analysis are briefly summarized for both models. Validation against experimental data is essential for all com- putational methods to demonstrate predictive accuracy. To that end, The Sandia Fracture Challenge and similar experimental data sets where both models could be benchmarked against are showcased. Emphasis is made to converge on common metrics for the evaluation of these two fracture mod- eling approaches. Both PD and PF models are assessed in terms of their computational effort and predictive capabilities with their relative advantages and challenges are summarized. A PREPRINT -JULY 13, 2021 1 Introduction and Wang [7] published a review on the usage of PD in geomaterials in 2021.
    [Show full text]