Level 4 — ApplicationSpacecraft Structures & Systems Engineering

Spacecraft Structures & Systems Engineering

60 minutes60 marksprintable — key stays hidden on paper

Level 4 — Application Examination

Time limit: 60 minutes Total marks: 60 Instructions: Answer all questions. Show all working. Use g=9.81 m/s2g = 9.81\ \text{m/s}^2 where needed. Symbols carry their usual meaning.


Question 1 — Launch load path & buckling (14 marks)

A cylindrical satellite adapter is a thin-walled aluminium tube of outer radius R=0.30 mR = 0.30\ \text{m}, wall thickness t=3.0 mmt = 3.0\ \text{mm}, and length L=1.2 mL = 1.2\ \text{m}. It supports a spacecraft of mass m=850 kgm = 850\ \text{kg} during a launch phase with steady axial acceleration of 6.0g6.0\,g plus a quasi-static bending moment from wind shear of M=4.0 kN⋅mM = 4.0\ \text{kN·m}. Aluminium: E=70 GPaE = 70\ \text{GPa}, yield stress σY=280 MPa\sigma_Y = 280\ \text{MPa}.

(a) Compute the axial compressive load and the axial membrane stress in the wall. (3)

(b) Compute the peak bending stress at the tube surface. (Treat as thin ring: IπR3tI \approx \pi R^3 t.) (3)

(c) Combine (a) and (b) to get the peak compressive stress, and state the margin of safety against yield if the required factor of safety is FOS=1.4\text{FOS} = 1.4. (3)

(d) Estimate the classical shell buckling stress σcr=EtR3(1ν2)\sigma_{cr} = \dfrac{E t}{R\sqrt{3(1-\nu^2)}} with ν=0.33\nu = 0.33, and comment on whether yield or buckling governs. (5)


Question 2 — Random vibration & modal survival (12 marks)

A component is mounted on a bracket whose fundamental mode is at fn=220 Hzf_n = 220\ \text{Hz} with quality factor Q=12Q = 12. The base random vibration input is a flat PSD of W0=0.04 g2/HzW_0 = 0.04\ \text{g}^2/\text{Hz} over the band containing fnf_n.

(a) Using the Miles' equation GRMS=π2fnQW0\text{GRMS} = \sqrt{\dfrac{\pi}{2} f_n\, Q\, W_0}, compute the RMS acceleration response of the component. (4)

(b) State the "3-sigma" peak acceleration used for design and convert it to a force on a component of mass 0.65 kg0.65\ \text{kg}. (4)

(c) The bracket natural frequency requirement is "fn200 Hzf_n \geq 200\ \text{Hz}." Rewrite this as a SMART requirement and state which verification method (analysis/test/inspection/demonstration) is most appropriate and why. (4)


Question 3 — Fatigue and fracture (12 marks)

A structural fitting experiences, per mission, the following cyclic loading blocks. Its S-N behaviour is N=1012σ3N = 10^{12}\,\sigma^{-3} (with σ\sigma in MPa, NN cycles to failure).

Block Stress amplitude σ\sigma (MPa) Cycles per mission nin_i
A 200 3000
B 100 40000

(a) Using Miner's rule, compute the cumulative damage per mission and the number of missions to failure. (6)

(b) The fitting contains a surface crack of length a=1.5 mma = 1.5\ \text{mm}. The stress intensity factor is K=YσπaK = Y\sigma\sqrt{\pi a} with Y=1.12Y = 1.12. Fracture toughness KIC=30 MPamK_{IC} = 30\ \text{MPa}\sqrt{\text{m}}. Find the critical stress for fast fracture at this crack length, and state whether the 200 MPa block is safe against fracture. (6)


Question 4 — Power, thermal & mass budget (12 marks)

A LEO smallsat requires an orbit-average electrical load of 180 W180\ \text{W}. Orbit period 95 min95\ \text{min}, eclipse fraction 35%35\%.

(a) The battery must supply the load through eclipse. Compute the energy (Wh) drawn per eclipse, and the required battery capacity if the allowable depth of discharge is DoD=25%\text{DoD} = 25\% and battery-to-load efficiency is 0.900.90. (4)

(b) Solar array must power the load AND recharge the battery during sunlight. If charging efficiency is 0.850.85, compute the minimum solar array output power needed (assume load runs continuously). (4)

(c) The spacecraft dry mass is estimated at 210 kg210\ \text{kg} with a 20%20\% system mass margin. Propellant is 18 kg18\ \text{kg}. Compute the design dry mass (with margin) and wet mass. (4)


Question 5 — Reliability & redundancy (10 marks)

A subsystem uses two identical units, each with constant failure rate λ=8×106 h1\lambda = 8\times10^{-6}\ \text{h}^{-1}.

(a) Compute the MTTF of a single unit, and its reliability at t=3t = 3 years (use 87668766 h/yr). (4)

