initial commit / add all books

15 files changed, 304 insertions, 0 deletions

diff --git a/2732/CH7/EX7.10/Ex7_10.sce b/2732/CH7/EX7.10/Ex7_10.sce
new file mode 100755
index 000000000..8924deaba
--- /dev/null
+++ b/2732/CH7/EX7.10/Ex7_10.sce
@@ -0,0 +1,21 @@
+clc

+//initialization of variables

+clear

+Ys=17000 //kg/cm^2

+E=2*10^6 //kg/cm^2

+d1=1 //mm

+d=1 //cm

+//calculations: 1 cm

+R=E*d/(2*Ys)

+M=Ys*%pi*d^3/32

+// results

+printf('%d cm daimeter wire:',d)

+printf('\n Minimum radius = %.2f cm',R)

+printf('\n Bending Moment = %.2f kg-cm',M)

+// calculations: 1 mm

+R1=R/(d1*10)

+M1=M/(d1*1000)

+// results

+printf('\n %d mm daimeter wire:',d1)

+printf('\n Minimum radius = %.2f cm',R1)

+printf('\n Bending Moment = %.2f kg-cm',M1)

diff --git a/2732/CH7/EX7.11/Ex7_11.sce b/2732/CH7/EX7.11/Ex7_11.sce
new file mode 100755
index 000000000..4833682f0
--- /dev/null
+++ b/2732/CH7/EX7.11/Ex7_11.sce
@@ -0,0 +1,25 @@
+clc

+//initialization of variables

+clear

+t=0.5 //cm

+s=2 //m

+p1=7.5 //cm

+p2=10 //cm

+d=p2/2

+str=1650 //kg/cm^2

+// calculations

+// ab

+IxX=p1*t^3/12+t*p1*d^2

+// bc

+alpha=atan(3/4)

+Ixx=t*(p1+d)^3/12

+Iyy=0

+Ixy=0

+Iuu=Ixx*cos(alpha)^2+Iyy*sin(alpha)^2-Ixy*sin(2*alpha)

+Ixx=Iuu+IxX

+IXX=Ixx*100/(2*p1)

+Z=IXX/(d+t/2)

+w=str*Z*8/(s^2*100)

+w=w/1000

+//Results

+printf('w = %.1f tonne/m',w)

diff --git a/2732/CH7/EX7.12/Ex7_12.sce b/2732/CH7/EX7.12/Ex7_12.sce
new file mode 100755
index 000000000..47a42392d
--- /dev/null
+++ b/2732/CH7/EX7.12/Ex7_12.sce
@@ -0,0 +1,45 @@
+clc

+//initialization of variables

+clear

+wb=10 //cm

+wh=20 //cm

+sb=0.5 //cm

+sh=10 //cm

+s=6 //m

+fs=1650 //kg/cm^2

+fw=150 //kg/cm^2

+Es=2*10^6 //kg/cm^2

+Ew=12*10^4 //kg/cm^2

+

+//calculations

+// Method 1

+as=2*fs/(21*Es)

+aw=2*fw/(20*Ew)

+a=min(as,aw)

+ss=a*Ew*wh/2

+//Moment resistance of steel portion

+F=(fs+1573)/2*sb*sh

+k=sb/3*(fs+2*1573)/(fs+1573)

+Ms=2*F*(10.5-k)

+//Moment resistance of wooden portion

+F=ss*wb*wb/2

+Mw=2*(F*(wb-wb/3))

+M=Ms+Mw

+//Total udl supported

+W=M*8/(s*100)

+

+//Results

+printf('Using method 1')

+printf('\n W = %d kg',W)

+

+//Method 2

+nE=Es/Ew

+nf=fs/fw

+Is=2*(0+sb*sh*10.25^2)

+Iw=0.6*wh^3/12

+I=Is+Iw

+W=fs*I*8/(s*100*10.5)

+

+//Results

+printf('\n Using method 2')

+printf('\n W = %d kg',W)

