DOI:

https://doi.org/10.14483/23448350.357

Publicado:

11/30/2006

Número:

Núm. 9 (2007): Enero-diciembre

Sección:

Ingeniería y Tecnología

Diseño e implementación de unidades de conversión de energía eléctrica para uso en redes de generación a partir de fuentes no convencionales

Design and implernentation of units of conversion of electrical energy for using in generation networks, from non-conventional sources

Autores/as

  • Cesar Leonardo Trujillo Rodríguez Universidad Distrital Francisco José de Caldas
  • Rafael Antonio Peña Suesca Universidad Distrital Francisco José de Caldas
  • Javier Antonio Guacaneme Moreno Universidad Distrital Francisco José de Caldas

Palabras clave:

Conversión de potencia bidireccional, fuentes no convencionales, modulación de ancho de pulso sinusoidal (es).

Palabras clave:

Bidirectional power conversion, non-conventional sources, sinusoidal pulse width modulation. (en).

Referencias

V. Gabriunas, Apuntes de electrónica", Bogotá, Universidad Distrital Francisco José de Caldas 2000, p. 355.

F. Marder, Generación de onda sinoidal. Modulación por ancho de pulso corriente, Department of Electric Engineering, Faculty of Exact and Natural Sciences and urveying, UNNE. Corrientes, 1999, p. 4.

M. Rashid, Electrónica de potencia, 3rd ed., México, McGraw-Hill, 1995, p. 702.

J. Gualda, S. Martínez y P. Martínez, Electrónica industrial: técnicas de potencia. 2nd ed., Madrid, Alfaomega Marcombo, 1992, p. 477.

Y. Zhang, "Snubber Considerations for IGBT Applications Power". International Rectifier. Application Note, International Versión. p. 7.

J. Grainger y W. Stevenson, Análisis de sistemas de potencia, México, McGraw-Hill, 1996, p. 740.

Cómo citar

APA

Trujillo Rodríguez, C. L., Peña Suesca, R. A., y Guacaneme Moreno, J. A. (2006). Diseño e implementación de unidades de conversión de energía eléctrica para uso en redes de generación a partir de fuentes no convencionales. Revista Científica, (9), 169–186. https://doi.org/10.14483/23448350.357

ACM

[1]
Trujillo Rodríguez, C.L. et al. 2006. Diseño e implementación de unidades de conversión de energía eléctrica para uso en redes de generación a partir de fuentes no convencionales. Revista Científica. 9 (nov. 2006), 169–186. DOI:https://doi.org/10.14483/23448350.357.

ACS

(1)
Trujillo Rodríguez, C. L.; Peña Suesca, R. A.; Guacaneme Moreno, J. A. Diseño e implementación de unidades de conversión de energía eléctrica para uso en redes de generación a partir de fuentes no convencionales. Rev. Cient. 2006, 169-186.

ABNT

TRUJILLO RODRÍGUEZ, Cesar Leonardo; PEÑA SUESCA, Rafael Antonio; GUACANEME MORENO, Javier Antonio. Diseño e implementación de unidades de conversión de energía eléctrica para uso en redes de generación a partir de fuentes no convencionales. Revista Científica, [S. l.], n. 9, p. 169–186, 2006. DOI: 10.14483/23448350.357. Disponível em: https://revistas.udistrital.edu.co/index.php/revcie/article/view/357. Acesso em: 16 abr. 2024.

Chicago

Trujillo Rodríguez, Cesar Leonardo, Rafael Antonio Peña Suesca, y Javier Antonio Guacaneme Moreno. 2006. «Diseño e implementación de unidades de conversión de energía eléctrica para uso en redes de generación a partir de fuentes no convencionales». Revista Científica, n.º 9 (noviembre):169-86. https://doi.org/10.14483/23448350.357.

Harvard

Trujillo Rodríguez, C. L., Peña Suesca, R. A. y Guacaneme Moreno, J. A. (2006) «Diseño e implementación de unidades de conversión de energía eléctrica para uso en redes de generación a partir de fuentes no convencionales», Revista Científica, (9), pp. 169–186. doi: 10.14483/23448350.357.

IEEE

[1]
C. L. Trujillo Rodríguez, R. A. Peña Suesca, y J. A. Guacaneme Moreno, «Diseño e implementación de unidades de conversión de energía eléctrica para uso en redes de generación a partir de fuentes no convencionales», Rev. Cient., n.º 9, pp. 169–186, nov. 2006.

