IV Conferencia Panamericana de END Buenos Aires Octubre 2007

Anlisis de las corrientes inducidas producidas por bobinas rectangulares planas

Resumen Se presentan los resultados tericos de los clculos de las corrientes inducidas por bobinas planas rectangulares sobre dos tipos de conductores planos no ferromagnticos: uno semi-infinito y otro de espesor finito, en condiciones del ensayo no destructivo por corrientes inducidas. Las expresiones que permiten los clculos de las corrientes forman parte de los resultados de este trabajo y son del tipo de transformada doble de Fourier. Se encontr que las mayores intensidades de las corrientes inducidas se presentan justo debajo de cada bobina, debido al pequeo lift-off que se tom para el clculo. Los mximos de esas intensidades responden a los tipos de arrollamientos que forman las bobinas, presentando una de las mismas dos mximos. Mediante el anlisis de la atenuacin de las corrientes al penetrar el material conductor, se pudo comprobar que se sigue una forma exponencial como la que propone el anlisis tradicional del efecto pelicular.

1. Introduccin
El uso de bobinas planas en ensayos de corrientes inducidas se muestra promisorio debido a sus ventajas para inspeccionar superficies planas: alta sensibilidad sobre este tipo de superficies, reduccin de las seales de acercamiento a los bordes en inspecciones para detectar fallas cercanas a esos bordes en las piezas inspeccionadas. Otras ventajas son su fcil adaptacin a superficies curvas debido a su flexibilidad o la reproducibilidad en cuanto a la fabricacin de las mismas si se hace mediante la tcnica de circuitos impresos [1]. Es por esto que, desde hace ya varios aos, nuestro grupo de trabajo le dedica especial atencin a este tipo de bobinas, estudiando la implementacin de las mismas en diferentes situaciones y analizando su respuesta desde el punto de vista experimental y terico [1-6]. En este camino es que resulta de especial inters conocer los campos electromagnticos (EM) producidos por este tipo de bobinas, los cuales nos permitan modelar las respuestas de las mismas; varios trabajos se pueden encontrar en este sentido [3, 4, 7, 8].

Las expresiones de los campos EM producidos por las bobinas que se estudian aqu se presentaron en trabajos anteriores [3, 4]. En ellos no fueron incluidos los clculos de las corrientes inducidas en el material, sino que los mismos estn centrados en el anlisis de los diagramas de impedancia. En este trabajo extendemos los resultados de [3, 4] para concluir con el clculo y el anlisis de las corrientes inducidas. En primer lugar se introduce una breve descripcin de las bobinas en estudio, para pasar luego a desarrollar las expresiones de las densidades de corrientes inducidas que se utilizan para el clculo, y por ltimo se presentan los resultados y el anlisis de los mismos.

2. Bobinas utilizadas
Las bobinas para las cuales se hicieron los clculos son las mismas que se presentaron en trabajos anteriores [1, 3, 4], no obstante lo cual haremos aqu una breve descripcin de las mismas para no perder continuidad en el desarrollo del trabajo. Se trata de bobinas espiraladas planas de seccin rectangular finita. Las bobinas fueron construidas con la tcnica de fabricacin de circuitos impresos y constan de caminos de cobre depositados sobre un substrato aislante de fibra de vidrio. La bobina identificada como bobina 1, est formada por dos arrollamientos de 6 vueltas cada uno, colocados en ambas caras del substrato aislante (constituyendo una bobina de 2 capas). La bobina que se llama bobina 2 es similar a la 1, pero formada por espirales de 10 vueltas. La otra bobina, que se identifica como bobina 3, consta de dos arrollamientos de 7 vueltas cada uno, pero que corren paralelos en la misma cara del substrato aislante (bobina de 1 capa). Ambos tipos de bobinas se muestran en la Figura 1. El ancho de los caminos de cobre y la separacin entre los mismos es de 0.2 mm y su espesor es aproximadamente de 0.065 mm. En todos los casos las bobinas estn recubiertas por una pelcula antisoldante de aproximadamente 0.02 mm. La forma y dimensiones reales de estas bobinas se tuvieron en cuenta para el clculo, y se aproximaron las espirales rectangulares por espiras concntricas.

