cz1

WSZYSTKIE TEKSTY SA TAKEI SAME JAK W KAZDEJ INNEJ PRACY - TAKZ E JEŚLI PRZEPISZESZ TO ZYWCEM TO JEST 100% PRAWDOPODOBNE ZE ZACZNIE CIE PYTAC O WSZYSTKO JAK I SKAD OBLICZYLES. Staraj się żeby praca nie była podobna, zmieniaj tekst, zdania, rysunki, układ strony, ale generalnie warto cos wiedzcie jak jest to policzone bo Ci ludzie zawsze i każdego potrafią zapytać ukradkiem o jakąś głupotę, im mniej człowiek wie tym go bardziej pytają.

Część I - zadanie 1 - wyprowadzenie równań równowagi płynu.

Rozważamy dowolną przestrzeń
wypełnioną płynem idealnym

- wnętrze obszaru

- powierzchnia ograniczająca obszar

ds - elementarny wycinek powierzchni dV - elementarny wycinek wnętrza obszaru

ΔFp - jednostkowa siła powierzchniowa

ΔFv - jednostkowa siła objętościowa

- wektor normalny do wycinka ds.

- gęstość sił objętościowych

- gęstość sił powierzchniowych

Siły powierzchniowe zależą od usytuowania rozpatrywanej powierzchni, nie są zatem wielkością stałą i niezależną. Niezależnie od doboru powierzchni opisywał będzie stan cieczy tensor napręzeń- charakteryzujący naprężenia w danym punkcie. Na jego wartości nie ma wpływu sposób doboru powierzchni.

- tensor naprężeń

Ze względu na równości odpowiednich naprężeń stycznych, macierz ta jest macierzą symetryczną, czyli:
, ponieważ
.

Rysunki przedstawiają elementarny wycinek obszaru Ω- czyli d Ω, oraz jego naprężenia i współrzędne wektora

Są to makroskopowe ujęcia punktu obszaru. Można spostrzec ze na jego powierzchnię działają siły powierzchniowe zależne od doboru powierzchni. Ściana powierzchni opisywana jest przez wektor

Ogólny stan naprężenia w danym punkcie przedstawia tensor naprężeń σ Siła powierzchniowa będzie zatem równa:

- jest to wektor [N/m] = [Pa]

(R³) (R^3x3) (R³)

Znalezienie siły wypadkowej wymaga scałkowania po powierzchni:

; [N]

Analizując budowę tensora σ zauważymy, że:

Tensor σ tensor kulisty tensor dewiatorowy - σ_D

Gdzie:
; oraz
- ślad macierzy σ

Tensor kulisty to

Czyli ostatecznie:

Tensor kulisty odpowiedzialny jest za opis matematyczny wszechstronnego ściskania/rozciągania. Dewiator jest wyrażeniem związanym ze zmianą postaci ciala - ścinaniem. Ponieważ przy rozpatrywaniu płynu nie możemy mowić o rozciąganiu, a σ₀ ma kierunek naprężeń rozciągających, to ma on zawsze wartość ujemna. Wartość tę nazywamy ciśnieniem.

;
;

Po uwzględnieniu założenia, że siła powierzchniowa w każdym kierunku jest jednakowa i skierowana prostopadle do powierzchni, to dla danej chwili:

Jest to całka wektoru

Rozważając kolejno siły objętościowe:

[N/m³] - gęstość siły objętościowej.

Są to siły wewnętrzne, wzajemne. Obliczamy ich wypadkową dla danej chwili:

Siła wypadkowa w całym obszarze to
:

Na mocy Twierdzenia Gaussa-Ostrogradskij-ego możemy zastąpić całkę powierzchniową całka objętosciową w poniższy sposób:

Aby zachowana była równowaga to:
dla każdego obszaru Ω czyli:

Musi tak być dla każdego obszaru Ω 0 (ponieważ równie dobrze wymagać można zachowania równowagi i tego obszaru). Zależność
= 0 nazywamy równaniem równowagi ośrodka ciągłego, lub równaniem Eulera. Obowiązuje ono dla płynu nielekkiego.

, a wtedy:

Możemy zatem napisać równanie równowagi dla płynu:

Ponieważ:
, oraz
to ostatecznie:

Dla zachowania pełnej równowagi spełnione musi być też równanie momentów.

Po przekształceniach analogicznych do tych z warunku na równowagę sił, możemy zauważyć że jeżeli siłą wypadkowa będzie równa zeru, to niezależnie od tego na jakim ramieniu będzie działała, również da moment równy zeru. A zatem warunkiem wystarczającym do spełnienia aby ośrodek był w równowadze, jest zależność:

Jest to warunek równowagi płynu doskonałego.

Zadanie 2 - cysterna

α%=	11	%		α=	6,3	°
L=	6,4	[m]
D=	1,9	[m]
d=	1	[m]
ρcieczy=	1000	[kg/m^3]
a=	-0,24	g		Współrzędne klapy:
g=	9,81	[ms^2]		Xk=	6,4	klapa tylna
Pa=	101325	[Pa]		Yk=	0,95