MLA

Trujillo Rodríguez, Cesar Leonardo, et al. «Diseño e implementación de unidades de conversión de energía eléctrica para uso en redes de generación a partir de fuentes no convencionales». Revista Científica, n.º 9, noviembre de 2006, pp. 169-86, doi:10.14483/23448350.357.

Turabian

Trujillo Rodríguez, Cesar Leonardo, Rafael Antonio Peña Suesca, y Javier Antonio Guacaneme Moreno. «Diseño e implementación de unidades de conversión de energía eléctrica para uso en redes de generación a partir de fuentes no convencionales». Revista Científica, no. 9 (noviembre 30, 2006): 169–186. Accedido abril 16, 2024. https://revistas.udistrital.edu.co/index.php/revcie/article/view/357.

Vancouver

1.
Trujillo Rodríguez CL, Peña Suesca RA, Guacaneme Moreno JA. Diseño e implementación de unidades de conversión de energía eléctrica para uso en redes de generación a partir de fuentes no convencionales. Rev. Cient. [Internet]. 30 de noviembre de 2006 [citado 16 de abril de 2024];(9):169-86. Disponible en: https://revistas.udistrital.edu.co/index.php/revcie/article/view/357

Descargar cita

Visitas

544

Dimensions


PlumX


Descargas

Los datos de descargas todavía no están disponibles.

Ingeniería y Tecnología

Revista Científica, 2007-08-00 nro:9 pág:169-186

Diseño e implementación de unidades de conversión de energía eléctrica para uso en redes de generación a partir de fuentes no convencionales

Design and implernentation of units of conversion of electrical energy for using in generation networks,from non-conventional sources

Cesar Leonardo Trujillo

Ingeniero electrónico Universidad Distrital Francisco José de Caldas. Magíster en Ingeniería EléctricaUniversidad Nacional de Colombia. Profesor Universidad Distrital Francisco José de Caldas. Universidad Distrital Francisco José de Caldas — Línea de Investigación en Fuentes Alternas de Energía LIFAE cltrujillo@udistrital.edu.co

Rodríguez Rafael Antonio

Ingeniero electrónico Universidad Distrital Francisco José de Caldas. Magíster en Ingeniería Eléctrica Universidad Nacional de Colombia. Profesor Universidad Distrital Fr;ncisco José de Caldas. rpena@udistrital.edu.co

Peña Suesca Javier Antonio Guacaneme Moreno

Ingeniero electrónico Universidad Distrital Francisco José de Caldas. Magíster en Ingeniería Eléctrica Universidad Nacional de Colombia. Profesor Universidad Distrital Francisco José de Caldas. jguacaneme@udistrital.edu.co

Resumen

En muchos lugares del mundo existen zonas adonde es realmente complejo transportar la energía eléctrica. Para estos casos, una solución es proveerse de sistemas propios de generación, entre los cuales puede considerarse el aprovechamiento de fuentes de energía, posiblemente no convencionales, presentes en estas zonas. Una alternativa es el diseño de un sistema que transforme y almacene energía en un banco de baterías y, a partir de este, convertirla en valores estándar A, C (por ejemplo, 120 V/60 Hz). Dichos procesos de conversión deben ser muy eficientes, porque en general se trata de energía disponible durante ciertos periodos, en muchos casos escasa y costosa.

En el presente artículo se analiza la etapa de conversión a partir de un bus DC el cual provee un medio de almacenamiento de la energía que se recolecta a través de algún tipo de fuente. Así mismo se analiza lo relacionado con los principios de operación, dimensionamiento de los componentes y características de diseño para la implementación de unidades electrónicas de conversión DC/AC de baja potencia, capaces de manejar potencia activa y reactiva en ambas direcciones (desde las baterías hacia la red y viceversa) y con posibilidad de interconexión a otros sistemas con generación. Además se presenta un análisis práctico del desempeño de la unidad electrónica.

Palabras clave
Conversión de potencia bidireccional, fuentes no convencionales, modulación de ancho de pulso sinusoidal.

Abstract

