◀ Back to Note

Tony's Intertwingled Memex • Reference Manual

Aerospace Fuel Tank Center of Mass and Slosh Baffle (Aerospace & Mechanical Engineering) Solution

Formulas & Definitions

An aerospace structural dynamics engineer at the UNLV Howard R. Hughes College of Engineering is designing a lightweight titanium fuel cell integrated into the wingtip structure of an autonomous reconnaissance aircraft.

The interior reservoir of the fuel cell is bounded by the three coordinate planes x=0x = 0x=0, y=0y = 0y=0, z=0z = 0z=0 and the upper inclined wing-skin aerodynamic fairing plane:

x+y+z=1x + y + z = 1x+y+z=1
with mass density ρ(x,y,z)=120z kg/m3\rho(x, y, z) = 120 z\text{ kg/m}^3ρ(x,y,z)=120z kg/m3.

Step-by-Step Solution:

  1. Extract Type I Regional Bounds:
    State the Type I Spatial Region Definition (defined in Triple Integrals over a General Region):

    E={(x,y,z)∣0≤x≤1, 0≤y≤1−x, 0≤z≤1−x−y}E = \{(x, y, z) \mid 0 \le x \le 1, \, 0 \le y \le 1 - x, \, 0 \le z \le 1 - x - y\}E={(x,y,z)∣0≤x≤1,0≤y≤1−x,0≤z≤1−x−y}

  2. Evaluate Enclosed Volume:
    State the Volume via Triple Integration Formula:

    V(E)=∭E1 dV=∫01∫01−x∫01−x−y1 dz dy dxV(E) = \iiint_{E} 1 \, dV = \int_{0}^{1} \int_{0}^{1-x} \int_{0}^{1-x-y} 1 \, dz \, dy \, dxV(E)=∭E​1dV=∫01​∫01−x​∫01−x−y​1dzdydx

    Inner zzz-integral:

    ∫01−x−y1 dz=1−x−y\int_{0}^{1-x-y} 1 \, dz = 1 - x - y∫01−x−y​1dz=1−x−y

    Middle yyy-integral:

    ∫01−x(1−x−y) dy=[(1−x)y−y22]01−x=(1−x)2−(1−x)22=12(1−x)2\int_{0}^{1-x} (1 - x - y) \, dy = \left[ (1 - x)y - \frac{y^2}{2} \right]_{0}^{1-x} = (1 - x)^2 - \frac{(1 - x)^2}{2} = \frac{1}{2}(1 - x)^2∫01−x​(1−x−y)dy=[(1−x)y−2y2​]01−x​=(1−x)2−2(1−x)2​=21​(1−x)2

    Outer xxx-integral:

    V(E)=∫0112(1−x)2 dx=[−16(1−x)3]01=0−(−16(1)3)=16 m3V(E) = \int_{0}^{1} \frac{1}{2}(1 - x)^2 \, dx = \left[ -\frac{1}{6}(1 - x)^3 \right]_{0}^{1} = 0 - \left(-\frac{1}{6}(1)^3\right) = \frac{1}{6}\text{ m}^3V(E)=∫01​21​(1−x)2dx=[−61​(1−x)3]01​=0−(−61​(1)3)=61​ m3

  3. Evaluate Total Mass:
    State the Total Mass Formula with Variable Density:

    M=∭Eρ(x,y,z) dV=∫01∫01−x∫01−x−y120z dz dy dxM = \iiint_{E} \rho(x, y, z) \, dV = \int_{0}^{1} \int_{0}^{1-x} \int_{0}^{1-x-y} 120 z \, dz \, dy \, dxM=∭E​ρ(x,y,z)dV=∫01​∫01−x​∫01−x−y​120zdzdydx

    Inner zzz-integral:

    ∫01−x−y120z dz=[60z2]01−x−y=60(1−x−y)2\int_{0}^{1-x-y} 120 z \, dz = \left[ 60 z^2 \right]_{0}^{1-x-y} = 60(1 - x - y)^2∫01−x−y​120zdz=[60z2]01−x−y​=60(1−x−y)2

    Middle yyy-integral:

    ∫01−x60(1−x−y)2 dy=[−20(1−x−y)3]01−x=0−(−20(1−x)3)=20(1−x)3\int_{0}^{1-x} 60(1 - x - y)^2 \, dy = \left[ -20(1 - x - y)^3 \right]_{0}^{1-x} = 0 - \left(-20(1 - x)^3\right) = 20(1 - x)^3∫01−x​60(1−x−y)2dy=[−20(1−x−y)3]01−x​=0−(−20(1−x)3)=20(1−x)3

    Outer xxx-integral:

    M=∫0120(1−x)3 dx=[−5(1−x)4]01=0−(−5(1)4)=5 kgM = \int_{0}^{1} 20(1 - x)^3 \, dx = \left[ -5(1 - x)^4 \right]_{0}^{1} = 0 - (-5(1)^4) = 5\text{ kg}M=∫01​20(1−x)3dx=[−5(1−x)4]01​=0−(−5(1)4)=5 kg

  4. Compute First Moment and Center of Mass Coordinate zˉ\bar{z}zˉ:
    State the First Moment of Mass Formula about the xyxyxy-Plane:

    Mxy=∭Ezρ(x,y,z) dV=∫01∫01−x∫01−x−y120z2 dz dy dxM_{xy} = \iiint_{E} z \rho(x, y, z) \, dV = \int_{0}^{1} \int_{0}^{1-x} \int_{0}^{1-x-y} 120 z^2 \, dz \, dy \, dxMxy​=∭E​zρ(x,y,z)dV=∫01​∫01−x​∫01−x−y​120z2dzdydx

    Inner zzz-integral:

    ∫01−x−y120z2 dz=[40z3]01−x−y=40(1−x−y)3\int_{0}^{1-x-y} 120 z^2 \, dz = \left[ 40 z^3 \right]_{0}^{1-x-y} = 40(1 - x - y)^3∫01−x−y​120z2dz=[40z3]01−x−y​=40(1−x−y)3

    Middle yyy-integral:

    ∫01−x40(1−x−y)3 dy=[−10(1−x−y)4]01−x=10(1−x)4\int_{0}^{1-x} 40(1 - x - y)^3 \, dy = \left[ -10(1 - x - y)^4 \right]_{0}^{1-x} = 10(1 - x)^4∫01−x​40(1−x−y)3dy=[−10(1−x−y)4]01−x​=10(1−x)4

    Outer xxx-integral:

    Mxy=∫0110(1−x)4 dx=[−2(1−x)5]01=2 kg⋅mM_{xy} = \int_{0}^{1} 10(1 - x)^4 \, dx = \left[ -2(1 - x)^5 \right]_{0}^{1} = 2\text{ kg}\cdot\text{m}Mxy​=∫01​10(1−x)4dx=[−2(1−x)5]01​=2 kg⋅m

    State the Center of Mass Coordinate Formula:

    zˉ=MxyM=2 kg⋅m5 kg=0.4 m\bar{z} = \frac{M_{xy}}{M} = \frac{2\text{ kg}\cdot\text{m}}{5\text{ kg}} = 0.4\text{ m}zˉ=MMxy​​=5 kg2 kg⋅m​=0.4 m

