ANALIZA UKŁADU STERUJĄCEGO PARAMETRAMI SPAWANIA LASEROWEGO STALOWYCH ELEMENTÓW O ZMIENNEJ GRUBOŚCI

Pełen tekst

(1)MODELOWANIE INŻYNIERSKIE nr 52, ISSN 1896-771X. ANALIZA UKŁADU STERUJĄCEGO PARAMETRAMI SPAWANIA LASEROWEGO STALOWYCH ELEMENTÓW O ZMIENNEJ GRUBOŚCI Adam Kulawik1a, Joanna Wróbel1b, Michael Nikolaevich Seidurov2c 1. Instytut Informatyki Teoretycznej i Stosowanej, Politechnika Częstochowska Department of Small Business in Welding Production, I.I. Polzunov Altai State Technical University a adam.kulawik@icis.pcz.pl, bjoanna.wrobel@icis.pcz.pl, cseidurov@mail.ru 2. Streszczenie Zastosowano sztuczne sieci neuronowe do wspomagania pracy laserowego automatu spawalniczego do spawania długich prostych elementów stalowych o zmiennej grubości. Nauczona sieć neuronowa pozwoli na zwiększenie poziomu automatyzacji procesu i autonomiczną reakcję maszyny na zmieniające się warunki spawania. Zastosowanie sieci neuronowej spowoduje odpowiednią kontrolę głębokości spoiny przy założonych wartościach temperatury. Do nauczenia sieci neuronowej wykorzystano dane otrzymane z obliczeń wykonanych przez autorski program komputerowy. W oprogramowaniu tym zaimplementowano numeryczne rozwiązanie równania przewodzenia ciepła z członem konwekcyjnym metodą elementów skończonych w sformułowaniu Petrova-Galerkina. W modelu uwzględniono ciepło topienia i krzepnięcia stali. Słowa kluczowe: układ sterujący, spawanie, sztuczne sieci neuronowe, model numeryczny. THE ANALYSIS OF THE CONTROL SYSTEM OF THE LASER WELDING PARAMETERS FOR STEEL ELEMENTS OF VARIABLE THICKNESS Summary The use of artificial neural networks to support the work of welding machine to weld of long straight elements of variable thickness has been presented. The use of artificial neural network, after the appropriate learning, allows to increase the level of automation of process and autonomous reaction the machines on changing welding conditions. Application of the artificial neural network allows the appropriate control of the depth of fusion for the assumed values of the temperature. To train the neural network the data obtained from calculations performed on the authorial software has been used. The numerical model of solution of conductivity equation with the convection term by using the Petrov-Galerkin finite element methods has been implemented in this software. In the model also the heat of melting and solidification have been taken into account. Keywords: control system, welding process, artificial neural network, numerical model. 1. WSTĘP Złożoność omawianej technologii spawania laserowego wymaga ciągłego monitorowania parametrów procesu oraz reakcji na zmienność warunków takiego procesu.. W celu zapewnienia dobrego jakościowo złącza spawanego, na całej długości ściegu, częstym wymogiem jest modyfikacja parametrów procesu. Pozwala to uzyskać. 112.