In many places of the world there are zones where it is really complex to transport the electrical energy. In these cases, a solution is to provide itseif with own systems of generation in which the advantage of non-conventional power sources presented in these zones can be considered. An alternative is the design of a system that transforms and stores energy in a bank of batteries and, from this one, to convert it to values standard AC (e.g. 120V/60Hz). These conversions must be highly efficient, since generally it is energy available during certain periods, in many cases little and expensive one.

In this paper is analyzed the convertion stage from a DC bus which provides a storage media of the energy that is collected through some type of source. They are also analyzed the operation principies, sizing of the components and characteristics of design for the implementation of DC/AC electronic units of conversion of low power able to handle active power and reactive power in both directions (from the batteries towards the network and vice versa) and with possibility of interconnection to others generation systems. Moreover, a practical analysis of the performance of the electronic unit is described.

Key words:
Bidirectional power conversion, non-conventional sources, sinusoidal pulse width modulation.


INTRODUCCIÓN

Debido al incremento en la demanda de energía eléctrica, a que las fuentes tradicionales de generación no son suficientes y que existen regiones desatendidas por este servicio, ha sido necesario investigar en nuevas fuentes de energía, como el viento y el sol.

Debido a que la energía eléctrica presenta ventajas como fácil transporte, distribución y conversión en otros tipos de energía, la energía colectada de las fuentes no convencionales es transformada en energía eléctrica DC que presenta ventajas para la interconexión y, en algunos casos, es almacenada en baterías para uso posterior. Es necesario disponer de sistemas que conviertan esta energía DC en energía AC; que para el caso de bajas potencias, se requiere a 208/120 Vr./60 HZ. Aparte de realizar dicha transformación estos sistemas, deben ser altamente eficientes y capaces de transferir energía hacia la red de energía alterna, dentro del principio de cogeneración.

Un sistema de generación y distribución de energía a partir de fuentes no convencionales tiene la característica de utilizar fuentes que no están disponibles las 24 horas del día, como en los casos del viento' y del sol. Lo anterior es ventajoso cuando el sistema está interconéctado con una red de generación convencional que suple los requerimientos de energía en los momentos que no se tiene generación de las fuentes no convencionales. Los sistemas de generación no convencional aislados, que deban garantizar servicio continuo, requieren sistemas adicionales de generación o sistemas de almacenamiento de energía; por ejemplo, bancos de baterías.

Teniendo en cuenta lo anterior y las ventajas de las interconexiones DC, estos sistemas requieren conversores AC/DC y DC/AC. A continuación se analizará una unidad de conversión prototipo pequeña que permite la transferencia de energía de las fuentes DC a la red AC, y viceversa.

UNIDAD DE CONVERSIÓN

Conversión DC/AC
Se implementó mediante configuración inversora de puente completo [1], como se muestra en la figura 1.

La configuración anterior permite obtener sobre la carga voltajes Vo, positivos, negativos o cero, según lo determine el circuito de control de los interruptores.

Funcionamiento bipolar

Si en el circuito de la figura 1, en el tiempo Ti están cerrados los interruptores Si y S4, y abiertos S2 y S3, entonces Vab = E. Si posteriormente en el tiempo T2 se cierran los interruptores S2 y S3, y se abren Si y S4, entonces Vab = -E. La forma de onda de Vab se muestra en la figura 2.

Para señales periódicas con T = T1+T2, el valor medio de Vab es

    (1)

Definiendo T 1/T como el ciclo útil (δ) de la señal, VabDC se puede expresar como

 :   (2)

Teniendo en cuenta que L es cortocircuito para DC, entonces

    (3)

es decir,

    (4)

Si el filtro pasabajos LC elimina las componentes AC, entonces el voltaje Vo es solo la componente DC. En la tabla 1 se muestran los valores de Vo para diferentes valores de δ.

Modulación SPWM

En la práctica, la conmutación de los interruptores se realiza a alta frecuencia (de 1 kHz a 20 kHz) y el filtro LC pasabajos se escoge para dejar pasar bajas frecuencias (50 Hz o 60 Hz) y bloquear las altas frecuencias. En estos casos el voltaje Vo puede variar a baja frecuencia y la ecuación (4) se puede escribir

(5)

Para obtener un voltaje sinusoidal de la forma

(6)

se requiere

(7)

que se obtiene despejando δ ( w t) de las ecuaciones (5) y (6).

