Investigación sobre modelado de redes térmicas equivalentes para especies raras.

Scientific Reports volumen 12, número de artículo: 18088 (2022) Citar este artículo

699 Accesos

1 Citas

Detalles de métricas

De crucial importancia para el diseño de transductores magnetoestrictivos gigantes (GMT) es analizar de forma rápida y precisa la distribución de temperatura. Con las ventajas de un bajo coste de cálculo y una alta precisión, se ha desarrollado el modelado de redes térmicas para el análisis térmico de GMT. Sin embargo, los modelos térmicos existentes tienen sus limitaciones para describir estos complicados comportamientos térmicos en los GMT: la mayoría de las investigaciones se centraron en el estado estacionario, que es incapaz de capturar las variaciones de temperatura; Generalmente se supone que la distribución de temperatura de las varillas magnetoestrictivas gigantes (GMM) es uniforme, mientras que el gradiente de temperatura en la varilla GMM es notable debido a su mala conductividad térmica; La distribución no uniforme de las pérdidas del GMM rara vez se introduce en el modelo térmico. Por lo tanto, en este artículo se establece un modelo de red térmica equivalente transitoria (TETN) de GMT, considerando los tres aspectos antes mencionados. En primer lugar, basándose en la estructura y el principio de funcionamiento de un GMT de vibración longitudinal, se llevó a cabo un análisis térmico. A continuación, de acuerdo con el proceso de transferencia de calor de GMT, se estableció el modelo TETN y se calcularon los parámetros correspondientes del modelo. Finalmente, la precisión del modelo TETN para el análisis temporal y espacial de la temperatura del transductor se verifica mediante simulación y experimentación.

El material magnetoestrictivo gigante (GMM), concretamente Terfenol-D, tiene las ventajas de una gran magnetoestricción y una alta densidad de energía. Estas características únicas pueden aprovecharse para permitir el desarrollo del transductor magnetoestrictivo gigante (GMT) que puede usarse en una amplia gama de aplicaciones, como transductores acústicos submarinos, micromotores, actuadores lineales, etc.1,2.

De particular preocupación es el posible sobrecalentamiento de los GMT submarinos, que generan un calor considerable debido a su alta densidad de potencia disipada cuando funcionan a máxima potencia y con un tiempo de excitación prolongado3,4. Además, las características de salida del GMT están estrechamente relacionadas con la temperatura debido al gran coeficiente de expansión térmica y su alta sensibilidad a la temperatura externa5,6,7,8. Revisando las publicaciones técnicas, los métodos para afrontar el análisis térmico GMT se pueden dividir en dos categorías principales9: métodos numéricos y de parámetros agrupados. El método de elementos finitos (FEM) es uno de los métodos de análisis numérico más utilizados. Xie et al.10 utilizaron FEM para modelar la distribución de la fuente de calor del actuador magnetoestrictivo gigante y realizaron el control de temperatura del actuador y el diseño del sistema de enfriamiento. Zhao et al.11 crearon una simulación FEM de campo de flujo turbulento acoplado y campo de temperatura y construyeron un dispositivo de control de temperatura de componentes inteligente GMM basado en los resultados de la simulación FEM. Sin embargo, FEM es muy exigente en términos de configuración del modelo y tiempo de cálculo. Por esta razón, FEM se considera un valioso apoyo para los cálculos fuera de línea, generalmente durante la etapa de diseño del transductor.

Un método de parámetros agrupados, a menudo denominado modelo de red térmica, se utiliza ampliamente en el análisis dinámico térmico en virtud de su forma matemática simple y su rápida velocidad de cálculo12,13,14. Este método ha jugado un papel esencial en la solución del problema de limitación térmica de los motores15,16,17. Mellor18 utilizó por primera vez un circuito térmico equivalente a T mejorado para simular el proceso de transferencia de calor de un motor. Verez et al.19 establecieron un modelo de red térmica tridimensional para máquinas síncronas de imanes permanentes de flujo axial. Boglietti et al.20 propusieron cuatro modelos de redes térmicas de diferentes complejidades para la predicción de transitorios térmicos de corto plazo en el devanado del estator. Finalmente, Wang et al.21 establecieron los circuitos térmicos equivalentes detallados para cada componente de la máquina síncrona de imanes permanentes y resumieron las ecuaciones de resistencia térmica. El error se puede controlar dentro del 5% en condiciones nominales.

En la década de 1990, el modelo de red térmica comenzó a aplicarse a transductores de baja frecuencia y alta potencia. Dubus et al.22 construyeron un modelo de red térmica para describir la transferencia de calor en un vibrador longitudinal de doble extremo y un transductor flextensional de clase IV en estado estacionario. Anjanappa et al.23 utilizaron el modelo de red térmica para realizar un análisis térmico bidimensional en estado estacionario de un miniactuador magnetoestrictivo. Zhu et al.24 establecieron el modelo de cálculo de desplazamiento y resistencia térmica equivalente en estado estacionario de un GMT para estudiar la relación entre la deformación térmica del Terfenol-D y los parámetros del GMT.