(b) The two units are in active (parallel) redundancy (subsystem works if at least one works). Write the subsystem reliability Rsys(t)R_{sys}(t) and evaluate it at t=3t = 3 years. (4)

(c) In one sentence each, distinguish cold standby from hot standby redundancy. (2)

Answer keyMark scheme & solutions

Question 1

(a) Axial load: F=ma=850×6.0×9.81=50,013 N50.0 kNF = m\,a = 850 \times 6.0 \times 9.81 = 50{,}013\ \text{N} \approx 50.0\ \text{kN}. (1) Wall area: A=2πRt=2π(0.30)(0.003)=5.655×103 m2A = 2\pi R t = 2\pi (0.30)(0.003) = 5.655\times10^{-3}\ \text{m}^2. (1) σaxial=F/A=50013/5.655×103=8.84×106 Pa=8.84 MPa\sigma_{axial} = F/A = 50013 / 5.655\times10^{-3} = 8.84\times10^{6}\ \text{Pa} = 8.84\ \text{MPa}. (1)

(b) I=πR3t=π(0.30)3(0.003)=2.545×104 m4I = \pi R^3 t = \pi (0.30)^3 (0.003) = 2.545\times10^{-4}\ \text{m}^4. (1) σbend=MR/I=(4000)(0.30)/2.545×104=4.71×106 Pa=4.71 MPa\sigma_{bend} = M R / I = (4000)(0.30)/2.545\times10^{-4} = 4.71\times10^{6}\ \text{Pa} = 4.71\ \text{MPa}. (2)

(c) Peak compressive stress =8.84+4.71=13.56 MPa= 8.84 + 4.71 = 13.56\ \text{MPa}. (1) Allowable with FOS: σallow=σY/FOS=280/1.4=200 MPa\sigma_{allow} = \sigma_Y/\text{FOS} = 280/1.4 = 200\ \text{MPa}. (1) Margin of Safety =σallowσpeak1=20013.561=13.7= \dfrac{\sigma_{allow}}{\sigma_{peak}} - 1 = \dfrac{200}{13.56} - 1 = 13.7 — very large positive margin (yield not a concern). (1)

(d) σcr=EtR3(1ν2)=70×109×0.0030.303(10.332)\sigma_{cr} = \dfrac{E t}{R\sqrt{3(1-\nu^2)}} = \dfrac{70\times10^9 \times 0.003}{0.30\sqrt{3(1-0.33^2)}}. (2) Denominator: 3(10.1089)=2.673=1.635\sqrt{3(1-0.1089)} = \sqrt{2.673} = 1.635; ×0.30=0.4905\times 0.30 = 0.4905. (1) σcr=2.1×108/0.4905=4.28×108 Pa=428 MPa\sigma_{cr} = 2.1\times10^{8} / 0.4905 = 4.28\times10^{8}\ \text{Pa} = 428\ \text{MPa} (classical). (1) Comparison: applied peak 13.6 MPaσcr13.6\ \text{MPa} \ll \sigma_{cr} classical, and even with a knockdown factor of ~0.3 (real shells) 128 MPa\sim128\ \text{MPa} still far above applied stress. Both yield and buckling have huge margins; the driver is elsewhere (or design is heavily over-conservative). (1)


Question 2

(a) GRMS=π2×220×12×0.04=π2×105.6=165.9=12.88 g\text{GRMS} = \sqrt{\dfrac{\pi}{2} \times 220 \times 12 \times 0.04} = \sqrt{\dfrac{\pi}{2}\times 105.6} = \sqrt{165.9} = 12.88\ \text{g}. (4)

(b) 3-sigma peak =3×12.88=38.6 g= 3 \times 12.88 = 38.6\ \text{g}. (2) Force =m×(3σ)g=0.65×38.6×9.81=246 N= m \times (3\sigma)\, g = 0.65 \times 38.6 \times 9.81 = 246\ \text{N}. (2)

(c) SMART version, e.g.: "The bracket's first natural frequency shall be 200 Hz\geq 200\ \text{Hz}, verified by modal test/analysis at ambient conditions." — Specific (frequency), Measurable (Hz value), Achievable, Relevant (avoids low-frequency coupling with launcher), Testable. (2) Verification: Test (modal/sine-sweep) is most appropriate — natural frequency is directly measurable on hardware; analysis (FEM) may support but a test provides authoritative verification of the as-built structure. (2)


Question 3

