Arxiv:2108.02336V1 [Cs.CE] 5 Aug 2021 Approximation to Capture the True Material Response with High Accuracy Efficiently

Total Page:16

File Type:pdf, Size:1020Kb

Arxiv:2108.02336V1 [Cs.CE] 5 Aug 2021 Approximation to Capture the True Material Response with High Accuracy Efficiently AFRACTURE MULTISCALE MODEL FOR PERIDYNAMIC ENRICHMENT WITHIN THE PARTITION OF UNITY METHOD: PART I APREPRINT Matthias Birner Institute for Numerical Simulation, University of Bonn, Bonn, Germany Fraunhofer SCAI, Sankt Augustin, Germany Patrick Diehl Center for Computation & Technology, Louisiana State University, Baton Rouge, LA Robert Lipton Department of Mathematics, Louisiana State University, Baton Rouge, LA Center for Computation & Technology, Louisiana State University, Baton Rouge, LA Marc Alexander Schweitzer Institute for Numerical Simulation, University of Bonn, Bonn, Germany Fraunhofer SCAI, Sankt Augustin, Germany August 6, 2021 ABSTRACT Partition of unity methods (PUM) are of domain decomposition type and provide the opportunity for multiscale and multiphysics numerical modeling. Different physical models can exist within a PUM scheme for handling problems with zones of linear elasticity and zones where fractures occur. Here, the peridynamic (PD) model is used in regions of fracture and smooth PUM is used in the surrounding linear elastic media. The method is a so-called global-local enrichment strategy, see [18, 1]. The elastic fields of the undamaged media provide appropriate boundary data for the localized PD simulations. The first steps for a combined PD/PUM simulator are presented. In part I of this series, we show that the local PD approximation can be utilized to enrich the global PUM arXiv:2108.02336v1 [cs.CE] 5 Aug 2021 approximation to capture the true material response with high accuracy efficiently. Test problems are provided demonstrating the validity and potential of this numerical approach. A PREPRINT -AUGUST 6, 2021 1 Introduction fracture growth) to validate the fundamental assumptions of our approach. As this paper provides a proof of con- cept, the reference cases considered are not too complex Peridynamics (PD) is a non-local generalization of clas- in geometry and loading. sical continuum mechanics [58] and naturally allows for the direct modelling of fracture phenomena and damage The paper is structured as follows: Section 2 introduces within its original formulation [13]. It can give good the preliminaries of the PUM and PD theory. Section 3 agreement with experimental data [14] in complex defor- sketches the approach to use peridynamics to construct mation scenarios. However, the computational expense appropriate multiscale enrichments on a local sub-domain incurred by exclusive use of peridynamics for large prob- for the global PUM. In Section 4 numerical results are pre- lems where damage and fracture only appear across small sented, which validate isolated steps of the overall multi- regions and most of the body behaves linearly elastic, can scale enrichment scheme. Specifically, we validate that be needlessly high. Therefor adaptive approaches have PD and PUM simulate a similar enough physical problem been proposed for such problems, e.g. force-based blend- despite different model formulations, in order to be able ing [52, 53, 56], the Arlequin approach [22], the variable to combine them. The issues here primarily are the differ- horizon method [57, 41], or discrete coupling [37, 55, 61, ent treatment of boundary conditions and the underlying 38, 30, 20, 62]. These methods directly couple PD and fi- material models. Moreover, we show examples of infor- nite element (FE) models to overcome the computational mation transfer in both directions: using the PD computed expense. These so-called concurrent coupling approaches crack path in a PUM simulation and using a PUM dis- introduce a non-physical interface (or overlapping hand- placement as boundary conditions on a PD problem. The shake zone) where the two computational models, PD and results clearly show the validity and potential of the pre- FE, must agree on the computed deformation behavior in sented approach. Finally, Section 5 concludes the paper a delicate way, and it is not obvious how computational with some remarks on future research steps. artifacts at this interface can be avoided in general. For example in [42] an out-of-balance force in the coupling region due to a lack of balance between the local and non- 2 Preliminaries local tractions is reported. For more general details we refer to the review paper [6]. 2.1 Classical continuum mechanics An alternative to the concurrent coupling approaches are Before we introduce the two models, we will look briefly so-called hierarchical coupling techniques where scale into the basics of classical continuum mechanics (CCM). separation is assumed, and a damage model together with Figure 1 shows a continuum in the reference configura- 3 an elastic model are utilized successively, exchanging in- tion W0 ⊂ R where a material point is identified with its 3 formation in a global and accumulated fashion, see e.g. position x 2 R . The reference configuration W0 refers to [19] and references therein. For instance, in multiscale the shape of the continuum at rest with no internal forces. 2 FEM like FE methods two FE calculations are carried out The deformation f : [0;T] × R3 ! R3 of a material point in a nested manner, one at the macroscale and the other at x 2 W0 to the so-called current configuration W(t) at time the microscale (in the vicinity of every quadrature point t 2 [0;T] is given as of the macroscale model). Obviously, such approaches (t;x) := id(x) + u(t;x) = x(t;x) (1) yield a tremendous computational cost and are typically f not trivial to implement which prohibits their widespread where u : [0;T] × R3 ! R3 refers to the displacement use. A number of multiscale PD concepts have been pro- posed [36, 54, 60, 3, 45, 60, 44, 32], however, many open u(t;x) := x(t;x) − x. (2) questions remain. The stretch s : [0;T] × R3 × R3 ! R3 between two ma- 0 We will pursue yet another approach to overcome the terial points x and x in the reference configuration W0 is computational expense of PD based on the partition of defined by unity (PU) approach [47]. The overall process we envi- s(t;x;x0) := f(t;x0) − f(t;x) = x(t;x) − x(t;x0). (3) sion is to construct a combined PD/PU simulator which automatically determines the region where PD should be employed so that the resulting local PD approximation 2.2 Partition of unity method can be utilized to construct a respective multiscale en- richment function for the global partition of unity method Partition of unity methods are a class of methods to (PUM) approximation to incorporate the true material re- numerically solve partial differential equations (PDEs). sponse including fracture growth with high accuracy effi- They were introduced in [39] to overcome limitations in ciently. In this paper, we present the first steps towards the the choice of basis functions of classical finite element realization of this process. Namely, we study the compati- methods (FEM). The key concept of a PUM is the use of bility of the different material models and the methods re- a compactly supported PU, that covers the computational quired for combining them. To accomplish this we seek to domain W. To each PU function ji a local approxima- outline the transfer of information between them. To this tion space Vi is attached, which yields a global, finite- PU end, we study several reference cases (with and without dimensional space V = ∑i jiVi, that is then used in a 2 A PREPRINT -AUGUST 6, 2021 3 of dimension ni, where i are spaces of polynomials of W0 f : W0 ! R W(t) P P degree d i . The spaces Ei denote so-called enrichment spaces of dimension dEi . The latter are arbitrary functions locally incorporated into the simulation, which we can x choose with respect to the problem at hand. We can ei- x(t;x) ther obtain the Ei from a-priori analytic knowledge about the structure of the solution to a problem. This is e.g. done in 2D fracture mechanics problems, where we have Figure 1: The continuum in the reference configura- 3 an analytic expansion of the solution around a crack tip tion W0 and after the deformation f : W0 ! R with available. Or the enrichments themselves are results of det(grad f) > 0 in the current configuration W(t) at time other simulations. In the PUM the global approximation t. space then reads as PU V := ∑jiVi = ∑jiPi + jiEi: (8) Galerkin approach. An advantage over classical FEMs is i i the ability to incorporate arbitrary basis functions. Here Observe, that we do not assume the discrete functions the intent is to use only a few problem specific basis func- dPi dEi ! tions, thereby requiring less degrees of freedom (DOF) PU s s t+dPi u = ∑ji ∑ ui yi + ∑ ui h; (9) to attain the desired global accuracy. Well-known in- i s=1 t=1 stances of FE based PUMs are the GFEM by Duarte and i.e. the employed basis functions, to be interpolatory. That Babuskaˇ [17] and the XFEM by Moes,¨ Dolbow and Be- j lytschko [40]. The particular meshfree PUM employed is, the coefficients ui do not necessarily correspond to in this paper was introduced in [47] and is referred to function values at specific points. We therefore use the as the PUM in the following. All computations in this direct splitting of the local spaces Vi presented in [48] study were carried out using the PUMA software frame- to enforce Dirichlet boundary conditions. By the above- work [49] developed at Fraunhofer SCAI. In the follow- mentioned flat-top property of the PU, it is possible to ing we present only a very short review of the PUM and construct a basis transformation that guarantees numerical refer the reader to [47] for details. stability and linear independence of the basis. The latter is of special importance when using enrichment functions Given a computational domain W, we assume to have a that are computed on-the-fly.
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]
  • 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.
    [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]