En comparación con la aplicación del motor, la estimación de la temperatura GMT es más complicada. La mayoría de los componentes del motor considerados a la misma temperatura generalmente se han simplificado a un solo nodo debido a la excelente conductividad térmica y magnética de los materiales utilizados13,19. Sin embargo, debido a la mala conductividad térmica del GMM, la suposición de una distribución uniforme de la temperatura ya no es válida. Además, el GMM tiene una permeabilidad muy baja, por lo que la generación de calor resultante de las pérdidas magnéticas normalmente no es uniforme a lo largo de la varilla del GMM. Además, la mayoría de las investigaciones se centraron en modelos de estado estacionario que son incapaces de capturar las variaciones de temperatura durante el funcionamiento de los GMT.

Para hacer frente a los tres desafíos técnicos antes mencionados, este artículo toma como objeto de investigación la vibración longitudinal GMT y modela con precisión diferentes componentes del transductor, especialmente la varilla GMM. Se establece el modelo completo de red térmica equivalente transitoria (TETN) del GMT. Se construyó un modelo FEM y una plataforma experimental para verificar la precisión y efectividad del modelo TETN para el análisis temporal y espacial de la temperatura del transductor.

Las dimensiones de la estructura y geometría de la vibración longitudinal GMT se muestran en las figuras 1a y b, respectivamente.

Vibración longitudinal GMT: (a) diagrama de configuración, (b) diagrama de dimensiones.

Los componentes principales incluyen la varilla GMM, la bobina de excitación, los imanes permanentes (PM), los yugos magnéticos, las placas de cubierta, el manguito y el resorte de disco. La bobina de excitación y los PM proporcionan campos magnéticos alternos y polarizados en CC para la varilla GMM, respectivamente. Los yugos magnéticos y la carcasa compuesta por placas de cubierta y manguito están hechos de hierro dulce, DT4, con una alta conductividad magnética. Junto con la varilla GMM y los PM, se forma un circuito magnético cerrado. La varilla de salida y la placa de presión están hechas de acero inoxidable 304 no magnético. Con el resorte de disco se puede aplicar una tensión previa estable a la varilla. Cuando pasa una corriente alterna a través de la bobina de excitación, la varilla del GMM vibraría en consecuencia.

En la Fig. 2 se ilustra el proceso de transferencia de calor dentro del GMT. La varilla del GMM y la bobina de excitación son las dos principales fuentes de calor del GMT. La bobina transfiere su calor a la carcasa por convección a través del aire interno y hacia la placa de cubierta por conducción. La varilla GMM generaría pérdidas magnéticas bajo la acción del campo magnético alterno y transferiría el calor a la carcasa por convección a través del aire interno y a los PM y yugos magnéticos por conducción. El calor transferido a la carcasa se disipa luego hacia el exterior mediante convección y radiación. Cuando el calor generado es igual al calor transferido, la temperatura de cada componente del GMT alcanza un estado estable.

Proceso de transferencia de calor en la vibración longitudinal GMT: (a) Diagrama de flujo de calor, (b) Rutas principales para la transferencia de calor.

Excepto el calor generado en la bobina de excitación y la varilla del GMM, todos los componentes dentro del circuito magnético cerrado también sufrirán pérdidas magnéticas. Por lo tanto, los PM, los yugos magnéticos, las placas de cubierta y la funda están todos laminados para suprimir las pérdidas magnéticas del GMT.

Los pasos principales para construir un modelo TETN para el análisis térmico del GMT son los siguientes25: en primer lugar, agrupar componentes que tengan temperaturas similares y representar a cada uno como un único nodo en la red; luego, estos nodos se acoplan con expresiones de transferencia de calor apropiadas, que representan la transferencia de calor por conducción y convección entre nodos. Mientras tanto, la fuente de calor correspondiente y la capacidad calorífica de cada componente se conectan en paralelo entre el nodo y el voltaje cero de tierra común para construir el modelo de red de calor equivalente. El siguiente paso es calcular los parámetros de la red térmica de cada componente del modelo, incluidas las resistencias térmicas, las capacidades caloríficas y las pérdidas de potencia. Finalmente, el modelo TETN se implementa en SPICE para su simulación. Y se puede obtener la distribución de temperatura en cada componente del GMT y su variación en el dominio temporal.

Para facilitar el modelado y el cálculo, es necesario simplificar el modelo térmico e ignorar las condiciones de contorno que tienen poco impacto en los resultados18,26. El modelo TETN propuesto en este artículo se establece sobre la base de los siguientes supuestos:

Considere únicamente las pérdidas de varilla y bobina, e ignore la pérdida de otros componentes.

Se ignora la influencia de la temperatura sobre los parámetros térmicos de diferentes materiales24,27.

Se ignora la transferencia de calor por radiación28. Sólo se consideran la conducción y la transferencia de calor por convección.

Se ignora la resistencia térmica de contacto entre diferentes componentes.

En GMT que tienen devanados devanados aleatoriamente, es imposible o indeseable modelar la posición de cada conductor individual. En el pasado, se han desarrollado varias estrategias de modelado para modelar la transferencia de calor y la distribución de temperatura dentro de un devanado: (1) conductividad térmica compuesta; (2) ecuaciones directas basadas en geometrías de conductores; (3) Circuito térmico equivalente a T29.

La conductividad térmica compuesta y las ecuaciones directas pueden considerarse soluciones más precisas que el circuito equivalente en T, pero dependen de varios factores, como los materiales, las geometrías de los conductores y la cantidad de aire residual en los devanados, que son difíciles de determinar29. Por el contrario, el circuito térmico T-equivalente, aunque como modelo aproximado, resulta más conveniente30. Se puede aplicar a la bobina de excitación en la vibración longitudinal GMT.

