A Fourier Transform 3D-Model for Wave Propagation in layered orthotropic Media
|
|
- Aubrey Harrell
- 5 years ago
- Views:
Transcription
1 A Fourier Transform 3D-Model for Wave Propagation in layered orthotropic Media F. Kramer*, H. Waubke Abstract In order to assess and prevent noise emmission in the environment of traffic routes, prognosis models for wave propagation are of interest for many applications. Guided by experience in railway traffic simulation, this contribution introduces a new model for wave propagation and simulation in layered orthotropic media. The model uses a component-wise Fourier transform to obtain the general solution within one layer; by a variation of constants, this yields the solution within a loaded layer. Then, balances at layer boundaries are used to derive the solution at the layer boundaries. Thus, a grid of stress respectively displacement vectors is obtained, and back transformation is performed. The model has been tested for physical accuracy concerning the splitting of layers and for numerical accuracy. The results are presented. INTRODUCTION The simulation of vibrations induced by machinery and traffic has become of increasing concern within the last decades, due to frequent construction of high speed trains and mass rapid transport systems. Theoretical studies on the propagation of waves in anisotropic layered media have been an important field of research for many years, and have been of particular use for applications in geophysics and seismology [1, 2]. Therefore, much research has been done on wave propagation in anisotropic layered media [3], layered orthotropic media however, has rarely been focused upon. Waubke [4] has suggested extensions and developments for this model in earlier papers. This paper draws these ideas together, and offers some concrete numerical results. In the first section, a model of wave propagation is derived by a Fourier transform of the classical wave equation. The balance of stresses and displacements at the layer boundary delivers a system of linear equations, yielding the stress and displacement level at the layer boundaries as functions ûk x,ky,z,ω and ˆσk x,k y,z,ω.
2 This model applied to a system of layers gives a transport matrix for the layer boundaries, yielding a system whose solution permits the calculation of displacement and stress levels at arbitrary points within the media. Due to the numerical stabililty of the approach, this is possible for large distances from the point of excitation. At first, the focus lies on one layer and the equations describing the process are derived. The strain vector ǫ is a vector of local derivatives of the displacement vector u, denoted by ǫ := D x, y, z u. The soil is assumed to consist of n orthotropic layers. For each layer, the following parameters are needed: the shear modulus G, the Young s modulus E, the Poisson ratios ν, the density ρ and width d. Additionally, homogenous behavior in the horizontal directions is assumed: E y = E x, G yz = G zx, and ν yz = ν zx hold. So, the elasticity matrix takes block diagonal form with blocks E 11 := 1 E x νxy E x νxy E x 1 E x νzx νzx νzx νzx 1 1, E 22 := diag2g xy, 2G zx, 2G zx 1 As can be found in [1], the wave equation takes the special form Bσ + ρω 2 û = ˆBE ˆDû + ρω 2 û =: Âû = ˆb. 2 for a differential operator B, where b is some external force. We denote the Fourier transform ˆf of a function fα by F α,kα fk α, and the Fourier back transform by F 1 k α,α.the displacement vector u := ux,y,z,t is Fourier transformed in all directions, yielding ûk x,k y,k z,ω. Thus, the differential operator D is substituted by the polynomial operator ˆD. Simple equations are obtained, expressing the link of the transformed displacement vector and the transformed stress vector: ˆσ = Eˆǫ = E ˆDû, where ˆD = ˆDk x,k y,k z is the transformed differential operator D. Note that E ˆD is symmetric. The transformed wave equation then takes the form ˆBE ˆDû + ρω 2 û =: Âû = ˆb. 3 In order to cope with the boundary conditions, it is useful to implement a grid in k x,k y,ω. Therefore k x,k y,ω are treated as constants, and are omitted in some cases. The general form of a solution DERIVING THE EQUATIONS The condition detâ = 0 is necessary for solutions of Âû = 0, û 0 to exist. The roots of the polynomial detâk z called eigenvalues are denoted by κ, and the vectors of the corresponding kernel are denoted as Ψ called eigenvectors. Note that
3 ICSV13, July 2-6, 2006, Vienna, Austria due to the symmetry of Â, it holds that detâk z = detâ k z. The homogeneous solution u h and the corresponding stress function σ h take the form û h k z = c h,i Ψ i δk z κ i, ˆσ h k z = c h,i E ˆDκ i Ψ i δk z κ i 4 for arbitrary coefficients c h,i for 1 i 6. We define Hk x,k y,κ i,ω := E ˆDk x,k y,κ i,ω for better readability, and write Hκ i =: H i for the sake of simplicity. As the model assumes the media to be layered with respect to the z-direction, û is therefore transformed back over the z direction. Equations for ũk x,k y,z,ω are derived by performing the transformation. This yields ũ h = c h,i Ψ i e jκiz, σ h = c h,i Hκ i Ψ i e jκiz, 5 By a variation of constants, the general form of the solution can be expressed as ũ = c h,i + c p,i zψ i e jκiz, σ = c h,i + c p,i zh i Ψ i e jκ iz 6 for z-dependent coefficients c p,i z from a particular solution. Calculation of the stress and displacement vector at layer boundaries In order to derive the equations for the particular solution, a discrete impact ata depth d inside a layer is assumed. Concretely, the impact at a point x 0,y 0,z 0,t 0 with constant direction f is assumed in the form px,y,z,t = δx x 0 δy y 0 δz z 0 δt t 0 f. The Fourier transform is applied on the pz. This yields ˆpk x,k y,k z,ω = expjx 0 k x + y 0 k y + z 0 k z + t 0 ωf. 7 As the model is invariant to translation, the origin of the coordinates x, y, z, ω is set to the point of impact, so that ˆp = f. Note that if the origin is assumed in the middle of the layer, a transformation factor expz p is yielded see Figure 1. A Fourier backward transform gives ˆp as right hand side of Eq. 3. It is to be taken into account that the balance of stresses only holds for the components that act at the interface of the layer boundary. the vector σ is reduced to dimension 3 by dropping the other components; the matrix H r is similarly reduced. A layer with width d is considered, loaded with an impact vector ˆp in depth d p. The layer is split at the depth of the impact d p, and two layers with identical material parameters are derived. For the upper virtual layer let ũ u k x,k y,z = d p,ω and ũ m k x,k y,z = 0,ω be the transformed displacement at the upper and lower boundary, respectively. For the lower virtual layer, let ũ n
4 Figure 1: Treatment of a loaded layer: the layer is split by a virtual layer boundary, at the depth of the impact. and ũ n be the displacement vector at the upper boundary and ũ d the one at the lower boundary. Obviously, ũ m = ũ n holds. For the virtual layer we establish the solution for the upper ũ z and the lower ũ z part: ũ z := a i Ψ i e jκiz, ũ z := b i Ψ i e jκiz, 8 such define ũ u := ũ d p,ũ := ũ 1 0, as well as ũ n := ũ 0,ũ d := ũ d d p, The equations for displacements and stress at the point of impact can be established now. The stress vectors σ r,m and σ r,n at the virtual layer boundary differ by ˆp. The coefficient vector a and b consist of a homogenoeous and particular part, a = a h +a p, b = b h + b p where a h = b h. The particular part b p = 0 as the lower virtual layer is treated as an unloaded layer. The particular part a p has to satisfy ũ p,m = a p,i Ψ i = 0, σ p,m = This gives a p as the solution of a linear equation Ψi Ma p = ˆp 0 with M := H r,i Ψ i a p,i H r,i Ψ i = ˆp. 9 1 i 6, ˆp 0 := 0 ˆp. 10 The coefficient of the displacements and stresses of the unified layers are denoted by c = c h + c p. The particular component of the solution is given by c p = a p = M 1ˆp 0. As the origin is in the middle of the layer, correction factors expjκ i z p for the partial components are derived. For an arbitrary z within the layer boundaries, this yields ũz = c h,i +c p,i e jκ iz p Ψ i e jκ iz d 2, σz = c h,i +c p,i e jκ iz p H r,i Ψ i e jκ iz d 2. 11
5 ICSV13, July 2-6, 2006, Vienna, Austria The factors exp jκ i z d/2 are derived from the backward Fourier transform of the i-th component for the upper layer boundary compare with Eq. 5; as a consequence, exp± 1jκ 2 id arise componentwise for the layer boundaries, and expjκ i z p for the particular solution. The matrices Θ and Ξ, as well as the vector p, are defined as follows, in order to express the displacements and stresses at the layer boundaries: This yields Θ := ũ =: Θc h + Ψ i e 1 2 jκ id Ψ i e 1 2 jκ id 1 i 6, Ξ r := H r,i Ψ i e 1 2 jκ id H r,i Ψ i e 1 2 jκ id 1 i 6 12 p := c p,i e jκ iz p 1 i Θu p 0 = ũ h + ũ p, σ r =: Ξ r c h + Ξr,u p 0 = σ r,h + σ r,p. 14 where Θ u, Ξ r,u consist of the upper three rows of Θ resp Ξ r. Substituting c h = Θ 1 ũ ũ p into the stress equations, and defining K := Ξ r Θ 1, it holds that σ r = Ξ r Θ 1 ũ ũ p + σ r,p = Kũ Kũ p + σ r,p, 15 for a loaded layer. The case for an unloaded layer can easily be deduced from the general case by setting the partial component to zero. BALANCE AT LAYER BOUNDARIES Now, a system of n layers is considered. For the i-th layer, let the matrices c i, Θ i, Ξ i, H i, K i denote the matrices defined in the section before, evaluated for the i-th layer. Let the index X u respectively X d denote the upper respectively lower half of a matrix X. The layers are numbered in the canonical way, such that ũ i is the displacement at the upper boundary of the i-th layer. The stress vectors are sometimes denoted alternatively as σ i,u = σ i. The index r is omitted in σ r, Ξ r etc. For unloaded layers, the balance of stresses yields Boundaries of loaded layers K i,d ũi ũ i+1 = K i+1,u ũi+1 ũ i+2 The layer k is assumed to be loaded. Let the indices h and p denote the homogeneous and the particular part of the solution. Eq.14 yields the homogeneous and the particular component of the stress vector as σ p,k = Ξk,u p k 0, σ h,k = Ξ k c k = Kk,uu K k,ud K k,du K k,dd 16 ũ h,k. 17
6 The stress vectors at the layer boundaries equate in the same way, up to the partial component of the solution. With σ k,u = σ k 1,d it holds that σ k 1,d = σ k,u = σ h,k,u + σ p,k,u, and as the stress vector at the upper boundary takes the form σ k,u = Ξ k,u c k + p k see Eq. 14, the balance of stresses at the upper loaded boundary takes the form Ξ k 1,d c k 1 = Ξ k,u c k + p k 18 and leads to K k 1,d ũk 1 ũ k = K k,u ũk ũ k+1 K k,uu ũ p,k + σ p,k,u. 19 after substitution c i = Θ 1 i ũ i ũ p,i for i = k 1,k. The vector ũ p,k is equivalent to ũ p from Eq. 14, and σ p,k,u is equivalent to σ r,p from Eq.14, for the k-th layer. For the lower boundary, since σ k,d = σ k+1,u, this yields K k,d ũk ũ k+1 + K k,du ũ p = K k+1,u ũk+1 ũ k Top boundary and bottom boundary For the displacement ũ 1, which occurs at the boundary between air and soil, the boundary conditions are assumed to be of von Neumann type. For the bottom halfspace, boundary conditions are derived by assuming that for a monofrequent wave with angular frequency Ω, the wave front κ i z + Ωt = 0 travels downward with time: z z Ω t = Re = Re t κ i > 0 ReΩκ i This eliminates three degrees of freedom. Thus, the matrices Θ n, Ξ n and K n shrink to size 3 3. At the bottom layer, this yields K n 1,d ũn 1 ũ n = K n ũ n. 22 Synthesis Again, a system of n layers is considered, including the k-th one, which is assumed to be loaded. Putting together the results, the equation to be solved is established defining the supervector U := ũ 1,...,ũ n T. So, it delivers a linear equation LU = V, where L is an n n block band matrix with 3 3 blocks.
7 ICSV13, July 2-6, 2006, Vienna, Austria Figure 2: Left: The relative error of U for a four layer system; Right: Difference between a four layer system and the same system split within the first layer. Obtaining the grids With the solution U it is possible to compute the stress vector as well as the displacement vector at arbitrary points. Extending the definition of Θ, Ξ and ˆp yields Θ s z := Ψ s,i e jκ s,iz ds 2, Ξ s z := H s,i Ψ s,i e jκ s,iz ds 2 23, 1 i 6 1 i 6 ˆp s z := { cp,i e jκ id p z 1 i 6 if z d p 0 z > d p. 24 for the layer with index s and for 0 z d s. The formulas for the displacement vector uk x,k y,z,ω and the stress vector σk x,k y,z,ω take the form u s k x,k y,z = Θ s zc h,s + p s z, σ s k x,k y,z = Ξ s zc h,s + p s z. 25 The coefficients c h,s and ĉ s are c h,s = Θ 1 s û h from Eq. 14. By implementing a grid in the z-direction, the four dimensional grid in k x,k y,z,ω can easily be derived by evaluating in z. A Fourier back transform in k x and k y finally yields a grid in the coordinates x,y,z,ω, which was the intended result. TESTING THE MODEL Numerical Accuracy: The numerical stability is checked by running a symbolic computation software MAPLE 9.5 with the data, and comparing the results to the ones of the C++ program. Results are computed on various points of the k x,k y [ 10, 10] [ 10, 10] plane for ω = 10.0 khz, accuracy up to the fifth decimal position is detected. See Figure 2 for the results. Physical consistency: The second test is checking physical consistency. For a system of layers, the solution is computed. By introducing a virtual layer boundary at
8 a random depth, one layer is split in two with identical material parameters. The results must coincide at the layer boundaries. Results are calculated for ω = 10.0 khz, k x,k y [ 10, 10] [ 10, 10]. See Figure 2 for the results. CONCLUSIONS The model is designed for a layer structure of n layers consisting of orthotropic media. A general solution within one layer is derived with a Fourier approach, and back Fourier transformed is performed to handle the boundary conditions. Some additional considerations allow the case of the loaded layer to be handled. The balances at the layer boundaries deliver equations for the displacement vectors at the layer boundaries. Thus, the stress and displacement vectors can be calculated at arbitrary points in the soil. The model presents a more realistic approach, as it incorporates the orthotropic behavior of soils. It showed high accuracy in algebraic as well as in physical tests. The problems the model has to deal with are merely caused by the fact that the soil parameters are in general not determined well enough. Empiric tests as well as future research will focus on that, promising attempts have already been made[5] and can probably be applied. ACKNOWLEDGEMENTS This work is financed by the Austrian Federal Ministry for Transport, Innovation and Technology. Thanks to Wolfgang Kreuzer for technical assistance. References [1] W.M. Ewing, W.S. Jardetzky and F. Press, Elastic Waves in Layered Media, McGraw-Hill, New York, [2] F. Ziegler and P. Borejko, Seismic Waves in Layered Soil: The Generalized Ray Theory, in Structural Dynamics-Recent Advances, edited by G. I. Schueller, Springer, Berlin, 1991, pp [3] I. Tsvankin, Seismic Signatures and Analysis of Reflection data in Anisotropic media, volume 29 of Handbook of Geophysical Exploration: Seismic Exploration, Elsevier science, Amsterdam, [4] H. Waubke, Waves in a Layered Orthotropic Medium with and without Random Properties, in 10th International Congress on Sound and Vibration, Stockholm, [5] H. Waubke, P. Balazs, B. Jilge and W. Kreuzer, Waves in Random Layers with Arbitrary Covariance Functions, in 12th International Congress on Sound and Vibration, Lisbon,2005.
LASER GENERATED THERMOELASTIC WAVES IN AN ANISOTROPIC INFINITE PLATE
LASER GENERATED THERMOELASTIC WAVES IN AN ANISOTROPIC INFINITE PLATE H. M. Al-Qahtani and S. K. Datta University of Colorado Boulder CO 839-7 ABSTRACT. An analysis of the propagation of thermoelastic waves
More informationBasic Equations of Elasticity
A Basic Equations of Elasticity A.1 STRESS The state of stress at any point in a loaded bo is defined completely in terms of the nine components of stress: σ xx,σ yy,σ zz,σ xy,σ yx,σ yz,σ zy,σ zx,andσ
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationFinite and Boundary Element Methods in Acoustics
Finite and Boundary Element Methods in Acoustics W. Kreuzer, Z. Chen, H. Waubke Austrian Academy of Sciences, Acoustics Research Institute ARI meets NuHAG Kreuzer, Chen, Waubke (ARI) FEM-BEM-FMM ARI meets
More informationElastic wavefield separation for VTI media
CWP-598 Elastic wavefield separation for VTI media Jia Yan and Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT The separation of wave modes from isotropic elastic wavefields is typically
More informationAvailable online at ScienceDirect. Procedia Engineering 144 (2016 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 144 (016 ) 170 177 1th International Conference on Vibration Problems, ICOVP 015 Propagation of Torsional Surface Wave in Anisotropic
More informationChapter 2 Governing Equations
Chapter Governing Equations Abstract In this chapter fundamental governing equations for propagation of a harmonic disturbance on the surface of an elastic half-space is presented. The elastic media is
More informationReduction of Finite Element Models of Complex Mechanical Components
Reduction of Finite Element Models of Complex Mechanical Components Håkan Jakobsson Research Assistant hakan.jakobsson@math.umu.se Mats G. Larson Professor Applied Mathematics mats.larson@math.umu.se Department
More informationISSN: X (p); (e)
TORSIONA SURFACE WAVE IN SEF-REINFORCED AYER SANDWICHED BETWEEN TWO VISCO-EASIC HAF-SPACES UNDER THE INITIA STRESS Nidhi Dewangan, Sanjeev A. Sahu and Soniya Chaudhary S.G.G. Govt P.G. College, Kurud,
More informationWe briefly discuss two examples for solving wave propagation type problems with finite differences, the acoustic and the seismic problem.
Excerpt from GEOL557 Numerical Modeling of Earth Systems by Becker and Kaus 2016 1 Wave propagation Figure 1: Finite difference discretization of the 2D acoustic problem. We briefly discuss two examples
More informationwhere d is the vibration direction of the displacement and c is the wave velocity. For a fixed time t,
3 Plane waves 3.1 Plane waves in unbounded solid Consider a plane wave propagating in the direction with the unit vector p. The displacement of the plane wave is assumed to have the form ( u i (x, t) =
More informationA MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC PROBLEMS IN SEMI-INFINITE MEDIA
8 th GRACM International Congress on Computational Mechanics Volos, 2 July 5 July 205 A MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC
More informationOptimal thickness of a cylindrical shell under dynamical loading
Optimal thickness of a cylindrical shell under dynamical loading Paul Ziemann Institute of Mathematics and Computer Science, E.-M.-A. University Greifswald, Germany e-mail paul.ziemann@uni-greifswald.de
More informationMECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso
MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Thermodynamics Derivation Hooke s Law: Anisotropic Elasticity
More informationTHE REFLECTION PHENOMENA OF SV-WAVES IN A GENERALIZED THERMOELASTIC MEDIUM
Internat. J. Math. & Math. Sci. Vol., No. 8 () 59 56 S67 Hindawi Publishing Corp. THE REFLECTION PHENOMENA OF SV-WAVES IN A GENERALIZED THERMOELASTIC MEDIUM ABO-EL-NOUR N. ABD-ALLA and AMIRA A. S. AL-DAWY
More informationAvailable online at ScienceDirect. Procedia Engineering 144 (2016 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 44 (06 ) 46 467 th International Conference on Vibration Problems, ICOVP 05 Propagation of Love waves in composite layered structures
More informationReflection and Transmission coefficient for VTI media
Reflection and Transmission coefficient for VTI media R. K. Sharma and R. J. Ferguson ABSTRACT VTI, R & T coefficients Presently, we obtain the reflection (R and transmission (T coefficients of plane waves
More informationA NOTE ON RELATIONSHIP BETWEEN FIXED-POLE AND MOVING-POLE APPROACHES IN STATIC AND DYNAMIC ANALYSIS OF NON-LINEAR SPATIAL BEAM STRUCTURES
European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS 212) J. Eberhardsteiner et.al. (eds.) Vienna, Austria, September 1-14, 212 A NOTE ON RELATIONSHIP BETWEEN FIXED-POLE
More information(TRAVELLING) 1D WAVES. 1. Transversal & Longitudinal Waves
(TRAVELLING) 1D WAVES 1. Transversal & Longitudinal Waves Objectives After studying this chapter you should be able to: Derive 1D wave equation for transversal and longitudinal Relate propagation speed
More informationProceedings, International Snow Science Workshop, Banff, 2014
SNOW SLAB FAILURE DUE TO BIOT-TYPE ACOUSTIC WAVE PROPAGATION Rolf Sidler 2 Simon Fraser University, 8888 University Drive, Burnaby V5A1S6, Canada ABSTRACT: Even though seismic methods are among the most
More informationSound Propagation through Media. Nachiketa Tiwari Indian Institute of Technology Kanpur
Sound Propagation through Media Nachiketa Tiwari Indian Institute of Technology Kanpur Lecture 14 Waves Propagation in Solids Overview Terminology: phasor, wave front, wave number, phase velocity, and
More informationUNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES
UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More informationEVALUATION OF THE EFFECTS OF TEMPERATURE ON RAILPAD PROPERTIES, RAIL DECAY RATES AND NOISE RADIATION
EVALUATION OF THE EFFECTS OF TEMPERATURE ON RAILPAD PROPERTIES, RAIL DECAY RATES AND NOISE RADIATION R. A. Broadbent, D. J. Thompson and C. J. C. Jones ISVR, University of Southampton, Southampton, SO7
More informationGround vibrations from accelerating and braking road vehicles
Loughborough University Institutional epository Ground vibrations from accelerating and braking road vehicles This item was submitted to Loughborough University's Institutional epository by the/an author.
More informationLinearized theory of elasticity
Linearized theory of elasticity Arie Verhoeven averhoev@win.tue.nl CASA Seminar, May 24, 2006 Seminar: Continuum mechanics 1 Stress and stress principles Bart Nowak March 8 2 Strain and deformation Mark
More informationEffects of tool eccentricity on wave dispersion properties in borehole acoustic logging while drilling
Effects of tool eccentricity on wave dispersion properties in borehole acoustic logging while drilling Yibing Zheng and M. Nafi Toksöz Earth Resources Laboratory Dept. of Earth, Atmospheric, and Planetary
More informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationRocking behaviour of a rigid foundation with an arbitrary embedment
Rocking behaviour of a rigid foundation with an arbitrary embedment H. Moghaddasi Islamic Azad University, Qazvin, Iran. M. Kalantari University of Tehran, Tehran, Iran N. Chouw The University of Auckland,
More informationEFFECT OF COUPLE-STRESS ON THE REFLECTION AND TRANSMISSION OF PLANE WAVES AT AN INTERFACE
International Journal of Modern Engineering Research (IJMER) www.ijmer.com Vol., Issue, pp-5-8 ISSN: 9-665 EFFET OF OUPLE-STRESS ON THE REFLETION AND TRANSMISSION OF PLANE WAVES AT AN INTERFAE Mahabir
More informationPEAT SEISMOLOGY Lecture 3: The elastic wave equation
PEAT8002 - SEISMOLOGY Lecture 3: The elastic wave equation Nick Rawlinson Research School of Earth Sciences Australian National University Equation of motion The equation of motion can be derived by considering
More informationChapter 4 Analysis of a cantilever
Chapter 4 Analysis of a cantilever Before a complex structure is studied performing a seismic analysis, the behaviour of simpler ones should be fully understood. To achieve this knowledge we will start
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationInternational Journal of Advanced Engineering Technology E-ISSN
Research Article INTEGRATED FORCE METHOD FOR FIBER REINFORCED COMPOSITE PLATE BENDING PROBLEMS Doiphode G. S., Patodi S. C.* Address for Correspondence Assistant Professor, Applied Mechanics Department,
More informationDevelopment and application of time-lapse ultrasonic tomography for laboratory characterisation of localised deformation in hard soils / soft rocks
Development and application of time-lapse ultrasonic tomography for laboratory characterisation of localised deformation in hard soils / soft rocks Erika Tudisco Research Group: Stephen A. Hall Philippe
More informationICSV14 Cairns Australia 9-12 July, 2007
ICSV14 Cairns Australia 9-1 July, 7 STUDY ON THE CHARACTERISTICS OF VIBRATION AND ACOUSTIC RADIATION OF DAMAGED STIFFENED PANELS Hong Ming 1, Guo Xin-yi, Liu Lian-hai 3 and Li Gen-tian 4 1 Department of
More informationFree vibration analysis of beams by using a third-order shear deformation theory
Sādhanā Vol. 32, Part 3, June 2007, pp. 167 179. Printed in India Free vibration analysis of beams by using a third-order shear deformation theory MESUT ŞİMŞEK and TURGUT KOCTÜRK Department of Civil Engineering,
More information2 Basic Equations in Generalized Plane Strain
Boundary integral equations for plane orthotropic bodies and exterior regions G. Szeidl and J. Dudra University of Miskolc, Department of Mechanics 3515 Miskolc-Egyetemváros, Hungary Abstract Assuming
More informationAnalytical Mechanics: Elastic Deformation
Analytical Mechanics: Elastic Deformation Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 1 / 60 Agenda Agenda
More informationMulti Degrees of Freedom Systems
Multi Degrees of Freedom Systems MDOF s http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March 9, 07 Outline, a System
More informationEigenvalues and Eigenvectors: An Introduction
Eigenvalues and Eigenvectors: An Introduction The eigenvalue problem is a problem of considerable theoretical interest and wide-ranging application. For example, this problem is crucial in solving systems
More informationMacroscopic theory Rock as 'elastic continuum'
Elasticity and Seismic Waves Macroscopic theory Rock as 'elastic continuum' Elastic body is deformed in response to stress Two types of deformation: Change in volume and shape Equations of motion Wave
More informationDiophantine Equations
Diophantine Equations Michael E. Pohst Institut für Mathematik Technische Universität Berlin May 24, 2013 Mordell s equation y 2 = x 3 + κ is one of the classical diophantine equations. In his famous book
More informationMEC-E8001 FINITE ELEMENT ANALYSIS
MEC-E800 FINIE EEMEN ANAYSIS 07 - WHY FINIE EEMENS AND IS HEORY? Design of machines and structures: Solution to stress or displacement by analytical method is often impossible due to complex geometry,
More informationDIMENSIONAL ANALYSIS OF A HIGH POWER ULTRASONIC SYSTEM USED IN ELECTRODISCHARGE MACHINING
Nonconventional Technologies Review Romania, March, 2012 2012 Romanian Association of Nonconventional Technologies DIMENSIONAL ANALYSIS OF A HIGH POWER ULTRASONIC SYSTEM USED IN ELECTRODISCHARGE MACHINING
More informationStress-induced transverse isotropy in rocks
Stanford Exploration Project, Report 80, May 15, 2001, pages 1 318 Stress-induced transverse isotropy in rocks Lawrence M. Schwartz, 1 William F. Murphy, III, 1 and James G. Berryman 1 ABSTRACT The application
More informationEffect of rail unevenness correlation on the prediction of ground-borne vibration from railways
Effect of rail unevenness correlation on the prediction of ground-borne vibration from railways Evangelos Ntotsios; David Thompson Institute of Sound and Vibration Research, University of Southampton,
More informationEffects of initial stresses on guided waves in unidirectional plates
Arch. Mech., 65,, pp. 3 6, Warszawa 3 Effects of initial stresses on guided waves in unidirectional plates X. M. ZHANG, J. G. YU School of Mechanical and Power Engineering Henan Polytechnic University
More informationEULERIAN DERIVATIONS OF NON-INERTIAL NAVIER-STOKES EQUATIONS
EULERIAN DERIVATIONS OF NON-INERTIAL NAVIER-STOKES EQUATIONS ML Combrinck, LN Dala Flamengro, a div of Armscor SOC Ltd & University of Pretoria, Council of Scientific and Industrial Research & University
More informationParameter Design of High Speed Coupling for 6 MW Wind Turbine Considering Torsional Vibration
Parameter Design of High Speed Coupling for 6 MW Wind Turbine Considering Torsional Vibration JongHun Kang 1, Junwoo Bae 2, Seungkeun Jeong 3, SooKeun Park 4 and Hyoung Woo Lee 1 # 1 Department of Mechatronics
More informationClassifying and building DLMs
Chapter 3 Classifying and building DLMs In this chapter we will consider a special form of DLM for which the transition matrix F k and the observation matrix H k are constant in time, ie, F k = F and H
More informationThe Basic Properties of Surface Waves
The Basic Properties of Surface Waves Lapo Boschi lapo@erdw.ethz.ch April 24, 202 Love and Rayleigh Waves Whenever an elastic medium is bounded by a free surface, coherent waves arise that travel along
More informationBending of Simply Supported Isotropic and Composite Laminate Plates
Bending of Simply Supported Isotropic and Composite Laminate Plates Ernesto Gutierrez-Miravete 1 Isotropic Plates Consider simply a supported rectangular plate of isotropic material (length a, width b,
More informationModelling vibration from surface and underground railways as an evolutionary random process
icccbe 010 Nottingham University Press Proceedings of the International Conference on Computing in Civil and Building Engineering W Tizani (Editor) Modelling vibration from surface and underground railways
More informationThe Finite Element Method
The Finite Element Method 3D Problems Heat Transfer and Elasticity Read: Chapter 14 CONTENTS Finite element models of 3-D Heat Transfer Finite element model of 3-D Elasticity Typical 3-D Finite Elements
More informationMeasurement of deformation. Measurement of elastic force. Constitutive law. Finite element method
Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationChapter 5 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
More informationStudy on elevated light rail induced vibration attenuation along the surrounding ground
Study on elevated light rail induced vibration attenuation along the surrounding ground Changqing Liu ; Yude Zhou ; Ying Tu 3 ; Weimin Xu 4 Shanghai Academy of Environmental Sciences 508 Qinzhou Rd, 0033
More informationIntroduction to structural dynamics
Introduction to structural dynamics p n m n u n p n-1 p 3... m n-1 m 3... u n-1 u 3 k 1 c 1 u 1 u 2 k 2 m p 1 1 c 2 m2 p 2 k n c n m n u n p n m 2 p 2 u 2 m 1 p 1 u 1 Static vs dynamic analysis Static
More informationMATHEMATICAL MODELLING OF HUMAN LIMB IN MATLAB
MATHEMATICAL MODELLING OF HUMAN LIMB IN MATLAB J. Daněk Department of Mathematics, University of West Bohemia Abstract The contribution deals with mathematical modelling of human limb and numerical simulation
More informationNon-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises
Non-linear and time-dependent material models in Mentat & MARC Tutorial with Background and Exercises Eindhoven University of Technology Department of Mechanical Engineering Piet Schreurs July 7, 2009
More informationRobert W. Vestrum and R. James Brown
Group versus phase velocity in laboratory measurements on orthorhombic samples Robert W. Vestrum and R. James Brown ABSTRACT In laboratory measurements of traveltimes across anisotropic materials, there
More informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More information= (G T G) 1 G T d. m L2
The importance of the Vp/Vs ratio in determining the error propagation and the resolution in linear AVA inversion M. Aleardi, A. Mazzotti Earth Sciences Department, University of Pisa, Italy Introduction.
More informationElastic wave-equation migration for laterally varying isotropic and HTI media. Richard A. Bale and Gary F. Margrave
Elastic wave-equation migration for laterally varying isotropic and HTI media Richard A. Bale and Gary F. Margrave a Outline Introduction Theory Elastic wavefield extrapolation Extension to laterally heterogeneous
More informationInfluence of Irregularity and Rigidity on the Propagation of Torsional Wave
Applied Mathematical Sciences, Vol. 4, 00, no. 7, 805-86 Influence of Irregularity and Rigidity on the Propagation of Torsional Wave S. Gupta, A. Chattopadhyay and S. Kundu Department of Applied Mathematics
More informationCodal Provisions IS 1893 (Part 1) 2002
Abstract Codal Provisions IS 1893 (Part 1) 00 Paresh V. Patel Assistant Professor, Civil Engineering Department, Nirma Institute of Technology, Ahmedabad 38481 In this article codal provisions of IS 1893
More informationReceiver. Johana Brokešová Charles University in Prague
Propagation of seismic waves - theoretical background Receiver Johana Brokešová Charles University in Prague Seismic waves = waves in elastic continuum a model of the medium through which the waves propagate
More informationLamb Wave Propagation Modeling Using Cellular Automata
6th European Workshop on Structural Health Monitoring - We.3.E.2 More info about this article: http://www.ndt.net/?id=14171 Lamb Wave Propagation Modeling Using Cellular Automata P. KLUSKA 1, W. J. STASZEWSKI
More informationGaussian vectors and central limit theorem
Gaussian vectors and central limit theorem Samy Tindel Purdue University Probability Theory 2 - MA 539 Samy T. Gaussian vectors & CLT Probability Theory 1 / 86 Outline 1 Real Gaussian random variables
More informationMohr's Circle for 2-D Stress Analysis
Mohr's Circle for 2-D Stress Analysis If you want to know the principal stresses and maximum shear stresses, you can simply make it through 2-D or 3-D Mohr's cirlcles! You can know about the theory of
More informationON THE NUMERICAL ANALYSIS OF COMPOSITE MATERIAL
The 4th International Conference Advanced Composite Materials Engineering COMAT 8- October, Brasov, Romania O THE UMERICAL AALYSIS OF COMPOSITE MATERIAL D.D. icoara Transilvania University, City Brasov,
More informationNDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.
CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationROTATING RING. Volume of small element = Rdθbt if weight density of ring = ρ weight of small element = ρrbtdθ. Figure 1 Rotating ring
ROTATIONAL STRESSES INTRODUCTION High centrifugal forces are developed in machine components rotating at a high angular speed of the order of 100 to 500 revolutions per second (rps). High centrifugal force
More informationFinite Elements for Elastic Shell Models in
Elastic s in Advisor: Matthias Heinkenschloss Computational and Applied Mathematics Rice University 13 April 2007 Outline Elasticity in Differential Geometry of Shell Geometry and Equations The Plate Model
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationResearch Article Dispersion of Love Waves in a Composite Layer Resting on Monoclinic Half-Space
Applied Mathematics Volume 011, Article ID 71349, 9 pages doi:10.1155/011/71349 Research Article Dispersion of Love Waves in a Composite Layer Resting on Monoclinic Half-Space Sukumar Saha BAS Division,
More informationGeneration and Propagation of vibrations induced by high-speed railways
Generation and Propagation of vibrations induced by high-speed railways João Manso Abstract In the years to come, Portugal has several challenges to overcome and one is to try to modernize its train network.
More informationEmbedded Foundation with Different Parameters under Dynamic Excitations
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 2287 Embedded Foundation with Different Parameters under Dynamic Excitations Jaya K P 1 and Meher Prasad
More informationApplication of Algebraic Inverse Method to Surface Wave Testing of Pavements Using Free Plate Solution
Application of Algebraic Inverse Method to urface Wave Testing of avements Using Free late olution T. Akhlaghi Abstract The use of surface waves of the Rayleigh type enables the properties of the component
More informationConstitutive Equations
Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module
More informationFig. 1. Circular fiber and interphase between the fiber and the matrix.
Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In
More information2C9 Design for seismic and climate changes. Jiří Máca
2C9 Design for seismic and climate changes Jiří Máca List of lectures 1. Elements of seismology and seismicity I 2. Elements of seismology and seismicity II 3. Dynamic analysis of single-degree-of-freedom
More informationComputation Time Assessment of a Galerkin Finite Volume Method (GFVM) for Solving Time Solid Mechanics Problems under Dynamic Loads
Proceedings of the International Conference on Civil, Structural and Transportation Engineering Ottawa, Ontario, Canada, May 4 5, 215 Paper o. 31 Computation Time Assessment of a Galerkin Finite Volume
More informationMixing Quantum and Classical Mechanics: A Partially Miscible Solution
Mixing Quantum and Classical Mechanics: A Partially Miscible Solution R. Kapral S. Nielsen A. Sergi D. Mac Kernan G. Ciccotti quantum dynamics in a classical condensed phase environment how to simulate
More informationJournal of Solid Mechanics and Materials Engineering
and Materials Engineering Analysis of In-Plane Problems for an Isotropic Elastic Medium with Many Circular Holes or Rigid Inclusions Mutsumi MIYAGAWA, Jyo SHIMURA, Takuo SUZUKI and Takanobu TAMIYA Dept.of
More informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More informationMP463 QUANTUM MECHANICS
MP463 QUANTUM MECHANICS Introduction Quantum theory of angular momentum Quantum theory of a particle in a central potential - Hydrogen atom - Three-dimensional isotropic harmonic oscillator (a model of
More informationVirtual distortions applied to structural modelling and sensitivity analysis. Damage identification testing example
AMAS Workshop on Smart Materials and Structures SMART 03 (pp.313 324) Jadwisin, September 2-5, 2003 Virtual distortions applied to structural modelling and sensitivity analysis. Damage identification testing
More informationELASTIC WAVE PROPAGATION IN THREE-DIMENSIONAL PERIODIC COMPOSITE MATERIALS
ELASTIC WAVE PROPAGATION IN THREE-DIMENSIONAL PERIODIC COMPOSITE MATERIALS B. Auld, Y. Shui, Y. Wang To cite this version: B. Auld, Y. Shui, Y. Wang. ELASTIC WAVE PROPAGATION IN THREE-DIMENSIONAL PERI-
More informationACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY THE SUPERPOSITION METHOD
Journal of Sound and Vibration (1999) 219(2), 265 277 Article No. jsvi.1998.1874, available online at http://www.idealibrary.com.on ACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY
More information7.2.1 Seismic waves. Waves in a mass- spring system
7..1 Seismic waves Waves in a mass- spring system Acoustic waves in a liquid or gas Seismic waves in a solid Surface waves Wavefronts, rays and geometrical attenuation Amplitude and energy Waves in a mass-
More informationUNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich
UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST
More informationA HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS
A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,
More information16. GAUGE THEORY AND THE CREATION OF PHOTONS
6. GAUGE THEORY AD THE CREATIO OF PHOTOS In the previous chapter the existence of a gauge theory allowed the electromagnetic field to be described in an invariant manner. Although the existence of this
More informationELASTICITY (MDM 10203)
LASTICITY (MDM 10203) Lecture Module 5: 3D Constitutive Relations Dr. Waluyo Adi Siswanto University Tun Hussein Onn Malaysia Generalised Hooke's Law In one dimensional system: = (basic Hooke's law) Considering
More informationA Continuum Theory for the Natural Vibrations of Spherical Viral Capsids
Viral Capsids A Continuum Theory for the Natural Vibrations of Spherical Viral Capsids Francesco Bonaldi Master s Thesis in Mathematical Engineering Advisors Prof. Paolo Podio-Guidugli Università di Roma
More information