z 0 > 0 z = 0 h; (x, 0)

Similar documents
P321(b), Assignement 1

Chapter 6. Differentially Flat Systems

Lecture 38: Equations of Rigid-Body Motion

Lecture 38: Equations of Rigid-Body Motion

Phys 7221 Homework # 8

S p e c i a l M a t r i c e s. a l g o r i t h m. intert y msofdiou blystoc

Solar Sailing near a collinear point

Two dimensional oscillator and central forces

Path Intergal. 1 Introduction. 2 Derivation From Schrödinger Equation. Shoichi Midorikawa

PHY 5246: Theoretical Dynamics, Fall November 16 th, 2015 Assignment # 11, Solutions. p θ = L θ = mr2 θ, p φ = L θ = mr2 sin 2 θ φ.

Written Examination. Antennas and Propagation (AA ) June 22, 2018.

Vibrations: Two Degrees of Freedom Systems - Wilberforce Pendulum and Bode Plots

Lecture Notes Multibody Dynamics B, wb1413

PHY 5246: Theoretical Dynamics, Fall Assignment # 9, Solutions. y CM (θ = 0) = 2 ρ m

The distance of the object from the equilibrium position is m.

MCE 366 System Dynamics, Spring Problem Set 2. Solutions to Set 2

Question 1: A particle starts at rest and moves along a cycloid whose equation is. 2ay y a

Statistical Mechanics Solution Set #1 Instructor: Rigoberto Hernandez MoSE 2100L, , (Dated: September 4, 2014)

University of Illinois at Chicago Department of Physics

Mathematical Modeling and response analysis of mechanical systems are the subjects of this chapter.

Chemistry 432 Problem Set 1 Spring 2018 Solutions

Distance travelled time taken and if the particle is a distance s(t) along the x-axis, then its instantaneous speed is:

Water is sloshing back and forth between two infinite vertical walls separated by a distance L: h(x,t) Water L

Homework 3. 1 Goldstein Part (a) Theoretical Dynamics September 24, The Hamiltonian is given by

Class Notes for Advanced Dynamics (MEAM535)

Nonlinear Dynamic Systems Homework 1

ECE 3209 Electromagnetic Fields Final Exam Example. University of Virginia Solutions

Pendulum with a vibrating base

the EL equation for the x coordinate is easily seen to be (exercise)

Hertz potentials in curvilinear coordinates

Variation Principle in Mechanics

Magnetostatics III Magnetic Vector Potential (Griffiths Chapter 5: Section 4)

MATH 189, MATHEMATICAL METHODS IN CLASSICAL AND QUANTUM MECHANICS. HOMEWORK 2. DUE WEDNESDAY OCTOBER 8TH,

Asynchronous Training in Wireless Sensor Networks

