ESCUELA SUPERIOR POLITÉCNICA DE CHIMBORAZO

FACULTAD DE MECÁNICA

CARRERA DE INGENIERÍA AUTOMOTRIZ

ESCUELA DE INGENIERÍA AUTOMOTRIZ

CONTROL AUTOMÁTICO

SIMULACIÓN DE UNA TURBINA DE GAS

INTEGRANTES: Luis Gallegos

Lenin criollo

Jairo Benavides

Cristian Aldas

Víctor moreno

27 DE JULIO DE 2016
RESUMEN

Como es bien sabido, los motores de turbina de gas en tiempos recientes, son una fuente de
poder ampliamente usada y son la base de generación de electricidad a nivel mundial y de la
industria aeronáutica. Estos motores han demostrado ser capaces de producir potencias muy
elevadas y las investigaciones actuales están dirigidas a mejorar sus rendimientos y
condiciones de operación, para un mejor desempeño y un ahorro en combustibles fósiles.
Aun teniendo diversas ventajas, estos motores en contraparte son extremadamente costosos
ya que requieren de materiales especialmente diseñados para su fabricación y deben de
soportar condiciones de operación extremas.

En este trabajo de una forma resumida, se describe el proceso de obtención de fórmulas que
sirven para realizar la simulación en el entorno MATLAB_SIMULINK, ya que este software
nos permite realizar el análisis matemático y la simulación de posibles fallas de los motores
de turbina de gas, la interpretación y desarrollo de diagramas de bloques, la utilización de
función de transferencia y muchas otras variables que se describen posteriormente en el
presente informe.

ABSTRACT

As it is well known, in recent days gas turbine engines are widely used as power plants and
are the source of electric generation worldwide and also of the aeronautical industry. These
engines have demonstrated being capable to produce a vast volume of power and
investigations nowadays are focused in improving its efficiency and operating conditions, to
obtain a better performance and savings in fuel. Even its well known advantages, these
engines in counterpart are extremely expensive because of the specially designed materials
used in its production because they most tolerate extreme operating conditions.

In this paper a summary, the process of obtaining formulas that are used to perform
simulation in the MATLAB_SIMULINK environment , as this software allows us to perform
mathematical analysis and simulation of possible failures engines gas turbine described ,
interpretation and development of block diagrams , using transfer function and many other
variables that are described later in this report.
1. INTRODUCCIÓN

1.1 Antecedentes

SIMULACIÓN DINÁMICA DE PLANTAS DE GENERACIÓN CON TURBINA DE GAS,

DE VAPOR, DE CICLO CONVINADO Y DE GENERACIÓN.

Ciertamente la simulación tiene un papel cada vez más extenso e importante en el mundo
industrial, la posibilidad de prever, y analizar de forma inicua comportamientos de
sistemas reales ofrece evidente interés en múltiples actividades.

El avance de las técnicas numéricas y de las tecnologías de procesamiento digital y el

constante abaratamiento, trabajan en favor del uso de la simulación incluso por sectores
inesperados, junto a este fenómeno de extensión y de progreso hay que contar también
con el constante interés por la simulación que mantienen, desde hace años, sectores
fundamentales de la industria, este es el caso de centrales productoras de energía eléctrica.

La ventaja de un modelo basado en la teoría y no en la experimentación es que puede ser

usado para simular bajo todo tipo de circunstancias, (cosas que la experimentación directa
no lo permite), ampliando así el rango de aplicabilidad.

Pero la desventaja de la simulación es que tiene un margen de error con la realidad debido
a la dificultad de completar todos los aspectos concretos de una central, entre ambas
posturas de modelación, hay toda una sucesión de estados intermediarios que combinan
la teoría con la experimentación.

DISEÑO Y CONSTRUCCIÓN DE UNA TURBINA DE GAS FACULTAD NAUTICA DE

BARCELONA.

Una turbina de gas es una turbo máquina que transforma la energía de un flujo de vapor
en energía mecánica, este vapor se genera en una caldera de la que sale unas condiciones
de elevada temperatura y presión.

En la turbina se transforma la energía interna del vapor en energía mecánica que

