Four Bar Linkages in Two Dimensions. A link has fixed length and is joined to other links and also possibly to a fixed point.

Similar documents
So far: simple (planar) geometries

SCHOOL OF COMPUTING, ENGINEERING AND MATHEMATICS SEMESTER 2 EXAMINATIONS 2011/2012 DYNAMICS ME247 DR. N.D.D. MICHÉ

Week 9 Chapter 10 Section 1-5

Chapter 3. r r. Position, Velocity, and Acceleration Revisited

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Spring 2002 Lecture #13

ENGN 40 Dynamics and Vibrations Homework # 7 Due: Friday, April 15

CHAPTER 10 ROTATIONAL MOTION

Section 8.3 Polar Form of Complex Numbers

Physics 4B. Question and 3 tie (clockwise), then 2 and 5 tie (zero), then 4 and 6 tie (counterclockwise) B i. ( T / s) = 1.74 V.

1 Matrix representations of canonical matrices

PHYS 1443 Section 003 Lecture #17

(c) (cos θ + i sin θ) 5 = cos 5 θ + 5 cos 4 θ (i sin θ) + 10 cos 3 θ(i sin θ) cos 2 θ(i sin θ) 3 + 5cos θ (i sin θ) 4 + (i sin θ) 5 (A1)

10/9/2003 PHY Lecture 11 1

Important Instructions to the Examiners:

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

Physics 111: Mechanics Lecture 11

Angular Momentum and Fixed Axis Rotation. 8.01t Nov 10, 2004

Section 8.1 Exercises

ˆ (0.10 m) E ( N m /C ) 36 ˆj ( j C m)

EN40: Dynamics and Vibrations. Homework 7: Rigid Body Kinematics

Chapter 11 Angular Momentum

Part C Dynamics and Statics of Rigid Body. Chapter 5 Rotation of a Rigid Body About a Fixed Axis

Physics 2A Chapter 3 HW Solutions

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project

Difference Equations

MAE140 - Linear Circuits - Fall 13 Midterm, October 31

Gravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11)

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

12. The Hamilton-Jacobi Equation Michael Fowler

Week 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product

= 1.23 m/s 2 [W] Required: t. Solution:!t = = 17 m/s [W]! m/s [W] (two extra digits carried) = 2.1 m/s [W]

General Tips on How to Do Well in Physics Exams. 1. Establish a good habit in keeping track of your steps. For example, when you use the equation

Chapter 10 Rotational motion

Rotational Dynamics. Physics 1425 Lecture 19. Michael Fowler, UVa

PHYSICS 203-NYA-05 MECHANICS

Please initial the statement below to show that you have read it

Indeterminate pin-jointed frames (trusses)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

International Mathematical Olympiad. Preliminary Selection Contest 2012 Hong Kong. Outline of Solutions

EN40: Dynamics and Vibrations. Homework 4: Work, Energy and Linear Momentum Due Friday March 1 st

Example: (13320, 22140) =? Solution #1: The divisors of are 1, 2, 3, 4, 5, 6, 9, 10, 12, 15, 18, 20, 27, 30, 36, 41,

MAE140 - Linear Circuits - Winter 16 Midterm, February 5

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.

The classical spin-rotation coupling

Mathematical Preparations

Work is the change in energy of a system (neglecting heat transfer). To examine what could

Mechanics Cycle 3 Chapter 9++ Chapter 9++

SCALARS AND VECTORS All physical quantities in engineering mechanics are measured using either scalars or vectors.

I 2 V V. = 0 write 1 loop equation for each loop with a voltage not in the current set of equations. or I using Ohm s Law V 1 5.

Linear Feature Engineering 11

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

Lecture 12: Discrete Laplacian

Slide. King Saud University College of Science Physics & Astronomy Dept. PHYS 103 (GENERAL PHYSICS) CHAPTER 5: MOTION IN 1-D (PART 2) LECTURE NO.

Physics 114 Exam 3 Spring Name:

Modeling of Dynamic Systems

Causal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim

Chapter 3 Differentiation and Integration

