7. Modele matematyczne elementów mikrohydraulicznych
7.3. Opracowanie modelu matematycznego mikrozaworu dławiącego
ROZWIĄZANIE MODELU
7.3.1. MODEL MATEMATYCZNY WYBRANEJ KONSTRUKCJI MIKROZAWORU DŁAWIĄCEGO
Niektóre odmiany konstrukcyjne zaworów hydraulicznych są wyposażone w gniazdo oraz element domykający w kształcie stożka – tzw. grzybek. Para grzybek– gniazdo (rys. 7.8), idealnie nadaje się do mikrozaworów dławiących. Elementy mają w tym wypadku proste kształty geometryczne, co ułatwia ich wykonanie. Zawory o tej konstrukcji odznaczają się dużą szczelnością w fazie zamknięcia, co w mikroprzepły-wach wydaje się istotne. Ponadto przemieszczenia elementu domykającego są bardzo małe, co przy zaworach dławiących umożliwia idealne kojarzenie ich z piezoelek-trycznymi elementami sterującymi. Elementy piezoelektryczne cechują się małymi wymiarami, dużą szybkością działania oraz małym poborem mocy, w związku z czym idealnie nadają się do automatyzacji sterowań dławieniowych układów mikrohydrau-licznych [7.7, 7.8].
Opracowanie modelu matematycznego dla nastawnego zaworu dławiącego spro-wadza się do określenia zależności między trzema wielkościami: natężeniem przepły-wu cieczy płynącej przez zawór Q, różnicą ciśnień po obu stronach zaworu dławiące-go ∆p oraz parametrem wynikającym z nastawy zaworu dławiącedławiące-go. W przypadku układu grzybek–gniazdo takim parametrem będzie przemieszczenie elementu domy-kającego względem gniazda z.
Związek między natężeniem przepływu cieczy przez zawór dławiący a różnicą ciśnień można opisać wzorem [7.5]:
2 1 2
p c Q c Q
∆ = + (7.37)
gdzie c1, c2 – współczynniki zależne od przemieszczenia elementu domykającego. Dla przepływów laminarnych, dla których liczba Reynoldsa jest mała, człon ze współczynnikiem c1 traci znaczenie i zależność jest liniowa, natomiast dla przepły- wów burzliwych człon ze współczynnikiem c2 traci znaczenie i zależność jest parabo-liczna. Ponieważ współczynniki c1 i c2 są zależne nie tylko od geometrii układu, ale również od przemieszczenia grzybka, postać zależności określona wzorem (7.37) nie jest zbyt użyteczna.
Liczba Reynoldsa, która określa charakter przepływu w szczelinie między gniaz-dem a stożkiem, dla układu z rysunku 7.8 jest definiowana jako:
v 2 π( sin cos ) h d Q Re v D z α α v = = − (7.38)
Rys. 7.8. Układ gniazdo–grzybek w kształcie stożka
gdzie: v – lepkość kinematyczną cieczy, m2/s, dh – średnica hydrauliczna najmniejsze-go przekroju przepływowenajmniejsze-go, m, v – prędkość cieczy w tym przekroju, m/s. D oraz
α mają wymiary jak na rysunku 7.8. Wzór (7.38) został szerzej opisany w rozdziale 5.
Warto zaznaczyć, że dla układu z rysunku 7.8 liczba Reynoldsa w wypadku stałej lepkości cieczy jest silnie zależna od natężenia przepływu, natomiast w bardzo małym stopniu zależna od przemieszczenia grzybka. W związku z tym w mikrohydraulice, gdzie są małe przepływy, liczba Reynoldsa przyjmuje również małe wartości i należy się spodziewać przepływów laminarnych bądź przejściowych.
Dla elementów klasycznej hydrauliki, gdzie duże są przepływy, a co za tym idzie również liczby Reynoldsa, bardzo popularna jest następująca zależność określająca natężenie przepływu w szczelinie dławiącej:
2 ( )
Q=µA z ρ ∆p (7.39)
gdzie: µ – bezwymiarowy współczynnik przepływu, ρ – gęstość cieczy, kg/m3, A(z) – pole najmniejszego przekroju, m2, przez który płynie ciecz i które można przed-stawić w postaci równania:
( ) π( sin cos ) sin
A z = D−z α α z α (7.40)
Współczynnik przepływu µ jest przyjmowany często w dużym uproszczeniu jako wartość stała. Założenie to jest przeważnie prawdziwe dla elementów klasycznej hy-drauliki, w której w oporach miejscowych występują przepływy burzliwe. Nie-uwzględnienie zmienności współczynnika przepływu dla małej liczby Reynoldsa po-wodować będzie jednak duże rozbieżności między przebiegiem modelowym a doświadczalnym.
W rozdziale 5 zaprezentowany został jeden ze sposobów wprowadzenia zmiennego współczynnika przepływu. Założone zostały dwa obszary zjawisk przepływowych dla
układu grzybek–gniazdo. W obszarze liczb Reynoldsa mniejszych od krytycznej, współczynnik przepływu był następującą funkcją liczby Reynoldsa:
(Re) a Re
µ
= (7.41)gdzie a – zależny współczynnik jedynie od kąta 2α (dla gniazda bez fazki). Taka po-stać wzoru uzasadniona jest analogią do opisu matematycznego przepływu laminarne-go w przewodach hydraulicznych.
W obszarze liczb Reynoldsa większych od krytycznej, współczynnik przepływu był wartością stałą, zależną również jedynie od kąta 2α (dla gniazda bez fazki).
W obszarze, w którym współczynnik przepływu można określić funkcją (7.41), w bardzo łatwy sposób uzyskuje się analityczny wzór łączący natężenie przepływu z różnicą ciśnień oraz przemieszczeniem grzybka. Aby to zrobić, należy dokonać od-powiednich przekształceń z użyciem wzorów (7.38), (7.39), (7.40) oraz (7.41). Zależ-ność ta wygląda następująco:
2 2 2 2 2 π sin ( cos sin ) a pz Q D z α ρ α α υ ∆ = − (7.42)
Nieco inne podejście do problemu jest przedstawione w pracy [7.6]. Zmienność współczynnika przepływu można określić następującą funkcją wykładniczą:
(Re) c(1 exp( b Re)
µ
=µ
− − (7.43)gdzie:
µ
c – stały współczynnik przepływu, występujący przy wartościach liczby Rey-noldsa dużo większe od krytycznej, b – parametr. Obydwa parametry są zależne od kąta 2α (z ostrokrawędziową formą gniazda).Zależność (7.41) opisuje zarówno zmienność współczynnika przepływu przy ni-skich wartościach, jak i jego stałość dla dużych wartości liczby Reynoldsa.
Niestety, analityczny zapis natężenia przepływu jako funkcji ciśnienia oraz prze-mieszczenia grzybka w sytuacji, gdy korzystamy z równania (7.43), jest niemożliwy ze względu na nieprzekształcalność wzorów (7.38), (7.39), (7.40) oraz (7.43).
7.3.2. ROZWIĄZANIE NUMERYCZNE MODELU MATEMATYCZNEGO
Rozwiązanie modelu matematycznego zostało zrealizowane dla zaworu dławiące-go o następujących parametrach: D = 2 mm (gniazdo niesfazowanane), 2α = 60°, v = 64·10–6 m2/s oraz ρ = 880 kg/m3. Powyższe własności odpowiadają własnościom oleju L-HL68 w temperaturze 42 °C. Założono ponadto, że zawór pracował będzie przy natężeniach przepływu 0–1,4 l/min i pod ciśnieniem 0–20 MPa.
W takim zminiaturyzowanym układzie grzybek–gniazdo przyjęcie stałego współ-czynnika przepływu jest zbyt dużym uproszczeniem. Biorąc po uwagę badania zaworu
modelowego (por. rozdz. 5), dla tak sprecyzowanego układu grzybek–gniazdo współ-czynnik przepływu przyjmuje wartość stałą dopiero dla liczb Reynoldsa większych od 300. Tymczasem maksymalna liczba Reynoldsa obliczona dla natężenia przepływu 1,4 l/min wynosi 118. Dlatego też badany układ grzybek–gniazdo będzie pracował przy przepływie laminarnym oraz przejściowym, gdy współczynnik przepływu będzie zmieniał się w całym zakresie przewidywanych wartości liczby Reynoldsa. Współ-czynnik przepływu należy więc aproksymować funkcją F1 (7.41) lub funkcją F2 (7.43).
Z badań zaworu modelowego wynika, że dla podanych parametrów zaworu dła-wiącego, dla funkcji F1 a = 0,063. Z pozycji [7.6] wynika natomiast, że dla funkcji F2 b = 0,1258, a
µ
c= 0,82.Mając określone parametry zaworu oraz równań (7.41) i (7.43), można przystąpić do rozwiązania modelu matematycznego zaworu. Rozwiązanie to zostało przedsta-wione na rysunkach 7.9, 7.10 oraz 7.11 w postaci:
• przykładowych charakterystyk oporów przepływu Q = f(∆p) dla dwóch różnych przemieszczeń grzybka z, gdzie współczynnik przepływu aproksymowany był funkcją F1 oraz F2,
• przykładowych charakterystyk sterowania Q = f(z) dla dwóch różnic ciśnień ∆p, gdzie współczynnik przepływu aproksymowany był funkcją F1 oraz F2,
• zbiorczego wykresu trójwymiarowego, określającego zależność Q = f(z,∆p), gdzie współczynnik przepływu aproksymowany był funkcją F1 (jest to ilustracja funkcji dwóch zmiennych oznaczona wzorem (7.42)).
Rys. 7.9. Charakterystyki oporów przepływu zaworu dławiącego aproksymowane dwiema funkcjami dla przemieszczeń grzybka
Z rysunków 7.10 oraz 7.11 wynika, że aproksymacja współczynnika przepływu funkcjami F1 oraz F2 daje zbliżone wartości w sytuacji, gdy natężenie przepływu wynosi około 0,8 l/min. Poniżej lub powyżej tej wartości występuje pewna rozbież-ność między funkcjami. W podrozdziale 10.4 wyniki badań doświadczalnych zostały porównane z wyżej opisanymi aproksymacjami.
Rys. 7.10. Charakterystyki sterowania zaworu dławiącego aproksymowane dwiema funkcjami dla różnic ciśnień
∆p = 2 MPa oraz ∆p = 5 MPa
Rys. 7.11. Ogólna charakterystyka przepływowa zaworu dławiącego zapisana jako funkcja dwóch zmiennych (zgodnie z wzorem (7.42))
7.4. OPRACOWANIE MODELU MATEMATYCZNEGO