symbolic - lengthy = sqrt(XX(i 1 XX(i)^2(YY(i 1-YY(i)^2%ds...

Info iconThis preview shows pages 1–2. Sign up to view the full content.

View Full Document Right Arrow Icon
clear c %everything is in millimeters % a = .019285; h = 500; m = -h/2; n = h/2; n th = 10; ts = 10; t c = a*(h/2)^2; c ss = 50; s %spline x = sym ('x'); eqn = sym (sqrt(1+(2*a*x)^2)); length = double (int (eqn, m, n)); %length l eqn2 = sym (-a*x^2 + c); %ezplot (eqn2, [-260, -240]); area = double (int (eqn2, m, n)); %area a centp = (3/5)*h; %centroid of the parabola c areap = ts*length; %area of parabola a %web areaw = h*th; centw = th/2; c %spar areas = 5000; cents = th; c %centroid of everything centroid = (areap*centp+areaw*centw+areas*cents)/(areap+areaw+areas) c %second moment of area %spline X=[-h/2,0,h/2]; Y=[0,c,0]; XX = -h/2:1:h/2; YY = spline(X, Y, XX); plot(XX,YY); p Ixspline=0; I for i = 1:1:size(XX) length1 = sqrt((XX(i+1) - XX(i))^2+(YY(i+1)-YY(i))^2); area1 = length1*ts; Ixspline = Ixspline + area1*((YY(i+1)-YY(i))/2)^2; end Ixspline; I Iyspline=0; I
Background image of page 1

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
for i = 1:1:size(XX)
Background image of page 2
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: lengthy = sqrt((XX(i+1) - XX(i))^2+(YY(i+1)-YY(i))^2); %ds areay = lengthy*ts; Iyspline = Iyspline + areay*((XX(i+1)-XX(i))/2-centp)^2; end Iyspline; I %web Ixweb = h*th^3/12; Iyweb = th*h^3/12; I %spar Ixspar = ss^4/12; Iyspar = Ixspar; I Ix = (Ixspline)+(Ixweb)+(2*(Ixspar+areas*(h/2)^2)) Iy = (Iyspline+areap*(centp-centroid)^2)+(Iyweb+areaw*(centw-centroid)^2)+ (Iyspar+areas*(cents-centroid)^2) ( %shear center %chart X = .33*area/h X sigma = 0; phi0 = 0; ty = 1; t for i = 1:1:size(XX,2)-1 sigma = sigma + (abs(XX(i))*ts*sqrt((XX(i+1) - XX(i))^2+(YY(i+1)-YY(i))^2)); end e H = [0, 0]; I = [h/2, -h/2]; JJ = -h/2:1:h/2; J WW = spline(I, H, JJ); W for i = 1:1:size(JJ,2)-1 sigma = sigma + abs(JJ(i))*th; end e phi1 = -(ty/Ix)*sigma + phi0; p tau1 = phi1/10; t...
View Full Document

This note was uploaded on 09/18/2011 for the course EC 101 taught by Professor Staff during the Spring '09 term at INSA Toulouse.

Page1 / 2

symbolic - lengthy = sqrt(XX(i 1 XX(i)^2(YY(i 1-YY(i)^2%ds...

This preview shows document pages 1 - 2. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online