Category Archives: Pisanie prac

Pisanie prac inżynierskich

Prawo Pascala, wykorzystywanie płynów w układzie zamkniętym

Hydraulika  –  jest  to nauka badająca  prawa  stanu  równowagi  i  ruchu cieczy.  Zagadnienia dotyczą między innymi urządzeń i budowli hydrotechniczych maszyn hydraulicznych i mechanizmów.

Żyjący w latach 1623-1662 – francuski matematyk, fizyk, pisarz i filozof Blaise Pascal badał między innymi zjawiska dotyczące ciśnienia atmosferycznego  i hydrostatyki.

W roku 1651 sformułował prawo przyjęte jako prawo Pascala mówiące, że: „Ciśnienie wywarte w jednym punkcie cieczy znajdującej się w równowadze (mechanicznej i cieplnej) rozchodzi się jednakowo we wszystkich kierunkach”.

Ciśnienie cieczy i gazu mierzymy wartością siły z jaką ciecz naciska na jednostkę powierzchni stykającego się z nią ciała.

Rys. 1   Ciecz w naczyniu zamkniętym poddana naciskowi tłoka.

rys1

Na pewną masę cieczy wywiera się za pomocą tłoka nacisk siłą F. Jeżeli pole przekroju poprzecznego tłoka wynosi A, to ciśnienie P jakie wywierane jest na jednostkę powierzchni wynosi       P = F / A

Gdzie:

F – siła w N

A – pole przekroju w m2

P – ciśnienie cieczy w Pa (paskalach)

Jednostką ciśnienia i naprężenia mechanicznego w układzie SI jest pascal – Pa

1 Pa = 1 N / m2

Współczynniki zamiany jednostek

1 N / m2 = 1 Pa

1 bar = 105 Pa = 0,1 Mpa

1 kG / m2 = 9,8066 Pa

1 kG / cm2 = 98,066 Pa

1 atm = 1 kG / cm2 = 98,066 Pa

Napęd hydrauliczny  – jest to napęd wywołany ruchem cieczy pod ciśnieniem, oparty na prawie Pascala, natomiast sterowanie hydrauliczne jest procesem kierowania pracą maszyny lub urządzenia, za pomocą cieczy pod ciśnieniem.

Każdy rodzaj napędu posiada swe zalety i wady. Oto niektóre z nich.

Więcej informacji na temat tej pracy znajdziecie na stronie:

Hydrauliczny sprzęt ratowniczy na wyposażeniu Straży Pożarnej

Ręczne i nożne pompy hydrauliczne

Występują w postaci jedno, dwu i trzystopniowej o mocnej zwartej konstrukcji i dużej sprawności. Pompy posiadają automatyczne zawory przełączające na poszczególne stopnie sprężania, zawory spustowo-ciśnieniowe ustawione na wymagane maksymalne ciśnienie i automatyczne odpowietrzenie.

Praca ich odbywać się może w pomieszczeniach zamkniętych i strefach zagrożonych wybuchem.

Rys. 11  Pompy mechaniczne, od lewej – dwustopniowa pompa nożna, trójstopniowa pompa ręczna.

r11Źródło: Archiwum autora

Serwer mediów – definicja, funkcjonalność i architektura

Niniejszy rozdział zawiera definicję bloku funkcjonalnego serwera mediów (Media Server, MS). Przedstawiono zakres jego funkcjonalności, opisano architekturę funkcjonalną oraz interfejsy do innych elementów sieciowych. Zaproponowano wprowadzenie dwóch typów serwera w zależności od roli pełnionej w sieci i zbioru wspieranych funkcji.

Definicja pojęcia serwera mediów

Przyjmijmy, że:

Serwer mediów (Media Server, MS) to element funkcjonalny pracujący w sieci pakietowej typu VOIP [1], który dokonuje operacji na strumieniach audio i video oraz integruje funkcjonalności niezbędne do realizacji różnego typu usług multimedialnych. Źródłem strumieni mogą być inne elementy sieciowe jak również lokalne zasoby serwera mediów (np. pliki multimedialne, syntezatory mowy etc).

Odnosząc się do przedstawionej w rozdziale 4.1 definicji, serwer mediów wspiera realizację usługi multimedialnej poprzez udostępnienie swojej funkcjonalności. Inne elementy sieciowe mogą przesyłać do MS żądanie dokonania pewnej operacji na strumieniu multimedialnym.

MS to bierny moduł, który nie zawiera żadnej logiki usługowej, nie ma więc możliwości samodzielnego świadczenia usług. Jego zasoby są dostępne dla wielu aplikacji usługowych jednocześnie, wykorzystywane zgodnie ze zdefiniowaną w aplikacjach logiką. Innymi słowy, MS to pewnego rodzaju zaawansowany „procesor sygnałowy” w sieci VOIP.

Umiejscowienie w danej sieci VOIP modułu MS pozwala na zlokalizowanie wszelkich operacji na sygnałach w jednym elemencie, do którego mogą odwoływać się istniejące i nowo wprowadzane aplikacje usługowe. Eliminuje to konieczność dodawania osobnych, dedykowanych modułów.


[1]   Należy podkreślić, iż definicja nie określa styku pomiędzy MS a sieciami z komutacją kanałów. Ten styk jest realizowany poprzez bramy medialne (Media Gateways), które są modułami o innej funkcjonalności niż MS. MS pracuje tylko w sieci z komutacją pakietów.

Rozszerzenia protokołu w kontekście usług multimedialnych

Protokół SIP dzięki prostemu mechanizmowi działania oraz przejrzystości składni jest szeroko wykorzystywany w ramach rozmaitych usług telekomunikacyjnych i internetowych. Nie pozwala jednak na realizację wyrafinowanych usług multimedialnych.

W odpowiedzi na ograniczenia protokołu, w ramach IETF podjęto prace, których rezultatem są liczne dokumenty opisujące rozszerzenia „podstawowej” specyfikacji SIP zawartej w [1]. Ponadto nieustannie opisywane są nowe możliwości zastosowania protokołu SIP.

W kontekście zastosowań multimedialnych, spośród opublikowanych draftów i norm RFC można wyróżnić m.in.:

  • [6] – definicja wiadomości SIP REFER
  • [7] – definicja wiadomości SIP INFO
  • [14] – definiuje architekturę usługi konferencji opartej na SIP
  • [24] – specyfikuje mechanizm reprezentowania pewnych usług multimedialnych w postaci adresów agentów SIP
  • [33] – opisuje sposób sterowania rozproszonym systemem konferencyjnym z wykorzystaniem protokołu SIP
  • [34] – specyfikuje mechanizm wskazywania skryptów VXML w wiadomości SIP przesyłanej do elementu sieciowego zdolnego do interpretacji tych skryptów
  • [36] – definicja zbioru funkcjonalności podstawowej interakcji z użytkownikiem w oparciu o protokół
  • [37] – definicja zbioru podstawowych funkcjonalności konferencyjnych

Wykorzystanie protokołu SIP do realizacji różnych usług multimedialnych to temat aktualnie szeroko rozważany w ramach organizacji standaryzacyjnych. Niniejsza praca dyplomowa pokazuje, jak wiele kwestii w tej dziedzinie wymaga jeszcze doprecyzowania.

[1] IETF RFC 3261 “SIP: Session Initiation Protocol”, Czerwiec 2002

Podsumowanie pracy magisterskiej

Celem pracy było zobrazowanie potencjału tkwiącego w systemach CMS na podstawie różnych przykładów oraz zastosowanie Joomla! do stworzenia strony nadleśnictwa.