Como se observa en la ecuación (7), el ciclo útil (ancho de pulso) varía de acuerdo con una función sinusoidal. Esta técnica se denomina SPWM (Sinusoidal Pulse Width Modulation) [2]. Es decir se modula el ancho del pulso de una portadora de alta frecuencia (20 kHz) con una señal sinusoidal moduladora de baja frecuencia (60 Hz).

Para la implementación de un sistema SPWM se debe tener en cuenta:

• Las frecuencias de la portadora y de la moduladora.

• Número de pulsos por ciclo de moduladora.

• Indice de modulación m = 2* δ m, 0<m<1.

• Sistema de filtrado.

• Dispositivos utilizados para la implementación.

En el prototipo que se implementó, la generación de la señal SPWM se realizó digitalmente usando un microcontrolador. El proceso consistió en dividir el periodo de la señal sinusoidal en p intervalos. Se discretizó la señal sinusoidal en 340 particiones para poder describirla como una sucesión de anchos de pulso variables en el tiempo.

Descripción del prototipo

La unidad de conversión prototipo es monofásica, de baja potencia (400 W) y presenta las siguientes características:

• Conversión 48 VDC/120 Vrms sinusoidal regulada a 60 Hz, y viceversa.

• Control de corriente AC constante cuando se presenta sobrecarga.

•Posibilidad de interconectarse con redes AC convencionales y con otros módulos semejantes.

Dimensionamiento de los dispositivos

Para la etapa de potencia de la figura 1, que funciona con modulación SPWM, los voltajes y corrientes en los dispositivos S1, S2, S3 y S4, idealmente, tienen las formas de onda que se muestran normalizadas en la figura 3a y 3b.

En las ecuaciones (8) se presentan las expresiones de voltaje y corriente sobre los dispositivos Si, S2, S3 y S4 [3].

E representa el voltaje de las baterías, Ip la corriente pico sobre el dispositivo y 340 es el número de particiones tomadas en un periodo.

La unidad prototipo funciona a partir de un banco de baterías de 48 VDC y utiliza un transformador de relación 1 a 6.5 que suministra una potencia de 400 W a una tensión de 120 Vrms. En los dispositivos se obtuvieron los valores de voltajes y corrientes que se muestran en la tabla 2.

Para la topología utilizada, con bajas tensiones y mayores corrientes, las pérdidas en los dispositivos aumentan y la eficiencia se reduce aproximadamente al 70% [4].

Los dispositivos se conmutan a frecuencias altas haciendo más evidente el efecto de las capacitancias e inductancias, originando sobretensiones. Por tanto, y para efectos de diseño, es conveniente escoger dispositivos que so porten picos varias veces superiores a la tensión de las baterías.

Filtro de salida

La señal SPWM de frecuencia portadora (fc) de 20 kHz y frecuencia moduladora (fm) de 60 Hz tiene el e spécíro que se muestra én la figura 4. La energía se concentra principalmente en 60 Hz, y sus armónicos en 20 kHz, 40 kHz y 60 kHz. De acuerdo con lo anterior, se escogió' Para el diseño del filtro una frecuencia de corte de aproximadaniente 180 Hz. Una simulación de la característica se observa en la figura 5.

La topología del filtro es LC, al cual se le han incluido la resistencia de pérdidas de la bobina (RB) y la resistencia de carga RL como se muestra en la figura 6.

El cálculo se hizo estableciendo su función de transferencia, igualándola a 1/√2 (valor de media potencia) y calculando para s= jw). La función de transferencia se presenta en la ecuación (9):

y se obtuvo el valor de L = 601 µH para RB = 0.06Ω, C = 25 µF, f= 200 Hz y RL =0.73Ω para Vo. 24 Vrms.

Transformador de salida

Sin considerar las pérdidas, idealmente el voltaje que puede generar el puente alimentado con el voltaje nominal de las baterías es 48/√2 Vrms. Considerando las pérdidas de voltaje en las baterías, el conversor DC/AC y el filtro, el voltaje a la salida de la unidad es 24 Vrms. Teniendo en cuenta que se requiere un transformador para elevar el voltaje a 120 V. y que para efectos de control es conveniente trabajar con índices de modulación inferiores a uno. Se escogió un transformador con relación de transformación de 6.5.

Protecciones

Debido a los sobrepicos de tensión por conmutación que se presentan en los dispositivos semiconductores, se requieren snubbers que limiten los voltajes máximo (Vpk) [5] (figura 7).

