Ingeniería Mecánica E-ISSN: Instituto Superior Politécnico José Antonio Echeverría.


Save this PDF as:
 WORD  PNG  TXT  JPG

Tamaño: px
Comenzar la demostración a partir de la página:

Download "Ingeniería Mecánica E-ISSN: Instituto Superior Politécnico José Antonio Echeverría."

Transcripción

1 Ingeniería Mecánica E-ISSN: Instituto Superior Politécnico José Antonio Echeverría Cuba Carro-Hernández, Gustavo Eugenio; Cruz-Inclán, Carlos Manuel; Arteche-Díaz, Josué; Leyva-Fabelo, Antonio Simulación numérica de un nuevo sensor de anisotropía magnética Ingeniería Mecánica, vol. 16, núm., mayo-agosto, 013, pp Instituto Superior Politécnico José Antonio Echeverría Ciudad de La Habana, Cuba Disponible en: Cómo citar el artículo Número completo Más información del artículo Página de la revista en redalyc.org Sistema de Información Científica Red de Revistas Científicas de América Latina, el Caribe, España y Portugal Proyecto académico sin fines de lucro, desarrollado bajo la iniciativa de acceso abierto

2 Ingeniería Mecánica. Vol. 16. No., mayo-agosto, 013, p ISSN Artículo Original Simulación numérica de un nuevo sensor de anisotropía magnética Numeric Simulation of a new magnetic anisotropy sensor Gustavo Eugenio Carro Hernández, Carlos Manuel Cruz Inclán, Josué Arteche Díaz, Antonio Leyva Fabelo Centro de Aplicación Tecnológica y Desarrollo Nuclear. CEADEN. La Habana. Cuba. Correo electrónico: Recibido: de noviembre de 01 Aceptado: 5 de marzo de 013 Resumen El objetivo del trabajo es evaluar la aplicabilidad de la metodología empleada para la determinación de las tensiones residuales con el sensor de anisotropía magnética SAM0 a una nueva versión más avanzada SAM. Se realizó la caracterización de la función de respuesta del nuevo sensor utilizando la simulación numérica, observándose que la misma presenta una periodicidad con respecto a las posiciones angulares en que se rotó el mismo para un valor dado de anisotropía magnética. Se demostró la existencia de una dependencia lineal entre los valores de la FEM medidas con el nuevo sensor y los valores de anisotropía magnética de la muestra. Las FEM medidas en los tres circuitos del SAM son equivalentes a las de tres sensores del primer modelo orientados a 10 entre sí. Se demuestra la aplicabilidad de la función de respuesta del SAM0 para describir las lecturas de su nueva versión. Palabras claves: sensor de anisotropía magnética, tensiones residuales, método de elementos finitos. Abstract The aim of this study is to evaluate the applicability of the residual stresses determination methodology, previously used with magnetic anisotropy sensor SAM0, with the new more advanced sensor SAM. The characterization of the new sensor response function was carried out using the numeric simulation. A periodic behavior with respect to the rotated angular positions for a given value of magnetic anisotropy was observed. A lineal dependence between the EMF values measured with the new sensor and the sample magnetic anisotropy was demonstrated. The EMF measured in the three SAM circuits are equivalent to those obtained for three SAM0 sensors, oriented 10 between them. The applicability of the SAM0 response function to describe the measurements made with the new sensor was established. Key words: magnetic anisotropy sensor, residual stress, finite element method. Sitio web: 113

