Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd Logo
Upload

Welcome to Scribd!

  • 10 views

    %plate With Stiffener Code

    Uploaded by

    ShashankGupta
    Assignment AE608A

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd

    %plate With Stiffener Code

    Uploaded by

    ShashankGupta
    10 views4 pages
    Assignment AE608A

    Original Title

    upadhay2

    Copyright

    © © All Rights Reserved

    Available Formats

    PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    Assignment AE608A

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    10 views4 pages

    %plate With Stiffener Code

    Uploaded by

    ShashankGupta
    Assignment AE608A

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    of 4

    %Plate with Stiffener Code

    clear all;
    %Calculation of Neutral Point
    t=0.005;
    n=2;
    E=[70e9 70e9];
    a=[t 0.01];
    b=[2 0.005];
    y1=[0 0];
    z1=[0.015 0.005];
    NY=0;
    NZ=0;
    D=0;
    for i=1:n
    NY=NY+E(i)*a(i)*b(i)*y1(i);
    NZ=NZ+E(i)*a(i)*b(i)*z1(i);
    D=D+E(i)*a(i)*b(i);
    i=i+1;
    end
    ynp=NY/D;
    znp=NZ/D;
    %Calculation of Area Moment of Inertias
    E1=0;
    A1=0;
    for i=1:n
    E1=E1+E(i)*a(i)*b(i);
    A1=A1+a(i)*b(i);
    i=i+1;
    end
    Estar=E1/A1;
    Iyy=[0 0];
    Izz=[0 0];
    Iyz=[0 0];
    for i=1:n
    Iyy(i)=Iyy(i)+a(i)*b(i)^3/12+a(i)*b(i)*(z1(i)-znp)^2;
    Izz(i)=Izz(i)+b(i)*a(i)^3/12+a(i)*b(i)*(y1(i)-ynp)^2;
    Iyz(i)=Iyz(i)+a(i)*b(i)*(z1(i)-znp)*(y1(i)-ynp);
    i=i+1;
    end
    Nyy=0;
    Nzz=0;
    Nyz=0;
    for i=1:n
    Nyy=Nyy+E(i)*Iyy(i);
    Nzz=Nzz+E(i)*Izz(i);
    Nyz=Nyz+E(i)*Iyz(i);
    i=i+1;
    end
    Iyystar=Nyy/Estar;
    Izzstar=Nzz/Estar;
    Iyzstar=Nyz/Estar;
    eps=1e-20; %Minimum Error

    1
    if Iyystar<eps
    Iyystar=0;
    end
    if Izzstar<eps
    Izzstar=0;
    end
    if Iyzstar<eps
    Iyzstar=0;
    end
    pp=pi;
    N=2%
    Iyycs=Iyystar;
    Zcs=znp;
    As=1%
    q=1%
    nu=1%
    D1=E*t^3/(12);
    nn=N^2;
    G=0.5*E/(1+nu);
    D2=G*t^3/3;
    p=1;
    for i=1:N
    q=1;
    for j=1:N
    K=
    0.5*((D1*((i*pp)^4+(j*pp/2)^4)+2*D2*(i*j*pp^2/2)^2)+E*(Iyycs
    +Zcs^2)*As*(i*pp)^4*(sin(j*pp/2)^2));
    F=((2*q)/(i*j*pp*pp))*(1-cos(i*pp))*(1-cos(j*pp));
    Al(p,q)=F/K(1);
    q=q+1;
    end
    p=p+1;
    end

    x=linspace(0,1);
    y=linspace(0,2);
    [X,Y]=meshgrid(x,y);
    u=0;
    for i=1:N
    for j=1:N
    u=u+Al(i,j)*sin(i*pp.*X).*sin(0.5*j*pp.*Y);
    end
    end

    Um=max(u(:));
    exx=0;
    eyy=0;
    exy=0;
    for i= 1:N
    for j=1:N
    exx=exx+0.5*t*Al(i,j)*(i*pp)^2*sin(i*pp.*X).*sin(0.5*j*pp.*Y);
    eyy=eyy
    +0.5*t*Al(i,j)*(0.5*j*pp)^2*sin(i*pp.*X).*sin(0.5*j*pp.*Y);

    2
    exy=exy-0.5*t*2*Al(i,j)*(0.5*j*pp)*(i*pp)*cos(i*pp.*X).*cos(0.5*j*pp.*Y);
    end
    end
    Sxx=(exx-nu*eyy)*70e9;
    Syy=(eyy-nu*exx)*70e9;
    Sxy=exy*(0.5*70e9/(1+nu));
    surf(X,Y,u);
    S1=max(Sxx(:));
    S2=max(Syy(:));
    S3=0; % since the Values are almost close to
    Zero
    S=[S1 S3;S3 S2];
    Eig=eig(S);
    R=rms(Eig);

    N =

    As =

    q =

    nu =

    3
    Published with MATLAB R2016a

    You might also like