diff --git a/2732/CH7/EX7.13/Ex7_13.sce b/2732/CH7/EX7.13/Ex7_13.sce
new file mode 100755
index 000000000..2df7eecc2
--- /dev/null
+++ b/2732/CH7/EX7.13/Ex7_13.sce
@@ -0,0 +1,18 @@
+clc

+//initialization of variables

+clear

+p=6 //mm

+Ixx=2375 //cm^4

+Es=2*10^6 //kg/cm^2

+EAl=0.667*10^6 //kg/cm^2

+d1=10.6 //cm

+d2=10 //cm

+// calculations

+I1=2*(0+p/10*10*10.3^2)

+I2=Ixx*EAl/Es

+I=I1+I2

+n=I/I2

+// results

+printf('stiffness ratio = %.2f ',n)

+n1=Es*d1/(d2*EAl)

+printf('\n Stress ration = %.2f ',n1)

diff --git a/2732/CH7/EX7.14/Ex7_14.sce b/2732/CH7/EX7.14/Ex7_14.sce
new file mode 100755
index 000000000..4d326d688
--- /dev/null
+++ b/2732/CH7/EX7.14/Ex7_14.sce
@@ -0,0 +1,17 @@
+clc

+//initilization of new variables

+clear

+wt=0.8 //cm

+ft=1.4 //cm

+w=10 //cm

+y=20 //cm

+// Sigma_y: yield stress is not given explicitly

+k1=wt*(40-2*ft)/2

+Zp=(14*19.3+k1*9.3)*2

+If=2*(w*ft^3/12+w*ft*19.3^2)

+Iw=wt*(40-2*ft)^3/12

+I=Iw+If

+Z=I/y

+sf=Zp/Z

+//Results

+printf('shape factor = %.2f ',sf)

diff --git a/2732/CH7/EX7.15/Ex7_15.sce b/2732/CH7/EX7.15/Ex7_15.sce
new file mode 100755
index 000000000..6b4d6cc43
--- /dev/null
+++ b/2732/CH7/EX7.15/Ex7_15.sce
@@ -0,0 +1,14 @@
+clc

+//initilization of new variables

+clear

+wt=0.8 //cm

+ft=1.4 //cm

+w=10 //cm

+y=20 //cm

+T=750 //T==750*sigma_y

+// calculations

+MpF=ft*w*(40-2*ft)

+c1=((40-2*ft)/2)^2-(T-MpF)/wt

+c=sqrt(3*c1)

+// results

+printf('Elastic core of %.1f cm depth is present',2*c)

diff --git a/2732/CH7/EX7.17/Ex7_17.sce b/2732/CH7/EX7.17/Ex7_17.sce
new file mode 100755
index 000000000..fb52f8f3f
--- /dev/null
+++ b/2732/CH7/EX7.17/Ex7_17.sce
@@ -0,0 +1,24 @@
+clc

+//initialization of new variables

+clear

+P=2000 //kg

+a=4 //cm

+b=1 //cm

+d=7 //cm

+r=3 //cm

+// calculations

+A=(a+b)/2*d

+xbar=(a+b*2)*d/(r*(a+b))

+rbar=r+xbar

+I=b*d^3/12+r*d^3/12

+Ixx=I-A*2.8^2

+e=Ixx/(rbar*A)

+f1=P*5.8*(xbar-0.62)/(A*0.62*r)

+f2=P*5.8*(-d+2.18)/(A*0.62*(5.18+d-2.18))

+str=P/A

+Str_i=f1+str

+Str_o=-f2-str

+//Results

+printf('stress at the inner side of the hook = %.1f kg/cm^2 (tensile)',Str_i)

+printf('\n stress at the outer side of the hook = %.1f kg/cm^2 (compressive)',Str_o)

+// approximations involved in the text