En la Fig. 3 se muestra un componente cilíndrico hueco general utilizado para representar la bobina de excitación y su circuito térmico equivalente en T que se deriva de las soluciones de las ecuaciones de conducción de calor. Se supone que el calor fluye dentro de la bobina de excitación en dirección radial y las direcciones axiales son independientes. Se ignora el flujo de calor circunferencial. En cada circuito equivalente a T, dos de los terminales representan las temperaturas superficiales apropiadas del componente y el tercer terminal T6 representa la temperatura promedio del componente. Las pérdidas P6 del componente se inyectan como fuente puntual en el nodo de temperatura media que se calculan en “Cálculo de pérdidas de calor de la bobina de excitación”. En el caso del modelado transitorio, la capacidad calorífica C6, que viene dada por la ecuación. (1), también se suma al nodo de temperatura media.

donde cec, ρec y Vec representan la capacidad calorífica específica, la densidad y el volumen de la bobina de excitación.

Las resistencias térmicas para el circuito térmico equivalente en T de la bobina de excitación con la longitud lec, la conductividad térmica λec, el radio exterior rec1 y el radio interior rec2 se presentan en la Tabla 1.

Bobina de excitación y su circuito térmico equivalente en T: (a) Componente cilíndrico hueco general, (b) Circuitos térmicos equivalentes en T axiales y radiales independientes.

También se ha demostrado que el circuito equivalente en T es preciso para otras fuentes de calor cilíndricas13. Como principal fuente de calor del GMT, la distribución de temperatura de la varilla GMM no es uniforme debido a su baja conductividad térmica y la discrepancia es particularmente pronunciada a lo largo del eje de la varilla. Por el contrario, la falta de uniformidad radial puede ignorarse porque el flujo de calor radial de la varilla GMM es mucho menor que el flujo de calor radial31.

Para representar con precisión el nivel de discretización axial de la varilla y obtener la temperatura más alta, la varilla GMM está representada por n nodos igualmente espaciados en la dirección axial, y el número de nodos n para el modelado de varillas GMM debe ser igual a un número impar. . El circuito térmico con n circuitos térmicos axiales equivalentes a T se presenta en la Fig. 4.

Modelo de red térmica con circuito de red equivalente a n T de varilla GMM.

Para determinar el número de nodos n para el modelado de varillas GMM, el resultado FEM se ilustra en la Fig. 5 como punto de referencia. El número de nodos n se ajusta en el circuito térmico de la varilla GMM como se muestra en la Fig. 4. Cada nodo se puede modelar como un circuito equivalente en T. Comparado con el resultado del FEM, se puede ver en la Fig. 5 que uno o tres nodos no pueden reflejar con precisión la distribución de temperatura de la varilla GMM (una longitud de alrededor de 50 mm) en el GMT. Cuando n aumenta a cinco, los resultados de la simulación mejorarán mucho y se acercarán al FEM. Un mayor aumento de n también puede producir resultados aún mejores, pero ocurrirá a expensas de un mayor tiempo de cálculo. Por lo tanto, en este artículo se seleccionan cinco nodos para modelar la varilla GMM.

Distribución de temperatura de la varilla GMM en FEM y modelo de red térmica con diferente complejidad.

Según el análisis comparativo anterior, en la Fig. 6 se presenta el circuito térmico preciso de la varilla GMM. T1 ~ T5 representan las temperaturas medias de las cinco secciones (Secciones 1 ~ 5) de la varilla. P1 ~ P5 representan las potencias térmicas totales en diferentes regiones de la varilla respectivamente, que se analizan en detalle en el siguiente capítulo. C1 ~ C5 representa las capacidades caloríficas de diferentes regiones, que se pueden calcular mediante

donde crod, ρrod y Vrod representan la capacidad calorífica específica, la densidad y el volumen de la varilla GMM.

Circuito térmico equivalente con 5 nodos para la varilla GMM en el GMT.

Aplicando el mismo método que para la bobina de excitación, las resistencias térmicas de transferencia de la varilla GMM en la Fig. 6 se pueden calcular mediante

donde lrod, rrod y λrod representan la longitud, el radio y la conductividad térmica de la varilla GMM.

Para la vibración longitudinal GMT estudiada en este artículo, los componentes en reposo y el aire interno se pueden modelar con configuraciones de un solo nodo.

Se puede considerar que estas regiones constan de uno o más cilindros. Los acoplamientos de transferencia de calor por conducción pura en el componente cilíndrico se determinan a partir de la ley de conducción de calor de Fourier como

donde λnhs es la conductividad térmica del material, lnhs es una longitud axial y rnhs1 y rnhs2 son los radios exterior e interior del componente de transferencia de calor, respectivamente.

La ecuación (5) se utiliza para calcular las resistencias térmicas radiales de estas regiones representadas por RR4 a RR12 en la Fig. 7. Mientras tanto, la ecuación (5) se utiliza para calcular las resistencias térmicas radiales de estas regiones representadas por RR4 a RR12 en la Fig. 7. (6) se utiliza para calcular las resistencias térmicas axiales representadas por RA15 a RA33 en la Fig. 7.

Modelo de red térmica equivalente de la vibración longitudinal GMT.

