Integraph

RC column design

Worked example: check a reinforced concrete column under axial load and biaxial bending against the N-M interaction surface, with hand verification of the squash load, the balanced point and the Bresler reciprocal.

Problem statement

An interior column in a multi-storey office building carries combined axial compression and biaxial bending from frame action. The column is 400 ×\times 400 mm and 3.5 m tall between floor levels. The exposure classification is A1 (interior, non-aggressive).

Check the column section capacity under the critical ULS load combination to AS 3600:2018.

Given data

ParameterValueUnitsSource
Column width (bb)400mmGiven
Column depth (DD)400mmGiven
Column height (LL)3500mmGiven
Concrete gradeN50AS 3600
fcf'_c50MPaN50 grade
Rebar gradeD500NAS/NZS 4671
fsyf_{sy}500MPaD500N
Cover35mmChosen — see the note below
Fitment sizeN10Assumed
Fitment spacing300mmAssumed
Reinforcement8-N248 bars evenly around the perimeter

On the cover. AS 3600:2018 Table 4.10.3.2 gives 20 mm for exposure classification A1 at every strength grade, so 35 mm is not that table’s number. It is a deliberate choice above the minimum: Cl 4.10.2 requires cover no less than the bar diameter, which is 24 mm here, and 35 mm leaves margin for the fitment cage and placing tolerance. The distinction matters because the cover sets the bar inset, and the bar inset moves every capacity on this page.

Design actions (critical ULS combination)

LoadValueUnits
NN^*2500kN
MxM^*_x120kN.m
MyM^*_y80kN.m

NN^* is positive in compression, MxM^*_x positive puts the top face in compression and MyM^*_y positive puts the left face in compression, so the two together put the peak compression in the top-left quadrant. See Platform — Sign and axis conventions.

These actions include second-order effects (moment magnification per AS 3600 Cl. 10.4 has been applied externally). The eccentricities they imply are worth writing down now, because one of the three biaxial methods is built on them:

ex=MxN=1202500=48 mm,ey=MyN=802500=32 mme_x = \frac{M^*_x}{N^*} = \frac{120}{2500} = 48 \text{ mm}, \qquad e_y = \frac{M^*_y}{N^*} = \frac{80}{2500} = 32 \text{ mm}

Reinforcement details

  • 8-N24 bars: Ast=8×452.4=3619A_{st} = 8 \times 452.4 = 3619 mm2^2
  • Reinforcement ratio: ρ=3619/(400×400)=2.26%\rho = 3619 / (400 \times 400) = 2.26\% — inside the Cl 10.7.1 range of 1%—4%
  • Bar inset from each face: 35+10+24/2=5735 + 10 + 24/2 = 57 mm, so the bars sit at ±143\pm 143 mm and 00 about the section centroid
  • Depth to the outermost bar layer: d=40057=343d = 400 - 57 = 343 mm

Step-by-step solution

Step 1: Define section geometry

Apply the Rectangular template with b=400b = 400 mm, D=400D = 400 mm.

Step 2: Set materials and member type

