Przykład 1

Dla podanej belki wyznaczyć częstości drgań wasnych i narvsować ich postacie. Sprawdzić warunek ortogonalności. Dane:

single-task-hero-img

Rozwiązanie

Stopień statycznej niewyznaczalności gdzie :
r - reakcje
p - przeguby
a - obwody zamknięte


Liczba stopni swobody dynamicznej

przyjmując
m =100 kg
𝑚1 =𝑚,𝑚2 =𝑚,𝑚3 =2.5𝑚

Ramę rozwiążemy METODĄ SIŁ.
Układ podstawowy metody sił

Wspólny stan jednostkowy

Reakcje
ΣMA=01212VB5=0 VB=0𝛿11=1EI(13222+252+13222)=7631EI Teraz równolegle robimy stan P dla wszystkich przypadków obciążania, obliczamy całkowanie δ1p osobno dla każdego przypadku, obliczamy x1 (rozwiązujemy równanie kanoniczne) i rysujemy wykresy ostateczne.


Delty dynamiczne z twierdzenia redukcyjnego 𝛿11=MOst1MP1EIdS𝛿11=1EI[56(21233811538)+13111]=77761EI 𝛿22=Most2MP2EIdS𝛿22=1EI[36(21.269951.2919)+26(21.269951.2919)]=93951EI 𝛿33=Most3MP3EIdS𝛿33=1EI[13222+26(222+236194+36192+42)+56(24361944019)]=1136571EI 𝛿12=Most1MP2EIdS𝛿12=1EI[36(21.239190+1.21538)26(21.239190+1.22338)]=791901EI𝛿21=𝛿12 𝛿13=Most3MP1EIdS𝛿13=1EI[56(21401913619)]=110571EI𝛿31=𝛿13𝛿23=MOst3MP2EIdS𝛿23=1EI[36(21.24895+123619)26(21.24895+1.24019)]=68951EI𝛿32=𝛿23 przyjmuję m =100 kg  m1 =1 m  m2 =1 m  m3 =2.5 m

Szukamy pierwiastków równania

m1𝛿111𝜔2 m2𝛿12 m3𝛿13 m1𝛿21 m2𝛿221𝜔2 m3𝛿23 m1𝛿31 m2𝛿32 m3𝛿331𝜔2=0 (mnożymy obustronnie przez El/m ) podstawiamy x =1𝜔2EIm

17776𝑥1(79190)52110571(79190)19395𝑥52(6895)1110571(6895)52113657𝑥 =0
7776x7919027557791909395x3419110576895284057x =0
Następnie należy metodą Sarrusa obliczyć wyznacznik z powyższej macierzy 3x3 i uprościć go do postaci wielomianu 3-go stopnia (jak poniżej). Najlepiej jeśli mamy do dyspozycji jakiś program, który może nam to uprościć, np. Mathcad.
1 𝑥3 +51.8167 𝑥2 89.4833 𝑥 +33.2632 =0
Pierwiastki równania x1=0.537x2=1.238x3=50.042 ponieważ x =1𝜔2 EI m 𝜔 =1x EIm
Częstości drgań własnych
𝜔1=1𝑥3=0.14136 EI m𝜔2=1𝑥2=0.899𝐸𝐼 m𝜔3=1𝑥1=1.3647𝐸𝐼 m Sztywność na zginanie  Dla I=2140108 m4 i modułu Younga dla stali E=200109 PaEI=EI=4280000.000Nm2 m=100 kg𝜔1=0.14136EIm=29.245rads𝜔2=0.899EIm=185.939rads𝜔3=1.365EIm=282.33rads
Wyznaczenie postaci drgań (m1𝛿11x)A11+m2𝛿12A21+m3𝛿13A31=0 m1𝛿12A11+(m2𝛿22x)A21+m3𝛿23A31=0 m1𝛿31A11+m2𝛿32A21+(m3𝛿33x)A31=0
Taki układ równań jak powyżej rozwiążemy trzy razy, albowiem raz za "x" podstawimy "x1", raz "x2" i raz "x3".
Po podstawieniu x, mamy układ 3 równań z trzema niewiadomymi A11, A21, A31. Ten układ równań zawsze ma nieskończenie wiele rozwiązań,
musimy więc pokazać jedno przykładowe rozwiązanie. To znaczy założymy za jedną z niewiadomych wartość równą 1, i zostanie nam układ 3 równań z dwoma niewiadomymi, gdzie jak rozwiążemy układ dwóch pierwszych równań to obliczymy pozostałe dwie niewiadome.

dla x =x1 =0.537 𝑧𝑎𝑘ł𝑎𝑑𝑎𝑚 A11 =1
0.476+0.416 A21+4.825A31=00.416+0.442 A211.7895A31=01.92980.716A21+49.288A31=0 Rozwiązuję układ dwóch pierwszych równań
A21=0.831 A31=0.027 I postać drgań


dla x =x2 =1.238 𝑧𝑎𝑘ł𝑎𝑑𝑎𝑚 A12 =1
0.2250.416A22+4.825A32=00.4160.259A221.7895A32=01.92980.716A22+48.587A32=0 Rozwiązuję układ dwóch pierwszych równań
A22=1.208 A32=0.0575 II postać drgań



dla x =x3 =50.042 𝑧𝑎𝑘ł𝑎𝑑𝑎𝑚 A13 =1
49.02860.416A23+4.825A33=00.41649.063A231.7895A33=01.92980.716A230.217A33=0 Rozwiązuję układ dwóch pierwszych równań
A23=0.378 A33=10.130 III postać drgań


Sprawdzenie warunków ortogonalności
m1A11A12+m2A21A22+m3A31A32=0 m1A12A13+m2A22A23+m3A32A33=0 m1A11A13+m2A21A23+m3A31A33=0
11m+1 m0.831(1.208)+2.5 m(0.027)(0.0575)0.00003 m011m+1 m0.831(1.208)+2.5 m(0.027)(0.0575)0.00003 m011m+1 m0.831(0.378)+2.5 m(0.027)10.12970.00213 m0
sprawdzenie daje pozytywny wynik