Scribd Logo
Upload

Welcome to Scribd!

  • 13 views

    Belka

    Uploaded by

    Paweł Czechowski

    Copyright:

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

    Belka

    Uploaded by

    Paweł Czechowski
    13 views5 pages

    Copyright

    © © All Rights Reserved

    Available Formats

    DOCX, PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    Copyright:

    © All Rights Reserved
    Download as DOCX, PDF, TXT or read online from Scribd
    Download as docx, pdf, or txt
    13 views5 pages

    Belka

    Uploaded by

    Paweł Czechowski

    Copyright:

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

    % MES_belka

    E=2e11; % moduł Younga [N/m2]

    mb=55; % masa belki [kg/mb]

    I0=15e-5; % porównawczy moment bezwładności [m4]

    le=3; % liczba elementów

    lw=4; % liczba wezłów

    N=8; % liczba stopni swobody

    ien=[1 2;2 3;3 4]' % numery wezłów w elementach

    Ee=[E,E,E]; % moduły Younga materiału elementów

    Ie=[I0,I0,I0]; % momenty bezwładności elementów

    Le=[2,2,2] % długości elementów

    mbe=[mb,mb,mb] % masa elementów [kg/mb]

    % Tablica alokacji

    lm=tab_aloc(2,2,lw,ien,le)

    % Budowa macierzy sztywności elementów

    k1=elstif_beam(1,Ee,Ie,Le)

    k2=elstif_beam(2,Ee,Ie,Le)

    k3=elstif_beam(3,Ee,Ie,Le)

    % Agregacja globalnej macierzy sztywności konstrukcji k

    k=zeros(N,N);

    k=agre_matrix(k,1,k1,2,2,le,lw,lm);

    k=agre_matrix(k,2,k2,2,2,le,lw,lm)

    k=agre_matrix(k,3,k2,2,2,le,lw,lm)

    % Budowa macierzy bezwładności elementów

    m1=elmas_beam(1,Le,mbe)

    m2=elmas_beam(2,Le,mbe)

    m3=elmas_beam(3,Le,mbe)

    %
    % Agregacja globalnej macierzy bezwładności konstrukcji m

    m=zeros(N,N);

    m=agre_matrix(m,1,m1,2,2,le,lw,lm);

    m=agre_matrix(m,2,m2,2,2,le,lw,lm)

    m=agre_matrix(m,3,m3,2,2,le,lw,lm)

    % Uwzględnienie warunków brzegowych

    % w macierzy sztywności konstrukcji

    k(7,:)=[];k(:,7)=[]; % usunięto 7. wiersz i kolumnę macierzy k(8x8)

    k(2,:)=[];k(:,2)=[]; % usunięto 2. wiersz i kolumnę macierzy k(7x7)

    k(1,:)=[];k(:,1)=[]; % usunięto 1. wiersz i kolumnę macierzy k(6x6)

    % w macierzy bezwładności konstrukcji

    m(7,:)=[];m(:,7)=[]; % usunięto 7. wiersz i kolumnę macierzy k(8x8)

    m(2,:)=[];m(:,2)=[]; % usunięto 2. wiersz i kolumnę macierzy k(7x7)

    m(1,:)=[];m(:,1)=[]; % usunięto 1. wiersz i kolumnę macierzy k(6x6)

    % Rozwiązanie zagadnienia własnego det(k-omega2*m)=0

    [A,Omega]=eig(k,m)

    % Wyznaczenie częstości kołowych omi [rad/sec] i fizycznych fi [Hz]

    om1=sqrt(Omega(1,1))

    f1=om1/2/pi

    om2=sqrt(Omega(2,2))

    f2=om2/2/pi

    om3=sqrt(Omega(3,3))

    f3=om3/2/pi

    om4=sqrt(Omega(4,4))
    f4=om4/2/pi

    om5=sqrt(Omega(5,5))

    f5=om5/2/pi

    % Sprawdzenie ortogonalności drgań i normalizacji wektorów własnych

    M=A'*m*A

    % Otrzymano macierz jednostkową

    ((3*2*10^(11)*15*10^(-5))/(0.2*55*6^4))^(0.5)*2/pi

    Sprawdzenie
    50,583−50,441=0,14200
    0,14200
    ∗100 %=0,28%
    50,583
    Błąd wyniku na poziomie 0,28% .
    Stwierdza się prawidłowość wykonania zadania.

    form1=A(:,1);

    form1=form1/A(1,1)

    form2=A(:,2);

    form2=form2/A(1,2)

    form3=A(:,3);

    form3=form3/A(1,3)

    You might also like