Elektromagnetno polje v bližini napake v materialu UDK: 621.191:620.179.6 ASM/SLA: S13h, X8k Božidar Brudar UVOD Osnovni princip odkrivanja napak s pomočjo elektromagnetnega polja je v tem, da preiskovani predmet namagnetimo in potem merimo, kako se magnetno polje v njem porazdeljuje1. Na porazdelitev polja namreč pomembno vplivajo tudi napake v materialu. Vpliv radialnih razpok v kovinskih palicah in ceveh na izmenično magnetno polje je določil dr. F. Forster že leta 1954 s poskusi na modelih2. Rezultate njegovih eksperimentov so več kot 20 let navajali v najrazličnejših člankih s tega področja3. Tudi v jeseniški železarni pregledujemo površino kovinskih palic s tem, da merimo spremembe v impedanci tuljave, skozi katero preiskovane palice potujejo (Defectomat 2.189, izdelan v Inštitutu dr. Forster). Praksa pa je pokazala1, da pri tem delu lahko nekatere napake zgrešimo, če se faze signalov zelo močno razlikujejo od faze, ki jo je povzročila napaka, ki jo ravno še dopuščamo. To je še posebno problematično pri feromagnetnih materialih, kjer skušamo z močnim dodatnim enosmernim magnetnim poljem doseči magnetno nasičenje. Meritve in teoretični izračuni so pokazali, da nikdar ne vemo, s kolikšno vrednostjo relativne diferencialne permeabilnosti imamo opra- Slika 1 Splošna oblika površinske razpoke v palici Fig. 1 General čase of a surface crack in a bar. viti. Zato je skoraj nemogoče napovedovati, kolikšna je fazna razlika med signali, ki ustrezajo zelo tankim radialnim razpokam in napakam drugih oblik. Vsi ti problemi so nas silili, da smo sami začeli natančneje raziskovati princip delovanja aparatov, kjer ugotavljamo defekte s pomočjo elektromagnetnih polj. Takšna raziskava je koristna predvsem za praktično delo, saj omogoča ne le, da kritično preverimo, kaj aparat zmore registrirati, ampak tudi česa ne more. Tega podatka navadno v trgovskih prospektih ni. Pri svojem delu nismo eksperimentirali z živo-srebrnimi modeli, ampak smo simulirali napake s pomočjo matematičnih modelov. Poiskali smo splošno metodo, po kateri lahko simuliramo različne napake v palicah ali ceveh in pridemo do podobnih rezultatov, kot jih je izmeril dr. F. Forster. Študirali pa smo tudi vpliv enosmernega magnetnega polja na možnost odkrivanja napak v feromagnetnem materialu z znanimi magnetnimi karakteristikami. V prvem delu tega sestavka so opisane izhodne enačbe za izračunavanje polja, v nadaljevanju pa bo opisan postopeik za izračunavanje polja v palicah z radialno razpoko in rezultati pri simuliranju podpovršinskih in površinskih razpok v ceveh. SEZNAM UPORABLJENIH SIMBOLOV a polmer palice a, notranji polmer cevi A, A, A amplituda vektorskega potenciala B, B, B amplituda gostote magnetnega pretoka b oddaljenost od osi, kjer je vpliv napake zanemarljiv E električna poljska jakost f, F frekvenca H, H, H magnetna poljska jakost i imaginarna enota I jakost električnega toka v tuljavi j. j gostota toka L« induktivnost prazne tuljave n število ovojev na enoto dolžine tuljave P magnetni pretok r, r polmer t čas V vektorski potencial W amplituda vektorskega potenciala x, y, z pravokotne koordinate Z, Z impedanca (t električna prevodnost Ho permeabilnost praznega prostora (i relativna permeabilnost = 1 dB relativna diferencialna Vo dH permeabilnost to = 2-rtf koordinatni sistem. Enačbi (4) in (5) nadomestimo s sistemom enačb za posamezne komponente Ar, A,, A,, oziroma Ar*, A9*, A,*: d2 Ar J_ dA, 1 _Ar d2Ar dr2 + r ' dr +r2 ' dtp2 r2 dz2 7 d A --Z. . -f F . Ar* = 0, r2 d tp d2A„ 1 dA 1 d2 A, A d2 A H—' *--r-: • —-—r ——— H—— T (6) dr2 dr 2 d Ar dz2 d + F.A/ = 0, (7) (P d2Az + 1dA1+J_i^Ai+d^+F ^ = : dr r2 dtp2 dz2 dr2 (8) d2Ar* + J. + i d2Ar* Ar* | d2Ar* dr2 koordinata v smeri ep II. MATEMATIČNA FORMULACIJA A. Osnovne enačbe Magnetno polje v materialu in zunaj njega lahko opišemo na več načinov. Največkrat uporabljamo izhodne enačbe za vektorski potencial, skalami potencial ali enačbe za gostoto magnetnega pretoka. Izbira je odvisna od narave problema. 1. Vektorski potencial Izhajamo iz enačbe: dV d2 V V2V — tloJAlT—--HlloEEo —r = W dt dt2 Predpostavljamo, da gre za sinusno nihanje polja: V = W . eiwt. Namesto enačbe (1) zapišemo: V2W — iwniK)0-W + HPOEEO"2W = 0. (2) Pri neporušnih preiskavah delamo navadno z razmeroma nizkimi frekvencami, tako da tretji člen v enačbi (2) lahko zanemarimo. Iščemo rešitev enačbe: . «r n in zapišemo enačbo (3) za realno in imaginarno komponento: y2 A + w ^ ^ A* = 0, (4) V2A* —wnyo2V I 2 dz2 r2 d V d

neka izbrana vrednost vektorskega potenciala in F brezdimenzijska frekvenca. - Gostoto magnetnega pretoka pa izračunamo iz relacije: B = rotW V cilindričnih koordinatah to zapišemo za realno komponento na primer takole: 1 dK dA_ B„ = B„ = r dtp d z dAr_ d K dz dr r dr, r d (20) (21) 3. Robni pogoji S sistemom diferencialnih enačb (6) do (11) opisujemo vektorski potencial V, oziroma gostoto magnetnega pretoka B v snovi, kjer je električna prevodnost končna. V praznem prostoru veljajo podobne enačbe. Upoštevati je treba le, da je v zraku na primer Problem je rotacijsko simetričen. Napaka take vrste bi bila na primer nenadna sprememba v premeru palice in bo natančneje opisana v nadaljevanju. 2.1. Pogoji za vektorski potencial Namesto enačb (25) do (30) zdaj pišemo: Av _ . _ h L> dr r d A* (33) dr dAr + A * 12- =0, r dAz (43) = 0. d A* dz dr = 0, dz dBr dBz dz dr = 0, (45) — • — (r B ) = 0. r dr V oddaljenosti r = b je Bz = Ba edina komponenta gostote magnetnega pretoka, ki je različna od nič. Zato je tam Bz = 1 in Bz* = 0. Enačba (45) povezuje komponenti Br in Bz. Komponenta B,, pa z Br in Bz ni povezana. Zato bi v tem primeru morali reševati enačbi za Br in Bz. Na mejnih ploskvah r = konst., oziroma z = konst. velja robni pogoj (45), oziroma: V dz dr Jout Velja pa tudi izraz (23), oziroma: dr r dz (46) Podobno bi lahko zapisali tudi za imaginarne komponente Br* in Bz*. 3. Poseben primer napake, kjer je — =0. dz Napaka take oblike se enakomerno razteza vzdolž palice. 3.1. Pogoji za vektorski potencial Namesto enačb (25) do (30) dobimo: d\ + dr 1 dAr _ r r d d

a. 5.2. Pogoji za gostoto magnetnega pretoka V področju a ^ r ^ b veljajo enaki pogoji kot pri palici. Tudi v področju r ^ velja Bz = konst. Robni pogoj pri r = a, pa lahko poiščemo takole: Na notranjem robu prereza velja: rot E = — ^ (53) d t Enačbo (53) integriramo na obeh straneh po notranjem prerezu cevi: ^rot (B + i B*) d S = [ p0 p i no ueAoii AAHHe npyTKa. YpaBiieiiHe aasi rvcTOTti MarHHTHoro noAa HaAo bu^hcahti, toalko aasi oahoh KOMnOHeHTM. OniicaHHtie cnocoGu, pa3BHTHe Koxopux BbinoAHeHO b HayHHO-HCCAeAOBaTeABHOM HHCTHTyTe MeTaAAyprHHecKoro 3aB0Aa >KeAe3apna EceHHijc 6 u ah nphmehehbi TaKJKe ycnenmo Ha cnenHaAbHbix npn-Mepax (AByxMepHbie H TpexMepnMe MaTeMaTHiecKHe MOAeAH no- BepXHOCTHiIX h nOAIJOBepXHOCTHbIX nOpOKOB). B o6Aacra kohtpoah 6e3 pa3pymeHHJi c nyAbCHpyioniHM h MHoroiacTHUM SAeKTpoMarHHTHtiM noAeM npeACTaBAaeT onHcaHHe MaTeMaTH^ecKHX mctoaob coBceM hobyk> 6a3y aas 3KcnepHMeHTaAJb-hbix paSoT. moaeahpobahhe ycaobhh npn nomoiuh cienhka ro-pa3AO AeuieBAe h 6bicTpee ot npaKTmeeKoro 9KcnepH.vieHTHpoBaHHH. n. Eeprep/Ile.Vie