a913a3d0-c8b4-4f4a-a1f8-4453af39f831 30 4 5 decimal

&[DATE] &[TIME] - &[FILENAME]

Calculation of biaxial demand (P, Mx, My) based on assumed responce

biaxial(α#,na#)
α# - orientation of neutral axis in degrees
na# - location of neutral axis as % rotated height in compression

variables c_outer_XY# (concrete outline) and c_rebar_XY# (rebar)
shall be defined prior to function call α# na# biaxial formulas used P# A Area Calculation Closing the polygon P#c P# P# 1 row stack eval : n# P#c rows eval : x#.i P#c i# 1 el eval : x#.i1 P#c i# 1 + 1 el eval : y#.i P#c i# 2 el eval : y#.i1 P#c i# 1 + 2 el eval : 12 / - (y#.i y#.i1 + (x#.i1 x#.i - (* i# 1 n# 1 - sum * eval 91 line : P# M.x 1st Moment of Area Closing the polygon P#c P# P# 1 row stack eval : n# P#c rows eval : x#.i P#c i# 1 el eval : x#.i1 P#c i# 1 + 1 el eval : y#.i P#c i# 2 el eval : y#.i1 P#c i# 1 + 2 el eval : 16 / - (x#.i1 x#.i - (y#.i 2^y#.i y#.i1 * + y#.i1 2^+ (* i# 1 n# 1 - sum * eval 91 line : P# M.y 1st Moment of Area Closing the polygon P#c P# P# 1 row stack eval : n# P#c rows eval : x#.i P#c i# 1 el eval : x#.i1 P#c i# 1 + 1 el eval : y#.i P#c i# 2 el eval : y#.i1 P#c i# 1 + 2 el eval : 16 / (y#.i1 y#.i - (x#.i 2^x#.i x#.i1 * + x#.i1 2^+ (* i# 1 n# 1 - sum * eval 91 line : θ Rot θ cos θ sin θ sin - θ cos 22 mat : 41 line check correctness of input na# na# 1000 random 100000 / Max : Locate Neutral Axis:Locate Neutral Axis rotate the section: c_outer_XY'# c_outer_XY# α# Rot * eval : rotate the rebar: c_rebar_temp# c_rebar_XY# 1 c_rebar_XY# rows range 12 range el eval : c_rebar_temp# c_rebar_temp# α# Rot * eval : c_rebar_XY'# c_rebar_temp# c_rebar_XY# 1 c_rebar_XY# rows range 3 c_rebar_XY# cols range el augment eval : locate Y' coordinate of NA: Y'.MIN# c_outer_XY'# 1 c_outer_XY'# rows range 2 el Min eval : Y'.MAX# c_outer_XY'# 1 c_outer_XY'# rows range 2 el Max eval : Y'.NA# Y'.MAX# na# Y'.MAX# Y'.MIN# - (* - eval : locate Y' coordinate of comp. extreme fiber: Y'.ext# Y'.MAX# eval : locate Y' coordinate of comp. block start: Y'.a# Y'.MAX# β.1 Y'.MAX# Y'.NA# - (* - eval : define rectangular stress block portion of the concrete section: #pgon c_outer_XY'# eval : #TLV Y'.a# eval : r #pgon rows eval : j 1 : k 1 r 1 - (range #pgon k 2 el #TLV ≥ #chunk j 1 el #pgon k 1 el eval : #chunk j 2 el #pgon k 2 el eval : j j 1 + : 31 line 0 if #pgon k 2 el #TLV - (#pgon k 1 + 2 el #TLV - (* 0 < #chunk j 1 el #pgon k 1 + 1 el #pgon k 1 el - (#TLV #pgon k 2 el - (* #pgon k 1 + 2 el #pgon k 2 el - / #pgon k 1 el + eval : #chunk j 2 el #TLV eval : j j 1 + : 31 line 0 if 21 line for c_outer_comp_XY'# #chunk : Calculate stress block properties (XY coordinates): c_outer_comp_XY# c_outer_comp_XY'# α# - Rot * eval : A.c_outer_comp# c_outer_comp_XY# A eval : xbar.c_outer_comp# c_outer_comp_XY# M.y A.c_outer_comp# / eval : ybar.c_outer_comp# c_outer_comp_XY# M.x A.c_outer_comp# / eval : summary of rebar force and moment: c_concr_ΣP# A.c_outer_comp# α.1 * φ.c * f'.c * eval : c_concr_ΣMx# c_concr_ΣP# ybar.c_outer_comp# * eval : c_concr_ΣMy# c_concr_ΣP# xbar.c_outer_comp# * eval : Calculate rebar properties: rebar coordinates: c_rebar_X# c_rebar_XY# 1 col eval : c_rebar_Y# c_rebar_XY# 2 col eval : rebar strain and stress: c_rebar_Y'# c_rebar_XY'# 2 col eval : c_rebar_ε c_rebar_Y'# Y'.NA# - Y'.ext# Y'.NA# - / ε.cu * eval : c_rebar_σ# φ.s f.y * c_rebar_ε sign * eval c_rebar_ε abs ε.sy ≥ φ.s c_rebar_ε * E.s * eval cases vectorize eval : rebar force: c_rebar_A# c_rebar_XY# 4 col eval : c_rebar_A# c_rebar_A# 1 n.plot /^mm 2^* vectorize eval : c_rebar_P# c_rebar_A# c_rebar_σ# * vectorize eval : rebar moment: c_rebar_Mx# c_rebar_P# c_rebar_Y# * vectorize eval : c_rebar_My# c_rebar_P# c_rebar_X# * vectorize eval : summary of rebar force and moment: c_rebar_ΣMx# c_rebar_Mx# sum eval : c_rebar_ΣMy# c_rebar_My# sum eval : c_rebar_ΣP# c_rebar_P# sum eval : Design Resistance - Pr, Mrx and Mry: P.rz# c_concr_ΣP# c_rebar_ΣP# + eval - : M.rx# c_concr_ΣMx# c_rebar_ΣMx# + eval - : M.ry# c_concr_ΣMy# c_rebar_ΣMy# + eval - : Output in kN and m; 1st column is α#, 2nd column is na#, 3rd is Pz, 4th is Mx, 5th is My α# ° / na# % / P.rz# kN / M.rx# kN m * / M.ry# kN m * / 15 mat eval 601 line :

function inputs (hidden):

]]>