(2) Adam Kulawik, Joanna Wróbel, Michael Nikolaevich Seidurov wymaganą głębokość wtopienia oraz szerokość i wysokość lica. Rozwój technologii spawania pozwala na znaczą poprawę jakości spawania i zwiększenie wydajności takiego procesu. W technologii spawania można wyróżnić następujące czynniki wpływające na redukcję kosztów: analizę zużycia materiałów i akcesoriów spawalniczych, unikanie przewymiarowania spoin, odpowiedni dobór średnicy elektrody, wykorzystanie dwustronnego ukosowania oraz wdrażanie mechanizacji, automatyzacji i robotyzacji procesów spawania [8]. Wiele z tych czynników może być uwzględnianych poprzez odpowiedni dobór parametrów procesu, gdy bazuje się na wynikach symulacji numerycznych. Zastosowanie obliczeń numerycznych pozwala poprawić technologię i zminimalizować występujące odkształcenia spawalnicze, zmniejszyć kosztochłonność materiałową, uzyskać dokładniejszą ocenę parametrów procesu spawania, wspomóc pracę robotów spawalniczych itp., dlatego często analizuje się procesy w trakcie których występują wady produkcyjne. Rozwiązaniem może być dokładniejsza kontrola przebiegu samego procesu z zastosowaniem modeli numerycznych opierających się na metodzie elementów skończonych lub metodzie różnic skończonych [6,7]. Jednak analiza numeryczna zagadnień niestacjonarnych o dużej zmienności parametrów jest wymagająca ze względów czasowych oraz sprzętowych. Stosowane współcześnie modele numeryczne pozwalają zarówno na wyznaczenie parametrów procesu a priori, jak i na kontrolę procesu poprzez aplikacje czasu rzeczywistego. Modele numeryczne wykorzystywane do poszukiwania optymalnych rozwiązań, często opatrzonych pewną niepewnością, coraz częściej bazują na metodach sztucznej inteligencji. Zastosowanie algorytmów genetycznych do optymalizacji procesów, ze względu na dużą złożoność obliczeniową, nie pozwala na reagowanie układu sterującego w czasie rzeczywistym. Natomiast metody oparte na zastosowaniu w układach sterujących sztucznych sieci neuronowych są jednymi z najszybszych [4,5].. Rys. 1. Schemat modelu. W prezentowanym modelu założono, że sieć neuronowa jest nauczona na podstawie danych otrzymanych z rozwiązania równania ustalonego przewodzenia ciepła (rys. 1). Wykonano 2000 symulacji dla różnych kombinacji prędkości przesuwu źródła oraz grubości spawanego elementu. Zakresy losowanych wartości były następujące dla: grubości [0.0026, 0.0053] m, prędkości [0.0012, 0.0075] m/s. Ze zbioru uczącego oraz testowego zostały wykluczone przypadki, w których maksymalna temperatura na górnej powierzchni nagrzewanego elementu była mniejsza od 1200 K lub większa niż 5000 K. Wartości wylosowanych prędkości oraz grubości przedstawione na rysunkach 2 i 3 wskazują na poprawną liczebność danych w rozważanych przedziałach. Pozostałe parametry źródła ciepła, takie jak moc (Q=1600 W) i promień (r0=0,001 m), przyjęto stałe. Własności materiałowe przyjęto jak dla stali średniowęglowej (C45).. 2. PODSTAWOWE ZAŁOŻENIA MODELU Podczas procesu spawania długich elementów stalowych o zmiennej grubości, gdy wykorzystuje się np. technologię spawania laserowego, mogą wystąpić wady połączenia polegające na zbyt małej głębokości wtopienia lub przetopienia i przepalenia miejsca połączenia. Zastosowanie sieci neuronowej pozwoli na zwiększenie poziomu automatyzacji procesu i autonomiczną reakcję maszyny spawającej na zmieniające się warunki spawania.. Rys. 2. Uzyskane poziomy grubości spawanego elementu w zbiorze uczącym. 113.