La capacidad calorífica para el circuito térmico de un solo nodo de las regiones antes mencionadas, incluidos C7 a C15 en la Fig. 7, puede estar dada por

donde ρnhs, cnhs y Vnhs representan la longitud, la densidad de capacidad calorífica específica y el volumen, respectivamente.

Las transferencias de calor por convección a través del aire interno dentro del GMT, así como entre la superficie de la carcasa y el ambiente, se modelan con una única resistencia conductora térmica de la siguiente manera:

donde A representa la superficie de contacto y h es el coeficiente de transferencia de calor. Algunos h típicos utilizados en los sistemas térmicos se enumeran en la Tabla 232. Según la Tabla 2, los coeficientes de transferencia de calor de las resistencias térmicas RH8–RH10 y RH14–RH18 en la Fig. 7, que representaron la convección entre el GMT y el ambiente, toman la constante valor de 25 W/(m2·K). El resto de coeficientes de transferencia de calor se establece en 10 W/(m2·K).

Según el proceso de transferencia de calor interno que se muestra en la Fig. 2, el modelo TETN completo del transductor se muestra en la Fig. 7.

Como puede verse en la Fig. 7, la vibración longitudinal GMT se dividió en 16 nodos, que se muestran con puntos rojos. Los nodos de temperatura representados en este modelo corresponden a las temperaturas promedio de sus respectivos componentes. La temperatura ambiente es T0, la temperatura de la varilla del GMM T1 ~ T5, la temperatura de la bobina de excitación T6, la temperatura de los PM T7 y T8, los yugos magnéticos T9 ~ T10, la temperatura de la carcasa T11 ~ T12 y T14, la temperatura del aire interno T13, y la temperatura de la varilla de salida es T15. Además, cada nodo está conectado al potencial térmico de tierra a través de C1 ~ C15, que representan respectivamente la capacidad térmica de cada región. P1 ~ P6 representan las potencias térmicas totales de la varilla GMM y la bobina de excitación respectivamente. Además, se utilizan 54 resistencias térmicas para representar las resistencias de transferencia de calor por conducción y convección entre nodos contiguos, que se han calculado en las secciones anteriores. La Tabla 3 muestra las diferentes propiedades térmicas de los materiales que componen el transductor.

La evaluación precisa de la cantidad de pérdidas y su distribución es de gran importancia para realizar un modelado térmico confiable. Las pérdidas de calor generadas en el GMT se pueden dividir en pérdida magnética de la varilla GMM, pérdida en julios de la bobina de excitación, pérdida mecánica y pérdida adicional. Las pérdidas adicionales y las pérdidas mecánicas contabilizadas son relativamente pequeñas y pueden ignorarse36.

La resistencia de CA de la bobina de excitación incluye: resistencia de CC Rdc y resistencia superficial Rs. Para una corriente alterna sinusoidal con un valor eficaz (RMS) de I, la potencia térmica total P6 de la bobina viene dada por35

La resistencia CC Rdc se puede calcular mediante

La resistencia Rs debida al efecto piel viene dada por

donde f y N son la frecuencia de la corriente de excitación y el número de vueltas. lCu y rCu representan el radio interior y exterior de la bobina, la longitud de la bobina y el radio del cable magnético de cobre que se define por su número AWG (American Wire Gauge). ρCu es la resistividad de su núcleo. μCu es la permeabilidad de su núcleo.

El campo magnético real dentro de la bobina de excitación, un solenoide, no es uniforme en toda la dirección longitudinal de la varilla. La discrepancia es particularmente pronunciada debido a la baja conductividad magnética de la varilla GMM y el PM. Pero es longitudinalmente simétrico. La distribución del campo magnético determina directamente la distribución de la pérdida magnética de la varilla GMM. Por lo tanto, para reflejar la distribución de pérdidas real, se toman para medir tres secciones de la varilla que se muestra en la Fig. 8.

Posiciones de las bobinas de inducción en la varilla GMM.

La pérdida magnética se puede obtener midiendo el bucle de histéresis dinámica37. Los bucles de histéresis dinámica de tres secciones se midieron basándose en la plataforma experimental que se muestra en la Fig. 11. Bajo la condición de que la temperatura de la varilla GMM sea estable por debajo de 50 °C, una fuente de alimentación de CA programable (Chroma 61512) impulsa la bobina de excitación en un rango de frecuencias con una corriente de prueba para producir un campo magnético y la densidad de flujo resultante se deriva integrando el voltaje inducido en las bobinas de inducción conectadas a la varilla GMM como se muestra en la Fig. 8. Los datos sin procesar se descargaron del hicorder de memoria (Daily MR8875 -30) y procesado en el software MATLAB para derivar los bucles de histéresis dinámica medidos que se muestran en la Fig. 9.

Bucles de histéresis dinámica medidos: (a) Sección 1/5: Bm = 0,044735 T, (b) Sección 1/5: fm = 1000 Hz; (c) Sección 2/4: Bm = 0.05955 T, (d) Sección 2/4: fm = 1000 Hz, (e) Sección 3: Bm = 0.07228 T, (f) Sección 3: fm = 1000 Hz.

Según la literatura37, la pérdida magnética total por unidad de volumen Pv de la varilla GMM se puede calcular a partir de:

