From 88b34823a23d9bef380e08c32d071b5e8fd010e9 Mon Sep 17 00:00:00 2001
From: Leahy <cpl11001@ad.engr.uconn.edu>
Date: Thu, 7 Dec 2017 17:02:52 -0500
Subject: [PATCH] 1

---
 ME3227_FinalProject.asv | 62 +++++++++++++++++++++++++++++++++++++++++
 ME3227_FinalProject.m   | 61 ++++++++++++++++++++++++++++++++++++++++
 2 files changed, 123 insertions(+)
 create mode 100644 ME3227_FinalProject.asv
 create mode 100644 ME3227_FinalProject.m

diff --git a/ME3227_FinalProject.asv b/ME3227_FinalProject.asv
new file mode 100644
index 0000000..14b6148
--- /dev/null
+++ b/ME3227_FinalProject.asv
@@ -0,0 +1,62 @@
+% Given Parameters: 
+% Shaft Material - AISI 4130 Q&T CD
+E = 29*10^6; %psi (Young's Modulus)
+sut = 118*10^3; %psi (ultimate strength)
+sy = 102*10^3; %psi (yield strength)
+p = 0.3; %(poisson's ratio)
+ga_w = 4; %lbf - make as input value
+ga_t = 50; %teeth number - input
+ga_a = 20; %deg - input
+ga_d = 4; %diametric pitch - input
+gc_w = 2; %lbf - make as input value
+gc_t = 25; %teeth number - input
+gc_a = 20; %deg - input
+gc_d = 4; %diametric pitch - input
+ss = 3000; %rpm (shaft speed)
+T = 3500; %lbf-in (torque)
+% Reliability for fatigue life = 0.99
+% nd, Safety factor for fatigue life = 2
+% nd, safety factor for all other criteria = 1
+% r/d = 0.02 - woodruff key
+gear_def = 0.01; % deflection at gears - from table 7-2
+bearing_slope = 0.001; %rad
+gear_slope = 0.005; %rad 
+d_a = ga_t/ga_d; % diameter of gear a
+d_c = gc_t/gc_d; % diameter of gear c
+Fct = T/(d_c/2); % tangential force on gear c
+Fat = T/(d_a/2); % tangential force on gear a
+Fcn = Fct * tand(gc_a); 
+Fan = Fat * tand(ga_a);
+
+% Find Reaction Forces - XY plane
+Rby = ((Fcn*3)+(Fan*9))/6;
+Rdy = Fan + Fcn - Rby; 
+
+% Calculate C1 and C2 - XY Plane
+c1 = (117*Fan-36*Rby+4.5*Fcn)/6;
+c2 = 4.5*Fan - 3*c1; 
+
+% Singularity Functions - x = 0 - 9 inches
+for i = 0:length(i)
+    if (0<=i) && (i<=3) 
+       EIdy_x = -Fan/2*(i)^2+c1;
+       EIy_x = -Fan/6*(i)^3+c1*i+c2;
+       
+   elseif (3<i) && (i<=6)
+       EIdy_x = -Fan/2*(i)^2+Rby/2*(i-3)^2+c1;
+       EIy_x = -Fan/6*(i)^3+Rby/6*(i-3)^2+c1*i+c2;
+   elseif (6<i) && (i<=9) 
+       EIdy_x = -Fan/2*(i)^2+Rby/2*(i-3)^2-Fcn/2*(i-6)^2+c1;
+       EIy_x = -Fan/6*(i)^3+Rby/6*(i-3)^2-Fcn/6*(i-6)^3+c1*i+c2;
+   end
+ end 
+ 
+
+
+
+
+
+% Diameter of the shaft due to slope at the bearings - Points B,D
+
+
+
diff --git a/ME3227_FinalProject.m b/ME3227_FinalProject.m
new file mode 100644
index 0000000..8d93779
--- /dev/null
+++ b/ME3227_FinalProject.m
@@ -0,0 +1,61 @@
+% Given Parameters: 
+% Shaft Material - AISI 4130 Q&T CD
+E = 29*10^6; %psi (Young's Modulus)
+sut = 118*10^3; %psi (ultimate stress)
+sy = 102*10^3; %psi (yield stress)
+p = 0.3; %(poisson's ratio)
+ga_w = 4; %lbf - make as input value
+ga_t = 50; %teeth number - input
+ga_a = 20; %deg - input
+ga_d = 4; %diametric pitch - input
+gc_w = 2; %lbf - make as input value
+gc_t = 25; %teeth number - input
+gc_a = 20; %deg - input
+gc_d = 4; %diametric pitch - input
+ss = 3000; %rpm (shaft speed)
+T = 3500; %lbf-in (torque)
+% Reliability for fatigue life = 0.99
+% nd, Safety factor for fatigue life = 2
+% nd, safety factor for all other criteria = 1
+% r/d = 0.02 - woodruff key
+gear_def = 0.01; % deflection at gears - from table 7-2
+bearing_slope = 0.001; %rad
+gear_slope = 0.005; %rad 
+d_a = ga_t/ga_d; % diameter of gear a
+d_c = gc_t/gc_d; % diameter of gear c
+Fct = T/(d_c/2); % tangential force on gear c
+Fat = T/(d_a/2); % tangential force on gear a
+Fcn = Fct * tand(gc_a); 
+Fan = Fat * tand(ga_a);
+
+% Find Reaction Forces - XY plane
+Rby = ((Fcn*3)+(Fan*9))/6;
+Rdy = Fan + Fcn - Rby; 
+
+% Calculate C1 and C2 - XY Plane
+c1 = (117*Fan-36*Rby+4.5*Fcn)/6;
+c2 = 4.5*Fan - 3*c1; 
+
+% Singularity Functions - x = 0 - 9 inches
+for x = linspace(0,9,10)
+   if (0<=x) && (x<=3) 
+       EIdy_x = -Fan/2*(x)^2+c1;
+       EIy_x = -Fan/6*(x)^3+c1*x+c2;
+   elseif (3<x) && (x<=6)
+       EIdy_x = -Fan/2*(x)^2+Rby/2*(x-3)^2+c1;
+       EIy_x = -Fan/6*(x)^3+Rby/6*(x-3)^2+c1*x+c2;
+   elseif (6<x) && (x<=9) 
+       EIdy_x = -Fan/2*(x)^2+Rby/2*(x-3)^2-Fcn/2*(x-6)^2+c1;
+       EIy_x = -Fan/6*(x)^3+Rby/6*(x-3)^2-Fcn/6*(x-6)^3+c1*x+c2;
+   end
+ end 
+ 
+
+
+
+
+
+% Diameter of the shaft due to slope at the bearings - Points B,D
+
+
+