O modelo de Gurson para dano dúctil : estratégia computacional e aplicações

Este trabalho trata do emprego do modelo de dano dúctil de Gurson e de alguns aspectos relativos a sua implementação computacional. Emprega-se como ponto de partida o programa Metafor, código de elementos finitos voltado para a simulação de problemas de conformação mecânica, contemplando portanto grandes deformações plásticas. Estudam-se alguns casos simples procurando identificar a influência de cada parâmetro do modelo de Gurson na resposta. É discutida a aplicação dos modelos de nucleação de vazios usualmente empregados em conjunto com a superfície de escoamento de Gurson em situações onde há reversão de solicitação. São propostas alterações nos modelos usuais de nucleação de forma a minorar resultados incoerentes verificados. Apresenta-se um algoritmo alternativo ao método de Newton-Raphson para a solução do sistema de equações associado à forma fraca das equações de equilíbrio, em situações onde a rigidez fica muito pequena. Tal algoritmo é implementado e testado. Um algoritmo de integração de tensões é apresentado e implementado, mostrando-se sua vantagem do ponto de vista de robustez, ou seja, obtenção da resposta correta a despeito do tamanho de incremento empregado. Também é discutido um procedimento para desconsiderar elementos que apresentem ruptura local Estudam-se problemas envolvendo geração de calor por deformação plástica em conjunto com a formulação do modelo de Gurson, mostrando que a consideração de acoplamento termo-mecânico com dano traz vantagens em algumas simulações. É discutida a implementação do modelo de Gurson sobre a versão tridimensional do Metafor, que contempla a formulação lagrangiana-euleriana arbitrária (LEA). Ao final, é feita uma simulação inicial do comportamento de espumas metálicas.

Implementação do conceito de marketing em jornais impressos : possibilidades de conciliação entre os parâmetros de desempenho da administração e os valores tradicionais do jornalismo

O principal objetivo desta tese é investigar a implementação do conceito de marketing em empresas jornalísticas e seus efeitos sobre a qualidade do conteúdo editorial e o desempenho financeiro dessas orga nizações. Há controvérsias sobre esse tema, dividindo os pesquisadores da área em duas correntes opostas: a pessimista e a otimista. Os pessimistas consideram incompatíveis os parâmetros de desempenho das ciências gerenciais, como sucesso financeiro e mercadológico, com a função maior dos meios de comunicação de massa — difundir as informações necessárias para a manutenção da sociedade democrática. Já os otimistas sustentam que são exatamente os veículos que adotam o conceito de marketing aqueles que enfatizam mais a qualidade do conteúdo editorial, sendo possível, portanto, conciliar os parâmetros de desempenho financeiro e mercadológico com os valores tradicionais do jornalismo. Foi desenvolvido um modelo teórico, baseado na perspectiva otimista, que foi testado através de um survey cross-section com editores-chefe e diretores comerciais de 184 diários brasileiros. Através da modelagem de equações estruturais, confirmaram-se sete das nove hipóteses, endossando a maioria das premissas otimistas. Os resultados sustentam os efeitos positivos da orientação para o mercado leitor sobre a integração entre os departamentos editorial e comercial e sobre a ênfase na qualidade do conteúdo editorial. Confirmou-se também o efeito positivo da orientação para o mercado anunciante sobre a integração entre os departamentos. A ênfase na qualidade do conteúdo editorial afeta positivamente o desempenho das vendas de espaço publicitário e o desempenho financeiro. E o desempenho das vendas de espaço publicitário afeta positivamente o desempenho financeiro. Não foi encontrado efeito significativo da orientação para o mercado anunciante sobre as vendas de espaço publicitário e foi rejeitada a hipótese de que a integração entre os departamentos editorial e comercial influencia as vendas de espaço publicitário. Em resumo, os resultados do estudo sustentam a possibilidade de conciliação entre os parâmetros de desempenho da administração e os valores tradicionais do jornalismo.

Sincronização idêntica em redes de osciladores caóticos