donde ABH es el área medida dentro de la curva BH a la frecuencia del campo magnético fm que es igual a la frecuencia de la corriente de excitación f.

Según el método de separación de pérdidas de Bertotti38, la pérdida magnética por unidad de masa de la varilla GMM Pm se puede expresar como la suma de la pérdida por histéresis Ph, la pérdida por corrientes parásitas Pe y la pérdida anormal Pa (13):

Desde la perspectiva del análisis de ingeniería38, la pérdida anormal y la pérdida por corrientes parásitas se pueden combinar en un término llamado pérdida total por corrientes parásitas. Por tanto, la fórmula de cálculo de pérdidas se puede simplificar a:

En las ecuaciones. (13)–(14), Bm es la amplitud de la densidad magnética del campo magnético de excitación. kh y kc son el coeficiente de pérdida por histéresis y el coeficiente de pérdida total por corrientes parásitas.

Según la literatura39, basándose en una fórmula simplificada de separación de pérdidas y datos de pruebas de pérdidas, se adopta el método de ajuste de curva polinomial para obtener matemáticamente el coeficiente de pérdida por histéresis y el coeficiente de pérdida total por corrientes parásitas, que se enumeran como (12) y (13):

donde a0, a1, a2, b0, b1, b2 y b3 son los parámetros relevantes de pérdida de material. Según los datos experimentales ajustados, los siete parámetros relacionados con las pérdidas en diferentes secciones se obtuvieron para diferentes secciones (Fig. 6), como se muestra en la Tabla 4.

La Tabla 5 muestra los valores de pérdida magnética P1 ~ P5 en diferentes secciones de la varilla GMM según la fórmula de separación de pérdidas ajustada en diferentes niveles de corriente.

Se puede ver en la Tabla 4 que la pérdida magnética en los extremos de la varilla del GMM es mucho menor que en el centro, lo que demuestra la evidente falta de uniformidad de la pérdida magnética a lo largo de la dirección axial de la varilla. Mientras tanto, la falta de uniformidad se amplifica continuamente con el aumento de la corriente de excitación. La distribución no uniforme de pérdidas de la varilla GMM afecta significativamente la precisión de cualquier enfoque de análisis térmico, especialmente para GMT de baja frecuencia y alta potencia.

Para verificar mejor la precisión del modelo TETN para la distribución de temperatura de la varilla GMM, se modeló la vibración longitudinal GMT empleando el software COMSOL Multiphysics que se muestra en la Fig. 10. Un sector de semisimetría con condiciones de contorno de simetría y un mallado libre de volumen generado por Sólo se utilizó el tipo de elementos piramidales (tetraédricos). La malla completa consta de 108.924 elementos de dominio, 23.976 elementos de límite y 2448 elementos de borde por completo.

Modelo de elementos finitos de vibración longitudinal GMT.

La temperatura de la varilla GMM T1 ~ T5 obtenida de la simulación FEM y el modelo TETN se comparan en la Fig. 11 bajo I = 3A yf = 1000 Hz.

Comparación de la temperatura de la varilla GMM en diferentes posiciones obtenida de la simulación FEM y el modelo TETN.

Como muestra la Fig. 11, la temperatura es más baja en el nodo 5. El modelo TETN produce una temperatura más pequeña que FEM con un error de alrededor de 0,89 °C en estado estable. Con el aumento de la temperatura del nodo, las diferencias entre el modelo TETN y FEM aumentaron a 2,31 °C y 2,93 °C respectivamente en el nodo 4 y el nodo 3. La Figura 11 muestra una buena concordancia entre los resultados de la simulación del modelo TETN y FEM con un pequeño error inferior al 2%, lo que indica que el modelo propuesto logra un resultado de simulación preciso para la distribución de temperatura de la varilla GMM.

Para validar la efectividad y precisión del modelo propuesto, se construyó una plataforma experimental de aumento de temperatura como se muestra en la Fig. 12a. Se elige una fuente de alimentación de CA programable (Chroma61512) para suministrar y controlar la corriente de excitación de la vibración longitudinal GMT. Los sensores de temperatura tipo K se utilizan para recopilar la temperatura en tiempo real en diferentes posiciones. Se utiliza un osciloscopio (Tektronix MDO34) para monitorear el voltaje y la corriente, mientras que una grabadora de memoria (Daily mr8875-30) para el almacenamiento de datos. La computadora obtiene datos de la memoria superior y se utiliza para procesar datos.

Plataforma experimental de aumento de temperatura para vibración longitudinal GMT: (a) Plataforma experimental, (b) Posiciones de prueba seleccionadas.

Comenzando con la vibración longitudinal GMT a temperatura ambiente, la prueba se realiza durante un tiempo de 20 minutos para evitar que la temperatura de la varilla del GMM alcance la mitad del punto Curie, lo que afectaría el funcionamiento confiable del GMT. Para fines de comparación, durante los ensayos experimentales se registran las temperaturas de la superficie central de la varilla, de la pared interior de la bobina de excitación y de la carcasa.

Para demostrar el rendimiento del modelo TETN para la estimación de temperatura, se midieron las temperaturas de la superficie central de la varilla, la pared interna de la bobina de excitación y la carcasa cuando oscila entre 1,5 y 3 A, como se muestra en la Fig. 12b.

