NONLINEAR OSCILLATIONS OF A FLOATING ELASTIC PLATE

Similar documents
Problem 1. Problem Engineering Dynamics Problem Set 9--Solution. Find the equation of motion for the system shown with respect to:

DETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.

Damped Vibration of a Non-prismatic Beam with a Rotational Spring

DYNAMIC ANALYSIS OF BEAM-LIKE STRUCTURES SUBJECT TO MOVING LOADS

METHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall

Fluid Physics 8.292J/12.330J % (1)

Principle Of Superposition

COMM 602: Digital Signal Processing

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS

Application of Homotopy Perturbation Method for the Large Angle period of Nonlinear Oscillator

2C09 Design for seismic and climate changes

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

Chapter 9: Numerical Differentiation

Numerical Study on MHD Flow And Heat Transfer With The Effect Of Microrotational Parameter In The Porous Medium

Streamfunction-Vorticity Formulation

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations

Stopping oscillations of a simple harmonic oscillator using an impulse force

DISPERSION OF THE SEISMIC RAYLEIGH WAVES IN ELASTIC CRUSTAL LAYER ON VISCOELASTIC MANTLE

Free Surface Hydrodynamics

Boundary layer problem on conveyor belt. Gabriella Bognár University of Miskolc 3515 Miskolc-Egyetemváros, Hungary

17 Phonons and conduction electrons in solids (Hiroshi Matsuoka)

ANALYSIS OF DAMPING EFFECT ON BEAM VIBRATION

( ) (( ) ) ANSWERS TO EXERCISES IN APPENDIX B. Section B.1 VECTORS AND SETS. Exercise B.1-1: Convex sets. are convex, , hence. and. (a) Let.

Time-Domain Representations of LTI Systems

FREE VIBRATIONS OF SIMPLY SUPPORTED BEAMS USING FOURIER SERIES

MULTIPLE TIME SCALES SOLUTION OF AN EQUATION WITH QUADRATIC AND CUBIC NONLINEARITIES HAVING FRAC- TIONAL-ORDER DERIVATIVE

Analysis of composites with multiple rigid-line reinforcements by the BEM

MATH 31B: MIDTERM 2 REVIEW

PHY4905: Nearly-Free Electron Model (NFE)

FREE VIBRATION RESPONSE OF A SYSTEM WITH COULOMB DAMPING

Added mass estimation of open flat membranes vibrating in still air

Vibratory Motion. Prof. Zheng-yi Feng NCHU SWC. National CHung Hsing University, Department of Soil and Water Conservation

Chapter 2 Feedback Control Theory Continued

Numerical Method for Blasius Equation on an infinite Interval

Where do eigenvalues/eigenvectors/eigenfunctions come from, and why are they important anyway?

THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE ABSTRACT

Measurement uncertainty of the sound absorption

Lecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods

NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE

A STUDY ON MHD BOUNDARY LAYER FLOW OVER A NONLINEAR STRETCHING SHEET USING IMPLICIT FINITE DIFFERENCE METHOD

On forward improvement iteration for stopping problems

mx bx kx F t. dt IR I LI V t, Q LQ RQ V t,

FIR Filter Design: Part II

Physics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.

CHAPTER 8 SYSTEMS OF PARTICLES

Dynamics of Structures 5th Edition Chopra SOLUTIONS MANUAL

The Transition Dipole Moment

Exact Solutions of the Generalized Benjamin Equation and (3 + 1)- Dimensional Gkp Equation by the Extended Tanh Method

L 5 & 6: RelHydro/Basel. f(x)= ( ) f( ) ( ) ( ) ( ) n! 1! 2! 3! If the TE of f(x)= sin(x) around x 0 is: sin(x) = x - 3! 5!

TMA4205 Numerical Linear Algebra. The Poisson problem in R 2 : diagonalization methods