diff --git a/2732/CH7/EX7.2/Ex7_2.sce b/2732/CH7/EX7.2/Ex7_2.sce
new file mode 100755
index 000000000..a974b4be4
--- /dev/null
+++ b/2732/CH7/EX7.2/Ex7_2.sce
@@ -0,0 +1,20 @@
+clc

+//initialization of variables

+clear

+s=3 //m

+n=60

+p=50 //kg

+// calculations

+W=n*p

+Rc=W*2/s

+Rb=W-Rc

+dx = 0.001;

+x = 0:dx:s

+n = s/dx +1;

+for i = 1:n

+    Sx(i) = -Rb + Rc*x(i)^2/6;

+    Mx(i) = Rb*x(i) - Rc*x(i)^3 /18;

+end

+//Results

+figure(1);plot(x,Sx);title("Shear force diagram");xlabel("X (in m)");ylabel("Shear force (in kg)");

+figure(2);plot(x,Mx);title("Bending Moment diagram");xlabel("X (in m)");ylabel("Bending Moment (in kg-m)");

diff --git a/2732/CH7/EX7.2/Ex7_2_a.jpeg b/2732/CH7/EX7.2/Ex7_2_a.jpeg
new file mode 100755
index 000000000..4c4c0274d
--- /dev/null
+++ b/2732/CH7/EX7.2/Ex7_2_a.jpeg
diff --git a/2732/CH7/EX7.2/Ex7_2_b.jpeg b/2732/CH7/EX7.2/Ex7_2_b.jpeg
new file mode 100755
index 000000000..6b3d19457
--- /dev/null
+++ b/2732/CH7/EX7.2/Ex7_2_b.jpeg
diff --git a/2732/CH7/EX7.20/Ex7_20.sce b/2732/CH7/EX7.20/Ex7_20.sce
new file mode 100755
index 000000000..08d951a0c
--- /dev/null
+++ b/2732/CH7/EX7.20/Ex7_20.sce
@@ -0,0 +1,42 @@
+clc

+//initialization of new variables

+clear

+t=1 //cm

+a=40 //cm

+A=236

+// calculations

+ybar=a*t*0.5+(50-1)*4*0.5/(a*t+(50-1)*4)

+y1bar=1.25*a-ybar

+IAA=a*t^3/3+(50-1)^3*4/12+(50-1)*4*25.5^2

+Io=IAA-A*ybar^2

+//part (1)

+r=y1bar/ybar

+// results

+printf('Ratio of maximum bending stress in the stem and flange')

+printf('\n Ratio = %.2f',r)

+//part(2)

+// calculations

+r=(2/3*388*29.56)-(2/3*160*20.44)-(228*20.44)

+r=r/(2*2/3*388*29.56)

+// results

+printf('\n Ratio of S.F in flange to total S.F')

+printf('\n Ratio = %.2f percent',r*100)

+// part (3)

+// calculations

+r=359*200/Io

+// results

+printf('\n Ratio of maximum shear stress in the flange to average sher stress in the stem')

+printf('\n Ratio = %.2f ',r)

+//part (4)

+// calculations

+s=10 //m

+r=r/0.922

+sigma=1650 //kg/cm^2

+shear=945 //kg/cm^2

+wsh=2*200*shear/(r*s)

+wsi=8*Io*sigma/(s^2*10*29.56)

+w=min(wsh,wsi)

+// results

+printf('\n Maximum u.d.l. = %d kg/m ',w)

+

+//wrong moment of Inertia (Io) in the text and hence part (3) and part (4) are wrong

diff --git a/2732/CH7/EX7.21/Ex7_21.sce b/2732/CH7/EX7.21/Ex7_21.sce
new file mode 100755
index 000000000..31e9abf0a
--- /dev/null
+++ b/2732/CH7/EX7.21/Ex7_21.sce
@@ -0,0 +1,17 @@
+clc

+//initialization of new variables

+clear

+a=30 //cm

+t=2.5 //cm

+S=15 //cm

+s=5 //Tonne

+// calculations

+I=a*a^3-25*25^3

+I=I/12