Przedstawione serwisy internetowe Wydziału Inżynierii Środowiska, Instytutu Badawczego Leśnictwa obrazują wykorzystanie CMS w dziedzinie edukacji. Z kolei projekt Walter przedstawia ideę wykorzystania CMS jako narzędzia, dzięki któremu możemy gromadzić wiedzę oraz materiały pochodzące z różnych instytucji a nawet z różnych krajów w jednym miejscu. Używanie przez wiele osób jednego narzędzia jest sporym plusem, gdyż nie zachodzi potrzeba dodatkowego przenoszenia i integrowania danych. Strony Krakowskiego Zarządu Komunalnego i Regionalnego Zarządu Gospodarki Wodnej pokazują przystosowanie systemu dla dwóch różnych instytucji, których strony pełnią rolę informacyjną, administracyjną i komunikacyjną. Strona Woda i Powietrze przedstawia komercyjny charakter wykorzystania CMS-u, jej twórcy skupili się na reklamie firm związanych z inżynierią środowiska. Strona organizacji Polskich Łowców Burz obrazuje, iż w sposób amatorski można stworzyć portal o ciekawej tematyce, który znajduje sporą liczbę odbiorców.

Pierwszym etapem mojej pracy było przygotowanie środowiska do instalacji Joomla! oraz zainstalowanie i skonfigurowanie programu. Następnie stosując ogólnodostępne moduły stworzyłem serwis WWW dla nadleśnictwa o przeznaczeniu zarówno edukacyjnym, promującym, informującym jak i komunikacyjnym. Zastosowałem komponenty, dzięki którym można umieszczać w prosty i szybki sposób aktualne informacje, przedstawiać oferty handlowe, mieć dostęp do struktury organizacyjnej pracowników nadleśnictwa, korzystać z Biuletynu Informacji Publicznej itd.

Serwis spełnił założone przeze mnie wymagania ustalone podczas rozmów z pracownikami nadleśnictwa, tzn. jest prosty w obsłudze nawet dla użytkowników mających słabe doświadczenie z komputerem. Koszt wykonania witryny www jest niski i nie wymaga zastosowania specjalistycznego oprogramowania. Informacje zawarte w Joomla! są poprawnie wyświetlane, a użycie dodatkowych modułów usprawnia pracę osobom zatrudnionym w nadleśnictwie.

Wydział Inżynierii Środowiska Politechniki Krakowskiej

Studentom Politechniki Krakowskiej jest dobrze znany portal Wydziału Inżynierii Środowiska (rys.6). Studenci na nim mogą pobrać zawsze aktualny plan zajęć, sprawdzić swoją średnią ocen (e-student), znaleźć materiały edukacyjne  albo zapisać się na moduły do wyboru. Dodatkowo znajdziemy w portalu całą administrację i wykaz pracowników Wydziału, także aktualności, np. bieżące informacje o wydziałowych wydarzeniach czy listę potrzebnych dokumentów na różne okazje, problemy, wnioski.

Rysunek 6. Strona internetowa Wydziału Inżynierii Środowiska Politechniki Krakowskiej

1

Źródło: wis.pk.edu.pl

Strona WIŚ stworzona CMS świetnie spisuje się w zakresie organizacyjn0-edukacyjnym. Jest przyjazna dla studentów, którzy w czytelny sposób mogą znaleźć na niej bieżące informacje. To jest duże ułatwienie biorąc pod uwagę, iż studiuje tutaj bardzo dużo osób. Sądzę, że ciekawym rozwiązaniem będzie umieszczenie w portalu  wyników z egzaminów – co w najbliższym czasie nastąpi…

Zbiór stanów użytkowania

Zbiór stanów użytkowania, które mogą wystąpić przy reali­zacji dowolnego zadania  Za  w dowolnych warunkach Ot, oznaczmy przez  s (s1,  s2, …,sw),  gdzie  w  jest liczbą moż­liwych stanów. Ze względu na oddziaływania zewnętrzne (głównie obciążenia dynamiczne) ważne są nie tylko liczby poszczegól­nych stanów, jakie mogą się pojawić przy wykonywaniu zadania, ale i rodzaje poprzedzających je stanów. Takie informacje za­wiera macierz prawdopodobieństw zmian stanów użytkowania. Praw­dopodobieństwo przejścia ze stanu j do stanu k oznaczmy przez qjk. W wielu przypadkach można założyć, że qjk nie zale­żą od czasu i są jednakowe dla każdego z zadań. Zbiorowi sta­nów  s  można przyporządkować następującą macierz przejść

q =                                                                              (13)

Wyrazy tej macierzy spełniają zależność

= 1            dla każdego  j.                                                                     (14)

Ważna ze względu na oddziaływania zewnętrzne liczba przejść njk ze stanu  j  do stanu  k,  przypadająca na jedno zadanie jest równa w przybliżeniu

njk = qjknj .                                                                                                        (15)

gdzie  nj  jest liczbą zdarzeń, polegających na pojawieniu się stanu  j  podczas wykonywania zadania. W pewnych przypadkach można przyjąć, że  nj nie zależą od czasu  i  od rodzaju za­dania. Z równości (15) wynika, że do wyznaczenia liczb njk potrzebna jest znajomość wyrazów macierzy  q  oraz wyrazów nj macierzy wierszowej  n.  Można je określić w badaniach eks­ploatacyjnych obiektów podobnych.

Ważną informacją o sposobie użytkowania niektórych obiek­tów może być znajomość czasów  tsj  (j = 1, 2,…, w)  trwania poszczególnych stanów użytkowania. W ogólnym przypadku można je opisać na przykład za pomocą zbioru dystrybuant

Fs1(ts1), Fs2(ts2),…, Fsw(tsw) .                                                                           (16)

Po zakończeniu zadania obiekt może przejść do stanu ocze­kiwania na nowe zadanie. Stan ten występuje na ogół tylko wów­czas, gdy zgłoszenie zadania wystąpiło dopiero po zakończeniu zadania poprzedniego. Można założyć, że proces zgłoszeń zadań ma cechy braku następstw, pojedynczości i stacjonarności. Wówczas może być on traktowany jako proces Poissona [5], a czas między zgłoszeniami opisuje zmienna losowa o rozkła­dzie wykładniczym. Gęstość prawdopodobieństwa tej zmiennej wy­raża się wzorem

fx(t) =                                                         (17)

gdzie  lz   jest intensywnością zgłoszeń.

Inną grupą elementów sposobu użytkowania są elementy zwią­zane z jakością użytkowania. Duża część z nich to błędy użyt­kowania. Jeśli przypisać im wielkości B1 ,B2,…  itd.  (na przykład przy skośnym ciągnięciu ładunku przez żuraw, istotną wielkością opisującą ten błąd, może być kątowe odchylenie liny od pionu), to w ogólnym przypadku opisem tych błędów może być zbiór ich dystrybuant lub gęstości prawdopodobieństwa

FB1(B1), FB2(B2),…, FBb(Bb) .                                                                        (18)

gdzie  b  jest liczbą błędów uwzględnionych w modelu.

Wpływ jakości użytkowania na te błędy może być uwzględnio­ny w wartościach parametrów zmiennych B1 ,B2,…, Bb.

Od jakości użytkowania zależą też inne elementy sposobu użytkowania, np. liczba sterowanych ruchów nieobligatoryjnych przy wykonywaniu zadania przez maszynę. Matematyczny opis tej zależności może mieć na przykład następującą postać

n(IO) = k(IO)  n(1) .                                                                                        (19)