1 1 2 = show that: over variables x and y. [2 marks] Write down necessary conditions involving first and second-order partial derivatives for ( x0, y

10-701/ Machine Learning Mid-term Exam Solution

Types of Waves Transverse Shear. Waves. The Wave Equation

Harmonic balance approaches to the nonlinear oscillators in which the restoring force is inversely proportional to the dependent variable

FIXED-FREE AGAINST FREE-FREE BEAMS FOR DYNAMIC YOUNG S MODULUS OF WOOD By: Mehran Roohnia

Dynamical stability and parametrical vibrations of the laminated plates with complex shape

LECTURE 14. Non-linear transverse motion. Non-linear transverse motion

Chapter 4. Fourier Series

The Transition Dipole Moment

(c) Write, but do not evaluate, an integral expression for the volume of the solid generated when R is

Avoidance of numerical singularities in free vibration analysis of Euler-Bernoulli beams using Green functions

Engineering Mechanics Dynamics & Vibrations. Engineering Mechanics Dynamics & Vibrations Plane Motion of a Rigid Body: Equations of Motion

6.1 Analysis of frequency selective surfaces

Solution of Differential Equation from the Transform Technique

NUMERICAL METHOD FOR SINGULARLY PERTURBED DELAY PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS

: Transforms and Partial Differential Equations

Modified Decomposition Method by Adomian and. Rach for Solving Nonlinear Volterra Integro- Differential Equations

Sequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018

Inelastic spherical caps with defects

2.004 Dynamics and Control II Spring 2008

Lesson 03 Heat Equation with Different BCs

Office: JILA A709; Phone ;

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS

Chapter 7 z-transform

MATHEMATICAL MODELLING OF ARCH FORMATION IN GRANULAR MATERIALS

The axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.

Math 257: Finite difference methods

Vibrational Spectroscopy 1

Chapter 7: The z-transform. Chih-Wei Liu

SOLUTIONS TO EXAM 3. Solution: Note that this defines two convergent geometric series with respective radii r 1 = 2/5 < 1 and r 2 = 1/5 < 1.

Taylor expansion: Show that the TE of f(x)= sin(x) around. sin(x) = x - + 3! 5! L 7 & 8: MHD/ZAH

REGRESSION (Physics 1210 Notes, Partial Modified Appendix A)

PLEASE MARK YOUR ANSWERS WITH AN X, not a circle! 1. (a) (b) (c) (d) (e) 3. (a) (b) (c) (d) (e) 5. (a) (b) (c) (d) (e) 7. (a) (b) (c) (d) (e)

(4 pts.) (4 pts.) (4 pts.) b) y(x,t) = 1/(ax 2 +b) This function has no time dependence, so cannot be a wave.

MATHEMATICS. 61. The differential equation representing the family of curves where c is a positive parameter, is of

Research Article Frequency Equations for the In-Plane Vibration of Circular Annular Disks

The target reliability and design working life

A comparative study of a system of Lotka-Voltera type of PDEs through perturbation methods

Dynamic Instability of Taut Mooring Lines Subjected to Bi-frequency Parametric Excitation

SRC Technical Note June 17, Tight Thresholds for The Pure Literal Rule. Michael Mitzenmacher. d i g i t a l

Similarity between quantum mechanics and thermodynamics: Entropy, temperature, and Carnot cycle

Stability Analysis of the Euler Discretization for SIR Epidemic Model

Dynamic Response of Second Order Mechanical Systems with Viscous Dissipation forces

AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC MAPS

The Riemann Zeta Function

x a x a Lecture 2 Series (See Chapter 1 in Boas)

Physics 7440, Solutions to Problem Set # 8

Taylor polynomial solution of difference equation with constant coefficients via time scales calculus

ON THE EIGENFUNCTION EXPANSION METHOD FOR THE CALCULATION OF GREEN S FUNCTIONS

Transcription:

