Beam-in-Socket Analysis

×

Beam-in-Socket

Beam-in-Socket Analysis (Pin-in-Socket Analysis)

1. Introduction and Physical Basis

The beam-in-socket analysis, also called pin-in-socket analysis, determines the stress and load distribution in a cylindrical pin or beam that is engaged in a socket (hole) over a finite length. This analysis is critical for:

• Pin joints in aircraft structures (lugs, clevises, shear pins)
• Shoulder bolts with no clamping force
• Continuous pistons in cylinders
• Dowel pins and alignment pins

Key Assumptions:

• Continuous contact: The pin maintains contact with the socket over the entire engagement length L
• Tight fit: No gap exists between pin and socket (interference or line-to-line fit)
• Rigid socket: The socket walls are assumed rigid compared to the pin
• Linear bearing distribution: The bearing pressure varies linearly along the engagement length
• Elastic behavior: Both pin and socket remain in the elastic range

2. Problem Setup and Coordinate System

Geometry:

• L = engagement length (depth of socket)
• D = pin diameter
• x = distance measured from the faying surface (loaded end), where 0 ≤ x ≤ L

Applied Loads (at x = 0, the faying surface):

• V_o = applied shear force (transverse to pin axis)
• M_o = applied moment (bending moment about pin axis)

Sign Convention:

• x-axis points INTO the socket (from faying surface toward pin end)
• Positive shear V acts in +y direction
• Positive moment M causes compression on +y side
• Bearing pressure ω(x) acts radially inward on the pin (reaction from socket)

3. Bearing Stress Distribution Model

The bearing pressure (force per unit length) ω(x) is assumed to vary linearly along the engagement length:

ω(x) = ω₁ - (ω₁ - ω₂) · (x / L)

where:

• ω₁ = bearing load per unit length at x = 0 (faying surface)
• ω₂ = bearing load per unit length at x = L (pin end)

This can be rewritten as:

ω(x) = ω₁ · (1 - x/L) + ω₂ · (x/L)

The bearing pressure is maximum at the loaded end (x=0) and minimum at the far end (x=L).

4. Equilibrium Equations

4.1 Force Equilibrium (Vertical Direction)

The total bearing force must balance the applied shear force. The bearing forces oppose the applied shear, giving:

∫₀^L ω(x) dx = -V_o

Note: The negative sign indicates that the bearing forces oppose the applied shear force V_o.

Substituting the linear distribution:

∫₀^L [ω₁ - (ω₁ - ω₂) · (x/L)] dx = -V_o

Evaluating the integral:

ω₁ · L - (ω₁ - ω₂) · (L/2) = -V_o

Simplifying:

ω₁ · L - (ω₁/2) · L + (ω₂/2) · L = -V_o

(L/2) · (ω₁ + ω₂) = -V_o

Therefore:

Equation 1: ω₁ + ω₂ = -2 · V_o / L

4.2 Moment Equilibrium (About x = 0)

The moment of the bearing forces about the faying surface must balance the applied moment. With x pointing into the socket and M_o defined as positive clockwise, the bearing forces create counterclockwise moments, giving:

∫₀^L ω(x) · x · dx = -M_o

Note: The negative sign accounts for the sign convention where positive M_o (clockwise) is opposed by the bearing load distribution.

Substituting:

∫₀^L [ω₁ - (ω₁ - ω₂) · (x/L)] · x · dx = -M_o

Evaluating:

ω₁ · (L²/2) - (ω₁ - ω₂) · (L²/3) = -M_o

Simplifying:

ω₁ · (L²/2) - ω₁ · (L²/3) + ω₂ · (L²/3) = -M_o

ω₁ · (L²/6) + ω₂ · (L²/3) = -M_o

Multiply by 6/L²:

Equation 2: ω₁ + 2 · ω₂ = -6 · M_o / L²

5. Solution for ω₁ and ω₂

We now have two equations with two unknowns:

Equation 1: ω₁ + ω₂ = -2 · V_o / L

Equation 2: ω₁ + 2 · ω₂ = -6 · M_o / L²

5.1 Solve for ω₂

Subtract Equation 1 from Equation 2:

(ω₁ + 2ω₂) - (ω₁ + ω₂) = (-6M_o/L²) - (-2V_o/L)

ω₂ = -6M_o/L² + 2V_o/L

Rearranging:

ω₂ = 2V_o/L - 6M_o/L²

5.2 Solve for ω₁

From Equation 1:

