Mechanisms Course Notebook
Mechanisms Course Notebook
Course Notebook
Instructor:
Jeremy S. Daily, Ph.D., P.E.
Fall 2012
Contents
1. Syllabus 6
1.1. Course Bulletin Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.2. Objectives . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.3. Course Outline . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.4. Course Policies . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.4.1. Text Books . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.4.2. Grading Procedures . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.4.3. Exam Policy . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.4.4. Computer Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.4.5. Late Submission and Absences . . . . . . . . . . . . . . . . . . . . . . . 8
1.4.6. Class Conduct . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.4.7. Academic Misconduct . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.4.8. Center for Student Academic Support . . . . . . . . . . . . . . . . . . . 9
2. Introduction To Mechanisms 10
2.1. Mechanisms Vocabulary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.2. Common Mechanisms . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.3. Kinematic Pairs (a.k.a. Joints) . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
2.3.1. Low Order Pairs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
2.3.2. High Order Pairs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.4. Degrees of Freedom . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.4.1. Mobility . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.4.2. Kutzbach Criteria . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.5. Grashof’s Law for Four-bar Mechanisms . . . . . . . . . . . . . . . . . . . . . . 16
2.5.1. Special Cases . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.5.2. Inversions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.6. Homework Problem Set 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
3. Position Analysis 20
3.1. Loop Closure Equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
3.1.1. Derivation of the Law of Cosines . . . . . . . . . . . . . . . . . . . . . 21
3.1.2. Derivation of the Law of Sines . . . . . . . . . . . . . . . . . . . . . . . 22
3.1.3. Inverted Slider Crank . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
3.1.4. Offset Slider Crank . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
Contents 3
4. Mechanism Synthesis 59
4.1. Geometric Constraint Programming . . . . . . . . . . . . . . . . . . . . . . . . 59
4.2. Homework Problem Set 5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
5. Velocity Analysis 75
5.1. Vector Operations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
5.1.1. Dot Product . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
5.1.2. Cross Product . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
5.1.3. Derivatives of Vector Products . . . . . . . . . . . . . . . . . . . . . . . 75
5.2. Velocity with a Rotating Reference Frame . . . . . . . . . . . . . . . . . . . . . 76
5.3. Graphical Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
5.3.1. Inverted Slider Crank . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
5.3.2. Four-Bar Mechanism . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
5.4. Analytical Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
5.4.1. Inverted Slider Crank . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
5.4.2. Four Bar Mechanism . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
5.5. Homework Problem Set 6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
6. Acceleration Analysis 93
6.1. Accelerations in Four-bar Mechanisms . . . . . . . . . . . . . . . . . . . . . . . 93
6.2. Inverted Slider Crank . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94
6.3. Homework Problem Set 7 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
Contents 4
7. Cams 99
7.1. Types of Cam Followers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
7.1.1. Flat Faced Radial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
7.1.2. Offset Roller Follower . . . . . . . . . . . . . . . . . . . . . . . . . . . 101
7.1.3. Barrel Cam with Roller Follower . . . . . . . . . . . . . . . . . . . . . . 101
7.1.4. Heavy Truck Brake Cams (S-Cams) . . . . . . . . . . . . . . . . . . . . 102
7.2. Cam Follower Motion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 102
7.2.1. Displacement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 102
7.2.2. Velocity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 102
7.2.3. Acceleration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 102
7.2.4. Jerk . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 102
7.3. Cam Follower Profiles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103
7.3.1. Constant Acceleration (Parabolic) . . . . . . . . . . . . . . . . . . . . . 103
7.3.2. Harmonic Motion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 106
7.3.3. Cycloidal Motion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 109
7.4. Cam Design . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
7.5. Homework Problem Set 8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 123
8. Gears 124
8.1. Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 124
8.1.1. Cog and Lantern Gears . . . . . . . . . . . . . . . . . . . . . . . . . . . 124
8.1.2. Common Types of Gears . . . . . . . . . . . . . . . . . . . . . . . . . . 125
8.2. Fundamental Law of Gearing . . . . . . . . . . . . . . . . . . . . . . . . . . . . 126
8.3. Conjugate Profiles and Involutometry . . . . . . . . . . . . . . . . . . . . . . . 126
8.3.1. The Involute Curve . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 127
8.3.2. Cycloidal Profiles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 127
8.3.3. Gear Sizing and Terminology . . . . . . . . . . . . . . . . . . . . . . . 127
8.4. Homework Problem Set 9 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 129
8.5. Gear Train Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 130
8.5.1. Taxonomy of gear trains: . . . . . . . . . . . . . . . . . . . . . . . . . . 130
8.5.2. Gear Train Ratio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 131
8.5.3. Idler Gears . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 131
8.5.4. Compound Gear Trains . . . . . . . . . . . . . . . . . . . . . . . . . . . 131
8.6. Homework Problem Set 10 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 132
8.7. Reverted Gear Train Design . . . . . . . . . . . . . . . . . . . . . . . . . . . . 133
8.7.1. Design Considerations . . . . . . . . . . . . . . . . . . . . . . . . . . . 133
8.8. Planetary Gear Trains . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 134
8.8.1. Advantages . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 134
8.8.2. Vector Approach to Planetary Gear Train Analysis . . . . . . . . . . . . 134
8.8.3. Tabular Planetary Gear Train Analysis . . . . . . . . . . . . . . . . . . . 136
8.8.4. Compound Planetary Gear Trains . . . . . . . . . . . . . . . . . . . . . 138
8.8.5. Plotting Planetary Gears . . . . . . . . . . . . . . . . . . . . . . . . . . 140
Contents 5
1.2. Objectives
This section of the syllabus is designed to give the student a larger picture of the purpose of
this class. This course is designed to provide the students the analytical skill necessary for
1.3. Course Outline 7
analyzing mechanism motion analysis. In addition to analysis, student will learn mechanism
synthesis using geometric constraint programming techniques. It is important for the student to
learn to use modern computer tools to aid in mechanism design and analysis. The familiarity
with rigid body analysis in modern computer analysis programs will give students a competitive
advantage in the current marketplace. The course covers position, velocity, and acceleration
analysis of planar linkages and slider crank mechanisms. Grashof’s laws, Kutzbach criterion,
and the Newton-Raphson method. Cam profiles and cam follower motion is studied along with
gears and gearing including epicyclic gear trains.
Exam 1: 20%
Exam 2: 20%
Final Exam: 30%
90-100 = A, 80-89 = B, 70-79 = C, 60-69 = D, < 60 = F
The instructor reserves the right to lower the minimum requirements for each letter grade.
Grades will be kept on WebCT and will be updated on a regular basis.
Kinetics:
Mechanism:
Machine:
Linkage:
Link:
Joint:
Skeleton Diagram:
• Wrapping Pairs
– Example: Belt on Pulley (sheave)
–
–
2.4.1. Mobility
Calculating the number of degrees of freedom for a mechanism is determining its mobility.
y
A
θ2
B
O2 θ3 x
%
% %
2.5.2. Inversions
• Every mechanism has a “ground” or “base” or “frame” link that is fixed. See
Cleghorn’s
•
Section
• 1.7
Let’s do an example to find the inversions for a four-bar mechanism where s = 1, l = 8, p = 6,
and q = 6 using the grids shown in Fig. 2.3.
Check Grashof’s Criteria:
6 6
5 5
4 4
3 3
2 2
1 1
0 0
(a) Crank Rocker (b) Crank Rocker
6 6
5 5
4 4
3 3
2 2
1 1
0 0
(c) Double Crank (d) Double Rocker
where n. . .
Consider a loop of three fixed lengths: A, B, and C with angles α , β , and γ .
A
C
γ
α
B
3.1. Loop Closure Equations 21
Let’s define the x-axis to be along length B and the vectors defining the loop to go in a clockwise
direction.
A
C
y x
γ
α
B
Breaking these vectors into components gives the following six expressions. Keep in mind all
angles are defined from the positive x-axis and are positive when measure in a counterclockwise
direction.
Ax :
Ay :
Bx :
By :
Cx :
Cy :
Now all the components in the x direction can be summed to zero and all the y components can
be summed to zero:
x: Ax + Bx +Cx = 0
y: Ay + By +Cy = 0
or
x:
y:
1. Rewrite each component equation to move the unknown and unneeded angle as a single
term on one side:
x:
y:
2. Square each component equation and add them together. Invoking the relationship
sin2 θ + cos2 θ = 1
x:
y:
Sum:
3. Now the unknown angle γ has been eliminated and only the remaining length A is un-
known. Expand the binomial terms:
p
A= C2 + B2 − 2BC cos α
Similarly q
C= A2 + B2 − 2AB cos γ
and q
B= A2 +C2 − 2AC cos β
5. Recall that β = π − α − γ . Take the sine of both angles and recall a trigonometric identity:
sin β = sin(π − α − γ )
sin β = sin(π ) cos(α + γ ) − cos(π ) sin(α + γ )
sin β = 0 − (−1) sin(α + γ )
sin β = sin(α ) cos(γ ) + cos(α ) sin(γ )
6. Notice the right hand side of the equation matches the term from the loop closure equations.
Make the appropriate substitution:
B sin α = A sin β
This is the law of sines. A similar formulation will give the familiar ratios:
Figure 3.1.: A pneumatic cylinder is an example of an inverted slider-crank where the stroke and
angle can change.
y A
r
θ3
θ2 x
%C %B
d
Figure 3.2.: Inverted slider crank. Let a be the fixed length from A to C (the crank length).
3.1. Loop Closure Equations 25
Another way to set up the loop closure equations is to find two different paths to the same point
in the mechanism. This formulation may be a little easier depending on how the angles are
defined. Consider two different loop closure formulations shown in Fig. 3.3 on the next page.
Both formulations describe the same physical system, so they should ultimately produce the
same solution.
3.1. Loop Closure Equations 26
~a ~r
θ3
θ2
x
d~
(a) ~a +~r + d~ = 0
~a ~r
θ3
θ2
x
d~
(b) ~a = d~ +~r
Figure 3.3.: Loop closure formulations for the inverted slider crank shown in Fig. 3.2.
3.1. Loop Closure Equations 27
For the inverted slider crank shown in Fig. 3.2, the following variables are known: a = 0.15 m,
d = 0.20 m and θ2 = 35◦ . The goal of the position analysis is to find r and θ3 .
1. Break the vectors into components
a) For the loop closure equation from Fig. 3.3a on the preceding page
x:
y:
b) For the loop closure equation from Fig. 3.3b on the previous page
x:
y:
2. Rearrange the loop closure equations to have all terms of θ3 on one side. Notice that either
formulation will give the same results.
x:
y:
3. Square each equation and add them together to eliminate θ3 :
5. Solve for r:
7. Compute the arctangent and adjust for the correct quadrant. Use the atan2 function.
The above example is also implemented in Mathematica and the Mathematica notebook can be
downloaded from
http://personal.utulsa.edu/~jeremy-daily/ME3212/InvertedSliderCrank.
nb
3.1. Loop Closure Equations 28
Given: r2 =
r3 =
e=
θ2 =
Find: θ3 and xB
1. Write the Loop equations:
2. Break into components
a) x :
b) y :
3. Rearrange y equation to solve for θ3
r3
θ3 r4
A
y
r2
θ2 θ4
x
%O %O
2 4
r1
3. Determine β :
6. Solve for θ3
7. Solve for θ4
See course website for a Mathematica notebook with the solutions to this example.
http://personal.utulsa.edu/~jeremy-daily/ME3212/4barPositions.nb
3.2. Coupler Curves 32
Point P is. . .
Analysis Steps: See
1. Write down loop closure equations and solve for θ3 . Cleghorn’s
Section
2. Write the equations for the vector describing point P. 4.2
3.2. Coupler Curves 33
5 c l c %c l e a r s c r e e n
c l e a r a l l %c l e a r memory
7 c l o s e a l l %c l o s e a l l f i g u r e windows
9 %i n p u t known v a l u e s
r2 =0.5
11 r 3 = 1 . 2
r4 =1.3
13 e = 0 . 3
15 %d e f i n e t h e t a 2 f r o m 0 t o 360 i n 5 d e g r e e i n c r e m e n t s u s i n g r a d i a n s
t h e t a 2 =[0: pi / 3 6 : 2 ∗ pi ] ;
17
%s o l v e f o r t h e t a 3
19 t h e t a 3 = a s i n ( ( e+ r 2 ∗ s i n ( t h e t a 2 ) ) / r 3 ) ;
21 %s o l v e f o r xB
xB= r 2 ∗ c o s ( t h e t a 2 )+ r 3 ∗ c o s ( t h e t a 3 ) ;
23 yB= z e r o s ( s i z e ( xB ) ) ; %T h i s i s n e e d e d f o r t h e c o u p l e r c u r v e
25 %D e t e r m i n e C o u p l e r c u r v e
xp = r2 ∗ cos ( t h e t a 2 ) − r4 ∗ cos ( t h e t a 3 ) ;
27 yp = e + r 2 ∗ s i n ( t h e t a 2 ) + r 4 ∗ s i n ( t h e t a 3 ) ;
29 f o r n = 1 : l e n g t h ( t h e t a 2 ) %I t e r a t e t h r o u g h t h e v a l u e s o f t h e t a 2
31 %d e f i n e k e y p o i n t s on mechanism i n a
%f a s h i o n t o p l o t t h e c o n f i g u r a t i o n
33 p t s x = [ 0 , r 2 ∗ c o s ( t h e t a 2 ( n ) ) , xB ( n ) , xp ( n ) ] ;
p t s y =[ e , e+ r 2 ∗ s i n ( t h e t a 2 ( n ) ) , 0 , yp ( n ) ] ;
35
%p l o t t h e c o u p l e r c u r v e p o i n t s and one mechanism o r i e n t i a t i o n
3.2. Coupler Curves 34
39 p l o t ( xp , yp , ’ . ’ , p t s x , p t s y , ’−o ’ , xB , yB , ’ y ’ )
41 %m o d i f y t h e p l o t
axis equal
43 x l a b e l ( ’ x−p o s i t i o n (m) ’ )
y l a b e l ( ’ y−p o s i t i o n (m) ’ )
45 t i t l e ( [ ’ C o u p l e r Curve : S l i d e r c r a n k shown w i t h ’ . . .
num2str ( t h e t a 2 ( n ) ∗ 1 8 0 / p i ) ’ d e g r e e i n p u t ’ ] )
47 i f t h e t a 2 ( n )== p i / 3 %f i n d t h e p o s i t i o n when t h e t a 2 i s 60 d e g r e e s
s a v e a s ( F , [ ’ C o u p l e r C u r v e _ ’ num2str ( t h e t a 2 ( n ) ∗ 1 8 0 / p i ) ’ . p d f ’ ] )
49 end
51 %C a p t u r e t h e c u r r e n t p l o t t o make an a n i m a t i o n
s e t ( gca , ’ P l o t B o x A s p e c t R a t i o ’ , [ 5 1 2 384 7 2 0 ] )
53 M( n )= g e t f r a m e ( F , [ 0 0 512 3 8 4 ] ) ;
end
55 %C r e a t e a m o vi e f i l e :
m o v i e 2 a v i (M, ’ O f f s e t S l i d e r C r a n k . a v i ’ )
57
F= f i g u r e %open new f i g u r e window
59 p l o t ( t h e t a 2 ∗ 1 8 0 / pi , xB , ’− ’ )
x l a b e l ( ’ Crank Angle , \ t h e t a _ 2 ( deg ) ’ )
61 y l a b e l ( ’ S l i d e D i s t a n c e , x_B (m) ’ )
t i t l e ( ’ S l i d e p o s i t i o n as a f u n c t i o n of i n p u t angle ’ )
63 g r i d on
axis tight
65 s a v e a s ( F , ’ S l i d e P o s i t i o n . p d f ’ )
The output generated by the above code is shown in Fig. 3.5 on the following page.
3.2. Coupler Curves 35
Coupler Curve: Slider crank shown with 60 degree input Slide position as a function of input angle
1.6
1.5
1.5
1.4
1 1.3
1.2
0.5 1.1
0 0.9
0.8
−0.5
0.7
−1.5 −1 −0.5 0 0.5 1 1.5 0 50 100 150 200 250 300 350
x−position (m) Crank Angle, θ (deg)
2
(a) (b)
Figure 3.5.: Output graph from Matlab code for the coupler curves of an Offset slider crank
mechanism.
Open Closure
A B C D E F G
1 ME!3212:!Mechanisms Plot!Coordinates
2 Four Bar!Linkage!Analysis Upper!Leg
3 Parameter Value Units Vertex X Y
4 r1 10 in 1 0 0
5 r2 3 in 2 1.5000 2.5981
6 r3 12 in 3 4.98362767 12.6957
7 r4 8 in
8 Lower!Leg
9 Radians Degrees Vertex X Y
10 theta2 2.094 120.000 3* 6.6708 7.2744
11 theta3 1.000 57.296 4 10.0000 0.0000
12 theta4 2.000 114.592
13 Difference 1.6872E+00 5.4213E+00
14 residual 5.6778E+00
15
16
14
17
18
12
19
20
21 10
22
23 8
24 Upper!Leg
25 6 Lower!Leg
26
27
4
28
29
2
30
31
32 0
33 4 2 0 2 4 6 8 10 12
34
3.2. Coupler Curves 38
3.2. Coupler Curves 39
A B C D E F G
1 ME!3212:!Mechanisms Plot!Coordinates
2 Four Bar!Linkage!Analysis Upper!Leg
3 Parameter Value Units Vertex X Y
4 r1 10 in 1 0 0
5 r2 3 in 2 1.5000 2.5981
6 r3 12 in 3 9.233836209 7.9632
7 r4 8 in
8 Lower!Leg
9 Radians Degrees Vertex X Y
10 theta2 2.094 120.000 3* 9.2338 7.9632
11 theta3 0.464 26.557 4 10.0000 0.0000
12 theta4 1.667 95.496
13 Difference 3.7786E 06 7.8999E 06
14 residual 8.7571E 06
15
16
9
17
18 8
19
20 7
21
6
22
23 5
24 Upper!Leg
25 4 Lower!Leg
26
27 3
28 2
29
30 1
31
32 0
33 4 2 0 2 4 6 8 10 12
34
3.2. Coupler Curves 40
A B C D E F G
1 ME!3212:!Mechanisms Plot!Coordinates
2 Four Bar!Linkage!Analysis Upper!Leg
3 Parameter Value Units Vertex X Y
4 r1 10 in 1 0 0
5 r2 3 in 2 =B5*COS(B10) =B5*SIN(B10)
6 r3 12 in 3 =F5+B6*COS(B11) =G5+B6*SIN(B11)
7 r4 8 in
8 Lower!Leg
9 Radians Degrees Vertex X Y
10 theta2 =2*PI()/3 =DEGREES(B10) 3* =F11+B7*COS(B12) =B7*SIN(B12)
11 theta3 0.4635153062=DEGREES(B11) 4 =B4 0
12 theta4 1.6667142836=DEGREES(B12)
13 Difference =F6 F10 =G6 G10
14 residual =SQRT(F13^2+G13^2)
15
16
9
17
18 8
19
20 7
21
6
22
23 5
24 Upper!Leg
25 4 Lower!Leg
26
27 3
28
2
29
30 1
31
32 0
33 4 2 0 2 4 6 8 10 12
34
Cross Closure
3.2. Coupler Curves 41
A B C D E F G
1 ME!3212:!Mechanisms Plot!Coordinates
2 Four Bar!Linkage!Analysis Upper!Leg
3 Parameter Value Units Vertex X Y
4 r1 10 in 1 0 0
5 r2 3 in 2 1.5000 2.5981
6 r3 12 in 3 9.030990743 3.1550
7 r4 8 in
8 Lower!Leg
9 Radians Degrees Vertex X Y
10 theta2 2.094 120.000 3* 14.3224 6.7318
11 theta3 0.500 28.648 4 10.0000 0.0000
12 theta4 1.000 57.296
13 Difference 5.2914E+00 3.5767E+00
14 residual 6.3869E+00
15
16
4
17
18
19 2
20
21
22 0
23 5 0 5 10 15 20
24 Upper!Leg
2
25 Lower!Leg
26
27 4
28
29
30 6
31
32
33 8
34
3.2. Coupler Curves 42
A B C D E F G
1 ME!3212:!Mechanisms Plot!Coordinates
2 Four Bar!Linkage!Analysis Upper!Leg
3 Parameter Value Units Vertex X Y
4 r1 10 in 1 0 0
5 r2 3 in 2 1.5000 2.5981
6 r3 12 in 3 5.884873205 6.8604
7 r4 8 in
8 Lower!Leg
9 Radians Degrees Vertex X Y
10 theta2 2.094 120.000 3* 5.8849 6.8604
11 theta3 0.908 52.019 4 10.0000 0.0000
12 theta4 2.111 120.957
13 Difference 1.8606E 06 5.1288E 07
14 residual 1.9300E 06
15
16
4
17
18
19 2
20
21
22 0
23 4 2 0 2 4 6 8 10 12
24 Upper!Leg
2
25 Lower!Leg
26
27 4
28
29
30 6
31
32
33 8
34
3.2. Coupler Curves 43
Given: r2 , xo , R, θ2
Find: r and θ3
Write down the Loop closure equations:
x:
y:
Rearrange to eliminate r:
This equation has just 1 unknown θ3 but no analytical solution. Therefore, we’ll solve using the
Newton-Raphson method.
3.4. Newton-Raphson Method 47
where
n is
For this example,
f=
f′ =
The algorithm is as follows:
1. Initialize the criteria: let ∆θ3 = 1
(1)
2. Guess a value for θ3
3. Do while |∆θ3 | > 0.000175
(n) (n)
4. Calculate f (θ3 ) and f ′ (θ3 )
5. ∆θ3 = − f / f ′
(n+1) (n)
6. θ3 = θ 3 + ∆θ 3
7. Increment n by 1 and loop back to step 3.
8. Continue to complete the solution by solving for
a) r =
or
r=
b) x =
9. Increment θ2 and repeat.
3.4. Newton-Raphson Method 48
In Matlab, the algorithm for the example looks like the following:
1 %ME 3 2 1 2 : Mechanisms
%Newton−Rahpson Method f o r s o l v i n g l o o p c l o s u r e e q u a t i o n s
3 %Dr . Jer em y D a i l y
%
5 clc ; clear all ; close al l
7 %knowns :
r 2 = 1 . 2 %m e t e r s
9 x0 = 1 . 8
R= . 4
11
t h e t a 2 = [ 0 : p i / 1 0 1 : 4 . ∗ p i ] ; %Two c y c l e s o f t h e c r a n k
13
t h e t a 3 = 0 ; %g u e s s a s o l u t i o n f o r f i r s t i t e r a t i o n
15 f o r n = 1 : l e n g t h ( t h e t a 2 )
d e l t a t h e t a 3 = 1 ; %s e t t h e l o o p c o n d i t i o n a l
17 %g u e s s t h e same s o l u t i o n a s t h e t i m e b e f o r e
19 w h i l e norm ( d e l t a t h e t a 3 ) > 0 . 0 0 0 1 7 5 %0 . 0 1 d e g r e e s
f = ta n ( t h e t a 3 ( n ) ) . ∗ ( x0+R∗ t h e t a 3 ( n)− r 2 ∗ c o s ( t h e t a 2 ( n ) ) ) . . .
21 −r 2 ∗ s i n ( t h e t a 2 ( n ) ) ;
d f d t h e t a 3 = ( 1 . / ( c o s ( t h e t a 3 ( n ) ) . ^ 2 ) ) . ∗ ( x0+R∗ t h e t a 3 ( n ) . . .
23 −r 2 ∗ c o s ( t h e t a 2 ( n ) ) ) + R∗ c o s ( t h e t a 3 ( n ) ) ;
d e l t a t h e t a 3 =−f / d f d t h e t a 3 ;
25 %u p d a t e t h e s o l u t i o n u n t i l i t i s w i t h i n t o l e r a n c e
t h e t a 3 ( n )= t h e t a 3 ( n )+ d e l t a t h e t a 3 ;
27 end
t h e t a 3 ( n +1)= t h e t a 3 ( n ) ;
29
end
31
%S i n c e t h e l a s t command i n t h e l o o p c r e a t e d an e x t r a t h e t a 3 e n t r y
33 %we m u s t remove i t by a s s i g n i n g i t t o t h e em p t y s e t
theta3 (n +1)=[];
35
%S o l v e f o r t h e r e m a i n i n g unkowns
37 r _ s i n e s = r 2 ∗ s i n ( t h e t a 2 ) . / s i n ( t h e t a 3 ) ; %law o f s i n e s
r _ c o s i n e s = s q r t ( r 2 . ^ 2 + ( x0+R∗ t h e t a 3 ) . ^ 2 . . .
39 − 2∗ r 2 . ∗ ( x0+R . ∗ t h e t a 3 ) . ∗ c o s ( t h e t a 2 ) ) ; %law o f c o s i n e s
x=R∗ t h e t a 3 ;
3.4. Newton-Raphson Method 49
41 p l o t ( t h e t a 2 , t h e t a 3 , ’ : ’ , t h e t a 2 , r _ s i n e s , ’− ’ , t h e t a 2 , r _ c o s i n e s , . . .
’−− ’ , t h e t a 2 , x , ’ −. ’ )
43 x l a b e l ( ’ \ t h e t a _ 2 ( r a d i a n s ) ’ )
ylabel ( ’ Position variables ’ )
45 l e g e n d ( ’ \ t h e t a _ 3 ( r a d ) ’ , ’ r : s i n e s (m) ’ , ’ r : c o s i n e s (m) ’ , ’ x (m) ’ )
3.5
θ3 (rad)
r:sines (m)
3
r:cosines (m)
x (m)
2.5
2
Position variables
1.5
0.5
−0.5
−1
−1.5
0 2 4 6 8 10 12 14
θ2 (radians)
Figure 3.6.: Output of Newton-Raphson example. What are the spikes from?
3.4. Newton-Raphson Method 50
1 F= f i g u r e ( 1 ) ;
p l o t ( t h e t a 2 ∗ 1 8 0 / pi , t h e t a 3 ∗ 1 8 0 / pi , ’ . ’ , . . .
3 t h e t a 2 ∗ 1 8 0 / pi , t h e t a 3 _ o p e n ∗ 1 8 0 / pi , ’ s ’ )
%No t e : t h e t a 3 _ o p e n m u s t be d e f i n e d i n o r d e r f o r t h i s t o work
5 l e g e n d ( ’ Newton−Raphson ’ , ’ A n a l y t i c a l ’ , ’ L o c a t i o n ’ , ’ N o rt h Wes t ’ )
x l a b e l ( ’ I n p u t Crank Angle , \ t h e t a _ 2 [ deg ] ’ )
7 y l a b e l ( ’ C o u p l e r Angle , \ t h e t a _ 3 [ deg ] ’ )
axis tight
9 s e t ( gca , ’ XTick ’ , [ 0 : 3 0 : 3 6 0 ] )
g r i d on
11 t i t l e ( ’ C o u p l e r a n g l e i n t h e open c l o s u r e c o n f i g u r a t i o n ’ )
3.5. Homework Problem Set 3 52
r2
r R
e
θ2 O3 θ3 x
O2
% %
r1
3. Perform the “Advanced Design” SolidWorks Tutorial. Turn in a print of your completed part.
3.5. Homework Problem Set 3 53
3.5. Homework Problem Set 3 54
55
50
Coupler Angle, θ [deg]
45
3
40
35
30
25
20
−20 Newton−Raphson
Analytical
−25
−30
−35
Coupler Angle, θ3 [deg]
−40
−45
−50
−55
−60
−65
0 30 60 90 120 150 180 210 240 270 300 330 360
Input Crank Angle, θ [deg]
2
Figure 3.7.: Graphs for the Newton-Raphson and analytical solution for θ3
3.6. Multi Loop Mechanisms 55
There are four equations after breaking up the loop closure equations:
x1 :
y1 :
x2 :
y2 :
3.7. Toggle and Limit Positions 56
Bottom center:
Number Synthesis:
Dimensional Synthesis:
Classical Analysis:
5. Create 5 3-point arcs. They do not have to be the same size yet and placement is arbitrary
at this time. These arcs will represent the scoop.
4.1. Geometric Constraint Programming 63
6. Make all the arcs equal by selecting them while holding the Shift key. Then click on the
Equal button under the Add Relations Pane.
4.1. Geometric Constraint Programming 64
9. Position and dimension the arcs in appropriate locations. The position will describe loca-
tions along the path of travel. Therefore, the arc must be facing up at the bottom location
and tilted at the drop off location. The other points are used to guide the path of travel.
Start by making the right edge of one of the arcs coincident with the origin and fixing its
location.
4.1. Geometric Constraint Programming 66
10. Draw congruent triangles from each of the arc endpoints. To ensure congruency, make all
corresponding legs of the triangle equal. The tip of the triangle will represent the connec-
tion to one of the links.
4.1. Geometric Constraint Programming 67
11. Draw a perimeter circle through three of the five tips of the triangles.
12. Make the points coincident with the circle.
4.1. Geometric Constraint Programming 68
13. Drag the other two points close to the circle and make them coincident. All five points
should be coincident with the circle.
4.1. Geometric Constraint Programming 69
14. Create another set of congruent triangles to represent the other connecting point for the
coupler.
4.1. Geometric Constraint Programming 70
15. Draw the other circle and maneuver the locations of the arc and the triangles to get a com-
pact package. This may be difficult and you may have to delete some fixed constraints.
4.1. Geometric Constraint Programming 71
16. Once the position is where you like, Select All and make a Block (Tools ⊲ Make Block).
4.1. Geometric Constraint Programming 72
17. Draw the frame (link 1) by connecting the two circle centers. Fix the ends. The draw
moving links 2 and 4. Trace the coupler. Add dimensions to keep the lengths fixed.
4.1. Geometric Constraint Programming 73
All dimensions should be determined and the mechanism synthesis is completed. The mecha-
nism can be actuated by a motor for continuous rotation or by a cylinder or cam to rock back and
forth.
4.2. Homework Problem Set 5 74
4. Measure Magnitudes
x
0v reference
B
C
y r5
r6
r4
r3
θ2 θ4
x
O2%
O4%
r1
θ̇2
θ3
A
Figure 5.1.: Four-Bar Mechanism Example
Use the space allocated in Fig. 5.2 to construct the velocity polygon for this example. Recall
velocity equivalence:
~vB =~vA +~vB/A =~vO4 +~vB/O4
1. Determine the angles for the lines of action.
2. Draw v~A
5. Find the intersection of the vB construction line and the vB/A construction line. Draw the
vector representing ~vB/A and measure its length.
5.3. Graphical Analysis 81
7. Velocity of C:
x
0v reference
• The ratio...
A
y r
θ̇2 r2
θ2
θ4 O2 x
%O %
4
r1
where r1 =
r2 =
r=
θ2 =
θ4 =
θ̇2 =
1. Loop Equations:
a) x :
b) y :
2. Differentiate with respect to time:
a)
b)
5.4. Analytical Analysis 84
We can write a computer program to solve the above system for various angles of θ2 and a
constant angular velocity. An example in Matlab is as follows:
1 %ME 3 2 1 2 : Mechanisms
%I n v e r t e d S l i d e r Crank
3 %F i n d i n g t h e v e l o c i t i e s o f l i n k s
%Dr . Jer em y D a i l y
5
%c l o s e an i n i t i a l i z e s y s t e m
7 clc
5.4. Analytical Analysis 85
clear all
9 close all
11 %I n v e r t e d s l i d e r c r a n k mechanism
%knowns :
13 r 1 =45
r 2 =25
15 omega2 =20
17 t h e t a 2 = l i n s p a c e ( 0 , 2 ∗ pi , 1 5 0 ) ;
19 %p o s i t i o n s o l u t i o n :
r = s q r t ( r 1 ^2+ r 2 ^2−2∗ r 1 ∗ r 2 ∗ c o s ( pi −t h e t a 2 ) ) ;
21 t h e t a 4 = a c o s ( ( r . ^ 2 + r 1 . ^ 2 − r 2 . ^ 2 ) . / ( 2 ∗ r ∗ r 1 ) ) ;
t h e t a 4 = a ta n2 ( r 2 ∗ s i n ( t h e t a 2 ) , r 1 + r 2 ∗ c o s ( t h e t a 2 ) ) ;
23
%v e l o c i t y a n a l y s i s
25 f o r i = 1 : l e n g t h ( t h e t a 2 )
A=[ c o s ( t h e t a 4 ( i ) ) −r ∗ s i n ( t h e t a 4 ( i ) ) ;
27 sin ( theta4 ( i )) r ∗ cos ( t h e t a 4 ( i ) ) ] ;
C=[ −r 2 ∗ omega2 ∗ s i n ( t h e t a 2 ( i ) ) r 2 ∗ omega2 ∗ c o s ( t h e t a 2 ( i ) ) ] ’ ;
29 v=A \ C ;
r d o t ( i )= v ( 1 ) ;
31 omega4 ( i )= v ( 2 ) ;
end
33
%P l o t t h e r e s u l t s
35 F= f i g u r e ( 1 )
subplot (2 ,2 ,1)
37 p l o t ( t h e t a 2 ∗ 1 8 0 / pi , t h e t a 4 )
x l a b e l ( ’ \ t h e t a _ 2 ( deg ) ’ )
39 y l a b e l ( ’ \ t h e t a _ 4 ( r a d ) ’ )
g r i d on
41 a x i s t i g h t
%s e t s t h e c u r r e n t a x i s t o h a ve t i c k s e v e r y 45 deg
43 s e t ( gca , ’ X t i c k ’ , [ 0 : 4 5 : 3 6 0 ] )
45 s u b p l o t ( 2 , 2 , 2 )
p l o t ( t h e t a 2 ∗ 1 8 0 / pi , r )
47 x l a b e l ( ’ \ t h e t a _ 2 ( deg ) ’ )
y l a b e l ( ’ s l i d e r p o s i t i o n (mm) ’ )
49 g r i d on
5.4. Analytical Analysis 86
axis tight
51 s e t s t h e c u r r e n t a x i s t o h av e t i c k s e v e r y 45 deg
s e t ( gca , ’ X t i c k ’ , [ 0 : 4 5 : 3 6 0 ] )
53
subplot (2 ,2 ,3)
55 p l o t ( t h e t a 2 ∗ 1 8 0 / pi , omega4 )
x l a b e l ( ’ \ t h e t a _ 2 ( deg ) ’ )
57 y l a b e l ( ’ \ omega_4 ( r a d / s ) ’ )
g r i d on
59 a x i s t i g h t
s e t ( gca , ’ X t i c k ’ , [ 0 : 4 5 : 3 6 0 ] )
61
subplot (2 ,2 ,4)
63 p l o t ( t h e t a 2 ∗ 1 8 0 / pi , r d o t )
x l a b e l ( ’ \ t h e t a _ 2 ( deg ) ’ )
65 y l a b e l ( ’ s l i d e r v e l o c i t y (mm/ s ) ’ )
g r i d on
67 a x i s t i g h t
s e t ( gca , ’ X t i c k ’ , [ 0 : 4 5 : 3 6 0 ] )
69
saveas (F , ’ In v ert ed Sl i d erCran k V el o c i t y . pdf ’ )
5.4. Analytical Analysis 87
70
0.5
60
0
4
40
30
−0.5
0 45 90 135 180 225 270 315 360 0 45 90 135 180 225 270 315 360
θ (deg) θ (deg)
2 2
6 400
0 0
4
−2
−200
−4
−6 −400
0 45 90 135 180 225 270 315 360 0 45 90 135 180 225 270 315 360
θ2 (deg) θ2 (deg)
Figure 5.3.: Output of the velocity analysis program of an inverted slider crank.
A
r
r2
r3
θ2
x θ̇2 O2 θ3 O4
r1 B
r3
r5
θ3 φ
A
r6
r2
C
r4
y
θ̇2
θ2 θ4
O4 x
O2% %
r1
Ans: θ̇4 = 45.51 rad/s, |vc | = 9.025 m/s and ∠vc = 137.9◦
3. Determine the velocities of points B, C, and D along with the angular velocity of the con-
necting rod (△ABC) in a double slider crank where r2 = 2 in., r3 = 10 in, CA = 4 in,
CD = 8 in., θ2 = −120◦ , θ˙2 = ω2 = 42 rad/s cw, and O2 is at (−3, 0).
5.5. Homework Problem Set 6 92
10
0
Ans: |vc | = 67.9 in/s and ∠vc = 154◦
6. Acceleration Analysis
Graphical techniques are discussed in Cleghorn’s Chapter 3. This section is focused on analytical
techniques. Simply put, the acceleration is the time derivative of the velocity equations.
Accelerations are needed to determine forces: F ~ = m~a
y:
θ̈2 =
Substitute numbers:
6.2. Inverted Slider Crank 94
vc,y =
ac,y =
Coriolis acceleration:
6.3. Homework Problem Set 7 96
A
r
a2 θ2
a3
O4 x
θ̇2 O2
θ3
d B
where a2 = 75 mm, a3 = 400 mm, d = 125 mm, r = 193.6 mm, θ2 = 150◦ , θ3 = 11.17◦ ,
and θ̇2 = 60rad/sec (constant). The position equations are as follows:
a2 cos θ2 + r cos θ3 =d
a2 sin θ2 − r sin θ3 =0
a3
a5
θ3 φ
A
a6
a2
C
a4
y
θ̇2
θ2 θ4
O4 x
O2% %
d
where a2 = 150 mm, a3 = 250 mm, a4 = 250 mm, d = 100 mm, a5 = 150 mm, a6 = 300
mm, θ2 = 105◦ , θ3 = 20.13◦, θ4 = 67.47◦, φ = 29.93◦, θ̇2 = ω2 = 56 rad/sec (constant),
θ̇3 = 27.84 rad/s, and θ˙4 = 45.51 rad/s. The position equations are as follows:
a2 cos θ2 + a3 cos θ3 − a4 cos θ4 =d
a2 sin θ2 + a3 sin θ3 − a4 sin θ4 =0
and the velocity equations are:
−a2 θ̇2 sin θ2 − a3 θ̇3 sin θ3 + a4 θ̇4 sin θ4 =0
a2 θ̇2 cos θ2 + a3 θ̇3 cos θ3 − a4 θ̇4 cos θ4 =0
3. Perform the “SolidWorks Motion” tutorial. Turn in a print of the time history of the contact
force.
7. Cams
A cam transmits rotary motion to another link by direct contact.
Vocabulary:
• Cam:
• Follower:
• Cam Profile:
• Trace Point:
• Pitch Curve:
• Base Circle:
• Prime Circle:
• Pressure Angle:
• Rise:
• Dwell:
• Fall:
•
•
•
0
7.1. Types of Cam Followers 101
0
7.2. Cam Follower Motion 102
7.2.1. Displacement
y = f (θ )
7.2.2. Velocity
v = ẏ =
7.2.3. Acceleration
a = ÿ =
7.2.4. Jerk
...
jerk = y =
7.3. Cam Follower Profiles 103
Interval Initial Angle Final Angle Type of Motion Initial Displacement Final Displacement
Parabolic 0 A
Const. Velocity A B
Parabolic B C
Dwell C C
Parabolic C D
Parabolic D 0
Fundamental Equations: See
Cleghorn’s
Section
7.5.2 and
7.5.3
7.3. Cam Follower Profiles 104
f (θ )
60
40
20
θ
30 60 90 120 150 180 210 240 270 300 330 360
f ′ (θ )
f ′′ (θ )
Interval Initial Angle Final Angle Type of Motion Initial Displacement Final Displacement
Harmonic 0 A
Const. Velocity A B
Harmonic B C
Dwell C C
Harmonic C 0
Fundamental Equations: See
Cleghorn’s
Section
7.5.4
7.3. Cam Follower Profiles 107
f (θ )
60
40
20
θ
30 60 90 120 150 180 210 240 270 300 330 360
f ′ (θ )
f ′′ (θ )
Interval Initial Angle Final Angle Type of Motion Initial Displacement Final Displacement
Cycloidal 0 A
Dwell A A
Cycloidal A B
Dwell B B
Cycloidal B 0
Fundamental Equations: See
Cleghorn’s
Section
7.5.5
7.3. Cam Follower Profiles 110
f (θ )
60
40
20
θ
30 60 90 120 150 180 210 240 270 300 330 360
f ′ (θ )
f ′′ (θ )
6
X = 130 cm, Y = 100 cm, e = 50 cm, rb = 70 cm See
Cleghorn’s
Example
5 7.1
0.30
0.20
0.10
θ
30 60 90 120 150 180 210 240 270 300 330 360
7.4. Cam Design 113
b) yc :
b) y′c :
a) xc =
b) yc =
6. Determine the radial length and orientation
a) rc =
b) γc =
7. Determine the angle with respect to the disk cam:
7.4. Cam Design 115
100
150 30
50
180 0
210 330
240 300
270
Figure 7.1.: Cam profile in polar coordinates. The inner circle shows the base radius and the
outer circle shows the maximum dimensions.
7.4. Cam Design 116
10 c l c %c l e a r s t h e command window
c l e a r a l l %c l e a r s memory
12 c l o s e a l l %c l o s e s a l l f i g u r e windows
14 %C o n s t r u c t cam f o l l o w e r p r o f i l e
16 %Harmonic r i s e f r o m 0 t o 100 d e g r e e s ( e v e r y d e g r e e )
L=.4
18 b e t a =100∗ p i / 1 8 0
t h e t a = [ 0 : p i / 1 8 0 : b e t a ] ; %m u s t be i n r a d i a n s
20 f =L∗(1 − c o s ( p i ∗ t h e t a / b e t a ) ) / 2 ; %f o l l o w e r d i s p l a c e m e n t
f p r i m e = p i ∗L∗ s i n ( p i ∗ t h e t a / b e t a ) / ( 2 ∗ b e t a ) ; %f o l l o w e r m o t i o n s l o p e
22 f d o u b l e p r i m e = p i ^2∗L∗ c o s ( p i ∗ t h e t a / b e t a ) / 2 / b e t a ^ 2 ; %f o l l o w e r a c c e l
30 %b u i l d t h e f o l l o w e r m o t i o n p r o f i l e s by c o n c a t e n a t i n g f u n c t i o n s
f =[ f s s t a r ] ;
32 f p r i m e =[ f p r i m e f p r i m e s t a r ] ;
f d o u b l e p r i m e =[ f d o u b l e p r i m e z e r o s ( s i z e ( s s t a r ) ) ] ;
34
%C y c l o i d a l F a l l f r o m 160 t o 300 d e g r e e s
36 L2 = . 1
L=L−L2
38 b e t a =140∗ p i / 1 8 0
t h e t a s t a r = [ 0 : p i / 1 8 0 : b e t a ] ; %m u s t be i n r a d i a n s
40 s s t a r =L2+L∗(1 − t h e t a s t a r / b e t a + s i n ( 2 ∗ p i ∗ t h e t a s t a r / b e t a ) / ( 2 ∗ p i ) ) ;
7.4. Cam Design 117
f p r i m e s t a r = −L / b e t a ∗(1 − c o s ( 2 ∗ p i ∗ t h e t a s t a r / b e t a ) ) ;
42 f d o u b l e p r i m e s t a r =−2∗ p i ∗L∗ s i n ( 2 ∗ p i ∗ t h e t a s t a r / b e t a ) / b e t a ^ 2 ;
t h e t a =[ t h e t a t h e t a ( end )+ t h e t a s t a r ] ;
44
%b u i l d t h e f o l l o w e r m o t i o n p r o f i l e s by c o n c a t e n a t i n g f u n c t i o n s
46 f =[ f s s t a r ] ;
f p r i m e =[ f p r i m e f p r i m e s t a r ] ;
48 f d o u b l e p r i m e =[ f d o u b l e p r i m e f d o u b l e p r i m e s t a r ] ;
50 %P a r o b o l i c f a l l f r o m 300 t o 360 d e g r e e s
L=L2
52 b e t a =60∗ p i / 1 8 0
54 %p a r t 1
t h e t a s t a r = [ 0 : p i / 1 8 0 : b e t a / 2 ] ; %m u s t be i n r a d i a n s
56 s s t a r = L∗(1 −2∗( t h e t a s t a r / b e t a ) . ^ 2 ) ;
f p r i m e s t a r =−4∗L∗ t h e t a s t a r / b e t a ^ 2 ;
58 f d o u b l e p r i m e s t a r =−4∗L / b e t a ^2∗ o n e s ( s i z e ( t h e t a s t a r ) ) ;
60 %b u i l d t h e f o l l o w e r m o t i o n p r o f i l e s by c o n c a t e n a t i n g f u n c t i o n s
f =[ f s s t a r ] ;
62 f p r i m e =[ f p r i m e f p r i m e s t a r ] ;
f d o u b l e p r i m e =[ f d o u b l e p r i m e f d o u b l e p r i m e s t a r ] ;
64
%p a r t 2
66 t h e t a s t a r =[ b e t a / 2 : p i / 1 8 0 : b e t a ] ; %m u s t be i n r a d i a n s
s s t a r = L∗(2 −4∗( t h e t a s t a r / b e t a ) + 2 ∗ ( t h e t a s t a r / b e t a ) . ^ 2 ) ;
68 f p r i m e s t a r =−4∗L / b e t a + 4∗L∗ t h e t a s t a r / b e t a ^ 2 ;
f d o u b l e p r i m e s t a r =4∗L / b e t a ^2∗ o n e s ( s i z e ( t h e t a s t a r ) ) ;
70
%remove r e p e a t e d e l e m e n t t o e n s u r e t h e a r r a y s a r e t h e same s i z e
72 s s t a r ( 1 ) = [ ] ;
fprimestar (1)=[];
74 f d o u b l e p r i m e s t a r ( 1 ) = [ ] ;
76 %b u i l d t h e f o l l o w e r m o t i o n p r o f i l e s by c o n c a t e n a t i n g f u n c t i o n s
f =[ f s s t a r ] ;
78 f p r i m e =[ f p r i m e f p r i m e s t a r ] ;
f d o u b l e p r i m e =[ f d o u b l e p r i m e f d o u b l e p r i m e s t a r ] ;
80
t h e t a s t a r = [ 0 : pi / 1 8 0 : beta ] ;
82 t h e t a =[ t h e t a t h e t a s t a r + t h e t a ( end ) ] ;
7.4. Cam Design 118
84 %P l o t t h e cam f o l l o w e r p r o f i l e
F= f i g u r e ( 1 )
86 s u b p l o t ( 3 , 1 , 1 )
p l o t ( t h e t a ∗ 1 8 0 / pi , f )
88 a x i s t i g h t
g r i d on
90 x l a b e l ( ’Cam a n g l e \ t h e t a , deg ’ )
ylabel ( ’ f ( t h e t a ) , rad ’ )
92 t i t l e ( ’ K i n e m a t i c c o e f f i c i e n t s f o r Cam F o l l o w e r R o t a t i o n ’ )
subplot (3 ,1 ,2)
94 p l o t ( t h e t a ∗ 1 8 0 / pi , f p r i m e )
axis tight
96 g r i d on
x l a b e l ( ’Cam a n g l e \ t h e t a , deg ’ )
98 y l a b e l ( ’ f ^ \ p r i m e ( \ t h e t a ) ’ )
subplot (3 ,1 ,3)
100 p l o t ( t h e t a ∗ 1 8 0 / pi , f d o u b l e p r i m e )
axis tight
102 g r i d on
x l a b e l ( ’Cam a n g l e \ t h e t a , deg ’ )
104 y l a b e l ( ’ f ^ { \ p r i m e \ p r i m e } ( \ t h e t a ) ’ )
s a v e a s ( F , ’ Cam Fo l l o w erM o t i o n . e p s ’ )
106
%Now t h a t t h e cam f o l l o w e r k i n e m a t i c p r o f i l e s a r e known ,
108 %we can d e s i g n t h e d i s k p r o f i l e t h a t g e n e r a t e s t h a t m o t i o n .
126 r c = s q r t ( xc . ^ 2 + yc . ^ 2 ) ;
gammac= a ta n2 ( yc , xc ) ;
128
%d e t e r m i n e t h e l o c a t i o n o f t h e r a d i a l i n r e l a t i o n t o a
130 %r e f e r e n c e on t h e r o t a t i n g cam :
b e t a =gammac−t h e t a ;
132
134 b i g R a d i u s =max ( r c )
136 F= f i g u r e ( 2 )
p o l a r ( beta , r c )
138 ho l d on
p o l a r ( beta , o n e s ( s i z e ( b e t a ) ) ∗ rb , ’ g ’ )
140 p o l a r ( beta , o n e s ( s i z e ( b e t a ) ) ∗ b i g R a d i u s , ’ r ’ )
ho l d o f f
142 t i t l e ( ’Cam p r o f i l e i n p o l a r form ’ )
saveas (F , ’ CamProfile_polar . eps ’ )
144
F= f i g u r e ( 3 )
146 p l o t ( [ 0 r c . ∗ c o s ( b e t a ) ] , . . .
[ 0 r c . ∗ s i n ( b e t a ) ] , r b ∗ c o s ( b e t a ) , r b ∗ s i n ( b e t a ) , ’−− ’ , . . .
148 bigRadius ∗ cos ( beta ) , bigRadius ∗ s i n ( beta ) , ’ : ’ )
ho l d on
150 i =1
p l o t ( [ X X−e ∗ s i n ( a l p h a ( i ) ) xc ( i ) . . .
152 X−e ∗ s i n ( a l p h a ( i )) − d∗ c o s ( a l p h a ( i ) ) ] ,
[Y Y−e ∗ c o s ( a l p h a ( i ) ) yc ( i ) . . .
154 Y−e ∗ c o s ( a l p h a ( i ) ) + d∗ s i n ( a l p h a ( i ) ) ] , ’−ok ’ )
axis equal
156 ho l d o f f
x l a b e l ( ’ x ( cm ) ’ )
158 y l a b e l ( ’ y ( cm ) ’ )
saveas (F , ’ CamSystemConfiguration . eps ’ )
160
F= f i g u r e ( 4 )
162 a x i s ([ −100 160 −120 1 8 0 ] )
axis equal
164 v= a x i s
%P l o t t h e m o t i o n o f t h e cam
166 f o r i = 1 : l e n g t h ( t h e t a )
7.4. Cam Design 120
p l o t ( [ 0 r c . ∗ cos ( beta + t h e t a ( i ) ) ] , . . .
168 [0 r c . ∗ s i n ( beta + t h e t a ( i ) ) ] , . . .
[X X−e ∗ s i n ( a l p h a ( i ) ) xc ( i ) . . .
170 X−e ∗ s i n ( a l p h a ( i )) − d ∗ c o s ( a l p h a ( i ) ) ] , . . .
[Y Y−e ∗ c o s ( a l p h a ( i ) ) yc ( i ) . . .
172 Y−e ∗ c o s ( a l p h a ( i ) ) + d ∗ s i n ( a l p h a ( i ) ) ] , ’−ok ’ )
axis (v)
174 t i t l e ( ’ S i m u l a t e d Cam Motion ’ )
x l a b e l ( ’ x ( cm ) ’ )
176 y l a b e l ( ’ y ( cm ) ’ )
M( i )= g e t f r a m e ( ) ;
178 end
%make a m o vi e
180 m o v i e 2 a v i (M, ’ CamMotion . a v i ’ )
The graphs generated with this program are shown in Figs. 7.2-7.3
A movie of this cam system can be seen at:
http://personal.utulsa.edu/~jeremy-daily/ME3212/CamMotion.avi
7.4. Cam Design 121
0.3
f(theta), rad
0.2
0.1
0
0 50 100 150 200 250 300 350
Cam angle θ, deg
0.2
f′(θ)
−0.2
0 50 100 150 200 250 300 350
Cam angle θ, deg
0.5
f′′(θ)
−0.5
0 50 100 150 200 250 300 350
Cam angle θ, deg
100
80
60
40
20
y (cm)
−20
−40
−60
−80
−100
−100 −50 0 50 100
x (cm)
Figure 7.3.: Cam Follower Motion Profiles for the example. The dashed line represents the base
circle. The cam rotates counterclockwise.
7.5. Homework Problem Set 8 123
8.1. Introduction
b b
8.8.1. Advantages
1.
2.
3.
4.
5.
b b b
~rP :
~rA/P :
~rB :
~rB/P :
8.8. Planetary Gear Trains 136
b b
See
1. Cleghorn’s
Section
4.4.2
2.
3.
8.8. Planetary Gear Trains 137
Example:
Ns = 34, NR = 78, N p = 22, the sun gear is the input, the ring gear is fixed, and the planet carrier
is the output. Build and analysis table:
Component Gear 1 (Sun) ω1 Crank (Carrier) ωarm Gear 3 (Planet) ω3 Gear 4 (ring) ω4
All Turn at x
rpm
Crank is
Fixed, Sun
turns at y
Absolute
Speeds
Example:
Ns = 34, NR = 78, N p = 22, the sun gear is fixed, the ring gear is the output, and the planet carrier
is the input. Use the same analysis table as before.
8.8. Planetary Gear Trains 138
-1
-2
-3 See
-3 -2 -1 0 1 2 3 4 5 6 7 8 9 10 11 12 Cleghorn’s
Example: an electric winch gear train where NS1 = 34, NP1 = 22, NS2 = 26, NP2 = 30, and
Section
NR = 78.
6.5
Table 1: First Stage
Component Gear 1 (Sun) ω1 Crank (Carrier) ωarm Gear 3 (Planet) ω3 Gear 4 (ring) ω4
All Turn at x
rpm
Crank is
Fixed, Sun
turns at y
Absolute
Speeds
Component Gear 1 (Sun) ω1 Crank (Carrier) ωarm Gear 3 (Planet) ω3 Gear 4 (ring) ω4
All Turn at x
rpm
Crank is
Fixed, Sun
turns at y
Absolute
Speeds
8.8. Planetary Gear Trains 140
14
12
10
0
-4 -2 0 2 4 6 8 10
Figure 8.1.: Output from the code listing to plot the configuration of a planetary gear train.
8.9. Differential Gear Trains 142
See
Cleghorn’s
Section
6.5
8.9. Differential Gear Trains 143
Example:When a vehicle is in a turn, determine a) the speed of each wheel and b) the angular
velocity of the planet carrier, c) the driveshaft speed.
b
8.9. Differential Gear Trains 144
A. Mechanisms Design Project
A computer mechanism design project will be given in class. You will have the choice of desiging
one of the following:
• Moving batters tee
• Runner friendly stroller
• Air-bag module near-deployment mechanism
• Remote motorcycle braking system
• Scissor lift safety mechanism
Figure A.1.: Flanged Sleeve Bearing description from http://www.mcmaster.com. This is an oil
impregnated bronze bearing that does not need lubrication.
Produce a Bill of Materials that includes the following entries: Quantity, Description, Source,
Part number, Price. If something is not available from a vendor, then mark it as “in-house.” The
part numbers should correspond to the labels in your original hand drawings. For example
Qty Description Source Part Number Price Total
8 1”x1-1/2” Shoulder Bolts McMaster-Carr 91259A132 $16.82 $134.56
2 Rear Guide Links In-house r2 TBD
.. .. .. .. .. ..
. . . . . .
Grand Total
Produce full page engineering drawings to manufacture each part. Include all pertinent dimen-
sions. Drawings are not required for items available from a vendor, only items in the bill of
materials whose source is in-house need to be drawn.
Create a SolidWorks assembly to demonstrate the functionality of the mechanism. No parts
should bind or interfere. Apply a linear motor to the hydraulic cylinder and produce an animation
of the assembly.
Technical drawings of each part, an image of the assembly in the low position, an image of the
assembly in the up position, and a bill of materials are due on the last day of classes.
A.4. SolidWorks Parts and Assembly 148
Figure A.2.: Shoulder Screw description from http://www.mcmaster.com. There are also metric
equivalent sizes that can be specified.
A.4. SolidWorks Parts and Assembly 149
Figure A.3.: Rod end description from http://www.mcmaster.com. These can be found in both
right hand thread and left hand thread configurations. A jam nut should always be
used to secure the rod end.
A.4. SolidWorks Parts and Assembly 150
Figure A.4.: Rod End drawing to interface with the hydraulic cylinder.
A.4. SolidWorks Parts and Assembly 151
Figure A.5.: External retaining ring description and drawing from www.mcmaster.com. These
are also known as snap rings.
A.4. SolidWorks Parts and Assembly 152