Modelagem de planos de recomposição de sistemas elétricos como grafos CPM/PERT

SISTEMAS DE POTÊNCIA

Modelagem de planos de recomposição de sistemas elétricos como grafos CPM/PERT

Alexandre de Assis Mota; Lia Toledo Moreira Mota; André Luiz Morelato França

DSEE/FEEC/UNICAMP - Avenida Albert Einstein, 400 - Cidade Universitária - Campinas (SP) - Brasil - C.P. 6101 - CEP 13083-852, Fone: +55-0xx-19-37883708, Fax : +55-0xx-19-32891395. [e-mail], [e-mail], [e-mail]

RESUMO

Esse trabalho formaliza um modelo de grafos CPM/PERT para viabilizar a representação de planos de recomposição de sistemas elétricos e expõe uma metodologia para derivação automática do mesmo a partir da lista de atividades de um procedimento de recomposição. Essa técnica de montagem automática foi aplicada com sucesso a uma instrução operativa de recomposição pós-blecautes do sistema interligado nacional (SIN), referente a uma porção da rede elétrica do sul do Brasil, que compreende as usinas hidrelétricas de Capivara e Taquaruçu, o corredor de transmissão associado a essas usinas e as subestações de Assis e Londrina.

Palavras-chave: Blecautes, Métodos de Caminho Crítico, Programação Orientada a Objetos, Recomposição de Sistemas de Energia Elétrica.

ABSTRACT

This work proposes a computational model based on PERT/CPM graphs to represent power system restoration plans. This model can be automatically derived using the precedence list represented by the restoration actions. This method was successfully applied to a restoration plan developed for a South-Brazilian power network, generating satisfactory results.

Keywords: Blackouts, Critical Path Methods, Objectoriented programming, Power System Restoration.

1 INTRODUÇÃO

No contexto do planejamento da recomposição de sistemas de potência, torna-se bastante desejável a especificação de um modelo computacional para representar planos de recomposição que viabilize, na medida do possível, o tratamento dos problemas inerentes às atividades de reenergização, segundo uma abordagem lógico-matemática bem definida. Deve ser objetivada, na escolha dessa abordagem, a sua capacidade de representar a natureza bastante diversa das inúmeras tarefas que são executadas durante o processo de recomposição, envolvendo atividades de caráter técnico (como, por exemplo, a atuação em controles e a atenção aos limites dos equipamentos e linhas de transmissão), gerencial (como a coordenação de interligação de diferentes áreas geo-elétricas), logístico (como o deslocamento de viaturas e equipes em campo) e estratégico (como a representação de cargas prioritárias para reenergização) (Adibi, 2000). Mais ainda, dois aspectos muito importantes têm que ser levados em conta na definição de um modelo para os planos de recomposição: a possibilidade de qualificar as diferentes linhas de ação definidas para abordar a reenergização sob um mesmo cenário de falta de energia, contemplando os diferentes aspectos do processo, e a viabilidade para modelar o caráter impreciso das atividades de recomposição, permitindo uma abordagem estocástica na sua especificação.

Assim, esse trabalho propõe uma metodologia para construção de um modelo computacional para os planos de recomposição (Mota et. al., 2004), baseado em métodos de caminho crítico (usualmente definidos como CPM/PERT), para representar o seqüenciamento temporal e a dependência lógica entre as atividades de recomposição.

Os métodos de caminho crítico CPM ("Critical Path Method") e PERT ("Program Evaluation Review Technique") são métodos de planejamento e programação de processos desenvolvidos, durante a década de 1950, pela Marinha dos Estados Unidos e pela companhia DuPont, com o objetivo de otimizar os processos com múltiplas atividades encadeadas, tanto no aspecto temporal quanto de aplicação de recursos (Moder, 1970; Sivazlian, 1975). Esses métodos definem uma base formal e uma abordagem generalizada para o gerenciamento de projetos e abrangem as características desejadas de flexibilidade e tratamento de imprecisões, desejadas para o modelo de planos de recomposição. Ambas as técnicas presumem que qualquer projeto ou plano pode ser desmembrado em uma seqüência lógica de ações predefinidas (ou atividades) encadeadas. A realização dessas atividades conduz à conclusão do projeto. Outra hipótese importante adotada nesses métodos é que as atividades possuem uma execução lógica de natureza tecnológica, ou seja, apesar de poderem ser executadas de uma maneira relativamente independente, elas possuem, necessariamente, precedências tecnológicas entre si. Na sua formulação, assume-se que a relação entre duração e custos de cada tarefa do processo é bem conhecida, o que permite a estimação de diversos parâmetros de gerenciamento e desempenho, como a alocação ótima de recursos e o controle de custos do projeto.