Snubber Cd: corresponde al condensador colocado lo más cerca posible y en paralelo a la rama de los interruptores. Se calcula mediante la ecuación (10):

Cd =(Ls +Lc + LB) I2o /(Vpk — Vcc)2

•donde Ls: es la inductancia del dispositivo, Lc: es la inductancia de los conductores y LB es la inductancia de la bobina.

• Snubber de carga y descarga: se compone de un circuito RC y un diodo. El condensador y la resistencia se pueden calcular a partir de la ecuación (11):

donde LT es la inductancia parásita del circuito cuando el interruptor del lado inferior se encuentra encendido; y fs es la frecuencia de conmutación.

Conversión AC/DC

Para explicar el comportamiento de la unidad, es conveniente modelar el sistema como se muestra en la figura 8, donde VTH y ZTH es el equivalente Thevenin de la unidad de conversión, V y Z es el equivalente de la fuente externa. Para facilidad del análisis, se asume que tanto VTH como V son voltajes sinusoidales de 60 Hz.

El flujo de potencia activa y potencia reactiva entre V y VTH depende de los valores de voltaje y fase de las fuentes y de los valores de la impedancia ZTH + Para el caso de impedancia puramente reactiva, las potencias están dadas por la expresión (12), donde la transferencia de potencia activa depende principalmente de las diferencias angulares, La transferencia de potencia reactiva depende de la diferencia de las magnitudes de voltaje [6].

En el caso de impedancias puramente resistivas, la potencia activa depende principalmente de las diferencias de voltaje; la potencia reactiva, de las diferencias angulares. En los casos intermedios, en donde la parte resistiva y reactiva es igualmente importante, las potencias activas y reactivas dependen de las diferencias de voltaje y de las diferencias angulares.

Para el prototipo implementado se realizaron ensayos de transferencia de potencia desde la red a las baterías, variando en primer lugar la magnitud de VTH, en fase con V, y variando en segundo lugar la fase de VTH, con magnitud constante.

Las variaciones de magnitud y fase de VTH se realizaron modificando el índice de modulación m y el desfasamiento φ de la expresión (13) sobre la que se puede realizar el control.

Vth = m•a •E• sen (wt+φ)

a: relación del transformador

En las figura 9 y 10 se muestran los resultados de los ensayos realizados sobre la unidad.

En la figura 9 se observa cine en las corrientes máximas para las cuales corresponde una potencia máxima de 168 W, aproximadamente, el equipo presenta una eficiencia del 50%

De la figura 10 se observa que en las corrientes máximas para las cuales corresponde una potencia máxima de 324 W, aproximadamente, el equipo presenta una eficiencia del 44%.

De lo anterior se deduce que se puede transferir más potencia desde la red a las baterías variando el ángulo de desfasamieno φ.

Modos de operación
La unidad de conversión tiene básicamente tres modos de operación.

Modo voltaje AC constante

En este modo, la unidad de conversión mantiene el voltaje en 120 Vrms, siempre y cuando la potencia exigida no sea mayor de 400 W. En caso que el equipo esté sometido a sobrecarga, pasa automáticamente al modo corriente AC constante hasta cuando desaparezca la sobrecarga.

Modo corriente AC constante

Se presenta automáticamente en la unidad en caso de sobrecarga; es decir, esta función actúa como un limitador de corriente. Cuando la unidad de conversión se encuentra en modo corriente AC constante suministra una corriente cercana a los 4 Arms.

Modo red

Es el modo de operación de interconexión con la red pública. Está en capacidad de entregar y de recibir energía. También se puede colocar en estado de flotación, con flujo de energía nulo.

La transferencia de energía se logra variando de manera controlada parámetros propios de la unidad, como amplitud y fase, lo que lleva a que se presenten tres estados de transferencia de energía, los cuales se citan a continuación:

• Estado de flotación. Cuando el índice de modulación y el ángulo de desfasamiento φ son tales que no se presenta intercambio de energía. Para este caso, el voltaje de salida de la unidad es igual al voltaje de la red en magnitud y fase.

• Estado de generación. En este modo, la energía fluye desde las baterías hacia la red. El ángulo de desfasamiento φ de la unidad se adelanta con respecto al de la red entregando potacia activa.