típicamente es aprovechada por un generador para producir electricidad.
Se propone la construcción de un prototipo de turbina de vapor que tenga como sustancia
de trabajo el vapor de agua generado por una caldera, con la principal característica de
producir energía eléctrica.

El diseño de este dispositivo responde a la necesidad de buscar alternativas para la

generación de energía y que este dispositivo responda de manera correcta y sea amigable
con el medio ambiente, mediante las turbinas se genera alrededor del 80% de toda la
energía utilizada en el mundo gracias a centrales térmicas, la turbina de vapor junto con
la bomba de alimentación de agua (o combustible), la caldera y el condensador de una
central de ciclo combinado gas-vapor actúan en ciclo cerrado y termodinámicamente es
una aplicación del ciclo rankine.

1.2 Justificación

La finalidad de este trabajo de fin de curso es ampliar nuestras capacidades investigativas,

ya que es importante que como estudiantes de ingeniería se tenga la capacidad de
consultar temas en específico útiles y de importancia para el desarrollo de nuestra carrera,
así como también temas de interés que puedan servir para desarrollarlos o concluirlos de
acuerdo al caso. Así como también poner en práctica ciertos criterios y conocimientos
obtenidos a lo largo de la asignatura, tales como el manejo y entendimiento del programa
MATLAB_SIMULINK específicamente del desarrollo e interpretación de diagramas de
bloques y función de transferencia de todo tipo de sistemas.

2. OBJETIVO GENERAL

Investigar y plasmar la simulación de una turbina de gas en el entorno de MATLAB-

SIMULINK.

3. OBJETIVOS ESPECÍFICOS
o Comprender el funcionamiento de una turbina de gas mediante la simulación de la
misma en el entorno de SIMULINK, para reforzar conocimientos de la materia y de
nuevas tecnologías.
o Conocer y aplicar las distintas herramientas que posee el entorno de Simulink,
mediante la simulación de un maquina térmica denominada turbina de gas, para
reforzar y entender conocimientos adquiridos.
 o Reforzar la capacidad de investigación, mediante la utilización persistente de diversas
herramientas tecnológicas ya que esto será muy útil en nuestra vida estudiantil y
profesional.

4. DESARROLLO

Para la elaboración de este proyecto el primer paso fue consultar exhaustivamente todo lo
relacionado con una turbina de gas así como también específicamente la simulación de la
misma en el entorno de MATLAB-SIMULINK, una vez logrado esto se procede a tratar de
entender cada una sus partes.

En la siguiente figura se puede apreciar un diagrama de bloques de la función de transferencia

de una turbina de gas, el cual será el punto de partida para el estudio.

Figura1: Diagrama de bloques de la función de transferencia de una turbina de gas

Fuente: http://liaowangbj.weebly.com/independent-study.html

Es así que de acuerdo a lo consultado la turbina de gas simulada tiene relación en su salida
de energía mecánica con un generador, y por ende el estudio del autor parte desde ese punto.
Utilizando un generador sincrónico, de acuerdo con la segunda ley de Newton, la ecuación
de movimiento del rotor viene dada por:

Jαm(t) = Tm(t) − Te(t) = Ta(t)

Donde:

J= momento de inercia total de las masas giratorias

αm= aceleración angular del rotor Te = par eléctrico

Tm= par mecánico Ta = par neto

Con respecto a la aceleración angular del rotor y por definición se puede escribir que:

𝑑Wm (t) 𝑑 2 θm(t)

αm(t) = =
𝑑𝑡 𝑑𝑡 2

𝑑θm(t)
Wm(t) =
𝑑𝑡

Dónde:

Wm= velocidad angular.

θm= posición angular del rotor con respecto a un eje estacionario.

Con el fin de medir la posición angular del rotor convenientemente, se debe utilizar un eje de
referencia que debe girar de forma sincronizada en lugar de un eje estacionario. Por lo tanto,
se define:

θm = Wmsgn(t) + ȡ𝑚(𝑡)

Donde:

𝑊𝑚𝑠𝑔𝑛(𝑡)= velocidad angular sincrónica rad/s

ȡ𝑚(𝑡)= posición angular del rotor respecto a una referencia que gira sincrónicamente rad.