2 PLANOS DE RECOMPOSIÇÃO

No Brasil, o produto gerado ao final da etapa de planejamento da recomposição é um conjunto de textos, denominados "instruções operativas de recomposição", que faz parte dos procedimentos de rede do Sistema Interligado Nacional (SIN), disponibilizados pelo Operador Nacional do Sistema Elétrico (ONS).

Assim, uma instrução operativa de recomposição, ou simplesmente instrução operativa (IO), é um documento que descreve as ações que devem ser tomadas e as condições que devem ser satisfeitas na operação do SIN, após um evento de interrupção de grandes proporções no fornecimento de energia elétrica.

2.1 Formato Tabular de Instruções Operativas

É sugerida, pelo ONS, a representação das instruções operativas em um formato tabular, para facilitar o entendimento, pelos operadores e engenheiros, das atividades descritas no plano sob condições de recomposição. Esse formato torna explícitas características importantes do processo, tais como a ordem de uma determinada atividade na seqüência planejada de ações, o agente executor dessa atividade, os procedimentos específicos e as condições (ou controles) que devem ser atendidas antes que uma ação em questão possa se iniciar. Essas características ficam organizadas de uma maneira lógica, definindo o conjunto de procedimentos e ações em uma cadeia ordenada de atividades.

Como exemplo, a tabela I apresenta o conjunto de ações para a reenergização da área Capivara (área 16 da região Sudeste do Brasil), representada geograficamente na figura 1, que é formada pelas Usinas Capivara e Taquaruçu, pelas subestações de Assis e Londrina e pelo tronco de transmissão associado, extraído da instrução operativa correspondente (ONS, 2004).

2.2 Lista de Ações de Recomposição

No Brasil, como as ações estão delineadas em uma seqüência lógica nesse formato tabular, é possível determinar listas de atividades para cada instrução operativa, conforme detalhado, anteriormente, na análise geral de planos de recomposição. Essas listas devem explicitar a relação de precedência entre as diferentes ações de recomposição. A tabela II apresenta uma possível lista de atividades correspondente às ações seqüenciais de reenergização dos diferentes equipamentos elétricos da área Capivara, extraídas por inspeção da tabela I. Cada linha da tabela II identifica uma atividade com uma letra (A, B, C, etc.), apresenta brevemente a descrição da mesma e especifica suas atividades precedentes, ou seja, as ações que devem ser completadas para que a atividade de recomposição em questão possa se iniciar.

3 GRAFOS CPM/PERT

Os métodos CPM e PERT são baseados na identificação do caminho crítico de um processo, isto é, da seqüência das atividades que determinam a duração total do projeto. Para atingir esse objetivo, utilizam uma representação gráfica das relações entre as atividades do processo, na forma de um grafo orientado (denominado grafo CPM/PERT).

Assim, a abordagem de um problema (como a recomposição de sistemas elétricos) através de métodos de caminho crítico é realizada a partir da definição inicial da lista de atividades do processo sob análise. Definida a lista de atividades, o grafo CPM/PERT pode ser montado, para viabilizar a aplicação dos algoritmos de cálculo de custos e tempos envolvidos com o projeto de maneira eficiente. Nesse trabalho, são abordados os grafos do tipo "AoA" ("Activities-on-Arcs") ou redes orientadas por eventos, onde as atividades do processo são associadas às arestas do grafo CPM/PERT em construção. Os vértices do grafo são denominados etapas, que recebem usualmente uma identificação seqüencial. Essas etapas são os estágios do processo onde ocorre a convergência de atividades e, portanto, onde fica explícita a interdependência entre as mesmas. Nesse tipo de rede, pode existir apenas uma atividade entre cada par de etapas. Assim que o grafo CPM/PERT esteja disponível pode-se aplicar as técnicas de identificação de caminho crítico.

4 MODELOS ADOTADOS PARA AS ETAPAS E ATIVIDADES DA RECOMPOSIÇÃO

