Combined Anelasticity and Nonlinearity : the Longitudinal Resonances of a Bar H

Total Page:16

File Type:pdf, Size:1020Kb

Combined Anelasticity and Nonlinearity : the Longitudinal Resonances of a Bar H COMBINED ANELASTICITY AND NONLINEARITY : THE LONGITUDINAL RESONANCES OF A BAR H. Jeong, D. Beshers To cite this version: H. Jeong, D. Beshers. COMBINED ANELASTICITY AND NONLINEARITY : THE LONGITUDI- NAL RESONANCES OF A BAR. Journal de Physique Colloques, 1987, 48 (C8), pp.C8-317-C8-322. 10.1051/jphyscol:1987846. jpa-00227150 HAL Id: jpa-00227150 https://hal.archives-ouvertes.fr/jpa-00227150 Submitted on 1 Jan 1987 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. JOURNAL DE PHYSIQUE Colloque C8, supplément au n°12, Tome 48, décembre 1987 C8-317 COMBINED ANELASTICITY AND NONLINEARITY : THE LONGITUDINAL RESONANCE OF A BAR H.S. JEONG and D.N. BESHERS Henry Krumb School of Mines, Columbia University, New York, NY 10027, U.S.A. Résumé - Nous avons ajoute un terme quadratique peu important I l'équation du modèle du solide anélastique standard et nous avons étudie son effet dans l'approximation du premier ordre. La résonance fondamen­ tale ne change pas, en première approximation, mais uneharmonique de deuxième ordre apparaît. En fait, la. deuxième harmonique résulte de deux ondes spatiales à même fréquence. Comme l'amortissement linéaire tend vers zéro, la forme mathématique de la solution change, les deux ondes dégénèrent en une seule qui n'est pas exactement sinusoïdale. Abstract - The addition of a small quadratic term to the standard anelas- tic solid model has been investigated to first order approximation. The fundamental resonance does not change, to first order, but a second har­ monic wave appears. The second harmonic actually consists of two spatial waves with one frequency. As the linear damping goes to zero, the mathe­ matical form of the solution changes, the two waves degenerating into one that is not entirely sinusoidal. I - INTRODUCTION The longitudinal resonance of a bar has often been exploited for the measurement of internal friction. The analysis of the experiments has usually been con­ ducted in an elementary fashion, and often that suffices. However, when the internal friction is nonlinear, that is to say dependent on the amplitude of oscillation, then there is a need for more careful consideration of the problem, a need universally neglected for lack of a suitable theory. A first step towards such a theory, an approximate treatment for the case of a small nonli­ near effect co-existing with a linear damping, is presented in this paper. II - THE MODEL Because nonlinear internal friction usually occurs in materials that also exhi­ bit linear damping, it is important to treat the linear and nonlinear aspects together. The basic model for linear damping is the standard anelastic solid (SAS) /1/. The longitudinal resonance of a bar made from a SAS has been treated already /2/. We consider here a bar of homogeneous material obeying a consti^ tutive equation based on the standard anelastic solid, but augmented by one nonlinear term: * ~ 2 a + ba - ae + ce + de , (1 ) where a is the longitudinal stress, e the corresponding strain, and a, b, c, d are material constants; in the usual notation, a = 1/(JU+6J), c = t a> and b - J c. The coefficient d is assumed subsequently to be small. For a long, thin bar the equations of stress equilibrium reduce to Article published online by EDP Sciences and available at http://dx.doi.org/10.1051/jphyscol:1987846 JOURNAL DE PHYSIQUE where p denotes the density of the material and u = u(x,t) the displacement of a cross-section of the bar. Strain and displacement are related by We eliminate stress in favor of u by differentiating (1) with respect to x, then using (2), and replacing E with u by (3), obtaining the equation of motion to be solved: putt + bputtt = au + cu + 2duxuxx. (4) XX xxt The bar is considered to be driven at one end, x = 0, by a periodic stress of amplitude 20 and circular frequency w and to be free at the other end, x = L. The boundaryOconditions are then, choosing the phase appropriately, o(~,t)= ooeiWt + 0 e-iwt atx=O (5) 0 We seek a steady-state solution of this problem, expecting to find resonant modes, and to examine the behavior in the vicinity of resonance. By steady- state we mean that every term in the solution must repeat in time with the fun- damental frequency; that is we assume that harmonics of the fundamental will occur, but not subharmonics. This assumption is consistent with experimental results in this laboratory. The nonlinear term necessitates care in the use of complex exponential func- tions, which are otherwise the simplest representation with which to work. The real part of any complex number Z is given by but R, as an operator, does not commute with multiplication. We shall therefore in effect require that R operate first by writing the trial solution in the form m u(x,t) = z [un(x)einwt + u *(x) e-inwt]. (8) n=O In this way, all our expressions are real, although sometimes we will not bother to write out the conjugate part. I11 - THE SOLUTION Substitution of (8) in (4) gives a very lengthy equation. Because exp(inwt) and exp(-inwt) are linearly independent functions, we require that the coefficient of each be zero, which has the effect of converting our partial differential equation into an infinite set of coupled ordinary differential equations. Each of these ordlnary equations may be viewed as consisting of two parts. One part contains four terms which correspond to the four terms of the standard linear solid at the frequency nw and the other part is an infinite series which has the coefficient d of the nonlinear term in (1). These ordinary equations are there- fore inhomogeneous. Now the solution to an inhomogeneous equation is the gem- era1 solution of the corresponding homogeneous equation plus any particular solution of the inhomogeneous equation. We can obtain the solutions to the homogeneous equations immediately, and then approximate the infinite series in terms of these known solutions, after which particular solutions can be found. This strategy is facilitated by the fact that the terms of the infinite series, while nonlinear, become products of trigonometric functions of x and therefore expressible as linear combinations of trigonometric functions of other periodic- i ties: the equations become linearized. The homogeneous equations The reduced equation for any n has the form where U' indicates dU/dx, and another equation which is the complex conjugate of (9). Equation (9) has the well-known form O" + K 2~n0= 0 "n (10) where the complex propagation constant K = kn + ian is determined by The general solution of (10) is unO = C sinK x + DncosK x (12) where the constants C and D are to be determined eventually by the boundary conditions. For the Boundary conditions given by (5) and (6). the solution is and u(x,t) = ulo(x)eiwt + C.C. (14) where C.C. stands for complex conjugate. This is the solution given previously /2/ for the resonance of a bar made of a standard linear solid. The Inhomogeneous equations Returning to the inhomogeneous equation, we start the process of approximating to the nonlinear terms by assuming that and that all other Un are small of order d, at least, and proceed to sort the terms by the powers of d, in the end neglecting all terms containing d to a power higher than 1. In the process of approximation we are also helped by assuming that 6J/J is small so that the strain amplitude at resonance is orders of magnitude less than unity. The only equations that survive this process are those for n = 0,1,2. For n = 1, the nonlinear term vanishes in this approxima.- tion. Therefore the present treatment will not, in fact, extend to the case of nonlinear damping. Rather, the approximation is one in which the fundamental is unaffected to first order, but small harmonic generation occurs. Our approxirnw tions have then reduced the infinite set of equations envisioned in (8) to three only. We focus our attention on the equations for n = 2, putting aside the equation for n = 0 which will not affect the results. Thereby, we discuss here a system with only two states, a closed system not an open one. We have then to solve the equations for n = 2 0 and the complex conjugate of (16). With U, given by (131, and denoting the JOURNAL DE PHYSIQUE coefficient of cosK1!L-x) in (13) by A, the particular solution of (16) is by the method of variation of parameters where 2 2 D = [(2K1) - K2 1 ( a + ic2w). The boundary conditions, (5) and (6), imply that the stress corresponding to U2 shall be zero at both ends of the bar, forcing the other terms in U to be proportional to the coefficient of the particular solution. Let P senote the coefficient of sin2K1 (L-x) in (17). Then the general solution of (16) is U2 = C sinK (L-x) + C2cosK2(L-x) + Psin2K1(L-x) 12 (19) where Cl = (,-K2/2K1)P, and C2 = P(K~/~K,)[CO~~K~L- cosK2Ll/sinK2L. Equation (19) shows that our solution is more complicated than we might have expected.
Recommended publications
  • 1 Revision 1 Single-Crystal Elastic Properties of Minerals and Related
    Revision 1 Single-Crystal Elastic Properties of Minerals and Related Materials with Cubic Symmetry Thomas S. Duffy Department of Geosciences Princeton University Abstract The single-crystal elastic moduli of minerals and related materials with cubic symmetry have been collected and evaluated. The compiled dataset covers measurements made over an approximately seventy year period and consists of 206 compositions. More than 80% of the database is comprised of silicates, oxides, and halides, and approximately 90% of the entries correspond to one of six crystal structures (garnet, rocksalt, spinel, perovskite, sphalerite, and fluorite). Primary data recorded are the composition of each material, its crystal structure, density, and the three independent nonzero adiabatic elastic moduli (C11, C12, and C44). From these, a variety of additional elastic and acoustic properties are calculated and compiled, including polycrystalline aggregate elastic properties, sound velocities, and anisotropy factors. The database is used to evaluate trends in cubic mineral elasticity through consideration of normalized elastic moduli (Blackman diagrams) and the Cauchy pressure. The elastic anisotropy and auxetic behavior of these materials are also examined. Compilations of single-crystal elastic moduli provide a useful tool for investigation structure-property relationships of minerals. 1 Introduction The elastic moduli are among the most fundamental and important properties of minerals (Anderson et al. 1968). They are central to understanding mechanical behavior and have applications across many disciplines of the geosciences. They control the stress-strain relationship under elastic loading and are relevant to understanding strength, hardness, brittle/ductile behavior, damage tolerance, and mechanical stability. Elastic moduli govern the propagation of elastic waves and hence are essential to the interpretation of seismic data, including seismic anisotropy in the crust and mantle (Bass et al.
    [Show full text]
  • Entropy Principle and Recent Results in Non-Equilibrium Theories
    Entropy 2014, 16, 1756-1807; doi:10.3390/e16031756 OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Article Entropy Principle and Recent Results in Non-Equilibrium Theories Vito Antonio Cimmelli 1;*, David Jou 2, Tommaso Ruggeri 3 and Péter Ván 4;5;6 1 Department of Mathematics, Computer Science and Economics, University of Basilicata, Potenza 85100, Italy 2 Departament de Física, Universitat Autònoma de Barcelona, Bellaterra, and Institut d’Estudis Catalans, Barcelona, Catalonia 08193, Spain; E-Mail: [email protected] 3 Department of Mathematics and Research Center of Applied Mathematics, University of Bologna, Bologna 40123, Italy; E-Mail: [email protected] 4 Department of Theoretical Physics, Wigner Research Centre for Physics, Institute for Particle and Nuclear Physics, Budapest 1121, Hungary; E-Mail: [email protected] 5 Department of Energy Engineering, BME, Budapest 1111, Hungary 6 Montavid Thermodynamic Research Group, Budapest 1112 , Hungary * Author to whom correspondence should be addressed; E-Mail: [email protected]; Tel.: +39-097-120-5885; Fax: +39-097-120-5896. Received: 30 December 2013; in revised form: 27 January 2014 / Accepted: 7 March 2014 / Published: 24 March 2014 Abstract: We present the state of the art on the modern mathematical methods of exploiting the entropy principle in thermomechanics of continuous media. A survey of recent results and conceptual discussions of this topic in some well-known non-equilibrium theories (Classical irreversible thermodynamics CIT, Rational thermodynamics RT, Thermodynamics of irreversible processes TIP, Extended irreversible thermodynamics EIT, Rational Extended thermodynamics RET) is also summarized. Keywords: continuum thermodynamics; entropy principle; rational thermodynamics; classical irreversible thermodynamics; extended irreversible thermodynamics; rational extended thermodynamics; exploitation of second law PACS Classification: 03.50.-z, 05.70.Ln, 05.70.Ce Entropy 2014, 16 1757 1.
    [Show full text]
  • Physical Origins of Anelasticity and Attenuation in Rock I
    2.17 Properties of Rocks and Minerals – Physical Origins of Anelasticity and Attenuation in Rock I. Jackson, Australian National University, Canberra, ACT, Australia ª 2007 Elsevier B.V. All rights reserved. 2.17.1 Introduction 496 2.17.2 Theoretical Background 496 2.17.2.1 Phenomenological Description of Viscoelasticity 496 2.17.2.2 Intragranular Processes of Viscoelastic Relaxation 499 2.17.2.2.1 Stress-induced rearrangement of point defects 499 2.17.2.2.2 Stress-induced motion of dislocations 500 2.17.2.3 Intergranular Relaxation Processes 505 2.17.2.3.1 Effects of elastic and thermoelastic heterogeneity in polycrystals and composites 505 2.17.2.3.2 Grain-boundary sliding 506 2.17.2.4 Relaxation Mechanisms Associated with Phase Transformations 508 2.17.2.4.1 Stress-induced variation of the proportions of coexisting phases 508 2.17.2.4.2 Stress-induced migration of transformational twin boundaries 509 2.17.2.5 Anelastic Relaxation Associated with Stress-Induced Fluid Flow 509 2.17.3 Insights from Laboratory Studies of Geological Materials 512 2.17.3.1 Dislocation Relaxation 512 2.17.3.1.1 Linearity and recoverability 512 2.17.3.1.2 Laboratory measurements on single crystals and coarse-grained rocks 512 2.17.3.2 Stress-Induced Migration of Transformational Twin Boundaries in Ferroelastic Perovskites 513 2.17.3.3 Grain-Boundary Relaxation Processes 513 2.17.3.3.1 Grain-boundary migration in b.c.c. and f.c.c. Fe? 513 2.17.3.3.2 Grain-boundary sliding 514 2.17.3.4 Viscoelastic Relaxation in Cracked and Water-Saturated Crystalline Rocks
    [Show full text]
  • An Ultrasonic Study on Anelasticity in Metals
    AN ULTRASONIC STUDY ON ANELASTICITY IN METALS J. P. Fulton, S. Nath, and B. Wincheski Analytical Services and Materials, Inc. 107 Research Drive Hampton, VA 23666 M. Namkung Mail Stop 231 NASA, Langley Research Center Hampton, VA 23681 INTRODUCTION Ultrasonic waves are highly sensitive to microstructural variations in materials and have been used extensively to investigate anharmonic effects in various metals and alloys[1-3]. A major focus of these studies is on the higher order elastic constants and their relation to the microstructure of the material. Ultrasonic techniques have also proven quite useful for char­ acterizing the stress state of a material [4-6]. Recently, while using the magnetoacoustic (MAC) method to investigate the residual stress in various steel samples, a time dependent change in the results was observed. It became apparent that the measurements were exhibit­ ing anelastic effects due to some intrinsic properties of the samples. Anelasticity is a phenomena that can occur whenever there is a rapid change in the exter­ nal stress applied to a crystalline solid. The change in stress will cause the induced strain to adopt a new equilibrium position. If the changes in the resulting strain does not take place instantaneously, the material is said to exhibit anelastic behavior. Investigations into anelas­ tic relaxation in crystalline solids have been carried out for many years and provide impor­ tant information relating to the microstructure of the material [7-13]. The time dependence of our MAC test results led us to propose a new approach for studying the causes of anelasticity in metals. Typical studies on anelasticity use strain gauges to monitor the recovery process.
    [Show full text]
  • The Mechanics of Ice
    COLD REGIONS SCIENCE AND ENGINEERING Monograph ll-C2b THE MECHANICS OF ICE John W. Glen December 1975 GB 2401 CORPS OF ENGINEERS, U.S. ARMY .U58m COLD REGIONS RESEARCH AND ENGINEERING LABORATORY no.ll-C2b HANOVER, NEW HAMPSHIRE 1975 Approved for public release; distribution unlimited. 5 fl rolomis The findings in this report are not to be construed as an official Department of the Army position unless so designated by other authorized docume* s. 0 Unclassified BUREAU OF RECLAMATION pENVER UBRARY C/ 92099625 ^ Y REPORT DOCUMENTATION PAGE ...............Q ? n Q 9 f i ? 5 1. REPORT NUMBER 2. GOVT ACCESSION NO. 3. RECIPt_______________ - ..------------ Monograph II-C2b -> 5. TYPE OF REPORT ft PERIOD COVERED 4. TITLE (end Subtitle) j m MECHANICS OF ICE 6. PERFORMING ORG. REPORT NUMBER 8. CONTRACT OR GRANT NUMBER*» 7. AUTHORf» European Research Office r J.W. Glen ,r Contract DAJA37-68-C*0208 £ ■ '« • ' 1 "• - TO. PROGRAM ELEMENT, PROJECT, TASK 9. PERFORMING ORGANIZATION NAME AND ADDRESS AREA ft WORK UNIT NJUMBERS Dr. John W. Glen Department of Physics ^ DA Project 1T062112A130 University of Birmingham ( ^ Task 01 ____ Birmingham. England ____; '____ ; _______ — 11. CONTROLLING OFFICE NAME AND ADDRESS 12. REPORT DATE Office, Chief of Engineers ^ December 1975 u 13. n u m b er o f p a g e s ' Washington, D.C. 47 14. MONITORING AGENCY NAME ft ADDRESS*?/ different from Controlling Office) 15. SECURITY CLASS, (of thie report) UvS ^Army Cold Regions Research and Engineering Laboratory-^ Unclassified Hanover, New Hampshire 03755 15«. DECLASSIFICATION/ DOWNGRADING SCHEDULE 16. DISTRIBUTION STATEMENT (of thla Report) Approved for public release; distribution unlimited.
    [Show full text]
  • The Anelasticity of the Mantle*
    CORE Metadata, citation and similar papers at core.ac.uk ProvidedGeophys. by Caltech i. AuthorsR. astr. - Main Soc. (1967) 14, 135-164. The Anelasticity of the Mantle* Don L. Anderson Summary The attenuation of seismic waves provides the most direct data regarding the non-elastic properties of the Earth. Recent experimental results from body waves, surface waves and free oscillations provide estimates of the anelasticity in various regions of the Earth. Results to date show that the upper mantle is more attenuating than the lower mantle, the maximum attenuation is in the vicinity of the low-velocity zone, a rapid increase in attenuation occurs in the vicinity of the C-region of the mantle and com­ pressional waves are less attenuated than shear waves. A frequency dependence of Q has not yet been discovered. Most laboratory measurements of attenuation have been performed at ultrasonic frequencies on pure specimens of metals, glasses, plastics and ceramics. A general feature of laboratory measurements is an exponential increase of attenuation with temperature on which are superimposed peaks which can be attributed to dislocation or other defect phenomena. Measurements on natural rocks at atmospheric pressure can be attributed to the presence of cracks. The intrinsic attenuation of rocks as a function of temperature and pressure is not known. However, on other materials grain boundary phenomena dominate at high temperature. This can be attributed to increased grain boundary mobility at high temperatures. High pressure would be expected to decrease this mobility. If attenuation in the mantle is due to an activated process it is probably controlled by the diffusion rate of defects at grain boundaries.
    [Show full text]
  • Anelasticity and Attenuation Part 1: Generalities Attenuation a Range Of
    Anelasticity and Attenuation Part 1: Generalities SIO 227C Bernard Minster April 2004 1 2 Attenuation A range of rheologies . Energy is lost to heat . Elasticity . Not explained by theory of elasticity . Inelasticity . Many mechanisms, most of which operate at . Plasticity the microscopic level . Viscosity . “slight” departure from elastic behavior . Viscoelasticity . Rheology based on simple considerations of . Visco-plasticity near-equilibrium thermodynamics. Elasto-plasticity . As much as possible, stick to linear theory . Anelasticity 3 4 1 Elasticity Viscosity . Linear elasticity (Green) . Linear viscosity Stress Stress . Hooke’s law σ = Kε σ = κε˙ . Nonlinear elasticity . Nonlinear viscosity Strain . Example. Strain rate σ = F(ε) σ ∝ε˙1 / N Stress Stress . Fully reversible deformation: . For instance, for water ice,N=3 Loading and unloading paths . Perhaps for the mantle as well are the same . Large stress leads to low effective . Equilibrium reached instantly Strain viscosity Strain rate . Deformation nonrecoverable 5 6 Plasticity Inelasticity . Instantaneous plasticity Stress . Nonlinear, . Rigid-plastic nonrecoverable . Time-dependent Strain . (example: material Stress Stress failure) . Elasto-plasticity Strain Strain rate 7 8 2 Anelasticity postulates Some Comments 1. For every stress there is a unique . Anelasticity is just like elasticity, plus the equilibrium value of strain and vice time dependence postulate versa . There can be an elastic component of the deformation in addition to the time-dependent 2. The equilibrium response is achieved component only after the passage of sufficient . Recovery is also time dependent time . Linearity is taken in the mathematical sense: 3. The stress-strain relationship is linear σ(αε1 + βε2) = ασ(ε1) + βσ(ε2) 9 10 Comparison of Rheologies Thermodynamic Substance Unique equilibrium Instan- Linear .
    [Show full text]
  • Anelastic Behavior of Small Dimensioned Aluminum
    metals Article Anelastic Behavior of Small Dimensioned Aluminum Enrico Gianfranco Campari 1 , Stefano Amadori 1, Ennio Bonetti 1, Raffaele Berti 1 and Roberto Montanari 2,* 1 Department of Physics and Astronomy, Bologna University, Viale Berti Pichat 6/2, I-40127 Bologna, Italy; [email protected] (E.G.C.); [email protected] (S.A.); [email protected] (E.B.); raff[email protected] (R.B.) 2 Department of Industrial Engineering, Rome University “Tor Vergata”, Via del Politecnico, 1-00133 Roma, Italy * Correspondence: [email protected]; Tel.: +39-0672597182 Received: 18 March 2019; Accepted: 9 May 2019; Published: 11 May 2019 Abstract: In the present research, results are presented regarding the anelasticity of 99.999% pure aluminum thin films, either deposited on silica substrates or as free-standing sheets obtained by cold rolling. Mechanical Spectroscopy (MS) tests, namely measurements of dynamic modulus and damping vs. temperature, were performed using a vibrating reed analyzer under vacuum. The damping vs. temperature curves of deposited films exhibit two peaks which tend to merge into a single peak as the specimen thickness increases above 0.2 µm. The thermally activated anelastic relaxation processes observed on free-standing films are strongly dependent on film thickness, and below a critical value of about 20 µm two anelastic relaxation peaks can be observed; both their activation energy and relaxation strength are affected by film thickness. These results, together with those observed on bulk specimens, are indicative of specific dislocation and grain boundary dynamics, constrained by the critical values of the ratio of film thickness to grain size.
    [Show full text]
  • Inference of Thermodynamic State in the Asthenosphere from Anelastic Properties, with Applications to North American Upper Mantle
    EarthArXiv Cover Page 2020-12-11 Inference of thermodynamic state in the asthenosphere from anelastic properties, with applications to North American upper mantle Christopher Havlin University of Illinois at Urbana-Champaign, havlin!illinois"ed$ Ben&amin Holtzman (amont-Doherty Earth *bservatory, Col$mbia University, benh!l#eo.col$mbia"ed$ Emily Hopper (amont-Doherty Earth *bservatory, Col$mbia University +his man$s ript is a non-peer revie,e# pre-print" It has been s$bmitte# for p$bli ation in Physics of the Earth and Planetary Interiors an# is c$rrently in revie," As s$ h, the content of the present man$s ript may change until the final version is p$blishe#, at whi h point a )*I lin. will dire t you to the final p$blishe# version" 1 Inference of thermodynamic state in the asthenosphere 2 from anelastic properties, with applications to North 3 American upper mantle 4 Christopher Havlin, Benjamin K. Holtzman*, Emily Hopper 5 Abstract 6 Inference of thermodynamic state and full-spectrum mechanical behavior of the litho- 7 sphere and asthenosphere is a central problem in geophysics, implicating our understand- 8 ing of the convection patterns, transient responses and chemical composition of the planet. 9 Anelasticity is responsible for significant relaxation of stress associated with seismic wave 10 propagation in the asthenosphere, while irreversible transient creep may be important in the 11 lithosphere. This paper focuses on the processes that may act at the time scales of seismic 12 wave propagation, and current questions in the effort to determine the dependence of these 13 effects on thermodynamic state.
    [Show full text]
  • Low-Frequency Measurements of Seismic Velocity and Attenuation in Antigorite Serpentinite Emmanuel C
    Low-Frequency Measurements of Seismic Velocity and Attenuation in Antigorite Serpentinite Emmanuel C. David1, Nicolas Brantut1, Lars N. Hansen2, and Ian Jackson3 Laboratory measurements of seismic velocity and atten- Reynard, 2013; David et al., 2018]. In addition, seismic at- uation in antigorite serpentinite at a confining pressure of tenuation (particularly shear attenuation) is very sensitive 2 kbar and temperatures up to 550 ◦C(i.e., in the antig- to many factors, including temperature, grain size, dislo- orite stability field) provide new results relevant to the in- cation density, melt fraction and redox conditions [see e.g. Jackson, 2015; Cline et al., 2018]. terpretation of geophysical data in subduction zones. A Among the serpentine group, antigorite is the mineral sta- polycrystalline antigorite specimen was tested via forced- ble over the widest depth range [Ulmer and Trommsdorff , oscillations at small strain amplitudes and seismic frequen- 1995]. The only existing, published attenuation data in ser- cies (mHz{Hz). The shear modulus has a temperature sensi- pentinites were reported on an anisotropic antigorite serpen- tivity, @G=@T , averaging −0:017 GPa K−1. Increasing tem- tinite by Kern et al. [1997], who measured the directional perature above 500 ◦C results in more intensive shear at- dependence of ultrasonic P- and S-wave velocities and their −1 associated attenuations up to 6 kbar confining pressure and tenuation (QG ) and associated modulus dispersion, with ◦ −1 700 C. More recently, Svitek et al. [2017] performed similar QG increasing monotonically with increasing oscillation pe- measurements for P-waves up to 4 kbar confining pressure. riod and temperature. This \background" relaxation is ade- However, such attenuation (and velocity) data were only ob- quately captured by a Burgers model for viscoelasticity and tained at the single MHz frequency of the ultrasonic pulse possibly results from intergranular mechanisms.
    [Show full text]
  • General Linear Response Analysis of Anelasticity
    Pram~na, Vol. 11, No. 4, October 1978, pp. 379-388, © printed in India General linear response analysis of anelasticity V BALAKRISHNAN Reactor Research Centre, Kalpakkam 603 102 MS received 1 May 1978 Abstract. Linear response theory is used to express the anelastic response (creep function and generalized compliance) of a system under an applied stress, in terms ot the equilibrium strain auto-correlation. These results extend an earlier analysis to cover inhomogeneous stresses and the tensor nature of the variables. For anelasti- city due to point defects, we express the strain compactly in terms of the elastic dipole tensor and the probability matrix governing dipole re-orientation and migration. We verify that re-orientations contribute to the deviatoric strain alone (Snoek, Zener, etc. effects), while the dilatory part arises solely from the long-range diffusion of the defects under a stress gradient (the Gorsky effect). Our formulas apply for arbitrary orien- tational multiplicity, specimen geometry, and stress inhomogeneity. The subse- quent development of the theory in any given situation then reduces to the modelling of the probability matrix referred to. In a companion paper, we apply our formalism to work out in detail the theory of the Gorsky effect (anelasticity due to long-range diffusion) for low interstitial concentrations, as an illustration of the advantages of our approach to the problem of anelastic relaxation. Keywords. Linear response theory; anelastic relaxation; compliance; elastic dipole; re-orientation; diffusion. 1. Introduction In an earlier paper (Balakrishnan et al 1978), referred to as BD¥ hereafter, a forma- lism was developed for mechanical response (specifically, anelasticity) based on the application of linear response theory (LRT).
    [Show full text]
  • Anelastic and Electromechanical Properties of Doped And&
    PROGRESS REPORT Ceria www.advmat.de Anelastic and Electromechanical Properties of Doped and Reduced Ceria Ellen Wachtel, Anatoly I. Frenkel, and Igor Lubomirsky* activity of ceria: the temporal, thermal, Room-temperature mechanical properties of thin films and ceramics of doped and compositional influence on the mate- and undoped ceria are reviewed with an emphasis on the anelastic behavior rial elasticity was little examined. How- of the material. Notably, the unrelaxed Young’s modulus of Gd-doped ceria ever, during the last decade, considerable ceramics measured by ultrasonic pulse-echo techniques is >200 GPa, while interest has developed in doped ceria as an essential component in micrometer-sized the relaxed biaxial modulus, calculated from the stress/strain ratio of thin fuel cells,[6] sensors, and even electrome- films, is≈ 10 times smaller. Oxygen-deficient ceria exhibits a number of chanical actuators.[7] These scaled-down anelastic effects, such as hysteresis of the lattice parameter, strain-dependent devices demand particularly precise Poisson’s ratio, room-temperature creep, and nonclassical electrostriction. mechanical engineering, e.g., because Methods of measuring these properties are discussed, as well as the thermal cycling of ceria films during device operation can be associated with applicability of Raman spectroscopy for evaluating strain in thin films of the generation of unpredictable mechan- Gd-doped ceria. Special attention is paid to detection of the time dependence ical strain. To give an example, the linear of anelastic effects. Both the practical advantages and disadvantages of thermal expansion coefficient of undoped anelasticity on the design and stability of microscopic devices dependent ceria is ≈11 × 10−6 K−1, which predicts that on ceria thin films are discussed, and methods of mitigating the latter are heating from 300 to 500 K will generate suggested, with the aim of providing a cautionary note for materials scientists strain 0.22%.
    [Show full text]