MAE140 - Linear Circuits - Winter 16 Final, March 16, 2016

Homework 2: Kinematics and Dynamics of Particles Due Friday Feb 7, 2014 Max Score 45 Points + 8 Extra Credit

Solutions to Selected Exercises

PHY2049 Exam 2 solutions Fall 2016 Solution:

11. Dynamics in Rotating Frames of Reference

Physics 207: Lecture 27. Announcements

Change. Flamenco Chuck Keyser. Updates 11/26/2017, 11/28/2017, 11/29/2017, 12/05/2017. Most Recent Update 12/22/2017

Math 261 Exercise sheet 2

Chapter 8. Potential Energy and Conservation of Energy

Physics 231. Topic 8: Rotational Motion. Alex Brown October MSU Physics 231 Fall

The Geometry of Logit and Probit

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.

Chapter 9. The Dot Product (Scalar Product) The Dot Product use (Scalar Product) The Dot Product (Scalar Product) The Cross Product.

From Biot-Savart Law to Divergence of B (1)

Physics 40 HW #4 Chapter 4 Key NEATNESS COUNTS! Solve but do not turn in the following problems from Chapter 4 Knight

Review: Fit a line to N data points

Physics 2A Chapters 6 - Work & Energy Fall 2017

Rigid body simulation

MTH 263 Practice Test #1 Spring 1999

ISQS 6348 Final Open notes, no books. Points out of 100 in parentheses. Y 1 ε 2

Math1110 (Spring 2009) Prelim 3 - Solutions

Selected Student Solutions for Chapter 2

Quadratic speedup for unstructured search - Grover s Al-

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

Summary with Examples for Root finding Methods -Bisection -Newton Raphson -Secant

Solutions for Homework #9

Lecture 10 Support Vector Machines II

measurement and the charge the electric field will be calculated using E =. The direction of the field will be the . But, so

Unit 5: Quadratic Equations & Functions

Structure and Drive Paul A. Jensen Copyright July 20, 2003

AERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY

Important Dates: Post Test: Dec during recitations. If you have taken the post test, don t come to recitation!

Physics 207 Lecture 6

Spin-rotation coupling of the angularly accelerated rigid body

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

(2mn, m 2 n 2, m 2 + n 2 )

Anglo-Chinese Junior College H2 Mathematics JC 2 PRELIM PAPER 1 Solutions

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Transcription:

Four bar lnkages 1 Four Bar Lnkages n Two Dmensons lnk has fed length and s oned to other lnks and also possbly to a fed pont. The relatve velocty of end B wth regard to s gven by V B = ω r y v B B = +y = ω B k( + y) = ω B ( y + ) 1 Vector r = + y n two dmensons θ ω B Note ω B s the counter clockwse drecton angular velocty. The relatve velocty s always perpendcular to the lne onng the lnkage ends, Formally ω = k ω and ω r = 0 0 ω y 0 k but Equaton 1 can be seen by nspecton. Or note V B = ω B r. The component s V B snθ = V B y r = ω B y. Ths we have n Equaton 1. The component arses n a smlar way. The relatve acceleraton has two components. One component s centrpetal and the other azmuthal or, one mght say, crcumferental. There s no Corols acceleraton snce the "radus" or barlnk length s constant. The centrpetal acceleraton of B w. r. t. a Brad = ω (ω r) = ω B 2 ( + y) 2 The azmuthal acceleraton a Bazm = ρ = B ( y + ) 3 a B = a Brad + a Bazm a Bazm = B 2 k(+y) = B 2 (-y+) y a Brad = - ω B 2 (+y) B = +y B, angular acceleraton ω B, angular velocty = ω B 2 ( + y) + B ( y + ) = ( ω B 2 B y ) + ( ω B 2 y + B )