A presente dissertação versa sobre a sincronização idêntica em redes de osciladores caóticos. Uma perspectiva razoavelmente histórica sobre a literatura da área é apresentada . O conceito de caos é introduzido junto com outras idéias da dinâmica não-linear: sistemas dinâmicos, exemplos de sistemas, atratores, expoentes de Liapunov, etc. A integração numérica de equações diferenciais é largamente utilizada, principalmente, para o cálculo de expoentes e o desenho do diagrama de fases. A sincronização idêntica é definida, inicialmente, em redes que não passam de um par de osciladores. A variedade de sincronização (conjunto de pontos no espaço de fases no qual a solução do sistema é encontrada se há sincronização) é determinada. Diferentes variantes de acoplamentos lineares são enfocadas: acoplamento interno, externo, do tipo mestre-escravo e birecional, entre outras. Para detectar sincronização, usa-se o conceito de expoente de Liapunov transversal, uma extensão do conceito clássico de expoente de Liapunov que caracteriza a sincronização como a existência de um atrator na variedade de sincronização. A exposição é completada com exemplos e atinge relativo detalhe sobre o assunto, sem deixar de ser sintética com relação à ampla literatura existente. Um caso de sincronização em antifase que usa a mesma análise é incluído. A sincronização idêntica também é estudada em redes de osciladores idênticos com mais de dois osciladores. As possibilidades de sincronização completa e parcial são explanadas. As técnicas usadas para um par de osciladores são expandidas para cobrir este novo tipo de redes. A existência de variedades de sincronização invariantes é considerada como fator determinante para a sincronização. A sincronização parcial gera estruturas espaciais, analisadas sob a denominação de padrões. Algumas relações importantes entre as sincronizações são explicitadas, principalmente as degenerescências e a relação entre a sincronização parcial e a sincronização completa do respectivo estado sincronizado para alguns tipos de acoplamento. Ainda são objetos de interesse as redes formadas por grupos de osciladores idênticos que são diferentes dos osciladores dos outros grupos. A sincronização parcial na qual todos os grupos de osciladores têm seus elementos sincronizados é chamada de sincronização primária. A sincronização secundária é qualquer outro tipo de sincronização parcial. Ambas são exemplificadas e analisadas por meio dos expoentes transversais e novamente por meio da existência de invariantes de sincronização. Obtém-se, então, uma caracterização suficientemente ampla, completada por casos específicos.

Efeitos da aplicação de altas pressões sobre compostos com estrutura de escuterudita

Neste trabalho foi estudado o comportamento em altas pressões e temperatura ambiente de dois compostos com estrutura de escuterudita, CoSb3 e LaFe3CoSb12. A determinação dos parâmetros das equações de estado isotérmicas, bem como a verificação das diferenças relativas nas curvas de compressibilidade desses compostos, foram os objetivos principais deste trabalho. Ambos compostos foram sintetizados conforme rota proposta na literatura. O sucesso da síntese foi verificado por meio de difração de raios X, não sendo observada a presença de fases contaminantes. Como meio gerador de altas pressões foi utilizada uma câmara de bigornas de diamante (DAC). A evolução das estruturas destes compostos com a pressão foi acompanhada por meio de difração de raios X por dispersão em energia. Até o limite de pressão estática de cerca de 14 GPa (definido pelas condições experimentais previamente existentes no laboratório), não foram observadas anomalias nas curvas de volume versus pressão desses compostos ou mesmo evidências de transições de fase. Com o objetivo de estender este limite de pressão, foram feitas melhorias na célula de alta pressão existente no laboratório, possibilitando assim a geração de pressões de até 40 GPa. Três equações de estado isotérmicas foram ajustadas aos dados experimentais, utilizando como parâmetros no ajuste o módulo volumétrico isotrópico e sua derivada primeira com relação à pressão. Como resultado deste procedimento, verificou-se que a escuterudita com lantânio apresenta um módulo volumétrico ligeiramente menor que o CoSb3, sendo, portanto, mais compressível, apesar da presença do íon lantânio preenchendo as cavidades deste composto. Foi observada também uma mudança na compressibilidade da escuterudita CoSb3 acima de 20 GPa, devido, talvez, ao congelamento do meio transmissor de pressão ou ao pinçamento do rubi (calibrante de pressão) entre as bigornas da DAC. O CoSb3 foi comprimido até a pressão de 40 GPa, apresentando indícios de um estranho comportamento na seqüência decrescente de pressão. Após a descompressão, o parâmetro de rede desse composto, a pressão ambiente, mostrou-se maior que o da fase original. Este aumento do parâmetro de rede não foi acompanhado de alterações significativas do espectro de difração de raios X. Este efeito anômalo dever ser sujeito a verificação por meio de novos experimentos.

Geração automática de leiaute através de matriz de células- MARTELO