Meaning of the Answer & Real-Life Application:

  • Quantitative Interpretation: The cell encloses a total volume of 16 m3≈0.167 m3\frac{1}{6}\text{ m}^3 \approx 0.167\text{ m}^361​ m3≈0.167 m3 containing 5 kg5\text{ kg}5 kg of fuel. Because the density ρ(z)=120z\rho(z) = 120zρ(z)=120z increases toward the upper reservoir ceiling, the vertical center of mass zˉ=0.4 m\bar{z} = 0.4\text{ m}zˉ=0.4 m is located higher than that of a uniform-density tetrahedron (which would have its centroid at zˉuniform=14=0.25 m\bar{z}_{\text{uniform}} = \frac{1}{4} = 0.25\text{ m}zˉuniform​=41​=0.25 m).
  • Engineering Application: In unmanned combat air vehicle (UAV) design, wingtip mass distribution strongly affects the wing's torsional natural frequency and roll moment of inertia IxxI_{xx}Ixx​. Locating the center of gravity at zˉ=0.4 m\bar{z} = 0.4\text{ m}zˉ=0.4 m rather than 0.25 m0.25\text{ m}0.25 m alters the vertical moment arm to the wing chord neutral axis. Autopilot flight control laws must incorporate this vertical offset into aileron feedforward gains to prevent roll-yaw dynamic coupling during turbulent wind gusts.