Material Properties & Cover

]]>

Concrete:

]]>

f'.c 35 MPa * :

γ.c 2350 :

γ.c 9.814 1000 / * kN m 3^/ * kN m 3^/ 23.063

w.c 24 kN m 3^/ * :

E.c 3000 f'.c MPa / sqrt * 6900 + (γ.c 2300 / (1.5^* MPa * : GPa 25.4563

f.cr 0.4 f'.c MPa * sqrt * 3.2 MPa * Min : MPa 2.366

α.1 0.85 0.0015 f'.c * 1 MPa / * - : 0.798

β.1 0.97 0.0025 f'.c * 1 MPa / * - : 0.883

α.c_ΔT 10106 -^* :

coefficient of thermal expansion

Concrete cover to reinforcing steel unless noted otherwise:

]]>

Cover for concrete cast in contact with soil:

]]>

cover.0 60 mm * :

cover.1 100 mm * :

Concrete failure strain:

]]>

Steel yield strain:

]]>

ε.cu 0.0035 :

ε.sy 0.002 :

Steel:

]]>

f.y 400 MPa * :

E.s 200 GPa * :

f.pu 1860 MPa * :

E.p 200 GPa * :

A.s15 200 mm 2^* :

A.s20 300 mm 2^* :

A.s25 500 mm 2^* :

A.s30 700 mm 2^* :

A.s35 1000 mm 2^* :

A.s45 1500 mm 2^* :

A.s55 2500 mm 2^* :

d.b15 16 mm * :