Esta dissertação enquadra-se em um processo de busca de soluções para a geração do leiaute de circuitos integrados que permitam aumentar a qualidade da previsibilidade do comportamento de um circuito após a sua implementação. Isso é importante face ao crescimento dos problemas referentes aos efeitos elétricos adversos que surgem em nanocircuitos, tais como eletromigração, efeito antena, contatos mal formados e outros, assim como o aumento da variabilidade do processo de fabricação em tecnologias submicrônicas. O foco deste trabalho de pesquisa é a busca de soluções regulares através do uso de matrizes de portas lógicas. A experimentação efetuada realiza a geração de uma matriz de portas NAND que viabiliza a implementação de equações lógicas mapeadas para redes de portas NAND e inversores, admitindo-se a parametrização do fanout máximo. Foi desenvolvida uma ferramenta de CAD, o MARTELO, que permite efetuar a geração automática de matrizes de portas lógicas, sendo que a versão inicial está voltada para a geração de matrizes com portas NAND em tecnologia CMOS. Os experimentos efetuados revelam que esta técnica é promissora, sendo apresentados alguns dos resultados obtidos.

Uma solução híbrida para a equação de Fokker-Planck dependente da energia usando a técnica da transformada de Laplace e diferenças finitas

Neste trabalho é obtida uma solução híbrida para a equação de Fokker-Planck dependente da energia, muito utilizada em problemas de implantação iônica. A idéia consiste na aplicação da transformada de Laplace na variável de energia e aplicação de um esquema de diferenças finitas nas variáveis espacial e angular desta equação. Tal procedimento gera um problema matricial simbólico para a energia transformada. Para resolver este sistema, procede-se a inversão de Laplace da matriz (sI+A), onde s é um parâmetro complexo, I a matriz identidade e A uma matriz quadrada gerada pela discretização das variáveis espacial e angular. A matriz A não é diagonalizável, desta forma, contorna-se este problema decompondo esta matriz na soma de outras duas, onde uma delas é diagonalizável. É gerado então um método iterativo de inversão, semelhante ao método da fonte fixa associado ao método de diagonalização, do qual o resultado fornecido são os valores para o fluxo de partículas do sistema. A partir disto pode-se determinar a energia depositada no sistema eletrônico e nuclear do alvo. Para validar os resultados obtidos faz-se a simulação de implantação de íons de B em Si numa faixa energética de 1keV a 50MeV, comparam-se os resultados com simulação gerada numericamente pelo software SRIM2003.

Determinação do comportamento carga-recalque de sapatas em solos residuais a partir de ensaios SPT

O comportamento carga-recalque de sapatas assentes em solos residuais é perfeitamente representado a partir da interpretação de provas de carga. No entanto, a execução de provas de carga em fundações reais ou em placa é indesejada pela maioria dos construtores e contratantes por razões de cronograma e custo. Por outro lado, sondagens de simples reconhecimento (SPT) são de rápida execução, exigidas pela normalização vigente de projeto e execução de fundações (ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS, 1996) e largamente utilizados na geotecnia nacional, sendo tal ensaio normalizado pela da norma NBR 6.484/2001. Deste modo, nesta pesquisa, pretende-se estabelecer uma metodologia semi-empírica para estimar o comportamento de sapatas assentes em solos residuais a partir de correlação estatística de resultados de sondagens de simples reconhecimento (NSPT) com a tensão admissível e os recalques sob tensão de trabalho, obtidos a partir de resultados de provas de carga em placa ou fundações reais assentes em solos residuais. Deste modo, a partir de retro-análise de um conjunto de resultados de provas de carga executados em diferentes solos residuais, encontrados na literatura, e o uso de conceitos de estatística, obteve-se diversas correlações que possibilitam estimar o comportamento tensão-recalque deste tipo de solo através dos resultados médios dos ensaios de SPT Como resultados destas correlações foram obtidas equações para a previsão de recalques e tensão admissível de fundações assentes sobre solos residuais, bem como do módulo de elasticidade de tais solos. As equações de estimativa de recalques são dependentes, além dos resultados de SPT, da geometria da fundação e da tensão aplicada. Já as equações para a estimativa da tensão admissível e módulo de elasticidade são dependentes somente dos valores médios do ensaio de SPT. E finalmente, para testar as equações desenvolvidas, foi realizada uma análise “Classe A” no comportamento carga-recalque de três provas de carga realizadas em locais diferentes, mostrando a eficiência da metodologia desenvolvida para a determinação da tensão admissível e dos recalques em sapatas sob a aplicação de carga de trabalho assentes em solos residuais.

Determinação da resposta dinâmica de feixes de condutores em linhas aéreas de transmissão

