β—€ Back to Note

Tony's Intertwingled Memex β€’ Reference Manual

Hydraulic Excavator Bracket Stress Analysis (Mechanical Engineering) Solution

Formulas & Definitions

A mechanical design engineer analyzes the localized contact interface on a cast-steel hydraulic pivot bracket for a heavy-duty mining excavator. The curved load-bearing shoulder of the bracket is modeled mathematically by the logarithmic surface:

z=ln⁑(x2+y)z = \ln(x^2 + y)z=ln(x2+y)
where x,y,zx, y, zx,y,z represent spatial coordinates in inches (in\text{in}in).

During heavy bucket digging, a hydraulic actuator piston pushes a spherical pin directly against the bracket shoulder at the localized contact coordinate (x0,y0)=(1,eβˆ’1)Β in(x_0, y_0) = (1, e - 1)\text{ in}(x0​,y0​)=(1,eβˆ’1)Β in.

  1. Compute the baseline height z0=f(1,eβˆ’1)z_0 = f(1, e - 1)z0​=f(1,eβˆ’1) on the bracket shoulder.
  2. Apply the multivariable chain rule to compute the partial derivatives fx(x,y)f_x(x, y)fx​(x,y) and fy(x,y)f_y(x, y)fy​(x,y).
  3. Evaluate the partial derivatives at (1,eβˆ’1)(1, e - 1)(1,eβˆ’1) and formulate the equation of the tangent plane to the bracket contact zone in both explicit (z=…z = \dotsz=…) and general (ax+by+cz=dax + by + cz = dax+by+cz=d) forms.
  4. If the hydraulic piston applies a normal compressive force of magnitude FN=50Β kipsF_N = 50\text{ kips}FN​=50Β kips, explain how finite element analysis (FEA) software uses the tangent plane's normal vector to resolve the force into axial and shear stresses along machine frame axes.

Step-by-Step Solution:

  1. Compute Baseline Coordinate Height:

    z0=f(1,eβˆ’1)=ln⁑(12+eβˆ’1)=ln⁑(1+eβˆ’1)=ln⁑(e)=1Β inz_0 = f(1, e - 1) = \ln\left(1^2 + e - 1\right) = \ln(1 + e - 1) = \ln(e) = 1\text{ in}z0​=f(1,eβˆ’1)=ln(12+eβˆ’1)=ln(1+eβˆ’1)=ln(e)=1Β in
    The pin contact point is P0(1,eβˆ’1,1)P_0(1, e - 1, 1)P0​(1,eβˆ’1,1).

  2. Compute General Partial Derivatives:

    • Let u(x,y)=x2+yu(x, y) = x^2 + yu(x,y)=x2+y. State the Logarithmic Chain Rule:
      βˆ‚βˆ‚x[ln⁑(u)]=1uβˆ‚uβˆ‚x,βˆ‚βˆ‚y[ln⁑(u)]=1uβˆ‚uβˆ‚y\frac{\partial}{\partial x}[\ln(u)] = \frac{1}{u} \frac{\partial u}{\partial x}, \quad \frac{\partial}{\partial y}[\ln(u)] = \frac{1}{u} \frac{\partial u}{\partial y}βˆ‚xβˆ‚β€‹[ln(u)]=u1β€‹βˆ‚xβˆ‚u​,βˆ‚yβˆ‚β€‹[ln(u)]=u1β€‹βˆ‚yβˆ‚u​
    • Partial derivatives:
      fx(x,y)=1x2+yβ‹…(2x)=2xx2+yf_x(x, y) = \frac{1}{x^2 + y} \cdot (2x) = \frac{2x}{x^2 + y}fx​(x,y)=x2+y1​⋅(2x)=x2+y2x​
      fy(x,y)=1x2+yβ‹…(1)=1x2+yf_y(x, y) = \frac{1}{x^2 + y} \cdot (1) = \frac{1}{x^2 + y}fy​(x,y)=x2+y1​⋅(1)=x2+y1​
  3. Evaluate at (1,eβˆ’1)(1, e - 1)(1,eβˆ’1) and Formulate Tangent Plane:

    • Note that x2+y=12+(eβˆ’1)=ex^2 + y = 1^2 + (e - 1) = ex2+y=12+(eβˆ’1)=e.
    • Partial derivative evaluations:
      fx(1,eβˆ’1)=2(1)e=2ef_x(1, e - 1) = \frac{2(1)}{e} = \frac{2}{e}fx​(1,eβˆ’1)=e2(1)​=e2​
      fy(1,eβˆ’1)=1ef_y(1, e - 1) = \frac{1}{e}fy​(1,eβˆ’1)=e1​
    • State the Tangent Plane Equation (defined in Tangent Planes):
      z=z0+fx(x0,y0)(xβˆ’x0)+fy(x0,y0)(yβˆ’y0)z = z_0 + f_x(x_0, y_0)(x - x_0) + f_y(x_0, y_0)(y - y_0)z=z0​+fx​(x0​,y0​)(xβˆ’x0​)+fy​(x0​,y0​)(yβˆ’y0​)
    • Substitute parameters:
      z=1+2e(xβˆ’1)+1e(yβˆ’(eβˆ’1))z = 1 + \frac{2}{e}(x - 1) + \frac{1}{e}(y - (e - 1))z=1+e2​(xβˆ’1)+e1​(yβˆ’(eβˆ’1))
      z=1+2exβˆ’2e+1eyβˆ’1+1ez = 1 + \frac{2}{e}x - \frac{2}{e} + \frac{1}{e}y - 1 + \frac{1}{e}z=1+e2​xβˆ’e2​+e1​yβˆ’1+e1​
      z=2ex+1eyβˆ’1ez = \frac{2}{e}x + \frac{1}{e}y - \frac{1}{e}z=e2​x+e1​yβˆ’e1​
    • General Form (clearing the denominator eee):
      2x+yβˆ’ez=12x + y - ez = 12x+yβˆ’ez=1
  4. FEA Contact Mechanics Application:

    • The surface normal vector at the contact point is nβƒ—=⟨2,1,βˆ’e⟩\vec{n} = \langle 2, 1, -e \ranglen=⟨2,1,βˆ’e⟩.
    • The unit normal vector is:
      n^=⟨2,1,βˆ’e⟩22+12+(βˆ’e)2=⟨2,1,βˆ’e⟩5+e2\hat{n} = \frac{\langle 2, 1, -e \rangle}{\sqrt{2^2 + 1^2 + (-e)^2}} = \frac{\langle 2, 1, -e \rangle}{\sqrt{5 + e^2}}n^=22+12+(βˆ’e)2β€‹βŸ¨2,1,βˆ’eβŸ©β€‹=5+e2β€‹βŸ¨2,1,βˆ’eβŸ©β€‹
    • In FEA boundary conditions, the 50Β kip50\text{ kip}50Β kip force is applied along Fβƒ—=βˆ’50n^Β kips\vec{F} = -50 \hat{n}\text{ kips}F=βˆ’50n^Β kips. FEA solvers project this force vector onto the local element coordinate axes to calculate von Mises yield stresses and verify that peak contact stresses do not exceed the steel's yield strength (Οƒy=70Β ksi\sigma_y = 70\text{ ksi}Οƒy​=70Β ksi).