3 Introducción Simulación numérica de un nuevo sensor de anisotropía magnética Las Tensiones Residuales (TR) de primer género, o macrotensiones, se hallan presentes prácticamente en todo material manufacturado o de estructura ensamblada. En los procesos de soldadura, recargue, tratamientos térmicos, proyección térmica, recubrimientos galvánicos y termodifusivos, o como resultado de procesos donde la deformación plástica esté presente, es importante conocer el comportamiento de las TR que se encuentran presentes en las piezas a producir o recuperar [1]. Existen varios métodos utilizados para la determinación de las TR [], encontrándose entre los más comunes los basados en la medición de las propiedades mecánicas, que son métodos destructivos, y los métodos basados en la difracción de radiaciones, que aunque son no destructivos, sí requieren equipos de alta complejidad y elevados costos. En cambio, los métodos magnéticos (en especial los basados en la anisotropía de la permeabilidad magnética) son más sencillos, económicos, portátiles, además de ser no destructivos. Se conoce de la utilización de Sensores de Anisotropía Magnética (SAM) para la medición de TR desde hace más de dos décadas. Por ejemplo, en 1987 Yamada [3] describe la teoría en la que se basa un SAM, con el cual es posible la medición de TR de flexión sin necesidad de contacto entre la superficie de la muestra y el sensor. En 1998 Govindaraju [4] emplea por primera vez el Método de Elementos Finitos (MEF) [5] mediante el software ANSYS TM EMAG como herramienta para el diseño de un nuevo SAM. Este autor relaciona la frecuencia de excitación del sensor con el poder de penetración de las líneas del campo magnético que pasan a través del material de estudio con el fin de evaluar las TR que se encuentran a diferentes profundidades. En [6] (004) Sakai extiende la aplicación de sensores SAM a superficies cilíndricas realizando la medición de tensiones en líneas de tuberías. Liu en el 011 [7] realiza un análisis teórico exhaustivo de la señal de salida de un SAM llegando a la conclusión de que el voltaje de salida del sensor es proporcional a la diferencia de las tensiones principales y está relacionada con sus direcciones. Esta serie de trabajos muestra el interés internacional por el desarrollo de los SAM como herramienta para la medición de tensiones en diferentes situaciones prácticas. Debido a la importancia del seguimiento y control en la industria cubana de las TR se han desarrollado en el Centro de Aplicaciones Tecnológicas y Desarrollo Nuclear, CEADEN, de Cuba, diferentes versiones de SAM. La primera de ellas se denominó SAM0 [8], cuya función de respuesta toma en consideración que la frecuencia del campo magnético excitador tiende a cero; es decir, el campo se considera cuasi estático y que las secciones transversales de los circuitos de excitación y de lectura se desprecian en comparación con la distancia entre sus puntas. En el 006 se consolida la creación del equipo medidor de TR denominado TENSO y se realizan mejoras al sensor SAM0 sobre la base de utilizar materiales de mayor calidad Permalloy C [9]. Este sensor es empleado en el 011 para la evaluación de TR en barras cilíndricas de acero sometidas a rodadura en frío [10, 11]. La figura 1 muestra los principales componentes del sensor SAM0, que consta de un circuito magnético excitador, formado a su vez por una bobina de cobre enrollada a un núcleo de hierro, y un circuito magnético de lectura con las mismas características que el de excitación pero orientado a 90 con respecto al primero. El circuito de lectura es el encargado de medir las variaciones de la Fuerza Electromotriz (FEM) generadas por la anisotropía magnética (Δµ), inducidas a su vez por las tensiones mecánicas (Δσ) existentes en la muestra. Fig. 1. Representación esquemática de los componentes del sensor SAM0 Ingeniería Mecánica. Vol. 16. No., mayo-agosto, 013. ISSN