gdzie  n(IO)   jest zdefiniowaną już poprzednio macierzą n, określoną dla klasy operatora o numerze

IO  , n( l )  jest taką samą macierzą określoną dla najwyższej jakości użytkowania, a  k(IO)  jest odpowiednio dobranym współczynnikiem. Czym wyż­sza jakość użytkowania, tym wartość współczynnika  k(IO)  jest mniejsza.

Wszystkie przedstawione powyżej elementy sposobu eksploa­tacji U dotyczyły użytkowania obiektu. Przedstawmy teraz dru­gą grupę elementów sposobu eksploatacji, charakteryzujących proces obsługiwania. Wartości wielu z nich narzucone są zwykle przez instrukcje obsługiwania. Chodzi tu głównie o okresy mię­dzy różnego rodzaju przeglądami i odnowami profilaktycznymi.

Użytkowanie obiektu może też być przerywane w sposób nie­zaplanowany przez pojawiające się niesprawności tego obiektu i Innych zespołów maszyny, której częścią jest badany obiekt. Tworzony model niezawodnościowy obiektu pozwala na wyznacze­nie chwil występowania ewentualnych jego niesprawności, nato­miast innych zespołów maszyny nie. Bardzo trudne jest też uzy­skanie wiarygodnych informacji o przyszłym odnawianiu poawaryjnym, zwłaszcza innych zespołów maszyny, np. informacji o czasach oczekiwania na naprawę. W takich przypadkach propo­nuje się pominąć czasy trwania odnów poawaryjnych. Oznacza to zubożenie niezawodnościowego modelu obiektu l rezygnację z ba­dania m.in. wpływu organizacji procesu obsługiwania na poziom niezawodności. Jednak przy zasadniczych celach, stawianych przed badaniami niezawodnościowymi w ramach proponowanego sy­stemu racjonalnego oddziaływania na poziom niezawodności (rozdz. 2), wprowadzenie takiego uproszczenia na ogół nie zmniejsza znacząco wartości modelu.

Łatwiej jest natomiast w tworzonym modelu eksploatacji obiektu uwzględnić efekty odnawiania poawaryjnego, podobnie jak efekty planowych przeglądów i odnów profilaktycznych. Efekty te wynikają z jakości tych zabiegów obsługiwania. Przy opisie matematycznym jakości obsługiwania można na przykład dokonać podziału jakości obsługiwania na kilka klas.  Każdej z klas przypisać można odpowiednie wartości wielkości określa­jących jakość zabiegów obsługiwania. Czym wyższa na przykład jakość przeprowadzania przeglądów, tym większe wartości mogą mieć prawdopodobieństwa

Pwyk1(Iprz), Pwyk2(Iprz),…, Pwykm(Iprz)                                                              (20)

wykrycia nieodpowiedniego stanu technicznego każdego z PK podlegających przeglądom. Symbolem  Iprz Pwyk1(Iprz), oznaczona została klasa jakości przeglądów, a symbolem m – liczba PK. Dla PK, które nie podlegają przeglądom, prawdopodobieństwa (31) nie są określane lub na podstawie umowy są przyjmowane jako równe zeru.

Jakość odnowy (profilaktycznej i poawaryjnej) jest bardzo dobra, gdy w wyniku tej odnowy przywraca się początkowy stan techniczny PK, a zła, gdy nie poprawia ona aktualnego nieod­powiedniego stanu technicznego. Jednak w ogólnym przypadku stan techniczny PK po odnowie może przyjmować dowolny poziom, będący rezultatem m.in. oddziaływań różnych czynników losowych na zabieg odnowy. Wówczas o jakości odnowy  PK może świadczyć na przykład średnie odchylenie kwadratowe  sod wielkości określającej stan techniczny po odnowie eOd od wartości oczekiwanej tej wielkości w chwili  t = 0  (na początku eksploatacji), czyli od EeO  Odchylenie sod jest określone za pomocą wzoru

sod (Iod) =  .                                                                 (21)

gdzie  Iod  jest symbolem klasy jakości odnowy. Czym wyższa jakość przeprowadzania odnowy, tym mniejsze są wartości  sod. Łatwo uzasadnić, że

sod (Iod) =  .                                            (22)

Jak widać, miara jakości odnowy  sod uwzględnia zarówno przesunięcie całego rozkładu (różnica systematyczna), jak i .różnice rozproszeń losowych wielkości określającej stan tech­niczny po odnowie eOd i w chwili  t = 0, czyli eO w  szcze­gólnym przypadku, gdy stan techniczny PK na początku eksploa­tacji i po panowie jest opisywany przez wielkości zdetermino­wane, odchylenie  sod jest równe

sod (Iod) = eOd(Iod) – eO                                                                                    (23)

Do odwzorowania stanu technicznego PK po odnowie informa­cja o  sod (Iod)   jest jednak zbyt skąpa. Konieczna jest do te­go dodatkowo znajomość na przykład postaci matematycznej roz­kładu zmiennej losowej  eOd oraz  EeOd, które na podstawie umowy można w wielu przypadkach przyjąć jako niezależne od ja­kości odnowy.

Przedstawione powyżej zasady tworzenia modelu procesu eks­ploatacji mogą być wykorzystane do modelowania dowolnego obiektu mechanicznego.

c) Model oddziaływań zewnętrznych na obiekt jest to zbiór informacji lub relacji matematycznych określających zbiór od­działywań roboczych i oddziaływań otoczenia (rozdz. 3). Naj­bardziej istotnymi oddziaływaniami zewnętrznymi są zazwyczaj obciążenia zewnętrzne. W przypadku, gdy rozpatrywanym obiektem jest maszyna, na obciążenia ta składają się z jednej strony obciążenia czynne (od silników), np. moment obrotowy na wałku silnika, a z drugiej strony – obciążenia bierne, wynikające przede wszystkim z charakteru pracy użytecznej wykonywanej przez tę maszynę w ramach danego stanu eksploatacji (p. punkt b). W przypadkach innych obiektów ich zewnętrzne obciążenia wynikają ze współpracy z sąsiednimi zespołami lub elementami maszyny.

Dla chwil odpowiednio odległych od początku trwania danego stanu eksploatacji obciążenia zewnętrzne są zwykle stałe lub stacjonarne. Prócz tego w ramach tego stanu występują obciąże­nia niestacjonarne, związane ze zmianą warunków pracy maszyny. Chodzi tu głównie o obciążenia dynamiczne pojawiające się przy zmianie stanów użytkowania. Zależą one od rodzajów stanów po­przedniego i rozpatrywanego. Obciążenia zewnętrzne obiektu po­jawiające się w stanie użytkowania  k  pod warunkiem, że sta­nem poprzednim był stan  j,  oznaczmy więc symbolem  Ok/j(t), gdzie  t może być tu traktowane jako czas użytkowania, liczo­ny na przykład dla każdego stanu  k  oddzielnie. Na ogół można przyjąć, że dla stanu innego niż stan użytkowania obciążenia zewnętrzne są równe zeru.

Jeśli wszystkich stanów użytkowania jest  w,  to zbiór ob­ciążeń Ok/j(t)  można przedstawić w sposób poglądowy za pomo­cą macierzy

º O(t)                                                     (24)

Ponieważ nie wszystkie przejścia ze stanu do stanu są moż­liwe, niektóre z wyrazów tej macierzy można traktować na pod­stawie umowy za równe zeru. Macierz (24) związana jest ściśle z macierzą (13) określającą prawdopodobieństwa zmiany stanów.

