Beam on elastic foundation
A reinforced-concrete grade beam (or continuous footing) of finite length resting on the soil, analysed as a beam on a Winkler foundation and designed to ACI 318-19 with ASCE 7-22 loads. The soil is a bed of independent springs: the contact pressure at a point is the modulus of subgrade reaction ks times the settlement there, so a stiff beam spreads its loads and a flexible one concentrates the pressure under them. The calculation gives the settlement, soil pressure, shear and moment along the beam for every load combination, checks the peak service pressure against the allowable bearing pressure and the settlement against your limit, and designs the concrete section for flexure and shear.
Inputs of the beam on elastic foundation #
- Beam: length L, width b in contact with the soil, depth h, and whether the self-weight is added to the dead load. The beam is free at both ends: the soil alone holds it up (lengths up to 80 ft, 24 m).
- Loads, unfactored and per load type (D, L, Lr, S, R, W, E): point loads P (downward positive) with an optional couple M (clockwise positive) at any position x from the left end, and line loads, uniform or linearly varying, over the full length or any part of it (From = To = 0 loads the full length). The ASCE 7-22 2.4.1 combinations are used for the soil pressure and the settlement, 2.3.1 for the concrete; wind and earthquake act in both directions.
- Modulus of subgrade reaction ks. It is not a soil constant: it depends on the loaded width and shape and on the soil, so it is your input, with the value the geotechnical engineer gives for this beam. When you only have a 1 ft (0.3 m) plate-load test value ks1, choose "From a 1 ft plate-load test": the calculation converts it to the beam width with Terzaghi's (1955) expressions, ks = ks1 ((B + 1)/2B)² for granular soil and ks = ks1 (1/B)(m + 0.5)/(1.5 m) for cohesive soil (B = b in ft, m = L/B), and shows the steps. Nothing is converted unless you choose it.
- Soil springs: compression only (the default: the soil cannot pull the beam down, so the springs that would be in tension are removed and the beam re-analysed until the contact zone settles, per combination) or linear (the classical Hetényi solution; any tension in the springs is then reported as a failed contact check).
- Allowable bearing pressure (gross, from the geotechnical report) and the settlement limit (your limit, commonly 1 in. total for buildings).
- Beam stiffness: Ec by ACI 318-19 19.2.2.1(b) (57,000√f'c), 19.2.2.1(a) (wc^1.5 × 33√f'c) or entered; the moment of inertia gross, a fraction of Ig, or entered. A stiffer beam attracts more moment and spreads the pressure; a cracked (softer) beam concentrates the pressure under the loads: Table 6.6.3.1.1(a) gives 0.35Ig for beams at factored loads and 6.6.3.2.2 lets you use 1.4 times that (0.49Ig) for service loads. Consider both; the calculation tells you when the service moment exceeds the cracking moment.
- Concrete and reinforcement: f'c, fy, λ, clear covers (3 in. cast against earth at the bottom), top and bottom bars (size and number, continuous over the length), stirrups (size, legs, spacing) or none.
Auto-design (editor toolbar) varies the width (4-in. steps) and depth (2-in. steps), bottom and top bars (size and count), stirrup size, spacing and legs; the length, loads, soil and limits are kept. Objectives: lowest cost index (concrete, formwork, reinforcement); least concrete, then cost. Every other input is kept; see Auto-design.
Checks performed #
| Check | Basis |
|---|---|
| Analysis | Beam finite elements with consistent Winkler springs, nodes at every load, at least 60 elements and λ·le ≤ 0.05; λ = (k/4EI)^¼ with k = ks·b; λL and Hetényi's classification (short and rigid below π/4, long above π) reported |
| Soil bearing pressure | Largest service pressure ks·y under the 2.4.1 combinations ≤ the allowable (gross) |
| Contact | Linear springs: no tension in the soil; compression-only springs: the length that lifts off is reported; a net upward load is not feasible |
| Settlement | Largest service settlement ≤ your limit |
| Flexure, top and bottom bars | 22.2 with φ by εt (Table 21.2.2); εt ≥ 0.004 (9.3.3.1); As,min 9.6.1.2 or 4/3 of the required area (9.6.1.3) |
| Bar spacing | Clear spacing 25.2.1 (3/4 in. aggregate); crack control 24.3.2 with fs = 2/3 fy |
| One-way shear | 22.5 (Table 22.5.5.1, λs; Vs = Av fyt d/s), section limit 22.5.1.2, Av,min 9.6.3.1, stirrup spacing along and across 9.7.6.2.2; the largest factored shear anywhere along the beam is used (no reduction at d from the loads) |
| Advisories | Development length of the top and bottom bars from the section of largest moment to the nearer end (25.4.2.4); service moment against the cracking moment (stiffness choice); skin reinforcement when h > 36 in. (9.7.2.3) |
The drawings (screen, report, DXF and PDF sheet) show the beam elevation on its soil springs with the loads by type (point loads keyed P1, P2, ... with their values in the sheet notes when there are more than three positions), the load positions and the length, the service soil pressure envelope (with the settlement at the labelled peaks) and the factored shear and moment envelopes under it, and the typical section with the bars, stirrups and cover. Bars are drawn continuous over the length as designed; cut-offs, hooks, splices and column dowels are not part of the calculation and the sheet notes say so.
Limitations #
- Winkler model: the springs are independent, so the soil next to a loaded area does not settle with it (no shear coupling, no elastic half-space); the results are only as good as ks. Vary ks (for example half and twice the value) when it is uncertain.
- Constant section and constant ks along the beam; no axial load, no torsion, no supports other than the soil (piles or walls under the beam are not modelled); the length is limited to 80 ft (24 m).
- Shear is checked at its largest value anywhere, conservatively; punching and bearing under the columns are not checked (use the spread or combined footing calculations where they govern).
- The plate-load conversion follows Terzaghi as given in the textbooks (Bowles, Foundation Analysis and Design, 9-6); the original paper is not on file and the value should be confirmed by the geotechnical engineer.
Verification examples #
Verification status: 20 values compared in 7 examples (5 published, 2 hand calculations from the cited equations): 18 within tolerance, 2 recorded deviations. Code basis: partial. Only some provisions were read in the target edition; the rest rests on older editions, examples or summaries.
backend/tests/test_verify_elastic_foundation_beam.py; the finite elements are also compared with Hetényi's closed forms for the infinite beam (force and couple), the semi-infinite beam and the finite beam with free ends in backend/tests/test_elastic_foundation_beam.py (within 0.1 %).
| Example | What it checks | Published | ILX Calcs |
|---|---|---|---|
| Hetényi (1946) §18: 120 in. wood beam, free ends, 5000 lb and a partial uniform load | Deflection at A / C / B | 0.03036 / 0.05193 / 0.00628 in. | 0.03035 / 0.05197 / 0.00632 |
| Hetényi §18 | Foundation pressure at A / C / B | 6.07 / 10.39 / 1.26 psi | 6.07 / 10.39 / 1.26 |
| Hetényi §18 | Moment and shear at the middle C | 9623.5 lb-in. / 75.3 lb | 9574 / 74.4 (the publication sums interpolated table values) |
| Hetényi §12, Ex. 1: equal loads spaced π/2λ on a long beam | Moment under a load / (P/4λ): two loads; middle and outer of three | 0.7921; 0.5842 and 0.7489 | 0.792; 0.584 and 0.749 |
| CodeCogs rail example (10 tons per rail, k = 400 lb/in.²) | 1/λ; half the depressed length; maximum moment | 54.8 in.; 129 in.; 137 ton-in. | 54.8; 129.1; 136.9 |
| Structville: 10 m × 0.4 m × 0.6 m concrete beam, 300 kN at 3 m | λ | 0.2828 per m | 0.2828 |
| Structville | End deflection; moment under the load (deviation: the hand calculation rounds Roark's constants) | 12.23 mm; 228.4 kN-m | 12.39; 231.0 (exact, as an independent closed-form solution) |
| Default design (hand) | φMn of 6 #8 top bars; As,min by 9.6.1.2 | 679.1 kip-ft; 3.96 in.² | 679.1; 3.96 |
| Plate-load conversion (hand) | ks from ks1 = 250 pci on sand, B = 3 ft | 111.1 pci | 111.1 |