Four bar lnkages 2 Eample 260 mm B 240 mm 100 mm 100 mm y ω ο O 1 C 60 mm O s rotatng at constant angular velocty ω o counterclockwse. Fnd all veloctes and acceleratons. Veloctes Take O. In eq. 1 60, y 80 The velocty of w. r. t. O s v O = ω o ( 80 +( 60) ) = ω o (80 + 60) 4 Smlarly v B = ω B ( 100 + 240 ) 5 where ω B s the ( at present unknown) counterclockwse angular velocty of B nd v CB = ω CB [ ( 180) + 0 )] = ω CB 180 6 where ω CB s the ( at present unknown) counterclockwse angular velocty of CB. Thus we can epress the velocty of C whch s zero as 0 = v CO = v O + v B + v CB = ω o (80 + 60) + ω B ( -100 + 240 ) + ω CB 180 7 Compare components 0 = ω o 60 + ω B 240 ω B = ω o 4 Compare components; 0 = ω o 80 ω B 100 + ω CB 180 ω CB = 7 ω o 12 8 9 Fnally v B = ω CB ( 180) = 105 ω o n -drecton as n the fgure 10

Four bar lnkages 3 cceleratons 260 mm B 240 mm 100 mm 100 mm y ω ο O C 1 60 mm Usng the formulae: Frst a O = ω o 2 ( 60 + 80) 11 Then usng eqq. 2 & 3 a B = ω B 2 ( 240+100) + B ( 100+ 240) 12 and a BC = ω B 2 ( + y) + B ( y + ) Note that we have ω B = ω o 4 ; B s yet to be determned. nd a CB = ω CB 2 ( 180) + CB [ 0 (- 180)] = 180 CB + 180 ω CB 2 13 Thus the acceleraton of C, whch s zero s gven by 0 = a CO = a O + a B + a CB = ω o 2 ( 60 80) ω B 2 ( 240+100) + B ( 100 + 240) + 180 ω CB 2 + 180 CB Compare coeffcents: 0 = 80 ω o 2 100ω B 2 + 240 B + 180 ω CB 2 14 whence 240 B = 80 ω o 2 + 100ω B 2 80 ω CB 2 = ω o 2 80 + 100 49 16 180 144 and B = 5 48 ω o 2 15 Compare coeffcents: 0 = 60ω o 2 240ω B 2 100 B + 180 CB 0 = ω o 2 60 240 49 100 69 16 +180 144 48 16 + 180 CB

Four bar lnkages 4 Thus CB = 83 432 ω o 2 16 From the angular acceleratons the lnear acceleratons follow farly easly. To recap the data: ω B = ω o 4 ; ω CB = 7 ω o 12 B = 5 48 ω o 2 ; CB = 83 432 ω o 2 From eq. 11 a O = ( 60 80)ω o 2 and we may wrte a BC = a BO = 180 CB 180 ω CB 2 = 415 12 245 4 = 34 58 61 25 260 mm B 100 mm β ω ο O y 60 mm 1 C

Four bar lnkages 5 lternatve approach We may work drectly wth the dagrams for veloctes and acceleratons eamnng the geometry of each. Velocty Dagram -125 g v B (v) v CB O -80-60 -40-20 () -10 v B -20 () v β -30-40 -50 e () d -60 The angles and β are taken from the geometry of the system. sn = 5 12 ; cos = 13 13 ;tan = 5 12 snβ = 4 5 ; cosβ = 3 5 () We may enter v at the angle β as shown and wth ω o = 1 usng mm s 1 as unts. Or take v wth component 80 vertcally and 60 horzontally. () v B s at an angle as shown and snce v B has only a horzontal component the length c g = 60 () Now tan = 5 12 then c-d = 5 12 60 = 25 (v) Whence v CB = 105 We have v B = 60 cos = 60 13 12 = 65 mm s 1 Snce the length of B s 260 mm ω B = 65 260 = 0 25 rad s 1 Smlarly the angular velocty of BC ω BC = 105 180 = 7 12 rad s 1.