Los resultados finales del modelo propuesto, FEM y la medición experimental en varias posiciones se presentan en la Tabla 6. Las diferencias máximas de temperatura final observadas entre los resultados de las simulaciones TETN y los datos medidos son 3,46 °C en la superficie central de la varilla GMM, 2,03 °C C en la pared interior de la bobina de excitación y 1,44 °C en la pared interior de la carcasa, respectivamente.

La Figura 13 muestra una comparación basada en el tiempo entre los resultados medidos y simulados en tres posiciones de prueba seleccionadas. Los resultados simulados para la varilla GMM son comparables a los resultados medidos con una ligera desviación dentro de 1,5 °C cuando la corriente de excitación es de 1,5 A y 2 A en la Fig. 13a. Con el aumento de la temperatura, las diferencias entre los resultados medidos y simulados que ocurren durante el período transitorio inicial aumentan a 5,01 °C y 15,92 °C en I = 2,5 A y 3 A, lo que probablemente se atribuye al valor menos preciso de la capacidad calorífica para Varilla GMM. Además, en realidad las pérdidas de potencia no son constantes a altas temperaturas, que se supone que no cambian según se miden a temperatura ambiente39.

Comparación de los valores medidos y los resultados del modelo de aumento de temperatura en diferentes posiciones de GMT: (a) La superficie central de la varilla GMM, (b) La pared interna de la bobina de excitación, (c) La pared interna del manguito.

Como se muestra en la Fig. 13b, el modelo propuesto es más preciso en la estimación de la temperatura de la bobina de excitación. Los errores máximos del modelo y experimento TETN están dentro de los 3 °C. Sin embargo, es evidente que la tendencia de aumento de temperatura de la bobina para FEM es bastante diferente a la del experimento. Probablemente la razón principal es que la bobina de excitación transfiere su calor por convección a través del aire interno. Sin embargo, COMSOL no puede describir con precisión el flujo de calor del aire interno. En la Fig. 13c, hay una buena concordancia entre los resultados medidos y calculados de la temperatura de la pared interior del armazón. Mientras que, con el aumento de la corriente de excitación, la temperatura medida aumenta más rápido que la temperatura calculada, ya que la radiación superficial de la carcasa, que en el modelo propuesto se supone ignorada, se vuelve cada vez más importante en la realidad.

En resumen, el método de modelado térmico propuesto se puede aplicar con éxito al GMT de baja frecuencia y alta potencia. Debería ser necesaria la optimización del modelo TETN para tener en cuenta la transferencia de calor radial si la distribución de pérdidas radiales de la varilla GMM se vuelve muy pronunciada a una frecuencia de operación más alta debido al mayor efecto de las corrientes parásitas.

En este artículo se ha establecido y aplicado un modelo TETN para estimar la distribución de temperatura de una vibración longitudinal GMT. En particular, se analizan y consideran en profundidad las características distintivas de la varilla GMM como fuente de calor y, además, la influencia de la temperatura y la distribución de la generación de calor de la varilla GMM. En este artículo, se informó el procedimiento completo para los circuitos térmicos equivalentes de cada parte del GMT y se calcularon en detalle los parámetros del modelo correspondiente. La precisión de este modelo en la estimación de la temperatura del transductor se verificó mediante simulación y experimentos FEM. Además, se analiza el análisis de errores en todas las posiciones de prueba seleccionadas, lo que debería arrojar luz sobre la mejora adicional del modelo TETN.

Los conjuntos de datos generados y/o analizados durante el estudio actual están disponibles del autor correspondiente a solicitud razonable.

Rong, C. y col. Modelado dinámico y análisis de un actuador magnetoestrictivo gigante de pila. Sens. Actuadores A. 276, 205–218 (2018).

Artículo CAS Google Scholar

Yadav, S. & Mehta, RK Modelado de vibración magnetoestrictiva y acústica en transformador convertidor. IET Electr. Aplicacion de energia 15, 332–347. https://doi.org/10.1049/elp2.12025 (2021).

Artículo de Google Scholar

Xue, G. y col. Estudios de simulación de la inyección en forma de bota de un inyector magnetoestrictivo gigante. Ciencia. Rep. 11, 22999. https://doi.org/10.1038/s41598-021-02529-z (2021).

Artículo ADS CAS PubMed PubMed Central Google Scholar

Butler, SC y Montgomery, RE A Análisis térmico de elementos transductores de anillo de alta potencia. División del Centro de Guerra Submarina Naval Newport Ri (2005).

Blottman, JB y cols. Diseño térmico de transductores activos de alta potencia con el código de elementos finitos atila. Proc. Inst. Acústico. 27(1), 87–94 (2005).

Google Académico

Debus, JC Análisis térmico en materiales piezoeléctricos y magnetoestrictivos utilizando aplicaciones ATILA del software ATILA FEM para materiales inteligentes 230–282 (Woodhead Publishing, 2013).

Google Académico

Zhao, Z. & Sui, X. Diseño y experimento de compensación de temperatura para un actuador magnetoestrictivo gigante. Ciencia. Rep. 11, 251. https://doi.org/10.1038/s41598-020-80460-5 (2021).

Artículo CAS PubMed PubMed Central Google Scholar

Nersessian, N., Or, SW y Carman, GP Caracterización magnetotermomecánica de compuestos de Terfenol-D unidos por polímeros de 1 a 3 tipos. J. Magn. Magn. Madre. 263(1–2), 101–112 (2003).

