3.6) Find the velocities of all the bodies of a slider crank mechanism, as shown in Example 3.

1, in
terms of  2 , the angular velocity of the crank shaft.

Solution:

There are 3 bodies in this system if the ground is ignored. The coordinates of the bodies are

T
for the crank shaft: q 2   Rx2 Ry2  2 
T
for the connecting rod: q3   Rx3 Ry3  3 
T
for the slider: q 4   Rx4 Ry4  4 

There are two simple ways to find the solution to this problem. One is using the methods of
multi-body dynamics, and the other is through simple geometric reasoning.

Due to the simplicity of this problem, simple geometric reasoning can produce the solution very
rapidly. Note from the diagram of the mechanism that the following relationships can be
deduced.

l 2 sin  2  l 3 sin  3  0 (EQAB)

l 2 cos  2  l 3 cos  3  Rx4 (EQAC)

l2
Rx2  cos  2 (EQAD)
2

l2
Ry2  sin  2 (EQAE)
2

l3
Rx3  l 2 cos  2  cos  3 (EQAF)
2

Ry2  Ry3 (EQAG)

These equations alone can be used to find all velocity relationships of this mechanism.

Differentiating Eq. (EQAB) and solving for  3 , one gets

l 2 cos  2 2
 3    (EQAH)
l 3 cos  3

Differentiating Eq. (EQAC) and substituting Eq. (EQAH), one gets

 l 2 cos  2 sin  3  2
 Rx4   l 2 sin  2    (EQAW)
 cos  3 
From Eqs. (EQAD) and (EQAE), one gets

l2
 Rx2   sin  2 2 (EQAI)
2

l2
 Ry2  cos  2 2 (EQAJ)
2

From Eqs. (EQAF) and (EQAH), one gets

 2 l 2 cos  2 sin  3  2
 R   l sin  
3 2
 
2cos  3
x
 

From Eqs. (EQAG) and (EQAJ), one gets

l2
 R  cos  2 2
3
y
2

Therefore, in matrix form, using geometric relationships leads directly to

 l2 
  sin  2 
 2 
 l2 
 Rx2   cos  2

2
 2  
 Ry   1 
  2   l cos  sin 
2 2 3
 3   l 2 sin  2  
 Rx   2 cos  3 
 R 3    l 2   2 (EQAV)
 3 
y
cos  2

    2 
 R 4   l cos 
2 2 
 x   3 
 Ry4   l cos  3 
 4  2 l 2 cos  2 sin  3 
    l sin  2  
 cos  3 
 0 
 
 0 

On the other hand, the methods of multi-body dynamics can also be used to find the same
solution, though written in a different form. In this method, the Jacobian matrix of the constraint
equations is required. First, the constraint equations are presented.

The constraint equations at Point O are

R 2  A 2uo2  0 (EQAK)

where
 T
cos  2  sin  2   l2 
A 
2
 and u 2
  0
 sin  cos  2 
2 o
 2 

Expanding Eq. (EQAK) produces

 l2 
 Rx2  cos  2  sin  2     0 
 2    2 
 Ry   sin 
2
cos  2   0 
 0 

which is

l2
Rx2  cos  2  0
2 (EQAL)
l2
Ry2  sin  2  0
2

The constraint equations at Point A are

R 2  A 2u 2A  R 3  A 3u 3A  0 (EQAM)

T T
l2   l3 
where u  2
A 0  and u 3A    0
2   2 

which is

l2 l3
Rx2  cos  2  Rx3  cos  3  0
2 2 (EQAN)
2
l l3
Ry2  sin  2  Ry3  sin  3  0
2 2

Expanding Eq. (EQAM) produces

l2   l3 
 Rx2  cos  2  sin  2     Rx  cos  3
3
 sin  3     0 
 2    2   3    2 
 Ry   sin 
2
cos  2     Ry   sin  3 cos  3   0 
 0   0 

The constraint equations for the revolute joint at Point B are

R 3  A3u 3B  R 4  0 (EQAO)

Note that Point B is the center of mass of Body 4. Therefore, u 4B  0