Até os dias de hoje, a abordagem estática tem sido a usualmente empregada para a avaliação da resposta de cabos de Linhas de Transmissão (LTs), apesar da resposta dinâmica ser, em muitos casos, reconhecidamente importante na avaliação do desempenho dos cabos. Para uma análise completa, não se pode desconhecer a natureza dinâmica da maioria dos fenômenos a que as LTs estão submetidas, sendo exemplos típicos as excitações mecânicas causadas pela ação do vento (carregamentos de baixa e alta freqüência) e a ruptura de cabos. Cabe salientar que mesmo a análise estática empregada tem sido simplificada, considerando normalmente apenas casos contemplados com soluções analíticas. A justificativa para emprego da análise estática baseia-se no menor esforço numérico exigido. Entretanto, essa justificativa já não se sustenta dado os grandes avanços na área computacional, possibilitando o estudo do desempenho de LTs na ocorrência de fenômenos que provocam carregamentos dinâmicos. Adicionalmente, há uma crescente demanda para que se obtenha um melhor entendimento de muitas questões relativas ao comportamento dinâmico das LTs, que possam explicar desempenhos observados e/ou sustentar o desenvolvimento de novas alternativas mais arrojadas de projeto e construção Entre os diversos métodos utilizados na engenharia estrutural, a integração direta das equações do movimento, através de métodos numéricos como as diferenças finitas centrais, se constitui numa poderosa ferramenta de cálculo. Esta ferramenta possibilita o tratamento de problemas envolvendo não linearidades geométricas e do material, como são os casos onde se avalia a resposta de cabos suspensos submetidos a carregamentos variáveis no tempo. Esta pesquisa objetiva a aplicação do método da integração direta das equações do movimento na análise de feixes de cabos de LT, quando submetidos à ação de carregamentos mecânicos variáveis no tempo, principalmente à excitação de ventos oriundos de fenômenos com natureza complexa (tormentas elétricas, por exemplo). É apresentado um método para determinar a resposta dinâmica de feixes que considera a interação entre o vento incidente e o movimento do condutor. A solução é obtida por integração numérica, no domínio do tempo, das equações de movimento de um modelo tridimensional não-linear discreto de feixe, as quais definem as forças nos espaçadores e cabos através de coeficientes aerodinâmicos obtidos experimentalmente. Também são apresentados modelos dos fenômenos meteorológicos mais comuns em nosso país (Brasil): tormenta extratropical (EPS) e tormenta elétrica (TS) Como ilustração, são apresentados exemplos de modelagem de vãos de LTs com condutor singelo e com feixes de condutores. Os exemplos demonstram a capacidade de avaliação de Estados Limites relacionados à distância relativa entre subcondutores, à estabilidade do feixe, à representação das suas propriedades e ao comportamento dinâmico, bem como aos carregamentos transmitidos às estruturas. A análise emprega conhecidas relações constitutivas para representar o comportamento tensão-deformação dos cabos. O enfoque utilizado possibilita a avaliação mais precisa de casos reais que ainda não podiam ser convenientemente tratados, além de permitir a extensão para estudos bem mais complexos, tais como feixes com disposições assimétricas.

"MGLTSmn" : aproximação angular "multigrid" em uma placa plana

Neste trabalho, desenvolve-se um método “multigrid” para a aproximação angular da solução da equação de transporte de partículas em uma placa plana, baseado na formulação LTSN com dependência contínua na variável angular. Para tanto, aplica-se a formulação LTSN sobre o conjunto de equações SN para determinar o fluxo angular de partículas nas N direções discretas referentes a uma malha grossa (N pequeno) e em seguida, usando os fluxos conhecidos, aplica-se a formulação LTSN com dependência angular contínua, para avaliar o fluxo angular nas M direções discretas referentes a uma malha fina (M grande). São apresentadas simulações numéricas que ilustram a capacidade desse método, denotado por MGLTSMN , no que diz respeito à redução do esforço computacional na aproximação da solução para problemas que requerem elevadas ordens de quadratura e alto grau de anisotropia.

Aerodinâmica de automóveis baseado no método dos contornos virtuais