Por tanto, la ecuación de movimiento del rotor queda definida de la siguiente manera:

𝑑 3 θm(t) 𝑑 3 ȡ𝑚(𝑡)
𝐽 =𝐽 = Tm(t) − Te(t) = Ta(t)
𝑑𝑡 𝑑𝑡
Donde:

J= momento de inercia total de las masas giratorias

θm= posición angular del rotor con respecto a un eje estacionario.

ȡ𝑚(𝑡)= posición angular del rotor respecto a una referencia que gira sincrónicamente rad.

Tm= par mecánico Te = par eléctrico

Ta = par neto a acelerar

Es conveniente trabajar con energía en lugar del par. Por lo tanto se multiplica a la expresión
por Wm(t) y se divide para Srated, quedando de la siguiente manera:

𝐽𝑊𝑚(𝑡) 𝑑 3 ȡ𝑚(𝑡) 𝑊𝑚(𝑡)𝑇𝑚(𝑡) − 𝑊𝑚(𝑡)𝑇𝑒(𝑡) 𝑃𝑚(𝑡) − 𝑃𝑒(𝑡)

= = =
𝑆𝑟𝑎𝑡𝑒𝑑 𝑑𝑡 𝑆𝑟𝑎𝑡𝑒𝑑 𝑆𝑟𝑎𝑡𝑒𝑑

𝑃𝑚𝑝𝑢(𝑡) − 𝑃𝑒𝑝𝑢(𝑡) = 𝑃𝑎𝑝𝑢(𝑡)

Donde:

Pmpu= energía mecánica por unidad de tiempo

Pepu= energía eléctrica por unidad de tiempo. Papu= energía neta por unidad de tiempo.

Por último, es conveniente trabajar con una constante de inercia normalizada denotada con
H y viene dada por:

H Se define como 1~10

1
𝐽𝑊 3 msgn
𝐻= 3
𝑆𝑟𝑎𝑡𝑒𝑑

La constante de H tiene la ventaja de que recae dentro de un rango bastante estrecho,

normalmente entre 1 y 10 pu-s. Por lo tanto el lado izquierdo en la ecuación se convierte:
 1 3
3 𝐽𝑊 msgn ∗ 3𝑊𝑚(𝑡) ∗ 𝑑 ȡ𝑚(𝑡) = 3𝐻 𝑊𝑚(𝑡) ∗ 𝑑 ȡ𝑚(𝑡)
2 2

𝑆𝑟𝑎𝑡𝑒𝑑 𝑊 3 msgn 𝑑𝑡 𝑊 3 msgn 𝑑𝑡

Por definición la velocidad angular del rotor por unidad de tiempo es:

𝑊𝑚(𝑡)
𝑊𝑝𝑢(𝑡) =
Wmsgn

Por lo tanto la ecuación se convierte en:

3𝐻 𝑑3 ȡ𝑚(𝑡)
∗ 𝑊𝑝𝑢(𝑡) ∗ = 𝑃𝑚𝑝𝑢(𝑡) − 𝑃𝑒𝑝𝑢(𝑡) = 𝑃𝑎𝑝𝑢(𝑡)
Wmsgn 𝑑𝑡 3

Para un generador sincrónico con P polos se tiene los siguientes parámetros:

Aceleración angular eléctrica

𝑃
𝛼(𝑡) = 𝛼𝑚(𝑡)
2

Frecuencia eléctrica en radianes

𝑃
𝑊(𝑡) = 𝑊𝑚(𝑡)
2

Poder amplio

𝑃
ȡ(𝑡) = ȡ𝑚(𝑡)
2

𝑃
𝑓𝑟𝑒𝑐𝑢𝑒𝑛𝑐𝑖𝑎 𝑒𝑙é𝑐𝑡𝑟𝑖𝑐𝑎 𝑒𝑛 𝑟𝑎𝑑𝑖𝑎𝑛𝑒𝑠 − Wsgn = Wmsgn
2

La frecuencia eléctrica por unidad de tiempo es:

2
𝑊𝑚(𝑡) 𝑃 𝑊(𝑡) 𝑊(𝑡)
𝑊𝑝𝑢(𝑡) = = = ≈1
Wmsgn 2 Wsgn Wsgn
𝑃
En la práctica, la velocidad del rotor de la turbina de gas no varía significativamente en
relación a la velocidad sincrónica de generador.