Znajduję siły objętościowe działające na ciecz w poruszającej się cysternie:

Siła ciężkości: fg =ρ g

Sila bezwładności fb = - ρ a ( skierowane przeciwnie do przyspieszenia cysterny - ale ponieważ przyspieszenie ma wartość ujemną, to sila bezwładności ma zwrot taki jak na rysunku - przeciwnie do osi ox)

[N/m³] - jest to wektor wypadkowy sił objętościowych, obliczam jego współrzędne (podstawiam tylko wartości bezwzględne składowych, ich zwrot jest uwzględniony w znakach minus i plus):

fx=	-1282
fy=	-9751
fz=	0

Mając współrzędne możemy usytuować wektor w płaszczyźnie cysterny. Wektor wskazuje kierunek najszybszego wzrostu cisnienia, po przeciwnej stronie jest kierunek spadku. W rogu, wskazywanym przez kierunek spadku ciśnienia znajduje się pęcherzyk powietrza, panuje w nim ciśnienie atmosferyczne.

Linie prostopadłe do wektora wypadkowego to izobary - linie stałego ciśnienia.

Wiemy, że:

Jest tak ze względu na to ze pochodne czastkowe dp/dx; dp/dy i dp/dz maja stale wartości

Dzieki temu możemy zastąpić przyrosty dx dy i dz wartościami x y i z.

I wtedy: p(x,y,z) =

p(x,y,z) =

-1282

x +

-9751

y +

0

z + c

C - obliczamy z warunku brzegowego na ciśnienie atmosferyczne w punkcie gdzie znajduje się pęcherzyk powietrza.

Współrzędne punktu pęcherzyka to:

Xp=	6,4
Yp=	1,9

Czyli p( x_p, y_p) = 101325 = -1282 x_p - 9751y_p + c

Po obliczeniu:

c=

128056

[Pa]

Równanie ciśnienia ma następującą postać:

p(x,y,z) =

-1282

x +

-9751

y +

128056

Aby znaleźć siłę naporu wody na klapę należy obliczyc ciśnienie w punkcie środka ciężkości tej klapy, a następnie pomnożyc je przez pole powierzchni tej klapy:

Współrzędne klapy:
Xk=	6,4	klapa tylna	p(xk, yk)=	110589	[Pa]
Yk=	0,95
Aklapy=πd^2/4=		0,79	[m^2]
F=A*p(xk, yk)=		86259,13	[N]	-odpowiedź

Zadanie 2 - pława

H0=	1	[m]
m0=	190	[kg]
δ=	0,003	[m]
ρstali=	7400	[kg/m^3]
ρwody=	1000	[kg/m^3]
g=	9,81	[m/s^2]

Najpierw obliczam objętość zanurzoną, masę i siłę wyporu pontonu, i zestawu pontonów:

Równanie równowagi układu:

35041	a3	=	7573	a2+	1864
35041	a3	-	7573	a2-	1864	=0

Za pomocą Excela obliczam wartości dla kolejnych wartości parametru a i znajduję pierwsze dodatnie rozwiązanie. Dokładne rozwiązanie to a=0,459 m, zatem przyjmuję do dalszych obliczeń wartość a = 0,46 [m]

Następnie trzeba obliczyc położenie środka ciężkości pontonów, zestawu pontonów i płąwy oraz położenie srodka wyporu. Pontony są symetryczne względem płaszczyzny poziomej, więc ich środek ciężkości jest w połowie wysokości pontonu

xc pont=	-0,23
M=	163
xc plawy=Ho=	1
m0=	190

xc=

0,431

[m]

Środek wyporu pontonu obliczam jako środek ciężkości objętości zanurzonej pontonu, przy czym zgodnie z założeniem, ponton przyjmuję zanurzony do polowy swojej wysokości.

-0,336 [m]

Odległość pomiędzy środkiem wyporu i środkiem ciężkości:

lm=|xcxw|=

0,767

[m]

Obliczam momenty bezwładności przekroju względem jego głównych centralnych osi bezwładności. Przekrój zanurzeniowy pontonu jest prostokątem o bokach a i 2a

0,0298 [m⁴]

0,00746 [m⁴]

Obliczam momenty bezwładności układu pontonów względem głownej centralnej osi bezwładności (ze względu na symetrię promienista układu momenty względem innych osi przechodzących przez środek symetrii będą takie same)

Przyjmuję promień rozstawu pław r = 2 m ;

A - pole przekroju zanurzeniowego pontonu = 2a²

2(0,0298+0,00746+0,423*4) = 3,46 [m4]

Rysunek:

(narysuj oczywiście przekroje prostokątne - takie jak masz powyżej)

Możemy teraz obliczyć wysokość metacentryczną:

= 3,46 [m⁴]/ 3,57a³ - 0,767 [m] = 9,96 - 0,767 =9,193 [m] >>>2 [m] - warunek stateczności pławy postawiony w zadaniu jest spełniony

1

---