Figura 1. Bobinas utilizadas. La bobina 1 es similar a la 2, pero de 6 vueltas

3. Densidad de corriente inducida

IV Conferencia Panamericana de END

Buenos Aires Octubre 2007

J = - j A

(1)

Figura 2. Problemas de inters: (a) conductor semi-infinito, (b) conductor finito. En la Figura 2 vemos las dos configuraciones para las cuales nos interesa conocer las densidades de corriente dentro del conductor a inspeccionar: un conductor infinito, Figura 2 (a) y un conductor finito, Figura 2 (b). En esta figura Ra, Rb y Rc son las regiones en las que se puede dividir el espacio para resolver el problema analticamente, segn corresponda. Como se menciona en [3, 4, 8] consideramos el problema del clculo de los campos EM producidos por una distribucin de corriente colocada sobre un conductor plano. Este conductor se supone istropo, homogneo y no ferromagntico, con conductividad elctrica y permeabilidad magntica = or, siendo r la permeabilidad relativa. La distribucin es excitada por una corriente de forma sinusoidal J(x,t) = J0(x)ej t, con = 2f. Para la resolucin se trabaja dentro de la aproximacin cuasiestacionara [9], aproximacin vlida en el rango de frecuencias utilizadas en los ensayos de corrientes inducidas. En [3] y [4] se establece que trabajando en el gauge de Coulomb ( A = 0 ) el potencial vector se puede suponer derivado de un potencial vectorial de segundo orden W [10]:
A = W con

(2)

Wa,b,c = W1,a,b,c + W2,a,b,c . z z

(3)

z Aqu es el versor en la direccin y sentido del eje z positivo, W1 y W2 son los potenciales escalares de los cuales depende W y a, b y c denotan las regiones en las que se divide el espacio para la resolucin del problema. Para la configuracin de la Figura 2 (a) existen slo las regiones a y b.
En [8] se obtienen las expresiones integrales para los potenciales escalares W1 y W2 para el caso de un conductor infinito, considerando una distribucin sinusoidal cualquiera J y se hace el clculo especfico para una distribucin formada por una sola espira rectangular. Ese trabajo es el punto de partida para [3] en el cual se resuelve el problema para bobinas rectangulares de cualquier nmero de vueltas y capas. En [4] se resuelve el

IV Conferencia Panamericana de END

Buenos Aires Octubre 2007

problema general de una distribucin sinusoidal sobre una placa conductora de espesor finito d (Figura 2 (b)) y se especifican esos resultados tambin para los casos de bobinas rectangulares de cualquier nmero de vueltas y capas; que son las de inters tambin en este trabajo. En todos los trabajos citados en el pargrafo anterior se parte de las ecuaciones que cumple A en cada regin, que segn las aproximaciones tenidas en cuenta toman la forma:
2 Aa = 0 en la regin R a ,

(4)

2 Ab 2 Ab = 0 donde 2 = j y
2 Ac = 0

en la regin R b ,

(5)

en la regin R c .

(6)

En las regiones sin conductor la densidad de flujo magntico slo depende del potencial escalar W1a,c de la siguiente forma, [3, 8]:
B a ,c = ( W1a,c ). z

(7)

En la regin Ra el potencial vector est formado por las contribuciones del potencial primario Aap, producido por la corriente excitadora y del potencial secundario Aas producido por las corrientes inducidas en el material. Por lo tanto W1a est formado por las mismas contribuciones: W1a = W1ap + W1as; cumpliendo, cada uno de estos potenciales escalares, las siguientes ecuaciones:
2W1ap = 0
2W1as = 0.

(8) (9)