W pewnych przypadkach prócz oddziaływań w postaci obciążeń mogą być istotne i inne oddziaływania zewnętrzne, z powodu du­żego ich wpływu na oddziaływania wewnętrzne i stan wytężeń w poszczególnych PK. Tymi innymi oddziaływaniami mogą być od­działywania określane na przykład za pomocą takich wielkości jak: temperatura otoczenia, wilgotność powietrza, zapylenie (czasami zapiaszczenie) powietrza itd. Ich opis zawarty jest w modelu eksploatacji obiektu (punkt b).

W ogólnym przypadku oddziaływania zewnętrzne dla zbioru rozpatrywanych obiektów, a więc także wyrazy macierzy O(t), można traktować jako procesy stochastyczne. W praktyce matema­tyczny model oddziaływań zewnętrznych może być przedstawiony w różny sposób, na przykład w postaci zbioru losowych lub zde­terminowanych funkcji czasu, zbioru histogramów oddziaływań (m.in. obciążeń) itd. W niektórych przypadkach jego zasadnicza część, dotycząca obciążeń, może być rezultatem pewnych rozwa­żań teoretycznych, wykorzystujących prawa statyki i dynamiki. Zwykle konieczne jest wówczas zbudowanie modelu dynamicznego (p. punkt d) rozpatrywanego obiektu lub maszyny, której jest on częścią.

d) Model oddziaływań wewnętrznych obiektu jest to zbiór relacji matematycznych pozwalających na wyznaczenie tych od­działywań pod wpływem oddziaływań zewnętrznych. W przypadku obiektów mechanicznych jest to głównie zbiór relacji, który pozwala na wyznaczenie obciążeń elementów (i PK) obiektu pod wpływem obciążeń zewnętrznych w warunkach innych oddziaływań zewnętrznych (np. temperatury), czyli na wyznaczenie wyrazów Ok/j,i(t) macierzy Oi(t) obciążeń elementu z i-tym PK, gdzie  j, k = 1, 2,…, w,  a  i = 1, 2,…, m,  spowodowanych ob­ciążeniami zewnętrznymi opisanymi przez wyrazy Ok/j,i(t)   macie­rzy  O(t)  (24).  Zasadniczą częścią modelu oddziaływań we­wnętrznych jest więc model dynamiczny badanego obiektu. Pozo­stałe relacje wchodzące w skład modelu oddziaływań wewnętrz­nych to na przykład wyrażenia pozwalające określić obciążenia Oi(t)  na podstawie rozwiązań modelu dynamicznego. Gdy rozpa­trywany obiekt jest częścią maszyny, to model dynamiczny tej maszyny może być wykorzystywany również do określenia obciążeń zewnętrznych tego obiektu (p. punkt o). Tworzenie modelu ob­ciążeń wewnętrznych polega więc w ogólnym przypadku na tworze­niu przede wszystkim odpowiedniego modelu dynamicznego, które­go zasadniczą część stanowi zbiór równań na przykład o nastę­pującej postaci

(25)

wynikającej z przyjmowanych zwykle założeń upraszczających. Równania takie tworzy się dla tych stanów użytkowania  j  oraz k,  dla których wyrazy macierzy (13) zmiany stanów są niezerowe. Symbolami  Mk/j , Ck/j  i  Kk/j  oznaczone są w wyrażeniu (25) odpowiednio macierze bezwładności, tłumienia i sztywno­ści, a symbole xk/j , i Ok/j (t)  użyte są do oznaczenia ko­lumnowych macierzy uogólnionych przemieszczeń i obciążeń ze­wnętrznych. Podczas tworzenia takiej postaci modelu dynamicz­nego określa się również wspomniane macierze oraz warunki po­czątkowe.

Własności elementów i ich połączeń w obiektach mechanicz­nych (opisywane przez parametry Mk/j , Ck/j , Kk/j ) mają często duże rozrzuty losowe. Wówczas parametry modelu można trak­tować jako zmienne losowe, a ze względu, na zjawiska zużycia i zmienność oddziaływań zewnętrznych (głównie otoczenia) – nawet jako procesy losowe. Również i obciążenia zewnętrzne mogą być w ogólnym przypadku traktowane jako procesy losowe. Jednakże modele takie, modele stochastyczne, nawet w przypadku nie­skomplikowanych maszyn, są na tyle złożone, że często nie na­dają się do analitycznego rozwiązania. Zwykle korzysta się wtedy z dalszych uproszczeń i buduje model o charakterze de­terministycznym, pozostając przy stochastycznym wymuszeniu, choć i takie modele nie są łatwe do analizy.

W przypadku konkretnego obiektu matematyczna postać modelu dynamicznego zależy od przyjętego w niezawodnościowym modelu nominalnym, a także matematycznym, opisu procesu eksploatacji lub opisu rodzajów i poziomów zewnętrznych obciążeń mechanicz­nych, od struktury funkcjonalnej obiektu itd. W dużym stopniu postać modelu dynamicznego zależy również od celu, jakiemu ma on służyć. W przypadku badań niezawodności ważne jest wyzna­czenie, na przykład na podstawie przebiegów obciążeń elemen­tów  Oi(t)  uzyskanych w wyniku analizy modelu dynamicznego, nie tylko ekstremalnych obciążeń elementów, lecz także innych charakterystyk przebiegu obciążeń dynamicznych istotnych z punktu widzenia niezawodności, np. szybkości tłumienia drgań, rozkładu szczytowych obciążeń itd.

Wyznaczanie obciążeń wewnętrznych jest znacznie łatwiej­sze, jeśli w modelu obciążeń zewnętrznych obiektu przyjęto, że obciążenia te są quasi-statyczne. Wówczas często dobrym mode­lem obciążeń elementów jest zbiór funkcyjnych związków linio­wych między obciążeniami uogólnionymi elementów Ok/j, i(t) º Ok, i(t) a obciążeniami uogólnionymi zewnętrznymi Ok/j(t) º Ok(t)  typu

Ok, i(t) = bk, i Ok(t)                                                                                            (26)

gdzie  bk, i  jest stałym współczynnikiem.

Jak już napisano, modelem oddziaływań wewnętrznych są głów­nie relacje matematyczne służące do wyznaczenia obciążeń ele­mentów pod wpływem obciążeń zewnętrznych. W pewnych jednak przypadkach do modelu tego mogą wchodzić relacje pozwalające na wyznaczenie obciążeń elementów powstających pod wpływem in­nych oddziaływań zewnętrznych niż obciążenia (np. zmian tempe­ratury) lub relacje pozwalające wyznaczyć Inne oddziaływania wewnętrzne niż obciążenia (np. termodynamiczne oraz odkształ­cenia lub wzajemne przemieszczenia elementów).

e) Model uogólnionych wytężeń w PK jest to zbiór relacji matematycznych służący do wyznaczenia uogólnionych wytężeń w PK, powstałych pod wpływem oddziaływań, którym poddany jest element obiegu. W przypadku obiektów mechanicznych istotne są zwykle wytężenia mechaniczne, np. naprężenia lub naciski. Są jednak i takie obiekty mechaniczne, dla których istotne są również wytężenia innego rodzaju, np. termodynamiczne. Wspom­niany zbiór relacji służy więc przede wszystkim do wyznaczenia wytężeń (mechanicznych i innych) bk/j, i(t) wywołanych w i-tym PK przez zewnętrzne oddziaływania, głównie obciążenia Ok/j, i(t). Element z tym PK, określono w modelu oddziały­wań wewnętrznych (p. punkt d), dla wszystkich wyszczególnio­nych w tym modelu wartości i, j  oraz

