2 files changed, 109 insertions, 0 deletions
diff --git a/764/CH4/EX4.12.b/result4_12.txt b/764/CH4/EX4.12.b/result4_12.txt
new file mode 100755
index 000000000..2c972d98c
--- /dev/null
+++ b/764/CH4/EX4.12.b/result4_12.txt
@@ -0,0 +1,58 @@
+-->//(Design against Static Load) Example 4.12
+ 
+-->//Refer Fig.4.39
+ 
+-->//Force acting on the bracket P (kN)
+ 
+-->P = 5
+ P  =
+ 
+    5.  
+ 
+-->//Ultimate tensile strength of FG200 material Sut (N/mm2)
+ 
+-->Sut = 200
+ Sut  =
+ 
+    200.  
+ 
+-->//Factor of safety fs
+ 
+-->fs = 3.5
+ fs  =
+ 
+    3.5  
+ 
+-->//Angle made by the line of force with the vertical theta (radian)
+ 
+-->theta = (60 * (%pi/180))
+ theta  =
+ 
+    1.0471976  
+ 
+-->//Distance between the point load and the horizontal axis h (mm)
+ 
+-->h = 150
+ h  =
+ 
+    150.  
+ 
+-->//Distance between the point load and the rigid support r (mm)
+ 
+-->r = 300
+ r  =
+ 
+    300.  
+ 
+-->//Ratio of depth to width of the cross-section ratio
+ 
+-->ratio = 2
+ ratio  =
+ 
+    2.  
+ 
+
+Value of t = 35.000000 mm
+
+Area of cross-section = (35.000000 x 70.000000) mm2
+ 
+\ No newline at end of file
diff --git a/764/CH4/EX4.12.b/solution4_12.sce b/764/CH4/EX4.12.b/solution4_12.sce
new file mode 100755
index 000000000..b26cdec54
--- /dev/null
+++ b/764/CH4/EX4.12.b/solution4_12.sce
@@ -0,0 +1,51 @@
+
+//Function to round-up a value such that it is divisible by 5
+function[v] = round_five(w)
+    v = ceil(w)
+    rem = pmodulo(v,5)
+    if (rem ~= 0)
+        v = v + (5 - rem)
+    end
+endfunction
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_12.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_12.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the permissible tensile stress sigmat (N/mm2)
+sigmat = Sut/fs
+//Calculate the horizontal component of the force Ph (N)
+Ph = (P * 1000) * sin(theta)
+//Calculate the vertical component of the force Pv (N)
+Pv = (P * 1000) * cos(theta)
+//Calculate the maximum bending moment Mb (N-mm)
+Mb = (Ph * h) + (Pv * r)
+//Assume the value of t to be 1mm
+t = 1
+//Calculate the value of y (mm)
+y = t
+//Calculate the second moment of area I (mm4)
+I = (t * ((ratio * t)^3))/12
+//Calculate the bending stress sigmab (N/mm2)
+sigmab = (Mb * y)/I
+//Calculate the direct tensile stress D (N/mm2)
+D = Ph/(ratio * (t^2))
+//Coefficients of the resulting cubic equation
+p = [sigmat 0 (-1 * D) (-1 * sigmab)]
+//Calculate the roots to obtain the the true value of t
+r = roots(p)
+real_part = real(r)
+for i = 1:1:3
+    if(real_part(i)>0)
+        t = real_part(i)
+        break
+    end
+end
+t = round_five(t)
+//Print results
+printf('\nValue of t = %f mm\n',t)
+printf('\nArea of cross-section = (%f x %f) mm2\n',t,(ratio * t))