Temat zadania: Belkę wieloprzęsłową przegubową rozwiązać analitycznie. Wyznaczyć reakcje i sporządzić wykresy sił przekrojowych M, N, T. Metodą kinematyczną wyznaczyć wielkości R, Mα, Tα.
Sprawdzenie statycznej wyznaczalności i geometrycznej niezmienności układu
Warunek ilościowy statycznej wyznaczalności układu.
3t=e
Gdzie:
e – Ilość więzi elementarnych
t – Ilość tarcz
3*3=9
Warunek ilościowy jest spełniony.
Warunek jakościowy geometrycznej niezmienności układu.
Dla sprawdzenia geometrycznej niezmienności wszystkich tarcz skorzystano z twierdzenia o dwóch tarczach, tj. gdy dwie tarcze są połączone trzema nierównoległymi, niezbieżnymi więziami elementarnymi, to tworzą jedną sztywną tarczę.
Na podstawie powyższego twierdzenia najpierw rozpatrzono połączenie tarcz 1-0. Spełniają one warunek twierdzenia, więc traktujemy je jako całość. Następnie rozpatrzono połączenie tarczy 2 z układem (1-0). Układ również spełnia warunek twierdzenia. Ostatnie rozpatrzono połączenie tarczy 3 z układem (2-1-0). Połączenie to również spełnia warunek twierdzenia o dwóch tarczach. Zapisując kolejność sprawdzeń:
– z tw. o 2 tarczach;
2-(1-0) – z tw. o 2 tarczach;
3-(2-1-0) – z tw. o 2 tarczach.
Wniosek:
Układ jest statycznie wyznaczalny i geometrycznie niezmienny.
Rozwiązanie belki metodą bezpośrednią
Dane:
P1=8kN
P2=22kN
M=37kN
q=11kN/m
Siłę P2 oraz reakcję podpory C rozłożono na siły składowe o kierunkach zgodnych z kierunkami głównych osi kierunkowych. Wartości składowych obliczono z wykorzystaniem własności trygonometrycznych kąta.
Tg(45o) = 1;
Cos(45o) = Sin(45o) = $\frac{\sqrt{2}}{2}$;
Stąd otrzymano:
RCy=RCx=$\frac{R_{C}}{\sqrt{2}}$; P2y=P2x=$\frac{P_{2}}{\sqrt{2}}$=15,556kN
Siłę P1 położoną na wsporniku zredukowano do punktu A. Zgodnie z prawami redukcji otrzymano siłę skupioną równa co do wartości, kierunku i zwrotu sile P1 oraz moment prawoskrętny o wartości MP1=64kNm.
Obliczenie reakcji
Do obliczenia reakcji w podporach wykorzystano następujące równania:
∑M3(prawy) = 0 => Rd=?
∑M2(prawy) = 0 => RCy = RCx = $\frac{R_{C}}{\sqrt{2}}$ =?
∑MA = 0 => RB=?
∑Y = 0 => VA=?
∑X = 0 => HA=?
∑M3P= 0
4RD-2P2y= 0
4Rd= 31,112
Rd= 7,778 kN
∑M2P= 0
M-2Rcy-6P2y-8Rd= 0
2RC= 5,888
RC= 2,944 kN
∑MA= 0
MP1+50q-5RB+M-12RC-16P2y-18RD= 0
5RB= 506,780
RB= 101,356 kN
∑Y= 0
P1-VA-RB+10q-RC-P2y-RD= 0
VA= 5,922 kN
∑X= 0
-HA-RCx+P2x= 0
HA= 12,612 kN
Sprawdzenie:
∑M4= ?
MP1-16P1+16VA+11RB-110q+M+4RCy+2RD=
64-128+94,752+1114,916-1210+37+11,776+15,556= 0
Wniosek:
Sprawdzenie wykazało poprawność obliczeń.
Zestawienie wartości reakcji:
HA= 12,612 kN
VA= 5,922 kN
RB= 101,356 kN
RC= $\sqrt{2}*R_{\text{Cy}} = \ \sqrt{2}*2,944 =$ 4,163 kN
RD= 7,778 kN
Obliczenie sił przekrojowych
Siły osiowe (N):
(A – C)
x ϵ(0;12)
N(A-C)=HA=12,612 kN
(C – 4)
x ϵ(12;16)
N(C-4)=P2x=15,556 kN
(4 – D)
x ϵ(16;18)
N(4-D)=0
Siły tnące (T):
(A – B)
x ϵ(0;5)
T(A-B)=VA-P1-q*x= -2,078-11x [kN]
T (A) = -2,078 kN T(BL) = -2,078-55= -57,078 kN
(B – 2)
x ϵ(5;10)
T (B-2)=T(BL)+RB-q*(x-5)= -57,078+101,356-11(x-5)= 44,278-11(x-5) [kN]
T(BP)=44,278-11(5-5)= 44,278 kN T(2)=44,278-11(10-5)=44,278-55= -10,722 kN
(2 – C)
x ϵ(10;12)
T(2-C)=T(2)=-10,722 kN
(C – 4)
x ϵ(12;16)
T (C-4)=RD-P2y=7,778-15,556= -7,778 kN
(4 – D)
x ϵ(16;18)
T(4-D)=RD=7,778 kN
Momenty zginające (M):
(A – B)
x ϵ(0;5)
M(A-B)=MP1+VAx-P1x-qx2/2=64+5,922x-8x-11x2/2=64-2,078x-11x2/2 [kNm]
M(A)=64 kNm M(B)=64-10,39-137,5= -83,89 kNm
(B – 2)
x ϵ(5;10)
M(B-2)=M(A-B)+RB(x-5)=64-2,078x-11x2/2+101,356(x-5) [kNm]
M(B)= -83,39kNm M(2)=64-20,78-550+506,78= 0 kNm
(D - 4)
x ϵ(0;2)
M(D-4)= -RDx= -7,778x [kNm]
M(D)= 0 M(4)= -15,556 kNm
(4 – C)
x ϵ(2;6)
M(4-C)= M(D-4)+P2y(x-2)= -7,778x+15,556(x-2) [kNm]
M(4)= -15,556 kNm M(CP)= -46,668+62,224= 15,556 kNm
(C – 2)
x ϵ(6;8)
M(C-2)= M(4-C)-M+RC(x-6)= -7,778x+15,556(x-2)-37+2,944(x-6) [kNm]
M(CL)= -46,668+62,224-37= -21,444 kNm M(2) = -62,224+93,336-37+5,888=0 kNm
Obliczenie momentu maksymalnego
W przęśle (B – 2) siła tnąca zmienia znak pomiędzy punktami charakterystycznymi. Punkt przecięcia wykresu sił tnących z osią zerową wyznacza miejsce powstania maksymalnego momentu zginającego. Położenie punktu przecięcia na osi x można obliczyć przyrównując równanie T(x) do zera.
T(B-2)= 44,278-11(x-5) [kN]
T(B-2)=0 => 44,278-11(x-5)=0
11x=99,278
x=9,02527m
Wykorzystując obliczoną odległość x można obliczyć wartość maksymalną momentu zginającego.
M(B-2)= 64-2,078x-11x2/2+101,356(x-5) [kNm]
M(9,02527)= 64-18,75-448,01+407,99=5,23 kNm
Wykresy sił wewnętrznych
Wyznaczenie R, Mα, Tα za pomocą metody kinematycznej
Wyznaczenie R:
L = 0
110*5ϕ1-R*5ϕ1+8*8ϕ1-37*ϕ2+15,556*2ϕ3=0
Zgodność przesunięć:
10 ϕ1=2 ϕ2 2 ϕ2=4 ϕ3
ϕ2=5 ϕ1 ϕ3=0,5 ϕ2
ϕ3=2,5 ϕ1
5Rϕ1=550 ϕ1+64 ϕ1-185 ϕ1+77,78 ϕ1
5R ϕ1=506,78 ϕ1
R = 101,356 kN
Wyznaczenie Tα:
L = 0
55*2,5ϕ1-Tα*5ϕ1+55*2,5 ϕ2+8*8 ϕ1-37* ϕ3+15,556* ϕ4
Zgodność przesunięć:
ϕ2= ϕ1 2ϕ3=5 ϕ2 4 ϕ4=2 ϕ3
ϕ3=2,5 ϕ2 ϕ4=0,5 ϕ3
ϕ3=2,5 ϕ1 ϕ4=1,25 ϕ1
5Tα ϕ1=137,5 ϕ1+137,5 ϕ1+64 ϕ1-92,5 ϕ1+38,89 ϕ1
5Tα ϕ1=285,39 ϕ1
Tα=57,078 kN
Wyznaczenie Mα:
L = 0
Mα* ϕ1+55*2,5 ϕ1-37* ϕ2+15,556* 2ϕ3
Zgodność przesunięć:
2ϕ2=5 ϕ1 4 ϕ3=2 ϕ2
ϕ2=2,5 ϕ1 ϕ3=0,5 ϕ2
ϕ3=1,25 ϕ1
Mα ϕ1=137,5 ϕ1-92,5 ϕ-1+38,89 ϕ1
Mα ϕ1=83,89 ϕ1
Mα=83,89 kNm
Politechnika Wrocławska
Wydział Budownictwa Lądowego i Wodnego
Instytut Inżynierii Lądowej
Zakład Dynamiki Budowli
Rok Studiów: 2
Projekt nr 1
Z Podstaw Statyki Budowli
Belka wieloprzęsłowa przegubowa
Opracował: Michał Wolański
Nr Albumu 183328
Grupa PN 13:15
Prowadzący: Doc. Dr Inż. Marek Kopiński
Rok akademicki:
2011/2012
Zestawienie wartości:
P1=8 kN
P2=22 kN
P2x=P2y=15,556 kN
M=37 kNm
q=11 kN/m
HA= 12,612 kN
VA= 5,922 kN
RB= 101,356 kN
RC= $\sqrt{2}*R_{\text{Cy}} = \ \sqrt{2}*2,944 =$ 4,163 kN
RD= 7,778 kN