Iteratioal Applied Mechaics Vol 46 o March (Russia Origial Vol 46 o October ) OLIEAR OSCILLATIOS OF A FLOATIG ELASTIC PLATE ÀÅ Buatov ad ÀÀ Buatov The multi-scales method is used to derive third-order equatios of gravitatioal bedig oscillatios of a thi elastic plate floatig o the surface of a homogeeous perfect icompressible fluid of fiite depth The equatios icorporate the compressive force ad oliear acceleratio of vertical displacemets of the plate Based o these equatios the deflectio of the plate ad the velocity potetial of the fluid iduced by a travelig periodic wave of fiite amplitude are expaded ito asymptotic series to terms of the third order of smalless The depedece of the oscillatio characteristics o the elastic modulus ad thicess of the plate compressive force the iitial legth ad steepess of the wave is aalyzed Keywords: thi elastic plate fluid of fiite depth gravitatioal bedig oscillatios iitial wave legth wave steepess Itroductio The oscillatios of a floatig elastic plate were studied i a liear formulatio i [ 7 8 4] disregardig the compressive force ad i [ 6 4] taig the compressive force ito accout The oliear oscillatios of a floatig absolutely flexible plate were studied i [9] The fiite-amplitude bedig oscillatios of a floatig elastic plate were aalyzed i [ 4] disregardig the fact that the acceleratio of vertical displacemet is oliear We will study the oscillatios of a logitudially compressed elastic plate floatig o the surface of a homogeeous perfect icompressible fluid of costat depth i which a periodic wave of fiite amplitude propagates We will use the method of multiple-scale asymptotic expasios [5] ad tae ito accout the fact that the acceleratio of the vertical displacemets of the plate is oliear Problem Formulatio Goverig Equatios Cosider a thi elastic plate floatig o the surface of a homogeeous perfect icompressible fluid of costat depth The plate ad fluid are ubouded i the horizotal directios Let the oscillatios of the plate be oliear ad the motio of the fluid be potetial With dimesioless variables x x z z t g t * ( / g ) * where is the wave umber the problem is to solve the Laplace equatio xx zz ( x H z ) () for the velocity potetial ( x z t) with the followig boudary coditios: D 4 4 Q 4 p x x z x t p t x z () at the plate fluid iterface ( z ) ad Marie Hydrophysical Istitute atioal Academy of Scieces of Uraie Kapitasaya St Sevastopol Uraie 99 e-mail: ewislad@listru Traslated from Priladaya Mehaia Vol 46 o pp 6 7 October Origial article submitted Jue 5 9 6-795//46-9 Spriger Sciece+Busiess Media Ic 9

z () o the bottom ( z H) The iitial coditios (t =)are f( x) t (4) D I () (4) the followig otatio is used: D g D Eh Q Q ( ) g h Å h ad are the ormal elastic modulus thicess desity ad Poisso s ratio of the plate; Q is the logitudial compressive force per uit width of the plate; ( xtis ) the deflectio of the plate or the rise of the plate fluid iterface; is the desity of the fluid; g is the acceleratio of gravity The velocity potetial ad the deflectio of the plate are related by the iematic coditio t x x z (5) The term with the factor i the dyamic coditio () is the iertia of the vertical displacemets of the plate The first term i bracets of this expressio characterizes the oliearity of the vertical acceleratio of the plate which was eglected i [ 4] Equatios for oliear Approximatios To solve problem () (5) we will use the multiple-scales method [5] Let us itroduce two ew variables T t T t slowly varyig from t Ò where is a small yet fiite ad let the followig expasios hold: f f O( ) O( ) f f f f O( ) () Substitutig from () ito () ad () we obtai the followig equalities up to terms of the third order of smalless: z z z x z Let us cosider dyamic () iematic (5) ad iitial (4) coditios ad represet the velocity potetial of the plate fluid iterface z i the form ( xt ) ( xt) z( xt) zz( xt) () Substitutig f f ( x t ) ad z ( xt ) ito () ad (5) respectively expressig the partial derivative with respect to time by the chai rule as t T T T cosiderig the depedece of o x ad t i () ad collectig the coefficiets of lie powers of ad equatig them to zero we obtai the equatios x z ( x H z ) () 4

4 4 D Q x4 x F * zt T (z ) (4) T z L ( ) (5) z ( z H) (6) f ( x) T G t ( ) (7) to fid oliear approximatios I (4) (7): F F * F F F F L G ( ) T z T T z x z T z L T z T z x x z T F T z T T x x z z T z T z z x x z T z T z T z z T T x x 4 4 T z T z z T T z T 5 4 T z T z z 5 T T z zt L 6 z x x x x T T z x 6 x z x z F z xz x F 7 xz x x z x xz x z 4

