Introduction to Waves in Structures. Mike Brennan UNESP, Ilha Solteira São Paulo Brazil

Size: px

Start display at page:

Download "Introduction to Waves in Structures. Mike Brennan UNESP, Ilha Solteira São Paulo Brazil"

Peter Hunt
5 years ago
Views:

1 Introduction to Waves in Structures Mike Brennan UNESP, Ilha Solteira São Paulo Brazil

2 Waves in Structures Characteristics of wave motion Structural waves String Rod Beam Phase speed, group velocity Low and high frequency behaviour Wave transmission and impedance

3 Wave Motion

4 Wave Motion Wave motion Wave motion Particle motion Particle motion

5 Harmonic Wave Motion ( x, ) = A sin( ωt ± k ) or W ( x, t ) = A cos( ωt ± kx ) W t x

6 Harmonic Wave Motion ( x, ) = A sin( ωt ± k ) or W ( x, t ) = A cos( ωt ± kx ) W t x Fix position Fix time

7 Harmonic Wave Motion ( x, ) = A sin( ωt ± k ) W t x W W { } ( ) ( ), Im j ω t ± x t = A e kx ( t ) x, = Ae ω e ± j t jkx ( x ) = e jkx W A ( ) = e jkx W x A

8 The Wavenumber (propagating waves) time distance T ω = π T Temporal frequency k λ π = = λ ω c Spatial frequency (wavenumber) k is the phase change per unit distance Ae jkx

9 Superposition of Waves

10 Interference of Waves ( ) Wave motion A sin ωt kx Wave motion A sin ( ωt kx φ ) Wave motion W ( x, t ) ( x, t ) = sin( ωt x ) + sin( ω k φ ) W A k A t x φ W = A x ( x, t ) cos sin( ωt k )

11 Standing Waves Wave motion A sin( ωt kx ) Wave motion A sin( ω t + kx ) Standing wave W ( x, t ) ( x, t ) = sin( ωt x ) + A sin( ωt + kx ) W A k ( x, t ) sin( ωt ) cos( x ) W = A k

12 Beating Wave motion Wave motion A A sin sin ( ω t k x ) 1 1 ( ω t k x ) Wave motion W ( x, t ) ( x, t ) = sin( ω t k x ) + A sin( ω k ) W A t x 1 1 W ω ω k k ω + ω k + k A x x ( x, t ) = cos t sin t

13 Waves in solids

14 Waves motion in a string

15 Wave motion in a string y + y dx x dy θ dx y T String segment T x x + dx ( sinθ ) ( sinθ ) df = T T y x + dx x

16 Wave motion in a string Applying the Taylor s series expansion gives f f x dx f x dx x ( + ) = ( ) x ( T sin θ ) df = ( T sinθ ) + dx T x ( sinθ ) y x x If θ is small, then sin θ θ = y x, then y df = T dx y x Stiffness force

17 Wave motion in a string The inertia force in the y direction is given by df y = ρ dx L t where ρ L is the mass per unit length Applying Newton`s second law gives the wave equation y x y 1 y = c t 0 T where c = is the phase speed ρ L

18 Wave motion in a string c = T ρ L is a constant - the wave is non-dispersive Wave speed c

19 Free harmonic wave motion in a string Start with the wave equation y 1 y = 0 x c t Assume time harmonic motion, gives d y k y 0 dx This has a solution + = where k y() x = Ae + Be jkx jkx ω = c

20 Wave motion in a string - impedance y f y T x y = ( ) j t kx Ae ω boundary condition ( ) f = T sinθ Ty 0 y y y = = jkae xx y yɺ = = j ωae t ( ω ) j t kx ( ω ) j t kx Z f y = = yɺ ( 0) ρ c L which is real and hence damping-like

21 Wave motion in a rod

22 Wave motion in a rod u f x x + dx f f + x df bar element of area S and density ρ f S = E u x..(1) Applying Newton`s second law gives f f + dx f = ρsdx x t u ()

23 Wave motion in a rod Now from (1) f = ES x x which combines with () to give the wave equation u u 1 u = x c t 0 E where c = is the phase speed (constant with frequency) ρ

24 Wave motion in a rod - impedance f x u rod area S and density ρ u = ( ) j t kx Ae ω f = ESu ( 0) ( ) u = jkae ω uɺ = j ωae ω j t kx ( ) j t kx f Z = = S E uɺ ( 0) ρ which is real and hence damping-like

25 In-Plane Wave motion rod c = E ρ plate c = ρ E ( 1 ν ) ν = Poisson's ratio 3D solid c = E ( 1 ν ) ( 1+ )( 1 ) ρ ν ν

26 Shear Waves

27 Shear Waves The wave equation is given by where w 1 w = x c t 0 c = G ρ Shear modulus ( ) E = 1+ν G ν = Poisson's ratio E = Young's modulus G = shear modulus Note that torsional waves are also shear waves

28 Wave motion in a beam

29 Wave motion in a beam y w Beam segment of length dx crosssectional area S M M + dm V V + dv Summing forces V V V + dx = ρsdx x x t w (1)

30 V x Wave motion in a beam Summing moments about the any point on the right hand face From (1) and () w + ρs = t M M + dx M Vdx = x 0 V 0 () Bending stiffness M w = and M = EI x x Combining these equations gives the wave equation EI w w 4 + ρs = 0 4 x t