(3) ANALIZA UKŁADU STERUJĄCEGO PARAMETRAMI SPAWANIA LASEROWEGO... LASEROWEGO. Rys. 5.. Położenie czujników temperatury (zadanie stacjonarne). Rys. 3. Uzyskane poziomy prędkości źródła ciepła w zbiorze uczącym. Przyjęcie powyższych założeń dla prezentowanego zadania pozwala na wyznaczenie profili temperatury dla przesuwającego się źródła ciepła w różnych warunkach warunk brzegowych i początkowych. Wykorzystanie do określeokreśl nia zbiorów uczącego oraz testowego obliczeń z zadania stacjonarnego znacznie skraca czas obliczeń. oblic W modelu tym otrzymuje się rozkłady temperatury dla każdego przypadku, wykonując jedynie kilka iteracji. Wynika to z uwzględniania ciepła topienia i krzepnięcia dla obszaru jeziorka spawalniczego.. Rys. 6. Usytuowanie i położenie ołożenie czujników. Działanie nauczonej sieci neuronowej neuronow zostało zweryfikowane w modelu numerycznym spawania elementu stalowego o zmiennej grubości (zadanie niestacjonarne). nie Założono,, że prędkość przesuwu źródła ciepła jest wyznaczana na podstawie przyjętego profilu temperatury. Ponadto dokonano pomiaru grubości przed przesuwająprzesuwaj cym się źródłem ciepła w odległości 0.0025 m od środka źródła. Sieć ieć neuronowa poprzez bieżącą modyfikację prędkości i danych z pomiaru grubości utrzymywała założony profil temperatury.. 3. MODELOWANIE ZJAWISK CIEPLNYCH PROCESU CESU SPAWANIA Do modelowania zjawisk cieplnych dla ustalonego przepływu ciepła wykorzystano rozwiązanie równania różniczkowego z członem konwekcyjnym (współrzędne Eulera): ∙ ∙

(4)

(5) Rys. 4. Uzyskane maksymalne poziomy temperatury (wskazane przez czujniki) w zbiorze uczącym. (1). gdzie: T [K] - temperatura, λ

(6) λT [W/mK] - współczynnik przewodzenia ciepła, [J/m3K] - ciepło efektywne, [m/s] - wektor prędkości, q [W/m3] - objętościowe źródło ciepła. adanie rozwiązano metodą elementów skończonych Zadanie w sformułowaniu Petrova-Galerkina Galerkina. Dla zadania stacjonarnego przepływu ciepła układ równań MES ma postać [6, 7]. Z otrzymanych wyników zostały wyodrębnione warwa tości temperatury dla czterech grup czujników (rys. 5 i 6). Na podstawie wykonanych obliczeń dla różnych grup czujników przyjęto, rzyjęto, że parametrami wejściowymi do sieci będą temperatury z czujników o numerach 17-23 17 umieszczone na dolnej powierzchni spawanego elementu oraz jego grubość. Wybór trzeciego zestawu czujników temperatury jako danych wejściowych został dokonany ze względu na możliwość utrzymania założonej temperatury na dolnej powierzchni spawanego elementu (rys. 9). Uzyskane maksymalne wartości temperatury dla trzeciej grupy czujników przedstawiono na rys. rys 4. Na wyjściu sieci uzyskuje się prędkości źródła spawającego.. . K B T

(7) B T B q∗ Q q. (2). gdzie elementowe macierze: przewodności K , warunków brzegowych B oraz źródeł wewnętrznych Q są następujące K

(8) ౛ λw, ϕ, dΩ ρc ౛ w ϕ, V ϕ dΩ. 114.

(9) Adam Kulawik, Joanna Wróbel, Michael Nikolaevich Seidurov Q = ౛ w ϕ dΩ . (3). . B = B + B = ౛ w ϕ dΓ + ౛ α w ϕ dΓ . ్. ొ. gdzie: w są funkcjami wagowymi, ϕ są funkcjami interpolacyjnymi. Do weryfikacji działania nauczonej sieci neuronowej zastosowano model nieustalonego przewodzenia ciepła we współrzędnych Lagrange'a ∙ − .

(10)

(11) . =. (4). = −