Artículo ADS CAS Google Scholar

Dubus, B. Análisis del límite térmico de proyectores de sonar de alta potencia y baja frecuencia. en el 1er Congreso Europeo sobre Acústica Subacuática (1992).

Xie, T., Wang, C., Yu, C., Xu, B. y Shi, R. Análisis y experimento de las características del aumento de temperatura de GMA. AIP Avanzado. 10(10), 105126 (2020).

ADS del artículo Google Scholar

Zhao, Z. & Sui, X. Diseño y experimento de compensación de temperatura para un actuador magnetoestrictivo gigante. Ciencia. Rep. 11(1), 1-14 (2021).

Google Académico

Vahedi, P., Ganji, B. & Afjei, E. Modelo térmico concentrado para el motor de reluctancia conmutada multicapa. IET Electr. Aplicacion de energia 14, 1873–1885. https://doi.org/10.1049/iet-epa.2019.1025 (2020).

Artículo de Google Scholar

Rostami, N. y col. Modelo térmico de parámetros concentrados para máquinas de imanes permanentes de flujo axial. Traducción IEEE. Magn. 49(3), 1178–1184 (2012).

ADS del artículo Google Scholar

Gan, Y., Wang, J., Liang, J., Huang, Z. & Hu, M. Desarrollo de un modelo de circuito equivalente térmico de un sistema de gestión térmica basado en tubos de calor para un módulo de batería con celdas cilíndricas. Aplica. Termia. Ing. 164, 114523 (2020).

Artículo de Google Scholar

Kang, M., Shi, T., Guo, L., Gu, X. & Xia, C. Análisis térmico del sistema de enfriamiento con circulación entre los orificios del rotor de PMSM cerrados basado en modelos modificados. Aplica. Termia. Ing. 1, 118054 (2022).

Artículo de Google Scholar

Cavazzuti, M., Gaspari, G., Pasquale, S. & Stalio, E. Gestión térmica de un motor eléctrico de Fórmula E: análisis y optimización. Aplica. Termia. Ing. 157, 113733 (2019).

Artículo de Google Scholar

Liu, F., Hu, J., Li, Y. y Wang, Q. Modelo térmico mejorado de motores refrigerados por aire forzado considerando la transferencia de calor en el devanado de alambre bobinado y la región final. IET Electr. Aplicacion de energia 14, 943–950. https://doi.org/10.1049/iet-epa.2019.0780 (2020).

Artículo de Google Scholar

Mellor, PH, Roberts, D. & Turner, DR Modelo térmico de parámetros agrupados para máquinas eléctricas de diseño TEFC. Proceso EEI B eléctrico. Aplicacion de energia 138(5), 205–218 (1991).

Artículo de Google Scholar

Verez, G., Tiegna, H., Barakat, G. & Hoblos, G. Modelado térmico analítico de máquinas síncronas de imanes permanentes de flujo axial. en Proc. 2012 20° Int. Conf. eléctrico. Mach. ICEM, 2799–2805 (2012).

Boglietti, A., Cavagnino, A. & Staton, D. Determinación de parámetros críticos en modelos térmicos de máquinas eléctricas. Traducción IEEE. Aplica. Ind. 44(4), 1150–1159 (2008).

Artículo de Google Scholar

Wang, Análisis del modelo térmico de parámetros agrupados del WB para PMSM 16–33 (Universidad de Zhejiang, 2014).

Google Académico

Dubus, B. Una evaluación analítica del calentamiento de proyectores de sonar de baja frecuencia. J. acústico. Soc. Soy. 95(4), 1983-1990 (1994).

ADS del artículo Google Scholar

Anjanappa, M. & Bi, J. Un estudio teórico y experimental de miniactuadores magnetoestrictivos. Madre inteligente. Estructura. 3(2), 83–91 (1994).

ADS del artículo Google Scholar

Zhu, Y. & Ji, L. Investigaciones teóricas y experimentales de la temperatura y la deformación térmica de un actuador magnetoestrictivo gigante. Sens. Actuadores A. 218, 167–178 (2014).

Artículo CAS Google Scholar

Boglietti, A. et al. Evolución y enfoques modernos para el análisis térmico de máquinas eléctricas. Traducción IEEE. Indiana electrón. 56(3), 871–882 (2009).

Artículo de Google Scholar

Xu, Y., Ai, M. & Yang, Y. Investigación de características de transferencia de calor basada en el método de red térmica en motor sumergible. En t. Trans. eléctrico. Sistema de energía. 28(3), e2507 (2018).

Artículo de Google Scholar

Debus, JC, Blottman, J. III. & Butler, S. Comportamiento térmico de dispositivos activos de alta potencia con el código de elementos finitos ATILA (análisis de transductores mediante integración de ecuaciones de LAplace). J. acústico. Soc. Soy. 120(5), 3274–3274 (2006).

ADS del artículo Google Scholar

De-xing, Z., Weifang, C. y Miaomiao, L. Un modelo de red térmica optimizado para estimar el rendimiento térmico de un par de rodamientos de bolas de contacto angular con lubricación por aire y aceite. Aplica. Termia. Ing. 131, 328–339 (2018).

Artículo de Google Scholar