On the General tab:

  • Design Standard: AS 3600
  • Member Type: Column
  • Under MaterialsConcrete grade: N50 (fc=50f'_c = 50 MPa), Rebar grade: D500N (fsy=500f_{sy} = 500 MPa)

On the Rebar/PT tab, under Cover:

  • Cover: 35 mm all sides

On the Settings tab (the gear icon), under Analysis models:

  • ULS model: Rectangular — the AS 3600 Cl 8.1.3 block, which is what the interaction sweep integrates and what the hand checks at the end of this page derive

Member Type is not cosmetic. It selects which checks run and which AS 3600 Section 5 axis distance applies, and it is what makes the Cl 10.1.2 minimum-moment row appear at all. ACS refuses to analyse a section that has none saved rather than assuming “beam” on your behalf. At N=2500N^* = 2500 kN against the 0.1fcAg=8000.1 f'_c A_g = 800 kN threshold, the load regime agrees with the declaration, so no mismatch note renders.

Step 3: Place reinforcement

Use the Perimeter Pattern tool:

  • Bar diameter: 24 mm (N24)
  • Number of bars: 8 (3 per face with corner sharing)

This distributes 8 bars evenly around the perimeter at the cover + fitment inset — three at d=57d = 57 mm, two at d=200d = 200 mm, three at d=343d = 343 mm.

Then the fitments. This section is a convex rectangle, so the whole outline is a single cell and the Cell Tie tool opens its dialog with no cell to pick first. In the dialog:

  • Stirrup Diameter: 10 mm (N10)
  • Grade: D500N
  • Spacing (mm): 300
  • Hook ends: 135° both ends

There is no “legs” input. AsvA_{sv} is derived from the tie geometry you drew — a closed perimeter tie on a rectangle crosses any horizontal cut twice. Nothing on this page depends on it (the combination carries no shear), but the cage has to exist for the bars to be placed against.

Step 4: Enter the design combination

Loads live on the Actions tab, under the ULS combinations heading. Add one row:

CombinationLimit statekϕk_\phiNN^*MxM^*_xMyM^*_yVyV^*_yVxV^*_x
ULS-1ULS12/132500 kN120 kN.m80 kN.m00

The kϕk_\phi column is the one to read twice. AS 3600 Amd 2:2021 Table 2.2.2 has two separate entries for a column, and they do not carry the same ϕ\phi:

  • item (a)(ii), pure axial compression — a fixed ϕ=0.65\phi = 0.65, with no kϕk_\phi in it
  • item (d), bending with axial compression — a floor of ϕo=0.65kϕ\phi_o = 0.65 \, k_\phi

So kϕk_\phi sets ϕ\phi everywhere on the interaction curve except its squash endpoint. The grid offers 12/13 and 1.0, and 12/13 is both the default it renders and the value it saves into the design. 1.0 is available only where Cl 10.3 makes the column short and Q/G0.25Q/G \geq 0.25; taking it otherwise is unconservative. This example designs at 12/13, which gives ϕo=0.65×12/13=0.60\phi_o = 0.65 \times 12/13 = 0.60.

API and MCP callers get no default at all: an AS 3600 ULS combination with no kϕk_\phi saved is refused with k_phi_not_set rather than analysed at 0.60 silently.

Because the member type is Column, AS 3600 Cl. 10.1.2 imposes a minimum design moment about each principal axis:

Mmin=N×0.05D=2500×0.05×0.400=50 kN.mM^*_{\min} = N^* \times 0.05 D = 2500 \times 0.05 \times 0.400 = 50 \text{ kN.m}

Both entered moments exceed it, so the floor does not govern and the results below are for the moments as entered. Where a floor does govern, the panel says so explicitly (”MyM^*_y taken as … — minimum design moment per AS 3600 Cl. 10.1.2”) rather than substituting silently.

Step 5: Review the interaction results

Switch to the ULS results tab and find the Interaction Diagram section.

That section is a set of readouts, not a plot. The plotted envelope lives on the canvas — the Interaction Surface tab, whose 2D / 3D toggle sits at the bottom of the canvas. Step 6 covers both views.

What the panel reports

The panel shows two axial end-caps and the biaxial utilisation checks. It does not report a uniaxial NN-MxM_x utilisation: that row was removed deliberately, because MyM^*_y was never one of its arguments and it therefore read lower than the true demand whenever there was out-of-plane bending. Biaxial (rigorous) is the utilisation for this envelope.

ReadoutValueUnits
Nu0N_{u0} (design squash, ϕ=0.65\phi = 0.65)5496.2kN
NtN_t (design pure tension)-1538.1kN
Biaxial (rigorous)0.747
Bresler contour (αn\alpha_n = 1.473)0.533
Bresler reciprocal0.905

All three are below 1.0, so the section is adequate. The Bresler reciprocal governs at 0.905, and the order the three come in is the subject of the Discussion below — it is not the order you might expect.

The panel’s contour label rounds the exponent to one decimal, so for this section it renders “Bresler contour (α=1.5)”. That 1.5 is a coincidence of rounding, not the method’s constant — see the α_n note under Hand calculation verification below.

Features of the curve

These are not on screen in the editor. They come from the design report, section N-M Interaction Diagram, whose Key Points table lists each labelled point with the ϕ\phi it was factored by, above a chart plotting the nominal and design NN-MxM_x curves together. The canvas plots the design envelope only — there is no nominal surface in the engine, and the 2D curve marks no point but ϕNu0\phi N_{u0}.

The nominal column is the unfactored capacity; the design column carries ϕ\phi from Table 2.2.2 at each point.

Curve featureNominalDesign (kϕk_\phi = 12/13)Units
Squash, Nu0N_{u0}8455.75496.2 (ϕ=0.65\phi = 0.65)kN
Balanced point(2241.0, 466.6)(1346.3, 280.3)(kN, kN.m)
Pure bending, MuoM_{uo}283.6241.0 (ϕ=0.85\phi = 0.85)kN.m
Pure tension, NtN_t-1809.6-1538.1kN

Two things to notice in that table, both of them ϕ\phi rather than mechanics:

  • The squash load is the same under either kϕk_\phi class, because Table 2.2.2(a)(ii) pins ϕ=0.65\phi = 0.65 at the pure-axial endpoint and carries no kϕk_\phi.
  • Pure bending takes ϕ=0.85\phi = 0.85, not 0.60 — Table 2.2.2(b), because at N=0N = 0 there is no compression-controlled floor to apply. Its ductility parameter is kuo=k_{uo} = 0.218, comfortably under the 0.36 limit.

N=2500N^* = 2500 kN sits above the nominal balanced axial load of 2241 kN, so the section is compression-controlled: the concrete crushes before the tension steel yields, and adding reinforcement buys less than it would on a beam.

Step 6: View the interaction envelope on the canvas

Click the Interaction Surface tab on the canvas. The toggle at the bottom of the canvas switches between two views of the same envelope.

3D (NN-MxM_x-MyM_y surface). Left-drag to rotate the surface and see the capacity envelope from different angles; scroll to move in and out, and middle-drag (or Shift / Ctrl / ⌘ + left-drag) to pan. The MxM_x-MyM_y contour at the applied axial load level shows the remaining moment capacity in all directions — and it is where the difference between the three methods below becomes visible rather than numerical.

2D (one curve, one combination). The curve is swept in the resultant-moment plane at this combination’s own bending angle θ\theta, so the x-axis is MR=Mx2+My2M_R = \sqrt{M_x^2 + M_y^2} and not MxM_x. The design point’s radial position against the curve is the printed utilisation, by construction. One combination at a time is deliberate: each has its own θ\theta and therefore its own curve, and a point overlaid on a curve from another plane could read comfortably inside an envelope that is not its own.

The design-point marker is coloured by verdict, not by identity: green when the utilisation is at or below 1.0, red when it is above, and a neutral colour when no ratio was returned at all. This combination passes at 0.747, so the marker renders green — a red marker here would mean the section had failed, not that you had found the design point.

Results summary

CheckDemandCapacityUtilisationStatus
Biaxial (rigorous)(2500, 120, 80)3D surface0.747Pass
Bresler reciprocalN=2500N^* = 2500 kNNuN_{u} = 2761.6 kN0.905Pass
Bresler contour (αn=1.473\alpha_n = 1.473)(120, 80) kN.m247.6 kN.m each axis0.533Pass
Reinforcement ratio (Cl 10.7.1)2.26%1%—4%Pass

The section is adequate as detailed: 400 ×\times 400 N50, 8-N24, N10-300 fitments. The Bresler reciprocal governs at 0.905, which leaves about 10% reserve — materially tighter than the rigorous surface’s 0.747 suggests.

Discussion

The three methods do not bracket each other

Methodkϕk_\phi = 12/13kϕk_\phi = 1.0
Bresler contour0.5330.452
Biaxial (rigorous)0.7470.661
Bresler reciprocal0.9050.801

The ordering is contour < rigorous < reciprocal, and it holds under both kϕk_\phi classes. The Bresler methods are not “slightly conservative” relative to the rigorous surface: one of them is markedly less conservative and the other markedly more. They are different approximations to the same surface, built from different read-offs, and for a square section with symmetric reinforcement at this eccentricity ratio they land 70% apart.

That matters because a designer who checks the contour row, sees 0.533, and stops has read the most generous of the three. Design to the governing row — here the reciprocal at 0.905 — and use the rigorous surface as the physical answer the two approximations are trying to estimate.

What kϕk_\phi is worth

Every utilisation in the table above falls on moving from 12/13 to 1.0 — the two biaxial rows by about 11.5%, the Bresler contour row by 15.2% — because ϕo\phi_o rises from 0.60 to 0.65 across the whole compression-controlled region of the curve. That is not a rounding matter, and 1.0 is not available by default: Cl 10.3 has to make the column short and Q/GQ/G has to reach 0.25. If neither has been established, 12/13 is the classification, and the app saves it rather than leaving the question open.

Other observations

  • The reinforcement ratio of 2.26% is within the Cl 10.7.1 range for columns (1%—4%).
  • For a symmetric section with symmetric reinforcement the Bresler methods are cheap sanity checks; for asymmetric sections the rigorous 3D surface is the only one of the three that reflects the actual geometry.

If the utilisation were too high, in order of effectiveness for a compression-dominated section:

  • Increase the column size — Nu0N_{u0} scales with AgA_g and the lever arms grow with DD
  • Increase the concrete grade — Nu0N_{u0} scales close to directly with fcf'_c
  • Add reinforcement (up to the 4% Cl 10.7.1 maximum) — the weakest of the three here, because the section is above the balanced point and the extra steel is not yielding

Hand calculation verification

Squash load (ϕNu0\phi N_{u0})

AS 3600:2018 Cl 10.6.2.2. Note the uniform-compression intensity factor α1\alpha_1, which is a different quantity from the bending block’s α2\alpha_2:

α1=clamp(10.003fc, 0.72, 0.85)=clamp(0.85, 0.72, 0.85)=0.85\alpha_1 = \mathrm{clamp}(1 - 0.003 f'_c,\ 0.72,\ 0.85) = \mathrm{clamp}(0.85,\ 0.72,\ 0.85) = 0.85 σs=min(fsy, Es×0.0025)=min(500, 500)=500 MPa\sigma_s = \min(f_{sy},\ E_s \times 0.0025) = \min(500,\ 500) = 500 \text{ MPa} Ac=40023619=156,381 mm2A_c = 400^2 - 3619 = 156{,}381 \text{ mm}^2 Nu0=0.85×50×156,381+500×3619=6646+1810=8456 kNN_{u0} = 0.85 \times 50 \times 156{,}381 + 500 \times 3619 = 6646 + 1810 = 8456 \text{ kN} ϕNu0=0.65×8456=5496 kN\phi N_{u0} = 0.65 \times 8456 = 5496 \text{ kN}

Matches the reported 8455.7 / 5496.2 kN. The ϕ\phi here is Table 2.2.2(a)(ii)‘s fixed 0.65 — kϕk_\phi does not enter.

Balanced point

The balanced state has the extreme tension bar reaching εsy=0.0025\varepsilon_{sy} = 0.0025 as the extreme compression fibre reaches εcu=0.003\varepsilon_{cu} = 0.003:

kub=0.0030.003+0.0025=0.5455,c=0.5455×343=187 mmk_{ub} = \frac{0.003}{0.003 + 0.0025} = 0.5455, \qquad c = 0.5455 \times 343 = 187 \text{ mm}

With the Cl 8.1.3 block at fc=50f'_c = 50: α2=0.850.0015×50=0.775\alpha_2 = 0.85 - 0.0015 \times 50 = 0.775, γ=0.970.0025×50=0.845\gamma = 0.97 - 0.0025 \times 50 = 0.845, so a=γc=158a = \gamma c = 158 mm. Taking moments about the gross centroid, with the three bar rows at y=+143y = +143, 00 and 143-143 mm:

Cc   = 0.775 x 50 x 400 x 158.1 = 2451 kN   at arm 121 mm
row 1 (3 bars, y = +143, compression, net of displaced concrete)  = +514 kN at arm +143 mm
row 2 (2 bars, y = 0, elastic tension)                            =   -38 kN at arm 0
row 3 (3 bars, y = -143, yielded tension)                         =  -679 kN at arm -143 mm

N = 2451 + 514 - 38 - 679 = 2248 kN
M = 2451(0.121) + 514(0.143) + 679(0.143)  = 467 kN.m

Against the reported nominal balanced point of (2241.0, 466.6) that is 0.3% on both ordinates — two independent constructions of the same state.

Pure bending (ϕMuo\phi M_{uo})

At N=0N = 0 the solve returns c=74.8c = 74.8 mm (kuo=74.8/343=0.218k_{uo} = 74.8/343 = 0.218), so a=0.845×74.8=63.2a = 0.845 \times 74.8 = 63.2 mm. The neutral axis sits at y=20074.8=125y = 200 - 74.8 = 125 mm, which puts the y=+143y = +143 row in compression at a strain of 0.003×(143125)/74.8=0.000710.003 \times (143-125)/74.8 = 0.00071 — elastic, at σ=143\sigma = 143 MPa:

Cc                  = 0.775 x 50 x 400 x 63.2       = 980 kN  at arm +168 mm
row 1 (compression, net of displaced concrete)      = +141 kN  at arm +143 mm
row 2 (2 bars, y = 0, yielded tension)              = -452 kN  at arm 0
row 3 (3 bars, y = -143, yielded tension)           = -679 kN  at arm -143 mm

N = 980 + 141 - 452 - 679 = -10 kN     (~0 as required)
M = 980(0.168) + 141(0.143) + 679(0.143) = 282 kN.m

Against the reported Muo=283.6M_{uo} = 283.6 kN.m that is 0.5%. ϕMuo=0.85×283.6=241.0\phi M_{uo} = 0.85 \times 283.6 = 241.0 kN.m — and note again that the 0.85 is Table 2.2.2(b), because at N=0N = 0 the item-(d) floor does not apply.

Bresler reciprocal

The reciprocal is a constant-eccentricity construction, which is what makes it a hand check you can actually reproduce. NuxN_{ux} and NuyN_{uy} are the uniaxial capacities on the radial line through the design point in each plane — at ex=48e_x = 48 mm and ey=32e_y = 32 mm respectively — not the capacities read off the diagram at constant moment. Reading them at constant moment gives numbers 15%—20% high and an unconservative utilisation.

From the envelope at those two eccentricities, with the squash load from above:

TermValueUnits
NuxN_{ux} at ex=48e_x = 48 mm3510.1kN
NuyN_{uy} at ey=32e_y = 32 mm3858.6kN
Nu0N_{u0}5496.2kN
NuN_u (reciprocal)2761.6kN
1Nu=13510.1+13858.615496.2=2.849×104+2.592×1041.819×104=3.621×104\frac{1}{N_u} = \frac{1}{3510.1} + \frac{1}{3858.6} - \frac{1}{5496.2} = 2.849\times10^{-4} + 2.592\times10^{-4} - 1.819\times10^{-4} = 3.621\times10^{-4} Nu=2762 kNNNu=25002761.6=0.905N_u = 2762 \text{ kN} \qquad \Rightarrow \qquad \frac{N^*}{N_u} = \frac{2500}{2761.6} = 0.905

That confirms the reciprocal row — 0.905 — and nothing else. Comparing a reciprocal hand check against the rigorous surface’s 0.747 and calling the agreement confirmation would be a category error: they are different constructions and they are supposed to differ.

Bresler contour exponent

The contour exponent is computed per AS 3600:2018 Cl 10.6.4, not taken as a constant. The clause pins ϕ=0.65\phi = 0.65 inside its own definition, so αn\alpha_n does not move with kϕk_\phi:

αn=0.7+1.7NϕNu0=0.7+1.7×25000.65×8456=1.473,1αn2\alpha_n = 0.7 + \frac{1.7 N^*}{\phi N_{u0}} = 0.7 + \frac{1.7 \times 2500}{0.65 \times 8456} = 1.473, \qquad 1 \leq \alpha_n \leq 2 (120247.6)1.473+(80247.6)1.473=0.344+0.189=0.533\left(\frac{120}{247.6}\right)^{1.473} + \left(\frac{80}{247.6}\right)^{1.473} = 0.344 + 0.189 = 0.533

with the uniaxial moment capacities at NN^* both 247.6 kN.m (equal, because the section is square and symmetrically reinforced). Matches the reported 0.533.

A fixed α=1.5\alpha = 1.5 is the ACI 318 value. Under AS 3600 — and under EN 1992 — the exponent is a function of the axial utilisation, and it happens to round to 1.5 for this particular column at this particular NN^*. It will not for the next one.