Este trabalho desenvolve um método numérico para a solução de escoamentos bidimensionais em torno de geometrias automobilísticas utilizando o método de diferenças finitas. O código computacional resolve as equações de Navier-Stokes e de Euler para uma distribuição adequada dos pontos discretos na malha. O método de integração empregado baseia-se no esquema explícito de Runge-Kutta de 3 estágios para as equações da quantidade de movimento e no de sub-relaxações sucessivas para a pressão na base Gauss-Seidel. Utilizou-se a técnica dos contornos virtuais em coordenadas cartesianas para resolver o escoamento sobre uma geometria simplificada, com a superfície coincidente com a malha computacional, e uma geometria automobilística mais complexa (BMW). Para a certificação da técnica empregada, optou-se pela utilização da teoria do escoamento potencial e pela comparação com dados experimentais encontrados na literatura e outros coletados em túnel de vento em escala reduzida. Houve dificuldade nesta comparação devido à falta de artigos relativos às simulações numéricas de escoamentos sobre automóveis e na aplicação da técnica dos contornos virtuais em geometrias complexas. Os resultados foram satisfatórios, com boas perspectivas para trabalhos futuros, contribuindo assim para o desenvolvimento da área.

A base dinâmica na existência de autovalores duplos no modelo de Timoshenko para uma viga uniforme livre-livre

Considerando uma viga uniforme do tipo Timoshenko com condições de contorno livre-livre, Geist e McLaughlin em [8]apresentam uma condição necessária e suficiente que garante a existência de freqüências naturais duplas. Esta condição foi obtida usando a formulação espectral, método clássico encontrado na literatura, para as equações de quarta ordem desacopladas do modelo de Timoshenko. O método clássico requer a obtenção de um vetor constante com oito componentes para que a solução deste modelo seja conhecida. Segundo Claeyssen [2], [3], [4], [5], [6], a solução do modelo de Timoshenko pode ser obtida usando a base dinâmica gerada por uma resposta impulso-matricial fundamental. Este método permite encontrar a solução do modelo de Timoshenko usando as equações de segunda ordem acopladas. Além disso, para que a solução seja conhecida é necessário obter um vetor constante com quatro componentes. O objetivo deste trabalho é estudar a condição necessária e suficiente que garante a existência de freqüências naturais duplas, apresentada por Geist e McLaughlin, para uma viga uniforme do tipo timoshenko com condições de contorno livre-livre e verificar se é possível obter esta mesma condição quando é utilizada a base dinâmica para obter a solução deste modelo.

Movimento de partículas carregadas em fluidos ionizados : fundamentos matemáticos da teoria de eletroforese capilar

Neste trabalho, discutimos o movimento de uma macromolécula carregada em um fluido ionizado. A interação do campo elétrico é descrita pela equação de Poisson-Boltzmann acoplada às equações governantes para a dinâmica do fluido e às equações dinâmicas da partícula. Uma formulação fraca é introduzida no caso em que o domínio ocupado pelo fluido é finito e um teorema de existência de soluções fracas, local em tempo, é estabelecido. Dois modelos são considerados: fluxos não-estacionários e estacionários. No primeiro caso, a hidrodinâmica do sistema é governada pelas equações de Navier-Stokes, considerando-se um termo forçante relacionado ao potencial elétrico; no segundo caso, uma velocidade de deslizamento, a qual depende não linearmente sobre os potenciais, é introduzida como uma condição de contorno para um problema estacionário de Stokes. O caso de um fluido ocupando uma região infinita é também discutido supondo-se uma hipótese de aproximação sobre o campo elétrico.

Rafaella : um esquema para um novo paradigma de criptografia de chave pública

O presente trabalho apresenta um novo esquema de criptografia de chave pública baseado no emprego de funções para representar as mensagens original e cifrada. No esquema proposto – denominado Rafaella -, o processo de cifração consiste na aplicação de um deslocamento no argumento da função que representa a mensagem, de modo que se f(x) descreve a mensagem original, então f(x+z) representa a respectiva mensagem cifrada. O deslocamento z representa um número complexo que, no esquema proposto, representa a forma das chaves privadas dos participantes. A dificuldade da resolução do problema inversos concentra-se na obtenção das partes real e imaginária do deslocamento z, que pode ser efetuada através de método de força bruta, ou da resolução de um problema de contorno. A segunda alternativa envolve a resolução de equações diferenciais. Dentre os métodos disponíveis para a resolução de equações diferenciais, o emprego dos chamados grupos de Lie constitui, via de regra, a estratégia mais apropriada para a obtenção de soluções analíticas, que demandam menor tempo de processamento do que as formulações numéricas. Mesmo assim, a solução obtida através da utilização dos grupos de Lie requer elevado número de operações simbólicas.

Aproximações de Galerkin clássico e Galerkin/mínimos quadrados de escoamentos viscosos incompressíveis : aplicação à medição de fluxos

