9 files changed, 512 insertions, 0 deletions
diff --git a/3683/CH18/EX18.1/Ex18_1.sce b/3683/CH18/EX18.1/Ex18_1.sce
new file mode 100644
index 000000000..e64a26f30
--- /dev/null
+++ b/3683/CH18/EX18.1/Ex18_1.sce
@@ -0,0 +1,12 @@
+Pu=3000//in kN
+fck=20//in MPa
+fy=415//in MPa
+l=3//unsupported length, in m
+//assume 1% steel
+Ag=Pu*10^3/(0.4*fck*0.99+0.67*fy*0.01)//in sq mm
+L=sqrt(Ag)//assuming a square column
+L=530//in mm
+Asc=0.01*L^2//in sq mm
+emin=l*10^3/500+L/30//in mm
+ep=0.05*L//>emin, hence OK
+mprintf("Column size - %d x %d mm",L,L)
diff --git a/3683/CH18/EX18.2/Ex18_2.sce b/3683/CH18/EX18.2/Ex18_2.sce
new file mode 100644
index 000000000..3e3856ed2
--- /dev/null
+++ b/3683/CH18/EX18.2/Ex18_2.sce
@@ -0,0 +1,15 @@
+Pu=1500//in kN
+fck=15//in MPa
+fy=250//in MPa
+l=2.75//unsupported length, in m
+//assume 1% steel
+Ag=Pu*10^3/(0.4*fck*0.99+0.67*fy*0.01)//in sq mm
+L1=225//assuming a square column
+L2=Ag/L1//in mm
+L2=880//in mm
+Asc=0.01*L1*L2//in sq mm
+e1=l*10^3/500+L1/30//in mm
+e2=l*10^3/500+L2/30//in mm
+ep1=0.05*L1//<e1
+ep2=0.05*L2//>e2, hence Ok
+mprintf("The column is safe on long dimension side but not on short dimension side. As such, the column be checked for eccentricity in short direction.")
diff --git a/3683/CH18/EX18.3/Ex18_3.sce b/3683/CH18/EX18.3/Ex18_3.sce
new file mode 100644
index 000000000..3f73c91ac
--- /dev/null
+++ b/3683/CH18/EX18.3/Ex18_3.sce
@@ -0,0 +1,35 @@
+b=225//in mm
+D=500//in mm
+c=45//cover, in mm
+Asc=2463//in sq mm
+Ast=Asc
+fck=15//in MPa
+fy=250//in MPa
+fcc=0.446*fck//in MPa
+//(i)
+xu=1.1*D//in mm
+m=0.43*D//in mm
+esc1=0.002*(xu-c)/(xu-m)
+esc2=0.002*(xu-D+c)/(xu-m)
+//by interpolation
+fsc1=217.5//in MPa
+fsc2=217.5*esc2/0.0010875//in MPa
+//stress block parameters for xu / D = 1.1
+n=0.384
+l=0.443
+A=n*fck*D//area of stress block
+r=l*D//distance of c.g., in mm
+Pu=(A*b+Asc*(fsc1-fcc)+Ast*fsc2)/10^3
+Mu=(A*b*(D/2-r)+Asc*(fsc1-fcc)*(D/2-c)-Ast*fsc2*(D/2-c))/10^6
+mprintf("(i) For xu = 1.1 D\nP=%f kN\nMu=%f kN-m\n",Pu,Mu)
+//answer in textbook is incorrect
+//(ii)
+xu=330//in mm
+esc=0.0035*(xu-c)/xu
+est=0.0035*(D-c-xu)/xu
+//by interpolation
+fsc=217.5//in MPa
+fst=217.5//in MPa
+Pu=(0.36*fck*b*xu+Asc*(fsc-fcc)-Ast*fst)/10^3//in kN
+Mu=(0.36*fck*b*xu*(D/2-0.416*xu)+Asc*(fsc-fcc)*(D/2-c)+Ast*fst*(D/2-c))/10^6//in kN-m
+mprintf("(ii) For xu = 330 mm\nP=%f kN\nMu=%f kN-m",Pu,Mu)
diff --git a/3683/CH18/EX18.4/Ex18_4.sce b/3683/CH18/EX18.4/Ex18_4.sce
new file mode 100644
index 000000000..b2960c1f4
--- /dev/null
+++ b/3683/CH18/EX18.4/Ex18_4.sce
@@ -0,0 +1,33 @@
+b=300//in mm
+D=400//in mm
+c=30//cover, in mm
+Asc=452//in sq mm
+Ast=Asc
+fck=15//in MPa
+fy=415//in MPa
+fcc=0.446*fck//in MPa
+//(i)
+xu=1.4*D//in mm
+m=0.43*D//in mm
+esc1=0.002*(xu-c)/(xu-m)
+esc2=0.002*(xu-D+c)/(xu-m)
+//by interpolation
+fsc1=356.8//in MPa
+fsc2=238.68//in MPa
+//stress block parameters for xu / D = 1.4
+n=0.417
+l=0.475
+A=n*fck*D//area of stress block
+r=l*D//distance of c.g., in mm
+Pu=(A*b+Asc*(fsc1-fcc)+Ast*fsc2)/10^3//in kN
+Mu=(A*b*(D/2-r)+Asc*(fsc1-fcc)*(D/2-c)-Ast*fsc2*(D/2-c))/10^6//in kN-m
+mprintf("(i) For xu = 1.4 D\nP=%f kN\nMu=%f kN-m\n",Pu,Mu)
+//(ii)
+xu=370//in mm
+esc=0.0035*(xu-c)/xu
+est=0.0035*(D-c-xu)/xu
+//by interpolation
+fsc=355.8//in MPa
+Pu=(0.36*fck*b*xu+Asc*(fsc-fcc))/10^3//in kN
+Mu=(0.36*fck*b*xu*(D/2-0.416*xu)+Asc*(fsc-fcc)*(D/2-c))/10^6//in kN-m
+mprintf("(ii) For xu = 370 mm\nP=%f kN\nMu=%f kN-m",Pu,Mu)
diff --git a/3683/CH18/EX18.5/Ex18_5.sce b/3683/CH18/EX18.5/Ex18_5.sce
new file mode 100644
index 000000000..b6c66e42d
--- /dev/null
+++ b/3683/CH18/EX18.5/Ex18_5.sce
@@ -0,0 +1,35 @@
+b=225//in mm
+D=500//in mm
+c=50//cover, in mm
+Asc=1520//in sq mm
+Ast=Asc
+fck=20//in MPa
+fy=500//in MPa
+fcc=0.446*fck//in MPa
+//(i)
+xu=1.3*D//in mm
+m=0.43*D//in mm
+esc1=0.002*(xu-c)/(xu-m)
+esc2=0.002*(xu-D+c)/(xu-m)
+//by interpolation
+fsc1=412.515//in MPa
+fsc2=183.794//in MPa
+//stress block parameters for xu / D = 1.3
+n=0.409
+l=0.468
+A=n*fck*D//area of stress block
+r=l*D//distance of c.g., in mm
+Pu=(A*b+Asc*(fsc1-fcc)+Ast*fsc2)/10^3//in kN
+Mu=(A*b*(D/2-r)+Asc*(fsc1-fcc)*(D/2-c)-Ast*fsc2*(D/2-c))/10^6//in kN-m
+mprintf("(i) For xu = 1.3 D\nP=%f kN\nMu=%f kN-m\n",Pu,Mu)
+//(ii)
+xu=400//in mm
+esc=0.0035*(xu-c)/xu
+est=0.0035*(D-c-xu)/xu
+//by interpolation
+fsc=422.11//in MPa
+fst=87.45//in MPa
+Pu=(0.36*fck*b*xu+Asc*(fsc-fcc)-Ast*fst)/10^3//in kN
+Mu=(0.36*fck*b*xu*(D/2-0.416*xu)+Asc*(fsc-fcc)*(D/2-c)+Ast*fst*(D/2-c))/10^6//in kN-m
+mprintf("(ii) For xu = 400 mm\nP=%f kN\nMu=%f kN-m",Pu,Mu)
+//answer in textbook for Mu in (ii) is incorrect
diff --git a/3683/CH18/EX18.6/Ex18_6.jpeg b/3683/CH18/EX18.6/Ex18_6.jpeg
new file mode 100644
index 000000000..01971b3cb
--- /dev/null
+++ b/3683/CH18/EX18.6/Ex18_6.jpeg
diff --git a/3683/CH18/EX18.6/Ex18_6.sce b/3683/CH18/EX18.6/Ex18_6.sce
new file mode 100644
index 000000000..cf18fb256
--- /dev/null
+++ b/3683/CH18/EX18.6/Ex18_6.sce
@@ -0,0 +1,328 @@
+b=250//width, in mm
+D=450//depth, in mm
+c=50//cover, in mm
+Asc=1472//in sq mm
+Ast=Asc
+fck=15//in MPa
+fcc=0.446*fck//in MPa
+fy=250//in MPa
+Es=2*10^5//in MPa
+ey=0.87*fy/Es//strain in mild steel at yield point
+fs=0.87*fy//stress in mild steel at yield point, in MPa
+
+//xu=infinity
+Pu1=(0.446*fck*(b*D-Asc-Ast)+(Asc+Ast)*fs)/10^3//in kN
+Mu1=0//in kN-m
+
+//xu=1.5 D
+xu=1.5*D//in mm
+m=0.43*D//in mm
+esc1=0.002*(xu-c)/(xu-m)
+if(esc1<=ey)
+    fsc1=esc1/ey*fs
+else
+    fsc1=fs
+end
+esc2=0.002*(xu-D+c)/(xu-m)//>ey
+if(esc2<=ey)
+    fsc2=esc2/ey*fs
+else
+    fsc2=fs
+end
+//stress block parameters for xu / D = 1.5
+n=0.422
+l=0.48
+A=n*fck*D//area of stress block
+r=l*D//distance of c.g.
+Pu2=(A*b+Asc*fsc1+Ast*fsc2)/10^3//in kN
+Mu2=(A*b*(D/2-r)+Asc*fsc1*(D/2-c)-Ast*fsc2*(D/2-c))/10^6//in kN-m
+
+//xu=1.3 D
+xu=1.3*D//in mm
+m=0.43*D//in mm
+esc1=0.002*(xu-c)/(xu-m)
+if(esc1<=ey)
+    fsc1=esc1/ey*fs
+else
+    fsc1=fs
+end
+esc2=0.002*(xu-D+c)/(xu-m)//>ey
+if(esc2<=ey)
+    fsc2=esc2/ey*fs
+else
+    fsc2=fs
+end
+//stress block parameters for xu / D = 1.3
+n=0.409
+l=0.468
+A=n*fck*D//area of stress block
+r=l*D//distance of c.g.
+Pu3=(A*b+Asc*fsc1+Ast*fsc2)/10^3//in kN
+Mu3=(A*b*(D/2-r)+Asc*fsc1*(D/2-c)-Ast*fsc2*(D/2-c))/10^6//in kN-m
+
+//xu=1.2 D
+xu=1.2*D//in mm
+m=0.43*D//in mm
+esc1=0.002*(xu-c)/(xu-m)
+if(esc1<=ey)
+    fsc1=esc1/ey*fs
+else
+    fsc1=fs
+end
+esc2=0.002*(xu-D+c)/(xu-m)//>ey
+if(esc2<=ey)
+    fsc2=esc2/ey*fs
+else
+    fsc2=fs
+end
+//stress block parameters for xu / D = 1.2
+n=0.399
+l=0.458
+A=n*fck*D//area of stress block
+r=l*D//distance of c.g.
+Pu4=(A*b+Asc*fsc1+Ast*fsc2)/10^3//in kN
+Mu4=(A*b*(D/2-r)+Asc*fsc1*(D/2-c)-Ast*fsc2*(D/2-c))/10^6//in kN-m
+
+//xu=1.1 D
+xu=1.1*D//in mm
+m=0.43*D//in mm
+esc1=0.002*(xu-c)/(xu-m)
+if(esc1<=ey)
+    fsc1=esc1/ey*fs
+else
+    fsc1=fs
+end
+esc2=0.002*(xu-D+c)/(xu-m)//>ey
+if(esc2<=ey)
+    fsc2=esc2/ey*fs
+else
+    fsc2=fs
+end
+//stress block parameters for xu / D = 1.1
+n=0.384
+l=0.443
+A=n*fck*D//area of stress block
+r=l*D//distance of c.g.
+Pu5=(A*b+Asc*fsc1+Ast*fsc2)/10^3//in kN
+Mu5=(A*b*(D/2-r)+Asc*fsc1*(D/2-c)-Ast*fsc2*(D/2-c))/10^6//in kN-m
+
+//xu = D
+xu=D//in mm
+m=0.43*D//in mm
+esc1=0.002*(xu-c)/(xu-m)
+if(esc1<=ey)
+    fsc1=esc1/ey*fs
+else
+    fsc1=fs
+end
+esc2=0.002*(xu-D+c)/(xu-m)//>ey
+if(esc2<=ey)
+    fsc2=esc2/ey*fs
+else
+    fsc2=fs
+end
+//stress block parameters for xu / D = 1
+n=0.361
+l=0.416
+A=n*fck*D//area of stress block
+r=l*D//distance of c.g.
+Pu6=(A*b+Asc*fsc1+Ast*fsc2)/10^3//in kN
+Mu6=(A*b*(D/2-r)+Asc*fsc1*(D/2-c)-Ast*fsc2*(D/2-c))/10^6//in kN-m
+
+//xu=400 mm
+xu=400//in mm
+esc=0.0035*(xu-c)/xu
+if(esc<=ey)
+    fsc=esc/ey*fs
+else
+    fsc=fs
+end
+est=0.0035*(D-xu-c)/xu
+if(est<=ey)
+    fst=est/ey*fs
+else
+    fst=fs
+end
+Pu7=(0.36*fck*b*xu+Asc*fsc-Ast*fst)/10^3
+Mu7=(0.36*fck*b*xu*(D/2-0.416*xu)+Asc*fsc*(D/2-c)+Ast*fst*(D/2-c))/10^6
+
+//xu=375 mm
+xu=375//in mm
+esc=0.0035*(xu-c)/xu
+if(esc<=ey)
+    fsc=esc/ey*fs
+else
+    fsc=fs
+end
+est=0.0035*(D-xu-c)/xu
+if(est<=ey)
+    fst=est/ey*fs
+else
+    fst=fs
+end
+Pu8=(0.36*fck*b*xu+Asc*fsc-Ast*fst)/10^3
+Mu8=(0.36*fck*b*xu*(D/2-0.416*xu)+Asc*fsc*(D/2-c)+Ast*fst*(D/2-c))/10^6
+
+//xu=350 mm
+xu=350//in mm
+esc=0.0035*(xu-c)/xu
+if(esc<=ey)
+    fsc=esc/ey*fs
+else
+    fsc=fs
+end
+est=0.0035*(D-xu-c)/xu
+if(est<=ey)
+    fst=est/ey*fs
+else
+    fst=fs
+end
+Pu9=(0.36*fck*b*xu+Asc*fsc-Ast*fst)/10^3
+Mu9=(0.36*fck*b*xu*(D/2-0.416*xu)+Asc*fsc*(D/2-c)+Ast*fst*(D/2-c))/10^6
+
+//xu=325 mm
+xu=325//in mm
+esc=0.0035*(xu-c)/xu
+if(esc<=ey)
+    fsc=esc/ey*fs
+else
+    fsc=fs
+end
+est=0.0035*(D-xu-c)/xu
+if(est<=ey)
+    fst=est/ey*fs
+else
+    fst=fs
+end
+Pu10=(0.36*fck*b*xu+Asc*fsc-Ast*fst)/10^3
+Mu10=(0.36*fck*b*xu*(D/2-0.416*xu)+Asc*fsc*(D/2-c)+Ast*fst*(D/2-c))/10^6
+
+//xu=300 mm
+xu=300//in mm
+esc=0.0035*(xu-c)/xu
+if(esc<=ey)
+    fsc=esc/ey*fs
+else
+    fsc=fs
+end
+est=0.0035*(D-xu-c)/xu
+if(est<=ey)
+    fst=est/ey*fs
+else
+    fst=fs
+end
+Pu11=(0.36*fck*b*xu+Asc*fsc-Ast*fst)/10^3
+Mu11=(0.36*fck*b*xu*(D/2-0.416*xu)+Asc*fsc*(D/2-c)+Ast*fst*(D/2-c))/10^6
+
+//xu=275 mm
+xu=275//in mm
+esc=0.0035*(xu-c)/xu
+if(esc<=ey)
+    fsc=esc/ey*fs
+else
+    fsc=fs
+end
+est=0.0035*(D-xu-c)/xu
+if(est<=ey)
+    fst=est/ey*fs
+else
+    fst=fs
+end
+Pu12=(0.36*fck*b*xu+Asc*fsc-Ast*fst)/10^3
+Mu12=(0.36*fck*b*xu*(D/2-0.416*xu)+Asc*fsc*(D/2-c)+Ast*fst*(D/2-c))/10^6
+
+//xu=250 mm
+xu=250//in mm
+esc=0.0035*(xu-c)/xu
+if(esc<=ey)
+    fsc=esc/ey*fs
+else
+    fsc=fs
+end
+est=0.0035*(D-xu-c)/xu
+if(est<=ey)
+    fst=est/ey*fs
+else
+    fst=fs
+end
+Pu13=(0.36*fck*b*xu+Asc*fsc-Ast*fst)/10^3
+Mu13=(0.36*fck*b*xu*(D/2-0.416*xu)+Asc*fsc*(D/2-c)+Ast*fst*(D/2-c))/10^6
+
+//xu=225 mm
+xu=225//in mm
+esc=0.0035*(xu-c)/xu
+if(esc<=ey)
+    fsc=esc/ey*fs
+else
+    fsc=fs
+end
+est=0.0035*(D-xu-c)/xu
+if(est<=ey)
+    fst=est/ey*fs
+else
+    fst=fs
+end
+Pu14=(0.36*fck*b*xu+Asc*fsc-Ast*fst)/10^3
+Mu14=(0.36*fck*b*xu*(D/2-0.416*xu)+Asc*fsc*(D/2-c)+Ast*fst*(D/2-c))/10^6
+
+//xu=200 mm
+xu=200//in mm
+esc=0.0035*(xu-c)/xu
+if(esc<=ey)
+    fsc=esc/ey*fs
+else
+    fsc=fs
+end
+est=0.0035*(D-xu-c)/xu
+if(est<=ey)
+    fst=est/ey*fs
+else
+    fst=fs
+end
+Pu15=(0.36*fck*b*xu+Asc*fsc-Ast*fst)/10^3
+Mu15=(0.36*fck*b*xu*(D/2-0.416*xu)+Asc*fsc*(D/2-c)+Ast*fst*(D/2-c))/10^6
+
+//xu=150 mm
+xu=150//in mm
+esc=0.0035*(xu-c)/xu
+if(esc<=ey)
+    fsc=esc/ey*fs
+else
+    fsc=fs
+end
+est=0.0035*(D-xu-c)/xu
+if(est<=ey)
+    fst=est/ey*fs
+else
+    fst=fs
+end
+Pu16=(0.36*fck*b*xu+Asc*fsc-Ast*fst)/10^3
+Mu16=(0.36*fck*b*xu*(D/2-0.416*xu)+Asc*fsc*(D/2-c)+Ast*fst*(D/2-c))/10^6
+
+//xu=100 mm
+xu=100//in mm
+esc=0.0035*(xu-c)/xu
+if(esc<=ey)
+    fsc=esc/ey*fs
+else
+    fsc=fs
+end
+est=0.0035*(D-xu-c)/xu
+if(est<=ey)
+    fst=est/ey*fs
+else
+    fst=fs
+end
+Pu17=(0.36*fck*b*xu+Asc*fsc-Ast*fst)/10^3
+Mu17=(0.36*fck*b*xu*(D/2-0.416*xu)+Asc*fsc*(D/2-c)+Ast*fst*(D/2-c))/10^6
+
+//xu=0.531 d
+d=D-c
+xu=0.531*d
+Pu18=0//in kN
+Mu18=0.149*fck*b*d^2/10^6//in kN-m
+
+Pu=[Pu1 Pu2 Pu3 Pu4 Pu5 Pu6 Pu7 Pu8 Pu9 Pu10 Pu11 Pu12 Pu13 Pu14 Pu15 Pu16 Pu17 Pu18]
+Mu=[Mu1 Mu2 Mu3 Mu4 Mu5 Mu6 Mu7 Mu8 Mu9 Mu10 Mu11 Mu12 Mu13 Mu14 Mu15 Mu16 Mu17 Mu18]
+xtitle('P-M Interaction Diagram', 'Mu (kN-m)', 'Pu (kN)')
+plot(Mu,Pu)
diff --git a/3683/CH18/EX18.8/Ex18_8.sce b/3683/CH18/EX18.8/Ex18_8.sce
new file mode 100644
index 000000000..1fa76e868
--- /dev/null
+++ b/3683/CH18/EX18.8/Ex18_8.sce
@@ -0,0 +1,31 @@
+b=250//column width in mm
+D=450//column depth in mm
+Asc=2*1472//in sq mm
+fck=15//in MPa
+fy=250//in MPa
+ex=200//in mm
+ey=150//in mm
+//from interaction curve
+//for ex=200 mm on x-axis
+Pum1=610//in kN
+Muy1=120//in kN-m
+//for ey=150 mm on y-axis
+Pum2=720//in kN
+Mux1=106//in kN-m
+//(i)
+Pu=300//in kN
+Mux=Pu*ey/10^3//in kN-m
+Muy=Pu*ex/10^3//in kN-m
+Puz=(0.45*fck*(b*D-Asc)+0.75*fy*Asc)/10^3//in kN
+a=Pu/Puz
+an=1+1/0.6*(a-0.2)
+b=(Mux/Mux1)^an+(Muy/Muy1)^an//<1
+mprintf("The column can take a load of 300 kN with ex=200 mm and ey=150 mm\n")
+//(ii)
+Pu=500//in kN
+Mux=Pu*ey/10^3//in kN-m
+Muy=Pu*ex/10^3//in kN-m
+a=Pu/Puz
+an=1+1/0.6*(a-0.2)
+b=(Mux/Mux1)^an+(Muy/Muy1)^an//>1
+mprintf("The section is not suitable for a load of 500 kN with ex=200 mm and ey=150 mm\n")
diff --git a/3683/CH18/EX18.9/Ex18_9.sce b/3683/CH18/EX18.9/Ex18_9.sce
new file mode 100644
index 000000000..c8521f0e3
--- /dev/null
+++ b/3683/CH18/EX18.9/Ex18_9.sce
@@ -0,0 +1,23 @@
+b=250//column width, in mm
+D=500//column depth, in mm
+lex=4//in m
+ley=4//in m
+Pu=300//in kN
+Asc=1472//in sq mm
+Ast=1472//in sq mm
+fck=15//in MPa
+fy=250//in MPa
+c=50//cover, in mm
+Max=Pu*10^3*D/2000*(lex/(D/10^3))^2/10^6//in kN-m
+May=Pu*10^3*b/2000*(ley/(b/10^3))^2/10^6//in kN-m
+Puz=(0.45*fck*(b*D-(Asc+Ast))+0.75*fy*(Asc+Ast))/10^3//in kN
+//to find Pb
+xu=(D-c)/(1+0.002/0.0035)//in mm
+fsc=217.5//in MPa
+fst=217.5//in MPa
+Pb=(0.36*fck*b*xu+fsc*Asc-fst*Ast)/10^3//in kN
+k=(Puz-Pu)/(Puz-Pb)//>1
+k=1
+Max=k*Max//in kN-m
+May=k*May//in kN-m
+mprintf("Additional Moments are:\nMax=%f kN/m\nMay=%f kN-m", Max,May)