Goldstein Problem 2.17 (3 rd ed. # 2.18)

Forces, Energy and Equations of Motion

Dynamics of structures

M2A2 Problem Sheet 3 - Hamiltonian Mechanics

Professor Fearing EE C128 / ME C134 Problem Set 2 Solution Fall 2010 Jansen Sheng and Wenjie Chen, UC Berkeley

ME Machine Design I

CHAPTER 4 GENERAL MOTION OF A PARTICLE IN THREE DIMENSIONS

Hamiltonian. March 30, 2013

m d v dt = q E + q c v B

PHYSICS 110A : CLASSICAL MECHANICS MIDTERM EXAM #2

Aperture Antennas 1 Introduction

Transient Phenomena in Quantum Bound States Subjected to a Sudden Perturbation

Semi-Classical Theory of Radiative Transitions

We first write the integrand into partial fractions and then integrate. By EXAMPLE 27 we have the identity

7 Pendulum. Part II: More complicated situations

The Motion of A Test Particle in the Gravitational Field of A Collapsing Shell

Electrodynamics HW Problems 06 EM Waves

Introduction to Nonlinear Control Lecture # 4 Passivity

Truncation Errors Numerical Integration Multiple Support Excitation

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 5 Solutions

Classical Mechanics Comprehensive Exam Solution

Semester University of Sheffield. School of Mathematics and Statistics

MATH325 - QUANTUM MECHANICS - SOLUTION SHEET 11

EE Homework 3 Due Date: 03 / 30 / Spring 2015

Solutions for homework 5

Consider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be.

University of California, Berkeley Department of Mechanical Engineering ME 104, Fall Midterm Exam 1 Solutions

The CR Formulation: BE Plane Beam

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012

Classical mechanics of particles and fields

Errata Instructor s Solutions Manual Introduction to Electrodynamics, 3rd ed Author: David Griffiths Date: September 1, 2004

LMI Methods in Optimal and Robust Control

Physics 322 Midterm 2

Waves and the Schroedinger Equation

Examensarbete. Separation of variables for ordinary differential equations. LiTH - MAT - EX / SE

Lecture 6: The Orbits of Stars

Mechanics Physics 151


Examination paper for Solution: TMA4285 Time series models

Chapter 5 Oscillatory Motion

PHY 5246: Theoretical Dynamics, Fall September 28 th, 2015 Midterm Exam # 1

Chemistry 532 Practice Final Exam Fall 2012 Solutions

System Control Engineering 0

AA242B: MECHANICAL VIBRATIONS

I. HARTLE CHAPTER 8, PROBLEM 2 (8 POINTS) where here an overdot represents d/dλ, must satisfy the geodesic equation (see 3 on problem set 4)

Control Systems. Frequency domain analysis. L. Lanari

Part III. The Synchro-Betatron Hamiltonian

CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems

University of California, Berkeley Physics H7B Spring 1999 (Strovink) SOLUTION TO PROBLEM SET 12 Solutions by P. Pebler. ( 5 gauss) ẑ 1+[k(x + ct)] 2

Generalized projective synchronization between two chaotic gyros with nonlinear damping

Department of Physics, Korea University Page 1 of 8

Yell if you have any questions

NYU Physics Preliminary Examination in Electricity & Magnetism Fall 2011

Relationship between POTENTIAL ENERGY and FORCE

Mechanical System. Seoul National Univ. School of Mechanical and Aerospace Engineering. Spring 2008

Chapter 5 Fast Multipole Methods

MAS 315 Waves 1 of 7 Answers to Example Sheet 3. NB Questions 1 and 2 are relevant to resonance - see S2 Q7

221A Lecture Notes Electromagnetic Couplings

Physics 214 Midterm Exam Solutions Winter 2017

Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems

Tutorial Exercises: Central Forces

Theoretical Basis of Modal Analysis

CHAPTER 12 OSCILLATORY MOTION

Module I: Electromagnetic waves

Numerical Modeling of Collective Effects in Free Electron Laser

Assignment 2. Goldstein 2.3 Prove that the shortest distance between two points in space is a straight line.

Transcription:

n = (q 1,..., q n ) T (,, t) V (,, t) L(,, t) = T V d dt ( ) L q i L q i = 0, i = 1,..., n. l l 0 l l 0 l > l 0 {x} + = max(0, x) x = k{l l 0 } +ˆ, k > 0 ˆ

(x, z) x z (0, z 0 ) (0, z 0 ) z 0 > 0 x z = 0 h; l; θ; ϕ; k; m; l 0. h m z 0 (x, 0) z 0 = l cos θ + h cos ϕ x = l sin θ + h cos ϕ. (x, ϕ) θ l θ l = l(x, ϕ) = x 2 + z 2 0 + h 2 2h(x sin ϕ + z 0 cos ϕ) = (x h sin ϕ) 2 + (z 0 h cos ϕ) 2. V (x, ϕ) = k(l l 0 ) 2, T = 1 2 mẋ2.

L = T V x L = 1 2 mẋ2 k ( (x h sin ϕ)2 + (z 0 h cos ϕ) 2 l 0 ) 2. ( (x ) 2k(x h sin ϕ) h sin ϕ)2 + (z 0 h cos ϕ) 2 l 0 mẍ =. (x h sin ϕ)2 + (z 0 h cos ϕ) 2 ϕ t E 1 ( ) 2 2 mẋ2 + k (x h sin ϕ)2 + (z 0 h cos ϕ) 2 l 0 = E. ẋ 0 ϕ ϕ ϕ + x x ϕ ϕ = ϕ(t) x mẍ + V x = 0, V t V ϕ x t t mẋ = K + O( t) K ẋ O( t) 0 t y 0 t [y] [x] = O( t).

t 0 x [x] = 0 ϕ x ϕ = ϕ ± ϕ = ϕ x N j m = N k{ j l 0 } + j, 1 j = j j. j j l 0 + a l 0 a > 0 j = (l 0 + a)(cos 2πj 2πj, sin ), j = 1,..., N. N N X = (0, 0) (0, 0) l 0 + a = ϵ ϵ > 0 l = l 0 + a = ϵl(u, v) j l 0 a ϵl(u cos 2πj N + v sin 2πj N ) + O(ϵ2 ).

c j = cos 2πj N s j = sin 2πj N j = (c j, s j ) ϵ (a(u, v) (c j u + s j v)(c j, s j )) + O(ϵ 2 ). (ü, v) = Nk ( 1 + a ) (u, v). 2m l 0 + a θ j = r = (r, 0) = (l 0 + a)(cos θ, sin θ) = (0, 0) l 0 = ((l 0 + a) cos θ r) 2 + (l 0 + a) 2 sin 2 θ l 0 p = l 0 + a l 0 = p 2 + r 2 2rp cos θ l 0. OXZ ϕ x k( p 2 + r 2 2rp cos θ l 0 )( cos ϕ), x 2 = 2 + 2 2 cos ϕ 2r cos ϕ = 2 + 2 2 = 2r 2 2rp cos θ. cos ϕ = r p cos θ x ( ) r p cos θ k r p cos θ l 0. p2 + r 2 2rp cos θ θ 2 l 2 0

ψ = 0 ψ 0 ψ π p 2 + r 2 2pr cos ψ = l 2 0. r ( ) π r p cos θ = ˆ 2k r p cos θ l 0 dθ ψ p2 + r 2 2rp cos θ ˆ ( ) π r p cos θ = ˆ 2k (π ψ)r + p sin ψ l 0 p2 + r 2 2rp cos θ dθ. r ψ = 0 ψ

ϕ = ϕ = 0 ϕ = ϕ + > 0 x ϕ ϕ ϕ ẋ < 0 ϕ = ϕ 0 ϕ ϕ + ϕ = ϕ ± ϕ = ϕ x 0 ϕ ± < π 2 x = h sin ϕ ± > 0 (x, ẋ) (h sin ϕ ±, 0) h z 0 f a = 1 2k(z 0 l 0 ). 2π mz 0

. x 4 2 0-2 -4-0.08-0.04 0 0.04 0.08 0.12 x ϕ = ϕ = 0 ẋ < 0 ϕ = ϕ + = 30 ẋ > 0 (x, ẋ) = (0.001, 0) ϕ = 0 1 mm SI m = 0.2 10 3 kg, z 0 = 0.2 m, h = 0.025 m. l 0 = 1z 2 0 = 0.1 m E 2 GP a 2 10 9 Nm 2 3 µm 3 10 6 m k = EA/l 0

A k = 10EA = 20 10 9 9π 10 12 = 18π 10 2. f a = 1 k 1 2π m = 15 π 8.5. 8 9 Hz ϕ + 30 ϕ + 1 mm 10 cm 8.03 Hz ϕ x = h sin ϕ x (x h sin ϕ) (x h sin ϕ) 2 ẋ < 0 ϕ = ϕ = 0 (x, ẋ) = (x, 0) x > 0 ( x, 0) ẋ > 0 ( x, 0) (2h sin ϕ + + x, 0) (x, 0) (x + 2h sin ϕ +, 0) A A = 2h sin ϕ +. ϕ ẋ > 0 ϕ ẋ < 0 ẋ = 0 x = 0 x = h sin ϕ + ẋ = 0

ω > 0

web b l z y y x body m (a) (b) ( c) x ψ ζ ϕ = ψ + ζ z = 0 P = (b cos ψ, b sin ψ, 0), ψ = ωt b > 0 z = 0 2π y = 0 x > 0 ω 1 ω 2π ω 2π C θ 0 θ π ψ + ζ ζ x C (x, y) ϕ ψ + ζ δ = ζ ϕ b 0 C = (x, y, z) l x = b cos ψ + l sin θ cos(ψ + ζ) y = b sin ψ + l sin θ sin(ψ + ζ) z = l cos θ.

L(ζ, θ, t) = 1 2 m(ẋ2 + ẏ 2 + ż 2 ) + mgz (x, y, z) ψ = ωt l θ l(ω + ζ) 2 cos θ sin θ bω 2 cos θ cos ζ + g sin θ = 0, l ζ sin θ + 2l(ω + ζ) θ cos θ + bω 2 sin ζ = 0. b = 0 l 2 θ sin 2 θ = H θ θ z ϕ = ψ + ζ 1 ω 2π 1 g 2π l z θ ζ sin ζ = 0 ζ = 0 ζ = π ζ = 0 θ θ + (0, π) 2 B ζ = π θ ( π, π) 2 b < l [ ( ) ] 2 3 4 g b 3 ω > 1. l l g 9.8 ms 2 l 0.05 m g 14 l

y x z z y x z x y (x, y, z) 0 < t < 20 ω = 5 (θ, θ, ζ, ζ) = (0.04, 0.01, 0.01, 0.05) ω = 9 (θ, θ, ζ, ζ) = (0.04, 0.01, 0.01, 0.05) ω = 12 (θ, θ, ζ, ζ) = (0.02, 0.07, 0.05, 0.01) ω < 13 g = 9.8, l = 0.05, b = 1 2 l. ω 5 13 2 Hz ω = 5 ω = 9 ω z > 0 z = 0

( 3 ) ( 1 = (l 0 + a), 1, 2 2 1 = (l 0 + a) ) 3, 1, 2 2 3 = (l 0 + a)(0, 1). x = 0 m = 3 1 k{ j l 0 } + ( j ) j. m = 0.2 10 3, l 0 = 0.17, a = 0.03, k = 18π 10 2 0.05 m a x

0.06 (a) y 0 (b) y 0.01 0-0.04 x -0.02 0 0.02-0.01-0.0004 0 0.0004 (x, y) 0 < t < 3 (x, ẋ, y, ẏ) = (0.02, 0, 0.05, 0) (x, ẋ, y, ẏ) = (0, 0.03, 0.02, 0) r A = 2π k (1 l 0 ) > 0 m 2p r r θ 2 = Ar + O(r 2 ), r 2 θ = H H r 2 r = Ar + H2 r = ( ) 1 3 r 2 Ar2 + H2. 2r 2 V (r) = 1 2 Ar2 + H2 2r 2 r > 0 r 2 = H/ A 1 2π V (r) = 2π 1 2A x

H = 0 r = 0

ω ω

\\ \

r I = l 0 π ψ r p cos θ p2 + r 2 2rp cos θ dθ = l π 0 2r ψ 0 2r 2 2rp cos θ p2 + r 2 2rp cos θ dθ p 2 ( I = l π 0 p 2 r 2 ) π 2r ψ p2 + r 2 2rp cos θ dθ + p2 + r 2 2rp cos θ dθ. ψ F (ϕ, k) E(ϕ, k) ϕ dθ ϕ F (ϕ, k) = 1 k2 sin 2 θ, E(ϕ, k) = 1 k 2 sin 2 θ dθ. ϕ = π 2 F ( π 2, k) = π 2 (1+ 1 4 k2 + 9 64 k4 +O(k 6 )), E( π 2, k) = π 2 (1 1 4 k2 3 64 k4 +O(k 6 )). A > B > 0 0 x π 1 A B cos θ dθ = 2 A + B F (δ, R) 0 A B cos θ dθ = 2 A + BE(δ, R) 2B sin θ A B cos θ.

δ = sin 1 (A + B)(1 cos θ) 2B, R = 2(A B cos θ) A + B. [ I = l 0 r (p r)f (δ θ, s) (p + r)e(δ θ, s) + 4pr sin θ p2 + r 2 2pr cos θ θ δ δ s = 2 pr. p + r r ψ = 0 ( = ˆ 2k πr + l ) 0 r ((p r)f ( π, s) (p + r)e( π, s)). 2 2 r ( ) = ˆ 2k π l 0 r + O(r 2 ). 2p ] π ψ,