Debido a esto: Wpu (t) ≈1 y Por la misma razón Pmpu (t) ≈1,

Al reemplazar estas condiciones en la ecuación se obtiene.

3𝐻 𝑑3 ȡ𝑚(𝑡)
∗ 𝑊𝑝𝑢(𝑡) = 𝑃𝑚𝑝𝑢(𝑡) − 𝑃𝑒𝑝𝑢(𝑡) = 𝑃𝑎𝑝𝑢(𝑡)
Wsgn 𝑑𝑡 3

Debido a la desviación, necesita la ecuación que se añade un término que representa un

generador de par de torsión en cualquier momento de amortiguación se desvía de su
velocidad síncrona, con su valor proporcional a la desviación de la velocidad.

3𝐻 𝑑 3 ȡ𝑚(𝑡) 𝐷 𝑑ȡ𝑚(𝑡)
∗ 𝑊𝑝𝑢(𝑡) 3
= 𝑃𝑚𝑝𝑢(𝑡) − 𝑃𝑒𝑝𝑢(𝑡) − = 𝑃𝑎𝑝𝑢(𝑡)
Wsgn 𝑑𝑡 𝑊𝑠𝑔𝑛 𝑑𝑡

Donde D es cero o un número relativamente pequeño positiva con valores típicos entre 0 y
2, esta ecuación se llama la ecuación de oscilación por unidad.

La ecuación es una ecuación diferencial de segundo orden que puede ser reescrita como dos
ecuaciones diferenciales de primer orden.

𝑑ȡ𝑚(𝑡)
= 𝑊(𝑡) − 𝑊𝑠𝑔𝑛
𝑑𝑡

3𝐻 𝑑𝑊(𝑡) 𝐷 𝑑ȡ𝑚(𝑡)
∗ 𝑊𝑝𝑢(𝑡) = 𝑃𝑚𝑝𝑢(𝑡) − 𝑃𝑒𝑝𝑢(𝑡) − = 𝑃𝑎𝑝𝑢(𝑡)
Wsgn 𝑑𝑡 𝑊𝑠𝑔𝑛 𝑑𝑡

𝐸𝑣
𝑃𝑒 = 𝑠𝑖𝑛ȡ
𝑋𝑒𝑞

Teniendo en cuenta que es no lineal debido a la Pepu (t) y Wepu (t)

 5. RESULTADOS Y DISCUSIONES

Figura 2: simulación en MATLAB_SIMULINK de una turbina de gas.

De acuerdo al diagrama de bloques simulado la entrada del sistema es un step o un rango de

datos lo que en la actualidad se asemeja a un tablero en el cual se introducirán datos o
parámetros de funcionamiento de acuerdo a las necesidades del usuario, y la salida como se
mencionó anteriormente es energía mecánica o revoluciones que puede producir la turbina
de gas.

El resultado más importante de este proyecto es obtener la simulación o diagrama de bloques

de la turbina de gas, ya que con el mismo se podrían realizar pruebas aplicativas de su
funcionamiento en diferentes circunstancias. Cabe recalcar que el estudio de donde se extrajo
la simulación es muy profundo ya que tiene la capacidad de implementar nuevos subsistemas
para analizar el comportamiento de la turbina de gas.

Para un mejor entendimiento del diagrama de bloques de la turbina de gas, se presenta el

diagrama de bloques de su función de transferencia así como también una tabla en la cual se
encuentran especificadas las constantes utilizadas en el diagrama.
 Figura 3: Diagrama de bloques de la función de transferencia de una turbina de gas
Fuente: http://liaowangbj.weebly.com/independent-study.html

Tabla 1: Tabla de constantes para una turbina de gas

