22-Mec-B9 Advanced Engineering Structures · December 2014
Question 8 of 8: Torsion of a three-cell thin-walled wing box
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
Paper format. National Exams, December 2014 — 07-Mec-B9 Advanced Engineering Structures. Three hours, open book, any non-communicating calculator permitted. Eight questions of equal total value (20 marks each); the rubric states that any five constitute a complete paper. All eight are solved here, because the complete set is the study resource.
Reference texts.
T. H. G. Megson, Aircraft Structures for Engineering Students, 6th ed. — open and closed thin-walled sections, shear flow, shear centre, structural idealisation (Ch. 16, 17, 18, 20, 21, 23).
A. C. Ugural and S. K. Fenster, Advanced Strength and Applied Elasticity, 5th ed. — yield criteria, torsion of thin-walled members, statically indeterminate axial members (Ch. 2, 4, 6).
A. P. Boresi and R. J. Schmidt, Advanced Mechanics of Materials, 6th ed. — torsion of open and closed sections, multiply connected cells (Ch. 6).
N. E. Dowling, Mechanical Behavior of Materials, 4th ed. — linear elastic fracture mechanics, the Paris law, strain-life and Miner’s rule (Ch. 8, 11, 14).
R. C. Hibbeler, Mechanics of Materials, 10th ed. — thermal stress in restrained members and combined loading (Ch. 4, 8).
Axes and sign convention used throughout. A right-handed set is used for every thin-walled question: the span axis runs along the member, $Y$ is vertical (upward) in the cross-section and $Z$ is horizontal. Section constants are written $I_{ZZ}=\int Z^{2}\,dA$, $I_{YY}=\int Y^{2}\,dA$ and $I_{YZ}=\int ZY\,dA$ about the centroid, and the bending stress is carried as the linear field $\sigma = aZ + bY$. Shear flows are counted positive in the direction of the walk stated with each figure; a negative answer simply means the flow runs the other way round the circuit.
Question 8: Torsion of a three-cell thin-walled wing box (20 marks)
Given. A rectangular three-cell box in pure torsion. Each cell is a closed loop bounded by an upper panel, a lower panel and two vertical webs; the two interior webs are shared between adjacent cells.
Quantity
Symbol
Value
Cell widths (left to right)
$w_1, w_2, w_3$
160, 300, 160 mm
Box depth
$h$
200 mm
Upper panel thickness
$t_u$
2.0 mm
Lower panel thickness
$t_l$
1.5 mm
Vertical web thickness
$t_w$
1.0 mm
Shear modulus
$G$
15 GPa
Applied torque (clockwise, constant)
$T$
10 000 N·m
Find. (a) the three cell shear flows, and (b) the largest shear stress in the section together with the wall in which it occurs.
Figure Q8 — the three-cell box with the converged cell flows. The interior webs carry only the difference between the flows of the cells they separate, which is why they are lightly stressed despite being the thinnest walls.
Approach. Write the rate of twist of each cell in terms of the three unknown cell flows, impose that all three cells twist at the same rate because the section is rigid in its own plane, and close the system with torque equilibrium; four equations, four unknowns.
Cell areas and the torque equation. The enclosed areas are $A_1 = A_3 = 160\times200 = 3.2\times10^{4}$ and $A_2 = 300\times200 = 6.0\times10^{4}$ mm2. For a multi-cell section the applied torque is shared according to $$T = 2\,(A_1q_1 + A_2q_2 + A_3q_3)$$with $T = 10\,000$ N·m $= 1.0\times10^{7}$ N·mm. This is one equation in three unknowns; the other two come from compatibility.
Rate of twist of a single cell. For cell $i$, $$\frac{d\theta}{dx} = \frac{1}{2A_iG}\oint \frac{q}{t}\,ds$$where the loop integral runs round the cell boundary and each wall carries the flow appropriate to it: an exterior wall carries $q_i$, while an interior web shared with cell $j$ carries the difference $q_i - q_j$.
Part (a) — assemble the compatibility equations. Writing out the loop integrals with the given gauges, cell 1 has $160/2 + 160/1.5 + 200/1 + 200/1 = 586.7$ mm of $ds/t$ on its own account and shares $200/1$ with cell 2. Doing the same for cells 2 and 3 and setting all three rates equal to a common $d\theta/dx$ gives three equations which, with the torque equation, form a four-by-four linear system in $q_1$, $q_2$, $q_3$ and $d\theta/dx$.
Solve the system. The box is symmetric about its vertical centre-line, so $q_1 = q_3$ before any arithmetic is done, which halves the work. Solving, $$\boxed{q_1 = q_3 = 34.72\ \text{N/mm},\qquad q_2 = 46.30\ \text{N/mm}}$$with a common rate of twist $d\theta/dx = 1.157\times10^{-5}$ rad/mm, equivalent to $0.663$ degrees per metre and an effective torsion constant $J = T/(G\,d\theta/dx) = 5.76\times10^{7}$ mm4.
Check the solution before using it. Substituting back, $2(A_1q_1+A_2q_2+A_3q_3) = 2(1.111\times10^{6}+2.778\times10^{6}+1.111\times10^{6}) = 1.000\times10^{7}$ N·mm, the applied torque exactly; and recomputing each cell’s twist rate independently returns the same $1.157\times10^{-5}$ rad/mm for all three. Both checks are worth doing because a sign error on a shared web is the classic failure mode of this calculation.
Part (b) — wall-by-wall shear stresses. Dividing each wall’s flow by its own thickness: the four outer webs and covers of cells 1 and 3 give $34.72/1.0 = 34.72$, $34.72/1.5 = 23.15$ and $34.72/2.0 = 17.36$ MPa; the cell-2 covers give $46.30/1.5 = 30.86$ and $46.30/2.0 = 23.15$ MPa; and the two interior webs carry only the difference $46.30-34.72 = 11.57$ N/mm, hence 11.57 MPa.
The governing wall. The largest value is therefore $$\boxed{\tau_{\max} = 34.72\ \text{MPa, in the two outer vertical webs (t = 1.0 mm)}}$$The result is worth a sentence of interpretation: the most heavily loaded flow is in the wide middle cell, but the most heavily stressed wall is an outer web, because it combines the full cell flow with the thinnest gauge. The interior webs, made of the same 1 mm sheet, are only a third as highly stressed because they see the flow difference rather than the flow itself.
Quantity
Result
(a) Shear flow in cell 1
34.72 N/mm
(a) Shear flow in cell 2
46.30 N/mm
(a) Shear flow in cell 3
34.72 N/mm
Flow carried by each interior web
11.57 N/mm
Rate of twist
$1.157\times10^{-5}$ rad/mm (0.663 degrees per metre)