close all
clear all


Est = 29e6;   %psi  - Steel
Eal = 10e6;   %psi -  Al
Efrp = 2.8e6;   %psi  - G10


Gst = 11.458e6;  %psi - Steel
Gal = 5.8e6; %psi - Al  
Gfrp = 1.5e6;  %psi - G10 (estimated)

E = Eal;
G = Gal;

L = 2;  %meters

P = 800;  %lbs per square meter
P = P*.0254^2;  %psi
P = 0.5;

L = 2/.0254;  %inches
b = 5.5;  %inches
a = 26-5.5;  %inches

do = 2.0;  %square shaft outer diameter
di = do-2*(1/4); %square shaft inner diameter
Ishaft = (do^4-di^4)/12;

t = 0.5;  %Plate thickness

T = 5;  %N-m/m/deg
T = 60*T*L*.0254/(4.45*.0254);  %in-lbs at 60 degrees


disp(['E = ' num2str(E) ' psi   t = ' num2str(t) 'in']);

%Chordwise Bending
w = P;  %Loading for 1" wide section
disp('Chordwise Bending');

I = (t^3)/12;  %I for 1" wide slice of span.

dtip = w*a^4/(8*E*I)
M_max = w*a^2/2
tau_max = M_max*(t/2)/I


%Spanwise Bending
w = P*(a+b);  %Loading per unit length (lbs/in)
disp('Spanwise Bending');

I = (1/12)*(do^4-di^4); 

dmax = 5*w*L^4/(384*E*I)
M_max = w*L^2/8
tau_max = M_max*(do/2)/I


%Torsion
disp('Torsion Displacement')

J = 2.25*((do/2)^4-(di/2)^4);
%J = 0.5*0.5^4; %for 1" diameter rod
theta = T*L/(J*G);
theta = theta*180/pi
tau_max = T*(do/2)/J