4 Gustavo Eugenio Carro-Hernández, Carlos Manuel Cruz-Inclán, Josué Arteche-Díaz, Antonio Leyva-Fabelo El procedimiento para el uso del SAM0 consiste en colocar las puntas sobre la superficie de un material ferromagnético y manualmente hacer rotar al sensor para colocarlo en, al menos, tres posiciones angulares (0,10,40 ), tomando la lectura de la FEM en cada uno de estos puntos (ε(θ)). Con ello se obtienen tres valores de ε(θ), necesarios para resolver el sistema de ecuaciones generado a partir de la función (1) que permite determinar la amplitud de la función de respuesta A (A=ε Δµ), el ángulo α que forman los ejes de las tensiones principales con los ejes de referencia x y y, y el coeficiente ε 0, que es una constante y no depende de las tensiones aplicadas. Los valores de anisotropía magnética Δµ determinados siguiendo este procedimiento van a ser proporcionales a las TR existentes en la muestra [6]. ) sin(( )) (1) ( 0 Con el fin de simplificar este procedimiento y de evitar la ocurrencia de errores en la medición, fue propuesto un nuevo sensor de anisotropía magnética denominado SAM. La idea original de este sensor se basa en colocar tres sensores SAM0 en las tres posiciones angulares antes mencionadas (0, 10, 40 ). Posteriormente con el objetivo de disminuir el tamaño del sensor, se integraron las tres bobinas de excitación a un solo circuito magnético con una única bobina y se conservaron los tres circuitos magnéticos de lectura independientes. En la figura se muestran la distribución de los diferentes componentes del SAM y las dos muestras en las cuales se posicionó el sensor en la simulación (caso (A) y caso (B)). Fig.. Representación esquemática del Sensor SAM y las muestras utilizadas en la simulación El objetivo de este trabajo es demostrar mediante la simulación numérica que la ecuación 1 para la determinación de los parámetros ε Δµ y α como parte de la metodología utilizada en el SAM0 continúa siendo aplicable al nuevo sensor SAM para los casos de dos situaciones prácticas diferentes: (A) para muestras donde los bordes se encuentran alejados de las puntas de excitación, y (B) cuando los bordes se encuentran cercanos a éstas, situación en que se manifiesta el denominado efecto de borde. También se pretende demostrar que el sensor SAM se comporta como tres sensores SAM0 orientados a 10 uno con respecto al otro. Materiales y Métodos Para la resolución de las complejas ecuaciones diferenciales que definen el comportamiento de los campos magnéticos, en las distintas configuraciones del sistema SAM+Muestra se utilizó el MEF, empleando como herramienta el sistema de códigos COMSOL Multiphysics versión 4.0. Se calculó el comportamiento de los campos magnéticos realizando un estudio dependiente del tiempo perteneciente al módulo AC/DC del software. Este proceder permite resolver la ecuación utilizando como aproximación un modelo cuasi-estático: 1 P e ( ˆ P M ) J () t Ingeniería Mecánica. Vol. 16. No., mayo-agosto, 013. ISSN

5 donde Simulación numérica de un nuevo sensor de anisotropía magnética 1 es el inverso del tensor de la permeabilidad magnética, P es el potencial magnético, M es la magnetización, σ es la conductividad eléctrica, t es el tiempo y J e es la densidad de corriente externa, que debido a sus características se puede representar como se muestra en la ecuación 3. J e ie ( )sin(ft) (3) S donde i e (0,1 A) es la corriente que pasa por los cables que conforman la bobina de excitación, S es la sección transversal de la bobina de excitación y f (700 Hz) es la frecuencia de J e. En la tabla 1 se muestran las propiedades de los materiales que conforman el modelo utilizado en la simulación numérica, donde µ r es la permeabilidad relativa y ε r la permitividad relativa. Tabla 1. Materiales utilizados en el modelo Materiales Componentes µ r σ (S/m) ε r Hierro Puntas de lectura y de excitación x Cobre Bobinas de lectura y de excitación x Material aislante Material que rodea al sensor 1 1x Las condiciones de frontera del modelo se establecieron considerando que toda la superficie externa de la muestra y del material aislante se encuentra aislada magnéticamente. Con esta opción se asume que la componente normal del vector densidad del flujo magnético ( B ) en dicha superficie es igual a cero ( n B 0 ). Esta condición es válida cuando se considera que la superficie externa se encuentra rodeada de aire, por lo que las líneas de flujo no escapan del modelo asumido. En la bobina excitadora se asume la presencia de una densidad de corriente externa, cuyas componentes ecuaciones (4) y (5). i sinft y S J e x x y J e y i S sinft x x y e J x y e J y se describen mediante las En el modelo se empleó el elemento tipo tetraédrico, donde se varió su tamaño en tres regiones, dependiendo de la exactitud necesaria para determinar la función de respuesta. En la tabla se observan los tamaños utilizados en el modelo y en la figura 3 se muestra un ejemplo del mallado realizado para el caso (A). (4) (5) Tabla. Tamaños de los Elementos Componentes Tamaño Máximo del Elemento Tamaño Mínimo del Elemento Bobinas 3,7 mm 0,666 mm Circuitos 5,55 mm 1,04 mm Resto 7,03 mm 1,58 mm Ingeniería Mecánica. Vol. 16. No., mayo-agosto, 013. ISSN

6 Gustavo Eugenio Carro-Hernández, Carlos Manuel Cruz-Inclán, Josué Arteche-Díaz, Antonio Leyva-Fabelo Fig. 3. Ejemplo de mallado utilizado en el caso (A) El tiempo t empleado en la simulación se dividió en intervalos t T / 4 s, donde T es el período de la oscilación de la corriente, desde un t 0 s hasta t 50 t s. 0 En el caso (A) se introdujeron variaciones en los valores principales µ i del tensor de la permeabilidad magnética ˆ en el plano xy de la muestra. Los valores de µ 1 se aumentaron manteniendo µ fijo, creando de 1 esta forma una anisotropía de la permeabilidad magnética. También para cada valor de Δµ se hizo variar α cada 15nπ/180 grados (n=0; 1; ; 3; 4), lo que es similar a considerar que el sensor se está rotando en estas posiciones angulares. Estos cambios se realizaron conforme al tensor 6: 1 cos ( ) sin ( ) ˆ xy sin( ) 0 1 xy sin( ) cos ( ) sin 0 ( ) (6) Para el caso (B), donde se consideran los bordes cercanos a las puntas de excitación, se aumentó la anisotropía magnética en la dirección axial y para cada valor de anisotropía magnética se rotó el sensor en las mismas posiciones angulares que el caso (A). Los cambios de anisotropía magnética introducidos en las muestras en una determinada dirección se considera que fueron inducidos mediante la aplicación de tensiones mecánicas que se encuentran en la misma dirección. Esto es debido a la dependencia lineal existente entre estos parámetros para pequeños valores de anisotropía magnética [9]. db Posteriormente se extrajeron los valores de FEM generada a partir de la componente que se encuentra dt en dirección normal a la sección transversal de los núcleos magnéticos de lectura. Para el cálculo de la FEM en cada bobina de lectura del sensor SAM se utilizó la ecuación 7. FEM N l k k db db x y sin( ) s cos( ) s (7) i1 dt i1 dt donde k es el número de nodos que se encuentran en dicha sección y N l es el número de vueltas del enrollado de las bobinas de lectura. Ingeniería Mecánica. Vol. 16. No., mayo-agosto, 013. ISSN

7 Simulación numérica de un nuevo sensor de anisotropía magnética Análisis de los resultados Como resultado de la simulación se obtuvo el comportamiento de la FEM con respecto a la posición angular para diferentes valores de anisotropía de la permeabilidad magnética en las tres bobinas de lectura para los dos casos estudiados. Estos valores se ajustaron por el método de Mínimos Cuadrados [1] a la función de respuesta utilizada en el SAM0 mostrada en la ecuación (1). En las figuras 4 (a), (b) y (c) se representan las lecturas de las bobinas 1, y 3 respectivamente del caso (A), mientras que en las figuras 5 (d), (e) y (f) se muestran las correspondientes al caso (B). (a) (d) (b) (e) (c) (f) Fig. 4. Comportamiento de la FEM de lectura con respecto a α para el caso (A) gráficas (a), (b) y (c) y con respecto a θ para el caso (B), gráficas (d), (e) y (f) Para conocer el grado de correlación entre la función de respuesta ε(θ) utilizada en el SAM0 y los valores de la FEM obtenidos en la simulación, se realizó una regresión lineal entre ellos para determinar el coeficiente de correlación lineal R en cada valor de anisotropía magnética. En la figura 5 se presentan algunos ejemplos para los dos casos estudiados, y en la tabla 3 se muestra el coeficiente R obtenido para todos los valores de anisotropías magnéticas simulados. Ingeniería Mecánica. Vol. 16. No., mayo-agosto, 013. ISSN

8 Gustavo Eugenio Carro-Hernández, Carlos Manuel Cruz-Inclán, Josué Arteche-Díaz, Antonio Leyva-Fabelo (a) (d) (b) (e) (c) (f) Fig.5. Ejemplos de regresión lineal entre ε(θ) y los valores de la FEM obtenidos en la simulación donde (a), (b) y (c) pertenecen a las bobinas 1, y 3 respectivamente del caso (A) (Δµ=8000) y (d), (e) y (f) al caso (B) (Δµ=5000). Tabla 3. Coeficiente de correlación lineal (R ) existente entre la función de respuesta ε(θ) y los valores de la FEM obtenidas en la simulación para los casos (A) y (B) Caso (A) Δμ Bobina 1 Bobina Bobina 3 Δμ Bobina 1 Bobina Bobina ,9948 0,9945 0, ,943 0,9488 0,8843 Caso (B) ,9976 0,9967 0, ,9681 0,9708 0, ,9987 0,9973 0, ,9768 0,9784 0, ,998 0,9974 0, ,9806 0,9817 0, ,973 0,9746 0,9116 En los resultados mostrados se observa una buena correlación entre los valores de la FEM obtenidos en la simulación numérica y la ecuación 1. Se puede notar que en el caso (B) los valores de R resultan ligeramente menores, lo cual puede estar relacionado con el empleo de un modelo geométrico más complejo, no simétrico (ver Fig. ), que provoca la no-homogeneidad en la FEM con relación al caso (A) (modelo Ingeniería Mecánica. Vol. 16. No., mayo-agosto, 013. ISSN

9 Simulación numérica de un nuevo sensor de anisotropía magnética geométrico simétrico), y por lo tanto una pequeña desviación de su comportamiento con respecto a la función de respuesta. Específicamente este fenómeno es atribuible al efecto de borde. El desfasaje que se observa entre las lecturas de cada bobina en la figura 4 se debe a la distinta posición angular relativa en que se encuentran cada una de ellas. A partir de estos resultados se determinó el grado de dependencia lineal entre los valores de la amplitud de la función de respuesta (A) y los valores de anisotropía magnética introducidos. En la figura 6 se muestran dichas dependencias pudiéndose apreciar la existencia de una buena correspondencia. (a) (d) (b) (e) (c) (f) Fig. 6. Dependencia lineal entre el parámetro A y Δµ, donde (a), (b) y (c) corresponden a las bobinas 1, y 3 del caso (A) y (d), (e) y (f) a las del caso (B) Como se puede apreciar, para cada caso estudiado, los resultados que se obtienen son similares independientemente de la bobina que se analice y están en correspondencia con la ecuación (1). Esto nos permite confirmar que cada circuito de lectura del sensor SAM se comporta como un sensor SAM0. La dependencia de la dirección de la anisotropía magnética con respecto a los ejes de referencia x y y con respecto a los valores de Δµ introducidos en la simulación es presentada en las figuras 7 (a) y (b). Ingeniería Mecánica. Vol. 16. No., mayo-agosto, 013. ISSN

10 Gustavo Eugenio Carro-Hernández, Carlos Manuel Cruz-Inclán, Josué Arteche-Díaz, Antonio Leyva-Fabelo (a) (b) Fig. 7. Dependencia de α en función de los cambios de anisotropía magnética para el caso (A) (a) y para el caso (B) (b) Se puede observar que α se mantiene constante, lo cual indica que este parámetro es independiente de los cambios de anisotropía magnética y sólo depende de la dirección en que ésta se encuentra. Conclusiones Mediante la simulación numérica utilizando el Método de Elemento Finitos se verificó la aplicabilidad de la función de respuesta empleada en el sensor de anisotropía magnética SAM0 al nuevo sensor SAM. La relación lineal obtenida entre los parámetro A y Δµ en cada bobina de lectura para cada caso simulado, demuestra el carácter constante del parámetro ε, lo que permite establecer la existencia de una relación igualmente lineal entre la amplitud de la función de respuesta y la diferencias entre las tensiones principales en el plano xy. La simulación muestra que la dirección en que se encuentra la anisotropía magnética con respecto a los ejes de referencia, no depende de las variaciones de Δµ, y que existe una dependencia lineal entre la amplitud de la señal FEM en cada circuito de lectura con respecto a los valores de Δµ. Se pudo observar que el valor calculado de α se mantiene constante ante los cambios de la anisotropía magnética introducidos en la simulación, solo dependiendo de la dirección en que Δµ se encuentra. Se concluyó que el nuevo sensor SAM puede ser considerado como si estuviese integrado por tres sensores SAM0 orientados a 10, demostrándose que es posible el empleo de la metodología utilizada para determinación de las tensiones residuales con el sensor SAM0 en el nuevo sensor avanzado SAM para los dos casos estudiados. Referencias 1. Hauk, V. Structural and residual stress analysis by nondestructive methods. 1ra ed. E. U. A.: Elsevier. 1997, p ISBN Withers, P., Turski, M., Edwards, L., et al. Recent advances in residual stress measurement. International Journal of Pressure Vessels and Piping. 008, vol. 85, n o 3, p ISSN DOI: /j.ijpvp Yamada, H., Ushiyama, S., Takeuchi, K., et al. Non-contact measurement of bending stress using a magnetic anisotropy sensor. IEEE Transactions on Magnetics. 1987, vol. 3, n o 5, p ISSN DOI /TMAG Govindaraju, M., Katragadda, G. y Wallace, J. L. Finite element analysis of a magnetic sensor to detect permeability changes due to residual stresses in ferromagnetic materials. SPIE Proceedings. 1998, vol. 3396, n o 47, p ISSN X. DOI / Bathe, K. Finite Element Method. En: Wiley encyclopedia of computer science and engineering. USA: John Wiley & Sons Ed ISBN DOI / ecse Sakai, Y., Unishi, H. y Yahata, T. Non-destructive method of stress evaluation in line pipes using magnetic anisotropy sensor. JFE Technical Report. 004, n o 3, p ISSN Ingeniería Mecánica. Vol. 16. No., mayo-agosto, 013. ISSN

11 Simulación numérica de un nuevo sensor de anisotropía magnética 7. Liu, H., Dun, C., Dou, L., et al. Theoretical analysis of magnetic sensor output voltage. Journal of Magnetism and Magnetic Materials. 011, vol. 33, n o 1, p ISSN DOI /j.jmmm Cruz, C., Suarez, J., Fernández, L., et al. Magnetoelastic sensor. En: Fourth European Conference on Residual Stresses. Cluny, France. 1996, p Cruz, C., López, N., Martín, J., et al. Desarrollo y consolidación de una versión avanzada del equipo medidor de tensiones residuales TENSO. En: Congreso Iberoamericano de Metalurgia y Materiales. IBEROMET 006, Cuba. 006, p ISBN Herrera, V., Cruz, C., Sierra, J., et al. Use of X ray diffraction and magnetoelastic effect for assessment of microstructural parameters and residual stresses in plastically deformed steel bars. En: XIII Workshop on Nuclear Physics and VII International Symposium on Nuclear and Related Techniques WONP-NURT 011. Cuba. 011, p ISBN: Fernández, T., Rodríguez, I., Quintana, R., et al. Determinación de tensiones por rayos X del Acero AISI 1045 deformado por rodillo. Ingeniería Mecánica. 01, vol. 15, n o 1, p [Consultado el: 1 de Marzo de 01]. Disponible en: ISSN Montgomery, D., Runger, G. y Hubele, N. Engineering statistics. Fourth Ed. USA: John Wiley & Sons. Ed. 008, p ISBN , Ingeniería Mecánica. Vol. 16. No., mayo-agosto, 013. ISSN