summaryrefslogtreecommitdiff
path: root/764/CH7
diff options
context:
space:
mode:
authorpriyanka2015-06-24 15:03:17 +0530
committerpriyanka2015-06-24 15:03:17 +0530
commitb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
treeab291cffc65280e58ac82470ba63fbcca7805165 /764/CH7
downloadScilab-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/CH7')
-rwxr-xr-x764/CH7/EX7.1.a/data7_1.sci9
-rwxr-xr-x764/CH7/EX7.1.b/result7_1.txt28
-rwxr-xr-x764/CH7/EX7.1.b/solution7_1.sce29
-rwxr-xr-x764/CH7/EX7.10.a/data7_10.sci17
-rwxr-xr-x764/CH7/EX7.10.b/result7_10.txt56
-rwxr-xr-x764/CH7/EX7.10.b/solution7_10.sce20
-rwxr-xr-x764/CH7/EX7.11.a/data7_11.sci15
-rwxr-xr-x764/CH7/EX7.11.b/result7_11.txt49
-rwxr-xr-x764/CH7/EX7.11.b/solution7_11.sce18
-rwxr-xr-x764/CH7/EX7.12.a/data7_12.sci15
-rwxr-xr-x764/CH7/EX7.12.b/result7_12.txt49
-rwxr-xr-x764/CH7/EX7.12.b/solution7_12.sce18
-rwxr-xr-x764/CH7/EX7.13.a/data7_13.sci19
-rwxr-xr-x764/CH7/EX7.13.b/result7_13.txt63
-rwxr-xr-x764/CH7/EX7.13.b/solution7_13.sce38
-rwxr-xr-x764/CH7/EX7.14.a/data7_14.sci23
-rwxr-xr-x764/CH7/EX7.14.b/result7_14.txt77
-rwxr-xr-x764/CH7/EX7.14.b/solution7_14.sce44
-rwxr-xr-x764/CH7/EX7.15.a/data7_15.sci21
-rwxr-xr-x764/CH7/EX7.15.b/result7_15.txt70
-rwxr-xr-x764/CH7/EX7.15.b/solution7_15.sce36
-rwxr-xr-x764/CH7/EX7.16.a/data7_16.sci17
-rwxr-xr-x764/CH7/EX7.16.b/result7_16.txt56
-rwxr-xr-x764/CH7/EX7.16.b/solution7_16.sce29
-rwxr-xr-x764/CH7/EX7.17.a/data7_17.sci15
-rwxr-xr-x764/CH7/EX7.17.b/result7_17.txt49
-rwxr-xr-x764/CH7/EX7.17.b/solution7_17.sce19
-rwxr-xr-x764/CH7/EX7.18.a/data7_18.sci30
-rwxr-xr-x764/CH7/EX7.18.b/graph7_18.pngbin0 -> 13942 bytes
-rwxr-xr-x764/CH7/EX7.18.b/result7_18.txt100
-rwxr-xr-x764/CH7/EX7.18.b/solution7_18.sce60
-rwxr-xr-x764/CH7/EX7.18.c/functions7_18.sci48
-rwxr-xr-x764/CH7/EX7.19.a/data7_19.sci27
-rwxr-xr-x764/CH7/EX7.19.b/graph7_19.pngbin0 -> 14660 bytes
-rwxr-xr-x764/CH7/EX7.19.b/result7_19.txt91
-rwxr-xr-x764/CH7/EX7.19.b/solution7_19.sce59
-rwxr-xr-x764/CH7/EX7.19.c/functions7_19.sci48
-rwxr-xr-x764/CH7/EX7.2.a/data7_2.sci9
-rwxr-xr-x764/CH7/EX7.2.b/result7_2.txt28
-rwxr-xr-x764/CH7/EX7.2.b/solution7_2.sce30
-rwxr-xr-x764/CH7/EX7.20.a/data7_20.sci35
-rwxr-xr-x764/CH7/EX7.20.b/graph7_20.pngbin0 -> 14216 bytes
-rwxr-xr-x764/CH7/EX7.20.b/result7_20.txt119
-rwxr-xr-x764/CH7/EX7.20.b/solution7_20.sce72
-rwxr-xr-x764/CH7/EX7.20.c/functions7_20.sci48
-rwxr-xr-x764/CH7/EX7.21.a/data7_21.sci31
-rwxr-xr-x764/CH7/EX7.21.b/graph7_21.pngbin0 -> 14461 bytes
-rwxr-xr-x764/CH7/EX7.21.b/result7_21.txt105
-rwxr-xr-x764/CH7/EX7.21.b/solution7_21.sce58
-rwxr-xr-x764/CH7/EX7.21.c/functions7_21.sci48
-rwxr-xr-x764/CH7/EX7.22.a/data7_22.sci35
-rwxr-xr-x764/CH7/EX7.22.b/graph7_22.pngbin0 -> 15049 bytes
-rwxr-xr-x764/CH7/EX7.22.b/result7_22.txt121
-rwxr-xr-x764/CH7/EX7.22.b/solution7_22.sce68
-rwxr-xr-x764/CH7/EX7.23.a/data7_23.sci29
-rwxr-xr-x764/CH7/EX7.23.b/graph7_23.pngbin0 -> 13791 bytes
-rwxr-xr-x764/CH7/EX7.23.b/result7_23.txt98
-rwxr-xr-x764/CH7/EX7.23.b/solution7_23.sce54
-rwxr-xr-x764/CH7/EX7.3.a/data7_3.sci17
-rwxr-xr-x764/CH7/EX7.3.b/result7_3.txt56
-rwxr-xr-x764/CH7/EX7.3.b/solution7_3.sce47
-rwxr-xr-x764/CH7/EX7.4.a/data7_4.sci17
-rwxr-xr-x764/CH7/EX7.4.b/result7_4.txt56
-rwxr-xr-x764/CH7/EX7.4.b/solution7_4.sce34
-rwxr-xr-x764/CH7/EX7.5.a/data7_5.sci17
-rwxr-xr-x764/CH7/EX7.5.b/result7_5.txt56
-rwxr-xr-x764/CH7/EX7.5.b/solution7_5.sce33
-rwxr-xr-x764/CH7/EX7.6.a/data7_6.sci15
-rwxr-xr-x764/CH7/EX7.6.b/result7_6.txt49
-rwxr-xr-x764/CH7/EX7.6.b/solution7_6.sce30
-rwxr-xr-x764/CH7/EX7.7.a/data7_7.sci17
-rwxr-xr-x764/CH7/EX7.7.b/result7_7.txt56
-rwxr-xr-x764/CH7/EX7.7.b/solution7_7.sce22
-rwxr-xr-x764/CH7/EX7.8.a/data7_8.sci15
-rwxr-xr-x764/CH7/EX7.8.b/result7_8.txt49
-rwxr-xr-x764/CH7/EX7.8.b/solution7_8.sce24
-rwxr-xr-x764/CH7/EX7.9.a/data7_9.sci15
-rwxr-xr-x764/CH7/EX7.9.b/result7_9.txt49
-rwxr-xr-x764/CH7/EX7.9.b/solution7_9.sce24
79 files changed, 3048 insertions, 0 deletions
diff --git a/764/CH7/EX7.1.a/data7_1.sci b/764/CH7/EX7.1.a/data7_1.sci
new file mode 100755
index 000000000..ed5eb7415
--- /dev/null
+++ b/764/CH7/EX7.1.a/data7_1.sci
@@ -0,0 +1,9 @@
+
+//(Threaded Joints) Example 7.1
+//Refer Fig.7.13 on page 232
+//Weight of electric motor P (kN)
+P = 10
+//Yield tensile strength of 30C8 Syt (N/mm2)
+Syt = 400
+//Factor of safety fs
+fs = 6
diff --git a/764/CH7/EX7.1.b/result7_1.txt b/764/CH7/EX7.1.b/result7_1.txt
new file mode 100755
index 000000000..196e3188c
--- /dev/null
+++ b/764/CH7/EX7.1.b/result7_1.txt
@@ -0,0 +1,28 @@
+-->//(Threaded Joints) Example 7.1
+
+-->//Refer Fig.7.13 on page 232
+
+-->//Weight of electric motor P (kN)
+
+-->P = 10
+ P =
+
+ 10.
+
+-->//Yield tensile strength of 30C8 Syt (N/mm2)
+
+-->Syt = 400
+ Syt =
+
+ 400.
+
+-->//Factor of safety fs
+
+-->fs = 6
+ fs =
+
+ 6.
+
+
+Standard size of the bolt is M20
+ \ No newline at end of file
diff --git a/764/CH7/EX7.1.b/solution7_1.sce b/764/CH7/EX7.1.b/solution7_1.sce
new file mode 100755
index 000000000..aedd4a6b7
--- /dev/null
+++ b/764/CH7/EX7.1.b/solution7_1.sce
@@ -0,0 +1,29 @@
+
+//Function to standardise the given bolt-size
+function[v] = standard(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('solution7_1.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_1.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate permissible tensile stress sigmat (N/mm2)
+sigmat = Syt/fs
+//Calculate the core diameter of the eye-bolt dc (mm)
+dc = ((4 * P * 1000)/(%pi * sigmat))^(1/2)
+//use equation 7.4
+//Calculate the nominal diameter d (mm)
+d = dc/0.8
+//Standardise the bolt size from Table 7.1
+d = standard(d)
+//Print results
+printf('\nStandard size of the bolt is M%d\n',d)
diff --git a/764/CH7/EX7.10.a/data7_10.sci b/764/CH7/EX7.10.a/data7_10.sci
new file mode 100755
index 000000000..08e304e02
--- /dev/null
+++ b/764/CH7/EX7.10.a/data7_10.sci
@@ -0,0 +1,17 @@
+
+//(Threaded Joints) Example 7.10
+//Refer Fig.7.25 on page 241
+//Number of bolts N
+N = 4
+//Tension in the slack side of the belt Tslack (kN)
+Tslack = 5
+//Tension in the tight side of the belt Ttight (kN)
+Ttight = 10
+//Distance between boltA and pt. C l1 (mm)
+l1 = 50
+//Distance between boltB and pt. C l2 (mm)
+l2 = 150
+//Eccentricity value l (mm)
+l = 200
+//Maximum permissible tensile stress sigmaMax (N/mm2)
+sigmaMax = 60
diff --git a/764/CH7/EX7.10.b/result7_10.txt b/764/CH7/EX7.10.b/result7_10.txt
new file mode 100755
index 000000000..ddfbe46f2
--- /dev/null
+++ b/764/CH7/EX7.10.b/result7_10.txt
@@ -0,0 +1,56 @@
+-->//(Threaded Joints) Example 7.10
+
+-->//Refer Fig.7.25 on page 241
+
+-->//Number of bolts N
+
+-->N = 4
+ N =
+
+ 4.
+
+-->//Tension in the slack side of the belt Tslack (kN)
+
+-->Tslack = 5
+ Tslack =
+
+ 5.
+
+-->//Tension in the tight side of the belt Ttight (kN)
+
+-->Ttight = 10
+ Ttight =
+
+ 10.
+
+-->//Distance between boltA and pt. C l1 (mm)
+
+-->l1 = 50
+ l1 =
+
+ 50.
+
+-->//Distance between boltB and pt. C l2 (mm)
+
+-->l2 = 150
+ l2 =
+
+ 150.
+
+-->//Eccentricity value l (mm)
+
+-->l = 200
+ l =
+
+ 200.
+
+-->//Maximum permissible tensile stress sigmaMax (N/mm2)
+
+-->sigmaMax = 60
+ sigmaMax =
+
+ 60.
+
+
+Area at the core cross-section(A) = 150.000000 mm2
+ \ No newline at end of file
diff --git a/764/CH7/EX7.10.b/solution7_10.sce b/764/CH7/EX7.10.b/solution7_10.sce
new file mode 100755
index 000000000..953ed796a
--- /dev/null
+++ b/764/CH7/EX7.10.b/solution7_10.sce
@@ -0,0 +1,20 @@
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution7_10.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_10.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the resisting force in appropriate bolt P2 (N)
+if (l1 > l2) then
+ P2 = ((Ttight + Tslack)*1000 * l * l1)/(2*(l1^2 + l2^2))
+else
+ P2 = ((Ttight + Tslack)*1000 * l * l2)/(2*(l1^2 + l2^2))
+end
+//Calculate the core cross-section area of the bolt A (mm2)
+A = P2/sigmaMax
+//Choose proper diameter from Table 7.1
+//Print results
+printf('\nArea at the core cross-section(A) = %f mm2\n',A)
diff --git a/764/CH7/EX7.11.a/data7_11.sci b/764/CH7/EX7.11.a/data7_11.sci
new file mode 100755
index 000000000..6aa17b803
--- /dev/null
+++ b/764/CH7/EX7.11.a/data7_11.sci
@@ -0,0 +1,15 @@
+
+//(Threaded Joints) Example 7.11
+//Refer Fig.7.26 on page 243
+//Screw pitch circle diameter 2b (mm)
+b = 300/2
+//Flange diameter 2a (mm)
+a = 400/2
+//External force acting on the bearing P (kN)
+P = 25
+//Eccentricity value l (mm)
+l = 150
+//Maximum permissible tensile stress in the cap screw sigmaMax (N/mm2)
+sigmaMax = 50
+//Number of cap screws n
+n = 4
diff --git a/764/CH7/EX7.11.b/result7_11.txt b/764/CH7/EX7.11.b/result7_11.txt
new file mode 100755
index 000000000..05d6f959d
--- /dev/null
+++ b/764/CH7/EX7.11.b/result7_11.txt
@@ -0,0 +1,49 @@
+-->//(Threaded Joints) Example 7.11
+
+-->//Refer Fig.7.26 on page 243
+
+-->//Screw pitch circle diameter 2b (mm)
+
+-->b = 300/2
+ b =
+
+ 150.
+
+-->//Flange diameter 2a (mm)
+
+-->a = 400/2
+ a =
+
+ 200.
+
+-->//External force acting on the bearing P (kN)
+
+-->P = 25
+ P =
+
+ 25.
+
+-->//Eccentricity value l (mm)
+
+-->l = 150
+ l =
+
+ 150.
+
+-->//Maximum permissible tensile stress in the cap screw sigmaMax (N/mm2)
+
+-->sigmaMax = 50
+ sigmaMax =
+
+ 50.
+
+-->//Number of cap screws n
+
+-->n = 4
+ n =
+
+ 4.
+
+
+The nominal diameter of the cap screw(d) = 15.000000 mm
+ \ No newline at end of file
diff --git a/764/CH7/EX7.11.b/solution7_11.sce b/764/CH7/EX7.11.b/solution7_11.sce
new file mode 100755
index 000000000..f71362b43
--- /dev/null
+++ b/764/CH7/EX7.11.b/solution7_11.sce
@@ -0,0 +1,18 @@
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution7_11.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_11.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the tensile force in the cap screw P1 (N)
+P1 = (2 * P * 1000 * l * (a + (b * cosd(180/n))))/(n * ((2 * (a^2)) + (b^2)))
+//Calculate the core diameter of the cap screw dc (mm)
+dc = ((P1 * 4)/(%pi * sigmaMax))^(1/2)
+//Calculate the nominal diameter of the cap screw d (mm)
+d = dc/0.8
+d = ceil(d)
+//Print results
+printf('\nThe nominal diameter of the cap screw(d) = %f mm\n',d)
diff --git a/764/CH7/EX7.12.a/data7_12.sci b/764/CH7/EX7.12.a/data7_12.sci
new file mode 100755
index 000000000..fe91b795a
--- /dev/null
+++ b/764/CH7/EX7.12.a/data7_12.sci
@@ -0,0 +1,15 @@
+
+//(Threaded Joints) Example 7.12
+//Refer Fig.7.27 on page 244
+//Number of bolts n
+n = 16
+//Pitch circle diameter of the bolts 2b (mm)
+b = 2000/2
+//Diameter of the pillar flange 2a (mm)
+a = 2250/2
+//Load acting on the crane P (kN)
+P = 50
+//Radius r (mm)
+r = 7500
+//Maximum permissible tensile stress sigmaMax (N/mm2)
+sigmaMax = 75
diff --git a/764/CH7/EX7.12.b/result7_12.txt b/764/CH7/EX7.12.b/result7_12.txt
new file mode 100755
index 000000000..232d53208
--- /dev/null
+++ b/764/CH7/EX7.12.b/result7_12.txt
@@ -0,0 +1,49 @@
+-->//(Threaded Joints) Example 7.12
+
+-->//Refer Fig.7.27 on page 244
+
+-->//Number of bolts n
+
+-->n = 16
+ n =
+
+ 16.
+
+-->//Pitch circle diameter of the bolts 2b (mm)
+
+-->b = 2000/2
+ b =
+
+ 1000.
+
+-->//Diameter of the pillar flange 2a (mm)
+
+-->a = 2250/2
+ a =
+
+ 1125.
+
+-->//Load acting on the crane P (kN)
+
+-->P = 50
+ P =
+
+ 50.
+
+-->//Radius r (mm)
+
+-->r = 7500
+ r =
+
+ 7500.
+
+-->//Maximum permissible tensile stress sigmaMax (N/mm2)
+
+-->sigmaMax = 75
+ sigmaMax =
+
+ 75.
+
+
+Area at the core cross-section(A) = 319.690265 mm2
+ \ No newline at end of file
diff --git a/764/CH7/EX7.12.b/solution7_12.sce b/764/CH7/EX7.12.b/solution7_12.sce
new file mode 100755
index 000000000..556f5c345
--- /dev/null
+++ b/764/CH7/EX7.12.b/solution7_12.sce
@@ -0,0 +1,18 @@
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution7_12.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_12.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the moment arm for load P l (mm)
+l = (r - a)
+//Calculate the absolute maximum force Pmax (N)
+Pmax = (2 * P * 1000 * l * (a + b))/(n * ((2 * (a^2)) + (b^2)))
+//Calculate the core cross-section area of the cap screw A (mm2)
+A = Pmax/sigmaMax
+//Choose proper diameter from Table 7.1
+//Print results
+printf('\nArea at the core cross-section(A) = %f mm2\n',A)
diff --git a/764/CH7/EX7.13.a/data7_13.sci b/764/CH7/EX7.13.a/data7_13.sci
new file mode 100755
index 000000000..983bc55e2
--- /dev/null
+++ b/764/CH7/EX7.13.a/data7_13.sci
@@ -0,0 +1,19 @@
+
+//(Threaded Joints) Example 7.13
+//Refer Fig.7.28 on page 245
+//Maximum force in the tie rod P (kN)
+P = 5
+//Tie rod inclination with the horizontal theta (degree)
+theta = 30
+//Number of bolts N
+N = 4
+//Tensile yield strength of 30C8 Syt (N/mm2)
+Syt = 400
+//Factor of safety fs
+fs = 5
+//Distance between boltA and pt.C l1 (mm)
+l1 = 25
+//Distance between boltB and pt.C l2 (mm)
+l2 = 150 + 25
+//Eccentricity value e (mm)
+e = 50 + 25
diff --git a/764/CH7/EX7.13.b/result7_13.txt b/764/CH7/EX7.13.b/result7_13.txt
new file mode 100755
index 000000000..4f97ec61c
--- /dev/null
+++ b/764/CH7/EX7.13.b/result7_13.txt
@@ -0,0 +1,63 @@
+-->//(Threaded Joints) Example 7.13
+
+-->//Refer Fig.7.28 on page 245
+
+-->//Maximum force in the tie rod P (kN)
+
+-->P = 5
+ P =
+
+ 5.
+
+-->//Tie rod inclination with the horizontal theta (degree)
+
+-->theta = 30
+ theta =
+
+ 30.
+
+-->//Number of bolts N
+
+-->N = 4
+ N =
+
+ 4.
+
+-->//Tensile yield strength of 30C8 Syt (N/mm2)
+
+-->Syt = 400
+ Syt =
+
+ 400.
+
+-->//Factor of safety fs
+
+-->fs = 5
+ fs =
+
+ 5.
+
+-->//Distance between boltA and pt.C l1 (mm)
+
+-->l1 = 25
+ l1 =
+
+ 25.
+
+-->//Distance between boltB and pt.C l2 (mm)
+
+-->l2 = 150 + 25
+ l2 =
+
+ 175.
+
+-->//Eccentricity value e (mm)
+
+-->e = 50 + 25
+ e =
+
+ 75.
+
+
+Area at the core cross-section(A) = 33.170154 mm2
+ \ No newline at end of file
diff --git a/764/CH7/EX7.13.b/solution7_13.sce b/764/CH7/EX7.13.b/solution7_13.sce
new file mode 100755
index 000000000..5d7fc8894
--- /dev/null
+++ b/764/CH7/EX7.13.b/solution7_13.sce
@@ -0,0 +1,38 @@
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution7_13.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_13.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the permissible shear stress tauMax (N/mm2)
+tauMax = ((50/100)*Syt)/fs
+//Calculate vertical component of the axial force in the tie rod Pv (N)
+Pv = P * 1000 * sind(theta)
+//Calculate horizontal component of the axial force in the tie rod Ph (N)
+Ph = P * 1000 * cosd(theta)
+//Calculate the direct tensile force on the bolts Dtensile (N)
+Dtensile = Pv/N
+//Calculate the direct shear force on the bolts Sshear (N)
+Sshear = Ph/N
+//Calculate the tensile force on appropriate bolt due to bending moment Ftensile (N)
+if (l1 > l2) then
+ Ftensile = (Ph * e * l1)/(2*((l1^2) + (l2^2)))
+else
+ Ftensile = (Ph * e * l2)/(2*((l1^2) + (l2^2)))
+end
+//Assume the core cross-section area of the bolts to be 1mm2 A
+A = 1
+//Calculate the resultant tensile stress on the bolts res (N/mm2)
+res = (Dtensile + Ftensile)/A
+//Calculate the shear stress in bolts Stau (N/mm2)
+Stau = Sshear/A
+//Calculate the maximum shear stress in the bolts tau (N/mm2)
+tau = (((res/2)^2) + (Stau^2))^(1/2)
+//Calculate the actual core cross-section area of the bolts A (mm2)
+A = tau/tauMax
+//Choose proper diameter from Table 7.1
+//Print results
+printf('\nArea at the core cross-section(A) = %f mm2\n',A)
diff --git a/764/CH7/EX7.14.a/data7_14.sci b/764/CH7/EX7.14.a/data7_14.sci
new file mode 100755
index 000000000..7fcd6fd76
--- /dev/null
+++ b/764/CH7/EX7.14.a/data7_14.sci
@@ -0,0 +1,23 @@
+
+//(Threaded Joints) Example 7.14
+//Refer Fig.7.29 on page 246
+//Force acting on bracket P (kN)
+P = 5
+//Angle made by the force with the vertical theta (degree)
+theta = 60
+//Number of bolts N
+N = 4
+//Tensile yield strength of 30C8 Syt (N/mm2)
+Syt = 400
+//Factor of safety fs
+fs = 5
+//Distance of the centre of gravity from the lower edge cg (mm)
+cg = 60 + 100
+//Distance between the point force and lower edge f (mm)
+f = 200
+//Distance of the force point from the tilting edge t (mm)
+t = 240
+//Distance between bolt1 and pt.C l1 (mm)
+l1 = 200 + 60
+//Distance between bolt2 and pt.C l2 (mm)
+l2 = 60
diff --git a/764/CH7/EX7.14.b/result7_14.txt b/764/CH7/EX7.14.b/result7_14.txt
new file mode 100755
index 000000000..3b6ab9f51
--- /dev/null
+++ b/764/CH7/EX7.14.b/result7_14.txt
@@ -0,0 +1,77 @@
+-->//(Threaded Joints) Example 7.14
+
+-->//Refer Fig.7.29 on page 246
+
+-->//Force acting on bracket P (kN)
+
+-->P = 5
+ P =
+
+ 5.
+
+-->//Angle made by the force with the vertical theta (degree)
+
+-->theta = 60
+ theta =
+
+ 60.
+
+-->//Number of bolts N
+
+-->N = 4
+ N =
+
+ 4.
+
+-->//Tensile yield strength of 30C8 Syt (N/mm2)
+
+-->Syt = 400
+ Syt =
+
+ 400.
+
+-->//Factor of safety fs
+
+-->fs = 5
+ fs =
+
+ 5.
+
+-->//Distance of the centre of gravity from the lower edge cg (mm)
+
+-->cg = 60 + 100
+ cg =
+
+ 160.
+
+-->//Distance between the point force and lower edge f (mm)
+
+-->f = 200
+ f =
+
+ 200.
+
+-->//Distance of the force point from the tilting edge t (mm)
+
+-->t = 240
+ t =
+
+ 240.
+
+-->//Distance between bolt1 and pt.C l1 (mm)
+
+-->l1 = 200 + 60
+ l1 =
+
+ 260.
+
+-->//Distance between bolt2 and pt.C l2 (mm)
+
+-->l2 = 60
+ l2 =
+
+ 60.
+
+
+Area at the core cross-section(A) = 34.874654 mm2
+ \ No newline at end of file
diff --git a/764/CH7/EX7.14.b/solution7_14.sce b/764/CH7/EX7.14.b/solution7_14.sce
new file mode 100755
index 000000000..dca3fc713
--- /dev/null
+++ b/764/CH7/EX7.14.b/solution7_14.sce
@@ -0,0 +1,44 @@
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution7_14.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_14.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the permissible shear stress tauMax (N/mm2)
+tauMax = ((50/100)*Syt)/fs
+//Horizontal component of force P Ph (N)
+Ph = P * 1000 * sind(theta)
+//Vertical component of force P Pv (N)
+Pv = P * 1000 * cosd(theta)
+//Calculate the direct tensile force on each bolt Dtensile (N)
+Dtensile = Ph/N
+//Calculate the turning moment due to horizontal component Mh (N-mm)
+Mh = Ph * (f - cg)
+//Calculate the direct shear force on each bolt Sshear (N)
+Sshear = Pv/N
+//Calculate the turning moment due to vertical component Mv (N-mm)
+Mv = Pv * t
+//Calculate the tensile force on appropriate bolt due to bending moment Ftensile (N)
+if (l1 > l2) then
+ Ftensile = ((Mh + Mv)*l1)/(2*((l1^2) + (l2^2)))
+else
+ Ftensile = ((Mh + Mv)*l2)/(2*((l1^2) + (l2^2)))
+end
+//Calculate the total tensile force on each bolt Pt (N)
+Pt = Dtensile + Ftensile
+//Assume the core cross-section area of the bolts to be 1mm2 A
+A = 1
+//Calculate the resultant tensile stress in the bolt res (N/mm2)
+res = Pt/A
+//Calculate the shear stress in bolts Stau (N/mm2)
+Stau = Sshear/A
+//Calculate the maximum shear stress in the bolts tau (N/mm2)
+tau = (((res/2)^2) + (Stau^2))^(1/2)
+//Calculate the actual core cross-section area of the bolts A (mm2)
+A = tau/tauMax
+//Choose proper diameter from Table 7.1
+//Print results
+printf('\nArea at the core cross-section(A) = %f mm2\n',A)
diff --git a/764/CH7/EX7.15.a/data7_15.sci b/764/CH7/EX7.15.a/data7_15.sci
new file mode 100755
index 000000000..8c703c035
--- /dev/null
+++ b/764/CH7/EX7.15.a/data7_15.sci
@@ -0,0 +1,21 @@
+
+//(Threaded Joints) Example 7.15
+//Refer Fig.7.30 on page 247
+//Vertical force acting on the rigid bracket P (kN)
+P = 10
+//Maximum permissible shear stress in any bolt tauMax (N/mm2)
+tauMax = 50
+//Number of bolts N
+N = 4
+//Eccentricity value when secondary shear is considered e (mm)
+e = 250
+//Distance between two bolts dist (mm)
+dist = 100 + 100
+//Angle made by the secondary shear forces with the vertical theta (degree)
+theta = 45
+//Eccentricity value when bracket tilting is considered et (mm)
+et = 300
+//Distance between bolts 1 and 3 from edge CC l1 (mm)
+l1 = 50 + 200
+//Distance between bolts 2 and 4 from edge CC l2 (mm)
+l2 = 50
diff --git a/764/CH7/EX7.15.b/result7_15.txt b/764/CH7/EX7.15.b/result7_15.txt
new file mode 100755
index 000000000..257e970c2
--- /dev/null
+++ b/764/CH7/EX7.15.b/result7_15.txt
@@ -0,0 +1,70 @@
+-->//(Threaded Joints) Example 7.15
+
+-->//Refer Fig.7.30 on page 247
+
+-->//Vertical force acting on the rigid bracket P (kN)
+
+-->P = 10
+ P =
+
+ 10.
+
+-->//Maximum permissible shear stress in any bolt tauMax (N/mm2)
+
+-->tauMax = 50
+ tauMax =
+
+ 50.
+
+-->//Number of bolts N
+
+-->N = 4
+ N =
+
+ 4.
+
+-->//Eccentricity value when secondary shear is considered e (mm)
+
+-->e = 250
+ e =
+
+ 250.
+
+-->//Distance between two bolts dist (mm)
+
+-->dist = 100 + 100
+ dist =
+
+ 200.
+
+-->//Angle made by the secondary shear forces with the vertical theta (degree)
+
+-->theta = 45
+ theta =
+
+ 45.
+
+-->//Eccentricity value when bracket tilting is considered et (mm)
+
+-->et = 300
+ et =
+
+ 300.
+
+-->//Distance between bolts 1 and 3 from edge CC l1 (mm)
+
+-->l1 = 50 + 200
+ l1 =
+
+ 250.
+
+-->//Distance between bolts 2 and 4 from edge CC l2 (mm)
+
+-->l2 = 50
+ l2 =
+
+ 50.
+
+
+Area at the core cross-section(A) = 141.035110 mm2
+ \ No newline at end of file
diff --git a/764/CH7/EX7.15.b/solution7_15.sce b/764/CH7/EX7.15.b/solution7_15.sce
new file mode 100755
index 000000000..17ac33c34
--- /dev/null
+++ b/764/CH7/EX7.15.b/solution7_15.sce
@@ -0,0 +1,36 @@
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution7_15.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_15.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the primary shear force on bolts Pshear (N)
+Pshear = (P * 1000)/N
+//Distance between bolt and centre of gravity r (mm)
+r = (((dist/2)^2) + ((dist/2)^2))^(1/2)
+//Calculate the secondary shear force on bolts Sshear (N)
+Sshear = (P * 1000 * e * r)/(4 * (r^2))
+//Calculate the resultant shear force on bolts Ps (N)
+Ps = (((Pshear + (Sshear * sind(theta)))^2) + ((Sshear * cosd(90 - theta))^2))^(1/2)
+//Calculate the resisting force set up in appropriate bolts due to bracket tilting Pt (N)
+if (l1 > l2) then
+ Pt = (P * 1000 * et * l1)/(2*((l1^2) + (l2^2)))
+else
+ Pt = (P * 1000 * et * l2)/(2*((l1^2) + (l2^2)))
+end
+//Assume the core cross-section area of the bolts to be 1mm2 A
+A = 1
+//Calculate the tensile stress sigmat (N/mm2)
+sigmat = Pt/A
+//Calculate the shear stress Stau (N/mm2)
+Stau = Ps/A
+//Calculate the maximum shear stress in the bolts tau (N/mm2)
+tau = (((sigmat/2)^2) + (Stau^2))^(1/2)
+//Calculate the actual core cross-section area of the bolts A (mm2)
+A = tau/tauMax
+//Choose proper diameter from Table 7.1
+//Print results
+printf('\nArea at the core cross-section(A) = %f mm2\n',A)
diff --git a/764/CH7/EX7.16.a/data7_16.sci b/764/CH7/EX7.16.a/data7_16.sci
new file mode 100755
index 000000000..e66aee17f
--- /dev/null
+++ b/764/CH7/EX7.16.a/data7_16.sci
@@ -0,0 +1,17 @@
+
+//(Threaded Joints) Example 7.16
+//Refer Fig.7.35 on page 256
+//Yield tensile strength of 30C8 Syt (N/mm2)
+Syt = 380
+//Young's modulus of 30C8 E1 (N/mm2)
+E1 = 207000
+//Young's modulus of aluminium E2 (N/mm2)
+E2 = 71000
+//Initial pre-load in the bolt Pi (kN)
+Pi = 5
+//External force acting on the bolted joint P (kN)
+P = 10
+//Factor of safety fs
+fs = 2.5
+//Length of the joint l (mm)
+l = 25 + 25
diff --git a/764/CH7/EX7.16.b/result7_16.txt b/764/CH7/EX7.16.b/result7_16.txt
new file mode 100755
index 000000000..72d73e10c
--- /dev/null
+++ b/764/CH7/EX7.16.b/result7_16.txt
@@ -0,0 +1,56 @@
+-->//(Threaded Joints) Example 7.16
+
+-->//Refer Fig.7.35 on page 256
+
+-->//Yield tensile strength of 30C8 Syt (N/mm2)
+
+-->Syt = 380
+ Syt =
+
+ 380.
+
+-->//Young's modulus of 30C8 E1 (N/mm2)
+
+-->E1 = 207000
+ E1 =
+
+ 207000.
+
+-->//Young's modulus of aluminium E2 (N/mm2)
+
+-->E2 = 71000
+ E2 =
+
+ 71000.
+
+-->//Initial pre-load in the bolt Pi (kN)
+
+-->Pi = 5
+ Pi =
+
+ 5.
+
+-->//External force acting on the bolted joint P (kN)
+
+-->P = 10
+ P =
+
+ 10.
+
+-->//Factor of safety fs
+
+-->fs = 2.5
+ fs =
+
+ 2.5
+
+-->//Length of the joint l (mm)
+
+-->l = 25 + 25
+ l =
+
+ 50.
+
+
+Area at the core cross-section(A) = 65.319549 mm2
+ \ No newline at end of file
diff --git a/764/CH7/EX7.16.b/solution7_16.sce b/764/CH7/EX7.16.b/solution7_16.sce
new file mode 100755
index 000000000..957a40614
--- /dev/null
+++ b/764/CH7/EX7.16.b/solution7_16.sce
@@ -0,0 +1,29 @@
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution7_16.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_16.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the permissible tensile stress sigmat (N/mm2)
+sigmat = Syt/fs
+//Assume the inner diameter of the circular plate to be 1mm di
+di = 1
+//Calculate the outer diameter of the circular plate do (mm)
+do = 2 * di
+//Calculate the stiffness of the bolts kb (N/mm)
+kb = (%pi/4*(di^2))*(E1/l)
+//Calculate the area of the two plates Ac (mm2)
+Ac = (%pi/4)*((do^2) - (di^2))
+//Calculate the combined stiffness of the two plates kc (N/mm)
+kc = (Ac * E2)/l
+//Calculate the resultant load on the bolt Pb (N)
+deltaP = (P * 1000)*(kb/(kb + kc))
+Pb = (Pi * 1000) + deltaP
+//Calculate the core cross-section area of the bolt A (mm2)
+A = Pb/sigmat
+//Choose proper diameter from Table 7.1
+//Print results
+printf('\nArea at the core cross-section(A) = %f mm2\n',A)
diff --git a/764/CH7/EX7.17.a/data7_17.sci b/764/CH7/EX7.17.a/data7_17.sci
new file mode 100755
index 000000000..69b2cd9c2
--- /dev/null
+++ b/764/CH7/EX7.17.a/data7_17.sci
@@ -0,0 +1,15 @@
+
+//(Threaded Joints) Example 7.17
+//Refer Fig. 7.36 on page 257
+//Pre-load in the bolt Pi (kN)
+Pi = 2.5
+//External force acting on the bolt P (kN)
+P = 5
+//Yield tensile strength of 30C8 Syt (N/mm2)
+Syt = 400
+//Factor of safety fs
+fs = 2.5
+//Assume stiffness of bolts to be 1N/mm kb
+kb = 1
+//Stiffness of parts held together kc (N/mm)
+kc = 2.5 * kb
diff --git a/764/CH7/EX7.17.b/result7_17.txt b/764/CH7/EX7.17.b/result7_17.txt
new file mode 100755
index 000000000..dad2f1cab
--- /dev/null
+++ b/764/CH7/EX7.17.b/result7_17.txt
@@ -0,0 +1,49 @@
+-->//(Threaded Joints) Example 7.17
+
+-->//Refer Fig. 7.36 on page 257
+
+-->//Pre-load in the bolt Pi (kN)
+
+-->Pi = 2.5
+ Pi =
+
+ 2.5
+
+-->//External force acting on the bolt P (kN)
+
+-->P = 5
+ P =
+
+ 5.
+
+-->//Yield tensile strength of 30C8 Syt (N/mm2)
+
+-->Syt = 400
+ Syt =
+
+ 400.
+
+-->//Factor of safety fs
+
+-->fs = 2.5
+ fs =
+
+ 2.5
+
+-->//Assume stiffness of bolts to be 1N/mm kb
+
+-->kb = 1
+ kb =
+
+ 1.
+
+-->//Stiffness of parts held together kc (N/mm)
+
+-->kc = 2.5 * kb
+ kc =
+
+ 2.5
+
+
+Tensile stress area of the bolt(A) = 24.553571 mm2
+ \ No newline at end of file
diff --git a/764/CH7/EX7.17.b/solution7_17.sce b/764/CH7/EX7.17.b/solution7_17.sce
new file mode 100755
index 000000000..59006bddc
--- /dev/null
+++ b/764/CH7/EX7.17.b/solution7_17.sce
@@ -0,0 +1,19 @@
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution7_17.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_17.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the permissible tensile stress sigmat (N/mm2)
+sigmat = Syt/fs
+//Calculate the resultant load on bolt Pb (N)
+deltaP = (P * 1000)*(kb/(kb + kc))
+Pb = (Pi * 1000) + deltaP
+//Calculate the tensile stress area of the bolt A (mm2)
+A = Pb/sigmat
+//Choose proper diameter from Table 7.1
+//Print results
+printf('\nTensile stress area of the bolt(A) = %f mm2\n',A)
diff --git a/764/CH7/EX7.18.a/data7_18.sci b/764/CH7/EX7.18.a/data7_18.sci
new file mode 100755
index 000000000..ff13e3369
--- /dev/null
+++ b/764/CH7/EX7.18.a/data7_18.sci
@@ -0,0 +1,30 @@
+
+//(Threaded Joints) Example 7.18
+//Refer Fig.7.35 on page 256
+//Maximum force acting on the assembly Pmax (kN)
+Pmax = 10
+//Minimum force acting on the assembly Pmin (kN)
+Pmin = 0
+//Yield tensile strength of 45C8 Syt (N/mm2)
+Syt = 380
+//Ultimate tensile strength of 45C8 Sut (N/mm2)
+Sut = 630
+//Young's modulus of the plain carbon steel E1 (N/mm2)
+E1 = 207000
+//Young's modulus of aluminium E2 (N/mm2)
+E2 = 71000
+//Fatigue stress concentration factor Kf
+Kf = 2.2
+//Expected reliability (%)
+reliability = 90
+//Initial pre-load in the bolt Pi (kN)
+Pi = 5
+//Factor of safety fs
+fs = 2
+//Length of the bolt l (mm)
+l = 25 + 25
+//Assume diameter of the bolt to be less than 7.5mm for Kb to be 1
+d = 1
+//As Ka is incorporated into Kf, assume the assembly to be ground
+//This is just to obtain a random value of Ka
+op = 1 \ No newline at end of file
diff --git a/764/CH7/EX7.18.b/graph7_18.png b/764/CH7/EX7.18.b/graph7_18.png
new file mode 100755
index 000000000..e13ca6307
--- /dev/null
+++ b/764/CH7/EX7.18.b/graph7_18.png
Binary files differ
diff --git a/764/CH7/EX7.18.b/result7_18.txt b/764/CH7/EX7.18.b/result7_18.txt
new file mode 100755
index 000000000..4b35b6c3a
--- /dev/null
+++ b/764/CH7/EX7.18.b/result7_18.txt
@@ -0,0 +1,100 @@
+-->//(Threaded Joints) Example 7.18
+
+-->//Refer Fig.7.35 on page 256
+
+-->//Maximum force acting on the assembly Pmax (kN)
+
+-->Pmax = 10
+ Pmax =
+
+ 10.
+
+-->//Minimum force acting on the assembly Pmin (kN)
+
+-->Pmin = 0
+ Pmin =
+
+ 0.
+
+-->//Yield tensile strength of 45C8 Syt (N/mm2)
+
+-->Syt = 380
+ Syt =
+
+ 380.
+
+-->//Ultimate tensile strength of 45C8 Sut (N/mm2)
+
+-->Sut = 630
+ Sut =
+
+ 630.
+
+-->//Young's modulus of the plain carbon steel E1 (N/mm2)
+
+-->E1 = 207000
+ E1 =
+
+ 207000.
+
+-->//Young's modulus of aluminium E2 (N/mm2)
+
+-->E2 = 71000
+ E2 =
+
+ 71000.
+
+-->//Fatigue stress concentration factor Kf
+
+-->Kf = 2.2
+ Kf =
+
+ 2.2
+
+-->//Expected reliability (%)
+
+-->reliability = 90
+ reliability =
+
+ 90.
+
+-->//Initial pre-load in the bolt Pi (kN)
+
+-->Pi = 5
+ Pi =
+
+ 5.
+
+-->//Factor of safety fs
+
+-->fs = 2
+ fs =
+
+ 2.
+
+-->//Length of the bolt l (mm)
+
+-->l = 25 + 25
+ l =
+
+ 50.
+
+-->//Assume diameter of the bolt to be less than 7.5mm for Kb to be 1
+
+-->d = 1
+ d =
+
+ 1.
+
+-->//As Ka is incorporated into Kf, assume the assembly to be ground
+
+-->//This is just to obtain a random value of Ka
+
+-->op = 1
+ op =
+
+ 1.
+
+
+Core cross-section area of the bolt(A) = 54.133961 mm2
+ \ No newline at end of file
diff --git a/764/CH7/EX7.18.b/solution7_18.sce b/764/CH7/EX7.18.b/solution7_18.sce
new file mode 100755
index 000000000..09221472e
--- /dev/null
+++ b/764/CH7/EX7.18.b/solution7_18.sce
@@ -0,0 +1,60 @@
+
+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('solution7_18.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_18.sci'
+//Obtain path of function file
+funcpath = path + filesep() + 'functions7_18.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+exec(funcpath,[-1])
+//Calculate the endurance limit stress for bolt Sdash (N/mm2)
+Sdash = (50/100)*Sut
+//Calculate Ka, Kb and Kc
+[Ka, Kb, Kc] = fluctuate(op, d, reliability)
+//Calculate Kd
+Kd = 1/Kf
+//Calculate the corrected endurance limit stress Se (N/mm2)
+Se = Kb * Kc * Kd * Sdash
+//Calculate the stiffness of bolt kb (N/mm)
+kb = ((%pi/4)*(d^2))*(E1/l)
+//Calculate the area of the two plates Ac (mm2)
+Ac = (%pi/4)*(((2*d)^2) - (d^2))
+//Calculate the stiffness of the plates kc (N/mm)
+kc = (Ac * E2)/l
+//Calculate the maximum force in the bolt PMAX (N)
+PMAX = (Pi * 1000) + ((kb/(kb + kc))*(Pmax * 1000))
+//Calculate the minimum force in the bolt PMIN (N)
+PMIN = (Pi * 1000) + ((kb/(kb + kc))*(Pmin * 1000))
+//Calculate the mean force and force amplitude
+Pm = (PMAX + PMIN)/2
+Pa = (PMAX - PMIN)/2
+//Plot modified Goodman diagram
+//The common quadrilateral in the plot is the area of concern
+y1 = {Se 0}
+x1 = {0 Sut}
+y2 = {Syt 0}
+x2 = {0 Syt}
+plot(x1,y1,'--*')
+plot(x2,y2,'-*')
+plot_format()
+title('Modified Goodman diagram (Example 7.18)')
+xlabel('sigmaM (N/mm2)')
+ylabel('sigmaA (N/mm2)')
+//Calculate the actual core cross-section area of the bolt A (mm2)
+A = (Pa + ((Pi * 1000)/((1 + (Sut/Se)) * fs)))/(Sut/((1 + (Sut/Se))*fs))
+//Choose proper diameter from Table 7.1
+//Print results
+printf('\nCore cross-section area of the bolt(A) = %f mm2\n',A)
diff --git a/764/CH7/EX7.18.c/functions7_18.sci b/764/CH7/EX7.18.c/functions7_18.sci
new file mode 100755
index 000000000..b53d32df5
--- /dev/null
+++ b/764/CH7/EX7.18.c/functions7_18.sci
@@ -0,0 +1,48 @@
+
+//Function generating the values of Ka, Kb and Kc
+function [Ka, Kb, Kc] = fluctuate(s, d, r)
+ //Calculate Ka
+ //Nomenclature:
+ //1 - Ground
+ //2 - Machined or cold drawn
+ //3 - Hot-rolled
+ //4 - Forged
+surface = [1 2 3 4]
+Ksurfa = [1.58 4.51 57.7 272]
+Ksurfb = [-0.085 -0.265 -0.718 -0.995]
+for j = 1:1:4
+ if (s == surface(j))
+ a = Ksurfa(j)
+ b = Ksurfb(j)
+ break
+ end
+end
+//From equation 5.18 on page 157
+Ka = a * (Sut^b)
+if (Ka > 1) then
+ Ka = 1
+end
+
+ //Calculate Kb
+ //d (mm)
+ if (d <= 7.5) then
+ Kb = 1
+ elseif ((d > 7.5) & (d <= 50))
+ Kb = 0.85
+ elseif (d > 50)
+ Kb = 0.75
+ else
+ printf('Error in Kb')
+ end
+
+ //Calculate Kc
+ // r (%)
+rel = [50 90 95 99 99.9 99.99 99.999]
+Krel = [1 0.897 0.868 0.814 0.753 0.702 0.659]
+for i = 1:1:7
+ if (r == rel(i)) then
+ Kc = Krel(i)
+ break
+ end
+end
+endfunction
diff --git a/764/CH7/EX7.19.a/data7_19.sci b/764/CH7/EX7.19.a/data7_19.sci
new file mode 100755
index 000000000..917efb6b0
--- /dev/null
+++ b/764/CH7/EX7.19.a/data7_19.sci
@@ -0,0 +1,27 @@
+
+//(Threaded Joints) Example 7.19
+//Maximum external force Pmax (kN)
+Pmax = 10
+//Minimum external force Pmin (kN)
+Pmin = 0
+//Assume the stiffness of the bolts to be 1N/mm kb
+kb = 1
+//Calculate the stiffness of the parts kc
+kc = 3 * kb
+//Overload percentage load (%)
+load = 50
+//Yield tensile strength of 50C4 Syt (N/mm2)
+Syt = 460
+//Ultimate tensile strength of 50C4 Sut (N/mm2)
+Sut = 660
+//Fatigue stress concentration factor Kf
+Kf = 2.2
+//Expected reliability (%)
+reliability = 90
+//Factor of safety fs
+fs = 2
+//Assume diameter of the bolt to be less than 7.5mm for Kb to be 1
+d = 1
+//As Ka is incorporated into Kf, assume the assembly to be ground
+//This is just to obtain a random value of Ka
+op = 1
diff --git a/764/CH7/EX7.19.b/graph7_19.png b/764/CH7/EX7.19.b/graph7_19.png
new file mode 100755
index 000000000..e7b7f8c0b
--- /dev/null
+++ b/764/CH7/EX7.19.b/graph7_19.png
Binary files differ
diff --git a/764/CH7/EX7.19.b/result7_19.txt b/764/CH7/EX7.19.b/result7_19.txt
new file mode 100755
index 000000000..d9a304656
--- /dev/null
+++ b/764/CH7/EX7.19.b/result7_19.txt
@@ -0,0 +1,91 @@
+-->//(Threaded Joints) Example 7.19
+
+-->//Maximum external force Pmax (kN)
+
+-->Pmax = 10
+ Pmax =
+
+ 10.
+
+-->//Minimum external force Pmin (kN)
+
+-->Pmin = 0
+ Pmin =
+
+ 0.
+
+-->//Assume the stiffness of the bolts to be 1N/mm kb
+
+-->kb = 1
+ kb =
+
+ 1.
+
+-->//Calculate the stiffness of the parts kc
+
+-->kc = 3 * kb
+ kc =
+
+ 3.
+
+-->//Overload percentage load (%)
+
+-->load = 50
+ load =
+
+ 50.
+
+-->//Yield tensile strength of 50C4 Syt (N/mm2)
+
+-->Syt = 460
+ Syt =
+
+ 460.
+
+-->//Ultimate tensile strength of 50C4 Sut (N/mm2)
+
+-->Sut = 660
+ Sut =
+
+ 660.
+
+-->//Fatigue stress concentration factor Kf
+
+-->Kf = 2.2
+ Kf =
+
+ 2.2
+
+-->//Expected reliability (%)
+
+-->reliability = 90
+ reliability =
+
+ 90.
+
+-->//Factor of safety fs
+
+-->fs = 2
+ fs =
+
+ 2.
+
+-->//Assume diameter of the bolt to be less than 7.5mm for Kb to be 1
+
+-->d = 1
+ d =
+
+ 1.
+
+-->//As Ka is incorporated into Kf, assume the assembly to be ground
+
+-->//This is just to obtain a random value of Ka
+
+-->op = 1
+ op =
+
+ 1.
+
+
+Core cross-section area of the bolt(A) = 39.413787 mm2
+ \ No newline at end of file
diff --git a/764/CH7/EX7.19.b/solution7_19.sce b/764/CH7/EX7.19.b/solution7_19.sce
new file mode 100755
index 000000000..bd4260cf2
--- /dev/null
+++ b/764/CH7/EX7.19.b/solution7_19.sce
@@ -0,0 +1,59 @@
+
+
+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('solution7_19.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_19.sci'
+//Obtain path of function file
+funcpath = path + filesep() + 'functions7_19.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+exec(funcpath,[-1])
+//Calculate the endurance limit stress for bolt Sdash (N/mm2)
+Sdash = (50/100)*Sut
+//Calculate Ka, Kb and Kc
+[Ka, Kb, Kc] = fluctuate(op, d, reliability)
+//Calculate Kd
+Kd = 1/Kf
+//Calculate the corrected endurance limit stress Se (N/mm2)
+Se = Kb * Kc * Kd * Sdash
+//Plot modified Goodman diagram
+//The common quadrilateral in the plot is the area of concern
+y1 = {Se 0}
+x1 = {0 Sut}
+y2 = {Syt 0}
+x2 = {0 Syt}
+plot(x1,y1,'--*')
+plot(x2,y2,'-*')
+plot_format()
+title('Modified Goodman diagram (Example 7.19)')
+xlabel('sigmaM (N/mm2)')
+ylabel('sigmaA (N/mm2)')
+//Calculate the external force at overload condition PbMax (N)
+PbMax = (Pmax * 1000) + ((load/100)*Pmax * 1000)
+//Calculate the initial pre-load Pi (N)
+Pi = PbMax/((kb + kc)/kc)
+//Calculate the maximum force in the bolt PMAX (N)
+PMAX = Pi + ((kb/(kb + kc))*(Pmax * 1000))
+//Calculate the minimum force in the bolt PMIN (N)
+PMIN = Pi + ((kb/(kb + kc))*(Pmin * 1000))
+//Calculate the mean force and force amplitude
+Pm = (PMAX + PMIN)/2
+Pa = (PMAX - PMIN)/2
+//Calculate the actual core cross-section area of the bolt A (mm2)
+A = (Pa + (Pi/((1 + (Sut/Se)) * fs)))/(Sut/((1 + (Sut/Se))*fs))
+//Choose proper diameter from Table 7.1
+//Print results
+printf('\nCore cross-section area of the bolt(A) = %f mm2\n',A)
diff --git a/764/CH7/EX7.19.c/functions7_19.sci b/764/CH7/EX7.19.c/functions7_19.sci
new file mode 100755
index 000000000..b53d32df5
--- /dev/null
+++ b/764/CH7/EX7.19.c/functions7_19.sci
@@ -0,0 +1,48 @@
+
+//Function generating the values of Ka, Kb and Kc
+function [Ka, Kb, Kc] = fluctuate(s, d, r)
+ //Calculate Ka
+ //Nomenclature:
+ //1 - Ground
+ //2 - Machined or cold drawn
+ //3 - Hot-rolled
+ //4 - Forged
+surface = [1 2 3 4]
+Ksurfa = [1.58 4.51 57.7 272]
+Ksurfb = [-0.085 -0.265 -0.718 -0.995]
+for j = 1:1:4
+ if (s == surface(j))
+ a = Ksurfa(j)
+ b = Ksurfb(j)
+ break
+ end
+end
+//From equation 5.18 on page 157
+Ka = a * (Sut^b)
+if (Ka > 1) then
+ Ka = 1
+end
+
+ //Calculate Kb
+ //d (mm)
+ if (d <= 7.5) then
+ Kb = 1
+ elseif ((d > 7.5) & (d <= 50))
+ Kb = 0.85
+ elseif (d > 50)
+ Kb = 0.75
+ else
+ printf('Error in Kb')
+ end
+
+ //Calculate Kc
+ // r (%)
+rel = [50 90 95 99 99.9 99.99 99.999]
+Krel = [1 0.897 0.868 0.814 0.753 0.702 0.659]
+for i = 1:1:7
+ if (r == rel(i)) then
+ Kc = Krel(i)
+ break
+ end
+end
+endfunction
diff --git a/764/CH7/EX7.2.a/data7_2.sci b/764/CH7/EX7.2.a/data7_2.sci
new file mode 100755
index 000000000..47cd1ed40
--- /dev/null
+++ b/764/CH7/EX7.2.a/data7_2.sci
@@ -0,0 +1,9 @@
+
+//(Threaded Joints) Example 7.2
+//Refer Fig.7.14 on page 232
+//Yield tensile strength of 30C8 Syt (N/mm2)
+Syt = 400
+//Factor of safety fs
+fs = 5
+//Shear load acting on the bolts Pt (kN)
+Pt = 5
diff --git a/764/CH7/EX7.2.b/result7_2.txt b/764/CH7/EX7.2.b/result7_2.txt
new file mode 100755
index 000000000..33171a4ac
--- /dev/null
+++ b/764/CH7/EX7.2.b/result7_2.txt
@@ -0,0 +1,28 @@
+-->//(Threaded Joints) Example 7.2
+
+-->//Refer Fig.7.14 on page 232
+
+-->//Yield tensile strength of 30C8 Syt (N/mm2)
+
+-->Syt = 400
+ Syt =
+
+ 400.
+
+-->//Factor of safety fs
+
+-->fs = 5
+ fs =
+
+ 5.
+
+-->//Shear load acting on the bolts Pt (kN)
+
+-->Pt = 5
+ Pt =
+
+ 5.
+
+
+The standard size of the bolt is M10
+ \ No newline at end of file
diff --git a/764/CH7/EX7.2.b/solution7_2.sce b/764/CH7/EX7.2.b/solution7_2.sce
new file mode 100755
index 000000000..7079dc7c8
--- /dev/null
+++ b/764/CH7/EX7.2.b/solution7_2.sce
@@ -0,0 +1,30 @@
+
+//Function to standardise the given bolt-size
+function[v] = standard(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('solution7_2.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_2.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the yield shear strength Ssy (N/mm2)
+Ssy = (50/100)*Syt
+//Calculate the permissible shear stress tau (N/mm2)
+tau = Ssy/fs
+//Shear load acting on one bolt P (kN)
+P = Pt/2
+//Calculate the diameter of the bolt shank d (mm)
+d = ((4 * P * 1000)/(%pi * tau))^(1/2)
+//Standardise the bolt size from Table 7.1
+d = standard(d)
+//Print results
+printf('\nThe standard size of the bolt is M%d\n',d)
diff --git a/764/CH7/EX7.20.a/data7_20.sci b/764/CH7/EX7.20.a/data7_20.sci
new file mode 100755
index 000000000..104c28166
--- /dev/null
+++ b/764/CH7/EX7.20.a/data7_20.sci
@@ -0,0 +1,35 @@
+
+//(Threaded Joints) Example 7.20
+//Number of bolts N
+N = 2
+//Engine speed n (rpm)
+n = 2000
+//Length of stroke l (mm)
+l = 100
+//Length of connecting rod c (mm)
+c = 200
+//Mass of reciprocating parts m (kg)
+m = 5
+//Overload percentage load (%)
+load = 50
+//Assume the stiffness of the bolts to be 1N/mm kb
+kb = 1
+//Calculate the stiffness of the parts kc
+kc = 4 * kb
+//Yield tensile strength of chromium-molybdenum steel Syt (N/mm2)
+Syt = 450
+//Ultimate tensile strength of chromium-molybdenum steel Sut (N/mm2)
+Sut = 600
+//Fatigue stress concentration factor Kf
+Kf = 3.0
+//Expected reliability (%)
+reliability = 90
+//Factor of safety fs
+fs = 2
+//Assume diameter of the bolt to be less than 7.5mm for Kb to be 1
+d = 1
+//As Ka is incorporated into Kf, assume the assembly to be ground
+//This is just to obtain a random value of Ka
+op = 1
+//Inclination of connecting rod to the line of stroke theta (degree)
+theta = 0
diff --git a/764/CH7/EX7.20.b/graph7_20.png b/764/CH7/EX7.20.b/graph7_20.png
new file mode 100755
index 000000000..7ccdd426e
--- /dev/null
+++ b/764/CH7/EX7.20.b/graph7_20.png
Binary files differ
diff --git a/764/CH7/EX7.20.b/result7_20.txt b/764/CH7/EX7.20.b/result7_20.txt
new file mode 100755
index 000000000..e581160f8
--- /dev/null
+++ b/764/CH7/EX7.20.b/result7_20.txt
@@ -0,0 +1,119 @@
+-->//(Threaded Joints) Example 7.20
+
+-->//Number of bolts N
+
+-->N = 2
+ N =
+
+ 2.
+
+-->//Engine speed n (rpm)
+
+-->n = 2000
+ n =
+
+ 2000.
+
+-->//Length of stroke l (mm)
+
+-->l = 100
+ l =
+
+ 100.
+
+-->//Length of connecting rod c (mm)
+
+-->c = 200
+ c =
+
+ 200.
+
+-->//Mass of reciprocating parts m (kg)
+
+-->m = 5
+ m =
+
+ 5.
+
+-->//Overload percentage load (%)
+
+-->load = 50
+ load =
+
+ 50.
+
+-->//Assume the stiffness of the bolts to be 1N/mm kb
+
+-->kb = 1
+ kb =
+
+ 1.
+
+-->//Calculate the stiffness of the parts kc
+
+-->kc = 4 * kb
+ kc =
+
+ 4.
+
+-->//Yield tensile strength of chromium-molybdenum steel Syt (N/mm2)
+
+-->Syt = 450
+ Syt =
+
+ 450.
+
+-->//Ultimate tensile strength of chromium-molybdenum steel Sut (N/mm2)
+
+-->Sut = 600
+ Sut =
+
+ 600.
+
+-->//Fatigue stress concentration factor Kf
+
+-->Kf = 3.0
+ Kf =
+
+ 3.
+
+-->//Expected reliability (%)
+
+-->reliability = 90
+ reliability =
+
+ 90.
+
+-->//Factor of safety fs
+
+-->fs = 2
+ fs =
+
+ 2.
+
+-->//Assume diameter of the bolt to be less than 7.5mm for Kb to be 1
+
+-->d = 1
+ d =
+
+ 1.
+
+-->//As Ka is incorporated into Kf, assume the assembly to be ground
+
+-->//This is just to obtain a random value of Ka
+
+-->op = 1
+ op =
+
+ 1.
+
+-->//Inclination of connecting rod to the line of stroke theta (degree)
+
+-->theta = 0
+ theta =
+
+ 0.
+
+
+Core cross-section area of the bolt(A) = 38.128117 mm2
+ \ No newline at end of file
diff --git a/764/CH7/EX7.20.b/solution7_20.sce b/764/CH7/EX7.20.b/solution7_20.sce
new file mode 100755
index 000000000..68bf7267b
--- /dev/null
+++ b/764/CH7/EX7.20.b/solution7_20.sce
@@ -0,0 +1,72 @@
+
+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('solution7_20.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_20.sci'
+//Obtain path of function file
+funcpath = path + filesep() + 'functions7_20.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+exec(funcpath,[-1])
+//Calculate the endurance limit stress for bolt Sdash (N/mm2)
+Sdash = (50/100)*Sut
+//Calculate Ka, Kb and Kc
+[Ka, Kb, Kc] = fluctuate(op, d, reliability)
+//Calculate Kd
+Kd = 1/Kf
+//Calculate the corrected endurance limit stress Se (N/mm2)
+Se = Kb * Kc * Kd * Sdash
+//Plot modified Goodman diagram
+//The common quadrilateral in the plot is the area of concern
+y1 = {Se 0}
+x1 = {0 Sut}
+y2 = {Syt 0}
+x2 = {0 Syt}
+plot(x1,y1,'--*')
+plot(x2,y2,'-*')
+plot_format()
+title('Modified Goodman diagram (Example 7.20)')
+xlabel('sigmaM (N/mm2)')
+ylabel('sigmaA (N/mm2)')
+//Calculate the angular velocity w (rad/s)
+w = (2 * %pi * n)/60
+//Calculate the crank radius r (m)
+r = (0.5 * l)/1000
+//Calculate the ratio of length of connecting rod to crank radius n1
+n1 = (c/(r * 1000))
+//Calculate the inertia force on bolt at normal running condition I (N)
+I = m * r * (w^2) * (cosd(theta) + (cosd(2 * theta)/n1))
+//Calculate the engine speed at overspeed condition nmax (rpm)
+nmax = n + ((load/100)*n)
+//Calculate the inertia force on bolt at overspeed Imax (N)
+Imax = m * r * (((2 * %pi * nmax)/60)^2) * (cosd(theta) + (cosd(2 * theta)/n1))
+//Calculate the force acting on each bolt under normal running condition P (N)
+P = I/N
+//Calculate the force acting on each bolt under overspeed condition PbMax (N)
+PbMax = Imax/N
+//Calculate the initial pre-load on the bolts Pi (N)
+Pi = PbMax/((kb + kc)/kc)
+//Calculate the maximum force in the bolt PMAX (N)
+PMAX = Pi + ((kb/(kb + kc))* P)
+//Calculate the minimum force in the bolt PMIN (N)
+PMIN = Pi
+//Calculate the mean force and force amplitude
+Pm = (PMAX + PMIN)/2
+Pa = (PMAX - PMIN)/2
+//Calculate the actual core cross-section area of the bolt A (mm2)
+A = (Pa + (Pi/((1 + (Sut/Se)) * fs)))/(Sut/((1 + (Sut/Se))*fs))
+//Choose proper diameter from Table 7.1
+//Print results
+printf('\nCore cross-section area of the bolt(A) = %f mm2\n',A)
diff --git a/764/CH7/EX7.20.c/functions7_20.sci b/764/CH7/EX7.20.c/functions7_20.sci
new file mode 100755
index 000000000..b53d32df5
--- /dev/null
+++ b/764/CH7/EX7.20.c/functions7_20.sci
@@ -0,0 +1,48 @@
+
+//Function generating the values of Ka, Kb and Kc
+function [Ka, Kb, Kc] = fluctuate(s, d, r)
+ //Calculate Ka
+ //Nomenclature:
+ //1 - Ground
+ //2 - Machined or cold drawn
+ //3 - Hot-rolled
+ //4 - Forged
+surface = [1 2 3 4]
+Ksurfa = [1.58 4.51 57.7 272]
+Ksurfb = [-0.085 -0.265 -0.718 -0.995]
+for j = 1:1:4
+ if (s == surface(j))
+ a = Ksurfa(j)
+ b = Ksurfb(j)
+ break
+ end
+end
+//From equation 5.18 on page 157
+Ka = a * (Sut^b)
+if (Ka > 1) then
+ Ka = 1
+end
+
+ //Calculate Kb
+ //d (mm)
+ if (d <= 7.5) then
+ Kb = 1
+ elseif ((d > 7.5) & (d <= 50))
+ Kb = 0.85
+ elseif (d > 50)
+ Kb = 0.75
+ else
+ printf('Error in Kb')
+ end
+
+ //Calculate Kc
+ // r (%)
+rel = [50 90 95 99 99.9 99.99 99.999]
+Krel = [1 0.897 0.868 0.814 0.753 0.702 0.659]
+for i = 1:1:7
+ if (r == rel(i)) then
+ Kc = Krel(i)
+ break
+ end
+end
+endfunction
diff --git a/764/CH7/EX7.21.a/data7_21.sci b/764/CH7/EX7.21.a/data7_21.sci
new file mode 100755
index 000000000..141f55e6a
--- /dev/null
+++ b/764/CH7/EX7.21.a/data7_21.sci
@@ -0,0 +1,31 @@
+
+//(Threaded Joints) Example 7.21
+//Maximum internal pressure MAXP (MPa)
+MAXP = 2
+//Minimum internal pressure MINP (MPa)
+MINP = 0
+//Diameter of the circle on which pressure acts D (mm)
+D = 400
+//Assume the stiffness of the bolts to be 1N/mm kb
+kb = 1
+//Calculate the stiffness of the parts kc (N/mm)
+kc = 4 * kb
+//Ultimate tensile strength Sut (N/mm2)
+Sut = 900
+//Yield strength Syt (N/mm2)
+Syt = 700
+//Endurance limit in bending Sdash (N/mm2)
+Sdash = 300
+//Fatigue stress concentration factor Kf
+Kf = 2.2
+//Factor of safety fs
+fs = 1.5
+//Number of bolts N
+N = 8
+//Assume expected reliability (%)
+reliability = 90
+//Assume diameter of the bolt to be less than 7.5mm for Kb to be 1
+d = 1
+//As Ka is incorporated into Kf, assume the assembly to be ground
+//This is just to obtain a random value of Ka
+op = 1
diff --git a/764/CH7/EX7.21.b/graph7_21.png b/764/CH7/EX7.21.b/graph7_21.png
new file mode 100755
index 000000000..aca198568
--- /dev/null
+++ b/764/CH7/EX7.21.b/graph7_21.png
Binary files differ
diff --git a/764/CH7/EX7.21.b/result7_21.txt b/764/CH7/EX7.21.b/result7_21.txt
new file mode 100755
index 000000000..640126a35
--- /dev/null
+++ b/764/CH7/EX7.21.b/result7_21.txt
@@ -0,0 +1,105 @@
+-->//(Threaded Joints) Example 7.21
+
+-->//Maximum internal pressure MAXP (MPa)
+
+-->MAXP = 2
+ MAXP =
+
+ 2.
+
+-->//Minimum internal pressure MINP (MPa)
+
+-->MINP = 0
+ MINP =
+
+ 0.
+
+-->//Diameter of the circle on which pressure acts D (mm)
+
+-->D = 400
+ D =
+
+ 400.
+
+-->//Assume the stiffness of the bolts to be 1N/mm kb
+
+-->kb = 1
+ kb =
+
+ 1.
+
+-->//Calculate the stiffness of the parts kc (N/mm)
+
+-->kc = 4 * kb
+ kc =
+
+ 4.
+
+-->//Ultimate tensile strength Sut (N/mm2)
+
+-->Sut = 900
+ Sut =
+
+ 900.
+
+-->//Yield strength Syt (N/mm2)
+
+-->Syt = 700
+ Syt =
+
+ 700.
+
+-->//Endurance limit in bending Sdash (N/mm2)
+
+-->Sdash = 300
+ Sdash =
+
+ 300.
+
+-->//Fatigue stress concentration factor Kf
+
+-->Kf = 2.2
+ Kf =
+
+ 2.2
+
+-->//Factor of safety fs
+
+-->fs = 1.5
+ fs =
+
+ 1.5
+
+-->//Number of bolts N
+
+-->N = 8
+ N =
+
+ 8.
+
+-->//Assume expected reliability (%)
+
+-->reliability = 90
+ reliability =
+
+ 90.
+
+-->//Assume diameter of the bolt to be less than 7.5mm for Kb to be 1
+
+-->d = 1
+ d =
+
+ 1.
+
+-->//As Ka is incorporated into Kf, assume the assembly to be ground
+
+-->//This is just to obtain a random value of Ka
+
+-->op = 1
+ op =
+
+ 1.
+
+
+Core cross-section area of the bolt(A) = 85.172067 mm2
+ \ No newline at end of file
diff --git a/764/CH7/EX7.21.b/solution7_21.sce b/764/CH7/EX7.21.b/solution7_21.sce
new file mode 100755
index 000000000..acef2b3bf
--- /dev/null
+++ b/764/CH7/EX7.21.b/solution7_21.sce
@@ -0,0 +1,58 @@
+
+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('solution7_21.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_21.sci'
+//Obtain path of function file
+funcpath = path + filesep() + 'functions7_21.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+exec(funcpath,[-1])
+//Calculate Ka, Kb and Kc
+[Ka, Kb, Kc] = fluctuate(op, d, reliability)
+//Calculate Kd
+Kd = 1/Kf
+//Calculate the corrected endurance limit stress Se (N/mm2)
+Se = Kb * Kd * Sdash
+//Plot modified Goodman diagram
+//The common quadrilateral in the plot is the area of concern
+y1 = {Se 0}
+x1 = {0 Sut}
+y2 = {Syt 0}
+x2 = {0 Syt}
+plot(x1,y1,'--*')
+plot(x2,y2,'-*')
+plot_format()
+title('Modified Goodman diagram (Example 7.21)')
+xlabel('sigmaM (N/mm2)')
+ylabel('sigmaA (N/mm2)')
+//Calculate the maximum and minimum force of fluid on cover P (N)
+P1 = (%pi/4)*(D^2)*MAXP
+P2 = (%pi/4)*(D^2)*MINP
+//Calculate the initial pre-load on the bolts Pi (N)
+Pi = 1.3 * P1
+//Calculate the maximum and minimum forces in the bolt (N)
+Pmax = Pi + ((kb/(kb + kc))*P1)
+Pmin = Pi + ((kb/(kb + kc))*P2)
+//Calculate the mean force and force amplitude
+Pm = (Pmax + Pmin)/2
+Pa = (Pmax - Pmin)/2
+//Calculate the actual core cross-section area of the bolts Atotal (mm2)
+Atotal = (Pa + (Pi/((1 + (Sut/Se)) * fs)))/(Sut/((1 + (Sut/Se))*fs))
+//Calculate the cross-section area of each bolt A (mm2)
+A = Atotal/N
+//Choose proper diameter from Table 7.1
+//Print results
+printf('\nCore cross-section area of the bolt(A) = %f mm2\n',A)
diff --git a/764/CH7/EX7.21.c/functions7_21.sci b/764/CH7/EX7.21.c/functions7_21.sci
new file mode 100755
index 000000000..b53d32df5
--- /dev/null
+++ b/764/CH7/EX7.21.c/functions7_21.sci
@@ -0,0 +1,48 @@
+
+//Function generating the values of Ka, Kb and Kc
+function [Ka, Kb, Kc] = fluctuate(s, d, r)
+ //Calculate Ka
+ //Nomenclature:
+ //1 - Ground
+ //2 - Machined or cold drawn
+ //3 - Hot-rolled
+ //4 - Forged
+surface = [1 2 3 4]
+Ksurfa = [1.58 4.51 57.7 272]
+Ksurfb = [-0.085 -0.265 -0.718 -0.995]
+for j = 1:1:4
+ if (s == surface(j))
+ a = Ksurfa(j)
+ b = Ksurfb(j)
+ break
+ end
+end
+//From equation 5.18 on page 157
+Ka = a * (Sut^b)
+if (Ka > 1) then
+ Ka = 1
+end
+
+ //Calculate Kb
+ //d (mm)
+ if (d <= 7.5) then
+ Kb = 1
+ elseif ((d > 7.5) & (d <= 50))
+ Kb = 0.85
+ elseif (d > 50)
+ Kb = 0.75
+ else
+ printf('Error in Kb')
+ end
+
+ //Calculate Kc
+ // r (%)
+rel = [50 90 95 99 99.9 99.99 99.999]
+Krel = [1 0.897 0.868 0.814 0.753 0.702 0.659]
+for i = 1:1:7
+ if (r == rel(i)) then
+ Kc = Krel(i)
+ break
+ end
+end
+endfunction
diff --git a/764/CH7/EX7.22.a/data7_22.sci b/764/CH7/EX7.22.a/data7_22.sci
new file mode 100755
index 000000000..925874965
--- /dev/null
+++ b/764/CH7/EX7.22.a/data7_22.sci
@@ -0,0 +1,35 @@
+
+//(Threaded Joints) Example 7.22
+//Refer Fig.7.43 and 7.44 on page 265
+//Weight of machine W (kg)
+W = 200
+//Distance between the machine and the mearest point of support dist (m)
+dist = 1
+//Rotating unbalanced force created by the machine F (N)
+F = 2000
+//Speed of rotation n1 (rpm)
+n1 = 14
+//Weight of channel w (kg/m)
+w = 20
+//Number of bolts N
+N = 2
+//Distance between bolt1 and the nearest point of support l1 (mm)
+l1 = 235 + 35
+//Distance between bolt2 and the nearest point of support l2 (mm)
+l2 = 35
+//Ultimate tensile strength Sut (MPa)
+Sut = 960
+//Yield tensile strength Syt (MPa)
+Syt = 850
+//Endurance limit in bending Sdash (MPa)
+Sdash = 500
+//Fatigue stress concentration factor Kf
+Kf = 3.0
+//Factor of safety fs
+fs = 2
+//Initial pre-load in each bolt Pi (kN)
+Pi = 55
+//Assume the stiffness of the bolts to be 1N/mm kb
+kb = 1
+//Calculate the stiffness of the parts kc
+kc = 3 * kb
diff --git a/764/CH7/EX7.22.b/graph7_22.png b/764/CH7/EX7.22.b/graph7_22.png
new file mode 100755
index 000000000..1642a09ee
--- /dev/null
+++ b/764/CH7/EX7.22.b/graph7_22.png
Binary files differ
diff --git a/764/CH7/EX7.22.b/result7_22.txt b/764/CH7/EX7.22.b/result7_22.txt
new file mode 100755
index 000000000..a95cfbfb2
--- /dev/null
+++ b/764/CH7/EX7.22.b/result7_22.txt
@@ -0,0 +1,121 @@
+-->//(Threaded Joints) Example 7.22
+
+-->//Refer Fig.7.43 and 7.44 on page 265
+
+-->//Weight of machine W (kg)
+
+-->W = 200
+ W =
+
+ 200.
+
+-->//Distance between the machine and the mearest point of support dist (m)
+
+-->dist = 1
+ dist =
+
+ 1.
+
+-->//Rotating unbalanced force created by the machine F (N)
+
+-->F = 2000
+ F =
+
+ 2000.
+
+-->//Speed of rotation n1 (rpm)
+
+-->n1 = 14
+ n1 =
+
+ 14.
+
+-->//Weight of channel w (kg/m)
+
+-->w = 20
+ w =
+
+ 20.
+
+-->//Number of bolts N
+
+-->N = 2
+ N =
+
+ 2.
+
+-->//Distance between bolt1 and the nearest point of support l1 (mm)
+
+-->l1 = 235 + 35
+ l1 =
+
+ 270.
+
+-->//Distance between bolt2 and the nearest point of support l2 (mm)
+
+-->l2 = 35
+ l2 =
+
+ 35.
+
+-->//Ultimate tensile strength Sut (MPa)
+
+-->Sut = 960
+ Sut =
+
+ 960.
+
+-->//Yield tensile strength Syt (MPa)
+
+-->Syt = 850
+ Syt =
+
+ 850.
+
+-->//Endurance limit in bending Sdash (MPa)
+
+-->Sdash = 500
+ Sdash =
+
+ 500.
+
+-->//Fatigue stress concentration factor Kf
+
+-->Kf = 3.0
+ Kf =
+
+ 3.
+
+-->//Factor of safety fs
+
+-->fs = 2
+ fs =
+
+ 2.
+
+-->//Initial pre-load in each bolt Pi (kN)
+
+-->Pi = 55
+ Pi =
+
+ 55.
+
+-->//Assume the stiffness of the bolts to be 1N/mm kb
+
+-->kb = 1
+ kb =
+
+ 1.
+
+-->//Calculate the stiffness of the parts kc
+
+-->kc = 3 * kb
+ kc =
+
+ 3.
+
+
+Core cross-section area of the bolt(A) = 197.638259 mm2
+
+Design is safe
+ \ No newline at end of file
diff --git a/764/CH7/EX7.22.b/solution7_22.sce b/764/CH7/EX7.22.b/solution7_22.sce
new file mode 100755
index 000000000..c66132b5e
--- /dev/null
+++ b/764/CH7/EX7.22.b/solution7_22.sce
@@ -0,0 +1,68 @@
+
+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('solution7_22.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_22.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the value of Kd
+Kd = 1/Kf
+//Calculate the endurance limit stress for bolt Se (N/mm2)
+Se = Kd * Sdash
+//Case1: Rotating force acting downward
+//Calculate the force P1d (N)
+P1d = ((w * 9.81 * (dist/2) * 1000) + (((W * 9.81) + F)*(dist * 1000)))/((l1 + ((l2^2)/l1)))
+//Case2: Rotating force acting upward
+//Calculate the force P1u (N)
+P1u = ((w * 9.81 * (dist/2) * 1000) + (((W * 9.81) - F)*(dist * 1000)))/((l1 + ((l2^2)/l1)))
+//Plot modified Goodman diagram
+//The common quadrilateral in the plot is the area of concern
+y1 = {Se 0}
+x1 = {0 Sut}
+y2 = {Syt 0}
+x2 = {0 Syt}
+plot(x1,y1,'--*')
+plot(x2,y2,'-*')
+plot_format()
+title('Modified Goodman diagram (Example 7.22)')
+xlabel('sigmaM (N/mm2)')
+ylabel('sigmaA (N/mm2)')
+//Calculate the maximum and minimum forces on the bolt (N)
+if (P1d > P1u) then
+ Pmax = (Pi * 1000) + ((kb/(kb + kc))*P1d)
+ Pmin = (Pi * 1000) + ((kb/(kb + kc))*P1u)
+else
+ Pmin = (Pi * 1000) + ((kb/(kb + kc))*P1d)
+ Pmax = (Pi * 1000) + ((kb/(kb + kc))*P1u)
+end
+//Calculate the mean force and force amplitude (N)
+Pm = (Pmax + Pmin)/2
+Pa = (Pmax - Pmin)/2
+theta = atand(Pa/Pm)
+//Calculate the actual core cross-section area of the bolt A (mm2)
+A = (Pa + ((Pi*1000)/(((1/tand(theta)) + (Sut/Se)) * fs)))/(Sut/(((1/tand(theta)) + (Sut/Se))*fs))
+//Choose proper diameter from Table 7.1
+//Check for static design
+//Calculate maximum tensile stress sigmat (N/mm2)
+sigmat = Pmax/A
+//Calculate the factor of safety fsNew
+fsNew = Syt/sigmat
+//Print results
+printf('\nCore cross-section area of the bolt(A) = %f mm2\n',A)
+if (fsNew > fs) then
+ printf('\nDesign is safe\n')
+else
+ printf('\nDesign is not safe\n')
+end
diff --git a/764/CH7/EX7.23.a/data7_23.sci b/764/CH7/EX7.23.a/data7_23.sci
new file mode 100755
index 000000000..7ca327ea5
--- /dev/null
+++ b/764/CH7/EX7.23.a/data7_23.sci
@@ -0,0 +1,29 @@
+
+//(Threaded Joints) Example 7.23
+//Refer Fig.7.47 on page 267
+//Maximum pressure in the vessel Pmax (MPa)
+Pmax = 1
+//Minimum pressure in the vessel Pmin (MPa)
+Pmin = 0
+//Seating pressure for the gasket PSeat (MPa)
+PSeat = 5
+//Number of bolts N
+N = 8
+//Assume the stiffness of the bolts to be 1N/mm kb
+kb = 1
+//Calculate the stiffness of the parts kc
+kc = 4 * kb
+//Factor of safety fs
+fs = 2
+//Ultimate tensile strength of bolt material Sut (N/mm2)
+Sut = 780
+//Yield tensile strength of the bolt material Syt (N/mm2)
+Syt = 580
+//Endurance limit in bending Sdash (N/mm2)
+Sdash = 260
+//Fatigue stress concentration factor Kf
+Kf = 3
+//Inner diameter of the gasket Di (mm)
+Di = 300
+//Outer diameter of the gasket Do (mm)
+Do = 300 + (2 * 50)
diff --git a/764/CH7/EX7.23.b/graph7_23.png b/764/CH7/EX7.23.b/graph7_23.png
new file mode 100755
index 000000000..822633721
--- /dev/null
+++ b/764/CH7/EX7.23.b/graph7_23.png
Binary files differ
diff --git a/764/CH7/EX7.23.b/result7_23.txt b/764/CH7/EX7.23.b/result7_23.txt
new file mode 100755
index 000000000..706e583c2
--- /dev/null
+++ b/764/CH7/EX7.23.b/result7_23.txt
@@ -0,0 +1,98 @@
+-->//(Threaded Joints) Example 7.23
+
+-->//Refer Fig.7.47 on page 267
+
+-->//Maximum pressure in the vessel Pmax (MPa)
+
+-->Pmax = 1
+ Pmax =
+
+ 1.
+
+-->//Minimum pressure in the vessel Pmin (MPa)
+
+-->Pmin = 0
+ Pmin =
+
+ 0.
+
+-->//Seating pressure for the gasket PSeat (MPa)
+
+-->PSeat = 5
+ PSeat =
+
+ 5.
+
+-->//Number of bolts N
+
+-->N = 8
+ N =
+
+ 8.
+
+-->//Assume the stiffness of the bolts to be 1N/mm kb
+
+-->kb = 1
+ kb =
+
+ 1.
+
+-->//Calculate the stiffness of the parts kc
+
+-->kc = 4 * kb
+ kc =
+
+ 4.
+
+-->//Factor of safety fs
+
+-->fs = 2
+ fs =
+
+ 2.
+
+-->//Ultimate tensile strength of bolt material Sut (N/mm2)
+
+-->Sut = 780
+ Sut =
+
+ 780.
+
+-->//Yield tensile strength of the bolt material Syt (N/mm2)
+
+-->Syt = 580
+ Syt =
+
+ 580.
+
+-->//Endurance limit in bending Sdash (N/mm2)
+
+-->Sdash = 260
+ Sdash =
+
+ 260.
+
+-->//Fatigue stress concentration factor Kf
+
+-->Kf = 3
+ Kf =
+
+ 3.
+
+-->//Inner diameter of the gasket Di (mm)
+
+-->Di = 300
+ Di =
+
+ 300.
+
+-->//Outer diameter of the gasket Do (mm)
+
+-->Do = 300 + (2 * 50)
+ Do =
+
+ 400.
+
+
+Core cross-section area of the bolt(A) = 74.889729 mm2
+ \ No newline at end of file
diff --git a/764/CH7/EX7.23.b/solution7_23.sce b/764/CH7/EX7.23.b/solution7_23.sce
new file mode 100755
index 000000000..11467f4cd
--- /dev/null
+++ b/764/CH7/EX7.23.b/solution7_23.sce
@@ -0,0 +1,54 @@
+
+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('solution7_23.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_23.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the value of Kd
+Kd = 1/Kf
+//Calculate the endurance limit stress for bolt Se (N/mm2)
+Se = Kd * Sdash
+//Plot modified Goodman diagram
+//The common quadrilateral in the plot is the area of concern
+y1 = {Se 0}
+x1 = {0 Sut}
+y2 = {Syt 0}
+x2 = {0 Syt}
+plot(x1,y1,'--*')
+plot(x2,y2,'-*')
+plot_format()
+title('Modified Goodman diagram (Example 7.23)')
+xlabel('sigmaM (N/mm2)')
+ylabel('sigmaA (N/mm2)')
+//Calculate the initial pre-load in the bolts Pitotal (N)
+Pitotal = PSeat*((%pi/4)*((Do^2) - (Di^2)))
+//Calculate the pre-load per bolt Pi (N)
+Pi = Pitotal/N
+//Calculate the total external load Fmax per bolt(N)
+Fmax = ((%pi/4)*((Di + ((Do - Di)/2))^2)*Pmax)/N
+//Calculate the total external load Fmin per bolt(N)
+Fmin = ((%pi/4)*((Di + ((Do - Di)/2))^2)*Pmin)/N
+//Calculate the maximum and minimum forces on the bolt (N)
+PMAX = Pi + ((kb/(kb + kc)) * Fmax)
+PMIN = Pi + ((kb/(kb + kc)) * Fmin)
+//Calculate the mean force and force amplitude (N)
+Pm = (PMAX + PMIN)/2
+Pa = (PMAX - PMIN)/2
+//Calculate the actual core cross-section area of the bolt A (mm2)
+A = (Pa + (Pi/((1 + (Sut/Se)) * fs)))/(Sut/((1 + (Sut/Se))*fs))
+//Choose proper diameter from Table 7.1
+//Print results
+printf('\nCore cross-section area of the bolt(A) = %f mm2\n',A)
diff --git a/764/CH7/EX7.3.a/data7_3.sci b/764/CH7/EX7.3.a/data7_3.sci
new file mode 100755
index 000000000..b04ddc374
--- /dev/null
+++ b/764/CH7/EX7.3.a/data7_3.sci
@@ -0,0 +1,17 @@
+
+//(Threaded Joints) Example 7.3
+//Refer Fig.7.16 on page 233
+//Eccentric force acting on the structure P (kN)
+P = 10
+//Eccentricity of the force from the C.G. of the bolts e (mm)
+e = 500
+//Centre distance between bolts 1 and 2 dist1 (mm)
+dist1 = 200
+//Centre distance between bolts 1 and 3 dist2 (mm)
+dist2 = 150
+//Tensile yield strength of 30C8 Syt (N/mm2)
+Syt = 400
+//Factor of safety fs
+fs = 2.5
+//Number of bolts N
+N = 4
diff --git a/764/CH7/EX7.3.b/result7_3.txt b/764/CH7/EX7.3.b/result7_3.txt
new file mode 100755
index 000000000..1c7959344
--- /dev/null
+++ b/764/CH7/EX7.3.b/result7_3.txt
@@ -0,0 +1,56 @@
+-->//(Threaded Joints) Example 7.3
+
+-->//Refer Fig.7.16 on page 233
+
+-->//Eccentric force acting on the structure P (kN)
+
+-->P = 10
+ P =
+
+ 10.
+
+-->//Eccentricity of the force from the C.G. of the bolts e (mm)
+
+-->e = 500
+ e =
+
+ 500.
+
+-->//Centre distance between bolts 1 and 2 dist1 (mm)
+
+-->dist1 = 200
+ dist1 =
+
+ 200.
+
+-->//Centre distance between bolts 1 and 3 dist2 (mm)
+
+-->dist2 = 150
+ dist2 =
+
+ 150.
+
+-->//Tensile yield strength of 30C8 Syt (N/mm2)
+
+-->Syt = 400
+ Syt =
+
+ 400.
+
+-->//Factor of safety fs
+
+-->fs = 2.5
+ fs =
+
+ 2.5
+
+-->//Number of bolts N
+
+-->N = 4
+ N =
+
+ 4.
+
+
+Standard size of the bolts is M20
+ \ No newline at end of file
diff --git a/764/CH7/EX7.3.b/solution7_3.sce b/764/CH7/EX7.3.b/solution7_3.sce
new file mode 100755
index 000000000..ac216e404
--- /dev/null
+++ b/764/CH7/EX7.3.b/solution7_3.sce
@@ -0,0 +1,47 @@
+
+//Function to standardise the given bolt-size
+function[v] = standard(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('solution7_3.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_3.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the permissible shear stress tau (N/mm2)
+tau = ((50/100)*Syt)/fs
+//Calculate the distance between bolt1 and C.G. of all bolts res (mm)
+res = (((dist1/2)^2) + ((dist2/2)^2))^(1/2)
+//Calculate the primary shear force at bolt1 Pshear (N)
+Pshear = (P * 1000)/N
+//Calculate the secondary shear force at bolt1 Sshear (N)
+Sshear = (P * 1000 * e)/(N * res)
+//Calculate angle theta (degree)
+theta = atand(dist2/dist1)
+//Calculate the resultant force on bolt1 P1 (N)
+P1 = (((Sshear * cosd(theta) - Pshear)^2) + ((Sshear * sind(theta))^2))^(1/2)
+//Calculate the resultant force on bolt2 P2 (N)
+P2 = (((Sshear * cosd(theta) + Pshear)^2) + ((Sshear * sind(theta))^2))^(1/2)
+//Obtain the bolt subjected to maximum shear force Pmax (N)
+if (P1 > P2) then
+ Pmax = P1
+else
+ Pmax = P2
+end
+//Calculate the core diameter of the bolt dc (mm)
+dc = ((4 * Pmax)/(tau * %pi))^(1/2)
+//Calculate the nominal diameter of the bolt d (mm)
+d = dc/0.8
+//Standardise the bolt size
+d = standard(d)
+//Print results
+printf('\nStandard size of the bolts is M%d\n',d)
diff --git a/764/CH7/EX7.4.a/data7_4.sci b/764/CH7/EX7.4.a/data7_4.sci
new file mode 100755
index 000000000..8321fdb20
--- /dev/null
+++ b/764/CH7/EX7.4.a/data7_4.sci
@@ -0,0 +1,17 @@
+
+//(Threaded Joints) Example 7.4
+//Refer Fig.7.19 on page 235
+//Force acting on the steel plate P (kN)
+P = 5
+//Yield tensile strength of 45C8 Syt (N/mm2)
+Syt = 380
+//Factor of safety fs
+fs = 3
+//Number of bolts N
+N = 3
+//Eccentricity value (distance between C.G. of bolts and force) e (mm)
+e = 200 + 30 + 75
+//Distance between bolt1 and bolt2 r1 (mm)
+r1 = 75
+//Distance between bolt2 and bolt3 r3 (mm)
+r3 = 75
diff --git a/764/CH7/EX7.4.b/result7_4.txt b/764/CH7/EX7.4.b/result7_4.txt
new file mode 100755
index 000000000..24a9dd2a1
--- /dev/null
+++ b/764/CH7/EX7.4.b/result7_4.txt
@@ -0,0 +1,56 @@
+-->//(Threaded Joints) Example 7.4
+
+-->//Refer Fig.7.19 on page 235
+
+-->//Force acting on the steel plate P (kN)
+
+-->P = 5
+ P =
+
+ 5.
+
+-->//Yield tensile strength of 45C8 Syt (N/mm2)
+
+-->Syt = 380
+ Syt =
+
+ 380.
+
+-->//Factor of safety fs
+
+-->fs = 3
+ fs =
+
+ 3.
+
+-->//Number of bolts N
+
+-->N = 3
+ N =
+
+ 3.
+
+-->//Eccentricity value (distance between C.G. of bolts and force) e (mm)
+
+-->e = 200 + 30 + 75
+ e =
+
+ 305.
+
+-->//Distance between bolt1 and bolt2 r1 (mm)
+
+-->r1 = 75
+ r1 =
+
+ 75.
+
+-->//Distance between bolt2 and bolt3 r3 (mm)
+
+-->r3 = 75
+ r3 =
+
+ 75.
+
+
+The standard size of the bolts is M20
+ \ No newline at end of file
diff --git a/764/CH7/EX7.4.b/solution7_4.sce b/764/CH7/EX7.4.b/solution7_4.sce
new file mode 100755
index 000000000..d9a2ee3c3
--- /dev/null
+++ b/764/CH7/EX7.4.b/solution7_4.sce
@@ -0,0 +1,34 @@
+
+//Function to standardise the given bolt-size
+function[v] = standard(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('solution7_4.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_4.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the permissible shear stress tau (N/mm2)
+tau = ((50/100)*Syt)/fs
+//Calculate the primary shear force on bolt3 Pshear (N)
+Pshear = (P * 1000)/N
+//Calculate the secondary shear force on bolt3 Sshear (N)
+Sshear = (P * 1000 * e * r1)/((r1^2) + (r3^2))
+//Calculate the resultant force on bolt3 P3 (N)
+P3 = Pshear + Sshear
+//Calculate the core diameter of the bolt dc (mm)
+dc = ((4 * P3)/(%pi * tau))^(1/2)
+//Calculate the nominal diameter of the bolt d (mm)
+d = dc/0.8
+//Standardise the bolt size
+d = standard(d)
+//Print results
+printf('\nThe standard size of the bolts is M%d\n',d)
diff --git a/764/CH7/EX7.5.a/data7_5.sci b/764/CH7/EX7.5.a/data7_5.sci
new file mode 100755
index 000000000..5855205a7
--- /dev/null
+++ b/764/CH7/EX7.5.a/data7_5.sci
@@ -0,0 +1,17 @@
+
+//(Threaded Joints) Example 7.5
+//Refer Fig.7.20 on page 236
+//Load applied on the bracket P (kN)
+P = 25
+//Eccentricity value e (mm)
+e = 100
+//Height of bolt1 from the base l1 (mm)
+l1 = 150
+//Height of bolt2 from the base l2 (mm)
+l2 = 25
+//Yield tensile strength of 45C8 Syt (N/mm2)
+Syt = 380
+//Factor of safety fs
+fs = 2.5
+//Number of bolts N
+N = 4
diff --git a/764/CH7/EX7.5.b/result7_5.txt b/764/CH7/EX7.5.b/result7_5.txt
new file mode 100755
index 000000000..c6659c4b2
--- /dev/null
+++ b/764/CH7/EX7.5.b/result7_5.txt
@@ -0,0 +1,56 @@
+-->//(Threaded Joints) Example 7.5
+
+-->//Refer Fig.7.20 on page 236
+
+-->//Load applied on the bracket P (kN)
+
+-->P = 25
+ P =
+
+ 25.
+
+-->//Eccentricity value e (mm)
+
+-->e = 100
+ e =
+
+ 100.
+
+-->//Height of bolt1 from the base l1 (mm)
+
+-->l1 = 150
+ l1 =
+
+ 150.
+
+-->//Height of bolt2 from the base l2 (mm)
+
+-->l2 = 25
+ l2 =
+
+ 25.
+
+-->//Yield tensile strength of 45C8 Syt (N/mm2)
+
+-->Syt = 380
+ Syt =
+
+ 380.
+
+-->//Factor of safety fs
+
+-->fs = 2.5
+ fs =
+
+ 2.5
+
+-->//Number of bolts N
+
+-->N = 4
+ N =
+
+ 4.
+
+
+Area at the core cross-section(A) = 98.022213 mm2
+ \ No newline at end of file
diff --git a/764/CH7/EX7.5.b/solution7_5.sce b/764/CH7/EX7.5.b/solution7_5.sce
new file mode 100755
index 000000000..55ae0923c
--- /dev/null
+++ b/764/CH7/EX7.5.b/solution7_5.sce
@@ -0,0 +1,33 @@
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution7_5.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_5.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the permissible shear stress tau (N/mm2)
+tau = ((50/100)*Syt)/fs
+//Calculate the direct shear force in bolt Dshear (N)
+Dshear = (P * 1000)/N
+//Assume the area at core cross-section to be 1mm2 A
+A = 1
+//Calculate the direct shear stress in bolt Sshear (N/mm2)
+Sshear = Dshear/A
+//Calculate the tensile force on the bolts Ftensile (N)
+if (l1 > l2) then
+ Ftensile = (P * 1000 * e * l1)/(2*((l1^2) + (l2^2)))
+else
+ Ftensile = (P * 1000 * e * l2)/(2*((l1^2) + (l2^2)))
+end
+//Calculate tensile stress in appropriate bolt sigmat (N/mm2)
+sigmat = Ftensile/A
+//Calculate the principal shear stress in bolt tauMax (N/mm2)
+tauMax = (((sigmat/2)^2) + (Sshear^2))^(1/2)
+//Calculate the actual area of core of cross-section A (mm2)
+A = tauMax/tau
+//Choose proper diameter from Table 7.1
+//Print results
+printf('\nArea at the core cross-section(A) = %f mm2\n',A)
+
diff --git a/764/CH7/EX7.6.a/data7_6.sci b/764/CH7/EX7.6.a/data7_6.sci
new file mode 100755
index 000000000..c8264780c
--- /dev/null
+++ b/764/CH7/EX7.6.a/data7_6.sci
@@ -0,0 +1,15 @@
+
+//(Threaded Joints) Example 7.6
+//Refer Fig.7.21 on page 237
+//Number of bolts N
+N = 4
+//Allowable tensile stress in the bolts sigmaMax (N/mm2)
+sigmaMax = 35
+//Load acting on the bracket P (kN)
+P = 25
+//Eccentricity value e (mm)
+e = 500
+//Height of bolt1 from the base l1 (mm)
+l1 = 550
+//Height of bolt2 from the base l2 (mm)
+l2 = 50
diff --git a/764/CH7/EX7.6.b/result7_6.txt b/764/CH7/EX7.6.b/result7_6.txt
new file mode 100755
index 000000000..492bb7078
--- /dev/null
+++ b/764/CH7/EX7.6.b/result7_6.txt
@@ -0,0 +1,49 @@
+-->//(Threaded Joints) Example 7.6
+
+-->//Refer Fig.7.21 on page 237
+
+-->//Number of bolts N
+
+-->N = 4
+ N =
+
+ 4.
+
+-->//Allowable tensile stress in the bolts sigmaMax (N/mm2)
+
+-->sigmaMax = 35
+ sigmaMax =
+
+ 35.
+
+-->//Load acting on the bracket P (kN)
+
+-->P = 25
+ P =
+
+ 25.
+
+-->//Eccentricity value e (mm)
+
+-->e = 500
+ e =
+
+ 500.
+
+-->//Height of bolt1 from the base l1 (mm)
+
+-->l1 = 550
+ l1 =
+
+ 550.
+
+-->//Height of bolt2 from the base l2 (mm)
+
+-->l2 = 50
+ l2 =
+
+ 50.
+
+
+Area at the core cross-section(A) = 401.446244 mm2
+ \ No newline at end of file
diff --git a/764/CH7/EX7.6.b/solution7_6.sce b/764/CH7/EX7.6.b/solution7_6.sce
new file mode 100755
index 000000000..f502d4716
--- /dev/null
+++ b/764/CH7/EX7.6.b/solution7_6.sce
@@ -0,0 +1,30 @@
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution7_6.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_6.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the direct shear force on each bolt Dshear (N)
+Dshear = (P * 1000)/N
+//Assume the core cross-section area to be 1mm2 A
+A = 1
+//Calculate the direct shear stress in each bolt tau (N/mm2)
+tau = Dshear/A
+//Calculate the tensile force on appropriate bolt Ftensile (N)
+if (l1 > l2) then
+ Ftensile = (P * 1000 * e *l1)/(2*((l1^2) + (l2^2)))
+else
+ Ftensile = (P * 1000 * e *l2)/(2*((l1^2) + (l2^2)))
+end
+//Calculate the tensile stress in bolt1 sigmat (N/mm2)
+sigmat = Ftensile/A
+//Calculate the principal stress in bolt1 S (N/mm2)
+S = (sigmat/2) + (((sigmat/2)^2) + (tau^2))^(1/2)
+//Calculate the actual core cross-section area of the bolt A (mm2)
+A = S/sigmaMax
+//Choose proper diameter from Table 7.1
+//Print results
+printf('\nArea at the core cross-section(A) = %f mm2\n',A)
diff --git a/764/CH7/EX7.7.a/data7_7.sci b/764/CH7/EX7.7.a/data7_7.sci
new file mode 100755
index 000000000..3cf1b0910
--- /dev/null
+++ b/764/CH7/EX7.7.a/data7_7.sci
@@ -0,0 +1,17 @@
+
+//(Threaded Joints) Example 7.7
+//Refer Fig.7.22 on page 238
+//Number of bolts N
+N = 6
+//Load acting on the bracket P (kN)
+P = 50
+//Height of bolt1 from the base l1 (mm)
+l1 = 300
+//Height of bolt2 from the base l2 (mm)
+l2 = 200
+//Height of bolt3 from the base l3 (mm)
+l3 = 100
+//Eccentricity value l (mm)
+l = 250
+//Maximum permissible tensile stress in bolts sigmaMax (N/mm2)
+sigmaMax = 100
diff --git a/764/CH7/EX7.7.b/result7_7.txt b/764/CH7/EX7.7.b/result7_7.txt
new file mode 100755
index 000000000..8b185977d
--- /dev/null
+++ b/764/CH7/EX7.7.b/result7_7.txt
@@ -0,0 +1,56 @@
+-->//(Threaded Joints) Example 7.7
+
+-->//Refer Fig.7.22 on page 238
+
+-->//Number of bolts N
+
+-->N = 6
+ N =
+
+ 6.
+
+-->//Load acting on the bracket P (kN)
+
+-->P = 50
+ P =
+
+ 50.
+
+-->//Height of bolt1 from the base l1 (mm)
+
+-->l1 = 300
+ l1 =
+
+ 300.
+
+-->//Height of bolt2 from the base l2 (mm)
+
+-->l2 = 200
+ l2 =
+
+ 200.
+
+-->//Height of bolt3 from the base l3 (mm)
+
+-->l3 = 100
+ l3 =
+
+ 100.
+
+-->//Eccentricity value l (mm)
+
+-->l = 250
+ l =
+
+ 250.
+
+-->//Maximum permissible tensile stress in bolts sigmaMax (N/mm2)
+
+-->sigmaMax = 100
+ sigmaMax =
+
+ 100.
+
+
+Area at the core cross-section(A) = 133.928571 mm2
+ \ No newline at end of file
diff --git a/764/CH7/EX7.7.b/solution7_7.sce b/764/CH7/EX7.7.b/solution7_7.sce
new file mode 100755
index 000000000..7152777ed
--- /dev/null
+++ b/764/CH7/EX7.7.b/solution7_7.sce
@@ -0,0 +1,22 @@
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution7_7.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_7.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the tensile force induced in appropriate bolt P1 (N)
+if ((l1 > l2) & (l1 > l3)) then
+ P1 = (P * 1000 * l * l1)/(2*((l1^2) + (l2^2) + (l3^2)))
+elseif ((l2 > l1) & (l2 > l3))
+ P1 = (P * 1000 * l * l2)/(2*((l1^2) + (l2^2) + (l3^2)))
+else
+ P1 = (P * 1000 * l * l3)/(2*((l1^2) + (l2^2) + (l3^2)))
+end
+//Calculate the core cross-section area of the bolt A (mm2)
+A = P1/sigmaMax
+//Choose proper diameter from Table 7.1
+//Print results
+printf('\nArea at the core cross-section(A) = %f mm2\n',A)
diff --git a/764/CH7/EX7.8.a/data7_8.sci b/764/CH7/EX7.8.a/data7_8.sci
new file mode 100755
index 000000000..1d39288b6
--- /dev/null
+++ b/764/CH7/EX7.8.a/data7_8.sci
@@ -0,0 +1,15 @@
+
+//(Threaded Joints) Example 7.8
+//Refer Fig.7.23 on page 239
+//Number of bolts N
+N = 2
+//Permissible tensile stress in the bolts sigmaMax (N/mm2)
+sigmaMax = 75
+//Load acting on the bracket P (kN)
+P = 20
+//Eccentricity value e (mm)
+e = 550
+//Distance of bolt1 from the C l1 (mm)
+l1 = 450
+//Distance of bolt2 from the C l2 (mm)
+l2 = 50
diff --git a/764/CH7/EX7.8.b/result7_8.txt b/764/CH7/EX7.8.b/result7_8.txt
new file mode 100755
index 000000000..eef3a3d02
--- /dev/null
+++ b/764/CH7/EX7.8.b/result7_8.txt
@@ -0,0 +1,49 @@
+-->//(Threaded Joints) Example 7.8
+
+-->//Refer Fig.7.23 on page 239
+
+-->//Number of bolts N
+
+-->N = 2
+ N =
+
+ 2.
+
+-->//Permissible tensile stress in the bolts sigmaMax (N/mm2)
+
+-->sigmaMax = 75
+ sigmaMax =
+
+ 75.
+
+-->//Load acting on the bracket P (kN)
+
+-->P = 20
+ P =
+
+ 20.
+
+-->//Eccentricity value e (mm)
+
+-->e = 550
+ e =
+
+ 550.
+
+-->//Distance of bolt1 from the C l1 (mm)
+
+-->l1 = 450
+ l1 =
+
+ 450.
+
+-->//Distance of bolt2 from the C l2 (mm)
+
+-->l2 = 50
+ l2 =
+
+ 50.
+
+
+Area at the core cross-section(A) = 455.284553 mm2
+ \ No newline at end of file
diff --git a/764/CH7/EX7.8.b/solution7_8.sce b/764/CH7/EX7.8.b/solution7_8.sce
new file mode 100755
index 000000000..5f257bd02
--- /dev/null
+++ b/764/CH7/EX7.8.b/solution7_8.sce
@@ -0,0 +1,24 @@
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution7_8.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_8.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the direct tensile force Dtensile (N)
+Dtensile = (P * 1000)/N
+//Calculate the tensile force due to bracket tilting tendency Ftensile (N)
+if (l1 > l2) then
+ Ftensile = (P * 1000 * e * l1)/((l1^2) + (l2^2))
+else
+ Ftensile = (P * 1000 * e * l2)/((l1^2) + (l2^2))
+end
+//Calculate the resultant tensile force on appropriate bolt res (N)
+res = Dtensile + Ftensile
+//Calculate the core cross-section area of bolt A (mm2)
+A = res/sigmaMax
+//Choose proper diameter from Table 7.1
+//Print results
+printf('\nArea at the core cross-section(A) = %f mm2\n',A)
diff --git a/764/CH7/EX7.9.a/data7_9.sci b/764/CH7/EX7.9.a/data7_9.sci
new file mode 100755
index 000000000..2820dafc4
--- /dev/null
+++ b/764/CH7/EX7.9.a/data7_9.sci
@@ -0,0 +1,15 @@
+
+//(Threaded Joints)Example 7.9
+//Refer Fig.7.24 on page 240
+//Number of bolts N
+N = 4
+//Load supported by the bracket P (kN)
+P = 25
+//Distance of bolt1 from C l1 (mm)
+l1 = 50
+//Distance of bolt2 from C l2 (mm)
+l2 = 200
+//Eccentricity value l (mm)
+l = 400
+//Permissible tensile stress in the bolt sigmaMax (N/mm2)
+sigmaMax = 50
diff --git a/764/CH7/EX7.9.b/result7_9.txt b/764/CH7/EX7.9.b/result7_9.txt
new file mode 100755
index 000000000..e9e4fc763
--- /dev/null
+++ b/764/CH7/EX7.9.b/result7_9.txt
@@ -0,0 +1,49 @@
+-->//(Threaded Joints)Example 7.9
+
+-->//Refer Fig.7.24 on page 240
+
+-->//Number of bolts N
+
+-->N = 4
+ N =
+
+ 4.
+
+-->//Load supported by the bracket P (kN)
+
+-->P = 25
+ P =
+
+ 25.
+
+-->//Distance of bolt1 from C l1 (mm)
+
+-->l1 = 50
+ l1 =
+
+ 50.
+
+-->//Distance of bolt2 from C l2 (mm)
+
+-->l2 = 200
+ l2 =
+
+ 200.
+
+-->//Eccentricity value l (mm)
+
+-->l = 400
+ l =
+
+ 400.
+
+-->//Permissible tensile stress in the bolt sigmaMax (N/mm2)
+
+-->sigmaMax = 50
+ sigmaMax =
+
+ 50.
+
+
+Area at the core cross-section(A) = 595.588235 mm2
+ \ No newline at end of file
diff --git a/764/CH7/EX7.9.b/solution7_9.sce b/764/CH7/EX7.9.b/solution7_9.sce
new file mode 100755
index 000000000..e9c2819d8
--- /dev/null
+++ b/764/CH7/EX7.9.b/solution7_9.sce
@@ -0,0 +1,24 @@
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution7_9.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data7_9.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the direct tensile stress in bolts Dtensile (N)
+Dtensile = (P * 1000)/N
+//Calculate the tensile force on appropriate bolt due to tilting tendency of the bracket Ftensile (N)
+if (l1 > l2) then
+ Ftensile = (P * 1000 * l * l1)/(2*((l1^2) + (l2^2)))
+else
+ Ftensile = (P * 1000 * l * l2)/(2*((l1^2) + (l2^2)))
+end
+//Calculate the resultant tensile force res (N)
+res = Dtensile + Ftensile
+//Calculate the core cross-section area of bolts A (mm2)
+A = res/sigmaMax
+//Choose proper diameter from Table 7.1
+//Print results
+printf('\nArea at the core cross-section(A) = %f mm2\n',A)