Expanding Eq. (EQAO) produces

 3    2   4   
 Ry   sin 
3
cos  3     Ry  0 
 0 

which is

l3
Rx3  cos  3  Rx4  0
2 (EQAP)
l3
Ry3  sin  3  Ry4  0
2

The equations for the prismatic joint for Body 4 are

In order to build the Jacobian matrix of the constraints, the derivatives of the constraint equations
are required.

The virtual changes in the coordinates for the revolute joint at Point O are

l2
 Rx2  sin  2 2  0
2 (EQAQ)
l2
 Ry  cos  2 2  0
2

which, when written in matrix form, produce

 l2 
sin  2   Rx  0 
2
 1 0
 
 2   Ry2   0 
  
2
2
l 2
 0 1  cos      0 
2

The virtual changes in the coordinates for the revolute joint at Point A are

l2 l3
 Rx2  sin  2 2   Rx3  sin  3 3  0
2 2 (EQAR)
2
l l3
 Ry2  cos  2 2   Ry3  cos  3 3  0
2 2

which, when written in matrix form, produce

 Rx2 
 2
 l2 l3 3 
 Ry 
1 0  2 sin 
2
1 0  sin    2  0 

 2   3
 0 
 l2 3
 R
cos  3   x3  0 
l
 0 1 cos  2 0 1
2 2   Ry 
 3
  
The virtual changes in the coordinates for the revolute joint at Point B are

l3
 Rx3  sin  3 3   Rx4  0
2 (EQAS)
l3
 Ry3  cos  3 3   Ry4  0
2

which, when written in matrix form, produce

 Rx2 
 2
 l3   Ry  0
1 0  2 sin 
3
1 0 0   2   

   3   0
 l3  R  
 0 1 cos  3 0 1 0   x3  0 
2   R y 
 3
  

The virtual changes in the coordinates for the prismatic joint for Body 4 are

 Ry4  0 
 4  (EQAT)
   0 

The Jacobian matrix of the constraints, Cq , can now be created.

 l2 
1 0 sin  2 0 0 0 0 00
 2 
 l2 
0 1  cos  2 0 0 0 0 0 0
2
 
 l2 l3 
1 0  sin  2 1 0  sin  3 0 0 0
2 2
 
Cq   0 0 0 0 
2
l l3
 1 cos  2 0 1 cos  3
2 2
 
0 l3
0 0 1 0  sin  3 1 0 0 
 2 
 3 
0 0 0 0 1
l
cos  3 0 1 0 
 2 
0 0 0 0 0 0 0 1 0 

0 0 0 0 0 0 0 0 1 

Note that Cq is an 8  9 matrix. The Jacobian matrix of the constraints is a product of

Cq q  0
where

 Rx2 
 2
 Ry 
  2 
 3
 Rx 
 q   Ry3 
 3
  
 R 4 
 x
 Ry4 
 4
  

Cq q  0 can be separated into

Cq d  q d  Cq i  q i  0 (EQAU)

where

1 0 0 0 0 0 0
0  l2 
0   sin  2 
 1 0 0 0 0 0 0  22 
    Rx2
  l 2
l3
1 0 1 0  sin  3 0 0 0  2
 Ry    2 cos  
 2    2 
 l3   Rx3   l 2
0 1 0 1 cos  3 0 0 0  3
 Ry    2 sin  
Cq d  2  ,  qd   , Cq i   2 ,
 l3    3   l cos  2 
0 0 1 0  sin  3 1 0 0   4  2 
 2   Rx   
 l 3
  R 4   0 
0 0 0 1 cos  3 0 1 0   y  
2   4  0
   
0 0 0 0 0 0 1 0  0 
0   
 0 0 0 0 0 0 1  0 

and  qi   2 

Cqd is an invertible 8  8 matrix. It is just Cq with the third column removed. Cqi is just the
third column of Cq .

Solving Eq. (EQAU) produces

 q d  Cq1 Cq  qi
d i

which is the same relationship as found in Eq. (EQAV). This can be verified using symbolic
mathematical manipulation software such as Maple.