Question 4 of 6: Vertical Stress Increase by Superposition
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
EGBC / Engineers Canada National Examination — Structural Engineering (legacy), 07-Str-A3 Geotechnical Materials and Analysis, May 2013. Three hours; closed book; one Casio or Sharp approved calculator; drawing instruments required; all charts and equations are supplied at the back of the paper. Total value 100 marks over six compulsory questions (20 + 10 + 10 + 20 + 20 + 20). Every question and every sub-part is solved in full below.
Reference texts: Das, Principles of Geotechnical Engineering, 9th ed. (Ch. 2 grain-size distribution and classification, Ch. 6 compaction, Ch. 7 permeability, Ch. 8 seepage and flow nets, Ch. 9 in-situ stresses and capillary rise, Ch. 10 stresses in a soil mass, Ch. 11 consolidation, Ch. 12 shear strength); Knappett & Craig, Craig’s Soil Mechanics, 8th ed. (flow nets, effective stress, heave and piping); Holtz, Kovacs & Sheahan, An Introduction to Geotechnical Engineering, 2nd ed.; Canadian Geotechnical Society, Canadian Foundation Engineering Manual (CFEM), 4th ed. — the governing Canadian practice document for settlement, excavation dewatering and factors of safety; ASTM D2487 (USCS classification) and ASTM D6913 (grain-size analysis).
Question 4: Vertical Stress Increase by Superposition (20 marks)
Given. The plan of Figure 3 sits inside a 50 m × 50 m envelope with point A at its north-east corner; the two flexible foundations and the depth of interest are listed below.
Given data — Question 4
Quantity
Symbol
Value
Contact stress, multiplex M1
q1
50 kPa
Contact stress, multiplex M2
q2
60 kPa
Depth of the point of interest below A
z
5 m
M1 footprint (two rectangles)
—
20 m × 40 m west column, plus a 50 m × 10 m south strip
M2 footprint (two rectangles)
—
10 m × 38 m east column, plus a 10 m × 10 m arm
Newmark chart influence value (from the supplied chart)
IN
0.005 (200 elements)
Find. The total increase in vertical stress $\Delta\sigma_z = \Delta\sigma_{z1} + \Delta\sigma_{z2}$ at 5 m depth vertically below point A.
[Figure not reproduced: Figure 3 (redrawn to scale) — plan of the two multiplex footprints. Origin at the south-west corner, x east and y north; A is at (50, 50). See the official exam paper.]
Approach. Boussinesq’s solution is linear in the applied load, so any footprint can be built up as a signed sum of rectangles that each have one corner vertically above the point of interest. For M1 the four such rectangles are read off the m–n chart; for M2 the irregular L-shape is handled with the Newmark influence chart as instructed, and the count is checked against the equivalent corner-rectangle solution.
Fix a frame and express both footprints as corner rectangles. Take the origin at the south-west corner of the 50 m square, $x$ east and $y$ north, so $A = (50,\,50)$. Write each rectangle in terms of its distance west of A ($a = 50 - x$) and south of A ($b = 50 - y$), because a rectangle described that way automatically has a corner over A. The M1 footprint is the union of the south strip $a \in [0,50],\ b \in [40,50]$ and the west column $a \in [30,50],\ b \in [10,40]$; the M2 footprint is the union of the east column $a \in [0,10],\ b \in [0,38]$ and the arm $a \in [10,20],\ b \in [20,30]$.
Write the signed superposition for M1. Denoting by $R(a,b)$ the influence factor of a rectangle $a \times b$ with its corner over A, the strip is $R(50,50) - R(50,40)$ and the column is $R(50,40) - R(50,10) - R(30,40) + R(30,10)$. The $R(50,40)$ terms cancel, leaving $$I_{M1} = R(50,50) - R(50,10) - R(30,40) + R(30,10)$$
Convert each rectangle to m and n at z = 5 m. With $m = B/z$ and $n = L/z$ (interchangeable), the four rectangles and their chart values are tabulated below. The values are the standard Boussinesq corner factors and can be read directly off the m–n chart supplied with the paper.
m–n superposition for M1 (z = 5 m, q1 = 50 kPa)
Rectangle a × b (m)
Sign
m = a/z
n = b/z
I
Signed I
50 × 50
+
10.0
10.0
0.2498
+0.2498
50 × 10
−
10.0
2.0
0.2398
−0.2398
30 × 40
−
6.0
8.0
0.2494
−0.2494
30 × 10
+
6.0
2.0
0.2397
+0.2397
Σ IM1
+0.00028
Evaluate the M1 contribution. Multiplying the net influence factor by the contact stress, $$\Delta\sigma_{z1} = q_1\,I_{M1} = 50 \times 0.00028 = 0.014\ \text{kPa}$$ so that $$\boxed{\;\Delta\sigma_{z1} \approx 0.01\ \text{kPa} \approx 0\;}$$ The four chart factors are each close to the limiting value 0.25 because every one of these rectangles is large compared with $z = 5$ m; the physics is carried entirely by the differences between them, and those differences very nearly cancel.
Interpret that result before moving on. The nearest loaded element of M1 to point A is the north-east corner of the south strip, at $(50,\,10)$, which is 40 m away in plan. The Boussinesq point-load solution decays as $\sigma_z \propto z^3/R^5$, and the familiar 2:1 construction spreads load only about one horizontal unit for every two vertical, so at a depth of only 5 m the stress bulb of M1 has long since died out at A. The superposition is not wasted work — it is the calculation that demonstrates the answer — but the engineering conclusion is that M1 has no practical influence on the stress 5 m below A.
Set up the Newmark chart for M2. The chart supplied has an influence value $I_N = 0.005$, that is 200 influence elements. Draw the plan of M2 on tracing paper to the scale set by the chart’s depth line, $\overline{AB} = z = 5\ \text{m}$, lay the tracing over the chart with point A on the centre of the circles, and count the elements covered by the loaded area.
Count the elements. The 10 m × 38 m east column has one corner exactly at A and, being only two depth-units wide, covers a compact block of elements around the centre; the 10 m × 10 m arm lies between two and four depth-units away and covers a negligible fraction of an element. Counting gives $$N \approx 48\ \text{elements}$$
Check the count against the corner-rectangle solution. A Newmark count is only as good as the tracing, so verify it. The east column has its corner at A with $m = 10/5 = 2.0$ and $n = 38/5 = 7.6$, giving $I = 0.2398$; the arm contributes $R(20,30) - R(20,20) - R(10,30) + R(10,20) = 0.0003$. Hence $I_{M2} = 0.2401$ and $\Delta\sigma_{z2} = 60 \times 0.2401 = 14.4\ \text{kPa}$, which back-figures to $N = 14.4/(0.005 \times 60) = 48.0$ elements — the count is confirmed.
Superpose. Adding the two contributions, $$\Delta\sigma_z = \Delta\sigma_{z1} + \Delta\sigma_{z2} = 0.014 + 14.406$$ $$\boxed{\;\Delta\sigma_z \approx 14.4\ \text{kPa at 5 m below point A}\;}$$
Question 4 — increase in vertical stress 5 m below A
Contribution
Method
Influence factor
Δσz (kPa)
Multiplex M1 (q1 = 50 kPa)
m–n corner rectangles
ΣI = 0.00028
0.014 ≈ 0
Multiplex M2 (q2 = 60 kPa)
Newmark chart, N ≈ 48, IN = 0.005
equivalent I = 0.2401
14.4
Total by superposition
—
—
14.4
Check: two readings of Figure 3 are worth recording. First, the right-hand dimension chain reads 38 m + 10 m against an overall height of 50 m, which leaves a 2 m unloaded strip between the base of the M2 column and the top of the M1 strip; the figure has been carried as drawn. Closing that gap changes Δσz1 by less than 0.01 kPa and leaves the total unchanged at 14.4 kPa. Second, point A is taken at the north-east corner of the 50 m envelope, exactly on the corner of the M2 column, as the dot in the figure shows.