Foi adotada, aqui, uma abordagem orientada a objetos na especificação dos elementos do grafo CPM/PERT. A utilização de programação orientada a objetos para a definição de modelos de atividades e de etapas do grafo permite a adoção de uma representação física para os mesmos, mais próxima da sua natureza tecnológica/procedural. A despeito da maior utilização de memória computacional, essa abordagem permite que os elementos do grafo sejam instanciados de maneira independente, preservando características individuais das atividades. Além disso, torna fácil a extração de diferentes formas de representação matemática do grafo, como listas de arestas, listas de incidência ou matrizes de incidência.

4.1 Modelo de Atividades

As atividades são representadas, na rede "AoA", por arcos orientados do grafo CPM/PERT que se iniciam e terminam em etapas distintas. Esses arcos ficam relacionados a características como as diferentes estimativas de tempo de conclusão de cada atividade, recursos envolvidos e relações entre custo e tempo do processo. Utilizando uma abordagem baseada em objetos, o instanciamento das atividades de um plano pode ser diretamente realizado a partir da lista de atividades associada ao processo tecnológico em questão. A figura 2 apresenta o modelo genérico para atividades de um processo.

4.2 Modelo de Etapas

Os nós do grafo CPM/PERT representam as etapas do processo sob análise, e correspondem aos pontos temporais onde algumas atividades convergem para dar início a outras atividades. Neles, fica explícita, portanto, a dependência entre os diferentes componentes do processo. Elas também podem ser caracterizadas diretamente a partir da lista de atividades, de acordo com o modelo explicitado na figura 3.

5 METODOLOGIA DE MONTAGEM AUTOMÁTICA DOS GRAFOS CPM/PERT

Uma vez definidos os modelos de objetos para as etapas e atividades, é necessário estabelecer um procedimento para satisfazer as relações de conectividade entre os elementos, caracterizando conseqüentemente a estrutura topológica do grafo CPM/PERT. Os passos que se seguem definem uma metodologia computacional, proposta por Mota et. al. (2004), que foi utilizada neste trabalho para realizar essa tarefa.

Passo 1: Definição das arestas

Neste passo, são definidas as arestas orientadas dos grafos CPM/PERT utilizando como entrada a estrutura da lista de atividades. Para cada componente dessa lista, são atribuídos os parâmetros de identificação e descrição de atividades precedentes. Este passo está representado pelo algoritmo a seguir.

Seja A = conjunto de atividades

Para cada linha da lista de atividades i

a = nova atividade

Atribuir nome de a conforme a linha da lista de atividades

Atribuir lista de pré-requisitos de a conforme a linha da lista de atividades

Adicionar a ao conjunto de atividades A

Passo 2: Definição das etapas de início das atividades

Uma vez definidos os objetos associados às atividades, pode-se instanciar novos objetos, representando as etapas de início dos mesmos. Esses objetos do tipo nó do grafo são derivados pela inspeção das configurações de pré-requisitos definidas em cada objeto atividade, sendo-lhes atribuído um nome arbitrário e seqüencialmente numerado, conforme o algoritmo a seguir. Além disso, um teste para verificar a existência da etapa deve ser conduzido para cada atividade, uma vez que atividades diferentes podem se iniciar no mesmo nó.

Seja N = conjunto de nós k = 0; Para cada atividade i definida em A P = conjunto de atividades pré-requisito de i Se não existe nó contendo exatamente todos os destinos de P k=k+1 n = novo nó que contém exatamente todos os destinos de P nome de n = "Nó" + k Adicionar n ao conjunto de nós N k=k+1 n = novo nó que não contém nenhuma atividade nome de n = "Start" Adicionar n ao conjunto de nós N

Passo 3: Conexão das origens das atividades

Uma vez que o conjunto inicial de nós é derivado a partir das listas de requisitos das próprias atividades, a conexão das origens das atividades às etapas já definidas é trivial. O algoritmo a seguir sintetiza as ações a serem tomadas nesse passo. É importante salientar que, computacionalmente, esse passo e o anterior podem ser realizados praticamente de maneira simultânea, embora, neste trabalho eles tenham sido explicitados separadamente com o propósito de proporcionar um melhor entendimento do processo.

Para cada atividade i definida em A P = conjunto de atividades pré-requisito de i Para cada nó j definido em N Se j contém exatamente todos os destinos de P então Faça a origem de i igual a j

Passo 4: Conexão dos possíveis destinos

