From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 764/CH7/EX7.18.b/solution7_18.sce | 60 +++++++++++++++++++++++++++++++++++++++
 1 file changed, 60 insertions(+)
 create mode 100755 764/CH7/EX7.18.b/solution7_18.sce

(limited to '764/CH7/EX7.18.b/solution7_18.sce')

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)
-- 
cgit