d.b20 19.5 mm * :

d.b25 25.2 mm * :

d.b30 29.9 mm * :

d.b35 35.7 mm * :

d.b45 43.7 mm * :

d.b55 56.4 mm * :

Resistance Factors:

]]>

φ.c 0.75 :

φ.s 0.9 :

matrix [x1, x2 | y1, y2] - start and end coordinates
no - number of points
format - plotting type for 2D plot, say "."
size - font size of a point
color - color, say "red" xy12# no# format# size# color# pts.line pts# no# 5 matrix : x1# x2# y1# y2# 22 mat xy12# : dx# x2# x1# - no# 1 - / : dy# y2# y1# - no# 1 - / : xi# x1# : yi# y1# : i 1 no# range pts# i 1 el xi# : pts# i 2 el yi# : pts# i 3 el format# : pts# i 4 el size# : pts# i 5 el color# : xi# xi# dx# + : yi# yi# dy# + : 71 line for pts# 81 line :

matrice di rotazione; rotazione rispetto all'origine degli assi cartesiani.
 θ: angolo di rotazione antiorario. rotation matrix; rotatation about the origin of the Cartesian coordinate system.
 θ: anticlockwise rotation. θ Rot θ cos θ sin θ sin - θ cos 22 mat :

rettangolo.
 B: larghezza;
 H: altezza;
 θ: angolo di rotazione antiorario. rectangle.
 B: width;
 H: height;
 θ: anticlockwise rotation. B H θ rectangle 00 B 0 B H 0 H 0052 mat θ Rot * :

trasla un poligono.
 P: matrice del poligono;
 u.x: spostamento in x;
 u.y: spostamento in y. geometric translation of a polygon.
 P: polygon matrix;
 u.x: x translation;
 u.y: y translation. P u.x u.y Translate use homogeneous coordinates #P P 1 P rows 1 matrix + (augment : #P 100010 u.x u.y 133 mat * 1 P rows 12 submatrix 31 line :

ruota un poligono rispetto all'origine degli assi cartesiani.
 P: matrice del poligono;
 θ: angolo di rotazione antiorario. rotating a polygon about the origin of the Cartesian coordinate system.
 P: polygon matrix;
 θ: anticlockwise rotation. P θ Rotate use homogeneous coordinates #P P 1 P rows 1 matrix + (augment : #P θ cos θ sin - 0 θ sin θ cos 000133 mat * 1 P rows 12 submatrix 31 line :

Sign convention:

]]>

Cross-section of column is a z-plane. Positive x and y are to the right and upward. Positive moments cause compression at positive x locations and positive y locations. Compression is positive

]]>

Define section geometry assuming section centroid at the zero coordinates:

]]>

a# 2000 mm * :

b# 1350 mm * :

c# 80 mm * :

n.plot 0.4 :

(power of size scaling for plotting rebar)

]]>

c_outer_XY# a# b# 0 rectangle a# 2 / - b# 2 / - Translate 11 sys eval : m - 0.675 m * - m 0.675 m * - m 0.675 m * m - 0.675 m * m - 0.675 m * - 52 mat

c_rebar_XY# a# 2 / c# - (- a# 2 / c# - (b# 2 / c# - (b# 2 / c# - (22 mat 2 . A.s35 mm 2^/ (n.plot^red pts.line a# 2 / c# - (- a# 2 / c# - (b# 2 / c# - (- b# 2 / c# - (- 22 mat 2 . A.s35 mm 2^/ (n.plot^red pts.line 21 sys eval :

c_rebar_XY# c_rebar_XY# sys2mat eval :

n.rebar_entries# c_rebar_XY# rows eval :

i## 1 n.rebar_entries# range i## 1 ≡ c_rebar_temp# c_rebar_XY# i## el eval : c_rebar_temp# c_rebar_temp# c_rebar_XY# i## el stack eval : if for

c_rebar_XY# c_rebar_temp# eval :

0 ° * 7.3 % * biaxial 0 7.3 3158.51 - 2869.08 - 015 mat