7 G G xz x xz T T T ote that the terms F ad F o the right-had sides of the dyamic coditios (4) for the secod ( ) ad third ( ) approximatios are due to the oliearity of the acceleratio of the vertical displacemets of the plate Deflectio of the Plate ad Velocity Potetial of the Fluid Equatios () (7) represet the geeral case of usteady oscillatios of fiite amplitude Let us solve these equatios for travelig periodic waves specifyig the first approximatio ( ) of the deflectio of the plate (the rise of the plate fluid iterface) i the form of a wave cos x T ( T T ) () travelig i the egative directio of the x-axis The the iematic coditio (5) yields z si z () To satisfy the boudary coditio (6) o the bottom we represet i the form Substitutig () ito () we obtai b b +cosh( zh)si () +sih H Hece the followig equalities hold: b si b + sih H + cosh( zh) (4) The dispersio equatio below follows from the dyamic coditio (4) with () (4) for oscillatios i liear approximatio: 4 ( Q D )( tahh) tah H (5) The expressio for ( Ò Ò ) i () follows from the subsequet approximatios Substitutig ad ito the right-had side of Eqs (4) ad (5) for the secod approximatio solvig the problem for ad assumig that there is o pricipal harmoic we obtai a cos b si (6) a 4 tah H b + cosh ( z H) ( tahh cothh ) tahh 4 + cosh HtahH ( 5tahH cothh ) ( 4Q 6D ) 4 ( 4Q 6D ) tahh ( tah H) 4 It appears that the fuctio does ot deped o Ò because ( Ò ) The solutios for the first () (4) (5) ad the secod (6) approximatios determie the right-had sides of dyamic (4) ad iematic (5) statemets of the problem for the third approximatio ( ) Elimiatig the terms resposible for secularity we obtai formulas for ad : a cos b si a b + cosh ( z H) + cosh H where ( l l ) tah H l l l ( 9Q 8D ) 4 4

5 l l l l a( cothh 6cothH) cothhcoth H 8 l a 5 coth H cothh ( cothh 5cothH) cothh 8 5 l a cothhcothh cothhcothh cothh 8 ( 9Q 8D ) tahh ( tah H) 4 T l l4( H) coth l a( cothh cothh) cothhcoth H 8 5 l4 l4 l4 l4 a cothh cothh ( cothh cothh) cothh 8 5 l4 a cothhcothh cothhcothh cothh 4 Thus we have the followig third-order expressios for the deflectio ad velocity potetial : a cos b si xt ( ) a Passig to the dimesioal quatities ( * g / x x z / t t/ g a where a is the amplitude of the iitial harmoic) we get acos a a cos a a cos a g / b si a g b si a g b si x ( ) t a g (7) (hereafter the subscript beeath x z t ad the asteris o are omitted) 4 Aalysis of the Results Solutio (7) is valid beyod small eighborhoods of the resoat values of the wave umber that are the positive real roots ad of the equatios ad respectively I the deep-water approximatio ( Í ) these equatios become D ( 76) Q 4 D ( ) Q 4 4 Rejectig the terms represetig the iertia of the plate ( ) we arrive at If there is o compressio ( Q ) the Q Q 4D 4D / Q 9Q 56D 78D / 4