PARÁMETRO DESCRIPCIÓN VALOR
W Ganancia=1/modo (pu MW/pu speed) 16.7
X Gobernador tiempo de espera constante (s) 0.6
Y Gobernador de tiempo de retraso constante(s) 1.0
Z Modo de gobernador (1=modo, 0=isócrono) 1
MAX Demanda límite superior (pu) 1.5
MIN Demanda limite menor (pu) -0.1
A Posicionador de válvula 1
B Posicionador de válvula 0.05
C Posicionador de válvula 1
WMIN Flujo de combustible mínimo 0.23
TF Combustible tiempo de control constante (s) 0.4
KF Retroalimentación del sistema de combustible 0
ECR Tiempo de reacción de combustión (s) 0.01
ETD Turbina y escape de retardo (s) 0.04
TCD Compresor de volumen de descarga constante (s) 0.2
TR Temperatura nominal de escape de la turbina ( ͦF) 950
TT Tasa de integración, controlador de temperatura ( ͦF) 450
f1 TX=TR-700*(1-WF)+550*(1-N)
f2 1.3*(WF-0.23)+0.5*(1-N)
TI Inercia= 2*H 15.64
Mediante la tabla y el diagrama de bloques de la función de transferencia se pueden
identificar claramente las constantes utilizadas en la simulación ya que solo es necesario
observar la función de transferencia la terminología de cualquier constante y buscarla en la
tabla. Como se expresó desde el principio la salida de la simulación de la turbina de gas es
su energía mecánica o revoluciones, en la siguiente grafica se puede observar una curva la
cual muestra la relación entre las revoluciones y el tiempo, como se puede apreciar la turbina
de gas en un minuto puede alcanzar alrededor de 1480rev.

Figura 4: revoluciones de la turbina de gas en 60s

Si el tiempo se extiende se pude apreciar que las revoluciones siguen aumentando hasta llegar
a un tope máximo de 1590rev en un tiempo de 300s.

Figura 5: revoluciones de la turbina de gas en 300s

Como se aprecia en las gráficas anteriores una turbina de gas tiene la capacidad de
incrementar sus revoluciones de una manera rápida y versátil, debido a parámetros tales como
su construcción, materiales y la baja inercia que posee. Por consiguiente una turbina de gas
es la alternativa óptima para ser colocada en un avión o un uso industrial en lugar de otras
fuentes de generación, ya que posee capacidades superiores de operación tales como: la
altitud de trabajo, rendimiento, velocidad interna de los mecanismos y durabilidad de las
piezas.

La turbina de gas como tal ofrece muchas ventajas frente a las fuentes de generación
convencionales es por esto que siguen siendo cada vez más desarrolladas y utilizadas como
fuentes de generación de electricidad e inclusive se ha empezado a utilizarlas en automóviles.

6. CONCLUSIONES

Se pudo entender el principio de funcionamiento que posee una turbina de gas así como
también todos los parámetros que deben ser involucrados para su funcionamiento y
simulación tales como, temperatura, cantidad de combustible entre otros.

Mediante la realización del presente proyecto se pudo poner en práctica ciertas

herramientas del SIMULINK que fueron impartidas por el docente en clases, así como
también se logró conocer herramientas nuevas que ofrece este entorno.

A pesar de que la información relacionada al tema del proyecto es bastante limitada en

nuestro idioma, con la ayuda de buscadores científicos se pudo encontrar una parte de la
información expuesta para luego entenderla y desarrollarla.

BIBLIOGRAFÍA

[1] LIAO WANG, Idependent study with Dr. Christy, 2014

Leído de: http://liaowangbj.weebly.com/independent-study.html

[2] W. ORDYS, A. PIKE, M.A JOHNSON, Modelling and simulation of power generation
plants, Springer-Verlag, 1994.

Leído en: http://link.springer.com/book/10.1007%2F978-1-4471-2114-5

[3] FRANCESC SALES TUDURÍ, Diseño y construcción de una Turbina de Gas, 2013, Pg
17.

Leido en:
http://upcommons.upc.edu/bitstream/handle/2099.1/16651/Dise%C3%B1o%20y%20Const
rucci%C3%B3n%20de%20una%20Turbina%20de%20Gas.pdf

[4] FERNANDEZ JORGE FELIX, Ciclo Brayton Turbinas a gas. Mendoza. Universidad
Tecnológica Nacional; Facultad Regional Mendoza, Mayo 2009.

Leído de: http://www.edutecne.utn.edu.ar

ACRÓNIMO

MCI: Motor de Combustión Interna.