Four bar lnkages 6 cceleraton Dagram s before the angles and β are taken from the geometry of the system. sn = 5 12 ; cos = 13 13 (v) 415/12 f g 10 20 30 40 50 60 snβ = 4 5 ; cosβ = 3 5 The order of steps s shown by roman numerals on the dagram. () We may put n a straght off and enter the 80 dmenson g a ( settng ω o = 1) () & ()a e s the radal acceleraton of B towards and ths we get knowng ω B. a e = 260 (ω o /4) 2. Hence the dmensons a b and b c. (v) The length f-e s the radal acceleraton of B w. r. t. C We have ω CB = 7/12 and so f-e = 245/4 (v) So, knowng the dmensons on the rght we fnd d e = 25. (v)snce c-e s the transverse acceleraton to B w. r. t. and we have the angle and can compute the dmenson d-c. -20-40 -60-80 25 (v) β e a B d 125/12 (v) 245/4 (v) c a b 15 () 80 () 25/4 () (v) Whence f O follows. The dagram s solved. To get the angular acceleratons, BC = f- O BC = 415 12 180 = 83 432 rad s 2 clockwse : see dagram. B = c e 260 = 325 12 260 = 5 48 rad s 2 counter clockwse : see dagram. Ths process could be done appromately wth a ruler and pencl. It does not always happen that the numbers come n as convenent ntegers as they do here.

Four bar lnkages 7 Slder and Crank & Connectng Rod ω 0 c r o θ B The slder B moves only horzontally. It s as f B were on the end on an nfntely long rod whose fed end was at ±. Veloctes. v = ω o ( c snθ + c cosθ) = ω o c ( snθ + cosθ ) 1 v B = ω B ( r ( sn) + r cos ) = ω B r ( sn+ cos) 2 Thus v B = v + v B = ( ω ο c snθ + ω ΒΑ r sn) + ( ω o c cosθ + ω B r sn) 3 Snce the vertcal or component of v B = 0 we get ω B = ω o c cosθ r cos 4 and v B = ω o c( snθ cosθ tan) 5 Note tan = c snθ c cosθ so v B = ω o c snθ + c snθ cosθ c snθ cosθ c cosθ = ω o c snθ c cosθ 6 We could have arrved at the last epresson by usng the cosne rule and the geometry of the connectons: r 2 = c 2 + 2 2 c cosθ 7 Take the dervatve w. r. t. tme: or snce ω o = θ ' 0 = 0 + 2 ' 2 c ' cosθ + 2 c snθ θ ' ' = ω o c snθ c cosθ = v B 8 The velocty dagram s shown below.

Four bar lnkages 8 v v B = r ω B c ω 0 θ v B Usng the sne rule we get: v B sn(θ+) = cω o cos ; v b = c ω osn(θ + ) cos o 9 whch agrees wth equaton 5 above. Note that there are varous other ways of approachng ths analyss: look carefully at the problem. For eample we could use the sne rule on the orgnal geometry sn c = snθ r : sn = r c snθ Take the dervatve w. r. t. : cos c ' = cosθ r θ ' ; ' = c cosθ r cos θ' so ω B = ' = c cosθ r cos θ' = ω o c cosθ r cos : whch s eq. 4

Four bar lnkages 9 cceleraton ω 0 c r o θ B We suppose that the crank O s gong at constant speed so O = 0 a = ω o 2 ( c cosθ + c snθ) = ω o 2 c ( cosθ + snθ) nd a B = ω B 2 ( r cos r sn) + B ( ( r sn) + r cos) = ω B 2 r ( cos + sn) + B r ( sn + cos) and a B = a + a B = ( ω o 2 c cosθ ω B 2 r cos+ B r sn) + ( ω o 2 c snθ+ ω B 2 r sn + B r cos) Snce B moves only horzontally the component of acceleraton s zero and so B = 0 = ω o 2 c snθ+ ω B 2 r sn + B r cos 1 r cos ( ω o 2 c snθ ω B 2 r sn) Snce ω B = ω o c cosθ r cos B = ω o 2 c r 2 cos 3 ( r snθ cos2 c sn cos 2 θ) and ω a o 2 c B = r r (snθ tan cosθ) c cos2 θ cos 3 Ths can hardly be counted as a smple epresson. For partcular cases the solutons may be smpler. In any event the complety of the formula suggests ether to use a computer or to nsert numercal values early n the analyss.

Four bar lnkages 10