Nesse ponto do processo, as atividades já possuem suas origens conectadas a nós definidos, mas seus destinos ainda estão desconexos. Esse passo consiste em conectar os destinos, para cada atividade definida em A, aos nós pertencentes ao conjunto N já definido, quando apenas uma etapa depende da conclusão da atividade em questão. Entretanto, é possível que mais de uma etapa dependa da conclusão de uma única atividade. Essas atividades são marcadas como possuindo dependência múltipla e tratadas nos passos que se seguem. O algoritmo a seguir sintetiza esse passo.

n = novo nó que não contém nenhuma atividade nome de n = "End" Para cada atividade i definida em A ND = número de etapas dependentes de i = 0 Para cada nó j definido em N Se j contém o destino de i Marcar j como dependente de i ND = ND + 1 Se ND > 1 Marcar i como tendo dependência múltipla Senão Se ND = 1 Conectar destino de i ao nó dependente marcado j Se ND = 0 Conectar destino de i ao nó final n Adicionar n ao conjunto de nós N

Passo 5: Criação de etapas de destino adicionais

Nesse passo, representado no algoritmo que se segue, são atribuídas etapas de destino adicionais para as atividades com dependência múltipla. Essas etapas correspondem ao momento da conclusão da atividade em questão e são representadas por nós conectados ao destino das atividades que não puderam ser conectados no passo anterior. A criação de tais etapas é necessária para que o passo seguinte, onde a conectividade do grafo é concluída, possa ser executado.

Seja NA = conjunto de nós adicionais Para cada atividade definida i Se i possui dependência múltipla n = novo nó que contém i nome de n = "Nó adicional" + nome de i Adicionar n ao conjunto de nós NA

Passo 6: Criação de atividades-fantasma

Para a execução desse passo, é necessário instanciar objetos atividade com duração e recursos nulos, ou seja, as atividades-fantasma (ou "dummies"). Nesse passo, as etapas que ainda não tem todos seus pré-requisitos atendidos sofrem uma busca conduzida pelas etapas adicionais criadas no passo anterior; quando a etapa de destino da atividade (ou atividades) pré-requisito é identificada, uma atividade-fantasma é criada, iniciando-se no nó de destino adicional e terminando na etapa que iniciou a busca. O objetivo de uma atividade fantasma é, portanto, atender às restrições de precedência impostas pela lista de atividades e que não podem ser satisfeitas usando apenas o conjunto de atividades reais. O algoritmo a seguir representa esse passo.

Seja AF = conjunto de atividades-fantasma Para cada nó adicional na definido em NA i = atividade incidente em na Para cada nó j definido em N Se j contém i então af = nova atividade fantasma com duração nula Faça a origem de af igual à na Faça o destino de af igual a j Adicionar af ao conjunto de atividades-fantasma AF

Passo 7: Filtragem de atividades-fantasma desnecessárias

Esse passo corresponde ao final do processo de atribuição das conexões do grafo CPM/PERT, onde as atividades-fantasma que não são relevantes são filtradas e descartadas da estrutura. A presença de atividades-fantasma desnecessárias pode ocorrer quando, para uma atividade com dependência múltipla, existe uma etapa que depende apenas da conclusão da atividade em questão. Nesse caso, essa etapa corresponde à conclusão individual da atividade e, portanto, coincide com a etapa adicional que foi criada para destino da atividade no passo 5. Embora esse passo tenha sido explicitado neste trabalho separadamente, para uma compreensão melhor do processo de construção do grafo, ele pode ser incorporado iterativamente aos passos 5 e 6, visando a obtenção de uma maior eficiência computacional, assim como no caso do passo 3. O algoritmo que se segue descreve esse passo.

Para cada atividade-fantasma af definida em AF j = nó de destino de af Se j contém exatamente af então n = nó de origem de af transfira todas as atividades dependentes de j para n elimine j atualize N elimine af atualize AF

6 EXEMPLOS DE APLICAÇÃO

6.1 Aplicação em um Processo Genérico

Pelo fato de agregar uma maior flexibilidade na sua formulação, decorrente da abordagem orientada a objetos, a metodologia proposta pode ser utilizada para modelar diferentes processos tecnológicos. Para verificar essa característica, a técnica de montagem automática de grafos CPM/PERT foi aplicada a um exemplo extraído de Goldbarg (2000), correspondente ao processo de montagem de uma cabine de simulação, cuja lista de atividades é reproduzida na figura 4.

