Przeprowadzono serię symulacji MES z wykorzystaniem różnych kombinacji długości przemieszczenia łuku, prędkości posuwu łuku oraz wypadkowej szybkości doprowadzania energii. W zbiorze treningowym zakres wartości zastosowanych dla długości przemieszczenia łuku wynosił 40 mm, 50 mm, 60 mm, 70 mm i 80 mm. Zakres wartości zastosowanych dla prędkości posuwu łuku wynosił 2 mm/s, 2,5 mm/s, 3 mm/s, 3,5 mm/s i 4 mm/s. Zakres wartości zastosowanych dla wypadkowej szybkości doprowadzania energii wynosił 1000 W, 1500 W, 2000 W, 2500 W i 3000 W. W zbiorze testowym zakres wartości zastosowanych dla długości przemieszczenia łuku wynosił 40 mm, 50 mm, 60 mm, 70 mm i 80 mm. Zakres wartości zastosowanych dla prędkości posuwu łuku wynosił 2,25 mm/s, 2,75 mm/s, 3,25 mm/s i 3,75 mm/s. Zakres wartości zastosowanych dla wypadkowej szybkości doprowadzania energii wynosił 1250 W, 1750 W, 2250 W i 2750 W. Naprężenia wzdłużne wzdłuż linii BD dla tych 205 symulacji zostały zebrane i przedstawione na Rysunku 2. Wynik ten został wcześniej opublikowany w25. W szczególności, wyniki studium przypadku zostały wykorzystane w tamtej pracy do potwierdzenia skuteczności przepływu pracy. Niniejsze badanie wykracza poza wcześniejsze analizy, dostarczając szczegółowego wyjaśnienia konkretnych ustawień zastosowanych w symulacji spawania oraz szczegółowo opisując metodologię skryptowania w języku Python wykorzystaną do ogólnej symulacji spawania.
Opracowany model ANN zastosowano do przeprowadzenia predykcji na zbiorze testowym. Rozbieżność między predykcją modelu a symulacją FE przedstawiono na Rysunku 3. Spośród 1680 punktów testowych większość różnic (759 przypadków) między wartościami naprężeń symulowanymi metodą analizy elementów skończonych a wartościami przewidzianymi przez ANN mieści się w zakresie 0-2 MPa, co stanowi 45,2% danych. Częstotliwość maleje wraz ze wzrostem różnicy, a bardzo niewiele punktów danych wykazuje różnice większe niż 10 MPa. Na Rysunku 4 przedstawiono wykres pudełkowy obrazujący rozkład bezwzględnych rozbieżności między naprężeniem symulowanym a przewidzianym dla różnych poziomów naprężenia symulowanego. Każde pudełko reprezentuje rozstęp międzykwartylowy reszt w konkretnym przedziale poziomów naprężeń symulowanych, a linia wewnątrz pudełka wskazuje wartość mediany. Histogram na pomocniczej osi y pokazuje liczbę punktów danych treningowych w każdym przedziale poziomów naprężeń, dostarczając dodatkowych informacji o rozkładzie danych wykorzystanych do uczenia modelu. Jak pokazano na rysunku, przedziały poziomów naprężeń z mniejszą liczbą punktów danych treningowych mają tendencję do wykazywania wyższych maksymalnych rozbieżności w zbiorze testowym. Wydajność w tych obszarach mogłaby zostać potencjalnie poprawiona poprzez zwiększenie ilości dostępnych danych treningowych dla tych poziomów naprężeń.
Do ilościowego określenia zdolności prognostycznej wykorzystano średni błąd kwadratowy (MSE), a wzór na wskaźniki wydajności MSE jest następujący:

gdzie yi to wartość naprężenia przed odwróceniem standaryzacji dla i-tego punktu próbki,
to przewidywana wartość naprężenia przed odwróceniem standaryzacji dla i-tej próbki, a n to całkowita liczba próbek.
Przewidywania sieci ANN w zakresie naprężeń były bardzo zbliżone do wyników symulacji FE, przy MSE wynoszącym 0,0024. Przewidywania dla czterech wybranych przypadków testowych przedstawiono na Rysunku 5; doprowadzone ciepło dla tych czterech przypadków wynosi odpowiednio 333 J/mm, 538 J/mm, 692 J/mm oraz 1222 J/mm. W każdym z tych czterech przypadków przewidywania ANN wykazują ogólnie dobrą zgodność z wynikami symulacji FE. Jednakże w scenariuszach, gdy prędkość spawania wynosi 2,25 mm/s przy długości spoiny 50 mm i wypadkowej mocy doprowadzanej energii 2750 W, lub gdy prędkość spawania wynosi 3,75 mm/s przy długości spoiny 40 mm i wypadkowej mocy doprowadzanej energii 1250 W, obserwuje się stosunkowo duże rozbieżności u dołu linii BD (w pobliżu głębokości 17 mm). Rozbieżności te mogą wynikać ze zwiększonej złożoności reguł rządzących tym obszarem, co prowadzi do powstania bardziej skomplikowanej rozmaitości w tych konkretnych obszarach przestrzeni parametrów. Możebnym rozwiązaniem byłoby zwiększenie gęstości zbioru danych treningowych w tym obszarze.