(12) . gdzie: [W] moc wiązki laserowej, ) [m] promień wiązki, [m] bieżący promień, [m] bieżąca głębokość penetracji wiązki laserowej, [m] głębokość penetracji wiązki laserowej wyznaczana z zależności [2, 3] "ೌ #. ಴೔ ೔షభ ೔సభ ೔ . . *+ ೛ ಮ , ∑ర. (10). gdzie: / - współczynnik absorpcji, 0 [K] - temperatura parowania, [K] - temperatura otoczenia, - liczba Pecleta.. gdzie: t[s] jest czasem. Równanie różniczkowe (2) rozwiązano z zastosowaniem metody elementów skończonych w sformułowaniu Galerkina. Przyjmując, że parametr γ <0,1> decyduje o przyjętym schemacie całkowania po czasie, układ równań MES przyjmuje postać: γK + B + M T = M − 1 − γ K + B T. . +γ B T. . . ∗. − B q. +1 − γ B T. . +. + Q q. ∗. − B q. +. . + Q q. (5). 4. ZASTOSOWANIE SZTUCZNYCH SIECI NEURONOWYCH W ELEMENTACH STERUJĄCYCH. . gdzie elementowe macierze: przewodności K , pojem-. ności cieplnej M , źródeł wewnętrznych Q oraz warunków brzegowych B są następujące: K = ౛ λw, ϕ, dΩ M =. . ∆ ౛. (6). w ϕ dΩ. Q = ౛ w ϕ dΩ . . B = B + B = ౛ w ϕ dΓ + ౛ α w ϕ dΓ . ్. ొ. W modelu krzepnięcia ciepło przemiany fazowej uwzględniono za pomocą modyfikacji ciepła efektywnego (sformułowanie pojemnościowe) [1]: C = ρT cT − ρ L . ౏ ై . . Rys. 7. Rozkłady mocy źródła ciepła na dwóch poziomach: z=0, z=2 mm. Budowa sieci neuronowej (oparta na schemacie logicznym działania komórki nerwowej) sprawia, że jest to jeden z najszybciej działających układów obliczeniowych dla zadań silnie nieliniowych. Struktura ta pozwala na wyznaczenie poprawnych wartości parametrów także dla zadań podobnych, ale niewystępujących jako przykłady wzorcowe. Zastosowano sztuczną sieć neuronową wielowarstwową (rys. 8) zbudowaną z następujących warstw: jednej wejściowej (9 neurony), dwóch ukrytych (9 i 5 neuronów) oraz jednej wyjściowej (1 neuron).. (7). ౏. gdzie: [kg/m3] gęstość fazy stałej, [J/kg] ciepło krzepnięcia, udział fazy stałej, , [K] są odpowiednio temperaturą likwidusu oraz solidusu. Udział fazy stałej wyznaczano z zastosowaniem reguły dźwigni f = f T =. ై ై ౏. , T ∈ T! , T . (8). W modelach zjawisk cieplnych (przepływ ustalony i nieustalony) do symulacji procesu spawania zastosowano źródło objętościowe uwzględniające głębokości penetracji (, ) wiązki laserowej (rys. 7), tzn. [3]:

(13) , =. "ೌ # $%బమ &. % '. (. %బ. &. exp 1 − 1 − . (9) Rys. 8. Zastosowana sztuczna sieć neuronowa. 115.

