c9 Apendice Atp-emtp

Embed Size (px)

Citation preview

  • 7/28/2019 c9 Apendice Atp-emtp

    1/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez357

    Apndice A

    ATP-EMTP

    A.1 Qu es ATP?.

    El Programa de Transitorios Alternativo (Alternative Transients Program, ATP) esconsiderado un programa informtico universal para la simulacin de fenmenoselectromagnticos transitorios y de naturaleza electromecnica en sistemas elctricos de potencia.Complejas redes elctricas y sistemas de control pueden ser de esta manera simulados.

    El Programa de Transitorios Electromagntico (Electromagnetic Transients Program, EMTP)fue desarrollado en el dominio pblico de la Administracin de Energa Bonneville (Bonneville Power

    Administration, BPA) de Portland, Oregon, como paso previo a su comercializacin en 1984 por elGrupo de Coordinacin de Desarrollo del EMTP (EMTP Developmet Coordination Group) y elInstituto de Investigacin de Energa Elctrica (Electric Power Research Institute, EPRI) de Palo Alto,California. Varios expertos en todo el mundo, sin embargo, han contribuido a su posteriordesarrollo y mejora en aos posteriores.

    Protegido por las leyes de Estados Unidos, ATP puede ser adquirido bajo licencia,expedida y garantizada por los propietarios del material.

    A.2 Requerimientos de hardware y software para ATP.

    ATP est disponible para la mayora de plataformas PC basadas en Intel: DOS, Windows3.1/9x/NT, OS/2, Linux, etc., y para otras computadoras (Digital Unix y VMS, Apple Macs, etc.)

    Un Pentium PC con una configuracin mnima de 128 Mb de memoria RAM, 20 Mb deespacio libre en disco duro y VGA grficos es suficiente para ejecutar ATP bajo MS-Windows.

  • 7/28/2019 c9 Apendice Atp-emtp

    2/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez358

    Las versiones ms populares de compiladores son las siguientes:

    MS-Windows 9x/NT/2000/XP:

    32 bit GNU-Mingw32 y Watcom ATP.

    MS-DOS, MS-Windows 3x/95/98:

    32 bit Salford ATP (requiere DBOS/486).

    Linux:

    Versin GNU de ATP.

    A.3 Licencia de ATP.

    Para bajar el programa de internet ha de obtenerse el permiso o licencia del Grupo deUsuarios de EMTP Canad/Amrica (Canadian/American EMTP User Group) o un grupo deusuarios regional autorizado. Mayor informacin al respecto puede encontrarse en la pgina web

    www.emtp.org .

    A continuacin se muestra la localizacin geogrfica de los distintos grupos de usuarios deATP-EMTP y de detalles de contacto:

    Figura. Localizacin de grupos de usuarios de ATP-EMTP.

  • 7/28/2019 c9 Apendice Atp-emtp

    3/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez359

    Nombre E-mail y direccin wwwCanadian/American EMTP User Group [email protected] EMTP-ATP Users Group(EEUG)

    [email protected], www.eeug.org

    Japanese ATP User Group (JAUG) [email protected],www.jaug.jp/~atp/index-e.htmLatin American EMTP User Group(CLAUE)

    [email protected],www.furnas.gov.br/atp

    Argentinian EMTP User Group (CAUE) [email protected]://iitree.ing.unlp.edu.ar/estudios/caue/caue.html

    Australian EMTP User Group (AEUG) [email protected] EMTP User Group (KEUG) [email protected] of China EMTP User Group [email protected] EMTP User Group (IEUG) [email protected] African ATP User Group(SAAUG)

    [email protected],www.ee.wits.ac.za/~atp

    Tabla. Direcciones email y pginas web de grupos de usuarios de ATP-EMTP.

    A.4 Capacidades del programa.

    La siguiente tabla muestra los lmites mximos de algunos componentes que se puedenutilizar en la simulacin de un circuito de potencia:

    Lneas 6000Elementos lineales 10000Interruptores 1200Fuentes 900Elementos no lineales 2250Mquinas sncronas 90

    A.5 Descripcin del entorno de trabajo.

    A fin de facilitar su uso, ATP-EMTP permite la integracin de una serie de programascomplementarios: editores de texto, un editor grfico de circuitos elctricos y aplicacionesespecficas para la visualizacin grfica de los resultados obtenidos. Se accede a dichos programas atravs del gestor del entorno, ATPCC.

  • 7/28/2019 c9 Apendice Atp-emtp

    4/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez360

    De esta manera, la simulacin de un sistema elctrico o electromecnico puede iniciarsecon la creacin del modelo en el editor grfico (fichero de extensin .adp o .cir). La ejecucin (runATP) del modelo en el mismo editor genera el fichero .atp y los ficheros de salida .lis y.pl4, stosltimos registran los resultados de la simulacin efectuada. En un editor grfico puede visualizarsea partir del fichero de extensin .pl4 los resultados obtenidos de esta manera.

    Figura. Interaccin entre los programas de ATP-EMTP.

    A.5.1 Gestor del entorno: ATPCC.Como se ha mencionado, ATP Control Center (ATPCC) facilita el acceso al conjunto de

    programas que pueden formar parte del paquete ATP-EMTP (el editor grfico ATPDraw, eleditor de textos PFE32, los programas de representacin de resultados PCPlot, PlotXY yGTPPlot) o ejecutar directamente el propio ATP.

    Se ha de hacer notar que permite el acceso o integracin de hasta diez programasadicionales, as como la posibilidad de trabajar con dos versiones de ATP.

  • 7/28/2019 c9 Apendice Atp-emtp

    5/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez361

    A.5.1.1 Pantalla principal.

    Al ejecutar ATPCC aparece la pantalla principal del programa. Se muestra en la figurasiguiente:

    Figura. Pantalla principal del programa ATPCC.

    A.5.1.2 Configuracin de los programas principales.

    La configuracin de los programas principales (ATP, ATPDraw, PCPlot y el Editor) serealiza a travs de la opcin Main Programs Settings del submen Options.

    Se muestra de esta manera la siguiente ventana de dilogo, compuesta por dos pestaas. Laprimera de las cuales hace referencia a las dos versiones del programa ATP que posteriormente sepodrn ejecutar desde ATPCC:

  • 7/28/2019 c9 Apendice Atp-emtp

    6/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez362

    Figura. Ventana de dilogo para configuracin de programas principales: ATP.

    Las opciones por defecto indican los siguientes aspectos:

    Save on: salvar los resultados obtenidos en un fichero de extensin .lis (disk) o

    presentarlos tambin en la pantalla (both).

    Output Parameters:

    *:se marca esta casilla cuando se desea que el nombre del fichero desalida coincida con el nombre del fichero .atp utilizado como entrada. En casocontrario se especifica el nombre del fichero de salida (incluyendo su extensin .liso .out).

  • 7/28/2019 c9 Apendice Atp-emtp

    7/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez363

    -R:se selecciona esta opcin si se quiere que, cuando ya exista un fichero .liscorrespondiente al fichero .atp utilizado como entrada, el resultado de una nuevasimulacin se sobreescriba en ese mismo fichero.

    En la pestaa correspondiente al resto de programas principales se acta de un modo similar:

    Figura. Ventana de dilogo para la configuracin del resto de los programas principales.

    A.5.1.3 Configuracin de programas adicionales.

    Para aadir, modificar o eliminar programas adicionales se selecciona, en la lnea de mens,la opcin Options Additional Programs. Al hacerlo, se presenta una ventana de dilogocompuesta de tres pestaas.

  • 7/28/2019 c9 Apendice Atp-emtp

    8/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez364

    Figura. Incorporacin de programas adicionales.

    Las siguientes opciones adicionales equivalen a:

    Cut extension of input file: elimina la extensin de los ficheros de entrada, antes

    de ejecutarlos con el programa adicional. Default PLOT: se utiliza cuando se pretende que el programa adicional aadido

    (por ejemplo: GTPPLOT o PlotXY) acte como programa grfico por defecto. Eneste caso, la extensin de los ficheros de entrada debe estar definida como .pl4.

    A.5.1.4 Configuracin de parmetros generales.

    Desde la opcin Settings del submen Options pueden considerarse una serie deparmetros de carcter general.

  • 7/28/2019 c9 Apendice Atp-emtp

    9/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez365

    Figura. Configuracin de parmetros generales.

    ATP Control Center Working Directory: permite definir el directorio de trabajo delprograma ATPCC.

    Log file. Permite definir el nmero de eventos a recordar en la ventana de ltimoscomandos utilizados.

    Save last window position. Permite guardar la posicin y el tamao de la ventana delprograma para la prxima ocasin en que se utilice.

    A.6 Compilador: ATP.

    Este programa constituye el ncleo fundamental de la aplicacin y se trata del compiladorque permite procesar los datos del circuito objeto de estudio. Maneja ficheros de tipo de texto deextensin .dat o .atp.

    Las dos versiones de compiladores ATP ms habituales son:

    Watcom ATP. Admite nombres de ficheros de gran longitud, ficheros de datosde ms de 150.000 lneas y puede ejecutar simultneamente varios casos.

    GNU ATP. Permite utilizar nombres largos para ficheros, proporciona unarranque rpido y requiere menos memoria.

    En el esquema de la siguiente figura se muestran los mdulos de simulacin disponibles en ATP,las subrutinas de apoyo y la interaccin entre ellos. Con todos estos elementos se pueden crear

    modelos de menor a mayor complejidad, que representan los diferentes equipos y componentesencontrados en los sistemas elctricos.

  • 7/28/2019 c9 Apendice Atp-emtp

    10/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez366

    Figura. Componentes de ATP.

    A.7Editor de textos: PFE32.

    El desarrollo del procesador grfico ATPDraw ha simplificado enormemente la creacinde los ficheros que son compilados posteriormente con ATP. Sin embargo, para los usuariosexpertos sigue siendo necesario el empleo de un editor de textos que permita trabajar, entre otros,con los ficheros de entrada a ATP (.atp o .dat) o de salida del mismo (.lis).

    Conviene sealar que la estructura de los ficheros .atp es muy rgida, ya que toda lainformacin del sistema elctrico a simular debe ocupar una posicin especfica en cada fila delfichero. Para facilitar la escritura de los datos en el fichero se incluye una serie de plantillas(template). Adems, resulta indispensable tener a mano el manual de usuarioATP-EMTP RuleBook

    .

  • 7/28/2019 c9 Apendice Atp-emtp

    11/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez367

    A.7.1 Estructura general de un fichero .atp o .dat.

    El fichero comienza con una primera lnea denominada BEGIN NEW DATA CASE. Acontinuacin le siguen dos lneas MISCELLANEUS DATA CARDS, donde se permite especificarel paso de integracin, tiempo mximo de la simulacin, cada cunto tiempo almacena puntos parasu posterior representacin, etc.

    Tal y como se ha comentado, puede adjuntarse una plantilla al editor de textos PFE32.Ello puede hacerse mediante la opcin Template Attach File. Una vez adjuntada, se puedeutilizar sin ms insertndola en el fichero sobre el que se est trabajando, con la opcin Template Edit.

    Figura. Insercin de plantillas.

    Precediendo a las lneas MISCELLANEUS pueden utilizarse una serie de rdenes quepermiten la realizacin de clculos especiales, como, por ejemplo: flujos de carga, clculo dearmnicos, etc. stas pueden ser:

    - CABLE CONSTANT- DATA BASE MODULE- FREQUENCY SCAN- HARMONIC FREQUENCY SCAN

    - HYSTERESIS- LINE CONSTANT

  • 7/28/2019 c9 Apendice Atp-emtp

    12/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez368

    - SATURATION- XFORMER- ZNO FITTER- FIX SOURCE (esta opcin permite calcular un flujo de cargas)

    Una fila de comentarios empieza con la letra C.

    Por otra parte, en cualquier lnea del programa pueden especificarse las siguientes tarjetasespeciales:

    $OPEN $CLOSE $SWIDTH $DISABLE

    $ENABLE $NEW EPSLIN, EPS $MONITOR $LISTOFF$LISTON $VINATGE, M $STOP $WATCH$COMMENT $DEBUG $UNITS $UNITS$PUNCH $ERASE $SETUP, file name $INCLUDE

    La estructura general de un fichero .dat (o de un fichero .atp) es la siguiente:

    BEGIN NEW DATA CASE

    Miscellaneus Data Cards (2 lneas)

    /TACS (si existen) /MODELS (si existen)

    (Sistemas de control)

    BLANK TACS ENDMODELS/BRANCH

    (R, L, C, Z en general, transformadores, lneas)

    /SWITCH

    (Interruptores)

    /SOURCE

    (Fuentes de V/I ideales: rampa, escaln, sinusoidal,

    mquina sncrona, motores)

    /OUTPUT

    (Variables de salida a representar y que se incluirn

    en el fichero .pl4)

    BLANK BRANCHBLANK SWITCH

    BLANK SOURCE

    BLANK OUTPUT

    BEGIN NEW DATA CASE

    BLANK

    A continuacin se muestra un ejemplo, el cual ha sido generado por el editor grfico ATPDraw delcircuito elctrico que se muestra a continuacin:

  • 7/28/2019 c9 Apendice Atp-emtp

    13/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez369

    BEGIN NEW DATA CASE

    C --------------------------------------------------------

    C Generated by ATPDRAW marzo, jueves 11, 2010

    C A Bonneville Power Administration program

    C Programmed by H. K. Hidalen at SEfAS - NORWAY 1994-98

    C --------------------------------------------------------

    C Example 1

    C Your first circuit

    C Rectifier bridge

    C Miscellaneous Data Card ....

    C dT >< Tmax >< Xopt >< Copt >

    5.E-5 .05

    500 1 1 1 1 0 0 1 0

    C 1 2 3 4 5 6 7 8

    C 345678901234567890123456789012345678901234567890123456789012345678901234567890

    /BRANCH

    C < n 1>< n 2>< R >< L >< C >

    C < n 1>< n 2>< R >< A >< B >0

    VA XX0002 33. 1. 0

    XX0002 33. 1. 0

    XX0021VA 33. 1. 0

    XX0021 33. 1. 0

    XX0002POS .01 1

    POS XX0021 1.E3 3

    XX0021POS 20. 3

    VS XX0025 1. 0

    VS XX0025 300. 0

    /SWITCH

    C < n 1>< n 2>< Tclose >< Ie >< type >

    11VA XX0002 0

    11 XX0002 0

    11XX0021VA 0

    11XX0021 0

    XX0025VA MEASURING 1

    /SOURCE

    C < n 1>< Ampl. >< Freq. >< A1 >< T1 >< TSTART >< TSTOP >

    14VS 0 167.7 60. -90. -1. 1.

    BLANK BRANCH

    BLANK SWITCH

    BLANK SOURCE

    2POS 7.500000E+0001

    2XX0021-7.500000E+0001

    3POS XX0021 1.500000E+0002

    VSBLANK OUTPUT

    BLANK PLOT

    BEGIN NEW DATA CASE

    BLANK

    Figura. Circuito elctrico.

  • 7/28/2019 c9 Apendice Atp-emtp

    14/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez370

    A.7.2 Subrutinas de apoyo.

    De forma similar, la estructura tpica para las subrutinas de apoyo es la siguiente:

    BEGIN NEW DATA CASEPalabra especial que define la subrutina de apoyoDatos especficos del equipo$PUNCHBLANK CARDBEGIN NEW DATA CASEBLANK

    Ejemplo de palabras especiales son: XFORMER, BCTRAN, SATURA, HYSDAT, LINECONSTANTS, CABLE CONSTANTS, DATA BASE MODULE, etc.

    A continuacin se muestra un ejemplo de DATA MODULARIZATION. El elementocreado ser una resistencia.

    Figura. Esquema de una resistencia.

    Fichero 1R.DAT:

    BEGIN NEW DATA CASE --NOSORT--

    C Module for a resistor

    C USAGE: INCLUDE, 1R,FRNODE,TONODE,VALUER,H

    DATA BASE MODULE

    $ERASE

    ARG, FRNODE, TONODE,

    ARG, VALUER, H

    NUM, VALUER, H

    DUM,

    /BRANCH

    C

  • 7/28/2019 c9 Apendice Atp-emtp

    15/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez371

    BEGIN NEW DATA CASE

    C Just one comment card only, since the KASEND in startup file is 1.$PUNCH

    BEGIN NEW DATA CASE

    BLANK

    El comando $ERASE clarifica el buffer donde se guardar el fichero .pch obtenido. DUM denotael nombre de los nodos internos del componente (en este caso no hay).

    La compilacin del fichero anterior genera el siguiente fichero 1R.PCH:

    KARD 3 3 3 3

    KARG 1 2 3 4

    KBEG 3 9 27 80

    KEND 8 14 32 80

    KTEX 1 1 0 0

    $ERASE

    /BRANCH

    C

  • 7/28/2019 c9 Apendice Atp-emtp

    16/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez372

    A.7.3 Mdulos de simulacin integrada.

    A.7.3.1 TACS.

    El mdulo TACS (Transient Analysis of Control Systems) se puede usar para simular el controlde convertidores HVDC, sistemas de excitacin de mquinas sncronas, funciones de limitacin deintensidad de pararrayos, cebado de arcos en interruptores y, en general, aquellos dispositivos ofenmenos que no se pueden modelizar con los componentes elctricos existentes en ATP.

    La programacin con TACS se emplea para resolver las ecuaciones diferenciales yalgebraicas asociadas. Para ello, se utiliza en ATPDraw una representacin basada en diagramas debloques.

    Adicionalmente, cualquier seal obtenida con el mdulo TACS se puede utilizar dentro dela red elctrica simulada, incorporndola como una fuente de tensin o intensidad, o como unaseal que controle la interruptor, vlvula o diodo.

    A.7.3.2 MODELS.

    Es un lenguaje de programacin que se usa en ATP-EMTP para simular variablesdependientes del tiempo con caractersticas especiales. Los ficheros MODELS se crean de formaindependiente al programa principal y se pueden llamar desde cualquier programa cuantas veces sedesee. Su estructura es muy similar a la de otros lenguajes de programacin, como por ejemploFORTRAN, con sentencias del tipo FOR, IF, WHILE, etctera. Han de especificarse lassiguientes declaraciones:

    MODEL DEFAULT

    INPUT -- Name of input variables. Variable names separated by ',' or CR

    OUTPUT -- Name of output variables

    DATA -- Name of data variables

    VAR -- Name of local+output variables

    HISTORY -- Default values of variables and expressions {DFLT:n}

    INIT -- Initialization

    ENDINIT

    EXEC -- Execution

    ENDEXEC

    ENDMODEL

  • 7/28/2019 c9 Apendice Atp-emtp

    17/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez373

    Cada uso de un fichero MODEL se introduce con el comando USE, especificando las

    entradas y salidas correspondientes a esa llamada. Los nombres de los elementos definidos en unfichero MODEL no son visibles desde el exterior y se pueden elegir sin tener en cuenta lasposibles coincidencias con nombres usados fuera de este fichero MODEL.

    A.8 EDITOR GRFICO: ATPDraw.

    ATP-EMTP incluye una aplicacin grfica, denominada ATPDraw, capaz de crear y editarsistemas elctricos para poder simularlos de forma interactiva.

    La opcin Edit Commands (ver figura a continuacin) permite el uso y ejecucin de

    otros programas de ATP en el entorno grfico de ATPDraw.

    Figura. Edit Commands.

  • 7/28/2019 c9 Apendice Atp-emtp

    18/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez374

    Figura. Interaccin entre ATPDraw y otros programas de ATP.

    A.8.1 ENTORNO DE TRABAJO.

    La pantalla principal de este programa ofrece el aspecto habitual de las aplicacionesdesarrolladas en entorno Windows.

  • 7/28/2019 c9 Apendice Atp-emtp

    19/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez375

    Figura. Pantalla principal de ATPDraw con el menu de seleccin de componentes.

    A.8.3 CREAR UN CIRCUITO NUEVO.Seleccionando el comando New en el men de ficheros (Filemenu) o pinchando en el

    smbolo de pgina en blanco de la barra de herramientas, aparece una nueva ventana para laconstruccin de un circuito nuevo.

    Una vez construido el circuito, el paso final es dar nombre a los nodos del mismo.ATPDraw puede automticamente dar nombre a todos los nodos, sin embargo el usuario puedeasimismo hacerlo,ATP Make Names, o dar nombre solo a los nodos de especial inters. Estoltimo se realiza simplemente haciendo click con el botn derecho del ratn en el nodo anombrar.

  • 7/28/2019 c9 Apendice Atp-emtp

    20/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez376

    Antes de crear el fichero de entrada ATP o ejecutar la simulacin, se han de especificar los

    parmetros de la simulacin (parameters miscellaneus), que se encuentran en el men ATP Settings

    Figura. Opciones de simulacin.

    Donde:

    Delta T: Paso de tiempo en segundos. Tmax: Tiempo mximo de simulacin. Xopt=0: Incuctancias en mH. Copt=0: Condensadores en uF.

  • 7/28/2019 c9 Apendice Atp-emtp

    21/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez377

    Figura. Opciones de salida.

    En la figura anterior se muestra el contenido de la pestaa Integer, donde:

    Print freq: un valor de 500 significa que cada 500 pasos de simulacin, se escribiren el fichero .lis.

    Plot freq: un valor de 5 significa que cada 5 pasos de simulacin, se escribir en elfichero .pl4.

    Las variables declaradas en el lenguaje de programacin MODELS en la seccin VAR aparecen enla pestaa Record en la seccin Variable. Para grabar una variable (su evolucin podr observarseen el visualizador grfico PCPlot) ha de seleccionarse sta y presionar el botnAdd. Igualmente seprocede para su eliminacin, presionando en este caso el botn Remove.

  • 7/28/2019 c9 Apendice Atp-emtp

    22/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez378

    Figura. Pestaa Record.

    La opcin ATP Make File del men principal genera un fichero ATP del circuito. Podrguardarse en el directorio que deseemos. Finalmente se procede a su ejecucin mediante la opcinrun ATP del menATP.

    A.8.4 COMPONENTES DEFINIDOS POR EL USUARIO.

    Adems de todos los componentes estndar que contiene el editor grfico, el usuariopuede definir sus propios modelos mediante ficheros de texto. La estructura de estos ficherospuede ser mediante MODELS o DATA BASE MODULE. Estos modelos definidos por el

    usuario se guardan, por defecto, en los subdirectoriosModyUsp.

    El primer paso para la creacin de un nuevo modelo es la generacin del fichero de textodonde se describe el nuevo elemento que se quiere modelizar. Este fichero de texto se puededesarrollar mediante una de las siguientes opciones:

    Utilizando las plantillas que incluye el editor de textos PFE32.

    Escribiendo directamente desde el editor de texto, con la ayuda del ATP-EMTP RuleBook para definir correctamente la estructura del fichero.

    Utilizando la opcin Data Base Module (DBM), se quiere que tenga parmetrosvariables.

  • 7/28/2019 c9 Apendice Atp-emtp

    23/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez379

    El fichero de texto creado se tiene que compilar mediante alguna de las versiones de ATP,

    de modo que se obtenga otro fichero con la extensin .lib. Este ltimo fichero es el que sealmacena en el directorio Usp de ATPDraw.

    Por otro lado, tambin hay que definir un icono y una ventana que representen el nuevomodelo de componente creado mediante el fichero .lib. Para ello, con la opcin Objects UserSpecified New sup-file se crea un nuevo componente, cuya extensin es .sup.

    Figura. Opcin New sup-file.

    En la ventana que aparece a continuacin hay que indicar los datos que definen dichocomponente y el nmero de nudos disponibles para su conexin externa.

  • 7/28/2019 c9 Apendice Atp-emtp

    24/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez380

    Figura. Definicin de los datos y nudos del nuevo modelo.

    El nmero de nodos y datos debe coincidir con las declaraciones ARG y NUM del fichero

    escrito mediante Data Base Module. Permite la introduccin de hasta 36 datos, siendo el nmeromximo de nodos 12.

    Kind:

    Nodo MODELS:

    0: nodo de salida.1: nodo de entrada de intensidad.2: nodo de entrada de tensin.3: nodo de entrada de estado de un interruptor.

    4: nodo de entrada de una variable de una mquina.5: variable TACS.6: nodo de tensin de la parte imaginaria de estado-fijo (imssv).7: nodo de intensidad de la parte imaginaria de estado-fijo (imssi).8: salida de otro model.

    Nodo TACS:

    0: nodo de salida.1: nodo de entrada de suma positiva.2: nodo de entrada de suma negativa.

    3: nodo de entrada desconectado.

  • 7/28/2019 c9 Apendice Atp-emtp

    25/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez381

    En la ventana aparecen dos botones. Pulsando el botn donde aparece una lmpara

    elctrica, se abre un editor de iconos que permite disear el icono con el que representar el nuevocomponente creado. Mediante el botn en el que aparece un interrogante, se accede a un editor deayuda donde se puede escribir el texto que se desea mostrar, cuando se solicite ayuda desde la

    ventana de dilogo correspondiente a este componente.

    Tras definir todos los datos, se pulsa la opcin save as y se guarda el fichero con extensin.sup en el subdirectorio Usp, con lo que queda definido el nuevo componente. Para poder utilizareste nuevo modelo como un componente ms del circuito elctrico a disear, se hace uso de laopcin User specified Files y se selecciona el fichero .sup correspondiente al componente.

    Figura. Insercin de un elemento del usuario.

    Para la creacin de un modelo mediante el lenguaje de programacin MODELS se accedea la opcin New mod-file del men Objects (ver figura a continuacin), donde se muestra laestructura a seguir (presentada en un apartado anterior) para la implementacin del modelo. Elmanual MODELS IN ATP-Language ofrece al lector interesado detallada informacin acerca deeste lenguaje de programacin.

  • 7/28/2019 c9 Apendice Atp-emtp

    26/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez382

    Figura. Opcin para crear un modelo mediante MODELS.

    Ejemplo de fichero .mod:

    MODEL FLASH_1

    comment ****************************************************************

    * *

    * Function: set or cancel the gap firing control signal ** Inputs : voltage and current across ZnO resistor *

    * Output : the firing signal to the electrical ZnO component *

    * *

    ************************************************************* endcomment

    INPUT V1 -- Voltage on positive side of ZNO [V]

    V2 -- Voltage of negative side of ZNO [V]

    iczn -- ZNO current [Amps]

    DATA Pset -- power setting [Megajoules/msec]

    Eset -- energy setting [Megajoules]

    fdel -- firing delay [msec]

    fdur -- firing duration [msec]

    VAR power -- power into ZnO resistor [Watts]

    trip -- gap firing control signal [0 or 1]

    energy -- energy into ZnO resistor [Joules]

    tfire -- time at which the gap was last fired [sec]

    vcap -- voltage difference across series caps [Volts]

    OUTPUT trip

    HISTORY INTEGRAL(power) {DFLT:0}

    INIT trip:=0

    tfire:=0ENDINIT

  • 7/28/2019 c9 Apendice Atp-emtp

    27/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez383

    EXEC

    ------------------------------------------------------------------

    vcap:=V1-V2power:=vcap*iczn

    energy:=INTEGRAL(power)

    ------------------------------------------------------------------

    IF trip>0 -- is already firing

    AND t-tfire>fdur*1.e-3 -- has exceeded firing duration

    THEN

    trip:=0 -- cancel the firing signal

    tfire:=0 -- null the firing time

    ENDIF

    ------------------------------------------------------------------

    IF trip=0 -- is not signaling to fire

    AND tfire=0 -- firing condition not yet detected

    AND ( power >= Pset * 1.e9 -- power setting exceededOR energy >= Eset * 1.e6 ) -- energy setting exceeded

    THEN

    tfire:=t -- set the firing detection time

    ENDIF

    ------------------------------------------------------------------

    IF trip=0 -- is not signaling to fire

    AND tfire>0 -- firing condition has been detected

    AND t-tfire>=fdel*1.e-3 -- firing delay exceeded

    THEN

    trip:=1 -- set the firing signal

    ENDIF

    ------------------------------------------------------------------

    ENDEXEC

    ENDMODEL

    Igualmente ha de definirse un fichero .sup, asociado a dicho modelo, tal y como se mostranteriormente. El nombre de los nodos ha de coincidir con las variables especificadas en lassecciones INPUT y OUTPUT del fichero .mod. Ambos ficheros habrn de almacenarse en lacarpeta /MOD.

    A.8.5 COMPONENTES EN ATPDRAW.Para mostrar el men de seleccin de componentes basta con pinchar en un espacio vaco

    de la ventana del circuito con el botn derecho del ratn.

  • 7/28/2019 c9 Apendice Atp-emtp

    28/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez384

    Figura. Men de seleccin de componentes.

    A.8.5.1 PROBES & 3-PHASE.

    Figura. Men Probes & 3-phase.

  • 7/28/2019 c9 Apendice Atp-emtp

    29/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez385

    Probe Volt

    Mide la tensin a tierra de un nodo.

    Probe Branch volt.

    Mide la tensin entre dos nodos.

    Probe Curr

    Mide la intensidad entre dos nodos.

    Probe Tacs

    Mide una seal de tipo TACS.

    Splitter

    Permite la transformacin de un nodo monofsico a un nodo trifsico.

  • 7/28/2019 c9 Apendice Atp-emtp

    30/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez386

    A.8.5.2 BRANCH LINEAR.

    Figura. Men Branch Linear.

    Tabla. Branch linear.

  • 7/28/2019 c9 Apendice Atp-emtp

    31/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez387

    A.8.5.3 BRANCH NONLINEAR.

    Figura. Men Branch Nonlinear.

    Tabla. Branch Nonlinear.

  • 7/28/2019 c9 Apendice Atp-emtp

    32/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez388

    A.8.5.4 LINE LUMPED.

    Figura. Men Line Lumped.

    Tabla. Line Lumped.

  • 7/28/2019 c9 Apendice Atp-emtp

    33/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez389

    A.8.5.5 LINE DISTRIBUTED.

    Figura. Men Line Distributed.

    Tabla. Line Distributed.

  • 7/28/2019 c9 Apendice Atp-emtp

    34/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez390

    A.8.5.6 SWITCHES.

    Figura. Men Switches.

    Tabla. Switches.

  • 7/28/2019 c9 Apendice Atp-emtp

    35/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez391

    A.8.5.7 SOURCES.

    Figura. Men Sources.

    Tabla. Sources.

  • 7/28/2019 c9 Apendice Atp-emtp

    36/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez392

    A.8.5.8 MACHINES.

    Figura. Men machines.

    Tabla. Machines.

  • 7/28/2019 c9 Apendice Atp-emtp

    37/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez393

    A.8.5.9 TRANSFORMERS.

    Figure. Men Transformers.

    Tabla. Transformers.

  • 7/28/2019 c9 Apendice Atp-emtp

    38/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez394

    A.8.5.10 TACS.

    Coupling to circuit

    Figura. Acoplamiento de un nodo elctrico trifsico a TACS.

  • 7/28/2019 c9 Apendice Atp-emtp

    39/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez395

    TACS sources

    Tabla. TACS sources.

    TACS transfer function

    Tabla. TACS transfer function.

    Fortran statement

    Tabla. Fortran statement.

  • 7/28/2019 c9 Apendice Atp-emtp

    40/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez396

    TACS devices

    Tabla. TACS devices.

  • 7/28/2019 c9 Apendice Atp-emtp

    41/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez397

    A.8.5.11 Frequency comp.

    Figura. Men Frequency comp.

    Tabla. Frequency comp.

  • 7/28/2019 c9 Apendice Atp-emtp

    42/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez398

    A.9Visualizador grfico: PCPlot.

    El programa PCPlot permite la visualizacin interactiva de los resultados guardados en losficheros de extensin .pl4 que genera ATP.

    Al abrir un nuevo fichero aparece la ventana de dilogo de la figura, donde se eligen lasvariables a representar y el tipo de representacin deseado.

    Las variables escogidas se pueden dibujar en funcin del tiempo o en funcin de una deellas, segn se seleccione la opcin correspondiente. En el segundo caso, la primera de las variablesseleccionada constituye el eje x y aparece sealizada con una X a su izquierda.

    A.10Visualizador grfico: PlotXY.

    Permite representar hasta ocho curvas en la misma grfica, realizar la transformada deFourier de una seal, representar en la misma hoja curvas de tres ficheros diferentes, representarlas curvas en funcin del tiempo o hacer una representacin X-Y, hacer un escalado automtico deejes, acceder al valor instantneo de forma numrica, exportar datos, etc.

  • 7/28/2019 c9 Apendice Atp-emtp

    43/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez399

  • 7/28/2019 c9 Apendice Atp-emtp

    44/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez400

    A.11 Oscilaciones numricas.

    A.11.1 Causas de las oscilaciones numricas.

    El programa EMTP de transitorios electromagnticos usa el mtodo de integracinnumrica denominado regla trapezoidal. La regla trapezoidal es una tcnica de integracinnumrica de segundo orden que es simple de implementar, estable y rpido. Sin embargo, estambin susceptible de oscilaciones numricas cuando se diferencian tensiones o corrientesrespecto al tiempo de paso. Ello es comn cuando se considera el paso de corriente a travs deuna inductancia o la tensin entre los terminales de un condensador.

    Tensin en una inductancia:dt

    diLv =

    Corriente en un condensador:dt

    dvCi =

    Las oscilaciones numricas se hacen presentes cuando se da una variacin de corrienterespecto del tiempo demasiado elevada a travs de una inductancia o, igualmente, cuando apareceuna diferencia de tensin significativamente alta entre los terminales de un condensador. Una

    posible solucin consistira en la realizacin de un circuito ms exacto, incluyendo condensadoresparsitos y resistencias internas.

    A.11.2 Representacin matemtica.

    El problema puede ser representado considerando la ecuacin

    dt

    dyx=

    usando la regla trapezoidal. Podemos integrar ambos lados de la ecuacin,

    )()( ttyxdtty

    t

    tt

    +=

    A continuacin reescribimos la ecuacin usando la regla trapezoidal. Se ha de hacer notar queestamos buscando una aproximacin trapezoidal de la frontera formada por x, t y t-t.

    ))()((2

    )()( ttxtxt

    ttyty +

    +=

  • 7/28/2019 c9 Apendice Atp-emtp

    45/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez401

    Se resuelve por consiguiente esta ecuacin para x(t), que puede representar la intensidad a travsde un condensador o la diferencia de potencial entre los terminales de una inductancia.

    ))()((2

    )()( ttytyt

    ttxtx

    +=

    El resultado esperado habra de ser el mostrado en la siguiente figura:

    Figura. Representacin esperada de x(t) e y(t) para el caso de ejemplo.

    As, para el primer paso de tiempo se obtiene:

    ttx

    =

    2)(

  • 7/28/2019 c9 Apendice Atp-emtp

    46/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez402

    El siguiente paso, x(t-t), sera:

    tttttx

    =

    +

    =+

    2)0.10.1(

    22)(

    Un nuevo paso, sin embargo, resultara:

    tttx

    =+

    2)2(

    Se aprecia en este caso la existencia de una oscilacin, como aparece en la siguiente figura. El valormedio de x(t) corresponde a la forma de onda deseada.

    Figura. Actual x(t) con oscilaciones numricas.

  • 7/28/2019 c9 Apendice Atp-emtp

    47/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez403

    Una transicin de y(t) de dos pasos de tiempo en lugar de uno resolvera esta situacin, como se

    muestra a travs del siguiente conjunto de ecuaciones. Se asume en este caso y(t-t)=0, y(t)=0.5 ey(t+t)=1.0.

    ( )

    ( )

    ...

    0.0)3(

    0.00.10.12

    0.0)2(

    0.05.00.122

    5.0)(

    125.0)(

    =+

    =

    +=+

    =

    +

    =+

    =

    =

    ttx

    tttx

    ttttx

    tttx

    En la siguiente figura se ilustra los valores de x(t) e y(t) en consideracin de estas ecuaciones:

    Figura. Representacin de y(t) y x(t) con tiempo de subida ms lento.

  • 7/28/2019 c9 Apendice Atp-emtp

    48/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez404

    De todo ello, se deduce que reducir el paso de tiempo permite en algunos casos paliar o

    eliminar la aparicin de oscilaciones numricas. Se debe sealar sin embargo que la aparicin deoscilaciones debidas a operaciones de conmutacin no se vern reducidas de esta manera, ya que elinterruptor ideal siempre cambiar en un paso de tiempo.

    A.11.3 Ejemplos.

    A.11.3.1 Intensidad a travs de una inductancia.

    Figura. Circuito de ejemplo.

    Figura. Tensin en el nodo 1.

  • 7/28/2019 c9 Apendice Atp-emtp

    49/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez405

    Figura. Zoom de la tensin en el nodo 1.

    A.11.3.2 Tensin a travs de un condensador.

    Figura. Circuito de ejemplo.

  • 7/28/2019 c9 Apendice Atp-emtp

    50/55

  • 7/28/2019 c9 Apendice Atp-emtp

    51/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez407

    A.11.4 Opciones para eliminar las oscilaciones numricas.

    A.11.4.1Aadiendo elementos al circuito.

    El modelo del sistema elctrico puede carecer de ciertos detalles que den lugar aoscilaciones numricas. Por ejemplo, podra tenerse en cuenta la impedancia de arco de uninterruptor en su apertura. Sin embargo, modelar una impedancia de arco es claramente complejo.En este sentido, puede optarse por incluir un condensador parsito como se ilustra en la siguientefigura. En lugar de oscilaciones numricas tendrn lugar oscilaciones LC que podrn ocurrir en laprctica.

    Igualmente, la diferencia de potencial acaecida entre los terminales de un condensadorobvia la resistencia serie efectiva del condensador. Su consideracin introduce un trmino deamortiguamiento.

    Figura. Interruptor en serie sin y con un condensador parsito.

  • 7/28/2019 c9 Apendice Atp-emtp

    52/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez408

    Figura. Simulacin con un condensador parsito aadido.

    A.11.4.2 Rama snubber.

    Sin embargo, no siempre es prctico introducir estos elementos parsitos adicionales. Estees el caso de ciertos dispositivos electrnicos de potencia. En ellos, la rama paralela o snubber,proporciona un camino alternativo a la gran corriente inductiva generada por la carga, tenindose

    as un circuito RLC de amortiguamiento. La corriente decrece a cero cuando el condensador estcargado. De esta manera, adems, se eliminan o reducen los problemas relacionados conoscilaciones numricas.

    Interruptores electrnicos de potencias con grandes reas de seguridad (Safe Operating Area,SOA) tales como IGBTs y MOSFETs no necesitan snubbers. En estos casos el usuario delprograma puede aadir snubbers numricos. El valor del condensador puede estar comprendidoentre 1 y 2 nF, si no se considera la resistencia. Por el contrario, si se considera una rama RCsnubber, como se muestra en la figura siguiente, la constante de tiempo =RC habr de ser 2 3 veces el tiempo de paso elegido para la simulacin. Ntese, sin embargo, que los valores de R yC podrn variar segn el sistema elctrico a simular.

  • 7/28/2019 c9 Apendice Atp-emtp

    53/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez409

    Figura. Opciones snubber.

    Llegados a este punto se considera la simulacin del paso de corriente a travs de unainductancia presentado anteriormente teniendo en cuenta una rama snubber. En la siguientefigura se muestra los resultados de tensin e intensidad obtenidos, donde se aprecia la ausencia deoscilaciones.

    Figura. Ejemplo de la corriente a travs de la inductancia con snubber aadida.

  • 7/28/2019 c9 Apendice Atp-emtp

    54/55

    Energa Solar Fotovoltaica

    Trabajo Fin de Mster Juan Antonio Rodrguez Gonzlez410

    A.11.4.3 Disminuir el tiempo de paso.

    Un tiempo de paso pequeo puede eliminar o paliar en algunos casos las oscilacionesnumricas presentes en la simulacin. Ello depende de la topologa del circuito y de la presencia deresistencias en el mismo. Se ha de sealar que el tiempo de paso requerido puede estar lejos de unasimulacin prctica. Adems, debido a que la amplitud de la oscilacin vara con 1/t, un tiempode paso pequeo puede aumentar la misma.

    A.11.4.4Amortiguamiento.

    Alternativamente, puede considerarse la introduccin de resistencias de amortiguamiento

    en el circuito elctrico. En el caso de una inductancia, la resistencia habra de disponerse enparalelo con la misma, como se muestra en la siguiente figura, proporcionando un caminoalternativo a la intensidad cuando el interruptor se abre. La intensidad en el circuito RL decaer ascon el tiempo.

    Se ha de considerar con cuidado el valor de la resistencia. La constante de tiempo RL es:

    dR

    L=

    En una primera aproximacin, podemos elegir la constante de tiempo en relacin con eltiempo de paso t, como sigue:

    t

    LR

    d

    =

    donde es el factor de amortiguamiento. Un valor de 0.15 para puede suponer un buencompromiso entre amortiguamiento y exactitud.

    Figura. Inductancia con resistencia de amortiguamiento.

  • 7/28/2019 c9 Apendice Atp-emtp

    55/55

    Energa Solar Fotovoltaica

    Igualmente puede considerarse para el caso del condensador. En este caso la resistencia secoloca en serie con el condensador, como se muestra en la siguiente figura. Para la eleccin de laconstante de tiempo puede considerarse la siguiente expresin:

    C

    tKsRd

    2

    =

    Se recomienda Ks entre 0.15 y 13.34

    Figura. Condensador con resistencia de amortiguamiento.

    A.11.4.5 Interpolacin.

    El esquema de interpolacin detecta el cambio de paso de la intensidad y la tensin cuandose presentan oscilaciones. Entonces considera el tiempo anterior a la misma y realiza unainterpolacin lineal a la intensidad cero o en la operacin del interruptor. Las ecuacionesdiferenciales se resuelven en este punto y se comienza de nuevo en el siguiente tiempo de paso.

    Netomac fue probablemente el primero de los programas EMTP que incluyeroninterpolacin.

    Tambin se considera interpolacin para el uso de modelos de lneas de transmisin ypulsos de disparo en dispositivos electrnicos de potencia.

    A.11.4.6 Cambiando el mtodo de integracin.

    La versin de EMTP de DCG-EPRI (EMTP version 3.0, tambin conocida comoEMTP96), elimina el problema de las oscilaciones numricas ejecutando dos pasos de tiempomediante el mtodo de integracin Backward Euler.

    El usuario invoca el mtodo de integracin escribiendo CDA (Critical Damping Adjustment)a continuacin del BEGIN NEW DATA CASE