(a) Cycles to failure: NA=1012×2003=1012/8×106=1.25×105N_A = 10^{12}\times 200^{-3} = 10^{12}/8\times10^{6} = 1.25\times10^{5}. (1) NB=1012×1003=1012/106=1.0×106N_B = 10^{12}\times 100^{-3} = 10^{12}/10^{6} = 1.0\times10^{6}. (1) Damage per mission: D=nANA+nBNB=3000125000+400001000000=0.024+0.040=0.064D = \dfrac{n_A}{N_A} + \dfrac{n_B}{N_B} = \dfrac{3000}{125000} + \dfrac{40000}{1000000} = 0.024 + 0.040 = 0.064. (3) Missions to failure =1/D=1/0.064=15.615= 1/D = 1/0.064 = 15.6 \Rightarrow 15 complete missions. (1)

(b) Critical stress: σc=KICYπa\sigma_c = \dfrac{K_{IC}}{Y\sqrt{\pi a}}. (2) πa=π×0.0015=4.712×103=0.06865\sqrt{\pi a} = \sqrt{\pi \times 0.0015} = \sqrt{4.712\times10^{-3}} = 0.06865. (1) σc=301.12×0.06865=300.07689=390 MPa\sigma_c = \dfrac{30}{1.12\times 0.06865} = \dfrac{30}{0.07689} = 390\ \text{MPa}. (2) Since 390 MPa>200 MPa390\ \text{MPa} > 200\ \text{MPa} applied, the 200 MPa block is safe against fast fracture (margin ~1.95). (1)


Question 4

(a) Eclipse duration =0.35×95=33.25 min=0.5542 h= 0.35 \times 95 = 33.25\ \text{min} = 0.5542\ \text{h}. (1) Energy per eclipse =180×0.5542=99.75 Wh= 180 \times 0.5542 = 99.75\ \text{Wh}. (1) Accounting for efficiency, energy drawn from battery =99.75/0.90=110.8 Wh= 99.75 / 0.90 = 110.8\ \text{Wh}. (1) Required capacity =110.8/0.25=443.3 Wh= 110.8 / 0.25 = 443.3\ \text{Wh}. (1)

(b) Sunlight duration =0.65×95=61.75 min=1.0292 h= 0.65 \times 95 = 61.75\ \text{min} = 1.0292\ \text{h}. (1) Recharge energy needed (into battery) =110.8 Wh= 110.8\ \text{Wh}; energy the array must generate for charging =110.8/0.85=130.4 Wh= 110.8/0.85 = 130.4\ \text{Wh}. (1) Array power during sunlight == (load) ++ (charge power) =180+130.4/1.0292=180+126.7=306.7 W= 180 + 130.4/1.0292 = 180 + 126.7 = 306.7\ \text{W}. (2)

(c) Design dry mass =210×1.20=252 kg= 210 \times 1.20 = 252\ \text{kg}. (2) Wet mass =252+18=270 kg= 252 + 18 = 270\ \text{kg}. (2)


Question 5

(a) Single-unit MTTF =1/λ=1/(8×106)=1.25×105 h= 1/\lambda = 1/(8\times10^{-6}) = 1.25\times10^{5}\ \text{h}. (2) t=3×8766=26298 ht = 3\times 8766 = 26298\ \text{h}; λt=8×106×26298=0.2104\lambda t = 8\times10^{-6}\times 26298 = 0.2104. R(t)=e0.2104=0.8103R(t) = e^{-0.2104} = 0.8103. (2)

(b) Active parallel: Rsys=1(1R)2=2RR2R_{sys} = 1 - (1-R)^2 = 2R - R^2. (2) Rsys=2(0.8103)(0.8103)2=1.62060.6566=0.9640R_{sys} = 2(0.8103) - (0.8103)^2 = 1.6206 - 0.6566 = 0.9640. (2)

(c) Cold standby: backup unit is unpowered/off until the primary fails, then switched in (lower wear, needs fault detection & switch). Hot standby: backup runs continuously (powered) alongside primary, enabling instant takeover but accumulating wear/aging. (2)


[
  {"claim":"Q1a axial stress ~8.84 MPa","code":"F=850*6.0*9.81; A=2*pi*Rational(3,10)*0.003; s=F/A; result = abs(float(s)-8.84e6) < 5e4"},
  {"claim":"Q1d shell buckling ~428 MPa","code":"scr=(70e9*0.003)/(0.30*sqrt(3*(1-0.33**2))); result = abs(float(scr)-4.28e8) < 5e6"},
  {"claim":"Q2a GRMS ~12.88 g","code":"g=sqrt((pi/2)*220*12*0.04); result = abs(float(g)-12.88) < 0.05"},
  {"claim":"Q3a missions to failure ~15.6","code":"D=3000/125000 + 40000/1000000; result = abs(1/D-15.625) < 0.1"},
  {"claim":"Q3b critical stress ~390 MPa","code":"sc=30/(1.12*sqrt(pi*0.0015)); result = abs(float(sc)-390) < 3"},
  {"claim":"Q5b redundant reliability ~0.964","code":"R=exp(-8e-6*3*8766); Rsys=2*R-R**2; result = abs(float(Rsys)-0.9640) < 0.002"}
]