+tau_zx=s*1000*27.5*t*25/(4*35000*t)

+FA=S*t*tau_zx

+tau_xy=s*1000*a*t*27.5/(4*35000*t)

+FB=tau_xy*t*S

+//Results

+printf('case A \n F = %d kg',FA)

+printf('\n case B \n F= %d kg',FB)

diff --git a/2732/CH7/EX7.23/Ex7_23.sce b/2732/CH7/EX7.23/Ex7_23.sce
new file mode 100755
index 000000000..777aaf6cd
--- /dev/null
+++ b/2732/CH7/EX7.23/Ex7_23.sce
@@ -0,0 +1,12 @@
+clc

+//initialization of variables

+clear

+h=40 //cm

+b=10 //cm

+t1=1.4 //cm

+t2=0.8 //cm

+Ixx=13989.5 //cm^4

+//calculations

+e=b^2*h^2*t1/(4*Ixx)*(1-t1/h-t1/b+t1^2/(b*h))*(1-t1/h)

+//Results

+printf('Shear center: \n e = %.2f cm',e)

diff --git a/2732/CH7/EX7.33/Ex7_33.sce b/2732/CH7/EX7.33/Ex7_33.sce
new file mode 100755
index 000000000..726830a55
--- /dev/null
+++ b/2732/CH7/EX7.33/Ex7_33.sce
@@ -0,0 +1,15 @@
+clc

+//initialization of new variables

+clear

+L=50 //cm

+k=15 //cm

+I=200 //cm^4

+II=40 //cm^4

+d=30 //cm

+Pd=40 //cm

+E=0.6667*10^6 //kg/cm^2

+//calculations

+delta=(100*10/2*16.33+L*d*35+L*k/2*25+d*k/2*45)

+delta1=delta/E

+//Results

+printf('deflection = %.2f mm',delta1*10^1)

diff --git a/2732/CH7/EX7.8/Ex7_8.sce b/2732/CH7/EX7.8/Ex7_8.sce
new file mode 100755
index 000000000..e3f0ea130
--- /dev/null
+++ b/2732/CH7/EX7.8/Ex7_8.sce
@@ -0,0 +1,34 @@
+clc

+//initialization of new variables

+clear

+b=10075 //mm

+h=10 //mm

+p1=7.5

+p2=9

+//part (a)

+ybar=1*p1*0.5+1*p2*5.5

+ybar=ybar/16.5

+xbar=1*p1*0.5+1*p1*4.75

+xbar=xbar/16.5

+printf('part (a)')

+printf('\n Centroid coordinates (x,y) = (%.2f, %.2f) cm',xbar,ybar)

+

+//part (b)

+Ixx=p1*1^3/12+p1*1*(3.23-0.5)^2+1*p2^3/12+p2*1*(5.5-3.23)^2

+Iyy=1*p1^3/12+p1*1*(3.75-2.43)^2+p2*1^3/12+p2*1*(2.43-0.5)^2

+Ixy=p1*1.32*2.73+9*(-1.93)*(-2.27)

+printf('\n part (b)')

+printf('\n Moment of Areas: \n Ixx = %.1f cm^4 \n Iyy = %.1f cm^4 \n Ixy=%.1f cm^4',Ixx,Iyy,Ixy)

+

+//part (c)

+alpha=0.5*atan(2*Ixy/(Iyy-Ixx))

+alpha=alpha*180/%pi

+printf('\n part (c)')

+printf('\n Direction of principal axes:')

+printf('\n alpha = %.2f degrees',alpha)

+

+//part (d)

+Iuu=(Ixx+Iyy)/2+sqrt((Iyy-Ixx)^2/4+Ixy^2)

+Ivv=(Ixx+Iyy)/2-sqrt((Iyy-Ixx)^2/4+Ixy^2)

+printf('\n part (d)')

+printf('\n Iuu = %.2f cm^4 \n Ivv = %.2f cm^4',Iuu,Ivv)