Para la otra regin sin conductor, Rc, el potencial es producido exclusivamente por las corrientes inducidas en el material conductor y la ecuacin para el potencial escalar W1c que se debe resolver es
2W1c = 0.

(10)

Por ltimo en la regin Rb, aplicando las condiciones de continuidad de los campos B y H en las interfases (z = 0 para los casos de las Figuras 2 (a) y (b) y z = -d para el caso de la Figura 2 (b) solamente) y debido a que no hay fuentes dentro del conductor se puede concluir que el potencial escalar W2b debe ser cero dentro del conductor para ambos casos, referirse a [3] y [4] para ms detalles. Por lo tanto el potencial escalar que

IV Conferencia Panamericana de END

Buenos Aires Octubre 2007

interviene en el clculo de las corrientes inducidas es W1b, y la ecuacin correspondiente es: 2W1b 2W1b = 0. (11)

La ecuaciones (8), (9), (10), y (11) determinan las soluciones de los problemas planteados. Estas ecuaciones tienen soluciones generales que se pueden expresar mediante transformadas dobles de Fourier. La resolucin del problema comienza calculando W1ap para alguna distribucin de corriente conocida, [8]. En los casos que nos preocupan en este trabajo, las distribuciones conocidas para las cuales se calcula W1ap corresponden a las de las bobinas en estudio en ausencia de la placa conductora, estos resultados son los presentados en [3] y [4]. Llamando X1ap(, , z), X1as(, , z), X1b(, , z) y X1c(, , z) a las transformadas bidimensionales de Fourier de W1ap(x, y, z), W1as(x, y, z), W1b(x, y, z) y W1c(x, y, z) respectivamente, las soluciones para las ecuaciones (8) a (11) son:

X 1ap ( , , z) = K1ap ( , )e kz
X 1as ( , , z) = K1as ( , )e kz

(12) (13)

ambas son vlidas tanto para la configuracin de conductor infinito como finito;

X 1b ( , ,z) = K1b ( , ) e z
slo para el conductor infinito y
X 1c ( , , z) = K1c ( , )ekz

(14)

(15) (16)

X 1b ( , ,z ) = A( , ) ez + B ( , ) ez

