Carrillo, José

Fernández Merodo, José Antonio

Pastor Pérez, Manuel

Repositorio da Universidade da Coruña

Ferroviaria '98 415 DINACAT 1.0 - ANÁLISIS MEDIANTE ELEMENTOS FINITOS DEL ACOPLAMIENTO DINÁMICO CATENARIA-PANTÓGRAFO EN LÍNEAS FERROVIARIAS José Antonio Fernández Merodo; Manuel Pastor Pérez Centro de Técnicas Aplicadas. CEDEX, Ministerio de Fomento. José Carrillo INECO. RESUMEN El trabajo propuesto tiene como objeto analizar el comportamiento del sistema mecánico catenaria pantógrafo. Se presenta un mode-lo numérico que simula en el tiempo la interacción entre ambos sistemas. Se estudia de esta manera el número y el tiempo de per-didas de contacto entre línea de contacto y pantógrafo para garan-tizar, en el caso de circulaciones a alta velocidad, la correcta trans-misión de corriente, manteniendo, por otra parte, un desgaste acep-table tanto en la línea como en el pantógrafo. DINACAT 1.0, FINITE ELEMENT SIMULATION OFTHE DYNAMIC CATENARY-PANTOGRAPH INTERACTION. ABSTRACT The purpose of this work is to anal y se the behaviour of the mecanic system: catenary-pantograph. A numerical technique to model, in the time domain, the interaction between both systems is presented. In order to ensure the correct electric transmission, in high speed traffic, the losses of contact between the contact line and the pantograph are studied. The number of ocurren ces of contact losses and the time when these take place, are analyzed. The fact that the 416 Universidade da Coruiia - Colegio de Ingenieros de Caminos, Canales y Puerws degree to which the cables and the pantograph wear out should be limitecl to an acceptable leve], is also taken into account in the stucly. Palabras clave: catenaria, pantógrafo, elementos finitos, algoritmos no lineales, algoritmos dinámicos, algoritmos de contacto. 1. INTRODUCCIÓN Las condiciones de transportes actuales nos obligan a crear sistemas de transporte más rápidos y mejores. Uno de los problemas que se plantean a la hora de mejorar el ferroca-rril eléctrico es el problema de la transmisión de energía; la transmisión de corriente de una línea de contacto al pantógrafo de la locomotora en el caso de circulaciones a alta velocidad es uno de los temas al que hay que dedicar una particular atención. Para garantizar una captación satisfactoria, el número y el tiempo de perdidas de contac-to entre línea de contacto y pantógrafo no tienen que ser demasiado importantes. Toda perdida de contacto genera además la formación de un arco voltaico y si este fenómeno se produce regularmente en puntos determinados de la catenaria, el hilo pierde sus capa-cidades mecánicas y eléctricas y debe ser sustituido. El objetivo de este trabajo es el desarrollo de un modelo numérico que permita modelar el comportamiento mecánico del sistema catenaria-pantógrafo teniendo en cuenta la mayoría de factores que pueden influir en el comportamiento vibratorio de este último. Hay que crear un instrumento que pueda simular, para una estructura de catenaria y un pantógrafo determinados, el paso ele un tren y así poder obtener las apreciaciones y conclusiones necesarias sobre modificaciones constructivas sin tener que efectuar la has-ta ahora necesaria, simulación real. El modelo de los dos sistemas estudiados tiene que ser lo más general posible para acer-carse lo más posible a la realidad. En este sentido, los elementos finitos pueden presen-tarse como una herramienta ele simulación muy útil. Por ello se ha desarrollado el pro-grama WINCAT, en el que en la actualidad se están implementando aspectos como el efecto del viento, etc. 2. ESTUDIO ESTÁTICO DE LA CATENARIA En primer lugar, se describirá la catenaria estudiada CR220 cuyo comportamiento será a continuación modelado siendo este estudio esencial para poder continuar más tarde con el análisis dinámico y el análisis del contacto. El punto ele partida es la obtención del estado de equilibrio de la catenaria, a partir del cual se realizará el estudio dinámico. A-------Ferroviaria '9R 417 La catenaria CR220 es una catenaria simple poligonal atirantada, formada por un susten-tador suspendido y dos hilos de contacto, sin péndola en Y, con flecha inicial de los hilos de contacto y regulación automática de la tensión mecánica. 1 vano de la catenaria CR220 En una primera etapa se estudian los movimientos en un plano vertical. El descentramiento lateral de ± 20 cm en todos los apoyos (del cable sustentador y de los hilos de contacto) no ha sido considerado ya que su influencia en el comportamiento vibratorio vertical para una longitud de vano de 60 m es mínima. Los hilos de contacto y el cable sustentador pueden ser representados por 2 tipos de elemento: el elemento tipo barra pretensada y el elemento tipo viga pretensada. Las péndolas son de Cu extraflexible, se representan mediante: elementos de tipo cable. Un estudio no lineal es necesario para poder introducirlos ya que estos últimos tipos de elemento sólo trabajan a tracción, las péndolas pueden quedarse desconectadas. Condiciones de contorno: los postes situados cada 60 m, son considerados como puntos fijos del cable sustentador en el plano vertical. Si se quisiere tener en cuenta la rigidez de la consola, basta con introducir un muelle vertical que sostenga al cable sustentador. En cada extremo de los hilos de contacto y del cable sustentador se colocan las poleas de regulación de la tensión. Estas poleas son de momento modeladas por puntos tijos, ya que el pantógrafo no está en contacto con el último tramo del último vano. En realidad el estudio se hace en los vanos centrales del cantón ya que estos, como veremos, son los más representativos. Los elementos aislantes, las grifas, pueden ser introducidos como masas puntuales. El efecto de los brazos de atirantado también puede ser incluido fácilmente. El modelo matemático se desarrolla utilizando el método de los elementos finitos. Las ecuaciones de equilibrio estático son discretizadas en primer lugar considerando un modelo lineal clásico y a continuación. se introducen las hipótesis no lineales y se desarrolla un algoritmo no lineal. 418 Universidade da Corw1a - Colegio de Ingenieros de Caminos, Canales y Puertos Modelo lineal clásico: donde U es el vector con las tres componentes globales de los desplazamientos en los nodos (es de dimensión 3n sin es el número de nodos del elemento, 2 translaciones y una rotación), K es la matriz de rigidez obtenida por el ensamblado y Fes el vector de fuerzas aplicadas que incluye: i) el efecto de las fuerzas de masa, ii) el efecto de las fuerzas de superficie, iii) fuerzas puntuales aplicadas, iv) el efecto de las tensiones iniciales. Modelo no lineal: con donde K, es la matriz de rigidez tangente (función también de las tensiones en los cables) y F- F(i-l) representa el balance entre fuerzas externas e internas. Este modelo es un modelo iterativo que tiene que ser resuelto por un método como el de Newton-Rapson hasta obtener la convergencia de la solución deseada. El primer objetivo del análisis es la obtención de la posición de equilibrio a partir de la cual deben efectuarse todos los cálculos posteriores. Esta posición de equilibrio tiene que acercarse lo más posible a la posición real de la catenaria. Posición de equilibrio de 8 vanos de la catenaria CR220 Ferroviaria '91? 419 Los vanos centrales son los más representativos y en ellos la flecha es de 36 mm. Una de las pruebas de rutina realizadas en las catenarias en funcionamiento es la obten-ción de los llamados diagramas de flexibilidad. Consiste en el desplazamiento debido a una carga vertical dirigida hacia arriba que recorre los hilos de contacto. De esta forma se identifican tanto las propiedades mecánicas generales como aspectos singulares como por ejemplo. puntos duros en la misma. Se da a continuación un ejemplo que ilustra la t1exibilidad obtenida para el caso particular estudiado. Diagrama de flexibilidad de la catenaria CR220 con una tuerza de 1 DON. Desplazamientos en m (2 vanos). 3. ESTUDIO DINÁMICO DE LA CATENARIA Durante el estudio dinámico las fuerzas exteriores dependen del tiempo. y el equilibrio del sistema tiene que tener cuenta tanto de las fuerzas de inercia como del amortigua-miento que disipa energía durante las vibraciones. La ecuación de equilibrio dinámico puede escribirse ahora como: M ·U(t)+C.U(t)+K ·U(t)= F(t) donde U y U contiene las aceleraciones y las velocidades de los nodos considerados, M es la matriz de masas y C la matriz de amortiguamiento. El método de resolución escogido para resolver estas ecuaciones de segundo grado ha sido el método de Newmark, uno de los más utilizados por ser un método incondicional-mente estable, adecuado a problemas de propagación de ondas. Los algoritmos de resolución han sido implementados tanto para el caso lineal como para el no lineal. 420 Universidade da Coruña - Cole!?ÍO de ln/{enieros de Caminos, Canales y Puertos Para calcular la velocidad de propagación de una onda en la catenaria CR220 se parte del estado inicial correspondiente al equilibrio, imponiendo un desplazamiento inicial en los dos nodos céntricos de los hilos de contacto. A continuación se da el estado de la catenaria para diferentes tiempos. En el tiempo= 0.08 s En el tiempo = O. Bs En el tiempo = 2s Deformada de la catenaria CR220 en diferentes tiempos (factor desplazamientos x5) La velocidad de propagación de una onda en la catenaria CR220 calculada por el progra-ma WINCAT es de 120 m/s, siendo esta velocidad idéntica a la obtenida por la aproxima-ción analítica. v= ltensiones masa lineal 4. ESTUDIO DEL PANTÓGRAFO Los pantógrafos tienen que cumplir la misión de establecer y mantener el contacto con la catenaria. El diseño de estos pantógrafos tiene como objetivo mantener lo más constante posible la fuerza estática que estos realizan sobre la catenaria. Por eso suelen ser dividi-dos en dos partes, una parte inferior, bastante más pasiva, que compensa las A diferencias de altura de la catenaria lo largo de grandes distancias (paso a Ferroviaria '98 nivel, túneL) y una parte superior, más ligero, que tiene que amortiguar los mo-vimientos de pequeña amplitud pero rá-pidos, de Jos hilos de contacto (imperfec-ciones de la altura del cable, flecha, vi-braciones el techo de la locomotora, pro-ceso de despegues ... ). 421 ¡ h2 _ -_ hZ El pantógrafo estudiado, modelo FAIVELEY GPU, es un pantógrafo pasi-vo a dos pisos que se puede idealizar como un sistema con varias masas, mue-lles y amortiguadores: Pantógrafo modelo Faiveley GPU M penr . . W + CP""'· . W + K 1wnt. · W = FI"'"'· donde M pent. representa la matriz de masas, Cpenr. representa la matriz de amortiguamientos, Kl'""'representa la matriz de rigidez, FP""'· representa el vector de fuerzas exteriores. 5. ESTUDIO DE LA INTERACCIÓN CATENARIA-PANTÓGRAFO Después de haber analizado el comportamiento dinámico de cada uno de los sistemas estudiados: el pantógrafo y la catenaria, por separado, se estudiará en este capítulo la forma en la que estos dos sistemas interaccionan entre si en el dominio del tiempo. La fuerza de contacto que actúa entre el pantógrafo y los hilos de contacto permite aco-plar a los dos sistemas. En caso de contacto, Jos dos sistemas están acoplados y es equivalente estudiar cada sistema considerando la misma fuerza de contacto. Los desplazamientos de la cabeza de pantógrafo y de los puntos de la catenaria donde se localiza el contacto en ese instante tienen que ser iguales. En el caso de no contacto, Jos dos sistemas están desacoplados , vibran libremente y la fuerza de contacto es nula. Si en un momento determinado no hay contacto, y de la solución se obtiene que la coor-denada final de la cabeza del pantógrafo es superior a la coordenada del punto del hilo de contacto se produce recuperación de contacto. Si en un momento determinado hay con-422 Universidade da Corwia - Colegio de Ingenieros de Caminos. Canales y Puertos tacto, y de la solución se obtiene que la fuerza de contacto es negativa se produce perdida de contacto. La catenaria CR220 está compuesta por dos hilos de contacto. la condición de contacto tiene que ser respetada en los dos hilos de la catenaria. Para una velocidad del tren igual a 200 km/h, representamos la deformada de la catenaria en diferentes momentos. en el tiempo= 0.04 s en el tiempo = 1 s en el tiempo = 2 s en el tiempo = 3 s Deformada de la catenaria CR220 en diferentes tiempos P-------Ferroviaria '98 423 f-uerzas de contacto para cada hilo de contacto en N tunc10n del tiempo Fuerza de contacto total en N función del tiempo Gracias a las representaciones gráficas de los desplazamientos de la cabeza del pantógra-fo y de las fuerzas de contacto en función del tiempo. podemos analizar y de esta manera juzgar del buen o del mal comportamiento de un sistema de catenaria y pantógrafo deter-minados. Para que los despegues del pantógrafo con los hilos de contacto sean mínimos, los des-plazamientos verticales del pantógrafo y las fuerzas de contacto tienen que permanecer lo más constantes posibles a lo largo del tiempo (acercándose a curvas horizontales). El programa WINCAT puede presentarse de este modo como un instrumento muy útil a la hora de diseñar y optimizar futuras catenarias y pantógrafos. Para la catenaria CR220 y el pantógrafo FAIVELEY GPU, para una velocidad de funcio-namiento de 220 km/h, el comportamiento parece aceptable puesto que no se produce ningún despegue o de fuerza de contacto inferior a SON. La fuerza media ele contacto es de 158.87 N y la desviación estándar de 23.83 N. 424 Universidade da Corwia - Cole¡;io de ln¡;enieros de Caminos. Canales y Puertos 6. CONCLUSIONES GENERALES, VÍAS DE DESARROLLO El programa WINCAT es capaz de representar, de manera bastante sencilla. el comporta-miento estático y vibratorio del sistema mecánico catenaria-pantógrafo. La catenaria ha sido discretizada por el método de los elementos finitos y el pantógrafo representado por un conjunto de masas. muelles y amortiguadores, resolviéndose el pro-blema planteado con un algoritmo de contacto dinámico no lineal. Gracias a este programa es posible, de manera bastante precisa. analizar y así prever el buen o el mal funcionamiento del sistema predefinido. El programa WINCAT puede presentarse de este modo como un instrumento muy útil a la hora de diseñar y optimizar futuras catenarias y pantógrafos. Se ha validado el programa comparando sus resultados con otras simulaciones numéri-cas. Por otra parte, en el futuro se compararán estos resultados con otros obtenidos en simulaciones reales con otros tipos de pantógrafos y catenarias a varias velocidades de funcionamiento. Se prevé mejorar el modelo desarrollando el programa para el caso en tres dimensiones. Se podrá de esta manera tener cuenta de un máximo de hipótesis no consideradas hasta ahora, como incluir la int1uencia de las acciones climáticas como la inducida por el viento. 7. AGRADECIMIENTOS Los autores de este trabajo desean agradecer a todas aquellas personas sin cuyo apoyo, unas veces moral, otras técnico, este artículo no hubiera visto la luz. En primer lugar agradecemos a todo el equipo técnico de INECO las largas horas de trabajo en las que pacientemente nos introdujeron en el mundo de las catenarias. Por otra parte, a D. Julian Santos y D. Ignacio Alonso de la Dirección General de Ferrocarriles y Transportes por Carretera. por su apoyo incondicional. Finalmente, a D. Jesús Montesinos, de la Jefatura de Albacete, por los ensayos y el tiempo que dedicó a los autores de este trabajo. /-i-------r: Ferroviaria '98 425 BIBLIOGRAFÍA Un programa de elementos finitos que simula el comportamiento dinámico del sistema mecánico catenaria-pantógrafo. Manuel PASTOR PÉREZ, José Antonio FERNÁNDEZ MERO DO CEDEX, Informe. Junio 1997 Methodes de calcul du comportement vibratoire de la caténaire et du pantographe dans le cadre de la circulation des trains a grande vitesse. FISCHER Wolfgang. SNCF, Traduction no 67-82. 1975 Perfectionnement des pantographes et des caténaires et méthode de calcul de leurs interactions mutuelles aux grandes vitesses de marche. I.A. BELYAEV, VA. YOLOGINE, A.Y. FREIFELD. Rail Intemational. Janvier 1978 Les études recentes de pantographes a la SNCF. Jean DAFFOS, Michel GARDOU Revue Générale des Chemins de Fer. Octubre 1986 Histoire de la caténaire 25 kV. Joseph LUPPI, Jean Pierre LAMON. Revue Générale des Chemins de Fer. Marzo 1992 Proyecto:Estudio dinámico a velocidades elevadas de la catenaria tipo Renfe Moderniza-da y Compensada. Universidad de Navarra. 1982 Servo actuated railway pantograph for high-speed running with constant constant force. G. GALEOTTI, M. GALANTI, S MAGRINI, P TONI. Journal of Rail and Rapid Transit. 1993 Periodical Dynamic Stabilities of a Catenary-Pantograph System. K. MANABE. QR of RTRI, Vol.35, N"2, Mayo 1994 Elasticita de oscillazioni delle condutture di contatto delle linee ferroviarie elettrificate. L. PASCUCCI Ingeniería Ferroviaria. Marzo 1962. 426 Universidade da Corw!a - Colegio de Ingenieros de Caminos, Canales y Puertos Velocita crtiche perle condutture di contatto delle linee ferroviare elettrificate L. PASCUCCI lngenieria Ferroviaria. Febrero 1964. Movimenti delle condutture di contatto delle linee ferro vi are elettrificate, alle al te velocita. L. PASCUCC! Ingenieria Ferroviaria. Enero 1967. Comportamento dinamico del sistema pantografi catenaria. Cesare BIANCHI, Giuliano TACC!. Alfredo VAN DI La Tecnica Professionale. Julio 1990 Captazione di corrente ad Alta Velocita. C. BIANCHI, A. LUZI, E. TASS!NI. lngenieria Ferroviaria. Junio 1993. El método de los elementos finitos. Vol. 1-2 O.C. Zienkiewics y R.L. Taylor Me. Graw Hill The finite element method Thomas J.R. Hughes Prentice Hall Finite Element Procedures in Engineering Analysis K.J. Bathe Prentice Hall Non-linear Finite Element Analysis of Solids and Structures. Vol 1-2 M.A. CR!SFIELD. WILEY The ana1ysis of cable and catenary structures. P. BROUGHTON, P. NDUMBARO. Thomas Telford. 4-------

DINACAT 1.0. Análisis mediante elementos finitos del acoplamiento dinámico catenaria-pantógrafo en líneas ferroviarias

http://ruc.udc.es/dspace/bitstream/2183/10652/1/CC%2041%20art%2039.pdf

DINACAT 1.0. Análisis mediante elementos finitos del acoplamiento dinámico catenaria-pantógrafo en líneas ferroviarias

Abstract

Similar works

Full text

Available Versions

Repositorio da Universidade da Coruña