• Estado de carga. En este modo de óperación, la energía fluye desde la red hacia las baterías. El ángulo de desfasamiento φ de lá unidad se atrasa con respecto al de la red recibiendo potencia activa para la carga de las baterías.

Operación de unidades en paralelo

En sistemas aislados que requieran potencias inayoíes a la capacidad nominal de lá unidad de conversión, se pueden interconectar varias unidades en paralelo. La sincronización del sistema se puede realizar con un . equipo Maestro y uno o varios esclavos con capacidad de generar salidas en fase con la salida del equipo Maestro.

CARACTERÍSTICAS DE DESEMPEÑO DE LA UNIDAD

Característica en vacío del conversor para diferentes valores de alimentación DC

Esta prueba se llevó a cabo aplicando diferentes índices de modulación para tres valores de batería (E= 36.4 VDC, 47 VDC y 61 VDC). Se midieron los voltajes de salida sin carga, obteniendo los resultados que se muestran en la figura 11.

A partir de la gráfica se concluye que el efecto de variar E, se traduce en un cambio de pendiente en la característica de vacío de la unidad, que era de esperarse, porque el voltaje de salida es directamente proporcional a E.

Características en carga e impedancia equivalente de la unidad
Para determinar la magnitud de las impedancias del sistema, se procedió a elaborar una serie de mediciones de voltaje versus. corriente de carga para diferentes índices de modulación.

Los puntos sobre los cuales se tomaron medidas de voltaje fueron la salida del puente, la salida del filtro y la salida del transformador. Estos voltajes se midieron para diferentes corrientes de salida hasta llegar a la corriente nominal. A continuación se muestra la gráfica para m = 68.4%.

En la figura 12,1a magnitud de la impedancia del equipo es aproximadamente 17.6 Ω referidos al secundario del transformador.

Eficiencia de la unidad

Con carga resistiva pura de 400W y corriente de 3.33 A nns, la potencia de entrada Pi es

Pi = Po + I2o * R = 400 + (3.3332 *17.6) = 595.16W

(14)

Finalmente, la eficiencia es

(15)

Parámetros de la unidad de conversión

A continuación se presentan los principales parámetros de la unidad de conversión, así como la captura en el osciloscopio de la forma de onda de la tensión de salida, atenuada por 10.

Como se observa en la figura 13,1a tensión de salida es uña señal de 120VAC a 60Hz, con baja distorsión armónica.

CONCLUSIONES

Los conversores DC/AC con SPWM facilitan la interconexión de sistemas de generación a partir de fuentes de energía no convencionales a redes alimentadas con sistemas de generación tradicionales.

En redes aisladas, el servicio se puede mejorar, en horas de baja generación, con sistemas de almacenamiento de energía que en algunos casos pueden ser bancos de baterías.

La tecnología SPWM presenta la gran ventaja de generar ondas sinusoidales con un bajo índice de distorsión armónica, y puede implementarse de manera rápida y sencilla. El prototipo construido puede conectarse directamente cuando la red está energizada, debido a que posee circuitos de control que permiten sincronizarse de manera rápida y segura.

En desarrollos futuros, es conveniente trabajar con conversores trifásicos y aplicar otros tipos de tecnologías de modulación a fin de lograr mayores ventajas en aplicaciones relacionadas con redes de generación distribuida y redes de fuentes de energía no convencional.

BIBLIOGRAFÍA

  1. V. Gabriunas, Apuntes de electrónica, Bogotá, Universidad Distrital Francisco José de Caldas 2000, p. 355.
  2. F. Marder, Generación de onda sinoidal. Modulación por ancho de pulso corriente, Department of Electric Engineering, Faculty of Exact and Natural Sciences and Surveying, UNNE. Corrientes, 1999, p. 4.
  3. M. Rashid, Electrónica de potencia, 3rd ed., México, McGraw-Hill,1995, p. 702.
  4. J. Gualda, S. Martínez y P. Martínez, Electrónica industrial: técnicas de potencia. 2nd ed., Madrid, Alfaomega Marcombo, 1992, p. 477.
  5. Y. Zhang, "Snubber Considerations for IGBT Applications Power". International Rectifier. Application Note, International Versión. p. 7.
  6. J. Grainger y W. Stevenson, Análisis de sistemas de potencia, México,McGraw-Hill, 1996, p. 740.

Creation date:
Loading...