ω₁ = -2V_o/L - ω₂

Substitute ω₂:

ω₁ = -2V_o/L - (2V_o/L - 6M_o/L²)

ω₁ = -2V_o/L - 2V_o/L + 6M_o/L²

ω₁ = -4V_o/L + 6M_o/L²

Or equivalently:

ω₁ = -(4V_o/L + 6M_o/L²) when M_o < 0 (counterclockwise)

5.3 Final Expressions for Bearing Intensities

Peak bearing at faying surface (x = 0):

ω₁ = -4 · V_o / L - 6 · M_o / L²

Or written as:

ω₁ = -(4V_o/L + 6M_o/L²)

Peak bearing at pin end (x = L):

ω₂ = 2 · V_o / L - 6 · M_o / L²

Note on Signs: For a tight-fit pin joint with V_o positive (upward applied load) and M_o = 0:

• ω₁ < 0 indicates bearing force in the negative direction (downward) at the faying surface, opposing the upward applied shear V_o
• ω₂ > 0 indicates bearing force in the positive direction (upward) at the far end
• This creates an S-shaped deflection: pin pushed down at loaded end, pushed up at far end
• With a tight fit, the pin maintains continuous contact throughout the engagement length on alternating sides
• The magnitude |ω₁| represents the intensity of the bearing force FROM the socket ONTO the pin

6. Peak Bearing Stress

The bearing stress is the bearing force per unit length divided by the pin diameter. The peak bearing stress occurs at whichever end has the higher magnitude:

f_br = max(|ω₁|, |ω₂|) / D

where D is the pin diameter. The absolute value is used since we are interested in the magnitude of the bearing stress, regardless of direction.

Note: While ω₁ is often larger in magnitude for typical loading (making the faying surface the most critical location), the peak can occur at either end depending on the applied moment M_o.

7. Internal Pin Shear Distribution

The internal shear force V(x) in the pin varies along its length due to the distributed bearing load. From equilibrium of a pin segment from 0 to x:

V(x) = V_o + ∫₀^x ω(ξ) dξ

Note: The bearing load ω acts as a distributed load on the pin, so it integrates with a positive sign to build up internal shear.

Substituting the linear bearing distribution ω(x) = ω₁ + (ω₂ - ω₁)·(x/L):

V(x) = V_o + ∫₀^x [ω₁ + (ω₂ - ω₁) · (ξ/L)] dξ

Integrating:

V(x) = V_o + ω₁ · x + (ω₂ - ω₁) · x² / (2L)

7.1 Location of Critical Shear

To find where the shear reaches an extremum (minimum or maximum), take dV/dx = 0:

dV/dx = ω₁ + (ω₂ - ω₁) · x / L = 0

Solving for x:

x_V,crit = -ω₁ · L / (ω₂ - ω₁)

Substituting this critical location back into V(x):

V_interior = V_o + (L/2) · ω₁² / (ω₁ - ω₂)

Note: The sign correction from earlier versions changes the sign of the second term. This interior shear value may be larger or smaller than the boundary values depending on the loading.

7.2 Maximum Shear

The maximum shear must be determined by comparing the interior critical value with the boundary values:

V_max = max(|V_o|, |V_interior|)

where:

• V_o is the applied shear at x = 0 (faying surface)
• V_interior is the shear at the interior critical location x_V,crit

Important: The maximum shear is not always at the faying surface - it can occur at an interior location depending on the bearing distribution.

8. Internal Pin Bending Moment Distribution

The internal bending moment M(x) is found by integrating the shear force:

M(x) = M_o + ∫₀^x V(ξ) dξ

Substituting V(x):

M(x) = M_o + ∫₀^x [V_o + ω₁ · ξ + (ω₂ - ω₁) · ξ² / (2L)] dξ

Integrating:

M(x) = M_o + V_o · x + ω₁ · x² / 2 + (ω₂ - ω₁) · x³ / (6L)

Note: The signs are corrected to match the equilibrium conventions established earlier. The bearing forces create moments that add to M_o.

8.1 Locations of Extreme Moments

Extreme moments occur where dM/dx = V(x) = 0. Setting V(x) = 0 and solving the quadratic equation:

V(x) = V_o + ω₁ · x + (ω₂ - ω₁) · x² / (2L) = 0

Note: We rearrange to use (ω₂ - ω₁) for cleaner algebra.

Using the quadratic formula with a = (ω₂ - ω₁)/(2L), b = ω₁, c = V_o:

x = [-ω₁ ± (ω₁² - 2 · V_o · (ω₂ - ω₁) / L)¹⁄²] / [(ω₂ - ω₁) / L]

This gives two potential critical locations:

x₁ = [-ω₁ - (ω₁² - 2 · V_o · (ω₂ - ω₁) / L)¹⁄²] / [(ω₂ - ω₁) / L]

x₂ = [-ω₁ + (ω₁² - 2 · V_o · (ω₂ - ω₁) / L)¹⁄²] / [(ω₂ - ω₁) / L]

Important Notes:

• The quadratic formula gives two roots, representing potential local maximum and minimum
• Both roots may be in the valid range [0, L], or only one, or neither
• The discriminant must be non-negative: ω₁² - 2·V_o·(ω₂ - ω₁)/L ≥ 0
• If the discriminant is negative, no interior extrema exist and the maximum moment occurs at a boundary

8.2 Maximum Moment Value

The maximum absolute moment must be found by evaluating M(x) at all potential locations:

M(x_i) = M_o + V_o · x_i + ω₁ · x_i² / 2 + (ω₂ - ω₁) · x_i³ / (6L)

Evaluate this at:

• x = 0 (faying surface): M(0) = M_o
• x = x₁ (if 0 ≤ x₁ ≤ L): M(x₁)
• x = x₂ (if 0 ≤ x₂ ≤ L): M(x₂)

Then:

M_max = max(|M_o|, |M(x₁)|, |M(x₂)|)

Note: We use absolute values because we're interested in the maximum bending stress magnitude, regardless of whether it's tension or compression.

9. Summary of Key Equations

Bearing Load Intensities:

ω₁ = -4 · V_o / L - 6 · M_o / L²     [at faying surface, x = 0]

ω₂ = 2 · V_o / L + 6 · M_o / L²     [at pin end, x = L]

Peak Bearing Stress:

f_br = max(|ω₁|, |ω₂|) / D

Internal Shear Distribution:

V(x) = V_o + ω₁ · x + (ω₂ - ω₁) · x² / (2L)

x_V,crit = -ω₁ · L / (ω₂ - ω₁)     [location of extreme shear]

V_interior = V_o + (L/2) · ω₁² / (ω₁ - ω₂)     [interior extreme shear]

V_max = max(|V_o|, |V_interior|)     [maximum shear magnitude]

Internal Moment Distribution:

M(x) = M_o + V_o · x + ω₁ · x² / 2 + (ω₂ - ω₁) · x³ / (6L)

x₁ = [-ω₁ - (ω₁² - 2·V_o·(ω₂ - ω₁)/L)¹⁄²] / [(ω₂ - ω₁)/L]     [first critical location]

x₂ = [-ω₁ + (ω₁² - 2·V_o·(ω₂ - ω₁)/L)¹⁄²] / [(ω₂ - ω₁)/L]     [second critical location]

M_max = max(|M_o|, |M(x₁)|, |M(x₂)|)     [maximum moment magnitude]

10. Special Case: Shear Only (M_o = 0)

For a symmetrically loaded pin joint with no applied moment and positive upward shear V_o:

ω₁ = -4 · V_o / L     (bearing downward at loaded end, opposing V_o)

ω₂ = 2 · V_o / L     (bearing upward at far end)

Note: The negative ω₁ indicates the bearing force opposes the applied shear. With a tight fit, the pin contacts the bottom of the socket near the faying surface (pushed down) and the top of the socket near the far end (pushed up), creating an S-shaped deflection curve.

11. Limitations and Practical Considerations

• Tight fit assumption: This analysis assumes continuous contact between pin and socket. Negative ω₂ indicates the pin contacts the opposite side of the socket bore at x = L, forming an S-curve deflection pattern.
• Elastic limit: The analysis assumes elastic behavior. Yielding at high bearing stresses invalidates the linear distribution.
• Clearance fits: Gaps between pin and socket (loose fits) require different analysis methods.
• Bending stiffness: Very flexible pins may have different bearing distributions (requires beam-on-elastic-foundation analysis).
• Edge effects: Stress concentrations at the socket edges are not captured by this simplified model.

12. Applications

This analysis is used for:

• Aircraft lug and pin joints (MIL-HDBK-5, MMPDS)
• Clevis and tang joints in landing gear
• Control surface hinges
• Shoulder bolts in mechanical assemblies
• Piston-cylinder contacts in hydraulic systems

Close