22-Mec-B9 Advanced Engineering Structures · December 2017
Question 8 of 8: Shear centre and shear flow of a box with a semi-elliptical leading edge
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
Paper format. National Exams, December 2017 — 16-Mec-B9
Advanced Engineering Structures. Three hours, open book, any non-communicating
calculator permitted. The paper prints eight questions of equal value (20 marks each) and
states that any five constitute a complete exam paper. All eight are solved below,
because the set is a study resource rather than a sitting.
Reference texts.
T. H. G. Megson, Aircraft Structures for Engineering Students, 6th ed. —
unsymmetrical bending, shear flow in open and closed thin-walled sections, shear centre,
structural idealisation, single- and multi-cell torsion (Ch. 15, 16, 17, 20, 23).
A. C. Ugural and S. K. Fenster, Advanced Strength and Applied Elasticity, 5th ed.
— three-dimensional states of stress, principal stresses, the Tresca and von Mises yield
criteria (Ch. 1, 4).
A. P. Boresi and R. J. Schmidt, Advanced Mechanics of Materials, 6th ed. —
shear centre, torsion of multiply connected cells, column stability (Ch. 6, 12).
N. E. Dowling, Mechanical Behavior of Materials, 4th ed. — strain-life
fatigue, the Coffin–Manson relation, the Palmgren–Miner rule and Paris-law crack
growth (Ch. 9, 11, 14).
R. C. Hibbeler, Mechanics of Materials, 10th ed. — Euler buckling of
pin-ended columns, slenderness ratio and the transition slenderness (Ch. 13).
Sign convention used throughout. A single right-handed frame is used
for every thin-walled question: x runs along the span from the root to the free end,
Y is vertically upward and Z horizontally to the right, both measured from the
section centroid. Direct stress is written $\sigma_x = aZ + bY$, and the two coefficients follow
from $aI_{ZZ} + bI_{YZ} = M_Y$ and $aI_{YZ} + bI_{YY} = M_Z$, in which
$I_{ZZ}=\int Z^{2}\,dA$, $I_{YY}=\int Y^{2}\,dA$, $I_{YZ}=\int ZY\,dA$,
$M_Y \equiv \int \sigma_x Z\,dA$ and $M_Z \equiv \int \sigma_x Y\,dA$. This one pair of
equations replaces the memorised unsymmetrical-bending fraction and is self-checking: for a
cantilever carrying a tip load, the fibres on the side the load points toward must go into
compression.
Question 8: Shear centre and shear flow of a box with a semi-elliptical leading edge (20 marks)
Given. A four-boom single cell. Booms 1 and 4 sit on the rear
spar, 100 mm above and below the axis of symmetry; booms 2 and 3 sit on the front spar, 500 mm
forward of them and also 100 mm from the axis. The leading edge is a semi-ellipse whose minor
semi-axis is the 100 mm half-depth and whose major semi-axis is therefore 200 mm, so the nose
projects 200 mm forward of the 2–3 line.
Given data — four-boom nose box
Quantity
Value
Rear-spar booms 1 and 4
600 mm$^{2}$ each, at $y = \pm100$ mm
Front-spar booms 2 and 3
450 mm$^{2}$ each, at $y = \pm100$ mm
Spar-to-spar distance
500 mm
Semi-elliptical nose
major radius 200 mm, minor radius 100 mm
Wall thickness
constant (cancels throughout)
Applied shear force
10 000 N upward, 100 mm forward of the shear centre
Find. (a) the chordwise position of the shear centre; (b) the shear flow in
each of the four walls under the offset load.
The idealised nose box. The curved leading edge is handled exactly through its
swept area, $\oint p\,\mathrm{d}s = 2A_{\text{swept}}$, so there is no need to approximate it by
chords.
Approach. Evaluate $I$ from the boom areas, cut the rear spar and walk the
open flow boom by boom, close for zero twist using the wall lengths (the nose arc length is
needed, and follows from a numerical evaluation of the elliptic integral), locate the shear
centre by moments, then superpose the torque of the offset load.
Part (a) — evaluate the second moment. All four booms lie 100 mm from
the axis of symmetry, so
$$I_{ZZ} = \sum B_r y_r^{2} = \left(2 \times 600 + 2 \times 450\right)(100)^{2}
= 2.10\times10^{7}\ \text{mm}^{4}.$$
Measure the leading-edge arc. An elliptic arc has no closed-form length, so
the semi-perimeter of the ellipse with semi-axes 200 and 100 mm is integrated numerically,
giving $484.42$ mm. Ramanujan's approximation returns $484.42$ mm as well, agreeing to five
significant figures, which validates the numerical quadrature.
Cut the rear spar and walk the open shear flow. Stepping by
$-(S_y/I)B_r y_r$ at each boom on the walk $1\to2\to3\to4$,
$$q_{b,12} = -28.57, \qquad q_{b,23} = -50.00, \qquad q_{b,34} = -28.57, \qquad
q_{b,41} = 0 \ \ \text{N/mm}.$$
Close the cell for zero twist. With constant thickness the twist condition
is $\oint q\,\mathrm{d}s = 0$, and the four wall lengths are 500, 484.42, 500 and 200 mm, so
$$q_{s,0} = -\frac{\sum q_b \ell}{\sum \ell} = \frac{52\,792}{1684.4} = +31.34\ \text{N/mm},$$
giving zero-twist flows of $+2.77$, $-18.66$, $+2.77$ and $+31.34$ N/mm.
Take moments about the rear spar. Choosing the mid-point of the rear spar as
the reference kills the moment of the spar flow. For the straight walls twice the swept area is
$Z_1Y_2-Z_2Y_1$; for the nose the same integral evaluated round the ellipse gives
$162\,832$ mm$^{2}$, and the four contributions sum to $2A_{\text{cell}} = 262\,832$ mm$^{2}$,
matching $2\left[500(200)+\tfrac{\pi}{2}(200)(100)\right]$ exactly.
Locate the shear centre. The moment of the zero-twist flows about that
reference is $-2.761\times10^{6}\ \text{N}\cdot\text{mm}$, so
$$\boxed{z_{SC} = 276.1\ \text{mm forward of the rear spar}
\ \ (223.9\ \text{mm aft of the 2--3 line})}$$
The centre lies forward of mid-chord because the nose contributes almost a quarter of the
enclosed area while carrying comparatively little of the direct stress.
Part (b) — find the torque of the offset load. The 10 000 N force acts
100 mm forward of the shear centre, so about that centre it applies
$$T = 10\,000 \times (-100) = -1.00\times10^{6}\ \text{N}\cdot\text{mm},$$
the sign denoting a nose-down (clockwise) couple in the chosen frame.
Convert the torque into a circulating flow and superpose.
$q_T = T/2A_{\text{cell}} = -1.00\times10^{6}/262\,832 = -3.805$ N/mm, which adds uniformly to
every wall:
$$\boxed{q_{12} = -1.03,\quad q_{23} = -22.46,\quad q_{34} = -1.03,\quad
q_{41} = +27.54 \ \ \text{N/mm}}$$
Check the answer statically. Summing $q\,\Delta y$ over the four walls
returns $10\,000.0$ N upward, and the moment about the rear spar comes to
$-3.761\times10^{6}\ \text{N}\cdot\text{mm}$, exactly $10\,000 \times 376.1$ mm — the load
acting 376.1 mm forward of the spar, which is 276.1 + 100 mm as required.
Read the load path. The rear spar carries by far the largest flow, 27.54
N/mm, because the bending shear and the torsional circulation reinforce one another there, while
the two cover panels almost cancel to 1.03 N/mm. Moving the load only 100 mm forward of the shear
centre is enough to change the nose flow by 20 % — a reminder of how sensitive a
single-cell nose box is to the chordwise position of the load.
Final results — Question 8
Quantity
Value
Second moment of area
$I = 2.10\times10^{7}$ mm$^{4}$
Leading-edge arc length
484.42 mm
Enclosed area
$2A_{\text{cell}} = 2.6283\times10^{5}$ mm$^{2}$
(a) Shear centre
276.1 mm forward of the rear spar
(223.9 mm aft of the front spar)