Desde a antigüidade a medição do escoamento dos fluidos tem sido uma marca de nossa civilização, ajudando a predizer a fertilidade das terras e o consumo d’água em fontes e aquedutos. Nos nossos dias, a área de medição de fluxo está bem estabelecida e ainda desperta grande interesse nas linhas de pesquisa da mecânica dos fluidos experimental e computacional. Em particular, o estudo da medição de fluxo com elementos intrusivos, tais como placas de orifício, é de grande interesse dado o preço baixo do medidor, e sua boa precisão quando comparada à sua simplicidade de instalação e manutenção. Esta dissertação tem como objetivo o estudo da aplicação de elementos finitos aos escoamentos de fluidos viscosos - via aproximação clássica de Galerkin e Galerkin/mínimos-quadrados (GLS) – com particular ênfase na aproximação das equações de Navier-Stokes incompressível no escoamento newtoniano através de um canal obstruído por uma placa de orifício. Inicialmente, são apresentadas as dificuldades do método de Galerkin clássico na aproximação de escoamentos incompressíveis; ou seja, através da simulação de escoamentos viscosos bem conhecidos - como o escoamento no interior de uma cavidade e através de uma expansão súbita - fica evidenciada a restrição imposta pela condição de Babuška-Brezzi quando da escolha dos subespaços aproximantes do campo de velocidade e pressão Como alternativa às patologias do método de Galerkin clássico, esta dissertação emprega a metodologia de Galerkin/mínimos-quadrados na simulação acima mencionada da placa de orifício, a qual permite o uso de elementos de igual-ordem para aproximar velocidade e pressão e capturar de maneira estável escoamentos sujeitos a altos números de Reynolds. Os testes computacionais realizados se apresentaram fisicamente realistas quando comparados com a literatura e dados experimentais, sendo todos desenvolvidos no Laboratório de Mecânica dos Fluidos Aplicada e Computacional (LAMAC) do Departamento de Engenharia Mecânica da Universidade Federal do Rio Grande do Sul.

Planejamento de perturbações para a identificação de modelos dinâmicos

A identificação de modelos é determinante no sucesso das modernas técnicas de controle avançado de processos. Um modelo para o sistema pode ser obtido através de modelagem rigorosa, baseada em equações governantes do sistema ou através da modelagem empírica e estimação de parâmetros. Embora mais rápida e fácil, a modelagem empírica necessita de alguns testes de identificação nos quais as variáveis manipuladas são variadas de modo que resultem em variações nas variáveis controladas. Os testes de identificação podem apresentar custos muito elevados tendo em vista que o sistema pode sair de seu ponto normal de operação, gerando produtos com folga de especificação. Este fato ocorre porque usualmente as perturbações aplicadas nas variáveis manipuladas nas indústrias de processos são independentes umas das outras, aumentando a duração do teste de identificação. Desta forma, neste trabalho foi desenvolvida uma nova metodologia de projeto de perturbações simultâneas para a identificação de modelos dinâmicos baseada na direcionalidade do sistema, com o objetivo de fornecer dados mais ricos para se capturar corretamente o comportamento multivariável do sistema e manter o processo no ponto de operação normal. As perturbações são projetadas conforme as características de um modelo simplificado do processo, ou pré-modelo. Este modelo inicial é obtido essencialmente de dados históricos de planta, selecionados através de uma sistemática análise de correlação desenvolvida neste trabalho A metodologia proposta é composta de duas partes: a primeira parte diz respeito à análise dos dados históricos de planta para obtenção de informações prelimirares as quais são utilizadas no planejamento de perturbações, tais como amplitude do ruído de medida, correlação entre as variáveis de processo, constante de tempo do sistema e matriz de ganhos. E a segunda parte consiste no cálculo da amplitude das perturbações baseado nos resultados da primeira etapa do planejamento. Para sistemas mal-condicionados verificou-se que as perturbações planejadas pela metodologia removem menos a planta de seu ponto de operação gerando resultados mais consistentes em relação às perturbações tradicionais. Já para sistemas bem-condicionados, os resultados são semelhantes. A metodologia foi aplicada em uma unidade piloto experimental e numa unidade de destilação da PETROBRAS, cujos resultados apontam pouca remoção dos sistemas do ponto de operação e modelos consistentes. A validação dos modelos também foi contemplada na dissertação, uma vez que foi proposto um novo critério de validação que considera a derivada dos dados de planta e a do modelo e não apenas os dados de planta e os dados da simulação das saídas do modelo.