A aplicação dos passo 1-7 do método de montagem automática gerou os resultados descritos a seguir.

Passo 1: Instanciamento das arestas

As arestas são instanciadas diretamente por inspeção da lista de atividades, resultando em nove objetos, conforme a figura 5.

Passo 2: Instanciamento das etapas de início

Assim como no caso das arestas, as etapas também são instanciadas diretamente por inspeção da lista de atividades, resultando em sete objetos (seis etapas extraídas da lista e mais uma instanciada como nó de início do processo), conforme ilustrado na figura 6.

Passo 3: Conexão das origens

De posse da lista de arestas e da lista de vértices, os inícios das atividades são conectados às etapas correspondentes, conforme ilustrado na figura 7.

Passo 4: Conexão dos possíveis destinos

A primeira ação nesse passo é o instanciamento de uma atividade adicional para representar o nó final do processo (nesse caso, instanciado como o vértice número 7). A seguir, os destinos das atividades são conectados ao conjunto de nós disponível (quando possível), resultando no grafo ilustrado na figura 8.

Os passos 5 (instanciamento de destinos adicionais), 6 (instanciamento de atividades-fantasma) e 7 (filtragem de atividades-fantasma) não são necessários nesse exemplo, uma vez que todas as atividades já possuem suas etapas inicial e final definidas.

6.2 Aplicação em um Plano de Recomposição

Para ilustrar a aplicação dessa metodologia de montagem automática de redes à modelagem de planos de recomposição do SIN, a seqüência de sete passos do algoritmo de montagem automática foi aplicada à lista de atividades associada à instrução operativa da região Capivara, explicitada previamente na tabela II.

Observando essa lista de atividades, é possível explorar as relações de precedência entre as diferentes ações que compõem o plano para determinar, então, a estrutura topológica do grafo de caminho crítico correspondente à IO de recomposição da área Capivara. A figura 9 representa o grafo CPM/PERT obtido para esse caso.

6.3 Observações Adicionais

Com o grafo CPM/PERT de uma IO disponível, é possível conduzir uma série de análises que fornecem informações relevantes acerca do processo de recomposição em questão. Os procedimentos envolvidos em tais análises estão descritos em maior detalhe em Mota et. al. (2005) e em Moder (1970). Como exemplo, a tabela III apresenta os tempos e respectivos valores de desvio padrão, estimados para a conclusão das atividades de recomposição da área Capivara, representadas anteriormente na tabela II.

Os resultados do processamento dessa estrutura com os métodos convencionais de caminho crítico e com as técnicas citadas, forneceram um tempo total estimado para recomposição dessa área igual a 93,5 unidades de tempo, com 9 atividades críticas (A, B, C, D, E, F, G, O e P) e um custo estimado para interrupção de aproximadamente 650.000 dólares, considerando-se valores de custo unitário de interrupção por segmento de consumo descritos por Gomes et. al. (1997).

A execução dessas tarefas tomou menos que um segundo, em um microcomputador pessoal, para processar desde a aquisição das informações acerca da lista de atividades (disponibilizada em arquivos tipo texto), passando pela geração do grafo de caminho crítico e terminando com a execução das funções de análise de tempos e custos associados ao grafo CPM/PERT.

É importante salientar que essa estrutura conexa de interdependência entre as ações de recomposição (o grafo CPM/PERT) pode ser aproveitada para a construção de representações gráficas em interfaces de usuário, permitindo que os operadores acompanhem, visualmente, os procedimentos de reenergização, durante treinamentos e no uso de programas de apoio à recomposição.

Finalmente, pode-se destacar que se o sistema sob análise fosse de maior porte (maior número de barras e linhas) ou se a solução do problema estivesse mais distante eletricamente do local sob recomposição no início do grafo, o modelo resultante da aplicação da metodologia teria um maior número de nós e arcos, uma vez que o plano de recomposição correspondente teria um maior número de ações de recomposição; os procedimentos para modelagem, entretanto, podem ser aplicados exatamente da mesma maneira, fornecendo resultados com esforço computacional proporcionalmente maior (Mota & Mota, 2006).

7 CONCLUSÕES

