diff options
| author | priyanka | 2015-06-24 15:03:17 +0530 |
|---|---|---|
| committer | priyanka | 2015-06-24 15:03:17 +0530 |
| commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
| tree | ab291cffc65280e58ac82470ba63fbcca7805165 /764/CH4 | |
| download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip | |
initial commit / add all books
Diffstat (limited to '764/CH4')
67 files changed, 2585 insertions, 0 deletions
diff --git a/764/CH4/EX4.1.a/data4_1.sci b/764/CH4/EX4.1.a/data4_1.sci new file mode 100755 index 000000000..24a2b04c6 --- /dev/null +++ b/764/CH4/EX4.1.a/data4_1.sci @@ -0,0 +1,13 @@ +
+//(Design against Static Load) Example 4.1
+//Refer Fig.4.10 on page 85
+//Tensile force acting on two plates P (kN)
+P = 50
+//Tensile yield strength of the plates Syt (N/mm2)
+Syt = 250
+//Factor of safety fs
+fs = 2.5
+//Number of rivets n
+n = 3
+//Length of plate L (mm)
+L = 200
diff --git a/764/CH4/EX4.1.b/result4_1.txt b/764/CH4/EX4.1.b/result4_1.txt new file mode 100755 index 000000000..80fe0019d --- /dev/null +++ b/764/CH4/EX4.1.b/result4_1.txt @@ -0,0 +1,44 @@ +-->//(Design against Static Load) Example 4.1
+
+-->//Refer Fig.4.10 on page 85
+
+-->//Tensile force acting on two plates P (kN)
+
+-->P = 50
+ P =
+
+ 50.
+
+-->//Tensile yield strength of the plates Syt (N/mm2)
+
+-->Syt = 250
+ Syt =
+
+ 250.
+
+-->//Factor of safety fs
+
+-->fs = 2.5
+ fs =
+
+ 2.5
+
+-->//Number of rivets n
+
+-->n = 3
+ n =
+
+ 3.
+
+-->//Length of plate L (mm)
+
+-->L = 200
+ L =
+
+ 200.
+
+
+Diameter of the rivet(d) = 21.000000 mm
+
+Thickness of plates(t) = 4.000000 mm
+
\ No newline at end of file diff --git a/764/CH4/EX4.1.b/solution4_1.sce b/764/CH4/EX4.1.b/solution4_1.sce new file mode 100755 index 000000000..89d1b9949 --- /dev/null +++ b/764/CH4/EX4.1.b/solution4_1.sce @@ -0,0 +1,26 @@ +
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_1.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_1.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Yield strength in shear for the plates Ssy (N/mm2)
+Ssy = (50/100) * Syt
+//Permissible shear stress for rivets Tau (N/mm2)
+Tau = Ssy / fs
+//Diameter of rivets d (mm)
+d = sqrt((P * 1000 * 4)/(n * %pi * Tau))
+//Round up d
+d = ceil(d)
+//Permissible tensile stress for plates Sigma (N/mm2)
+Sigma = (Syt/fs)
+//Thickness of plates t (mm)
+t = (P * 1000)/(Sigma * (L - (n * d)))
+//Round up t
+t = ceil(t)
+//Print results
+printf('\nDiameter of the rivet(d) = %f mm\n',d)
+printf('\nThickness of plates(t) = %f mm\n',t)
diff --git a/764/CH4/EX4.10.a/data4_10.sci b/764/CH4/EX4.10.a/data4_10.sci new file mode 100755 index 000000000..3eb5a2289 --- /dev/null +++ b/764/CH4/EX4.10.a/data4_10.sci @@ -0,0 +1,13 @@ +
+//(Design against Static Load) Example 4.10
+//Refer Fig.4.29
+//Distance between the axis of the column and the load e (mm)
+e = 500
+//Tensile yield strength of FeE200 Syt (N/mm2)
+Syt = 200
+//Factor of safety fs
+fs = 4
+//Load supported by the column P (kN)
+P = 25
+//Ratio of outer diameter to inner diameter ratio
+ratio = 0.8
\ No newline at end of file diff --git a/764/CH4/EX4.10.b/result4_10.txt b/764/CH4/EX4.10.b/result4_10.txt new file mode 100755 index 000000000..9be912889 --- /dev/null +++ b/764/CH4/EX4.10.b/result4_10.txt @@ -0,0 +1,44 @@ +-->//(Design against Static Load) Example 4.10
+
+-->//Refer Fig.4.29
+
+-->//Distance between the axis of the column and the load e (mm)
+
+-->e = 500
+ e =
+
+ 500.
+
+-->//Tensile yield strength of FeE200 Syt (N/mm2)
+
+-->Syt = 200
+ Syt =
+
+ 200.
+
+-->//Factor of safety fs
+
+-->fs = 4
+ fs =
+
+ 4.
+
+-->//Load supported by the column P (kN)
+
+-->P = 25
+ P =
+
+ 25.
+
+-->//Ratio of outer diameter to inner diameter ratio
+
+-->ratio = 0.8
+ ratio =
+
+ 0.8
+
+
+Outer diameter(d0) = 160.000000 mm
+
+Inner diameter(di) = 128.000000 mm
+
\ No newline at end of file diff --git a/764/CH4/EX4.10.b/solution4_10.sce b/764/CH4/EX4.10.b/solution4_10.sce new file mode 100755 index 000000000..ec317545f --- /dev/null +++ b/764/CH4/EX4.10.b/solution4_10.sce @@ -0,0 +1,44 @@ +
+//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_10.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_10.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the permissible tensile stress sigmat (N/mm2)
+sigmat = Syt/fs
+//Let the outer diameter of the cross-section be 1mm d0
+d0 = 1
+//Calculate the direct compressive stress D (N/mm2)
+D = (P * 1000)/((%pi/4)*((1 - (ratio^2))*(d0^2)))
+//Calculate the value of y (mm)
+y = d0/2
+//Calculate the second moment of area I (mm4)
+I = (%pi/64)*((1 - (ratio^4))*(d0^4))
+//Calculate the tensile stress due to bending moment B (N/mm2)
+B = (P * 1000 * e * y)/I
+//Coefficients of the resultant cubic equation
+p = [sigmat 0 D (-1 * B)]
+r = roots(p)
+real_part = real(r)
+for i = 1:1:3
+ if(real_part(i)>0)
+ d0 = real_part(i)
+ break
+ end
+end
+d0 = round_five(d0)
+//Print results
+printf('\nOuter diameter(d0) = %f mm\n',d0)
+printf('\nInner diameter(di) = %f mm\n',(0.8 * d0))
diff --git a/764/CH4/EX4.11.a/data4_11.sci b/764/CH4/EX4.11.a/data4_11.sci new file mode 100755 index 000000000..f1afd7ba4 --- /dev/null +++ b/764/CH4/EX4.11.a/data4_11.sci @@ -0,0 +1,16 @@ +
+//(Design against Static Load) Example 4.11
+//Refer Fig.4.38
+//Tension in the wire rope P (kN)
+P = 5
+//Ultimate tensile strength of FG200 Sut (N/mm2)
+Sut = 200
+//Factor of safety fs
+fs = 2.5
+//Ratio of depth to width of cross-section ratio
+ratio = 2
+//Distance between the load and cantilever support A dist (mm)
+dist = 1500 + 500
+//Distance between the free end and pulley centre pdist (mm)
+pdist = 500
+
\ No newline at end of file diff --git a/764/CH4/EX4.11.b/graph4_11.png b/764/CH4/EX4.11.b/graph4_11.png Binary files differnew file mode 100755 index 000000000..518222afa --- /dev/null +++ b/764/CH4/EX4.11.b/graph4_11.png diff --git a/764/CH4/EX4.11.b/result4_11.txt b/764/CH4/EX4.11.b/result4_11.txt new file mode 100755 index 000000000..2b9fd99fc --- /dev/null +++ b/764/CH4/EX4.11.b/result4_11.txt @@ -0,0 +1,44 @@ +-->//(Design against Static Load) Example 4.11
+
+-->//Refer Fig.4.38
+
+-->//Tension in the wire rope P (kN)
+
+-->P = 5
+ P =
+
+ 5.
+
+-->//Ultimate tensile strength of FG200 Sut (N/mm2)
+
+-->Sut = 200
+ Sut =
+
+ 200.
+
+-->//Factor of safety fs
+
+-->fs = 2.5
+ fs =
+
+ 2.5
+
+-->//Ratio of depth to width of cross-section ratio
+
+-->ratio = 2
+ ratio =
+
+ 2.
+
+-->//Distance between the load and cantilever support A d (mm)
+
+-->d = 1500 + 500
+ d =
+
+ 2000.
+
+
+Width of the cross-section(w) = 60.000000 mm
+
+Depth of the cross-section(d) = 120.000000 mm
+
\ No newline at end of file diff --git a/764/CH4/EX4.11.b/solution4_11.sce b/764/CH4/EX4.11.b/solution4_11.sce new file mode 100755 index 000000000..fe83fab3e --- /dev/null +++ b/764/CH4/EX4.11.b/solution4_11.sce @@ -0,0 +1,63 @@ +
+//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
+
+function[] = plot_format()
+ //Get the handle of current axes
+ g = gca()
+ //Give labels and set label properties
+ g.labels_font_color=5
+ g.font_size=3
+ g.grid=[1,1]
+ g.box="off"
+endfunction
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_11.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_11.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate permissible bending stress sigmab (N/mm2)
+sigmab = Sut/fs
+//Distance of B l (mm)
+l = pdist
+//Increment length by 20mm and plot the bending moment variation over cantilever
+//Infinite for loop
+for i = 1:1:%inf
+ Mb(i) = (P * 1000)*(l)
+ len(i) = l
+ if(l == dist)
+ break
+ end
+ l = l + 20
+end
+plot(len,Mb)
+plot_format()
+title('Variation of bending moment against cantilever length')
+xlabel('Cantilever Length(mm)')
+ylabel('Bending Moment (Hogging) (N-mm)')
+//Calculate maximum bending moment at A Mbmax (N-mm) from the plot
+Mbmax = max(Mb)
+//Assume the width of the cross-section to be 1mm w
+w = 1
+//Calculate the value of y (mm)
+y = w
+//Calculate second moment of area I (mm4)
+I = (w * ((ratio * w)^3))/12
+//Calculate the width of the cross-section (mm)
+w = ((Mbmax * y)/(I * sigmab))^(1/3)
+w = round_five(w)
+//Calculate the depth of the cross-section d (mm)
+d = 2 * w
+//Print results
+printf('\nWidth of the cross-section(w) = %f mm\n',w)
+printf('\nDepth of the cross-section(d) = %f mm\n',d)
diff --git a/764/CH4/EX4.12.a/data4_12.sci b/764/CH4/EX4.12.a/data4_12.sci new file mode 100755 index 000000000..630a480e4 --- /dev/null +++ b/764/CH4/EX4.12.a/data4_12.sci @@ -0,0 +1,17 @@ +
+//(Design against Static Load) Example 4.12
+//Refer Fig.4.39
+//Force acting on the bracket P (kN)
+P = 5
+//Ultimate tensile strength of FG200 material Sut (N/mm2)
+Sut = 200
+//Factor of safety fs
+fs = 3.5
+//Angle made by the line of force with the vertical theta (radian)
+theta = (60 * (%pi/180))
+//Distance between the point load and the horizontal axis h (mm)
+h = 150
+//Distance between the point load and the rigid support r (mm)
+r = 300
+//Ratio of depth to width of the cross-section ratio
+ratio = 2
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))
diff --git a/764/CH4/EX4.13.a/data4_13.sci b/764/CH4/EX4.13.a/data4_13.sci new file mode 100755 index 000000000..f880eb8a1 --- /dev/null +++ b/764/CH4/EX4.13.a/data4_13.sci @@ -0,0 +1,13 @@ +
+//(Design against Satic Load) Example 4.13
+//Refer Fig.4.40
+//Force acting on the overhang crank P (kN)
+P = 1
+//Tensile yield strength of 45C8 material Syt (N/mm2)
+Syt = 380
+//Factor of safety fs
+fs = 2
+//Horizontal distance between the point load and A h (mm)
+h = 500
+//Vertical distance between the point load and A v (mm)
+v = 250
diff --git a/764/CH4/EX4.13.b/result4_13.txt b/764/CH4/EX4.13.b/result4_13.txt new file mode 100755 index 000000000..3948fa064 --- /dev/null +++ b/764/CH4/EX4.13.b/result4_13.txt @@ -0,0 +1,41 @@ +-->//(Design against Satic Load) Example 4.13
+
+-->//Refer Fig.4.40
+
+-->//Force acting on the overhang crank P (kN)
+
+-->P = 1
+ P =
+
+ 1.
+
+-->//Tensile yield strength of 45C8 material Syt (N/mm2)
+
+-->Syt = 380
+ Syt =
+
+ 380.
+
+-->//Factor of safety fs
+
+-->fs = 2
+ fs =
+
+ 2.
+
+-->//Horizontal distance between the point load and A h (mm)
+
+-->h = 500
+ h =
+
+ 500.
+
+-->//Vertical distance between the point load and A v (mm)
+
+-->v = 250
+ v =
+
+ 250.
+
+The value of diameter(d) = 31.061601 mm
+
\ No newline at end of file diff --git a/764/CH4/EX4.13.b/solution4_13.sce b/764/CH4/EX4.13.b/solution4_13.sce new file mode 100755 index 000000000..72a4a0fe2 --- /dev/null +++ b/764/CH4/EX4.13.b/solution4_13.sce @@ -0,0 +1,36 @@ +
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_13.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_13.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the yield strength in shear Ssy (N/mm2)
+Ssy = (50/100)*Syt
+//Calculate the permissible shear stress tau (N/mm2)
+tau = Ssy/fs
+//Calculate the bending moment Mb (N-mm)
+Mb = (P * 1000)*v
+//Calculate the torsional moment Mt (N-mm)
+Mt = (P * 1000)*h
+//Assume value of diameter to be 1mm d
+d = 1
+//Calculate the value of y (mm)
+y = d/2
+//Calculate the second moment of area I (mm4)
+I = (%pi/64)*(d^4)
+//Calculate the polar moment of area J (mm4)
+J = (%pi/32)*(d^4)
+//Calculate the bending shear stress B (N/mm2)
+B = (Mb * y)/I
+//Calculate the torsional shear stress T (N/mm2)
+T = (Mt * (d/2))/J
+//Refer Fig.4.40(c) for Mohr's Circle
+taumax = sqrt(((B/2)^2) + (T^2))
+//Calculate the true value of diameter d (mm)
+d = ((taumax/tau)^(1/3))
+//Print results
+printf('The value of diameter(d) = %f mm\n',d)
+
diff --git a/764/CH4/EX4.14.a/data4_14.sci b/764/CH4/EX4.14.a/data4_14.sci new file mode 100755 index 000000000..1b920daf8 --- /dev/null +++ b/764/CH4/EX4.14.a/data4_14.sci @@ -0,0 +1,13 @@ +
+//(Design against Static Load) Example 4.14
+//Refer Fig.4.41
+//Force acting on the overhang crank P (kN)
+P = 1
+//Tensile yield strength of 30C8 material Syt (N/mm2)
+Syt = 400
+//Factor of safety fs
+fs = 2
+//Horizontal distance between the point load and D h (mm)
+h = 500
+//Vertical distance between the point load and D v (mm)
+v = 50 + 25 + 100
diff --git a/764/CH4/EX4.14.b/result4_14.txt b/764/CH4/EX4.14.b/result4_14.txt new file mode 100755 index 000000000..3a5c7eabb --- /dev/null +++ b/764/CH4/EX4.14.b/result4_14.txt @@ -0,0 +1,41 @@ +-->//(Design against Static Load) Example 4.14
+
+-->//Refer Fig.4.41
+
+-->//Force acting on the overhang crank P (kN)
+
+-->P = 1
+ P =
+
+ 1.
+
+-->//Tensile yield strength of 30C8 material Syt (N/mm2)
+
+-->Syt = 400
+ Syt =
+
+ 400.
+
+-->//Factor of safety fs
+
+-->fs = 2
+ fs =
+
+ 2.
+
+-->//Horizontal distance between the point load and D h (mm)
+
+-->h = 500
+ h =
+
+ 500.
+
+-->//Vertical distance between the point load and D v (mm)
+
+-->v = 50 + 25 + 100
+ v =
+
+ 175.
+
+The value of diameter(d) = 29.992391 mm
+
\ No newline at end of file diff --git a/764/CH4/EX4.14.b/solution4_14.sce b/764/CH4/EX4.14.b/solution4_14.sce new file mode 100755 index 000000000..63858ab65 --- /dev/null +++ b/764/CH4/EX4.14.b/solution4_14.sce @@ -0,0 +1,36 @@ +
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_14.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_14.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the yield strength in shear Ssy (N/mm2)
+Ssy = (50/100)*Syt
+//Calculate the permissible shear stress tau (N/mm2)
+tau = Ssy/fs
+//Calculate the bending moment Mb (N-mm)
+Mb = (P * 1000)*v
+//Calculate the torsional moment Mt (N-mm)
+Mt = (P * 1000)*h
+//Assume value of diameter to be 1mm d
+d = 1
+//Calculate the value of y (mm)
+y = d/2
+//Calculate the second moment of area I (mm4)
+I = (%pi/64)*(d^4)
+//Calculate the polar moment of area J (mm4)
+J = (%pi/32)*(d^4)
+//Calculate the bending shear stress B (N/mm2)
+B = (Mb * y)/I
+//Calculate the torsional shear stress T (N/mm2)
+T = (Mt * (d/2))/J
+//Refer Fig.4.40(c) for Mohr's Circle
+taumax = sqrt(((B/2)^2) + (T^2))
+//Calculate the true value of diameter d (mm)
+d = ((taumax/tau)^(1/3))
+//Print results
+printf('The value of diameter(d) = %f mm\n',d)
+
diff --git a/764/CH4/EX4.15.a/data4_15.sci b/764/CH4/EX4.15.a/data4_15.sci new file mode 100755 index 000000000..b25032893 --- /dev/null +++ b/764/CH4/EX4.15.a/data4_15.sci @@ -0,0 +1,20 @@ +
+//(Design against Static Load) Example 4.15
+//Refer Fig.4.50 on page 121
+//Safety valve blow-off pressure p (MPa)
+p = 1.5
+//Effective diameter of the valve opening d (mm)
+de = 50
+//Distance between fulcrum and lever dead weights l1 (mm)
+l1 = 1000
+//Distance between the fulcrum and the pin connecting the valve spindle to the lever l2 (mm)
+l2 = 100
+//Tensile yield strength of 30C8 Syt (N/mm2)
+Syt = 400
+//Factor of safety fs
+fs = 5
+//Permissible bearing pressure at the pin in the lever Bp (N/mm2)
+Bp = 25
+//Ratio of width to thickness ratio
+ratio = 3
+
\ No newline at end of file diff --git a/764/CH4/EX4.15.b/result4_15.txt b/764/CH4/EX4.15.b/result4_15.txt new file mode 100755 index 000000000..44ecb8d05 --- /dev/null +++ b/764/CH4/EX4.15.b/result4_15.txt @@ -0,0 +1,74 @@ +-->//(Design against Static Load) Example 4.15
+
+-->//Refer Fig.4.50 on page 121
+
+-->//Safety valve blow-off pressure p (MPa)
+
+-->p = 1.5
+ p =
+
+ 1.5
+
+-->//Effective diameter of the valve opening d (mm)
+
+-->de = 50
+ de =
+
+ 50.
+
+-->//Distance between fulcrum and lever dead weights l1 (mm)
+
+-->l1 = 1000
+ l1 =
+
+ 1000.
+
+-->//Distance between the fulcrum and the pin connecting the valve spindle to the lever l2 (mm)
+
+-->l2 = 100
+ l2 =
+
+ 100.
+
+-->//Tensile yield strength of 30C8 Syt (N/mm2)
+
+-->Syt = 400
+ Syt =
+
+ 400.
+
+-->//Factor of safety fs
+
+-->fs = 5
+ fs =
+
+ 5.
+
+-->//Permissible bearing pressure at the pin in the lever Bp (N/mm2)
+
+-->Bp = 25
+ Bp =
+
+ 25.
+
+-->//Ratio of width to thickness ratio
+
+-->ratio = 3
+ ratio =
+
+ 3.
+
+
+
+The diamater of the pin(dp) = 11.000000 mm
+
+The length of the pin(lp) = 20.000000 mm
+
+The thickness of the lever cross-section(b) = 15.000000 mm
+
+The width of the lever cross-sectio(d) = 45.000000 mm
+
+The dead weight(P) = 294.524311 mm
+
+The lever design is safe
+
\ No newline at end of file diff --git a/764/CH4/EX4.15.b/solution4_15.sce b/764/CH4/EX4.15.b/solution4_15.sce new file mode 100755 index 000000000..49c6f60b5 --- /dev/null +++ b/764/CH4/EX4.15.b/solution4_15.sce @@ -0,0 +1,82 @@ +
+//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) then
+ v = v + (5 - rem)
+ end
+endfunction
+
+//Function to round-up a value such that it is divisible by 10
+function[v] = round_ten(w)
+ v = ceil(w)
+ rem = pmodulo(v,10)
+ if (rem ~= 0) then
+ v = v + (10 - rem)
+ end
+endfunction
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_15.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_15.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the permissible stresses for lever and pin sigmat (N/mm2)
+sigmat = Syt/fs
+//Calculate the yield strength in shear Ssy (N/mm2)
+Ssy = (50/100)*Syt
+//Calculate the permissible stress in shear tau (N/mm2)
+tau = Ssy/fs
+//Calculate the maximum steam load F (N)
+F = (%pi/4)*((de^2) * p)
+//Calculate the dead weight P (N)
+P = (F * l2)/l1
+//Calculate reaction at fulcrum R (N)
+R = F - P
+//Assume length and diameter of pin to be equal (lp = dp)
+dp = sqrt(F/Bp)
+dp = ceil(dp)
+lp = dp
+//Calculate the shear stress in pin tau1 (N/mm2)
+tau1 = F/(2 * (%pi/4) * (dp^2))
+//Calculate gunmetal bush thickness t (mm)
+t = 2
+//Calculate inside diameter of the boss di (mm)
+di = dp + (2 * t)
+//Calculate the outside diameter of the boss d0 (mm)
+d0 = 2 * di
+//Calculate the maximum banding moment Mb (N-mm)
+Mb = P*(l1 - l2)
+//Assume the thickness of the cross-section to be 1mm b
+b = 1
+//Calculate the width of the cross-section d (mm)
+d = ratio * b
+//Calculate the value of y (mm)
+y = d/2
+//Calculate the second moment of area I (mm4)
+I = (b * ((ratio * b)^3))/12
+//Calculate the true value of b (mm)
+b = ((Mb * y)/(sigmat * I))^(1/3)
+b = round_five(b)
+//Calculate the true value of d (mm)
+d = ratio * b
+lp = round_ten(lp)
+//For lever cross-section
+y1 = d/2
+I1 = ((b * (d^3)) + ((lp - b) * (d0^3)) - (lp * (di^3)))/12
+//Calculate the bending stress for the modified design sigmab (N/mm2)
+sigmab = (Mb * y1)/I1
+//Print results
+printf('\nThe diamater of the pin(dp) = %f mm\n',dp)
+printf('\nThe length of the pin(lp) = %f mm\n',lp)
+printf('\nThe thickness of the lever cross-section(b) = %f mm\n',b)
+printf('\nThe width of the lever cross-sectio(d) = %f mm\n',d)
+printf('\nThe dead weight(P) = %f mm\n',P)
+//Check for design safety
+if(tau1<tau & sigmab<sigmat)
+ printf('\nThe lever design is safe\n')
+end
diff --git a/764/CH4/EX4.16.a/data4_16.sci b/764/CH4/EX4.16.a/data4_16.sci new file mode 100755 index 000000000..353c23f28 --- /dev/null +++ b/764/CH4/EX4.16.a/data4_16.sci @@ -0,0 +1,19 @@ +
+//(Design against Static Load) Example 4.16
+//Refer Fig.4.53 on page 123
+//Load to be raised F (kN)
+F = 5
+//Length of short arm l1 (mm)
+l1 = 100
+//Length of long arm l2 (mm)
+l2 = 450
+//Tensile yield strength of 30C8 Syt (N/mm2)
+Syt = 400
+//Factor of safety fs
+fs = 5
+//Permissible bearing pressure at the pin in the lever p (N/mm2)
+p = 10
+//Ratio of width to thickness ratio
+ratio = 3
+//Ratio of length to diameter of the fulcrum pin rf
+rf = 1.25
diff --git a/764/CH4/EX4.16.b/result4_16.txt b/764/CH4/EX4.16.b/result4_16.txt new file mode 100755 index 000000000..4fa2729aa --- /dev/null +++ b/764/CH4/EX4.16.b/result4_16.txt @@ -0,0 +1,71 @@ +-->//(Design against Static Load) Example 4.16
+
+-->//Refer Fig.4.53 on page 123
+
+-->//Load to be raised F (kN)
+
+-->F = 5
+ F =
+
+ 5.
+
+-->//Length of short arm l1 (mm)
+
+-->l1 = 100
+ l1 =
+
+ 100.
+
+-->//Length of long arm l2 (mm)
+
+-->l2 = 450
+ l2 =
+
+ 450.
+
+-->//Tensile yield strength of 30C8 Syt (N/mm2)
+
+-->Syt = 400
+ Syt =
+
+ 400.
+
+-->//Factor of safety fs
+
+-->fs = 5
+ fs =
+
+ 5.
+
+-->//Permissible bearing pressure at the pin in the lever p (N/mm2)
+
+-->p = 10
+ p =
+
+ 10.
+
+-->//Ratio of width to thickness ratio
+
+-->ratio = 3
+ ratio =
+
+ 3.
+
+-->//Ratio of length to diameter of the fulcrum pin rf
+
+-->rf = 1.25
+ rf =
+
+ 1.25
+
+
+Diameter of the fulcrum pin(df) = 20.242469 mm
+
+Length of the fulcrum pin(lf) = 25.303086 mm
+
+Force to be applied at the long arm end(P) = 1111.111111 mm
+
+Thickness of the cross-section of the lever(b) = 16.091490 mm
+
+Width of the cross-section of the lever(d) = 48.274469 mm
+
\ No newline at end of file diff --git a/764/CH4/EX4.16.b/solution4_16.sce b/764/CH4/EX4.16.b/solution4_16.sce new file mode 100755 index 000000000..8408277ca --- /dev/null +++ b/764/CH4/EX4.16.b/solution4_16.sce @@ -0,0 +1,52 @@ +
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_16.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_16.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the permissible stresses for lever and pin sigmat (N/mm2)
+sigmat = Syt/fs
+//Calculate the yield strength in shear Ssy (N/mm2)
+Ssy = (50/100)*Syt
+//Calculate the permissible stress in shear tau (N/mm2)
+tau = Ssy/fs
+//Calculate force acting on the lever P (N)
+P = (F * 1000 * l1)/l2
+//Calculate the reaction at the fulcrum R (N)
+R = sqrt(((F * 1000)^2) + (P^2))
+//Calculate the diameter of the fulcrum pin df (mm)
+df = (R/(p * rf))^(1/2)
+//Calculate the length of the fulcrum pin lf (mm)
+lf = df * rf
+//Calculate the shear stress in pin tau1 (N/mm2)
+tau1 = R/(2 * (%pi/4) * (df^2))
+//Dimensions of the boss
+//Inner diameter di (mm)
+di = ceil(df)
+//Outer diameter d0 (mm)
+d0 = 2 * di
+//Length l (mm)
+l = ceil(lf)
+//Dimensions of cross-section of lever
+//Calculate the maximum bending moment Mb (N-mm)
+Mb = (F * 1000)*l1
+//Assume the thickness b to be 1mm
+b = 1
+d = ratio * b
+//Calculate the value of y (mm)
+y = d/2
+//Calculate second moment of area I (mm4)
+I = (b * (d^3))/12
+//Calculate the true value of b (mm)
+b = ((Mb * y)/(sigmat * I))^(1/3)
+//Calculate the width of lever cross-section d (mm)
+d = ratio * b
+//Print results
+printf('\nDiameter of the fulcrum pin(df) = %f mm\n',df)
+printf('\nLength of the fulcrum pin(lf) = %f mm\n',lf)
+printf('\nForce to be applied at the long arm end(P) = %f mm\n',P)
+printf('\nThickness of the cross-section of the lever(b) = %f mm\n',b)
+printf('\nWidth of the cross-section of the lever(d) = %f mm\n',d)
diff --git a/764/CH4/EX4.17.a/data4_17.sci b/764/CH4/EX4.17.a/data4_17.sci new file mode 100755 index 000000000..8706943e5 --- /dev/null +++ b/764/CH4/EX4.17.a/data4_17.sci @@ -0,0 +1,24 @@ +
+//(Design against Static Load) Example 4.17
+//Refer Fig.4.54 on page 124
+//For beam
+//Tensile yield strength of FeE250 Syt (N/mm2)
+Syt = 250
+//Factor of safety fs
+fs = 5
+//Ratio of width to thickness of the beam cross-section ratio
+ratio = 2
+//For vessel
+//Ultimate tensile strength of FG200 Sut(N/mm2)
+Sut = 200
+//Internal guage pressure p (MPa)
+p = 0.25
+//Diameter D (mm)
+D = 500
+//Refer table for ISO Metric threads-Coarse series
+//Area mm2
+area = [561 817 1120 1470]
+//Pitch (mm)
+pitch = [3.5 4 4.5 5]
+//Distance between vessel centre line and beam end l (mm)
+l = 325
diff --git a/764/CH4/EX4.17.b/result4_17.txt b/764/CH4/EX4.17.b/result4_17.txt new file mode 100755 index 000000000..ea87d2b4f --- /dev/null +++ b/764/CH4/EX4.17.b/result4_17.txt @@ -0,0 +1,92 @@ +-->//(Design against Static Load) Example 4.17
+
+-->//Refer Fig.4.54 on page 124
+
+-->//For beam
+
+-->//Tensile yield strength of FeE250 Syt (N/mm2)
+
+-->Syt = 250
+ Syt =
+
+ 250.
+
+-->//Factor of safety fs
+
+-->fs = 5
+ fs =
+
+ 5.
+
+-->//Ratio of width to thickness of the beam cross-section ratio
+
+-->ratio = 2
+ ratio =
+
+ 2.
+
+-->//For vessel
+
+-->//Ultimate tensile strength of FG200 Sut(N/mm2)
+
+-->Sut = 200
+ Sut =
+
+ 200.
+
+-->//Internal guage pressure p (MPa)
+
+-->p = 0.25
+ p =
+
+ 0.25
+
+-->//Diameter D (mm)
+
+-->D = 500
+ D =
+
+ 500.
+
+-->//Refer table for ISO Metric threads-Coarse series
+
+-->//Area mm2
+
+-->area = [561 817 1120 1470]
+ area =
+
+ 561. 817. 1120. 1470.
+
+-->//Pitch (mm)
+
+-->pitch = [3.5 4 4.5 5]
+ pitch =
+
+ 3.5 4. 4.5 5.
+
+-->//Distance between vessel centre line and beam end l (mm)
+
+-->l = 325
+ l =
+
+ 325.
+
+
+Select M42 screw
+
+The nominal diameter and pitch of the screw are 42.000000 mm and 4.500000 mm respectively
+
+The width of the beam cross-section(h) = 130.000000 mm
+
+The thickness of the beam cross-section(b) = 65.000000 mm
+
+The value of d0 = 107.000000 mm
+
+The diameter of the pins(d) = 25.000000 mm
+
+The diameter of links 1 and 2(d2) = 25.000000 mm
+
+The width of the vessel extension cross-section(h1) = 90.000000 mm
+
+The thickness of the vessel extension cross-section(b1) = 45.000000 mm
+
\ No newline at end of file diff --git a/764/CH4/EX4.17.b/solution4_17.sce b/764/CH4/EX4.17.b/solution4_17.sce new file mode 100755 index 000000000..9dedee473 --- /dev/null +++ b/764/CH4/EX4.17.b/solution4_17.sce @@ -0,0 +1,100 @@ +
+//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) then
+ v = v + (5 - rem)
+ end
+endfunction
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_17.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_17.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the permissible stresses for steel parts sigmat1 (N/mm2)
+sigmat1 = Syt/fs
+//Compressive stress sigmac1 (N/mm2)
+sigmac1 = sigmat1
+//Yield strength in shear Ssy (N/mm2)
+Ssy = (50/100)*Syt
+//Shear stress tau1 (N/mm2)
+tau1 = Ssy/fs
+//For cast iron parts
+//Permissible tensile stress sigmat2 (N/mm2)
+sigmat2 = Sut/fs
+//Calculate the force acting on the cover F (N)
+F = (%pi/4)*(p * (D^2))
+//Calculate the stressed area of the screw a (mm2)
+a = F/sigmac1
+//Select the standard screw from the table
+for i = 1:1:4
+ if(a < area(i))
+ break
+ end
+end
+//Print the size of the standard screw selected
+if(i==1)
+ printf('\nSelect M30 screw\n')
+ d1 = 30
+elseif(i==2)
+ printf('\nSelect M36 screw\n')
+ d1 = 36
+elseif(i==3)
+ printf('\nSelect M42 screw\n')
+ d1 = 42
+else
+ printf('\nSelect M48 screw\n')
+ d1 = 48
+end
+//Calculate the maximum bending moment on the beam Mb (N-mm)
+Mb = l * (F/2)
+//Assume the thickness of the beam to be 1mm b
+b = 1
+h = ratio * b
+//Calculate the second moment of area I (mm4)
+I = (b * ((ratio * b)^3))/12
+//Calculate the value of y (mm)
+y = h/2
+//Calculate the true value of b (mm)
+b = ((Mb * y)/(sigmat1 * I))^(1/3)
+b = round_five(b)
+//Calculate the true value of h (mm)
+h = ratio * b
+//Calculate the value of d0 (mm)
+d0 = d1 + (b/2) + (b/2)
+//Calculate the diameter of the pins d (mm)
+d = (((F/2) * 4)/(2 * %pi * tau1))^(1/2)
+//Calculate the diameter of links 1 and 2 d2 (mm)
+d2 = (((F/2) * 4)/(%pi * sigmat1))^(1/2)
+//Maximum length of the vessel extension e (mm)
+e = l - (D/2)
+//Calculate the maximum bending moment on the extension Mb1 (N-mm)
+Mb1 = (F/2)*e
+//Assume the width of extension to be 1mm b1
+b1 = 1
+h1 = ratio * b1
+//Calculate y1 (mm)
+y1 = h1/2
+//Calculate I1 (mm4)
+I1 = (b1 * (h1^3))/12
+//Calculate the true value of b1 (mm)
+b1 = ((Mb1 * y1)/(sigmat2 * I1))^(1/3)
+b1 = round_five(b1)
+//Calculate the true value of h1
+h1 = ratio * b1
+//Print results
+printf('\nThe nominal diameter and pitch of the screw are %f mm and %f mm respectively\n',d1,pitch(i))
+printf('\nThe width of the beam cross-section(h) = %f mm\n',h)
+printf('\nThe thickness of the beam cross-section(b) = %f mm\n',b)
+printf('\nThe value of d0 = %f mm\n',d0)
+printf('\nThe diameter of the pins(d) = %f mm\n',d)
+printf('\nThe diameter of links 1 and 2(d2) = %f mm\n',d2)
+printf('\nThe width of the vessel extension cross-section(h1) = %f mm\n',h1)
+printf('\nThe thickness of the vessel extension cross-section(b1) = %f mm\n',b1)
+
+
diff --git a/764/CH4/EX4.18.a/data4_18.sci b/764/CH4/EX4.18.a/data4_18.sci new file mode 100755 index 000000000..55ce19a95 --- /dev/null +++ b/764/CH4/EX4.18.a/data4_18.sci @@ -0,0 +1,25 @@ +
+//(Design against Static Load) Example 4.18
+//Refer Fig.4.57 on page 126
+//Tensile yield strength of FeE250 Syt (N/mm2)
+Syt = 250
+//Factor of safety fs
+fs = 5
+//Diameter of bars to be sheared D (mm)
+D = 6.25
+//Ultimate shear strength of the material Sus (N/mm2)
+Sus = 350
+//Permissible bearing pressure on the pins p (N/mm2)
+p = 10
+//Pin length to diameter ratio r1
+r1 = 1.25
+//Link cross-section width to thickness ratio r2
+r2 = 2
+//Distance between pinA and pinB l1 (mm)
+l4 = 400
+//Distance between bar to be sheared and pinA & force application point and pinC l2 (mm)
+l2 = 100
+//Distance between the force applied and pinD l3 (mm)
+l3 = 1000
+//Thickness of gunmetal bush over pin C t (mm)
+t = 2.5
diff --git a/764/CH4/EX4.18.b/result4_18.txt b/764/CH4/EX4.18.b/result4_18.txt new file mode 100755 index 000000000..2579a2ab0 --- /dev/null +++ b/764/CH4/EX4.18.b/result4_18.txt @@ -0,0 +1,94 @@ +-->//(Design against Static Load) Example 4.18
+
+-->//Refer Fig.4.57 on page 126
+
+-->//Tensile yield strength of FeE250 Syt (N/mm2)
+
+-->Syt = 250
+ Syt =
+
+ 250.
+
+-->//Factor of safety fs
+
+-->fs = 5
+ fs =
+
+ 5.
+
+-->//Diameter of bars to be sheared D (mm)
+
+-->D = 6.25
+ D =
+
+ 6.25
+
+-->//Ultimate shear strength of the material Sus (N/mm2)
+
+-->Sus = 350
+ Sus =
+
+ 350.
+
+-->//Permissible bearing pressure on the pins p (N/mm2)
+
+-->p = 10
+ p =
+
+ 10.
+
+-->//Pin length to diameter ratio r1
+
+-->r1 = 1.25
+ r1 =
+
+ 1.25
+
+-->//Link cross-section width to thickness ratio r2
+
+-->r2 = 2
+ r2 =
+
+ 2.
+
+-->//Distance between pinA and pinB l1 (mm)
+
+-->l4 = 400
+ l4 =
+
+ 400.
+
+-->//Distance between bar to be sheared and pinA & force application point and pinC l2 (mm)
+
+-->l2 = 100
+ l2 =
+
+ 100.
+
+-->//Distance between the force applied and pinD l3 (mm)
+
+-->l3 = 1000
+ l3 =
+
+ 1000.
+
+-->//Thickness of gunmetal bush over pin C t (mm)
+
+-->t = 2.5
+ t =
+
+ 2.5
+
+
+Diameter of pins(d1) = 15.000000 mm
+
+Length of pins(l1) = 20.000000 mm
+
+Diameter of link(d) = 10.000000 mm
+
+Width of lever cross-section(h) = 40.000000 mm
+
+Thickness of lever cross-section(b) = 20.000000 mm
+
+The design is safe
+
\ No newline at end of file diff --git a/764/CH4/EX4.18.b/solution4_18.sce b/764/CH4/EX4.18.b/solution4_18.sce new file mode 100755 index 000000000..aecb236f9 --- /dev/null +++ b/764/CH4/EX4.18.b/solution4_18.sce @@ -0,0 +1,81 @@ +
+//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) then
+ v = v + (5 - rem)
+ end
+endfunction
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_18.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_18.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate permissible stresses
+//Permissible tensile stress sigmat (N/mm2)
+sigmat = Syt/fs
+//Calculate the yield strength in shear Ssy (N/mm2)
+Ssy = (50/100)*Syt
+//Calculate the permissible stress in shear tau (N/mm2)
+tau = Ssy/fs
+//Maximum force required to shear the bar Ps (N)
+Ps = (%pi/4)*(D^2)*Sus
+//Calculate the force acting on each pin P1 (N)
+P1 = (Ps * l2)/l4
+//Calculate reaction at pinA RA (N)
+RA = Ps - P1
+//Calculate the force applied at lever P (N)
+P = (P1 * l2)/l3
+//Calculate reaction at pinD RD (N)
+RD = P1 - P
+//Calculate the diameter of pins d1 (mm)
+d1 = (P1/(p * r1))^(1/2)
+d1 = round_five(d1)
+//Calculate the length of pin l1 (mm)
+l1 = r1 * d1
+l1 = round_five(l1)
+//Calculate the shear stress in pin tau1 (N/mm2)
+tau1 = (P1 * 4)/(2 * %pi * (d1^2))
+//Calculate the diameter of the link d (mm)
+d = ((P1 * 4)/(%pi * sigmat))^(1/2)
+d = round_five(d)
+//Calculate the maximum bending moment on lever Mb (N-mm)
+Mb = P * (l3 - l2)
+//Assume the thickness of the cross-section to be 1mm b
+b = 1
+//Calculate the width of the cross-section h (mm)
+h = r2 * b
+//Calculate the value of y (mm)
+y = h/2
+//Calculate the second moment of area I (mm4)
+I = (b * ((r2 * b)^3))/12
+//Calculate the true value of b (mm)
+b = ((Mb * y)/(sigmat * I))^(1/3)
+b = round_five(b)
+//Calculate the true value of h (mm)
+h = r2 * b
+//Calculate the inner diameter of the boss di (mm)
+di = d1 + (2 * t)
+//Calculate the outer diameter of the boss d0 (mm)
+d0 = 2 * di
+//Calculate the second moment of area I1 (mm4)
+I1 = ((l1 * (d0^3)) - (l1 * (di^3)))/12
+//Calculate y1 (mm)
+y1 = d0/2
+//Calculate the bending stress B (N/mm2)
+B = (Mb * y1)/I1
+//Print results
+printf('\nDiameter of pins(d1) = %f mm\n',d1)
+printf('\nLength of pins(l1) = %f mm\n',l1)
+printf('\nDiameter of link(d) = %f mm\n',d)
+printf('\nWidth of lever cross-section(h) = %f mm\n',h)
+printf('\nThickness of lever cross-section(b) = %f mm\n',b)
+//Chech design
+if((tau1 < tau) & (B - sigmat)<2)
+ printf('\nThe design is safe\n')
+end
\ No newline at end of file diff --git a/764/CH4/EX4.19.a/data4_19.sci b/764/CH4/EX4.19.a/data4_19.sci new file mode 100755 index 000000000..87310fd51 --- /dev/null +++ b/764/CH4/EX4.19.a/data4_19.sci @@ -0,0 +1,14 @@ +
+//(Design against Static Load) Example 4.19
+//Refer Fig.4.65 and Fig.4.66
+//Tensile yield strength of 45C8 Syt (N/mm2)
+Syt = 380
+//Factor of safety fs
+fs = 3.5
+//All dimensions in mm
+bi = 90
+bo = 30
+h = 120
+Ro = 170
+Ri = 50
+
diff --git a/764/CH4/EX4.19.b/result4_19.txt b/764/CH4/EX4.19.b/result4_19.txt new file mode 100755 index 000000000..7c458f537 --- /dev/null +++ b/764/CH4/EX4.19.b/result4_19.txt @@ -0,0 +1,49 @@ +-->//(Design against Static Load) Example 4.19
+
+-->//Refer Fig.4.65 and Fig.4.66
+
+-->//Tensile yield strength of 45C8 Syt (N/mm2)
+
+-->Syt = 380
+ Syt =
+
+ 380.
+
+-->//Factor of safety fs
+
+-->fs = 3.5
+ fs =
+
+ 3.5
+
+-->//All dimensions in mm
+
+-->bi = 90
+ bi =
+
+ 90.
+
+-->bo = 30
+ bo =
+
+ 30.
+
+-->h = 120
+ h =
+
+ 120.
+
+-->Ro = 170
+ Ro =
+
+ 170.
+
+-->Ri = 50
+ Ri =
+
+ 50.
+
+
+
+The load carrying capacity(P) = 94828.181768 N
+
\ No newline at end of file diff --git a/764/CH4/EX4.19.b/solution4_19.sce b/764/CH4/EX4.19.b/solution4_19.sce new file mode 100755 index 000000000..d7fa73adc --- /dev/null +++ b/764/CH4/EX4.19.b/solution4_19.sce @@ -0,0 +1,29 @@ +
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_19.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_19.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the permissible stresses for steel parts sigmat (N/mm2)
+sigmat = Syt/fs
+//Calculate the eccentricity e (mm)
+Rn = (((bi + bo)/2)*h)/((((bi * Ro) - (bo * Ri))/h)*log(Ro/Ri)-(bi - bo))
+R = Ri + ((h*(bi + (2 * bo)))/(3 * (bi + bo)))
+e = R - Rn
+hi = Rn - Ri
+A = (h * (bi + bo))/2
+//Assume the value load P to be 1N
+P = 1
+//Calculate the bending moment Mb (N-mm)
+Mb = R * P
+//Calculate bending stress at the inner fibre B (N/mm2)
+B = (Mb * hi)/(A * e * Ri)
+//Calculate the direct tensile stress T (N/mm2)
+T = P/A
+//Calculate the actual load carrying capacity P (N)
+P = sigmat/(B + T)
+//Print results
+printf('\nThe load carrying capacity(P) = %f N\n',P)
diff --git a/764/CH4/EX4.2.a/data4_2.sci b/764/CH4/EX4.2.a/data4_2.sci new file mode 100755 index 000000000..ea5c0fb30 --- /dev/null +++ b/764/CH4/EX4.2.a/data4_2.sci @@ -0,0 +1,11 @@ +
+//(Design against Static Load) Example 4.2
+//Refer section 4.10 on page 90
+//Axial tensile force acting on acting on each rod P (kN)
+P = 50
+//Tensile yield strength of the material Syt (N/mm2)
+Syt = 400
+//Factor of safety for rod f1
+f1 = 6
+//Factor of safety for cotter f2
+f2 = 4
diff --git a/764/CH4/EX4.2.b/result4_2.txt b/764/CH4/EX4.2.b/result4_2.txt new file mode 100755 index 000000000..61640e334 --- /dev/null +++ b/764/CH4/EX4.2.b/result4_2.txt @@ -0,0 +1,51 @@ +-->//(Design against Static Load) Example 4.2
+
+-->//Refer section 4.10 on page 90
+
+-->//Axial tensile force acting on acting on each rod P (kN)
+
+-->P = 50
+ P =
+
+ 50.
+
+-->//Tensile yield strength of the material Syt (N/mm2)
+
+-->Syt = 400
+ Syt =
+
+ 400.
+
+-->//Factor of safety for rod f1
+
+-->f1 = 6
+ f1 =
+
+ 6.
+
+-->//Factor of safety for cotter f2
+
+-->f2 = 4
+ f2 =
+
+ 4.
+
+
+Diameter of the rods(d) = 31.000000 mm
+
+Thickness of cotter(t) = 10.000000 mm
+
+Diameter of spigot(d2) = 40.000000 mm
+
+Outer diameter of socket(d1) = 55.000000 mm
+
+Diameter of spigot collar(d3) = 47.000000 mm
+
+Diameter of socket collar(d4) = 80.000000 mm
+
+Width of cotter(b) = 50.000000 mm
+
+Thickness of spigot collar(t1) = 15.000000 mm
+
+Design of cotter joint is safe
+
\ No newline at end of file diff --git a/764/CH4/EX4.2.b/solution4_2.sce b/764/CH4/EX4.2.b/solution4_2.sce new file mode 100755 index 000000000..0f4a52c27 --- /dev/null +++ b/764/CH4/EX4.2.b/solution4_2.sce @@ -0,0 +1,105 @@ +
+//Function to calculate roots of a quadratic equation
+function[r] = quadratic(a,b,c)
+ //Calculate discriminant D
+ D = (b^2)-(4 * a * c)
+ r1 = ((-1 * b)+ sqrt(D))/(2 * a)
+ r2 = ((-1 * b)- sqrt(D))/(2 * a)
+ if(r1 > 0)
+ r = r1
+ else
+ r = r2
+ end
+endfunction
+
+//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)
+ v = v + (5 - rem)
+endfunction
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_2.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_2.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Refer Fig.4.17 on page 93
+//For rod
+//Permissible tensile stress sr1 (N/mm2)
+sr1 = Syt/f1
+//Permissible compressive stress sr2 (N/mm2)
+sr2 = 2 * sr1
+//Yield strength of the material in shear Ssy (N/mm2)
+Ssy = (50/100)*Syt
+//Permissible shear stress tau1 (N/mm2)
+tau1 = Ssy/f1
+//For cotter
+//Permissible tensile stress sc1 (N/mm2)
+sc1 = Syt/f2
+//Permissible shear stress tau2 (N/mm2)
+tau2 = Ssy/f2
+//Calculate diameter of the rods d (mm)
+d = sqrt((4 * P * 1000)/(%pi * sr1))
+//Round up d
+d = ceil(d)
+//Calculate thickness of cotter t (mm)
+t = 0.31 * d
+//Round up t
+t = ceil(t)
+//Calculate the diameter of spigot d2 (mm)
+//a,b,c are the coefficients of the resulting quadratic equation
+//ax^2 + bx + c = 0 x=d2
+a = (%pi/4)*sr1
+b = -1* t * sr1
+c = -1 * P * 1000
+//Call the declared functions
+d2 = quadratic(a, b, c)
+d2 = round_five(d2)
+//Calculate outer daimeter of the socket d1 (mm)
+//a,b,c are the coefficients of the resulting quadratic equation
+//ax^2 + bx + c = 0 x=d1
+a = (%pi/4)*(sr1)
+b = (-1)*(t * sr1)
+c = (((-1 * (%pi/4) * (d2^2))+(t * d2))*sr1)-(P * 1000)
+//Call the declared functions
+d1 = quadratic(a, b, c)
+d1 = round_five(d1)
+//Calculate diameters of spigot collar d3 and socket collar d4 (mm)
+d3 = 1.5 * d
+d3 = ceil(d3)
+d4 = 2.4 * d
+d4 = round_five(d4)
+//Calculate dimensions a1 and c1 (mm)
+a1 = 0.75 * d
+a1 = ceil(a1)
+c1 = 0.75 * d
+c1 = ceil(c1)
+//Calculate width of cotter b (mm)
+b = (P * 1000)/(2 * tau2 * t)
+//Calculate thickness of spigot collar t1 mm
+t1 = 0.45 * d
+t1 = round_five(t1)
+//Print results
+printf('\nDiameter of the rods(d) = %f mm\n',d)
+printf('\nThickness of cotter(t) = %f mm\n',t)
+printf('\nDiameter of spigot(d2) = %f mm\n',d2)
+printf('\nOuter diameter of socket(d1) = %f mm\n',d1)
+printf('\nDiameter of spigot collar(d3) = %f mm\n',d3)
+printf('\nDiameter of socket collar(d4) = %f mm\n',d4)
+printf('\nWidth of cotter(b) = %f mm\n',b)
+printf('\nThickness of spigot collar(t1) = %f mm\n',t1)
+//Check for crushing and shear stresses in spigot end
+sigc = (P * 1000)/(t * d2)
+tauc = (P * 1000)/(2 * a1 * d2)
+//Check for crushing and shear stresses in socket end
+sigs = (P * 1000)/((d4 - d2) * t)
+taus = (P * 1000)/(2 * (d4 - d2) * c1)
+if ((sigc < sr2) & (tauc < tau1) & (sigs < sr2) & (taus < tau1))
+ printf('\nDesign of cotter joint is safe\n')
+else
+ printf('\nDesign of cotter joint is not safe\n')
+end
diff --git a/764/CH4/EX4.20.a/data4_20.sci b/764/CH4/EX4.20.a/data4_20.sci new file mode 100755 index 000000000..1091703a4 --- /dev/null +++ b/764/CH4/EX4.20.a/data4_20.sci @@ -0,0 +1,9 @@ +
+//(Design against Static Load) Example 4.20
+//Refer Fig.4.67
+//Tensile yield strength of 30C8 Syt (N/mm2)
+Syt = 400
+//Factor of safety fs
+fs = 3.5
+//Load on the curved link P (kN)
+P = 1
diff --git a/764/CH4/EX4.20.b/result4_20.txt b/764/CH4/EX4.20.b/result4_20.txt new file mode 100755 index 000000000..413f733c4 --- /dev/null +++ b/764/CH4/EX4.20.b/result4_20.txt @@ -0,0 +1,28 @@ +-->//(Design against Static Load) Example 4.20
+
+-->//Refer Fig.4.67
+
+-->//Tensile yield strength of 30C8 Syt (N/mm2)
+
+-->Syt = 400
+ Syt =
+
+ 400.
+
+-->//Factor of safety fs
+
+-->fs = 3.5
+ fs =
+
+ 3.5
+
+-->//Load on the curved link P (kN)
+
+-->P = 1
+ P =
+
+ 1.
+
+
+The diameter of the link(D) = 20.105902 mm
+
\ No newline at end of file diff --git a/764/CH4/EX4.20.b/solution4_20.sce b/764/CH4/EX4.20.b/solution4_20.sce new file mode 100755 index 000000000..ddfa0d549 --- /dev/null +++ b/764/CH4/EX4.20.b/solution4_20.sce @@ -0,0 +1,31 @@ +
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_20.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_20.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the permissible stresses for steel parts sigmat (N/mm2)
+sigmat = Syt/fs
+//Assume the diameter of the link to be 1mm D
+D = 1
+R = 4 * D
+Ri = (4 * D)-(0.5 * D)
+Ro = (4 * D)+(0.5 * D)
+Rn = ((sqrt(Ro) + sqrt(Ri))^2)/4
+//Calculate the eccentricity e (mm)
+e = R - Rn
+hi = Rn - Ri
+A = (%pi/4)*(D^2)
+//Calculate the bending moment Mb (N-mm)
+Mb = (P * 1000)*(4 * D)
+//Calculate the bending stress B (N/mm2)
+B = (Mb * hi)/(A * e * Ri)
+//Calculate the direct tensile stress T (N/mm2)
+T = (P * 1000)/A
+//Calculate the true diameter of the link D (mm)
+D = ((B + T)/sigmat)^(1/2)
+//Print results
+printf('\nThe diameter of the link(D) = %f mm\n',D)
diff --git a/764/CH4/EX4.21.a/data4_21.sci b/764/CH4/EX4.21.a/data4_21.sci new file mode 100755 index 000000000..b238cca70 --- /dev/null +++ b/764/CH4/EX4.21.a/data4_21.sci @@ -0,0 +1,11 @@ +
+//(Design against Static Load) Example 4.21
+//Refer Fig.4.68 and 4.65
+//Load capacity of press P (kN)
+P = 100
+//Ultimate tensile strength of FG200 Sut (N/mm2)
+Sut = 200
+//Factor of safety fs
+fs = 3
+//Distance between force application point and C-Frame centre l (mm)
+l = 1000
diff --git a/764/CH4/EX4.21.b/result4_21.txt b/764/CH4/EX4.21.b/result4_21.txt new file mode 100755 index 000000000..db2ed2f0f --- /dev/null +++ b/764/CH4/EX4.21.b/result4_21.txt @@ -0,0 +1,35 @@ +-->//(Design against Static Load) Example 4.21
+
+-->//Refer Fig.4.68 and 4.65
+
+-->//Load capacity of press P (kN)
+
+-->P = 100
+ P =
+
+ 100.
+
+-->//Ultimate tensile strength of FG200 Sut (N/mm2)
+
+-->Sut = 200
+ Sut =
+
+ 200.
+
+-->//Factor of safety fs
+
+-->fs = 3
+ fs =
+
+ 3.
+
+-->//Distance between force application point and C-Frame centre l (mm)
+
+-->l = 1000
+ l =
+
+ 1000.
+
+
+The value of t = 100.000000 mm
+
\ No newline at end of file diff --git a/764/CH4/EX4.21.b/solution4_21.sce b/764/CH4/EX4.21.b/solution4_21.sce new file mode 100755 index 000000000..5181a08a8 --- /dev/null +++ b/764/CH4/EX4.21.b/solution4_21.sce @@ -0,0 +1,55 @@ +
+//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) then
+ v = v + (5 - rem)
+ end
+endfunction
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_21.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_21.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate permissible tensile stress sigma (N/mm2)
+sigma = Sut/fs
+//Assume value of t to be 1mm
+t = 1
+//All dimensions in mm
+bi = 3 * t
+h = 3 * t
+Ri = 2 * t
+Ro = 5 * t
+ti = t
+to = 0.75 * t
+Rn = ((ti * (bi - to)) + (to * h))/((bi - to)*log((Ri + ti)/Ri)+(to * log(Ro/Ri)))
+R = Ri + ((((1/2)*to*(h^2)) + ((1/2)*(ti^2)*(bi - to)))/((to*h) + (ti*(bi - to))))
+//Calculate eccentricity e (mm)
+e = R - Rn
+hi = Rn - Ri
+A = (bi * ti) + (to * Ri)
+//Calculate the direct tensile stress T (N/mm2)
+T = (P * 1000)/A
+//The polynomial obtained is as follows
+C = (((R * P * 1000 * hi)/(A * e * Ri)) + T)
+D = (P * 1000 * l * hi)/(A * e * Ri)
+p = [(-1 * sigma) 0 C D]
+r = roots(p)
+//Calculate the true value of t (mm)
+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('\nThe value of t = %f mm\n',t)
+
+
diff --git a/764/CH4/EX4.22.a/data4_22.sci b/764/CH4/EX4.22.a/data4_22.sci new file mode 100755 index 000000000..69717ce68 --- /dev/null +++ b/764/CH4/EX4.22.a/data4_22.sci @@ -0,0 +1,17 @@ +
+//(Design against Static Load) Example 4.22
+//Refer Fig.4.70 on page 136
+//Initial temperature of the steel tube Ti (degree)
+Ti = 25
+//Final temperature of the steel tube Tf (degree)
+Tf = 250
+//Length of the tube l (mm)
+l = 200
+//Area of cross-section of the tube A (mm2)
+A = 300
+//Joint separation j (mm)
+j = 0.15
+//Modulus of elasticity of steel E (N/mm2)
+E = 207000
+//Coefficient of thermal expansion of steel alpha (per degree celsius)
+alpha = 10.8 * (10^(-6))
diff --git a/764/CH4/EX4.22.b/result4_22.txt b/764/CH4/EX4.22.b/result4_22.txt new file mode 100755 index 000000000..5f5da5883 --- /dev/null +++ b/764/CH4/EX4.22.b/result4_22.txt @@ -0,0 +1,58 @@ +-->//(Design against Static Load) Example 4.22
+
+-->//Refer Fig.4.70 on page 136
+
+-->//Initial temperature of the steel tube Ti (degree)
+
+-->Ti = 25
+ Ti =
+
+ 25.
+
+-->//Final temperature of the steel tube Tf (degree)
+
+-->Tf = 250
+ Tf =
+
+ 250.
+
+-->//Length of the tube l (mm)
+
+-->l = 200
+ l =
+
+ 200.
+
+-->//Area of cross-section of the tube A (mm2)
+
+-->A = 300
+ A =
+
+ 300.
+
+-->//Joint separation j (mm)
+
+-->j = 0.15
+ j =
+
+ 0.15
+
+-->//Modulus of elasticity of steel E (N/mm2)
+
+-->E = 207000
+ E =
+
+ 207000.
+
+-->//Coefficient of thermal expansion of steel alpha (per degree celsius)
+
+-->alpha = 10.8 * (10^(-6))
+ alpha =
+
+ 0.0000108
+
+
+Force acting on the tube(P) = 104328.000000 N
+
+Resultant stress(sigma) = 347.760000 N/mm2
+
\ No newline at end of file diff --git a/764/CH4/EX4.22.b/solution4_22.sce b/764/CH4/EX4.22.b/solution4_22.sce new file mode 100755 index 000000000..208c1f522 --- /dev/null +++ b/764/CH4/EX4.22.b/solution4_22.sce @@ -0,0 +1,22 @@ +
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_22.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_22.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Temperature change Dt (degree)
+Dt = Tf - Ti
+//Calculate the expansion of tube delta (mm)
+delta = alpha * l * Dt
+//Net compression of tube c (mm)
+c = delta - j
+//Calculation of force P (N)
+P = (A * E * c)/l
+//Calculate the resultant stress sigma (N/mm2)
+sigma = P/A
+//Print results
+printf('\nForce acting on the tube(P) = %f N\n',P)
+printf('\nResultant stress(sigma) = %f N/mm2\n',sigma)
diff --git a/764/CH4/EX4.3.a/data4_3.sci b/764/CH4/EX4.3.a/data4_3.sci new file mode 100755 index 000000000..f2f580f44 --- /dev/null +++ b/764/CH4/EX4.3.a/data4_3.sci @@ -0,0 +1,12 @@ +
+//(Design against Static Load) Example 4.3
+//Tensile yield strength of steel 30C8 Syt (N/mm2)
+Syt = 400
+//Factor of safety fs
+fs = 4
+//Tensile force acting on the rods P (kN)
+P = 50
+//Inside diameter of the socket d2 (mm)
+d2 = 50
+//Outside diameter of the socket d4 (mm)
+d4 = 100
diff --git a/764/CH4/EX4.3.b/result4_3.txt b/764/CH4/EX4.3.b/result4_3.txt new file mode 100755 index 000000000..5d1cfdedd --- /dev/null +++ b/764/CH4/EX4.3.b/result4_3.txt @@ -0,0 +1,50 @@ +-->//(Design against Static Load) Example 4.3
+
+-->//Tensile yield strength of steel 30C8 Syt (N/mm2)
+
+-->Syt = 400
+ Syt =
+
+ 400.
+
+-->//Factor of safety fs
+
+-->fs = 4
+ fs =
+
+ 4.
+
+-->//Tensile force acting on the rods P (kN)
+
+-->P = 50
+ P =
+
+ 50.
+
+-->//Inside diameter of the socket d2 (mm)
+
+-->d2 = 50
+ d2 =
+
+ 50.
+
+-->//Outside diameter of the socket d4 (mm)
+
+-->d4 = 100
+ d4 =
+
+ 100.
+
+
+Shear failure criterion
+
+Thickness of cotter(t) = 10.000000 mm
+
+Width of cotter(b) = 50.000000 mm
+
+Bending failure criterion
+
+Thickness of cotter(t1) = 11.000000 mm
+
+Width of cotter(b1) = 55.000000 mm
+
\ No newline at end of file diff --git a/764/CH4/EX4.3.b/solution4_3.sce b/764/CH4/EX4.3.b/solution4_3.sce new file mode 100755 index 000000000..5f3fd332a --- /dev/null +++ b/764/CH4/EX4.3.b/solution4_3.sce @@ -0,0 +1,35 @@ +
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_3.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_3.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate permissible stresses for cotter (N/mm2)
+//Tensile stress sigma
+sigma = Syt/fs
+//Yield strength in shear for cotter Ssy (N/mm2)
+Ssy = (50/100)*Syt
+//Shear stress tau
+tau = Ssy/fs
+//Shear failure criterion
+//Calculate thickness of cotter t (mm)
+t = sqrt((P * 1000)/(2 * 5 * tau))
+//Calculate width of cotter b (mm)
+b = 5 * t
+//Bending failure criterion
+//Calculate thickness of cotter t1 (mm)
+t1 = ((((P * 1000)/2)*((d2/4) + ((d4 - d2)/6))*(5/2))/((100 * (5^3))/12))^(1/3)
+//Round up t1
+t1 = ceil(t1)
+//Calculate width of cotter b1 (mm)
+b1 = 5 * t1
+//Print results
+printf('\nShear failure criterion\n')
+printf('\nThickness of cotter(t) = %f mm\n',t)
+printf('\nWidth of cotter(b) = %f mm\n',b)
+printf('\nBending failure criterion\n')
+printf('\nThickness of cotter(t1) = %f mm\n',t1)
+printf('\nWidth of cotter(b1) = %f mm\n',b1)
\ No newline at end of file diff --git a/764/CH4/EX4.4.a/data4_4.sci b/764/CH4/EX4.4.a/data4_4.sci new file mode 100755 index 000000000..7bcc8d6ae --- /dev/null +++ b/764/CH4/EX4.4.a/data4_4.sci @@ -0,0 +1,10 @@ +
+//(Design against Static Load) Example 4.4
+//Tensile yield strength of plain carbon steel 40C8 Syt (N/mm2)
+Syt = 380
+//Diameter of each rod d (mm)
+d = 50
+//Thickness of cotter t (mm)
+t = 15
+//Factor of safety fs
+fs = 6
diff --git a/764/CH4/EX4.4.b/result4_4.txt b/764/CH4/EX4.4.b/result4_4.txt new file mode 100755 index 000000000..d288900e3 --- /dev/null +++ b/764/CH4/EX4.4.b/result4_4.txt @@ -0,0 +1,43 @@ +-->//(Design against Static Load) Example 4.4
+
+-->//Tensile yield strength of plain carbon steel 40C8 Syt (N/mm2)
+
+-->Syt = 380
+ Syt =
+
+ 380.
+
+-->//Diameter of each rod d (mm)
+
+-->d = 50
+ d =
+
+ 50.
+
+-->//Thickness of cotter t (mm)
+
+-->t = 15
+ t =
+
+ 15.
+
+-->//Factor of safety fs
+
+-->fs = 6
+ fs =
+
+ 6.
+
+
+Load acting on rods(P) = 124354.709205 N
+
+Inside diameter of socket(d2) = 65.000000 mm
+
+Outside diameter of socket(d1) = 85.000000 mm
+
+Diameter of socket collar(d4) = 135.000000 mm
+
+a = 35.000000 mm
+
+c = 30.000000 mm
+
\ No newline at end of file diff --git a/764/CH4/EX4.4.b/solution4_4.sce b/764/CH4/EX4.4.b/solution4_4.sce new file mode 100755 index 000000000..dd24a6e23 --- /dev/null +++ b/764/CH4/EX4.4.b/solution4_4.sce @@ -0,0 +1,73 @@ +
+//Function to calculate roots of a quadratic equation
+function[r] = quadratic(a,b,c)
+ //Calculate discriminant D
+ D = (b^2)-(4 * a * c)
+ r1 = ((-1 * b)+ sqrt(D))/(2 * a)
+ r2 = ((-1 * b)- sqrt(D))/(2 * a)
+ if(r1 > 0)
+ r = r1
+ else
+ r = r2
+ end
+endfunction
+
+//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_4.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_4.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Permissible stresses (N/mm2)
+//Compressive stress sigmac
+sigmac = (2 * Syt)/fs
+//Yield strength in shear of the material Ssy (N/mm2)
+Ssy = (50/100)*Syt
+//Shear stress tau
+tau = Ssy/fs
+//Tensile stress sigmat
+sigmat = Syt/fs
+//Calculate the load acting on the rods P (N)
+P = (%pi/4)*((d^2) * sigmat)
+//Calculate inside diameter of the socket d2 (mm)
+//a,b,c are the coefficients of the resulting quadratic equation
+//ax^2 + bx + c = 0 x=d2
+a = (%pi/4)*sigmat
+b = (-1)*(t * sigmat)
+c = (-1)*(P)
+d2 = quadratic(a, b, c)
+d2 = round_five(d2)
+//Calculate outside of socket d1 (mm)
+//a,b,c are the coefficients of the resulting quadratic equation
+//ax^2 + bx + c = 0 x=d1
+a = (%pi/4)*(sigmat)
+b = (-1)*(t * sigmat)
+c = (((-1 * (%pi/4) * (d2^2))+(t * d2))*sigmat)-(P)
+d1 = quadratic(a, b, c)
+d1 = round_five(d1)
+//Calculate daimeter of socket collar d4 (mm)
+d4 = (P/(sigmac * t)) + d2
+d4 = round_five(d4)
+//Calculate dimensions of a and c (mm)
+a = P/(2 * d2 * tau)
+a = round_five(a)
+c = P/(2 * (d4 - d2) * tau)
+c = round_five(c)
+//Print results
+printf('\nLoad acting on rods(P) = %f N\n',P)
+printf('\nInside diameter of socket(d2) = %f mm\n',d2)
+printf('\nOutside diameter of socket(d1) = %f mm\n',d1)
+printf('\nDiameter of socket collar(d4) = %f mm\n',d4)
+printf('\na = %f mm\n\nc = %f mm\n',a,c)
+
diff --git a/764/CH4/EX4.5.a/data4_5.sci b/764/CH4/EX4.5.a/data4_5.sci new file mode 100755 index 000000000..f40ea4781 --- /dev/null +++ b/764/CH4/EX4.5.a/data4_5.sci @@ -0,0 +1,8 @@ +
+//(Design against Static Load) Example 4.5
+//Tensile force acting on the rods P (N)
+P = 50000
+//Tensile yield strength for 30C8 material Syt (N/mm2)
+Syt = 400
+//Factor of safety fs
+fs = 5
diff --git a/764/CH4/EX4.5.b/result4_5.txt b/764/CH4/EX4.5.b/result4_5.txt new file mode 100755 index 000000000..5f7cffaa2 --- /dev/null +++ b/764/CH4/EX4.5.b/result4_5.txt @@ -0,0 +1,36 @@ +-->//(Design against Static Load) Example 4.5
+
+-->//Tensile force acting on the rods P (N)
+
+-->P = 50000
+ P =
+
+ 50000.
+
+-->//Tensile yield strength for 30C8 material Syt (N/mm2)
+
+-->Syt = 400
+ Syt =
+
+ 400.
+
+-->//Factor of safety fs
+
+-->fs = 5
+ fs =
+
+ 5.
+
+
+Diameter of rods(D) = 30.000000 mm
+
+Enlarged diameter of rods(D1) = 33.000000 mm
+
+Dimensions of a and b are 25.000000 mm and 40.000000 mm respectively
+
+Diameter of pin(d) = 40.000000 mm
+
+Dimensions of d0 and d1 are 80.000000 mm and 60.000000 mm respectively
+
+Design of knuckle joint is safe
+
\ No newline at end of file diff --git a/764/CH4/EX4.5.b/solution4_5.sce b/764/CH4/EX4.5.b/solution4_5.sce new file mode 100755 index 000000000..64726a76c --- /dev/null +++ b/764/CH4/EX4.5.b/solution4_5.sce @@ -0,0 +1,63 @@ +
+//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_5.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_5.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Refer Fig.4.24
+//Permissible stresses (N/mm2)
+//Tensile stress sigmat
+sigmat = Syt/fs
+//Compressive stress
+sigmac = Syt/fs
+//Yield strength in shear for the material Ssy (N/mm2)
+Ssy = (50/100)*Syt
+//Shear stress tau
+tau = Ssy/fs
+//Calculate diameter of the rods D (mm)
+D = sqrt((4 * P)/(%pi * sigmat))
+D = round_five(D)
+//Calculate enlarged diameter of rods D1 (mm)
+D1 = 1.1 * D
+//Calculate dimensions a and b (mm)
+a = 0.75 * D
+a = round_five(a)
+b = 1.25 * D
+b = round_five(b)
+//Calculate diameter of pin d
+d = ((32/(%pi * sigmat))*((P/2)*((b/4) + (a/3))))^(1/3)
+d = round_five(d)
+//Calculate dimensions d0 and d1
+d0 = 2 * d
+d1 = 1.5 * d
+//Print results
+printf('\nDiameter of rods(D) = %f mm\n',D)
+printf('\nEnlarged diameter of rods(D1) = %f mm\n',D1)
+printf('\nDimensions of a and b are %f mm and %f mm respectively\n',a,b)
+printf('\nDiameter of pin(d) = %f mm\n',d)
+printf('\nDimensions of d0 and d1 are %f mm and %f mm respectively\n',d0,d1)
+//Check for stresses in eye
+s1 = (P/(b * (d0 - d)))
+c1 = (P/(b * d))
+t1 = (P/(b * (d0 - d)))
+//Check for stresses in fork
+s2 = (P/(2 * a * (d0 - d)))
+c2 = (P/(2 * a * d))
+t2 = (P/(2 * a * (d0 - d)))
+if ((s1<sigmat) & (c1<sigmac) & (t1<tau) & (s2<sigmat) & (c2<sigmac) & (t2<tau))
+ printf('\nDesign of knuckle joint is safe\n')
+else
+ printf('\nDesign of knuckle joint is not safe\n')
+end
\ No newline at end of file diff --git a/764/CH4/EX4.6.a/data4_6.sci b/764/CH4/EX4.6.a/data4_6.sci new file mode 100755 index 000000000..4990b0207 --- /dev/null +++ b/764/CH4/EX4.6.a/data4_6.sci @@ -0,0 +1,13 @@ +
+//(Design against Static Load) Example 4.6
+//Refer Fig.4.25
+//Data according to standard table h (mm)
+h = {150,175,200,225,250}
+//Ixx (mm4)
+Ixx = {(688.2 * (10^4)), (1096.2 * (10^4)), (1696.6 * (10^4)), (2501.9 * (10^4)), (3717.8 * (10^4))}
+//Total weight of the bars W (kN)
+W = 75
+//Permissible bending stress sigmab (N/mm2)
+sigmab = 165
+//Length of cantilever beam l (mm)
+l = 2000
diff --git a/764/CH4/EX4.6.b/result4_6.txt b/764/CH4/EX4.6.b/result4_6.txt new file mode 100755 index 000000000..53ed0da1c --- /dev/null +++ b/764/CH4/EX4.6.b/result4_6.txt @@ -0,0 +1,42 @@ +-->//(Design against Static Load) Example 4.6
+
+-->//Refer Fig.4.25
+
+-->//Data according to standard table h (mm)
+
+-->h = {150,175,200,225,250}
+ h =
+
+ 150. 175. 200. 225. 250.
+
+-->//Ixx (mm4)
+
+-->Ixx = {(688.2 * (10^4)), (1096.2 * (10^4)), (1696.6 * (10^4)), (2501.9 * (10^4)), (3717.8 * (10^4))}
+ Ixx =
+
+ 6882000. 10962000. 16966000. 25019000. 37178000.
+
+-->//Total weight of the bars W (kN)
+
+-->W = 75
+ W =
+
+ 75.
+
+-->//Permissible bending stress sigmab (N/mm2)
+
+-->sigmab = 165
+ sigmab =
+
+ 165.
+
+-->//Length of cantilever beam l (mm)
+
+-->l = 2000
+ l =
+
+ 2000.
+
+
+ISLB200 is suitable for this application
+
\ No newline at end of file diff --git a/764/CH4/EX4.6.b/solution4_6.sce b/764/CH4/EX4.6.b/solution4_6.sce new file mode 100755 index 000000000..e25066af7 --- /dev/null +++ b/764/CH4/EX4.6.b/solution4_6.sce @@ -0,0 +1,44 @@ +
+//Function to print the designation
+function[] = desig(z)
+ printf('\nISLB%d is suitable for this application\n',z)
+endfunction
+
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_6.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_6.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Load supported by each beam P (kN)
+P = (W/2)
+//Distance of centre of gravity of the load from the rigid support d (mm)
+d = (l/3)
+//Calculate the bending moment Mb (N-mm)
+Mb = (P * 1000)*d
+//Calculation of ratio Ixx/y (mm3)
+ratio = Mb/sigmab
+//Selection of beam
+//There are 5 standard beams provided
+for i = 1:1:5
+ x = Ixx(i)/(h(i)/2)
+//Break the loop if value of x exceeds value of ratio
+ if (x > ratio)
+ break
+ end
+end
+//Print the suitable beam designation obtained
+if (i == 1)
+ desig(150)
+elseif (i == 2)
+ desig(175)
+elseif (i == 3)
+ desig(200)
+elseif (i == 4)
+ desig(250)
+else
+ desig(250)
+end
diff --git a/764/CH4/EX4.7.a/data4_7.sci b/764/CH4/EX4.7.a/data4_7.sci new file mode 100755 index 000000000..236870223 --- /dev/null +++ b/764/CH4/EX4.7.a/data4_7.sci @@ -0,0 +1,13 @@ +
+//(Design against Static Load) Example 4.7
+//Refer Fig.4.26
+//Initial tension in the hacksaw blade P (N)
+P = 300
+//Tensile yield strength of 30C8 Syt (N/mm2)
+Syt = 400
+//Factor of safety (fs)
+fs = 2.5
+//Ratio of depth to width of cross-section ratio
+ratio = 3
+//Length of hacksaw frame l (mm)
+l = 200
diff --git a/764/CH4/EX4.7.b/result4_7.txt b/764/CH4/EX4.7.b/result4_7.txt new file mode 100755 index 000000000..0d16b4371 --- /dev/null +++ b/764/CH4/EX4.7.b/result4_7.txt @@ -0,0 +1,44 @@ +-->//(Design against Static Load) Example 4.7
+
+-->//Refer Fig.4.26
+
+-->//Initial tension in the hacksaw blade P (N)
+
+-->P = 300
+ P =
+
+ 300.
+
+-->//Tensile yield strength of 30C8 Syt (N/mm2)
+
+-->Syt = 400
+ Syt =
+
+ 400.
+
+-->//Factor of safety (fs)
+
+-->fs = 2.5
+ fs =
+
+ 2.5
+
+-->//Ratio of depth to width of cross-section ratio
+
+-->ratio = 3
+ ratio =
+
+ 3.
+
+-->//Length of hacksaw frame l (mm)
+
+-->l = 200
+ l =
+
+ 200.
+
+
+Value of t = 6.266535 mm
+
+Area of cross-section = (6.266535 x 18.799604) mm2
+
\ No newline at end of file diff --git a/764/CH4/EX4.7.b/solution4_7.sce b/764/CH4/EX4.7.b/solution4_7.sce new file mode 100755 index 000000000..a67b1f7bb --- /dev/null +++ b/764/CH4/EX4.7.b/solution4_7.sce @@ -0,0 +1,47 @@ +
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_7.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_7.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate permissible tensile stress sigmat (N/mm2)
+sigmat = Syt/fs
+//Assume the wudth of the cross-section to be 1mm t
+t = 1
+//Calculate direct compressive stress sigmac (N/mm2)
+sigmac = P/(t * (ratio * t))
+//Calculate maximum bending moment Mb (N-mm)
+Mb = P * l
+//Calculate y
+y = 1.5 * t
+//Calculate the second moment of area I (mm4)
+I = (t * ((ratio * t)^3))/12
+//Calculate tensile bending stress at the lower fibre sigmab(N/mm2)
+sigmab = (Mb * y)/I
+//Finding the real value of width t (mm)
+//On superimposing the stress values, a cubic equation in t is obtained
+a = 0
+b = (sigmac/sigmat)
+c = (-1 * (sigmab/sigmat))
+//Define polynomial
+p = [1,a,b,c]
+//Calculate roots of this polynomial
+r = roots(p)
+real_part = real(r)
+for i = 1:1:3
+ if(real_part(i)>0)
+ t = real_part(i)
+ break
+ end
+end
+//Print results
+printf('\nValue of t = %f mm\n',t)
+printf('\nArea of cross-section = (%f x %f) mm2\n',t,(ratio * t))
+
+
+
+
+
diff --git a/764/CH4/EX4.8.a/data4_8.sci b/764/CH4/EX4.8.a/data4_8.sci new file mode 100755 index 000000000..fcf36dd27 --- /dev/null +++ b/764/CH4/EX4.8.a/data4_8.sci @@ -0,0 +1,14 @@ +
+//(Design against Static Load) Example 4.8
+//Refer Fig.4.27
+//Force acting on offset link P (kN)
+P = 25
+//Ultimate tensile strength of FG300 Sut (N/mm2)
+Sut = 300
+//Factor of safety (fs)
+fs = 3
+//Force offset distance e (mm)
+e = 10
+//Ratio of depth to width of the cross-section ratio
+ratio = 2
+
diff --git a/764/CH4/EX4.8.b/result4_8.txt b/764/CH4/EX4.8.b/result4_8.txt new file mode 100755 index 000000000..945b162fe --- /dev/null +++ b/764/CH4/EX4.8.b/result4_8.txt @@ -0,0 +1,45 @@ +-->//(Design against Static Load) Example 4.8
+
+-->//Refer Fig.4.27
+
+-->//Force acting on offset link P (kN)
+
+-->P = 25
+ P =
+
+ 25.
+
+-->//Ultimate tensile strength of FG300 Sut (N/mm2)
+
+-->Sut = 300
+ Sut =
+
+ 300.
+
+-->//Factor of safety (fs)
+
+-->fs = 3
+ fs =
+
+ 3.
+
+-->//Force offset distance e (mm)
+
+-->e = 10
+ e =
+
+ 10.
+
+-->//Ratio of depth to width of the cross-section ratio
+
+-->ratio = 2
+ ratio =
+
+ 2.
+
+
+
+Value of t = 25.443742 mm
+
+Area of cross-section = (25.443742 x 50.887483) mm2
+
\ No newline at end of file diff --git a/764/CH4/EX4.8.b/solution4_8.sce b/764/CH4/EX4.8.b/solution4_8.sce new file mode 100755 index 000000000..72db15254 --- /dev/null +++ b/764/CH4/EX4.8.b/solution4_8.sce @@ -0,0 +1,36 @@ +
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_8.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_8.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the permissible tensile stress for the link sigmat (N/mm2)
+sigmat = Sut/fs
+//Assume the value of t to be 1mm
+t = 1
+//Calculate the direct tensile stress D (N/mm2)
+D = (P * 1000)/(t * (ratio * t))
+//Calculate the value of y (mm)
+y = t
+//Calculate the second moment of area I (mm4)
+I = (t * ((ratio * t)^3))/12
+//Calculate the dimensions of cross-section t
+a = (-1 * D) + ((-1 * P * 1000)/I)
+b = (-1 * P * 1000 * e)/I
+p = [sigmat 0 a b]
+r = roots(p)
+real_part = real(r)
+for i = 1:1:3
+ if(real_part(i)>0)
+ t = real_part(i)
+ break
+ end
+end
+//Print results
+printf('\nValue of t = %f mm\n',t)
+printf('\nArea of cross-section = (%f x %f) mm2\n',t,(ratio * t))
+
+
diff --git a/764/CH4/EX4.9.a/data4_9.sci b/764/CH4/EX4.9.a/data4_9.sci new file mode 100755 index 000000000..568fcfea7 --- /dev/null +++ b/764/CH4/EX4.9.a/data4_9.sci @@ -0,0 +1,13 @@ +
+//(Design against Static Load) Example 4.9
+//Refer Fig.4.28
+//Maximum force acting on the frame P (kN)
+P = 20
+//Tensile yield strength of 45C8 material Syt (N/mm2)
+Syt = 380
+//Factor of safety fs
+fs = 2.5
+//Width of the cross-section of the frame w (mm)
+w = 150
+//Distance between the press tool and the frame d (mm)
+d = 200
diff --git a/764/CH4/EX4.9.b/result4_9.txt b/764/CH4/EX4.9.b/result4_9.txt new file mode 100755 index 000000000..de734bb80 --- /dev/null +++ b/764/CH4/EX4.9.b/result4_9.txt @@ -0,0 +1,42 @@ +-->//(Design against Static Load) Example 4.9
+
+-->//Refer Fig.4.28
+
+-->//Maximum force acting on the frame P (kN)
+
+-->P = 20
+ P =
+
+ 20.
+
+-->//Tensile yield strength of 45C8 material Syt (N/mm2)
+
+-->Syt = 380
+ Syt =
+
+ 380.
+
+-->//Factor of safety fs
+
+-->fs = 2.5
+ fs =
+
+ 2.5
+
+-->//Width of the cross-section of the frame w (mm)
+
+-->w = 150
+ w =
+
+ 150.
+
+-->//Distance between the press tool and the frame d (mm)
+
+-->d = 200
+ d =
+
+ 200.
+
+
+The thickness of the plate(t) = 5.263158 mm
+
\ No newline at end of file diff --git a/764/CH4/EX4.9.b/solution4_9.sce b/764/CH4/EX4.9.b/solution4_9.sce new file mode 100755 index 000000000..625afd38c --- /dev/null +++ b/764/CH4/EX4.9.b/solution4_9.sce @@ -0,0 +1,30 @@ +
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_9.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_9.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the permissible stress for the plates sigmat (N/mm2)
+sigmat = Syt/fs
+//Force acting on each plate P (kN)
+P = P/2
+//Calculate the bending moment at the frame Mb (N-mm)
+Mb = (P * 1000) * (d + (w/2))
+//Assume the plate thickness to be 1mm t
+t = 1
+//Calculate the direct tensile stress D (N/mm2)
+D = (P * 1000)/(w * t)
+//Calculate the value of y (mm)
+y = (w/2)
+//Calculate the second moment of area I (mm4)
+I = (t * (w^3))/12
+//Calculate the bending stress B (N/mm2)
+B = (Mb * y)/I
+//Calculate the plate thickness t (mm)
+t = (D + B)/sigmat
+//Print results
+printf('\nThe thickness of the plate(t) = %f mm\n',t)
+
\ No newline at end of file |