(14) ANALIZA UKŁADU STERUJĄCEGO PARAMETRAMI SPAWANIA LASEROWEGO... LASEROWEGO Do uczenia sieci wykorzystano algorytm wstecznej propagacji błędu oraz w celu zwiększenia prędkości uczenia uwzględniono momentową metodę wstecznej propagacji błędu [4,5].. 5. WYNIKI OBLICZEŃ Wykonano weryfikację działania sieci neuronowej dla elementu stalowego o zmiennej grubości. Zostały przeprz prowadzone obliczenia dla dwóch wymaganych profili temperatury (rys. 9-11) dla czujników temperatury 1717 24.. Rys. 12.. Zmiana prędkości źródła ciepła, zadanie niestacjonarne przypadek 1. Rys. 9. Założony profil il temperatury, przypadek 1 i 2. Rys. 13.. Zmiana prędkości źródła ciepła, zadanie zad niestacjonarne przypadek 2. Rys. 10.. Rozkład temperatury dla położenia źródła ciepła x=0.023 m, przypadek 1. Rys. 11.. Rozkład temperatury dla położenia źródła źródł ciepła x=0.023 m, przypadek 2. Przyjęto, że źródło ciepła znajdowało się za prawą krawędzią elementu spawanego i przesuwało się z prędkością początkową Vx=-0.005 0.005 m/s. Następnie na postapost wie zmiany grubości sieć neuronowa aktualizowała prędkość przesuwu (rys. 12 i 13), powodowało to zmiany temperatury (rys. 14 i 15). Dla obu przypadków głębokość penetracji wyznaczowyznacz na na podstawie wzoru (10) nie przekraczała grubości elementu stalowego (rys. 16).. Rys. 14.. Zmiany temperatury w wybranych czujnikach, zadanie zad niestacjonarne przypadek 1. Zmiana prędkości źródła ciepła pozwoliła na dosyć skuteczną kontrolę temperatury maksymalnej na pop wierzchni górnej i dolnej spawanego elementu (rys. 12 15). Sieć neuronowa dobrze reagowała na zmianę grubogrub ści spawanego elementu (rys. 10-18). 18).. 116.

(15) Adam Kulawik, Joanna Wróbel, Michael Nikolaevich Seidurov. Rys. 17. Położenie obszaru przetopionego, przypadek 1. Rys. 15. Zmiany temperatury w wybranych czujnikach, zadanie niestacjonarne przypadek 2 Rys. 18. Położenie obszaru przetopionego, przypadek 2. 6. WNIOSKI Zastosowanie sztucznej sieci neuronowej do kontroli parametrów procesu spawania laserowego elementów stalowych o zróżnicowanej grubości jest skuteczną techniką. Zastosowanie takiego układu wraz z czujnikiem grubości umożliwia pełna automatyzację procesu spawania. Należy jednak zaznaczyć, że kalibrację sztucznej sieci neuronowej należy wykonać dla każdej zmiany gatunku stali czy zakresu grubości elementów. Prezentowana metoda kontroli parametrów spawania laserowego jest skuteczna zarówno w przypadkach całkowitego jak i częściowego przetopienia elementu stalowego (rys. 17 i 18). Zastosowanie w urządzeniach spawalniczych elementów, bazujących na sztucznej inteligencji, takich jak sieci neuronowe, pozwoli na redukcję kosztów.. Rys. 16. Głębokość penetracji wiązki laserowej, przypadek 1 i 2. Zmiany głębokości przetopienia na granicach obszarów o różnej grubości spowodowane są tym, że pomiar grubości dokonywany jest tylko przed źródłem ciepła.. Literatura 1. 2. 3. 4. 5. 6. 7. 8.. Bokota A.: Modelowanie krzepnięcia i stygnięcia dwuskładnikowych stopów metali: pola temperatury, stężeń i naprężeń. Częstochowa: Pol. Częstochowska, 2001. Monografie nr 79. Lankalapalli K.: Model-based laser weld penetration depth estimation. „Welding In the World” 1997, 6, Vol. 39, p. 304-313. Piekarska W.: Analiza numeryczna zjawisk termomechanicznych procesu spawania laserowego: pole temperatury, przemiany fazowe i naprężenia. Częstochowa: Pol. Częstochowska, 2007. Monografie nr 135. Rojas R.: Neural network: a systematic introduction. Berlin: Springer, 1996. Rutkowski L.: Computational intelligence: methods and techniques. Berlin: Springer, 2005. Wait R., Mitchell A. R.: Finite element analysis and applications. Chichester: John Wiley & Sons, 1985. Zienkiewicz O. C., Taylor R. L.: The finite element method. Oxword: Butterworth-Heinemann, 2000. Reducing your welding costs: http://www.lincolnelectric.com/en-us/support/process-and-theory/pages/reducingwelding-costs-detail.aspx. 117.

(16)