Rycina 1: Schemat procedury. Zaproponowany schemat składa się z czterech kroków (oznaczonych niebieskimi blokami): budowy bazowego modelu FE, programowania, generowania danych oraz opracowywania modelu zastępczego. Pomarańczowe bloki reprezentują dane wejściowe każdego kroku, a zielone bloki reprezentują dane wyjściowe każdego kroku. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 2: Wygenerowany zbiór danych. Przedstawiono łącznie 205 zestawów danych dotyczących naprężeń wzdłużnych, z których każdy odpowiada unikalnej kombinacji długości przesuwu, wypadkowej szybkości doprowadzania energii i prędkości spawania. W każdym zestawie naniesiono 20 punktów danych dla różnych węzłów wzdłuż linii BD. Aby odróżnić poszczególne zestawy, każdy z nich został oznaczony innym kolorem, który odpowiada wypadkowej szybkości doprowadzenia ciepła. Wartość tę można obliczyć, dzieląc wypadkową szybkość doprowadzania energii przez prędkość spawania. Rycina została zmodyfikowana na podstawie opracowania25. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 3: Częstotliwość występowania rozbieżności między wynikami symulacji FE a przewidywaniami ANN. Różnice między wynikami symulacji FE a przewidywaniami ANN w zbiorze testowym zostały podzielone na 10 poziomów. Adnotacje na wykresie wskazują częstotliwość i procent występowania każdego poziomu różnicy w przypadkach testowych. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rysunek 4: Rozkład średnich reszt bezwzględnych i zbioru danych treningowych dla różnych poziomów naprężeń. Wykres pudełkowy przedstawia rozkład bezwzględnych rozbieżności pomiędzy naprężeniem symulowanym a przewidywanym dla różnych poziomów naprężeń symulowanych. Pomarańczowa linia wewnątrz każdego pudełka reprezentuje medianę, wskazując punkt środkowy rozbieżności bezwzględnych. Wysokość każdego pudełka odpowiada rozstępowi międzykwartylnemu, który obejmuje środkowe 50% danych (od 25. do 75. percentyla), odzwierciedlając zmienność rozbieżności. Wąsy rozciągają się od pudełka, pokazując minimalne i maksymalne wartości rozbieżności. Niebieski histogram przedstawia liczbę punktów danych treningowych odpowiadających każdemu poziomowi naprężeń w zbiorze treningowym. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 5: Porównanie wyników symulacji z wynikami przewidzianymi przez sieć ANN. Na rysunku przedstawiono cztery przypadki testowe: (A) prędkość spawania 3,75 mm/s, długość spawu 40 mm i wypadkowa szybkość doprowadzenia energii 1250 W; (B) prędkość spawania 3,25 mm/s, długość spawu 60 mm i wypadkowa szybkość doprowadzenia energii 1750 W; (C) prędkość spawania 3,25 mm/s, długość spawu 50 mm i wypadkowa szybkość doprowadzenia energii 2250 W; oraz (D) prędkość spawania 2,25 mm/s, długość spawu 50 mm i wypadkowa szybkość doprowadzenia energii 2750 W. Najbardziej znacząca rozbieżność we wszystkich przypadkach testowych została zaznaczona niebieską ramką, wykazując różnicę 19,4 MPa na głębokości 17 mm od górnej powierzchni w przypadku testowym (D). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Tabela 1: Informacje o geometrii próbki spoiny, obszaru zagęszczonej siatki oraz spoiny.W tabeli wymieniono długość, szerokość i głębokość prostokątnej próbki spoiny oraz obszaru zagęszczonej siatki. Wymiary spoiny zostały przedstawione na szkicu. W procesie modelowania szerokość należy podzielić przez dwa, ponieważ utworzono tylko połowę modelu. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 2: Właściwości materiałowe stali AISI 316LN14. Gęstość oraz współczynnik Poissona są traktowane jako niezależne od temperatury, natomiast pozostałe właściwości fizyczne są traktowane jako zależne od temperatury. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 3: Właściwości plastyczne AISI 316LN24. Właściwości mechaniczne stali nierdzewnej AISI 316LN w funkcji temperatury; granica plastyczności przy zerowym odkształceniu plastycznym; parametry utwardzania kinematycznego (C1, Gamma 1, C2, Gamma 2); parametr utwardzania izotropowego (Q-inf) oraz parametr utwardzania (b). Prosimy kliknąć tutaj, aby pobrać tę tabelę.
Uzupełniający plik z kodem 1: Create_input_files.py Kliknij tutaj, aby pobrać ten plik.
Uzupełniający plik z kodem 2: Run_thermal Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający z kodem 3: DFLUX.for Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający z kodem 4:Run_Mechanical Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający z kodem 5: extract_data.py Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający z kodem 6: ANN_development.py Kliknij tutaj, aby pobrać ten plik.