test function call:

]]>

PzMxMy NA rotation % height in compr Pz Mx My degree % kN kN.m kN.m 25 mat :

PzMxMy PzMxMy 10 ° * 5 % * biaxial stack (:

PzMxMy NA rotation % height in compr Pz Mx My degree % kN kN.m kN.m 105 811.34 392.47 - 557.59 - 35 mat

Now that we have a function that produces input loads based on crossection responce, we will use the function to solve for crossection responce based on input loads. For example lets solve for "% height in compression" na#, given known neutral axis orientation (10 degrees):

]]>

Equation solver using a for loop:

rightSide# 3000 - :

guess1# 20 % * :

guess2# 1005 % * :

precision# 7 :

xi guess1# eval :

xii guess2# eval :

fi 10 ° * xi biaxial 3 el rightSide# - eval : 1856.7586 -

fii 10 ° * xii biaxial 3 el rightSide# - eval : 54962.8125 -

df fii fi - xii xi - / eval fi UnitsOf xi UnitsOf / try : 5391.4776 -

precision# 7 :

i 110000 range xii xi fi df / - eval : fii 10 ° * xii biaxial 3 el rightSide# - eval : df fii fi - xii xi - / eval df try : fii rightSide# / abs 10 precision# -^< break if xi xii : fi fii : 61 line for

xii 0.1654

10 ° * xii biaxial 3 el 3000 -

success !!!

Lets automate solution using The Newton-Raphson Solver:

]]>

leftSide### rightSide### variable### guess1### guess2### precision### solveNR solveNRcode### :

solveNRcode### xi### guess1### eval : xii### guess2### eval : value eval\0040\ leftSide### rightSide### - variable### value ≡ at eval variable### value : leftSide### rightSide### - eval 21 line try 11 line : fi### xi### eval\0040\ : fii### xii### eval\0040\ : df### fii### fi### - xii### xi### - / eval fi### UnitsOf xi### UnitsOf / try : i 110000 range xii### xi### fi### df### / - eval : fii### xii### eval\0040\ : df### fii### fi### - xii### xi### - / eval df try : fii### rightSide### variable### xii### ≡ at / abs 10 precision### -^< break if xi### xii### : fi### fii### : 61 line for xii### 81 line :

testing solver:

]]>

h 3^5 + h 2^- h + 50 h 01 - 7 solveNR 3.8207922517468

h 3^5 + h 2^- h + h 3.8207922517468 ≡ at 50

solver for biaxial(alpha, na)

]]>

alpha# rightSide### guess1### guess2### precision### solvePcr leftSide### alpha# variable### biaxial 3 el : solveNRcode### 21 line :

trying to solve:

]]>

α 10 ° * :

na α 3000 - 0.01 0.1 10 solvePcr : 0.16539

α na biaxial 10 16.539 3000 - 2631.388 - 2556.91 - 15 mat

u 4500 - 5004000 - range : j 1 u rows range na j el α u j el 0.01 0.1 10 solvePcr eval : for 21 line

<= division by zero ?

]]>

Project j 1 na rows range sol j el α na j el biaxial eval : for Unest the stack v sol 1 el : j 2 sol rows range v v sol j el stack : for v 51 line :

na 0.1937 0.1846 0.1752 0.1654 0.155 0.1425 0.1287 0.1144 0.0986 91 mat

...

]]>

Project 10 19.4 4500 - 3448.7 - 2749.1 - 10 18.5 4000 - 3179.8 - 2696 - 10 17.5 3500 - 2907.4 - 2631.4 - 10 16.5 3000 - 2631.4 - 2556.9 - 10 15.5 2500 - 2351.9 - 2464.7 - 10 14.2 2000 - 2070.3 - 2281.7 - 10 12.9 1500 - 1783.8 - 2070.6 - 10 11.4 1000 - 1491.6 - 1815 - 10 9.9 500 - 1194.2 - 1524.1 - 95 mat