Staton, D., Boglietti, A. & Cavagnino, A. Resolver los aspectos más difíciles del análisis térmico de motores eléctricos en motores de inducción industriales pequeños y medianos. Traducción IEEE. Conversaciones de energía. 20(3), 620–628 (2005).

ADS del artículo Google Scholar

Perez, IJ & Kassakian, JG Un modelo térmico estacionario para máquinas eléctricas giratorias con entrehierro suave. eléctrico. Mach. Sistema de energía. 3(3–4), 285–303 (1979).

Artículo de Google Scholar

El-Refaie, AM, Harris, NC, Jahns, TM y Rahman, KM Análisis térmico de una máquina síncrona de PM interior multibarrera utilizando un modelo de parámetros agrupados. Traducción IEEE. Conversaciones de energía. 19(2), 303–309 (2004).

ADS del artículo Google Scholar

Bahman, AS, Ma, K. & Blaabjerg, F. Un modelo térmico agrupado que incluye acoplamiento térmico y condiciones de límite térmico para módulos IGBT de alta potencia. Traducción IEEE. Electron de potencia. 33(3), 2518–2530 (2017).

ADS del artículo Google Scholar

Roy, P. Modelado térmico de motores síncronos de imanes permanentes para aplicaciones en vehículos eléctricos (tesis doctoral, Universidad de Windsor). (2020).

Jung, YG, Choi, SC, Oh, CS y Paik, UG Tensión residual y propiedades térmicas de materiales clasificados funcionalmente de circonio/metal (níquel, acero inoxidable 304) fabricados mediante prensado en caliente. J. Mater. Ciencia. 32(14), 3841–3850 (1997).

Artículo ADS CAS Google Scholar

Angara, R. Sistema de conducción de campo magnético de alta frecuencia y alta amplitud para actuadores magnetostrictivos (Universidad de Maryland, 2009).

Google Académico

Grunwald, A. & Olabi, AG Diseño de un actuador magnetoestrictivo (MS). Sens. Actuadores A Phys. 144(1), 161-175 (2008).

Artículo CAS Google Scholar

Du, R. & Robertson, P. Modelado de bucles de histéresis magnética dinámica y pérdidas de potencia en acero laminado. en Proc. UKSim-AMSS 17º Int. Conf. Computadora. Modelo. Simulación, UKSim 2015, 432–437 (2016).

Guo, P. y col. Modelo de cálculo de pérdida de energía magnética de coeficiente variable para materiales magnetoestrictivos considerando esfuerzos de compresión. AIP Avanzado. 11(3), 035227 (2021).

ADS del artículo Google Scholar

Huang, W., Xia, Z. y Guo, P. Prueba característica de alta frecuencia y cálculo de pérdidas de la aleación TbDyFe a temperatura variable. Traducción IEEE. Magn. (2021).

Descargar referencias

Centro Nacional de Investigación de Tecnología de Ingeniería de Control y Conversión de Energía Eléctrica (Universidad de Hunan), Changsha, Hunan, China

Zhihe Zhang, Xin Yang y Yukai Chen

También puedes buscar este autor en PubMed Google Scholar.

También puedes buscar este autor en PubMed Google Scholar.

También puedes buscar este autor en PubMed Google Scholar.

ZZ e YC: análisis y simulación de datos; XY: concibió el manuscrito; ZZ: escribió el manuscrito original y dibujó las imágenes del manuscrito; XY y ZC: manuscrito revisado y verificado todos los datos.

Correspondencia a Xin Yang.

Los autores declaran no tener conflictos de intereses.

Springer Nature se mantiene neutral con respecto a reclamos jurisdiccionales en mapas publicados y afiliaciones institucionales.

Acceso Abierto Este artículo está bajo una Licencia Internacional Creative Commons Attribution 4.0, que permite el uso, compartir, adaptación, distribución y reproducción en cualquier medio o formato, siempre y cuando se dé el crédito apropiado a los autores originales y a la fuente. proporcione un enlace a la licencia Creative Commons e indique si se realizaron cambios. Las imágenes u otro material de terceros en este artículo están incluidos en la licencia Creative Commons del artículo, a menos que se indique lo contrario en una línea de crédito al material. Si el material no está incluido en la licencia Creative Commons del artículo y su uso previsto no está permitido por la normativa legal o excede el uso permitido, deberá obtener permiso directamente del titular de los derechos de autor. Para ver una copia de esta licencia, visite http://creativecommons.org/licenses/by/4.0/.

Reimpresiones y permisos

Zhang, Z., Yang, X. y Chen, Y. Investigación sobre modelado de redes térmicas equivalentes para transductores magnetoestrictivos gigantes de tierras raras. Representante científico 12, 18088 (2022). https://doi.org/10.1038/s41598-022-22959-7

Descargar cita

Recibido: 31 de mayo de 2022

Aceptado: 21 de octubre de 2022

Publicado: 27 de octubre de 2022

DOI: https://doi.org/10.1038/s41598-022-22959-7

Cualquier persona con la que comparta el siguiente enlace podrá leer este contenido:

Lo sentimos, actualmente no hay un enlace para compartir disponible para este artículo.

Proporcionado por la iniciativa de intercambio de contenidos Springer Nature SharedIt

Al enviar un comentario, acepta cumplir con nuestros Términos y pautas de la comunidad. Si encuentra algo abusivo o que no cumple con nuestros términos o pautas, márquelo como inapropiado.