Nesse trabalho, é explorado o modelo baseado em grafos dirigidos, do tipo CPM/PERT, para representar as seqüências de ações de reenergização de sistemas elétricos de potência após blecautes. Entretanto, é fundamental destacar que a operação de redes de energia no seu estado restaurativo envolve dois aspectos que tornam essa modelagem uma tarefa árdua, quando se considera a construção manual dos diagramas de caminho crítico:

(a) o planejamento da recomposição necessita da análise de um grande número de opções de reenergização para um mesmo cenário de interrupção e, assim, requer que um grande número de diferentes listas de atividades seja analisado;

(b) na utilização das instruções operativas como apoio à recomposição em tempo real, pode ser necessária a alteração de algumas das atividades que compõem o plano em decorrência de imprevistos na operação e, portanto, o diagrama do processo deve ser atualizado, o que, manualmente, é extremamente difícil nessa situação.

Nesse sentido, a metodologia de modelagem automática dos planos, conforme exposta previamente, permite o armazenamento das informações contidas nos planos de recomposição em estruturas matemáticas e computacionais bem definidas. Essas estruturas podem ser exploradas para fornecer um grande número de informações sobre os planos de recomposição, relacionadas a sua duração, custos e decisões sobre investimentos (Moder, 1970; Mota, 2005; Sivazlian, 1975) em um tempo bastante reduzido, o que vem ao encontro das necessidades específicas do problema de restabelecimento do serviço após distúrbios de grandes proporções.

AGRADECIMENTOS

Os autores gostariam de agradecer à Fundação de Amparo a Pesquisa do Estado de São Paulo (FAPESP) pelo suporte durante o desenvolvimento deste trabalho.

Artigo submetido em 01/11/2005

1a. Revisão em 12/06/2006

2a. Revisão em 09/10/06

Aceito sob recomendação do Editor Associado Prof. Carlos A. Castro

Referências bibliográficas

• Adibi, M.M. (2000). Power System Restoration - Methodologies and Implementation Strategies, New York: IEEE Press.
• Adibi, M.M., & D.P. Milanicz (1999). Estimating Restoration Duration, IEEE Transactions on Power Systems, vol. 14, n. 4, Nov.
• Goldbarg, M.C., & H.P.L. Luna (2000). Otimizaçăo Combinatória e Programaçăo Linear - Modelos e Algoritmos, Rio de Janeiro: Editora Campus.
• Gomes, P. & M. Th. Schilling (1997). Custo de Interrupçăo: Conceituaçăo, Metodologia de Avaliaçăo, Valores Existentes e Aplicaçőes, Anais do XIV SNPTEE - Seminário Nacional de Produçăo e Transmissăo de Energia Elétrica, Belém (PA).
• Moder, J.J. & C.R. Phillips (1970). Project Management with CPM and PERT, New York: Van Nostrand Reinhold Company.
• Mota, A.A., L.T.M. Mota e A. Morelato (2004). Metodologia Orientada a Objetos para a Construçăo Automática de Grafos CPM/PERT, Anais do XXXVI SBPO - Simpósio Brasileiro de Pesquisa Operacional, Săo Joăo Del Rei (MG).
• Mota, A.A., L.T.M. Mota e A. Morelato (2005). Avaliaçăo da Qualidade de Planos de Recomposiçăo de Sistemas de Energia Elétrica, Anais do XVIII SNPTEE - Seminário Nacional de Produçăo e Transmissăo de Energia Elétrica, Curitiba (PR).
• Mota, A.A. & L.T.M. Mota (2006). Análise dos Tempos de Recomposiçăo das Áreas Capivara e Marimbondo Relatório 02/2006, Operador Nacional do Sistema Elétrico - ONS, Brasília (DF).
• Morelato, A. & A. Monticelli (1989). Heuristic Search Approach to Distribution System Restoration", IEEE Transactions on Power Delivery, vol. 4, n.4.
• ONS - Operador Nacional do Sistema Elétrico (2004). Procedimentos de Rede: Módulo 10, Centro Nacional de Operaçăo do Sistema - CNOS, Brasília (DF).
• Sivazlian, B.D. & L.E. Stanfel (1975). Optimization Techniques in Operations Research, New Jersey: Prentice-Hall.

Datas de Publicação

• Publicação nesta coleção
  02 Abr 2007
• Data do Fascículo
  Dez 2006

Histórico

• Aceito
  09 Out 2006
• Revisado
  12 Jun 2006
• Recebido
  01 Nov 2005