summaryrefslogtreecommitdiff
path: root/3776/CH6/EX6.5/Ex6_5.sce
diff options
context:
space:
mode:
Diffstat (limited to '3776/CH6/EX6.5/Ex6_5.sce')
-rw-r--r--3776/CH6/EX6.5/Ex6_5.sce30
1 files changed, 30 insertions, 0 deletions
diff --git a/3776/CH6/EX6.5/Ex6_5.sce b/3776/CH6/EX6.5/Ex6_5.sce
new file mode 100644
index 000000000..551491f53
--- /dev/null
+++ b/3776/CH6/EX6.5/Ex6_5.sce
@@ -0,0 +1,30 @@
+clear
+//Given
+//We will divide this into three parts
+F = 8 //k - force applied
+d = 16 //inch -distance
+l_1 = 1 //in
+l_2 = 3 //in
+b_1 = 4 //in
+b_2 = 1 //in
+A_1 = l_1* b_1 //in2 - area of part_1
+y_1 = 0.5 //in com distance from ab
+A_2 =l_2*b_2 //in2 - area of part_1
+y_2 = 2.5 //in com distance from ab
+A_3 = l_2*b_2 //in2 - area of part_1
+y_3 = 2.5 //in com distance from ab
+
+y_net = (A_1*y_1 +A_2*y_2 + A_3*y_3)/(A_1+A_2+A_3) //in - The com of the whole system
+c_max = (4-y_net) //in - The maximum distace from com to end
+c_min = y_net //in - the minimum distance from com to end
+I_1 = b_1*(l_1**3)/12 + A_1*((y_1-y_net)**2) //Parallel axis theorm
+I_2 = b_2*(l_2**3)/12 + A_2*((y_2-y_net)**2)
+I_3 = b_2*(l_2**3)/12 + A_2*((y_2-y_net)**2)
+I_net = I_1 + I_2 + I_3 //in^4 - the total moment of inertia
+M_c = F*d*c_max
+stress_cmax = M_c/I_net //ksi - The maximum compressive stress
+
+M_t= F*d*c_min
+stress_tmax = M_t/I_net //ksi - The maximum tensile stress
+printf("\n The maximum tensile stress %0.3f ksi",stress_tmax )
+printf("\n The maximum compressive stress %0.1f ksi",stress_cmax)