Jak wiadomo, nawet przy zdeterminowanym charakterze od­działywań na element, wytężenia w PK tego elementu mogą mieć duże rozproszenia losowe. Na przykład w miejscach spiętrzenia naprężeń bardzo duży wpływ na ich maksymalne wartości mają wy­miary tych fragmentów elementu, które wywołują spiętrzenia (np. promienie zaokrągleń odsadzeń na wałkach, promienie za­okrągleń przy dnach rowków wpustowych, wymiary spoin, itd.). Ze względu na to, że wymiary te mają często stosunkowo duże rozproszenia losowe, gradienty naprężeń i maksymalne napręże­nia w miejscach spiętrzeń też mogą mieć duże rozproszenia lo­sowe.

Jeśli ten fakt uwzględniono już w modelu początkowego sta­nu technicznego obiektu, to model wytężeń mechanicznych może być modelem deterministycznym. W tym przypadku służy on do wy­znaczenia nominalnych wytężeń (bez uwzględnienia karbów) w PK, i wpływ wymiarów, kształtu itd. na spiętrzenie naprężeń uwzględnia się w modelu początkowego stanu technicznego, a objawia się on zmniejszeniem niszczących naprężeń nominalnych. Jeśli tego nie uwzględniono w modelu początkowego stanu tech­nicznego, powinno to być uwzględnione w modelu wytężeń mechanicznych i wówczas ten model ma charakter stochastyczny (pro­babilistyczny), a wytężenia bk/j, i(t) są procesami lub zmiennymi losowymi.

Przy budowie modeli uogólnionych wytężeń korzysta się z różnych teoretycznych metod analitycznych, z teoretycznych metod numerycznych (np. MES), z metod eksperymentalnych (ela-stooptyka, tensometria i inne).

Postać modelu uogólnionych wytężeń w każdym PK jest cha­rakterystyczna dla tego PK i elementu, w którym on występuje orał dla rodzaju i charakteru oddziaływań, którym poddany jest ten element. Zależy również od postaci matematycznego opisu cech zdatności (p. punkt f). Na przykład do opisu nie­których cech zdatności nie jest potrzebna znajomość przebiegu wytężeń bk/j, i(t)   w czasie eksploatacji, a jedynie pewne wielkości charakterystyczne dla tego przebiegu, takie jak: na­prężenia maksymalne, naprężenia maksymalne i średnie (lub od­powiednie rozkłady prawdopodobieństw tych wielkości ) itd.

Określone z pomocą omawianego modelu uogólnione wytężenia w PK mogą być podstawą do oceny stopnia pogorszenia istotnych własności elementów w tych PK w rozpatrywanym okresie eksploa­tacji i do oceny bieżących wartości cech zdatności (np. zapa­su wytrzymałości, przez porównanie wyznaczonych naprężeń ł na­prężeniami krytycznymi). Do ocen takich konieczne jest jednak istnienie modeli przedstawionych w dalszej części tekstu.

f) Model zmian stanu technicznego obiektu stanowi matema­tyczny opis tych zmian, głównie procesów degradacji stanu technicznego obiektu pod wpływem oddziaływań wewnętrznych (zwłaszcza obciążeń), a bezpośrednio – wskutek między innymi działania uogólnionych wytężeń (np. zmiennych naprężeń). W przypadku obiektu mechanicznego są to zwykle procesy zmęcze­nia objętościowego i powierzchniowego, zużycia ciernego, peł­zania itd. W zależności od celu badań używane modele tych pro­cesów są deterministyczne lub stochastyczne. Z punktu widzenia niezawodności szczególnie interesujące są modele stochastycz­ne. Najwygodniej jest odwzorować zmiany stanu technicznego nie za pomocą wiel­kości charakteryzujących ten stan (p. punkt a), lecz za pomocą wielkości przyjętych wcześniej cech zdatności.

Przedstawmy dla przykładu ogólną postać jednego z modeli zmiany stanu technicznego PK obiektu, spowodowanej zjawiskiem zmęczenia materiału. Za miarę zmian stanu technicznego można w tym przypadku przyjąć na przykład wielkość względnego uszko­dzenia zmęczeniowego  D.  Wielkość ta może też spełniać rolę cechy zdatności. W ogólnym przypadku, gdy naprężenie s(t) w PK jest traktowane jako proces stochastyczny, wielkość  D (a właściwie jej wartość oczekiwana) jest określona zależno­ścią

,                                                                           (27)

gdzie:

f(ss)  jest gęstością prawdopodobieństwa szczytowych war­tości naprężenia s(t),

t – czasem użytkowania,

nO – efektywną liczbą oscylacji naprężenia w jednostce czasu,

N(s) – trwałością zmęczeniową dla poziomu naprężeń  s wyznaczoną na podstawie odpowiedniego równania krzywej zmęcze­niowej,

Wd – dolnym ograniczeniem wytrzymałości zmęczeniowej  W.

Jeśli naprężenie  s(t)  jest stacjonarnym procesem normal­nym o wąskim widmie i zerowej wartości oczekiwanej.

dla ss > 0                                                        (28)

(29)

gdzie  D2s i  D2 s są wariancjami procesów naprężenia i jego pochodnej względem czasu.

Równanie krzywej zmęczeniowej przyjmijmy w postaci

(30)

gdzie:

h  jest wykładnikiem potęgi decydującym o pochyleniu krzy­wej zmęczeniowej,

WO  jest charakterystyczną wytrzymałością zmęczeniową,  a

NO  jest trwałością zmęczeniową PK  odpowiadającą pozio­mowi naprężeń  ss = WO.

Dla stali konstrukcyjnych WO  jest na ogół trwałą wytrzy­małością zmęczeniową.

Podstawiając zależności (28) – (30) do związku (27), otrzymuje się wyrażenie

(31)

gdzie:

G(x)  jest funkcją gamma,  a G(x1, x2)  jest niepełną funkcją gamma.

Miara uszkodzenia zmęczeniowego  D(t)  (31) jest więc li­niową funkcją czasu.

Wytrzymałość zmęczeniowa WO  jest traktowana w wyrażeniu (31) jako wielkość zdeterminowana. Ponieważ jednak ta własność materiału (i konstrukcji) ma duże rozrzuty losowe, w wielu przypadkach przyjmuje się, że WO jest zmienną losową. Opisu­je się ją za pomocą rozkładu Weibulla lub logarytmo-normalnego, a zazwyczaj za pomocą rozkładu normalnego (p. punkt a). Wówczas opia funkcji losowej D(t),  a nawet wyznaczenie nie­których jej parametrów, jest już zadaniem znacznie trudniej­szym.

Należy zaznaczyć, że w piśmiennictwie można znaleźć rów­nież i inne modele zmian stanu technicznego PK wskutek zmęcze­nia. Jest to zresztą charakterystyczne dla modelowania w ogóle. Na przykład przedstawiony powyżej proces pogarszania się włas­ności wytrzymałościowych można też opisać za pomocą modelu Serensena opartego na hipotezie stopniowego obniżania się trwałej wytrzymałości zmęczeniowej wskutek działania tzw. cyk­li przeciążeniowych, za pomocą modeli opartych na teorii łańcuchów Markowa, za pomocą modeli proponowa­nych przez mechanikę pękania.

Modelem zmian stanu technicznego obiektu jest zbiór rela­cji matematycznych, takich jak na przykład relacja (31), opi­sujących pogarszanie się stanu technicznego we wszystkich PK tego obiektu.

