A simple model for the small-strain behaviour of soils
|
|
- Trevor Park
- 6 years ago
- Views:
Transcription
1 A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo , São Paulo, Brazil, Keywors: constitutive moels; elasticity; soil mechanics; soil eformation. Abstract: The constitutive equation iscusse in this note eliminates some efects of linear elasticity in the escription of the small-strain behaviour of soils. It is capable of representing volume changes in pure shear an ifferent values of bulk moulus in compression an expansion. The new equation provies a simple escription of soil behaviour at small strains in that it prouces linear stress-strain relations that can approximate the initial part of experimental stress-strain curves. 1. Introuction Linear elasticity is not a goo moel for soils even at the small-strain level. onsier, for instance, ientical isotropic soil samples subject to ifferent stress paths in the triaxial cell. Although the beginning of each particular stress-strain curve can be approximate by a tangent line through the origin, experience shows that ifferent pairs of elastic parameters (say, bulk an shear moulus) woul be necessary to represent the small-strain response in the various stress paths. The simplest example of this fact occurs in the isotropic compression or expansion of a normally consoliate sample: the bulk moulus for a subsequent virgin loaing is smaller than the bulk moulus for unloaing. Another shortcoming of the linear elastic moel is that it cannot preict volume changes in pure shear (e.g., volume increase in ense sans, volume ecrease in loose sans) as well as values of Skempton s A parameter ifferent from 1/3 in the unraine compression of a saturate soil (see, e.g., Lambe an Whitman 1979, for experimental results). This note presents a constitutive moel for the small-strain behaviour of soils that, unlike linear elasticity, is able to preict the occurrence of volume changes in pure shear, the existence of ifferent values of bulk moulus in isotropic compression an expansion an ifferent values of Young s moulus an Poisson s ratio in axial compression an axial extension. On the other han, like linear elasticity, the
2 propose moel yiels linear stress-strain relations in straight stress-paths an so provies a relatively simple approximate representation of the behaviour at small strains. Throughout the text, we use the soil mechanics sign convention for stresses an strains (compressive stresses an strains are positive etc). 2. The new moel As a starting point, consier the isotropic linear elastic equation, as use in soil mechanics: T = KI ( + G. (1) 1 ) I 2 T, an I are, respectively, the effective stress increment tensor, the strain tensor an the ientity tensor; is the eviatoric part of, I ( ) is the first invariant of, K is the bulk moulus an G is the 1 shear moulus. Soils are always subject to initial effective stresses, whether in the fiel or in laboratory. q. 1 thus yiels the increment T that shoul be ae to the initial effective stress tensor to obtain the final effective stress tensor (applications of linear elasticity to soil mechanics problems is a classical subject; see, e.g., Terzaghi 1943; Davis an Selvaurai 1996). Now the first term on the right-han sie of q. 1 will be replace by a more general function of to prouce the new moel. The resulting equation reas T = α( ) I + 2G, (2) where 1 ( ) / 2 + ( K K ) I1( ) / 2 + ξg I 2 ( ) α ( ) = ( K + K ) I. I ) enotes the secon 2 ( invariant of. The parameters K (compression bulk moulus), K (expansion bulk moulus) an G (shear moulus) are positive whereas ξ (ilatancy coefficient) may be negative, positive or null. Their physical meaning will become clear in the examples below. Note that, for I ( ) 0, 1 K 1 ( ) + ξg I 2 ( ) α ( ) = I while, for I ( ) 0, 1 K 1 ( ) + ξg I 2 ( ) α ( ) = I. Of course, q. 2 reuces to q. 1 if K = K an ξ=0. In the next section the behaviour of the propose equation will be investigate in typical situations (isotropic compression an expansion, pure shear, axial compression an extension, unraine axial compression). For the analysis of cases in which T is given an must be foun, it is convenient to invert the stress-strain relation expresse by q. 2. As = I ( +, it suffices to write I ( ) an 1 ) I / 3 1
3 in terms of T. learly, = / 2G (in which T is the eviator stress increment tensor). T Further, it can be conclue that, if I ( ) 0, then I ( ) ) / K while, if I ( ) 0, K 1 1 T 1 I ( ) ) /, where ρ( T) = I ( T) / 3 ξ I ( ) / 2. Hence, the following equation hols: 1 T 1 2 T 1 = β( T) I + T, (3) 2G in which β( T) = (1/ 6K + 1/ 6K ) ρ( T) + (1/ 6K 1/ 6K ) ρ( T). Therefore, for ρ( T ) 0, β ( T) T) / 3K, whereas, for ρ( T ) 0, β ( T) T) / 3K. In any case, the volumetric strain is given by ε = I ( ) = 3β( ) v 1 T. 3. Moel behaviour in simple cases a. Isotropic compression an expansion The response of q. 2 to an isotropic strain ( = εi ) is an isotropic stress increment ( T = σi ). The relation between the volumetric strain ( ε v = 3ε ) an σ epens on the sign of ε. It is σ = K εv in case of compression ( ε > 0 ) an σ = K ε v in case of expansion ( ε < 0 ). Thus, if K K, the moel preicts a ifferent linear stress-strain relation in each case. b. Pure shear For a pure shear efine by the stress increment matrix [ T ] = 0 0 τ (4) 0 τ 0 ( τ > 0 ), q. 3 gives the following strain matrix: ε 0 0 [] = 0 ε γ / 2. (5) 0 γ / 2 ε
4 Thus the shear strain γ is relate to the shear stress τ through γ = τ / G, as in linear elasticity. By introucing ρ( T ) = ξτ / 2 in the expression for β ( T), the relation between the volumetric strain ( ε v = 3ε ) an τ can be obtaine. If ξ < 0, then ε v = ξτ / 2K, a positive number, inicating volume ecrease (contraction). On the other han, if ξ > 0, then ε v = ξτ / 2K, which is negative, inicating volume increase (expansion). Therefore, unlike linear elasticity, the present moel preicts volume changes in pure shear. c. Axial compression an extension Now the only non-zero stress increment component is, say, T = σ (axial compression: σ > 0 ; axial 33 extension: σ < 0 ). Accoring to q. 3, the corresponing strain matrix is iagonal, with 11 = 22. In compression ρ( T ) = mσ, with m = 1 3ξ / 2 3, whereas in extension ρ( T ) = nσ, with n = 1 + 3ξ / 2 3. Table 1 gives the obtaine expressions for the Young s moulus ( σ / 33 ) in compression ( ) an in extension ( ) in terms of the basic parameters, G, K, K an ξ. The corresponing values of the Poisson s ratio ( 22 / 33 ) in compression ( ν ) an in extension ( ν ) are given in all cases by ν = / 2G 1 an ν = / 2G 1. Note that there are two ifferent expressions to compute both an, each expression associate to an interval of values of ξ. In aition, Table 1 gives expressions for the volumetric strain an, in the last column, restrictions that must be impose, if any, on the basic parameters in orer that the values of Young s moulus an Poisson s ratio are positive, as expecte on physical grouns. The meaning of the symbols H, L, L, H appearing in the first column of Table 1will now be explaine. They constitute a classification base on the value of the ilatancy coefficient that helps unerstan the behaviour of the material efine by q. 2 in regar to volume changes. The symbols an stan for two main material classes relate to volume changes in pure shear. lass ( ξ < 0 ): materials that contract in pure shear. lass ( ξ 0 ): materials that expan ( ξ > 0 ) or suffer no volume change ( ξ = 0 ) in pure shear. lasses an are further ivie into two subclasses accoring to the behaviour in axial compression an extension. Subclass H ( ξ < 2 3 / 3; high contractibility): materials in class that contract both in axial compression an extension. Subclass L ( 2 3 / 3 ξ < 0 ; low contractibility;): materials in class that contract in axial compression an expan ( 2 3 / 3 < ξ ) or suffer no volume
5 change ( 2 3 / 3 = ξ ) in axial extension. Subclass L ( 0 ξ 2 3 / 3: low expansibility): materials in class that expan in axial extension an contract ( 0 ξ < 2 3 / 3) or suffer no volume change ( ξ = 2 3 / 3 ) in axial compression. Subclass H ( ξ > 2 3 / 3 ; high expansibility): materials in class that expan both in axial compression an extension. 1. The appropriate expressions for, an ε v in each subclass can be obtaine irectly from Table Table 1: Young s moulus an volumetric strain in axial compression an extension. Axial compression ( σ > 0 ; m = 1 3ξ / 2 3 ) ξ 2 3 / 3 ( m 0 ) H, L, L = 3 K + mg ε v = mσ / 3K (7) (6) ( ε v 0 ) > 0. ν > 0, if ξ > 2 3 / 3 3K / G. ξ > 2 3 / 3 ( m < 0 ) H = 3 K + mg ε v = mσ / 3K (9) (8) ( ε v < 0 ) > 0 an > 0, if ν ξ < 2 3G / K. Axial extension ( σ < 0 ; n = 1 + 3ξ / 2 3 ) ξ < 2 3 / 3 ( n < 0 ) H = 3 K + ng ε v = nσ / 3K (11) (10) ( ε v > 0 ) > 0 an > 0, if ν ξ > 2 3G / 3 2 3K. ξ 2 ( n 0 ) 3 / 3 L, L, H = 3 K + ng ε v = nσ / 3K (13) (12) ( ε v 0 ) > 0. ν > 0, if ξ < 2 3 / 3 + 3K / G.
6 . Unraine axial compression Here the strain matrix reas: ε / [] = 0 ε / 2 0, (14) 0 0 ε where ε > 0 is the axial compressive strain. From q. 2, the non-zero stress increment components rea (recall they are effective stresses): 3 T 11 = T22 = (1 ξ) Gε, (15) 2 3 T 33 = (2 + ξ) Gε. (16) 2 In the traitional unraine compression test the total lateral stress remains constant. Hence the pore pressure increment u is the negative of the lateral effective stress increment: 3 u = ( 1 ξ) Gε. (17) 2 It is positive if ξ < 2 3 / 3 (H, L, L). Only soils in the class H generate negative pore pressures. The corresponing Skempton s pore pressure parameter is A = 1/ 3 3ξ / Final remarks learly, avance constitutive moels can correct the mentione efects of linear elasticity an, in aition, are able to escribe soil behaviour in complex stress paths (for funamentals of elastoplastic moels see Davis e Selvaurai 2005; as for hypoplastic moels, Kolymbas 2000, Naer 2003, 2010). But, ue to its simplicity, the propose moel may be avantageous if we are intereste only in moelling small-strain behaviour in monotonic short stress-paths, far from failure. Besies, the moel may have a iactic value: it can be an intermeiate stage between linear elasticity an more complex theories.
7 References Davis, R.O. an Selvaurai, A.P.S. (1996). lasticity an Geomechanics. ambrige University Press. Davis, R.O. an Selvaurai, A.P.S. (2005). Plasticity an Geomechanics. ambrige University Press. Kolymbas, D. (2000). Introuction to Hypoplasticity. Balkema. Lambe, T.W. an Whitman, R.V. (1979). Soil Mechanics. Wiley. Naer, J.J. (2003). Analysis of a new hypoplastic equation. International Journal of Plasticity 19, Naer, J.J. (2010). A new constitutive theory extening Hypoplasticity. Mechanics Research ommunications 37, Terzaghi, K. (1943). Theoretical Soil Mechanics. Wiley.
EVALUATION OF LIQUEFACTION RESISTANCE AND LIQUEFACTION INDUCED SETTLEMENT FOR RECLAIMED SOIL
386 EVALUATION OF LIQUEFACTION RESISTANCE AND LIQUEFACTION INDUCED SETTLEMENT FOR RECLAIMED SOIL Lien-Kwei CHIEN 1, Yan-Nam OH 2 An Chih-Hsin CHANG 3 SUMMARY In this stuy, the fille material in Yun-Lin
More informationDISCRETE MODELING OF STRAIN ACCUMULATION IN GRANULAR SOILS UNDER CYCLIC LOADING
11th Worl Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Flui Dynamics (ECFD VI) E. Oñate, J. Oliver
More informationA note on the Mooney-Rivlin material model
A note on the Mooney-Rivlin material moel I-Shih Liu Instituto e Matemática Universiae Feeral o Rio e Janeiro 2945-97, Rio e Janeiro, Brasil Abstract In finite elasticity, the Mooney-Rivlin material moel
More informationStrength Analysis of CFRP Composite Material Considering Multiple Fracture Modes
5--XXXX Strength Analysis of CFRP Composite Material Consiering Multiple Fracture Moes Author, co-author (Do NOT enter this information. It will be pulle from participant tab in MyTechZone) Affiliation
More informationAn Anisotropic Hardening Model for Springback Prediction
An Anisotropic Harening Moel for Springback Preiction Danielle Zeng an Z. Ceric Xia Scientific Research Laboratories For Motor Company Dearborn, MI 48 Abstract. As more Avance High-Strength Steels (AHSS
More informationA PROCEDURE FOR DETERMINATION OF THE ALTERNAT MODEL PARAMETERS
1. Mohamme Y. FATTAH,. Omar. F. S. al DAMLUJI, 3. Yousif. J. al SHAKARCHI A PROCEDURE FOR DETERMINATION OF THE ALTERNAT MODEL PARAMETERS 1. BUILDING AND CONSTRUCTION DEPARTMENT, UNIVERSITY OF TECHNOLOGY,
More informationVUMAT for Fabric Reinforced Composites
VUMAT or Fabric Reinorce Composites. Introuction This ocument escribes a constitutive mo or abric reinorce composites that was introuce in Abaqus/Exicit 6.8. The mo has been imemente as a built-in VUMAT
More informationThe Phenomenon of Anomalous Rock Embrittlement
https://papers.acg.uwa.eu.au/p/574_29_tarasov/ B.G. Tarasov A.V. Dyskin School of Civil an Resource Engineering The University of Western Australia The paper analyses a phenomenon of rock behaviour - the
More informationResilient Modulus Prediction Model for Fine-Grained Soils in Ohio: Preliminary Study
Resilient Moulus Preiction Moel for Fine-Graine Soils in Ohio: Preliminary Stuy by Teruhisa Masaa: Associate Professor, Civil Engineering Department Ohio University, Athens, OH 4570 Tel: (740) 59-474 Fax:
More informationAdhesive Wear Theory of Micromechanical Surface Contact
International Journal Of Computational Engineering esearch ijceronline.com Vol. Issue. hesive Wear Theory of Micromechanical Surface Contact iswajit era Department of Mechanical Engineering National Institute
More informationState-Space Model for a Multi-Machine System
State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal
More informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More informationSemiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom
PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationCrack-tip stress evaluation of multi-scale Griffith crack subjected to
Crack-tip stress evaluation of multi-scale Griffith crack subjecte to tensile loaing by using periynamics Xiao-Wei Jiang, Hai Wang* School of Aeronautics an Astronautics, Shanghai Jiao Tong University,
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationChaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena
Chaos, Solitons an Fractals (7 64 73 Contents lists available at ScienceDirect Chaos, Solitons an Fractals onlinear Science, an onequilibrium an Complex Phenomena journal homepage: www.elsevier.com/locate/chaos
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More informationSurvey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013
Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing
More information3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects
IOP Conference Series: Materials Science an Engineering 3-D FEM Moeling of fiber/matrix interface eboning in UD composites incluing surface effects To cite this article: A Pupurs an J Varna 2012 IOP Conf.
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationExperimental Investigation on the Dynamic Shear Modulus and Damping Ratio of Aeolian Soils
2016 International Conference on Power Engineering & Energy, Environment (PEEE 2016) ISBN: 978-1-60595-376-2 Experimental Investigation on the Dynamic Shear Moulus an Damping Ratio of Aeolian Soils Shengyang
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson
JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More informationLQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR
5 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR Tatsunori Kaneko* an Yasuto Asano* * Department of Mechanical Engineering,
More informationAn Iterative Incremental Learning Algorithm for Complex-Valued Hopfield Associative Memory
An Iterative Incremental Learning Algorithm for Complex-Value Hopfiel Associative Memory aoki Masuyama (B) an Chu Kiong Loo Faculty of Computer Science an Information Technology, University of Malaya,
More informationTHE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE
Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek
More informationMULTISCALE FRICTION MODELING FOR SHEET METAL FORMING
MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING Authors J. HOL 1, M.V. CID ALFARO 2, M.B. DE ROOIJ 3 AND T. MEINDERS 4 1 Materials innovation institute (M2i) 2 Corus Research Centre 3 University of
More informationA Review of Multiple Try MCMC algorithms for Signal Processing
A Review of Multiple Try MCMC algorithms for Signal Processing Luca Martino Image Processing Lab., Universitat e València (Spain) Universia Carlos III e Mari, Leganes (Spain) Abstract Many applications
More informationInternational Conference on Advances in Energy, Environment and Chemical Engineering (AEECE-2015)
International Conference on Avances in Energy, Environment an Chemical Engineering (AEECE-2015) Stuy on Damage Characteristic of Unergroun Cavern Blasting Excavation Base on Dynamic Damage Constitutive
More information2Algebraic ONLINE PAGE PROOFS. foundations
Algebraic founations. Kick off with CAS. Algebraic skills.3 Pascal s triangle an binomial expansions.4 The binomial theorem.5 Sets of real numbers.6 Surs.7 Review . Kick off with CAS Playing lotto Using
More informationSYNTHESIS ON THE ASSESSMENT OF CHIPS CONTRACTION COEFFICIENT C d
SYNTHESIS ON THE ASSESSMENT OF HIPS ONTRATION OEFFIIENT Marius MILEA Technical University Gheorghe Asachi of Iași, Machine Manufacturing an Inustrial Management Faculty, B-ul D Mangeron, nr 59, A, 700050,
More informationLATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION
The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische
More informationFLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS. 1. Introduction
FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS ALINA BUCUR, CHANTAL DAVID, BROOKE FEIGON, MATILDE LALÍN 1 Introuction In this note, we stuy the fluctuations in the number
More informationPredictive Control of a Laboratory Time Delay Process Experiment
Print ISSN:3 6; Online ISSN: 367-5357 DOI:0478/itc-03-0005 Preictive Control of a aboratory ime Delay Process Experiment S Enev Key Wors: Moel preictive control; time elay process; experimental results
More informationAnalytic Scaling Formulas for Crossed Laser Acceleration in Vacuum
October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945
More informationHow the potentials in different gauges yield the same retarded electric and magnetic fields
How the potentials in ifferent gauges yiel the same retare electric an magnetic fiels José A. Heras a Departamento e Física, E. S. F. M., Instituto Politécnico Nacional, México D. F. México an Department
More information7. STRESS ANALYSIS AND STRESS PATHS
7-1 7. STRESS ANALYSIS AND STRESS PATHS 7.1 THE MOHR CIRCLE The discussions in Chapters and 5 were largely concerned with vertical stresses. A more detailed examination of soil behaviour requires a knowledge
More informationConservation laws a simple application to the telegraph equation
J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness
More informationApproximate reduction of dynamic systems
Systems & Control Letters 57 2008 538 545 www.elsevier.com/locate/sysconle Approximate reuction of ynamic systems Paulo Tabuaa a,, Aaron D. Ames b, Agung Julius c, George J. Pappas c a Department of Electrical
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More informationThe proper definition of the added mass for the water entry problem
The proper efinition of the ae mass for the water entry problem Leonaro Casetta lecasetta@ig.com.br Celso P. Pesce ceppesce@usp.br LIE&MO lui-structure Interaction an Offshore Mechanics Laboratory Mechanical
More informationIntroduction to Markov Processes
Introuction to Markov Processes Connexions moule m44014 Zzis law Gustav) Meglicki, Jr Office of the VP for Information Technology Iniana University RCS: Section-2.tex,v 1.24 2012/12/21 18:03:08 gustav
More informationFinite element analysis of the dynamic behavior of concrete-filled double skin steel tubes (CFDST) under lateral impact with fixed-simply support
Structures Uner Shock an Impact XIII 383 Finite element analysis of the ynamic behavior of concrete-fille ouble skin steel tubes (CFDST) uner lateral impact with fixe-simply support M. Liu & R. Wang Architecture
More informationTOWARDS THERMOELASTICITY OF FRACTAL MEDIA
ownloae By: [University of Illinois] At: 21:04 17 August 2007 Journal of Thermal Stresses, 30: 889 896, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0149-5739 print/1521-074x online OI: 10.1080/01495730701495618
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the
More informationPhysics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1
Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems
More informationDuctility and Failure Modes of Single Reinforced Concrete Columns. Hiromichi Yoshikawa 1 and Toshiaki Miyagi 2
Ductility an Failure Moes of Single Reinforce Concrete Columns Hiromichi Yoshikawa 1 an Toshiaki Miyagi 2 Key Wors: seismic capacity esign, reinforce concrete column, failure moes, eformational uctility,
More informationProblem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs
Problem Sheet 2: Eigenvalues an eigenvectors an their use in solving linear ODEs If you fin any typos/errors in this problem sheet please email jk28@icacuk The material in this problem sheet is not examinable
More informationDIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10
DIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10 5. Levi-Civita connection From now on we are intereste in connections on the tangent bunle T X of a Riemanninam manifol (X, g). Out main result will be a construction
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More information2-7. Fitting a Model to Data I. A Model of Direct Variation. Lesson. Mental Math
Lesson 2-7 Fitting a Moel to Data I BIG IDEA If you etermine from a particular set of ata that y varies irectly or inversely as, you can graph the ata to see what relationship is reasonable. Using that
More informationAssessment of the Buckling Behavior of Square Composite Plates with Circular Cutout Subjected to In-Plane Shear
Assessment of the Buckling Behavior of Square Composite Plates with Circular Cutout Sujecte to In-Plane Shear Husam Al Qalan 1)*, Hasan Katkhua 1) an Hazim Dwairi 1) 1) Assistant Professor, Civil Engineering
More informationLinear and quadratic approximation
Linear an quaratic approximation November 11, 2013 Definition: Suppose f is a function that is ifferentiable on an interval I containing the point a. The linear approximation to f at a is the linear function
More informationd dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1
Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of
More informationThe thin plate theory assumes the following deformational kinematics:
MEG6007 (Fall, 2017) - Solutions of Homework # 8 (ue on Tuesay, 29 November 2017) 8.0 Variational Derivation of Thin Plate Theory The thin plate theory assumes the following eformational kinematics: u
More informationTime-of-Arrival Estimation in Non-Line-Of-Sight Environments
2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor
More informationMATHEMATICS BONUS FILES for faculty and students
MATHMATI BONU FIL for faculty an stuents http://www.onu.eu/~mcaragiu1/bonus_files.html RIVD: May 15, 9 PUBLIHD: May 5, 9 toffel 1 Maxwell s quations through the Major Vector Theorems Joshua toffel Department
More informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
More informationConservation Laws. Chapter Conservation of Energy
20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More informationA Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges
Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationPeriods of quadratic twists of elliptic curves
Perios of quaratic twists of elliptic curves Vivek Pal with an appenix by Amo Agashe Abstract In this paper we prove a relation between the perio of an elliptic curve an the perio of its real an imaginary
More informationMarkov Chains in Continuous Time
Chapter 23 Markov Chains in Continuous Time Previously we looke at Markov chains, where the transitions betweenstatesoccurreatspecifietime- steps. That it, we mae time (a continuous variable) avance in
More information1 Introuction In the past few years there has been renewe interest in the nerson impurity moel. This moel was originally propose by nerson [2], for a
Theory of the nerson impurity moel: The Schrieer{Wol transformation re{examine Stefan K. Kehrein 1 an nreas Mielke 2 Institut fur Theoretische Physik, uprecht{karls{universitat, D{69120 Heielberg, F..
More informationINTERNATIONAL CONFERENCE ON ADVANCES IN MATERIALS AND PROCESSING TECHNOLOGIES (AMPT2010), October 2010, Paris, (France)
Provie by the author(s) an University College Dublin Library in accorance with publisher policies. Please cite the publishe version when available. Title Extrusion process by finite volume metho using
More informationCalculus of Variations
Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,
More informationMomentum and Energy. Chapter Conservation Principles
Chapter 2 Momentum an Energy In this chapter we present some funamental results of continuum mechanics. The formulation is base on the principles of conservation of mass, momentum, angular momentum, an
More informationRastislav Ďuriš, Justín Murín. A Solution for 1D Non-linear Problems for Finite Elements With Full Stiffness Matrices
Rastislav Ďuriš, Justín Murín A Solution for D Non-linear Problems for Finite lements With Full Stiffness Matrices Scientific Monographs in Automation an Computer Science ite by Prof. Dr. Peter Husar Ilmenau
More informationA. Incorrect! The letter t does not appear in the expression of the given integral
AP Physics C - Problem Drill 1: The Funamental Theorem of Calculus Question No. 1 of 1 Instruction: (1) Rea the problem statement an answer choices carefully () Work the problems on paper as neee (3) Question
More informationPoroelasticity and tidal loading
Chapter 6: Poroelasticity an tial loaing 63 Chapter 6: Poroelasticity an tial loaing 6. Introuction It is establishe above, in Chapter 4, that there are many known examples of the tial moulation of seafloor
More information3.2 Shot peening - modeling 3 PROCEEDINGS
3.2 Shot peening - moeling 3 PROCEEDINGS Computer assiste coverage simulation François-Xavier Abaie a, b a FROHN, Germany, fx.abaie@frohn.com. b PEENING ACCESSORIES, Switzerlan, info@peening.ch Keywors:
More informationApproximate Reduction of Dynamical Systems
Proceeings of the 4th IEEE Conference on Decision & Control Manchester Gran Hyatt Hotel San Diego, CA, USA, December 3-, 6 FrIP.7 Approximate Reuction of Dynamical Systems Paulo Tabuaa, Aaron D. Ames,
More informationOn the dynamics of the electron: Introduction, 1, 9
On the ynamics of the electron: Introuction,, 9 Henri oincaré Introuction. It seems at first that the aberration of light an relate optical an electrical phenomena will provie us with a means of etermining
More informationNonPAS: A Program for Nonlinear Analysis of Flexible Pavements
International Journal of Integrate Engineering, Vol. 7 No. (05) pp. -8 NonPAS: A Program for Nonlinear Analysis of Flexible Pavements Ali eza Ghanizaeh,*, Arash Ziaie, Department of Civil Engineering,
More informationarxiv: v3 [quant-ph] 16 Feb 2016
What is the probability of a thermoynamical transition? arxiv:504.00020v3 [quant-ph] 6 Feb 206 Álvaro M. Alhambra,, Jonathan Oppenheim,, 2, an Christopher Perry, Department of Physics an Astronomy, University
More informationThe dielectric constant of Na 0.4 K 0.6 Br and its large temperature variation
The ielectric constant of Na.4 K.6 Br an its large temperature variat V. Katsika-Tsigourakou*, P. Bekris Sect of Soli State Physics, Department of Physics, Natal an Kapoistrian University of Athens, Panepistimiopolis,
More informationDiophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations
Diophantine Approximations: Examining the Farey Process an its Metho on Proucing Best Approximations Kelly Bowen Introuction When a person hears the phrase irrational number, one oes not think of anything
More informationPDE Notes, Lecture #11
PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =
More information6. Friction and viscosity in gasses
IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner
More informationThe theoretical limits to the power output of a muscle tendon complex with inertial and gravitational loads
Receive 1 November Accepte 17 March 3 Publishe online 11 June 3 The theoretical limits to the power output of a muscle tenon complex with inertial an gravitational loas Apostolos Galantis * an Roger C.
More informationSYNCHRONOUS SEQUENTIAL CIRCUITS
CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents
More informationThe Press-Schechter mass function
The Press-Schechter mass function To state the obvious: It is important to relate our theories to what we can observe. We have looke at linear perturbation theory, an we have consiere a simple moel for
More informationThe total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More informationA USER-FRIENDLY TOOL FOR SIMULATING THE TIME-DEPENDENT FIELD DISTRIBUTION IN PE INSULATION ON THE BASIS OF A PHYSICAL APPROACH
A USER-FRIENLY TOOL FOR SIMULATING THE TIME-EPENENT FIEL ISTRIBUTION IN PE INSULATION ON THE BASIS OF A PHYSICAL APPROACH Séverine LE ROY, Laboratoire Plasma et Conversion Energie, (France), severine.leroy@laplace.univ-tlse.fr
More information18 EVEN MORE CALCULUS
8 EVEN MORE CALCULUS Chapter 8 Even More Calculus Objectives After stuing this chapter you shoul be able to ifferentiate an integrate basic trigonometric functions; unerstan how to calculate rates of change;
More informationarxiv:physics/ v4 [physics.class-ph] 9 Jul 1999
AIAA-99-2144 PROPULSION THROUGH ELECTROMAGNETIC SELF-SUSTAINED ACCELERATION arxiv:physics/9906059v4 [physics.class-ph] 9 Jul 1999 Abstract As is known the repulsion of the volume elements of an uniformly
More informationPRACTICE 4. CHARGING AND DISCHARGING A CAPACITOR
PRACTICE 4. CHARGING AND DISCHARGING A CAPACITOR. THE PARALLEL-PLATE CAPACITOR. The Parallel plate capacitor is a evice mae up by two conuctor parallel plates with total influence between them (the surface
More informationApplications of First Order Equations
Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more
More informationSmectic-C tilt under shear in smectic-a elastomers
Smectic-C tilt uner shear in smectic-a elastomers Olaf Stenull an T. C. Lubensky Department of Physics an Astronomy, University of Pennsylvania, Philaelphia, Pennsylvania 19104, USA J. M. Aams an Mark
More informationBoth the ASME B and the draft VDI/VDE 2617 have strengths and
Choosing Test Positions for Laser Tracker Evaluation an Future Stanars Development ala Muralikrishnan 1, Daniel Sawyer 1, Christopher lackburn 1, Steven Phillips 1, Craig Shakarji 1, E Morse 2, an Robert
More informationAcute sets in Euclidean spaces
Acute sets in Eucliean spaces Viktor Harangi April, 011 Abstract A finite set H in R is calle an acute set if any angle etermine by three points of H is acute. We examine the maximal carinality α() of
More informationNonlinear Dielectric Response of Periodic Composite Materials
onlinear Dielectric Response of Perioic Composite aterials A.G. KOLPAKOV 3, Bl.95, 9 th ovember str., ovosibirsk, 639 Russia the corresponing author e-mail: agk@neic.nsk.su, algk@ngs.ru A. K.TAGATSEV Ceramics
More information