para conductor finito, con k = 2 + 2 y 2 = k2 + 2. Como se mencion ms arriba, las frmulas para K1ap se conocen, ya que se tienen las soluciones de W1ap para las bobinas estudiadas en el espacio libre, las cuales son las transformadas inversas de Fourier de X1ap:
j x+y ) W1ap = F 1 X 1ap = K1ap ( , ) ekz e ( d d
+ +

(17)

Aplicando las condiciones de contorno para la continuidad de la componente normal de B y para la componente tangencial de H en z = 0 y en z = -d para las ecuaciones de los potenciales escalares, segn corresponda para el caso de conductor infinito o finito, se obtienen las expresiones buscadas para K1as, A, B y K1c.

IV Conferencia Panamericana de END

Buenos Aires Octubre 2007

Para el caso de conductor infinito, segn [8] se tiene:

r k r k +
2r k r k +

(18)

K1ap = K1ap

(19)

La resolucin de la configuracin con conductor de espesor finito da las siguientes expresiones, [4]:
K1as = K1ap

(1 e ) g h ,
2 d

g 2e2 d h 2

(20)

K1c = K1ap e(

k )d

g 2 h2 , g 2e2 d h 2 y

(21)

A = K1ap

h g 2 h2 2 g 2e2 d h 2

(22)

B = K1ap con g = 1 rk/ y h = 1 + rk/.

g 2 d g 2 h 2 e 2 g 2e2 d h 2

(23)

Volviendo al potencial escalar W1b, tenemos entonces las siguientes expresiones:

W1b, inf = K1ap ( , )

+ +

2r k z j(x+y) e e d d r k +

(24)

W1b, fin = A( , )e z + B ( , )e z e j(x+y)d d , para un conductor infinito y finito respectivamente.

(25)

Las densidades de corriente se obtienen combinando (24) y (25) con (1) y (2), con lo que se llega a:
J inf = K1ap

+ +

+ +

2r k z e ( x + y ) e j(x+y)d d r k +

(26)

J fin = Ae z + Be z ( x + y ) e j(x+y)d d .

(27)

IV Conferencia Panamericana de END

Buenos Aires Octubre 2007

Resolviendo la antitransformada de Fourier por algn mtodo conveniente obtenemos las distribuciones de corriente buscadas, y estos son los resultados presentados en este trabajo. Antes de seguir debemos dedicarnos brevemente al coeficiente K1ap. Este coeficiente es de suma importancia porque contiene la geometra de la bobina a utilizar en cada clculo, que para nosotros son las bobinas descriptas en la seccin 2. K1ap para los dos tipos de bobinas que se estudian aqu fueron calculados en [3].

4. Resultados
La resolucin numrica de las ecuaciones (26) y (27) nos permite computar y conocer las densidades de corrientes inducidas. Para todos los clculos, que se presentan a continuacin, se consider al material conductor como aluminio, tomando para el mismo: = 35.4 x 106 S/m y r = 1, y a una frecuencia de ensayo de 100 kHz. Estos parmetros arrojan para la profundidad de penetracin estndar , el valor de 2.67 x 10-4 m. Aqu es (vase ec. (29)).

(28)

Con el fin de dar algunos resultados numricos, se presentan las densidades de corrientes en z = 0 para las bobinas 1 y 3 (Figuras (3) y (5)) en el caso de conductor infinito; y para la bobina 1 sobre un conductor finito de espesor d = 1 = 2.67 x 10-4 m (Figura (7)). Para cada tipo de bobina, las formas que toman los caminos de las corrientes inducidas son muy parecidas ya sea para conductor finito o infinito, es por esto que hemos decidido mostrar en este trabajo slo una distribucin para cada una. Para todos los clculos la separacin en la direccin x (o y) entre dos puntos consecutivos en los que se calcula la corriente se tom igual al ancho de los caminos de cobre de las bobinas, esto es 0.2 mm. En todas estas figuras podemos ver que la densidad de corriente toma sus mayores valores en la regin del material conductor que est justo debajo de toda la seccin de las bobinas. Esto se debe a que, para el clculo, la separacin que se tom entre las bobinas y el material conductor en muy pequea: corresponde a 0.02 mm, que es aproximadamente el espesor de la mscara antisoldante que protege a las bobinas. En las Figuras (3), (5) y (7) se muestran como rectngulos de lnea llena los lmites interno y externo de las bobinas. Las lneas punteadas indican, de manera aproximada, la posicin de los mximos de J. Con el fin de mostrar mejor cmo vara la intensidad de las corrientes inducidas, se representa en las Figuras (4) y (6) la componente y de J en funcin de x, para los valores correspondientes a J a lo largo del eje x positivo. Como podemos apreciar de la Figura (6) (tambin marcado en la Figura (5)) la bobina 3 presenta dos mximos de la intensidad de J. Esta particularidad puede deberse al diseo de la bobina 3: dos espirales rectangulares paralelos sobre la misma cara.

IV Conferencia Panamericana de END

Buenos Aires Octubre 2007

Figura 3. Corrientes inducidas por la bobina 1 (2 capas de 6 vueltas cada una), en z = 0 sobre conductor semi-infinito.

Figura 4. Componente y de la distribucin de corriente inducida por la bobina 1 sobre un conductor de espesor semi-infinito a lo largo del eje x positivo, en z = 0.

IV Conferencia Panamericana de END

Buenos Aires Octubre 2007

Figura 5. Corrientes inducidas por la bobina 3 (1 capa con dos bobinados paralelos de 7 vueltas cada uno), en z = 0 sobre conductor infinito.

Figura 6. Componente y de la distribucin de corriente inducida por la bobina 3 sobre un conductor de espesor semi-infinito a lo largo del eje x positivo, en z = 0. Se ven los dos mximos.

IV Conferencia Panamericana de END

Buenos Aires Octubre 2007

Figura 7. Corrientes inducidas por la bobina 2 (2 capas de 10 vueltas cada una), en un conductor finito de espesor d = 1 = 2.67 x 10-4 m, para z = 0.
Por otro lado se analiz la atenuacin de las corrientes con la penetracin en el material conductor semi-infinito (Figuras (8) y (9)), comparando los resultados obtenidos con la relacin clsica del efecto pelicular [9], la cual se puede escribir de la siguiente manera [11]:

Jz = e cos (t - ) J0

(29)

con = z , siendo z la profundidad de penetracin de la seal en el conductor y Jz

simboliza la intensidad de J a la profundidad z. Para realizar este anlisis se calcul la densidad de corriente en funcin de la profundidad, correspondiente al punto en el plano xy en el cual se alcanza el mayor valor de J para z = 0. Si consideramos la atenuacin segn la relacin clsica debemos esperar una forma exponencial como el primer factor de la ecuacin (29), esto se grafica en lnea llena en la figura (8) como funcin de la profundidad (dada en unidades de ). Sin embargo, si se grafica el mdulo de J calculado, en funcin tambin de la profundidad, obtenemos lo que se muestra en la figura (9) con puntos. La diferencia, mencionada en el prrafo anterior, entre el mdulo de J calculado en funcin de la profundidad y el factor e se debe a que la corriente calculada por

IV Conferencia Panamericana de END

Buenos Aires Octubre 2007

10

cualquiera de las ecuaciones (26) y (27) es un vector complejo cuya parte real es la que se interpreta como la densidad de corriente J y es esa misma componente la que aparece en (29). Para interpretar los resultados obtenidos con la expresin (29) llamamos a la fase del vector complejo J calculado (utilizamos para no confundirla con la definicin de = z ). Tanto la componente x como la componente y de la parte real de J calculada tiene la forma de (29); pero considerando en lugar de la fase , esto nos dice que deberemos dividir a cada valor calculado de J por el correspondiente cos() para obtener el factor exponencial decreciente e- que se ajuste a e-, esto es lo que se representa con puntos en la figura (8), como se ve se logra el ajuste adecuado. De la misma manera el J calculado en funcin de la penetracin es tambin una expresin de la forma (29), por lo que la funcin que ajusta a los puntos que lo representan ser: e cos ( ) , esto se muestra en la figura (9).

Figura 8. |J|/cos() (puntos) y e- (lnea llena), en funcin de la profundidad de penetracin.

Figura 9. Mdulo de J (puntos) y e-cos() (lnea llena), en funcin de la profundidad de penetracin.

IV Conferencia Panamericana de END

Buenos Aires Octubre 2007

11

5. Conclusiones
El conocimiento de los campos EM producidos por bobinas rectangulares de seccin finita sobre conductores planos semi-infinitos o de espesor finito [3, 4], permiti obtener las expresiones analticas para el clculo de las corrientes inducidas en esas situaciones. La evaluacin de esas expresiones, que nos permiten el clculo de las corrientes, se realiz reproduciendo la geometra real de las bobinas que se utilizan en nuestro grupo de trabajo, as como el lift-off tpico de las mismas. Si bien para el clculo se fij una situacin particular: conductor de aluminio y frecuencia de inspeccin de 100 kHz, esto no rest generalidad a los resultados obtenidos. Para todas las bobinas analizadas, las mayores intensidades de las corrientes inducidas estn justo debajo del rea de las mismas, atenundose rpidamente al alejarse de ellas sobre el plano que las contiene. Las formas de las densidades de corrientes estn de acuerdo a lo que se esperaba segn la geometra de las bobinas, esto es una indicacin de que el clculo realizado fue correcto. Del anlisis de las intensidades de J en z = 0 se localizaron los mximos de las mismas en el plano xy para cada una de las bobinas: Las bobinas construidas con dos espirales paralelos en caras opuestas mostraron tener un solo mximo. Mientras que la bobina que consta de dos espirales en la misma cara del sustrato aislante presenta dos mximos para la intensidad de J. Como otro resultado del mismo clculo se estudi el efecto de la atenuacin de las corrientes con la penetracin en el material. Se concluye que se reproducen las caractersticas generales del efecto pelicular para un conductor plano, otro indicador de que es satisfactorio el clculo de las corrientes realizado para cada una de las situaciones planteadas. Este trabajo complementa otros que se vienen realizando con anterioridad [1, 3, 4] y que son de gran utilidad para nosotros y para todos aquellos que trabajen con bobinas planas rectangulares.

Agradecimientos
Los autores agradecen a Dr. Theodoros Theodoulidis, de University of West Macedonia, Grecia. Tambin agradecen a las instituciones CNEA y CONICET y la ayuda econmica de la agencia argentina ANPCYT a travs de los proyectos PICT 12113 y 14426.

Referencias
1. 2. J O Fava and M C Ruch, ' Design, Construction and Characterization of Ect Sensors with Rectangular Planar Coils', Insight, Vol 46, No 5, pp 268-274, 2004. M C Ruch, L S Lanzani, and J O Fava, 'Non Destructive Studies of Hydrogen in Zircaloy-4', in Technical Meeting on Behaviour of High Corrosion Resistance Zr Based Alloys, Buenos Aires: IAEA and CNEA, 24-28 October, 2005.

IV Conferencia Panamericana de END

Buenos Aires Octubre 2007

12

3.

4.

5.

6.

7.

8.

9. 10. 11.

J O Fava and M C Ruch, 'Calculation and Simulation of Impedance Diagrams of Planar Rectangular Spiral Coils for Eddy Current Testing', NDT & E International, Vol 39, No 5, pp 414-424, 2006. J O Fava and M C Ruch, 'Simulation of Eddy Currents Impedance Diagrams in the Inspection of a Conducting Plate with Planar Rectangular Coils', in NDT 2006 Conference & Exhibition, Stratford-upon-Avon: BINDT, 11-14 September, 2006. J O Fava and M C Ruch, 'Informe Tcnico in-13-E-Im 079/05: "Evaluacin Por Corrientes Inducidas Del Contenido De Hidrgeno Y De La Presencia De Ampollas En Tubos De Presin Candu Hidrurados"': ENDE-CNEA, 2005. G Domizzi, et al., 'Hydrogen Determination and Blister Detection in Zr-Alloys', in Technical Meeting on Behaviour of High Corrosion Resistance Zr Based Alloys, Buenos Aires: IAEA and CNEA, 24-28 October, 2005. N Tsopelas and N J Siakavellas, 'Performance of Circular and Square Coils in Electromagnetic-Thermal Non-Destructive Inspection', NDT & E International, Vol 40, No 1, pp 12-28, 2007. T P Theodoulidis and E E Kriezis, 'Impedance Evaluation of Rectangular Coils for Eddy Current Testing of Planar Media', NDT & E International, Vol 35, No 6, pp 407-414, 2002. H L Libby, 'Introduction to Electromagnetic Nondestructive Test Methods', New York: Wiley-Interscience, pp 128-129, 145, 1971. W R Smythe, 'Static and Dynamic Electricity', 3rd ed, New York: MacGraw-Hill Inc, pp 284-286, 1968. N Ida, 'Numerical Modeling for Electromagnetic Non-Destructive Evaluation', Engineering Non-Destructive Evaluation Series, London: Chapman & Hall, pp 82-84, 1995.

IV Conferencia Panamericana de END

Buenos Aires Octubre 2007

13