g) Model granicy obszaru zdatności obiektu jest to matema­tyczny opis granic cech zdatności poszczególnych PK oraz mate­matyczny opis struktury niezawodnościowej tego obiektu, przy­jętej przy tworzeniu nominalnego modelu niezawodnościowego (p. krok czwarty w podrozdz.2).

Ustalenie granicznej wartości cechy zdatności PK  odbywa się na podstawie sformułowanej w modelu nominalnym definicji niesprawności togo PK i jest łatwe w przypadku, gdy przekro­czenie tej wartości przez cechę (na skutek zmian stanu tech­nicznego obiektu) uniemożliwia teoretycznie dalsze funkcjono­wanie PK w obiekcie, czyli jest równoznaczne z jego niesprawnością fizyczną. Chodzi tu o zdarzenia, które z założenia po­winny prowadzić do takich postaci naruszenia geometrii, jak:

pęknięcia, trwałe odkształcenia wskutek utraty stateczności itd. W tych przypadkach istnieje wyraźna granica między sta­nem zdatności i stanem niezdatności PK.  Jest nią na przy­kład krytyczna wartość zapasu uogólnionej wytrzymałości  (ze względu na utratę stateczności, ze względu na doraźne pęknię­cie itd.)  Zgr   = O  lub krytyczna wartość względnego uszkodze­nia zmęczeniowego  Dgr   = 1  (lub

Dgr   = a  według modelu Serensena).

W innych przypadkach, gdy zmiana makro- lub mikrogeometrii zachodzi w sposób stopniowy, to ustalenie tej granicy nie jest łatwe. Na ogół są do tego wykorzystywane odpowiednie informa­cje uzyskane z eksploatacji już istniejących obiektów podob­nych. Na przykład w przypadku wysięgnika żurawia samojezdnego taką informacją jest dopuszczalna maksymalna wartość luzu  między ślizgiem a współpracującym z nim członem wysięgnika, wynikająca z praktyki eksploatacji takich żurawi i narzucona w instrukcji obsługi.

Graniczna wartość cechy zdatności  PK,  zmieniającej się stopniowo, może być też określona na podstawie analizy wpływu zmian tej cechy na efektywność eksploatowania (lub dobroć) obiektu. Stopniowa zmiana cechy zdatności w czasie eksploata­cji obiektu powoduje między innymi wzrost ryzyka powstania uszkodzenia, wywołującego duże straty; pogorszenie poprawności funkcjonowania obiektu (np. spadek wydajności) itd. Są to nie­które z czynników określających efektywność eksploatowania obiektu. Po pewnym czasie eksploatacji obiektu cecha zdatności osiąga wartość, przy której efektywność ta zaczyna gwałtownie spadać albo obniża się do niedopuszczalnego poziomu. Taka war­tość cechy zdatności może być traktowana jako jej wartość gra­niczna. Do przeprowadzenia jej wyboru konieczne jest przyjęcie miary efektywności i dokonanie matematycznego opisu zależności tej miary od rozpatrywanej cechy zdatności. W niniejszej pracy zagadnieniem tym nie zajmowano się.

Chwila, w której co najmniej jedna z cech zdatności Zn(t) i-tego PK osiąga wartość graniczną  Zngr  jest traktowana ja­ko chwila powstania niesprawności tego PK (p. rozdz. 3), a czas  Ti ,  jaki do tej chwili upłynął od początku eksploatacji (lub od ostatniej odnowy poawaryjnej), jest czasem bezawaryj­nej pracy i-tego PK.

Druga część modelu granicy obszaru zdatności obiektu – to matematyczny opis struktury niezawodnościowej obiektu, przy­jętej w modelu nominalnym. Przy tworzeniu tego opisu należy podjąć decyzję o przyjęciu lub nieprzyjęciu założenia upra­szczającego polegającego na pominięciu zależności stocha­stycznych między zmiennymi losowymi   Ti (i  = 1, 2,….m) (spowodowanych zależnościami stochastycznymi między cechami zdatności tych PK). W przypadku obiektów mechanicznych, zwła­szcza przenoszących napęd, te zależności mogą być silne mię­dzy innymi z powodu istnienia silnych zależności między ob­ciążeniami poszczególnych PK. Aby podjąć wspomnianą decyzję, należy ocenić, jak silne są te zależności stochastyczne. W pewnych przypadkach oceny tej można również dokonać na podstawie analizy tych zależności między cechami zdatności obiektu. Duży wpływ na siłę tych zależności ma cha­rakter zjawisk fizycznych uwzględnionych przez model nominalny (p. rozdz. 5).

Najczęściej przyjmowaną strukturą niezawodnościową dla obiektów mechanicznych, zwłaszcza przenoszących napęd, jest. struktura szeregowa (p. krok czwarty w p. rozdz. 4.2). Matema­tyczny opis takiej struktury można ująć następującym związkiem między czasem  T bezawaryjnej pracy obiektu i czasami T- bez­awaryjnej pracy jego PK

T = min (T1, T2,…, Tm)                                                                                  (32)

oraz zbiorem E  związków, które można przedstawić w następują­cej ogólnej postaci

ye (T1, T2,…, Tm) = 0     e Î E                                                                        (33)

uwzględniających stochastyczne zależności między zmiennymi losowymi T.

Załóżmy, że T  i , Ti,   (i = 1, 2,…, m)  są zmiennymi loso­wymi ciągłymi określonymi dla  T ³ 0  i Ti ³ 0  o skończo­nych co najmniej pierwszych dwóch momentach. W szczególnym przypadku, gdy korelacja między zmiennymi losowymi T   jest liniowa (lub w przybliżeniu liniowa), dobrymi miarami tych za­leżności stochastycznych są współczynniki korelacji rmh  mię­dzy zmiennymi Th  i Tm .  Wówczas związki (33)  można przed­stawić na przykład w postaci

rmh   =  bmh     m, h = 1, 2,…, m,                                                                     (34)

gdzie  bmh  jest wartością współczynnika rmh .  Wartości te można w wielu przypadkach wyznaczyć analitycznie. W tym celu, na podstawie zbudowanych niezawodnościowych modeli odpowied­nich PK, trzeba uprzednio wyznaczyć zależności wielkości Tm i Th od różnych czynników konstrukcyjnych, technologicznych i eksploatacyjnych. Z relacji (3) – (5) wynika, że zmienne Tm i Th  są funkcjonałami, które ogólnie można przedstawić w po­staci

Tm  = fm [eO, Ga(t)] ,   gdzie  0 £ t £ Tm                                                          (35)

Th = fh [eO, Ga(t)] ,    gdzie  0 £ t £ Th                                                         (36)

Przyporządkowują one realizacjom  eO  początkowego stanu technicznego i realizacjom Ga(t)   oddziaływań zewnętrznych, występujących aż do chwili pojawienia się niesprawności, rea­lizacje Tm i Th   czasów bezawaryjnej pracy PK o numerach m i h.  Postacie tych funkcjonałów wynikają z postaci nie­zawodnościowych modeli tych PK. Ponieważ zawierają one zbiory wspólnych zmiennych losowych i zbiory wspólnych procesów losowych, więc zmienne Tm i Th są zależne stochastycznie. Zwykle wielkości Tm i Th   udaje się przedstawić jako funkcje tylko zmiennych losowych (zwłaszcza, gdy oddziaływania zewnętrzne potraktuje się jako procesy sta­cjonarne). Wówczas relacje (35) i (36) można zastąpić relacja­mi

Tm = fm ( Y1, Y2,…,Yy; L1, L2,…, Ll )  ,                                                          (37)

Th = fh ( Y1, Y2,…,Yy; V1, V2,…, Vv )  ,                                                          (38)