31 Wave speed in a beam Note that the wave equation is fourth order rather than second order. What is the phase speed? Is it dependent on frequency? First consider the equation for a single degree of freedom mass-spring-damper system d u ω u 0 n dt + = (Ordinary differential equation - time) Now consider the wave equation for a rod u 1 u = x c t 0 (Partial differential equation - space and time)

32 Wave speed in a beam Assume harmonic motion the wave equation becomes where k d u k u 0 dx + = (Ordinary differential equation - space) ω = c Now consider the wave equation for a beam 4 w ρs w + = x EI t 4 0 (4 th order Partial differential equation - space and time)

33 Wave speed in a beam Assume harmonic motion the wave equation becomes 4 d w k w = dx (Ordinary differential equation - space) where k ω ρs = = c EI 1 4 ω 1 So c EI = ρs 1 4 ω 1 Note that the wave speed is dependent on frequency - it is dispersive

34 Group velocity Wave motion Wave motion A A sin sin ( ω t k x ) 1 1 ( ω t k x ) Carrier Velocity c g W ( x, t ) W ω ω k k ω + ω k + k A x x ( x, t ) = cos t sin t c g ω ω ω k k k 1 = = 1 In the limit c g d ω = dk

35 Group velocity flexural waves The wavenumber is given by k ρs = EI ω 1 EI ω = k ρs d ω EI = k c dk = ρs The group velocity is twice the phase velocity 1

36 Wave characteristics

37 Wave characteristics

38 Wave characteristics

39 Bending Waves motion in a beam Start with the wave equation EI w w 4 + ρs = 0 4 x t Assume time harmonic motion, gives 4 d w k w = dx This has a solution where w() x = A e + A e + A e + A e Nearfield waves kx kx jkx jkx ω k = c

40 Wave motion in a beam - impedance w x f w = A e e + A e e j ωt kx j ωt jkx 4 Beam of area S and density ρ 0 = EIw ( 0) Z f EI A = = + wɺ ( 0) 1/4 3/4 ( )( ρ ) 1/ ω (1 j ) f = EIw ( 0) which has a real part - damping-like and a positive imaginary part - mass-like

41 Bending Waves in a plate

42 Bending Waves in a plate Bending wave behaviour similar to the beam but in two dimensions The equation of motion is w w + h = t t 4 D ρ 0 where = + + x x y y Assume time harmonic motion, gives 4 4 w k w = 0 and D = where 11 k Eh ( ν ) = 3 ω c

43 Bending Wavenumber in a plate The bending number for a plate is given by k 1ρ1 1 ( ν ) 4 1 = ω Eh The bending number for a rectangular beam is given by k 1ρ = Eh 1 4 ω 1

44 Bending Waves frequency limits Bending wave behaviour stems from the Euler- Bernoulli beam equation Neglects rotary inertia of the beam Assumes infinite shear stiffness Assumptions are ok for low frequencies, when λ > 10 depth of the beam

45 Material Properties Material Young s modulus E (N m - ) Density ρ (kg m -3 ) Poisson s ratio, ν Loss factor, η Mild Steel e11 7.8e Aluminium Alloy 7.1e10.7e Brass 10e10 8.5e < Copper 1.5e10 8.9e Glass 6e10.4e Cork 1..4e Rubber 1e6-1e9 1e Plywood 5.4e9 6e 0.01 Perspex 5.6e9 1.e Light Concrete 3.8e9 1.3e Brick 1.6e10 e

46 Wavelength for 1 mm Euler-Bernoulli beam 10 3 Wavelength (mm) Frequency (Hz)

47 Timoshenko beam y w Beam segment of length dx crosssectional area S M ρi M + dm V V + dv Equation of motion x E x κg x t κg t I w E w I w I 1 S w + ρ + + ρ + ρ = 0 4 t 4

48 Timoshenko beam Assuming time harmonic motion Equation of motion the equation of motion becomes 4 S d w 4 S L 4 k k L B d w dx which leads to k k k + w = κ x κ S L B 0 - shear - Longitudinal - Bending k k k κ 4 κ 1 S 4 1 S = + k ± k + k L B L As ω 0, then k = k As ω, then k = k and k = k B S L

49 Wave motion in beams summary At low frequencies Torsional waves non-dispersive waves that propagate at the shear wave speed Longitudinal waves non-dispersive waves that propagate at the compressional wave speed Bending waves nearfield waves (do not propagate) and flexural waves that propagate at the bending wave speed At high frequencies Torsional waves, Longitudinal waves Flexural waves that propagate at (a) compressional wave speed (b) shear wave speed)

50 Reflection and Transmission of Waves Hard boundary Note the phase differences Soft boundary

51 Reflection and Transmission of Waves 1 High speed to low speed A Z Z r = A Z Z A t = + 1 A Z Z Low speed to high speed

52 Surface wave in a Solid Rayleigh wave Water wave

53 Summary Characteristics of wave motion Structural waves String Rod Beam Phase speed, group velocity Low and high frequency behaviour Wave transmission and impedance

COPYRIGHTED MATERIAL. Index

COPYRIGHTED MATERIAL. Index Index A Admissible function, 163 Amplification factor, 36 Amplitude, 1, 22 Amplitude-modulated carrier, 630 Amplitude ratio, 36 Antinodes, 612 Approximate analytical methods, 647 Assumed modes method,