Wave Dispersion and Basic Concepts of Peridynamics Compared To

Total Page:16

File Type:pdf, Size:1020Kb

Wave Dispersion and Basic Concepts of Peridynamics Compared To Wave Dispersion and Basic Concepts of Peridynamics Zdenek P. Bazant1 McCormick Institute Professor and W.P. Murphy Professor of Civil and Mechanical Engineering Compared to Classical Nonlocal and Materials Science, Northwestern University, Damage Models 2145 Sheridan Road, CEE/A135, The spectral approach is used to examine the wave dispersion in linearized bond-based Evanston, IL 60208 and state-based peridynamics in one and two dimensions, and comparisons with the clas- e-mail: [email protected] sical nonlocal models for damage are made. Similar to the classical nonlocal models, the peridynamic dispersion of elastic waves occurs for high frequencies. It is shown to be Wen Luo stronger in the state-based than in the bond-based version, with multiple wavelengths Northwestern University, giving a vanishing phase velocity, one of them longer than the horizon. In the bond-based 2145 Sheridan Road, and state-based, the nonlocality of elastic and inelastic behaviors is coupled, i.e., the dis- CEE/A135, persion of elastic and inelastic waves cannot be independently controlled. In conse- Evanston, IL 60208 quence, the difference between: (1) the nonlocality due to material characteristic length for softening damage, which ensures stability of softening damage and serves as the Viet T. Chau localization limiter, and (2) the nonlocality due to material heterogeneity cannot be dis- Northwestern University, tinguished. This coupling of both kinds of dispersion is unrealistic and similar to the orig- 2145 Sheridan Road, inal 1984 nonlocal model for damage which was in 1987 abandoned and improved to be CEE/A135, nondispersive or mildly dispersive for elasticity but strongly dispersive for damage. With Evanston, IL 60208 the same regular grid of nodes, the convergence rates for both the bond-based and state- based versions are found to be slower than for the finite difference methods. It is shown Miguel A. Bessa that there exists a limit case of peridynamics, with a micromodulus in the form of a Delta Northwestern University, function spiking at the horizon. This limit case is equivalent to the unstabilized imbricate 2145 Sheridan Road, continuum and exhibits zero-energy periodic modes of instability. Finally, it is empha- CEE/A135, sized that the node-skipping force interactions, a salient feature of peridynamics, are Evanston, IL 60208 physically unjustified (except on the atomic scale) because in reality the forces get trans- mitted to the second and farther neighboring particles (or nodes) through the displace- ments and rotations of the intermediate particles, rather than by some potential permeating particles as on the atomic scale. [DOI: 10.1115/1.4034319] 1 Introduction types of micromoduli. Zimmermann [25] compared the dispersion behavior of discrete bond-based peridynamic networks based on Peridynamics is a recently promulgated radically new nonlocal intermolecular forces with classical gradient and integral type theory of solid mechanics with an intriguing new name.2 It nonlocal continua. For integral-type nonlocal elasticity and for eschews Cauchy’s concept of stress and strain and is based on various gradient formulations, the dispersion relation was studied, central force interactions between material points at various finite e.g., by Jirasek [26]. However, a comprehensive investigation of distances3 [e.g., Refs. 2–18]. The novel idea of peridynamics [13] the various versions of peridynamics is still lacking. is to perform the nonlocal spatial (weighted) averaging over dis- This work analyzes the wave propagation behavior obtained placement differences, whereas in the classical nonlocal damage from the continuous and discrete descriptions for both bond- and model [19–23], the averaging is performed over inelastic strains state-based peridynamics. The effects originating from the discre- or over a scalar softening parameter that depends on the displace- tization of continuum peridynamics are isolated from the effects ment derivatives. The introduction of nonlocality into the contin- of the governing laws of the different peridynamic formulations. uum model is essential for preventing spurious localization of A comparison with the classical nonlocal continua is also pro- softening damage [19]. On the other hand, an artificial or exces- vided. The main objectives of this study are fourfold: sive nonlocality introduced into elasticity may raise challenging questions when compared to the nonlocal damage model, (1) Analytical examination of the wave dispersion properties especially with respect to wave dispersion, material property char- of linear peridynamics; acterization, and numerical discretization. (2) Comparison of peridynamic continuum to the earlier con- This study is focused on the wave dispersion properties. Silling tinuous nonlocal models for quasibrittle materials [13], and Weckner and Abeyaratne [24], obtained a general [19–23,27–33] (although comparisons to the extensive expression for the dispersion relation of the bond-based peridy- work on lattice particle models for quasibrittle materials namic continuum and evaluated the expression by using different [e.g., Refs. 34–37] are beyond the present scope); (3) Comparison of important metrics such as critical time step, 1Corresponding author. numerical error as well as rate of convergence of the dis- 2In Greek, “peri” means “about” and “dyna” means “force”. Although crete peridynamic formulations with classical numerical “dynamics” is normally understood as the opposite of “statics”, “peridynamics” is methods, such as finite difference schemes and the finite applied to statics, too. element method. 3In this regard, it bears partial similarity to the 1977 “network model” of Burt and Dougill [1] for concrete, which was, however, used only to simulate the strain- (4) Discussion of some problematic or improvable physical softening constitutive relation aspects of the current forms of peridynamics. Contributed by the Applied Mechanics Division of ASME for publication in the JOURNAL OF APPLIED MECHANICS. Manuscript received June 26, 2016; final manuscript The wave dispersion analysis will pursue an approach similar received July 23, 2016; published online August 30, 2016. Editor: Yonggang Huang. to that used earlier for spurious wave reflection in nonuniform Journal of Applied Mechanics Copyright VC 2016 by ASME NOVEMBER 2016, Vol. 83 / 111004-1 Downloaded From: http://appliedmechanics.asmedigitalcollection.asme.org/ on 08/30/2016 Terms of Use: http://www.asme.org/about-asme/terms-of-use finite element meshes [38,39] and for the original nonlocal models where uS ¼ uXðÞS; t and u ¼ uXðÞ; t ; of distributed damage [40] (see also p. 926 in Ref. [27]). uS þ XS À u À X nuðÞS À u; XS À X; t ¼ (4) 2 Overview of Peridynamic Continuum Models jjuS þ XS À u À Xjj 2.1 Bond-Based Peridynamics. The continuum limit of a jjuS þ XS À u À Xjj discrete (bond-based) peridynamic network [13] is a special kind sðÞuS À u; XS À X ¼ À 1 (5) jjX À Xjj of a nonlocal continuum defined at every point X by the following S equation of motion [9,14]: Here CðÞXS À X is a specified material property called micromo- dulus which has the same dimension as f (although these dimen- qðÞX uX€ðÞ; t ¼ FXðÞ; t þ bXðÞ; t (1) sions have been perfectly justified mathematically, it is strange ð from the physical viewpoint that the material stiffness needs to be where FXðÞ; t ¼ fuX½ðÞS; t À uXðÞ; t ; XS À X; t dV (2) fundamentally characterized by quantities of typical dimensions 6 V such as N/m in 3D). where X; X ¼ Lagrangian coordinate vectors of paired continuum S 2.2 State-Based Peridynamics. The state-based version of points, u; u ¼ displacement vectors of these points, t ¼ time, S peridynamics, developed by Silling et al. [15], brings about a sig- q ¼ mass density, V ¼ nonlocal averaging domain bounded by a nificant improvement in the material modeling capability. In this horizon radius, d (or briefly “horizon”), b ¼ distributed body version, an effective strain tensor is constructed from dyadic prod- force, FXðÞ; t ¼ elastic resisting force, and f ¼ force density trans- ucts of all the pairwise force vectors emanating from a node. The mitted by the nodal pair to point of coordinate vector X in the use of tensors circumvents the Poisson ratio restrictions of bond- undeformed state, called the “pairwise response function,” whose based peridynamics (1/4 in 3D, 1/3 in 2D) and allows using the dimension is N/m6, N/m5, or N/m4, respectively, for three- standard tensorial constitutive relations (as well as the microplane dimensional (3D), two-dimensional (2D), and one-dimensional models). The average of the displacement differences over all the (1D) continuum. Force f is calculated as [9,14] (Fig. 1) pairs is used to approximate the displacement derivative, which was shown in Ref. [18] to coincide, for uniform grids, with the so- called Lanczos derivative (a particular case of the reproducing fuðÞÀ u; X À X; t ¼ CðÞX À X sðÞu À u; X À X; t S S S S S kernel approximation [41]). Connection to other mesh-free meth- nuðÞS À u; XS À X; t (3) ods was elucidated by Bessa et al. [18]. In state-based peridynamics, with correspondence to classical constitutive laws, the conservation of linear momentum equation is formulated as [15] ð  À1 qðÞX uX€ðÞ; t ¼ aðÞXS À X rðÞXS; t Á K ðÞXS V à À1 þrðÞX; t Á K ðÞX Á XS dV þ bXðÞ; t (6) where aðÞXS À X is the weight function, typically used in the peri- dynamics literature as the window (or separator) function, and KXðÞis the so-called shape tensor ð KXðÞ¼ aðÞXS À X ðÞXS À X ðÞXS À X dV (7) V and where the stress at each point of the peridynamic continuum is obtained from the constitutive law for that material, rðÞX; t rðÞFdefðÞX; t , applied to the nonlocal deformation gradi- ent tensor, which is ð FdefðÞX; t ¼ aðÞXS À X ðÞXS À X ðÞXS À X dV V Á KÀ1ðÞX dV (8) This tensor is substituted into the aforementioned constitutive law giving r, which is then substituted into Eq.
Recommended publications
  • Chapter 3 Constitutive Relations
    Chapter 3 Constitutive Relations Maxwell’s equations define the fields that are generated by currents and charges. However, they do not describe how these currents and charges are generated. Thus, to find a self-consistent solution for the electromagnetic field, Maxwell’s equations must be supplemented by relations that describe the behavior of matter under the influence of fields. These material equations are known as constitutive relations. The constitutive relations express the secondary sources P and M in terms of the fields E and H, that is P = f[E] and M = f[H].1 According to Eq. (1.19) this is equivalent to D = f[E] and B = f[H]. If we expand these relations into power series ∂D 1 ∂2D D = D + E + E2 + .. (3.1) 0 ∂E 2 ∂E2 E=0 E=0 we find that the lowest-order term depends linearly on E. In most practical situ- ations the interaction of radiation with matter is weak and it suffices to truncate the power series after the linear term. The nonlinear terms come into play when the fields acting on matter become comparable to the atomic Coulomb potential. This is the territory of strong field physics. Here we will entirely focus on the linear properties of matter. 1In some exotic cases we can have P = f[E, H] and M = f[H, E], which are so-called bi- isotropic or bi-anisotropic materials. These are special cases and won’t be discussed here. 33 34 CHAPTER 3. CONSTITUTIVE RELATIONS 3.1 Linear Materials The most general linear relationship between D and E can be written as2 ∞ ∞ D(r, t) = ε ε˜(r r′, t t′) E(r′, t′) d3r′ dt′ , (3.2) 0 − − −∞Z−∞Z which states that the response D at the location r and at time t not only depends on the excitation E at r and t, but also on the excitation E in all other locations r′ and all other times t′.
    [Show full text]
  • 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]
  • Modeling Optical Metamaterials with Strong Spatial Dispersion
    Fakultät für Physik Institut für theoretische Festkörperphysik Modeling Optical Metamaterials with Strong Spatial Dispersion M. Sc. Karim Mnasri von der KIT-Fakultät für Physik des Karlsruher Instituts für Technologie (KIT) genehmigte Dissertation zur Erlangung des akademischen Grades eines DOKTORS DER NATURWISSENSCHAFTEN (Dr. rer. nat.) Tag der mündlichen Prüfung: 29. November 2019 Referent: Prof. Dr. Carsten Rockstuhl (Institut für theoretische Festkörperphysik) Korreferent: Prof. Dr. Michael Plum (Institut für Analysis) KIT – Die Forschungsuniversität in der Helmholtz-Gemeinschaft Erklärung zur Selbstständigkeit Ich versichere, dass ich diese Arbeit selbstständig verfasst habe und keine anderen als die angegebenen Quellen und Hilfsmittel benutzt habe, die wörtlich oder inhaltlich über- nommenen Stellen als solche kenntlich gemacht und die Satzung des KIT zur Sicherung guter wissenschaftlicher Praxis in der gültigen Fassung vom 24. Mai 2018 beachtet habe. Karlsruhe, den 21. Oktober 2019, Karim Mnasri Als Prüfungsexemplar genehmigt von Karlsruhe, den 28. Oktober 2019, Prof. Dr. Carsten Rockstuhl iv To Ouiem and Adam Thesis abstract Optical metamaterials are artificial media made from subwavelength inclusions with un- conventional properties at optical frequencies. While a response to the magnetic field of light in natural material is absent, metamaterials prompt to lift this limitation and to exhibit a response to both electric and magnetic fields at optical frequencies. Due tothe interplay of both the actual shape of the inclusions and the material from which they are made, but also from the specific details of their arrangement, the response canbe driven to one or multiple resonances within a desired frequency band. With such a high number of degrees of freedom, tedious trial-and-error simulations and costly experimen- tal essays are inefficient when considering optical metamaterials in the design of specific applications.
    [Show full text]
  • Transparent Metamaterials with a Negative Refractive Index Determined by Spatial Dispersion
    Transparent metamaterials with a negative refractive index determined by spatial dispersion V.V. Slabko Siberian Federal University, Krasnoyarsk 660041, Russia Institute of Physics of the Russian Academy of Sciences, Krasnoyarsk 660036, Russia Abstract: The paper considers an opportunity for the creation of an artificial two-component metamaterial with a negative refractive index within the radio and optical frequency band, which possesses a spatial dispersion. It is shown that there exists a spectral region where, under certain ratio of volume fraction unit occupied with components, the metamaterial appears transparent without attenuating the electromagnetic wave passing through it and possessing a phase and group velocity of the opposite sign. 1. Introduction J.B. Pendry in co-authorship papers suggested the idea of possible creation of artificial metamaterials with a negative refractive index in the gigahertz frequency band [1,2]. The idea proceeded from the work by V.G. Veselago where the processes of electromagnetic wave propagation in isotropic materials with simultaneously negative permittivity ( ( ) 0) and magnetic permeability ( ( ) 0) was theoretically considered [3]. The latter leads to negative refraction index ( n 0 ), and phase propagation (wave vector k ) and energy-flow (Poynting vector S ) appear counter-directed, with the left-handed triplet of the vectors of magnetic field, electrical field, and the wave vector. Those were called optical negative-index materials (NIM). Interest to such metamaterials is aroused both by possible observation of extremely unusual properties, e.g. a negative refraction at the materials boundary with positive and negative refractive index, Doppler’s and Vavilov-Cherenkov’s effect, nonlinear wave interactions; and by possibilities in the solution of a number of practical tasks, which can overcome the resolution threshold diffracting limit in optical devices and creation of perfect shielding or “cloaking” [4-7].
    [Show full text]
  • Spatial Dispersion Consideration in Dipoles
    Basics of averaging of the Maxwell equations A. Chipouline, C. Simovski, S. Tretyakov Abstract Volume or statistical averaging of the microscopic Maxwell equations (MEs), i.e. transition from microscopic MEs to their macroscopic counterparts, is one of the main steps in electrodynamics of materials. In spite of the fundamental importance of the averaging procedure, it is quite rarely properly discussed in university courses and respective books; up to now there is no established consensus about how the averaging procedure has to be performed. In this paper we show that there are some basic principles for the averaging procedure (irrespective to what type of material is studied) which have to be satisfied. Any homogenization model has to be consistent with the basic principles. In case of absence of this correlation of a particular model with the basic principles the model could not be accepted as a credible one. Another goal of this paper is to establish the averaging procedure for metamaterials, which is rather close to the case of compound materials but should include magnetic response of the inclusions and their clusters. We start from the consideration of bulk materials, which means in the vast majority of cases that we consider propagation of an electromagnetic wave far from the interfaces, where the eigenwave in the medium has been already formed and stabilized. A basic structure for discussion about boundary conditions and layered metamaterials is a subject of separate publication and will be done elsewhere. 1. Material equations
    [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]
  • Operator Approach to Effective Medium Theory to Overcome a Breakdown of Maxwell Garnett Approximation
    Downloaded from orbit.dtu.dk on: Oct 05, 2021 Operator approach to effective medium theory to overcome a breakdown of Maxwell Garnett approximation Popov, Vladislav; Lavrinenko, Andrei; Novitsky, Andrey Published in: Physical Review B Link to article, DOI: 10.1103/PhysRevB.94.085428 Publication date: 2016 Document Version Publisher's PDF, also known as Version of record Link back to DTU Orbit Citation (APA): Popov, V., Lavrinenko, A., & Novitsky, A. (2016). Operator approach to effective medium theory to overcome a breakdown of Maxwell Garnett approximation. Physical Review B, 94(8), [085428]. https://doi.org/10.1103/PhysRevB.94.085428 General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. Users may download and print one copy of any publication from the public portal for the purpose of private study or research. You may not further distribute the material or use it for any profit-making activity or commercial gain You may freely distribute the URL identifying the publication in the public portal If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim. PHYSICAL REVIEW B 94, 085428 (2016) Operator approach to effective medium theory to overcome a breakdown of Maxwell Garnett approximation Vladislav Popov* Department of Theoretical Physics and Astrophysics, Belarusian State University, Nezavisimosti avenue 4, 220030, Minsk, Republic of Belarus Andrei V.
    [Show full text]
  • Vector Circuit Theory for Spatially Dispersive Uniaxial Magneto-Dielectric Slabs 3
    Machine Copy for Proofreading, JEWA (Ikonen et al.) VECTOR CIRCUIT THEORY FOR SPATIALLY DISPERSIVE UNIAXIAL MAGNETO-DIELECTRIC SLABS Pekka Ikonen, Mikhail Lapine, Igor Nefedov, and Sergei Tretyakov Radio Laboratory/SMARAD, Helsinki University of Technology P.O. Box 3000, FI-02015 TKK, Finland. Abstract—We present a general dyadic vector circuit formalism, applicable for uniaxial magneto-dielectric slabs, with strong spatial dispersion explicitly taken into account. This formalism extends the vector circuit theory, previously introduced only for isotropic and chiral slabs. Here we assume that the problem geometry imposes strong spatial dispersion only in the plane, parallel to the slab interfaces. The difference arising from taking into account spatial dispersion along the normal to the interface is briefly discussed. We derive general dyadic impedance and admittance matrices, and calculate corresponding transmission and reflection coefficients for arbitrary plane wave incidence. As a practical example, we consider a metamaterial slab built of conducting wires and split-ring resonators, and show that neglecting spatial dispersion and uniaxial nature in this structure leads to dramatic errors in calculation of transmission characteristics. 1 Introduction 2 Transmission matrix 3 Impedance and admittance matrices 4 Reflection and transmission dyadics 5 Practically realizable slabs 6 Specific examples 7 Conclusions Acknowledgment References 1. INTRODUCTION arXiv:physics/0605243v1 [physics.class-ph] 29 May 2006 It is well known that many problems dealing with reflections from multilayered media can be solved using the transmission line analogy when the eigen-polarizations are studied separately (see e.g. [1]). In this case, the amplitudes of tangential electric and magnetic fields are treated as equivalent scalar voltages and currents in the equivalent transmission line section.
    [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]