gdzieL1, L2,…, Ll   oraz V1, V2,…, Vv są zbiorami różnych zmiennych losowych lub zdeterminowanych opisujących początko­wy stan techniczny eO i oddziaływania zewnętrzne Ga ,  a Y1, Y2,…,Yy jest zbiorem wspólnych zmiennych losowych opisu­jących eO i Ga.  Wykorzystując relacje (37) i (38). można wyznaczyć współczynnik korelacji  rmh na podstawie znanego wyrażenia

(39)

Zasadnicze uproszczenie, pozwalające wyznaczyć przybliżo­ną wartość  rmh ,  polega zwykle na rozwinięciu funkcji fm  i fh  oraz ich iloczynu fm fh  w szeregi Taylora w otoczeniu punktów określonych przez zbiory współrzędnych odpowiednio ZWm( EY1, EY2,…, EYy; EL1 ,EL2,…, ELl ),

ZWh( EY1, EY2,…, EYy; EV1 ,EV2,…, EVV ) oraz ZWmh( EY1, EY2,…, EYy; EL1 ,EL2,…, ELl ;

EV1 ,EV2,…, EVV )  i zachowaniu w tych szeregach pierwszych kilku wyrazów. Uproszczenie to ułatwia wyznaczenie momentów zmiennych losowych, występujących w formule (39). Mimo tego uproszczenia wyprowadzenie wzoru na współczynnik rmh , jest na ogół niełatwe. Jeżeli wspomniane szeregi są wolnozbieżne, to w celu zwiększenia dokładności wyznaczenia rmh  należy przy obliczaniu momentów zmiennych Tm , Th  i Tm Th  uwzględnić większą liczbę wyrazów rozwinięć. To jednak znacznie utrudnia obliczenia.

W trudniejszych przypadkach do wyznaczania współczynników korelacji może być użyta metoda symulacji. W tym przypadku po­lega ona na losowaniu realizacji zmiennych losowych Y1, Y2,…, Yy; L1 ,L2,…, Ll ; V1 ,V2,…, Voraz wyznaczeniu na podstawie relacji (37) i (38) realizacji zmiennych Tm i Th    Odpowiednio liczny zbiór par realizacji (Tm , Th)  umożliwia wy­znaczenie estymatorów momentów występujących w związku (39) oraz obliczenie wartości odpowiedniego współczynnika korela­cji. Metoda symulacji może ułatwić wyznaczanie tych współczyn­ników korelacji również i w tych przypadkach, gdy poszczegól­nym PK przyporządkowane są większe liczby cech zdatności niż 1.

Jeżeli zależności stochastyczne z założenia nie występu­ją, do opisu szeregowej struktury niezawodnościowej obiektu złożonego wystarcza relacja (32) (lub odpowiednio inne dla innych struktur).

Nie jest również potrzebne wyrażanie wspomnianych zależności stochastycznych na przykład za pomocą relacji (33), gdy do badań niezawodności, wykorzystujących niezawodnościowy mo­del obiektu, użyta jest metoda symulacji stanów obiektu i otoczenia oraz zjawisk fizycznych prowadzących do niesprawno­ści (p. podrozdz. 5.3 [tej pracy inżynierskiej]). Wówczas bowiem zależności te’ są uwzględniane pośrednio przez losowanie między innymi wartości zmiennych losowych, wspólnych dla różnych PK, zmiennych cha­rakteryzujących początkowy stan techniczny  eO i oddziaływa­nia zewnętrzne Ga  na pojedynczy egzemplarz obiektu. Decydu­ją one o stochastycznych zależnościach między cechami zdatności, np. między Zm(t)  i Zh(t),  co wynika z relacji (3) i (4)  lub między czasami bezawaryjnej pracy poszczególnych PK, np. między Tm  i Th,  co wynika z relacji (35) i (36). Taki­mi zmiennymi losowymi, wspólnymi dla PK o numerach m i h są wielkości Y1,Y2,…,Yy  występujące w relacjach (37) i (38).

Matematyczne opisy innych struktur niezawodnościowych niż szeregowa (jeśli były przyjęte w modelu nominalnym) mogą być dokonane w sposób podobny do przedstawionego powyżej.

Przedstawiony w niniejszym rozdziale w punktach a – g zbiór modeli częściowych tworzy niezawodnościowy model obiektu me­chanicznego, nadający się do teoretycznych badań niezawodności tego obiektu. W rozdziale tym zostały zaprezentowane zasady i kolejność postępowania przy tworzeniu pełnej postaci niezawod­nościowego modelu obiektu mechanicznego. W praktyce mogą być stosowane również i inne postacie niezawodnościowego modelu takiego obiektu, będące jednakże szczególnymi przypadkami mo­delu przedstawionego. Różnice między nimi wynikają przede wszystkim z jakości i ilości posiadanych informacji o opisa­nych powyżej procesach i zjawiskach prowadzących do niespraw­ności.

Te inne postacie modelu mogą się różnić od przedstawionej postaci pełnej przede wszystkim brakiem niektórych części mo­delu. W pewnych przypadkach na przykład wśród danych wykorzy­stywanych przy tworzeniu modelu są dane charakteryzujące bez­pośrednio obciążenia zewnętrzne obiektu. Pozwalają one zwykle zbudować model oddziaływań zewnętrznych, a część modelu niezawodnościowego dotycząca procesu eksploatacji staje się zbędna. Posiadanie takich danych umożliwia więc pewne uproszczenie niezawodnościowego modelu obiektu.