m m 4 4 4 h m 4 h m a b Fig H m H m 4 7 4 75 5 6 7 75 8 85 65 6 6 9 Q m 6 9 Q m 4 5 55 4 45 5 55 6 a Fig b / 4 4D / 4 9D I the shallow-water approximatio ( H ) the resoat values ad are the roots of the equatios solvig which yields 4HD 5D H Q 4 9HD D H Q (8) 4 5 8 H H H Q 6 ( ) 5D / 5 9H 9 H( H Q ) 5D / eglectig the iertia of the plate ( i (8)) leads to Q /( 5D ) Q /( D ) As the umerical result demostrate if the depth of the fluid remais costat the roots ad of the equatios ad icrease with decreasig elastic modulus of the plate ad icreasig compressive force A icrease i the depth of the fluid also leads to a icrease i ad However the upper limit of the depth at which this effect is sigificat decreases with decreasig elastic modulus but icreases with decreasig compressive force If Å ad Q remai costat a decrease i the thicess of the plate leads to a icrease i ad Figure illustrates the depedece of (dashed lies) ad (solid lies) o the thicess of a plate with / 87 4 (ice) for H = m Curves ad correspod to 9 ad 5 8 M/m forq i Fig a ad toq adq D for Å 9 /m i Fig b Figures a ad b shows the isolies of ad o the plae H Q for h = 5 m ad Q D m If the plate is absolutely flexible ((Å Q = ) the solutio (7) for a fluid of fiite depth is valid for ay values of the wave umber [9] 44

/ a / a x m 4 6 8 x m Fig Fig 4 The distributio (calculated by formulas (7)) of the vertical displacemet of the plate fluid iterface (plate deflectio) i the directio of wave propagatio shows that the structure of perturbatios depeds ot oly o the flexural rigidity of the plate ad the compressive force but also o the fluid depth ad the wave legth ( / ) ad steepess ( a ) of the iitial pricipal harmoic The effect of these parameters ca be judged from the chage i the cotributio of the higher harmoics to oscillatios ad hece the deflectio wave profile This is illustrated by Figs ad 4 plotted for Í m t 7 sec 5 m ad Í m t sec m respectively for the same ad / as i Fig The solid dashed ad dash-ad-dot lies correspod to Q = D 95 D for Å 9 /m h m à m The curves of ( x) demostrate the effect of the compressive force o the phase velocity of the flexural wave of the plate It decreases with icreasig compressive force Q for ay ad fixed elastic modulus E The phase lag (the wave propagates i the egative directio of the x-axis) is due also to the icrease i the elastic modulus at costatq ad If adq is fixed the phase velocity icreases with E Whe there is o compressio (Q = ) a icrease i E leads to a icrease i the phase velocity of the flexural wave for all The oliearity of the acceleratio of vertical displacemets causes the phase of oscillatios to shift i the directio of wave motio ad their amplitude to wealy icrease if the plate is absolutely flexible (E = ) If the plate is elastic the effect of the oliearity of the acceleratio o the phase is strog whe It is maifested as a phase lag o matter whether the compressive force is preset or abset This effect weaes with icreasig wavelegth of the iitial pricipal harmoic Coclusios Usig the multiple-scales method we have derived the equatios for three oliear approximatios of the solutio for the oscillatios of a floatig elastic plate with the compressive force ad the oliear acceleratio of vertical displacemets tae ito accout For a travelig periodic wave of fiite amplitude we have obtaied asymptotic expasios of the third order describig the bedig of the plate ad the velocity potetial of the fluid These expasios are valid beyod small eighborhoods of two resoat values of the wave umber The depedece of the bedig characteristics ad resoat wave umbers of a ice plate o its thicess ad elastic modulus the compressive force ad the fluid depth has bee aalyzed as a example It has bee show that these parameters have a effect o the phase velocity of the deflectio wave ad the cotributio of higher harmoics to oscillatios ad hece ad o the deflectio wave profile REFERECES A E Buatov Effect of logitudial compressio o trasiet vibratios of a elastic plate floatig o the surface of a liquid It Appl Mech 7 o 75 79 (98) A E Buatov ad L V Cheresov Trasiet vibratios of a elastic plate floatig o a liquid surface It Appl Mech 6 o 8 878 88 (97) O M Gladu ad V S Fedoseo oliear steady-state oscillatios of a elastic plate floatig o the surface of a fluid of fiite depth Izv A SSSR Meh Zhid Gaza o 46 54 (989) 4 R V Gol dshtei ad A V Marcheo Log waves i the ice fluid system subject to ice pressure i: Electrical ad Mechaical Properties of Ice [i Russia] Gidrometeoizdat Leigrad (989) pp 88 5 5 A H ayfeh Perturbatio Methods Wiley ew Yor (97) 45

6 V I Pozhuev ad P Polyaova Effect of a movig load o a ice cover: ostatioary problem Stroit Meh Rasch Sooruzh o 6 46 5 (99) 7 V A Tacheo ad V V Yaovlev Usteady bedig-gravitatioal waves i fluid-plate system It Appl Mech o 66 7 (984) 8 D E Kheisi Oscillatrios of a ifiite plate floatig o the surface of a perfect fluid: ostatioary problem Izv A SSSR Meh Mashiostr o 6 67 (96) 9 A E Buatov ad A A Buatov Propagatio of surface wave of fiite amplitude i a basi with floatig broe ice It J Offsh Polar Eg 9 o 6 66 (999) A E Buatov ad V V Zharov Formatio of the ice cover s flexural oscillatios by actio of surface ad iteral ship waves Part Surface waves It J Offsh Polar Eg o 7 o (997) D G Daffy The respose of floatig ice to a movig vibratig load Cold Reg Sci Tech 5 64 (99) A D Kerr The critical velocities of a movig o a floatig ice plate that is subjected to i-plae forces Cold Reg Sci Tech o 6 67 76 (98) R M S M Schules R J Hosig ad A D Seyd Waves due to a steadily movig source o a floatig ice plate Part J Fluid Mech 8 97 8 (987) 4 V A Squire R J Hosig A D Kerr ad P J Laghore Movig Loads o Ice Plates Kluwer Dordrecht (996) 46