Force Analysis - 4 Bar - Bill Pollard

Problem Statement

Given position, velocity, acceleration, and inertia properties such as mass and mass moment of inertia for each moving link of a four-bar linkage, we are able to perform force analysis for the linkage. Three static equilibrium equations, in terms of forces in the X and Y directions and moment about the center of gravity of the link, can be written for each link. The free body diagram is also available.

Schematic

flickr:5265940905

Assumptions and Parameters

Unit Type: SI US Customary

Link lengths (m):
r1: 0.50
r2: 0.15
r3: 0.5
r4: 0.1778

Link masses (kg):
m2: 0.36
m3: 1.0
m4: 0.64

Link moments of inertia (kg-m^2 or lbf-ft-sec^2):
ig2:0.00136
ig3:0.01345
ig4:0.00429

Position Vector for Center of Gravity

Magnitudes (m):
rg2:0.0508
rg3:0.1524
rg4:0.0889

Mode for all angles: Degrees
theta1: 0.0

Select and input the known angle (theta2, theta3, or theta4):
theta2: 150

Departures from Link Positions:
delta2: 0.0
delta3: 0.0
delta4: 0.0

External load (N-m or ft-lbf) on joint (1,4):
tl:5.0

Select and input the known angular velocity:
omega2:5 rad/s

Select and input the known angular acceleration:
alpha2: -5 rad/s

Force Analysis

For link 2, we get

F12x + F32x + Fg2x = 0 (1)

-m2g + F12y + F32y + Fg2y = 0 (2)

Ts + (-rg2) x F12 + (r2-rg2) x F32 + Tg2 = 0 (3)

where rg2=rg2*exp(i(theta2+delta2)) is the position vector from joint Ao to the center of gravity of link 2. F12 and F32 are the joint forces acting on link 2. Fg2 and Tg2 are the inertia force and inertia moment of link 2. m2 is the mass of link 2 and Ts is the driving torque. For link 3 we get

F23x + F43x + Fg3x = 0 (4)

-m3g + F23y + F43y + Fg3y = 0 (5)

(-rg3) x F23 + (r3-rg3) x F43 + Tg3 = 0 (6)

where rg3 = rg3*exp(i(theta3+delta3)) is the position vector from joint A to the center of gravity of link 3. F23 and F43 are the joint forces acting on link 3. Fg3 and Tg3 are the inertia force and inertia moment of link 3. m3 is the mass of link 3.

For link 4 we get

F34x + F14x + Fg4x = 0 (7)

-m4g + F34y + F14y + Fg4y = 0 (8)

(-rg4) x F14 + (r4-rg4) x F34 + Tg4 + Tl = 0 (9)

where rg4=rg4*exp(i(theta4+delta4)) is the position vector from joint Bo to the center of gravity of link 4. F14 and F34 are the joint forces acting on link 4. Fg4 and Tg4 are the inertia force and inertia moment of link 4. m4 is the mass of link 4 and Tl is the torque of external load.

Equations (3), (6), and (9) can be expressed as

Ts - rg2*cos(theta2+delta2)*F12y + rg2*sin(theta2+delta2)*F12x

[r2*cos(theta2) - rg2*cos(theta2+delta2)]*F32y - [r2*sin(theta2)

- rg2*cos(theta2+delta2)]*F32x + Tg2 = 0 (10)

-rg3*cos(theta3+delta3)*F23y + rg3*sin(theta3+delta3)*F23x

[r3*cos(theta3) - rg3*cos(theta3+delta3)]*F43y - [r3*sin(theta3)

- rg3*cos(theta3+delta3)]*F43x + Tg3 = 0 (11)

-rg4*cos(theta4+delta4)*F14y + rg4*sin(theta4+delta4)*F14x

[r4*cos(theta4) - rg4*cos(theta4+delta4)]*F34y - [r4*sin(theta4)

- rg4*cos(theta4+delta4)]*F34x + Tg4 + Tl = 0 (12)

Note that Fijx = -Fjix and Fijy = -Fjiy, equations (1-9) can be written as nine linear equations in terms of nine unknowns. They can be expressed in a symbolic form

Ax = b (13)

where x = the transpose of (F12x, F12y, F23x, F23y, F34x, F34y, F14x, F14y, Ts) and is a vector consisting of the unknown forces and input torque, b = the transpose of (Fg2x, Fg2y-m2g, Tg2, Fg3x, Fg3y-m3g, Tg3, Fg4x, Fg4y-m4g, Tg4 + Tl) and is a vector that contains extenal load plus inertia forces and inertia torques.

The interface below allows the user to find the joint forces and the required input torque given the link lengths, theta1, one additional position, one angular velocity, one angular acceleration, the external load and the inertia properties for the mechanism.

Fourbar Dynamics Analysis Results:

r1 = 0.500, r2 = 0.150, 3 = 0.500, r4 = 0.178;
m2 = 0.360, m3 = 1.000, m4 = 0.640;
ig2 = 0.001, ig3 = 0.013, ig4 = 0.004;
rg2 = 0.051, rg3 = 0.152, rg4 = 0.089;
delta2 = 0.000 radians (0.00 degrees),
delta3 = 0.000 radians (0.00 degrees),
delta4 = 0.000 radians (0.00 degrees);
tl = 5.000 N-m (6.779 lb-ft);
theta1 = 0.000 radians (0.00 degrees);
theta2 = 2.618 radians (150.00 degrees),
omega2 = 5.000 rad/sec (286.48 deg/sec),
alpha2 = -5.000 rad/sec^2 (-286.48 deg/sec^2)

Circuit 1: (Joint Reaction Forces and Input Torques:

F12x = 49.296 N,
F12y = 13.794 N,
F23x = 48.855 N,
F23y = 10.412 N,
F34x = 45.315 N,
F34y = 1.207 N,
F14x = -44.244
N, F14y = 5.331 N,
Ts = -5.183 N*m

Circuit 2: (Joint Reaction Forces and Input Torques:

F12x = -38.374 N,
F12y = 23.004 N,
F23x = -38.816 N,
F23y = 19.622 N,
F34x = -41.821 N,
F34y = 11.415 N,
F14x = 42.333 N,
F14y = -5.925 N,
Ts = 0.195 N*m