Inna szczególna postać przedstawionego modelu powstaje wówczas, gdy badania dotyczą tylko pojedynczego PK. Wśród po­trzebnych danych są wtedy zwykle dane charakteryzujące bezpo­średnio uogólnione wytężenia (np. naprężenia lub naciski wy­stępujące w badanym PK. W niezawodnościowym modelu tego PK zbędne stają się te części pełnej postaci modelu, które doty­czą: procesu eksploatacji obiektu zawierającego ten PK, proce­su oddziaływań zewnętrznych na ten obiekt, procesu oddziaływań wewnętrznych. Zdecydowana większość krajowych i światowych pu­blikacji z zakresu niezawodnościowego modelowania obiektów me­chanicznych dotyczy takich właśnie szczególnych postaci mode­li, tzn. modeli PK.

Porównanie metod kontroli przeciążenia

[1, 3, 10]

W rozdziale tym porównam trzy algorytmy kontroli przeciążenia. W celu porównania algorytmów skorzystam z symulacji przeprowadzonych przez niezależne organizacje. Symulacje zostały przeprowadzone za pomocą własnych pakietów symulacyjnych. Chcąc powtórzyć podobne symulacje lub chcąc rozszerzyć zakres symulacji należałoby albo napisać do tego celu program albo zakupić odpowiedni pakiet symulacyjny. Jednym z takich komercyjnych programów umożliwiających przeprowadzenie symulacji metod kontroli przeciążenia jest pakiet OPNET firmy MIL3. Inne znane dostępne pakiet nie pozwalają na zbyt dużą ingerencję w mechanizmy kontroli przeciążenia.

Do symulacji trzech algorytmów zostanie użyta konfiguracja sieci składająca się z pięciu przełączników i z siedmiu par: nadawca-odbiorca. Odległość pomiędzy dwoma sąsiednimi przełącznikami wynosi 1000km a między przełącznikiem a węzłem sieci 100m. Wszystkie łącza mają przepustowość równą 45Mb/s.

Uwzględniając odległość połączenia w symulacji możemy wyróżnić trzy rodzaj połączeń:

Długodystansowe: połączenie od nadawcy 0 przechodzi przez cztery przełączniki. Połączenie to oznaczymy jaki VC0.

Średniodystansowe: połączenia od nadawcy 1 i 4 składają się z dwóch przełączników. Połączenia oznaczamy odpowiedni VC1 i VC4.

Krótkodystansowe: połączenia od nadawcy 2, 3, 5 i 6 przechodzą tylko przez jeden przełącznik. Połączenia oznaczamy odpowiedni VC2, VC3, VC5  i VC6.

Konfigurację testową przedstawia Rysunek 12.

Symulacja dotyczyć się będzie jednego z ważniejszych aspektów  działania algorytmu, tzn. „sprawiedliwego” rozdziału dostępnego pasma pomiędzy poszczególne kanały wirtualne. Używając metody sugerowanej przez ATM Forum możemy obliczyć optymalne teoretyczne pasmo dla każdego połączenia (współczynnik fairshare). Obliczenia te zostały przedstawione w Tabela 3

VC0 VC1 VC2 VC3 VC4 VC5 VC6
3.75 3.75 7.5 3.75 3.75 7.5 11.25

Tabela 3 Oczekiwane pasma przepustowe dla każdego połączenia w Mb/s

Symulacje zostały przeprowadzone dla dwóch parametrów ICR (initial cell rate), czyli początkowej prędkości nadawania, równych PCR i PCR/20.

 Rysunek 12. Testowa konfiguracja sieci05

Niezawodność strukturalna szlifierki

Niezawodność strukturalna szlifierki jest określona wzorem:

08

R0(t) – niezawodność strukturalna (prawdopodobieństwo poprawnej pracy)

wiertarki

Ri(t) – niezawodność i-tego elementu

3a. Zebranie danych niezawodnościowych o elementach wiertarki.

Uzyskane z literatury dane o rozkładach [11, 17] i parametrach funkcji rozkładu [3, 4, 5] dla elementów wiertarki zestawiono w tabeli 2.

3b. Zebranie danych o uszkodzeniach powstałych w trakcie eksploatacji

Dane o uszkodzeniach lin stalowych obciążonych podobnie jak lina wiertarki zestawiono w tabeli 3.

 Zestawienie danych o uszkodzeniach lin stalowych (tabela 3)

Czas poprawnej pracy liny

ti

Liczba uszkodzeń w zadanym przedziale czasu

ni

Liczba uszkodzeń do czasu ti

Ni

Ni/åni
400

800

1200

1600

1880

2280

2800

3600

4000

4100

6200

6400

4

2

3

2

2

1

1

1

1

1

1

1

4

6

9

11

13

14

15

16

17

18

19

20

0,20

0,30

0,45

0,54

0,64

0,70

0,76

0,83

0,85

0,90

0,95

1,00

å ni 20
  1. Wyznaczenie modelu matematycznego rzeczywistego rozkładu uszkodzeń.

Do wyznaczenia modelu rozkładu zastosowano metodę graficzną. Sporządzoną siatkę funkcyjna rozkładu wykładniczego [5].

Dla liniowej interpolacji danych eksploatacyjnych punkty o współrzędnych (tab. 3) naniesiono na siatkę i poprowadzono linię prostą tak, aby odchylenia punktów od prostej były najmniejsze (rys. 11).

Mała odległość punktów od prostej stanowi podstawę do wstępnego stwierdzenia zgodności rzeczywistego rozkładu uszkodzeń z rozkładem wykładniczym.

  1. Wyznaczenie parametrów otrzymanego rozkładu

Intensywność uszkodzeń l dla rozkładu wykładniczego obliczono na podstawie wykresu

w funkcji t (rys. 11) przy wykorzystaniu tabel siatki prawdopodobieństwa rozkładu [5].

gdzie:

a – kąt zawarty między prostą interpolacyjną a osią odciętych (rys. 11)

L – długość odcinka odpowiadającego obszarowi zmienności zmiennej niezależnej

t [mm],

Dt – obszar zmienności zmiennej t, [h]

Kt – współczynnik skali na osi odciętych [mm/h].

Zestawienie danych i wyników obliczeń do wyznaczenia parametru l(t) (tabela 4)

Wielkości odczytane z rysunku Wyniki obliczeń
a[°] L[mm] tmax[h] tmin[h] Dt[h] tga Kt[mm/h] l[10-6h-1]
41°36’ 160 6400 0 6400 0,8875 0,025 500
  1. Weryfikacja przyjętego modelu matematycznego rzeczywistego rozkładu uszkodzeń elementów.

Zgodność rozkładu doświadczalnego z wykładniczym stwierdzono przy pomocy testu Kołmogorowa; polegającego na sprawdzeniu warunku zgodności:

gdzie:

D – maksymalna odległość punktów doświadczalnych od prostej interpolacyjnej, liczona według osi rzędnych siatki funkcyjnej,

n – liczba doświadczalnie otrzymanych punktów.

Jeżeli to istnieje zgodność rozkładu doświadczalnego z założonym rozkładem teoretycznym.

Z rysunku 11 odczytano D = 0,003, przy n = 20. a więc rozkład doświadczalny odpowiada rozkładowi wykładniczemu.

  1. Obliczenie wartości funkcji niezawodności poszczególnych elementów urządzenia:

Wartości funkcji niezawodności Ri(t) w przedziale czasowym 500-5000 [h] wyzna­czono z tabel 2.1; 3.1; 4.2 [23], opierając się na parametrach rozkładu uszkodzeń elemen­tów i zestawiono w tabeli 5.

Parametr a dla bębna, obliczono z zależności analitycznej [23]

Przebiegi funkcji niezawodności Ri(t) przedstawiono na rysunku 12.

  1. Obliczenie niezawodności strukturalnej wiertarki:

Niezawodność strukturalną wiertarki R0(t) określoną wzorem w punkcie 2, obliczono na podstawie wyników zestawionych w tabeli 5.

Wartości niezawodności strukturalnej w przedziale czasowym 500-5000 [h] zestawio­no w tabeli 6.

Niezawodność strukturalna wiertarki w przedziale czasowym 500 – 5000 h

 Czas pracy t [h] 500 1000 1500 3000 5000
Niezawodność strukturalna R0(t) 0,5686 0,2680 0,1203 0,0082 0

Obliczenia niezawodnościowe prowadzone w toku projektowania dają możliwość ko­rygowania projektu na poszczególnych jego etapach w aspekcie wymaganej niezawodności a więc zmniejszają ryzyko zaprojektowania urządzenia którego niezawodność jest nie­wystarczająca. Na schemacie rysunku 14 przedstawiono kolejność czynności projektowych i towarzyszące im obliczenia niezawodności. W procesie projektowania wyodrębniono trzy stopnie: pierwszy, obejmujący etapy I, II, III oraz drugi i trzeci odpowiadające eta­pom IV i V.

Każdy stopień zamyka porównanie niezawodności obliczonej R0(t) niezawodnością wymaganą Rw(t). W przypadku, gdy oszacowana niezawodność urządzenia jest nie mniej­sza od wymaganej, należy podjąć działania objęte następnym stopniem projektowania. Jeżeli natomiast niezawodność obliczona jest niewystarczająca, należy: w ramach nie­zmienionej konstrukcji zastosować elementy o wyższych wskaźnikach lub, o ile znajdzie potrzeba, wyeliminować słabe ogniwa. O wyborze kolejności wykorzystania środków zmierzających do podwyższenia niezawodności decyduje wartość wskaźnika kolejności przejścia k, przypisana liniom obrazującym następstwa czasowe działań projektowych.