PS 11 2017 art 70 PRZEGLĄD  SPAWALNICTWA Vol. 89 11/2017 Wykorzystanie szumu Barkhausena w ocenie twardości  warstw nawęglanych oraz hartowanych indukcyjnie   na stali AMS 6414 Using of Barkhausen noise in hardness evaluation of carburized and inductively tempered layers on AMS 6414 steel Mgr inż. Krzysztof Fryczowski – Politechnika Śląska; dr hab. inż. Maciej Roskosz, prof. AGH – AGH Akademia Górniczo-Hutnicza;  dr inż. Dominik Kukla – Instytut Podstawowych Problemów Techniki PAN; dr inż. Maciej Szwed– Politechnika Warszawska. Autor korespondencyjny/Corresponding author: krzysztof.fryczowski@polsl.pl Streszczenie Analizowano możliwości określenia zmian twardości warstw nawęglanych i hartowanych indukcyjnie na podsta- wie wielkości ilościowo opisujących szum Barkhausena. Badania przeprowadzono na zestawie trzech próbek o sta- łej grubości warstwy nawęglanej (1,2 mm) o różnych war- tościach twardości. Pomiary przeprowadzono dla trzech konfiguracji aparatury MEB-4C Mag-Lab s.c. różniących się częstotliwością prądu magnesującego. Analizowano ob- wiednię szumu Barkhausena i jej punkty charakterystyczne, energię, RMS, periodogramy FFT oraz rozkłady liczby zda- rzeń. Uzyskano interesujące korelacje diagnostyczne mogą- ce posłużyć w rozwiązywaniu zagadnień odwrotnych badań nieniszczących. Słowa kluczowe: twardość; szum Barkhausena; zagadnienie odwrotne Abstract The possibility of hardness changes determination of carburized and induction hardened layers using quanti- ties describing Barkhausen noise quantitatively was ana- lyzed. The studies were carried out on three samples of con- stant thickness of carburized layer (1.2 mm) with different hardness values. Measurements were made for three con- figurations of MEB4-C device provided by Mag-Lab s.c. that differed in magnetization frequency. The Barkhausen noise envelope and its characteristic parameters, energy, RMS, FFT periodograms and event distributions were analyzed. Interesting diagnostic correlations have been obtained which can be used to solve the non-destructive testing inverse problems. Keywords: hardness; Barkhausen noise; inverse problems Wprowadzenie W trakcie tradycyjnego pomiaru twardości wgłębnik twar- dościomierza odkształca plastycznie fragment powierzchni, na której dokonywany jest pomiar. W przypadkach niektó- rych elementów takie uszkodzenie powierzchni jest niedo- puszczalne. W innych przypadkach odkształcenie plastycz- ne powierzchni nie stanowi problemu eksploatacyjnego, natomiast uniemożliwia przeprowadzania pomiaru twardo- ści dokładnie w tym samym miejscu, co stanowi problem diagnostyczny. Rozwiązanie tych problemów jest możliwe przez zastosowanie tzw. zagadnień odwrotnych badań nieniszczących, w których właściwości mechaniczne okre- ślane są na podstawie innych wielkości fizycznych. Celem badań zaprezentowanych w artykule była analiza możliwo- ści opracowania zależności diagnostycznych pomiędzy wiel- kościami ilościowo opisującymi szum Barkhausena a twar- dością Vickersa HV0,5 zmierzoną na powierzchni zestawu Krzysztof Fryczowski, Maciej Roskosz, Dominik Kukla, Maciej Szwed przeglad Welding Technology Review trzech próbek, ze stali AMS 6414 o różnej twardości. Bada- nia dotyczące wykorzystania szumu Barkhausena w ocenie twardości prowadzone są obecnie w ośrodkach badawczych i obejmują analizy rozkładu twardości w funkcji odległości od powierzchni próbki [1,2], rozkładu twardości w spoinie i w strefie wpływu ciepła [3] oraz wpływu odkształcenia pla- stycznego oraz wpływu eksploatacji na zmiany twardości [4]. Szczegóły badań Badania przeprowadzono na zestawie trzech próbek, ze stali AMS 6414, poddanych procesowi nawęglania, w at- mosferze o różnym stężeniu węgla, a następnie hartowaniu. W wyniku tych procesów otrzymano próbki o różnej twar- dości i stałej grubości warstwy nawęglanej równej 1,2 mm. DOI:.http://dx .doi .org/10 .26628/ps .v89i11 .837 71PRZEGLĄD  SPAWALNICTWA Vol. 89 11/2017 Powierzchnia próbek została przygotowana do badań po- przez szlifowanie [5]. Rdzeń próbek ma strukturę martenzytyczną. W przypadku warstw widoczne jest coraz większe rozdrobnienie martenzy- tu wraz ze wzrostem twardości. Wynika to z faktu zwiększenia ilości miejsc zarodkowania martenzytu podczas hartowania wraz ze wzrostem zawartości węgla w warstwie nawęglonej, co w konsekwencji prowadzi do ograniczenia jego rozrostu. Dodatkowo w najtwardszej warstwie można odnotować obec- ność wydzieleń cementytu [5]. Pomiary szumu Barkhausena wykonywano przy użyciu systemu pomiaru, analizy i rejestracji impulsów napięcio- wych natężenia Barkhausena MEB-4C firmy Mag-Lab s.c. Pomiary przeprowadzono dla kierunku magnesowania rów- noległego i prostopadłego do kierunku śladów szlifowania dla trzech konfiguracji ustawień urządzenia pomiarowego – szczegółowe informacje zawarto w tablicy I. Wynikiem pomiaru jest czasowy rozkład napięcia szumu Barkhausena UBN. Na rysunku 1 przedstawiono jeden pełny cykl przemagnesowania. Linia niebieska to przebieg prądu magnesującego, kolorem szarym zaznaczono nieprzetworzo- ny, zarejestrowany w urządzeniu pomiarowym szum Barkhau- sena, linia koloru czarnego to wyznaczona obwiednia szumu Barkhausena. Zmiany natężenia prądu magnesującego mają przebieg trójkątny symetryczny. W pełnym cyklu zmian moż- na wyróżnić połowę cyklu, w którym wartość natężenia prądu maleje (oznaczenie I ) i połowę cyklu, w którym natężenie prą- du rośnie (oznaczenie I ). Tablica I. Szczegóły konfiguracji pomiarowej Table I. Details of the measurement configurations Tablica II. Wyniki pomiarów twardości HV0.5 Table II. Results of HV0.5 hardness measurements Rys. 1. Szczegółowy opis pojedynczego cyklu przemagnesowania Fig. 1. A detailed description of a single magnetization cycle Rys. 2. Sposób wyznaczania wielkości charakterystycznych obwiedni szumu Barkhausena Fig.  2. Method of determination of BN envelopes charateristic parameters oznaczenie konfiguracji Konfiguracja C1 Konfiguracja C2 Konfiguracja C3 częstotliwość próbkowania fp [kHz] 800 800 800 amplituda prądu magnesującego [mA] 200 200 200 częstotliwość prądu magnesującego [Hz] 2,02 5,31 10,4 wzmocnienie sygnału 50 dB 50 dB 25 dB liczba cykli przemagnesowania 5 15 25 Oznaczenie  próbki Wyniki pomiarów HV0,5 Wartość średnia  HV0,5 Odchylenie  standardowe 1 757 749 746 768 748 757 740 757 757 753 753,2 7,40 2 815 800 809 789 801 809 803 814 815 803 805,8 7,82 3 834 827 827 821 851 846 821 853 846 836 836,2 11,53 Badania mikrotwardości na powierzchni warstw nawęgla- nych próbek wykonano przy pomocy mikrotwardościomierza Zwick Roel a wyniki pomiarów pokazano w tablicy II. Analiza wyników Analiza obwiedni przebiegu czasowego MBN Obwiednię otrzymano w wyniku operacji uśredniania wygładzającego na wartościach bezwzględnych impulsów napięciowych UBNi. Wykorzystano do tego celu operację wie- lokrotnego filtrowania wygładzającego przy wykorzystaniu filtru Savitzkiego-Golaya. Przeprowadzono 10 iteracji wygładzania wielomianem 3-go stopnia o różnej szerokości okna dla poszczególnych kon- figuracji (C1-1500, C2-300 próbek, C3-100 próbek). Następnie uśredniono wartości obwiedni z kilku cykli przemagnesowa- nia (szczegóły w tabl. I) odpowiadające tym samym warto- ściom prądu magnesującego. Do analizy wybrano fragmenty obwiedni odpowiadające połowie cyklu przemagnesowania, w którym wartość natę- żenia prądu maleje I . Można w nich wyróżnić lokalne eks- trema, które można scharakteryzować poprzez ich wartości odpowiadające napięciu obwiedni PEAKVAL i miejsca ich wy- stępowania odpowiadającą natężeniu prądu magnesującego PEAKPOS co graficznie zostało przedstawione na rysunku 2. ↑ ↑ ↑ Im UEN 72 PRZEGLĄD  SPAWALNICTWA Vol. 89 11/2017 0 -20 -40 -60 -80 -100 0,8 0,6 0,4 0,2 0 U EN  [V ] Im [mA] Rys.  6.  Zależność pomiędzy twardością i PEAKPOS dla ekstremum E2 – konf. C3 Fig.  6. Dependence between hardness and PEAKPOS values for extremum E2 – conf. C3 Rys.  7.  Zależność pomiędzy twardością i PEAKPOS dla ekstremum E3 – konf. C3 Fig.  7. Dependence between hardness and PEAKPOS values for extremum E3 – conf. C1 Rys.  5.  Zależność pomiędzy twardością i PEAKPOS dla ekstremum E1 – konf. C3 Fig.  5. Dependence between hardness and PEAKPOS values for extremum E1 – conf. C1 Rys.  3.  Fragmenty obwiedni szumu Barkhausena reprezentatywne dla konfiguracji C1 i C2 Fig. 3. Fragments of BN envelopes representative for C1 and C2 con- figuration Rys. 4.  Fragmenty obwiedni szumu Barkhausena reprezentatywne dla konfiguracji C3 Fig. 4. Fragments of BN envelopes representative for C3 configuration W przypadku konfiguracji C1 i C2, obwiednie charaktery- zują się występowaniem dwóch lokalnych ekstremów, nato- miast dla konfiguracji C3 trzech. Miejsca ich występowania przykładowo dla konfiguracji C1 pokazano na rysunku 3, a dla konfiguracji C3 na rysunku 4. W przypadku konfiguracji C1 i C2 zmiany PEAKPOS dla eks- tremum E1 i E2 nie umożliwiają opracowania korelacji dia- gnostycznej. Dla konfiguracji C3 w przypadku ekstremum 1 (rys. 5) wzrostowi twardości towarzyszy przesuwanie się war- tości maksymalnej PEAKVAL w kierunku niższych wartości na- tężenia prądu magnesującego, natomiast dla ekstremum E2 występuje zależność odwrotna (rys. 6). Ekstremum E3 zo- stało zaobserwowane tylko dla konfiguracji C3. Wyższym wartościom twardości odpowiada większa wartość bez- względna natężenia prądu magnesującego PEAKPOS (rys. 7). W przypadku analizy wartości maksymalnej obwiedni PEAKVAL dla ekstremum E1 dla konfiguracji C1 i C2 większym wartościom twardości odpowiadają mniejsze wartości mak- symalne obwiedni szumu Barkhausena PEAKVAL, natomiast dla konfiguracji C3 brak jednolitego trendu zmian. Przykła- dową zależność pomiędzy twardością a PEAKVAL pokazano na rysunku 8. Wartość PEAKVAL dla ekstremum E2 umożliwia opracowanie korelacji jedynie dla konfiguracji C3 – przykład na rysunku 9. Na rysunku 10 pokazano zależność pomiędzy twardością i PEAKVAL dla ekstremum E3 obwiedni szumu Barkhausena. Większym wartościom twardości odpowiada- ją niższe wartości PEAKVAL. Analiza zmian energii szumu Barkhausena Energia szumu Barkhausena zwana również jego inten- sywnością definiowana jest jako suma iloczynów kwadratu kolejnych impulsów napięciowych UBNi i kroku czasowego, wynikającego z częstotliwości próbkowania fpr, wyrażoną w postaci (1): EBN= Σ [(UBNi)2•∆t], ∆t=1/fpr (1) gdzie: N – liczba próbek sygnału szumu Barkhausena zarejestro- wanego dla jednego cyklu przemagnesowania. Na rysunkach od 11 do 13 pokazano zależności pomię- dzy twardością i energią EBN obwiedni szumu Barkhausena. Dla wszystkich kierunków magnesowania i wszystkich kon- figuracji występuje jednakowy trend zmian – wyższym war- tościom twardości odpowiadają niższe wartości energii EBN. Otrzymane współczynniki determinacji korelacji liniowej R2 pomiędzy twardością i wartościami średnimi energii EBN miesz- czą się w przedziale od 0,898 do 0,999. Podobne, pozytywne wyniki otrzymano przy analizie wartości skutecznej URMS. N i=1 0 -20 -40 -60 -80 -100 2 1,6 1,2 0,8 0,4 0 0,2 0,1 0 2 1 0 -85 -90 -95 -100 -40 -50 -60 U EN  [V ] Im [mA] -48 -50 -52 -54 -56 -58 880 840 800 760 720 H V 0, 5 PEAKPOS [mA] 880 840 800 760 720 -86 -88 -90 -92 -94 -96 PEAKPOS [mA] 880 840 800 760 720 -92 -94 -96 -98 -100 PEAKPOS [mA] H V 0, 5 H V 0, 5 73PRZEGLĄD  SPAWALNICTWA Vol. 89 11/2017 Rys.  9.  Zależność pomiędzy twardością i PEAKVAL dla ekstremum E2 – konf. C3 Fig. 9. Dependence between hardness and PEAKVAL values for extremum E2 – conf. C3 Rys.  10.  Zależność pomiędzy twardością i PEAKVAL dla ekstremum E3 – konf. C3 Fig. 10. Dependence between hardness and PEAKVAL values for extremum E3 – conf. C3 Rys.  8.  Przykładowa zależność pomiędzy twardością i PEAKVAL dla ekstremum E1 Fig. 8. Exemplary dependence between hard- ness and PEAKVAL values for extremum E1 Rys.  12.  Zależność pomiędzy twardością i EBN – konf. C2 Fig.  12. Dependence between hardness and EBN values – conf. C2 Rys.  13.  Zależność pomiędzy twardością i EBN – konf. C3 Fig. 13. Dependence between hardness and EBN values – conf. C3 Rys.  11.  Zależność pomiędzy twardością i EBN – konf. C1 Fig.  11. Dependence between hardness and EBN values – conf. C1 Analiza liczby zdarzeń szumu Barkhausena Zmierzony sygnał szumu Barkhausena stanowi zbiór impulsów napięciowych, wśród których możemy wyróżnić tzw. zdarzenia. Dla ustalonego progu napięciowego, wystą- pienie zdarzenia określane jest na podstawie 3 kolejnych próbek sygnału, których wartości są wyższe od tego progu i w którym próbka środkowa ma wartość wyższą od próbek skrajnych. Uśrednione rozkłady liczby zdarzeń w funkcji napię- cia dyskryminacji Ug opracowano dla trzech konfiguracji, a w każdej z nich dla dwóch kierunków magnesowania. Dla poszczególnych twardości można wyróżnić obszary, w których liczba zdarzeń różni się między sobą. Podobnie jak dla analizy częstotliwościowej analizowano rozkład współczynnika R2 korelacji liniowej pomiędzy liczbą zda- rzeń i twardością w funkcji napięcia dyskryminacji Ug – przy- kład na rysunku 18. Dla każdej konfiguracji wybrano napięcie dyskryminacji dla którego wartości R2 dla obu kierunków magnesowania świadczą o dobrej korelacji z twardością HV. Zależności te przedstawiono na rysunkach od 19 do 21. Analiza częstotliwościowa szumu Barkhausena  z wykorzystaniem FFT Zadaniem analizy częstotliwościowej jest badanie para- metrów i właściwości sygnału poprzez analizę jego struktury częstotliwościowej, umożliwioną dzięki zastosowaniu trans- formacji Fouriera, w wyniku której otrzymujemy tzw. widmo częstotliwościowe sygnału, prezentujące miarę korelacji (po- dobieństwa) sygnału analizowanego do poszczególnych ze- spolonych funkcji harmonicznych. Uśrednione periodogramy FFT otrzymano dla trzech kon- figuracji, a w każdej z nich dla dwóch kierunków magnesowa- nia. Dla poszczególnych twardości można wyróżnić obszary, w których amplitudy składowych częstotliwościowych różnią się między sobą. Opracowano algorytm wyznaczania współ- czynnika determinacji liniowej R2 pomiędzy wartościami am- plitudy dla poszczególnych składowych częstotliwościowych periodogramów odpowiadających trzem wartościom twardo- ści a twardością. W wyniku tego otrzymano rozkład wartości R2 w funkcji częstotliwości. Umożliwia to wybór częstotliwo- ści bądź zakresu częstotliwości, w których korelacja linio- wa pomiędzy wartościami amplitudy widma a twardością jest najlepsza – przykład na rysunku 14. Dla każdej konfiguracji wybrano jedną częstotliwość, dla której wartości R2 dla obu kierunków magnesowania świad- czą o dobrej korelacji z twardością HV. Zależności te przed- stawiono na rysunkach od 15 do 17. 0,2 0,450,4 0,40,50,6 0,550,8 0,80,61 0,651,2 1,2 1,6 20,71,4 880 840 800 760 720 880 840 800 760 720 880 840 800 760 720 H V 0, 5 H V 0, 5 H V 0, 5 PEAKVAL [V] PEAKVAL [V] PEAKVAL [V] 80 90 100 150 0170 40110 160 20180 60120 190 80130 200 100210 880 840 800 760 720 880 840 800 760 720 880 840 800 760 720 H V 0, 5 H V 0, 5 H V 0, 5 EBN EBN EBN 74 PRZEGLĄD  SPAWALNICTWA Vol. 89 11/2017 Rys.  16.  Zależność pomiędzy twardością i amplitudą składowej 20,4 kHz – konf. C2 Fig. 16. Dependence between hardness and amplitude of 20.4 kHz frequency compo- nent – conf C2 Rys.  15.  Zależność pomiędzy twardością i amplitudą składowej 33,1 kHz – konf. C1 Fig. 15. Dependence between hardness and amplitude of 33.1 kHz frequency compo- nent – conf. C1 Rys.  17.  Zależność pomiędzy twardością i amplitudą składowej 4,2 kHz – konf. C3 Fig. 17. Dependence between hardness and amplitude of 4.2 kHz frequency component – conf C3 Rys.  18.  Przykładowe wyniki analizy liczby zdarzeń: a) uśrednione rozkłady całkowitej liczby zdarzeń dla poszczególnych twardości, b) rozkład współczynnika R2 Fig. 18. Exemplary results of number of events analysis: a) averaged number of total events distributions for individual hardness values, b) R2 coefficient distribution Rys. 14. Wyniki analizy częstotliwościowej: a) uśrednione periodogramy FFT dla poszczególnych twardości, b) rozkład współczynnika R2 Fig. 14. Exemplary results of frequency analysis: a) averaged FFT periodograms for individual hardness values, b) R2 coefficient distribution 00 1010 55 1515 2020 10 8 6 4 2 0 1 0,8 0,6 0,4 0,2 0 A FF T [ m V ] f [kHz] f [kHz] R2 1,2 3,4 21,4 3,8 61,3 3,6 41,5 41,6 4,2 8 880 840 800 760 720 H V 0, 5 AFFT [mV] AFFT [mV] AFFT [mV] 880 840 800 760 720 H V 0, 5 880 840 800 760 720 H V 0, 5 Ug [V] Ug [V] 3000 2000 1000 0 1 0,8 0,6 0,4 0,2 0 N oE TO T 0,6 0,60,8 0,81 11,2 1,21,4 1,41,6 1,6 R2 75PRZEGLĄD  SPAWALNICTWA Vol. 89 11/2017 Podsumowanie  W celu opracowania korelacji umożliwiających wyznaczenie twardości na podstawie pomiarów szumu Barkhasuena analizowano obwiednię szumu Barkhausena i jej punkty charakterystyczne, energię, RMS, periodogramy FFT oraz rozkła- dy liczby zdarzeń. Uzyskano interesujące korelacje diagnostyczne mogące posłużyć w rozwiązywaniu zagadnień odwrot- nych badań nieniszczących. Dobre i bardzo dobre korelacje uzyskano dla energii i wartości skutecznej szumu Barkhausena, całkowitej liczby zdarzeń dla wybranych wartości napięcia dyskryminacji. Umożliwia to wielowariantową analizę sygna- łu diagnostycznego a przez to zwiększa pewność wyznaczania twardości na podstawie pomiarów szumu Barkhasuena. W potencjalnych aplikacjach tych korelacji należy pamiętać o uwzględnieniu czynników wpływających na sygnał diagno- styczny, np. chropowatość i kształt powierzchni, wpływ krawędzi elementu. Przedstawione.w.artykule.wyniki.zostały.uzyskane.w.badaniach.finansowanych. przez.Narodowe.Centrum.Nauki.w.ramach.umowy.UMO-2014/15/B/ST8/00203. oraz.w.ramach.projektu.OPUS.o.numerze.rej ..2014/15/B/ST8/0436 . Literatura [1] Moorthy V., Shaw B.A., Evans J.T. : Evaluation of tempering induced chan- ges in the hardness profile of case-carburised EN36 steel using magnetic Barkhausen noise analysis, NDT&E International, 36(1), pp. 43-49, 2003. [2] Stupakov O., Perevertov O., Tomas I., Skrbek B.: Evaluation of surface de- carburization depth by magnetic Barkhausen noise technique, Journal of Magnetism and Magnetic Materials, 323(12), pp. 1692-1697, 2011. [3] Vourna P. , Ktena A., Tsakiridis P. E., Hristoforou E.: A novel approach of ac- curately evaluating residual stress and microstructure of welded electrical steels, NDT&E Int., Vol. 71, pp. 33-42, 2015. [4] Roskosz M., Fryczowski K.: Wykorzystanie szumu Barkhausena jako me- tody pomiaru twardości materiału ferromagnetycznego, Przegląd Spawal- nictwa, Vol. 88, No 10, 2016, s. 143-145. [5] Kukla D., Piotrowski L., Szwed M.: Ocena twardości warstw nawęglanych i hartowanych indukcyjnie na stali AMS 6414 z zastosowaniem magne- tycznych technik nieniszczących, Przegląd Spawalnictwa, Vol. 88, No 10, 2016, s. 46-50. Rys.  20.  Zależność pomiędzy twardością i całkowitą liczbą zdarzeń powyżej progu dyskryminacji Ug=0,97 V – konf. C2 Fig. 20. Dependence between hardness and total numer of events – treshold voltage Ug=0.97V – conf. C2 Rys.  19.  Zależność pomiędzy twardością i całkowitą liczbą zdarzeń powyżej progu dyskryminacji Ug=1 V – konf. C1 Fig. 19. Dependence between hardness and total numer of events – treshold voltage Ug=1V – conf. C1 Rys.  21.  Zależność pomiędzy twardością i całkowitą liczbą zdarzeń powyżej progu dyskryminacji Ug=0,2 V – konf. C3 Fig. 21. Dependence between hardness and total numer of events – treshold voltage Ug=0.2V – conf. C3 400 520 80560 120600 160640 200680 720600 800 1000 NoETOT NoETOT NoETOT 880 840 800 760 720 880 840 800 760 720 880 840 800 760 720 H V 0, 5 H V 0, 5 H V 0, 5