O tema deste trabalho versa sobre a produção de anéis de centragem de alto-falantes de uma empresa de alto-falantes do Rio Grande do Sul. O objetivo geral é melhorar o processo de fabricação de anéis de centragem de alto-falantes.Como objetivo específico, quer-se conhecer melhor, através de estudos de Desdobramento da Função Qualidade (QFD), Projeto deExperimentos e Análise de Regressão, os parâmetros de processo que afetam a qualidade dos mesmos,e apontar soluções para amelhoria do processo, reduzindo refugos e retrabalhos. Em função dos resultados, foi elaborado um plano de ação demelhorias do processo, partindodas conclusões do experimento e visando-se melhorar a situaçãoatual do processo na empresa.Na apresntação dos resultados, testa-se o trabalho no intuito de melhorar o processo na empresa em questão.

A maioria dos métodos de síntese e sintonia de controladores, bem como métodos de otimização e análise de processos necessitam de um modelo do processo em estudo. A identificação de processos é portanto uma área de grande importância para a engenharia em geral pois permite a obtenção de modelos empíricos dos processos com que nos deparamos de uma forma simples e rápida. Mesmo não utilizando leis da natureza, os modelos empíricos são úteis pois descrevem o comportamento específico de determinado processo. Com o rápido desenvolvimento dos computadores digitais e sua larga aplicação nos sistemas de controle em geral, a identificação de modelos discretos foi amplamente desenvolvida e empregada, entretanto, modelos discretos não são de fácil interpretação como os modelos contínuos pois a maioria dos sistema com que lidamos são de representação contínua. A identificação de modelos contínuos é portanto útil na medida que gera modelos de compreensão mais simples. A presente dissertação estuda a identificação de modelos lineares contínuos a partir de dados amostrados discretamente. O método estudado é o chamado método dos momentos de Poisson. Este método se baseia em uma transformação linear que quando aplicada a uma equação diferencial ordinária linear a transforma em uma equação algébrica evitando com isso a necessidade do cálculo das derivadas do sinais de entrada e saída Além da análise detalhada desse método, onde demonstramos o efeito de cada parâmetro do método de Poisson sobre o desempenho desse, foi realizado também um estudo dos problemas decorrentes da discretização de sinais contínuos, como por exemplo o efeito aliasing decorrente da utilização de tempos de amostragem muito grandes e de problemas numéricos da identificação de modelos discretos utilizando dados com tempos de amostragem muito pequenos de forma a destacar as vantagens da identificação contínua sobre a identificação discreta Também foi estudado um método para compensar a presença de offsets nos sinais de entrada e saída, método esse inédito quando se trata do método dos momentos de Poisson. Esse trabalho também comprova a equivalência entre o método dos momentos de Poisson e uma metodologia apresentada por Rolf Johansson em um artigo de 1994. Na parte final desse trabalho são apresentados métodos para a compensação de erros de modelagem devido à presença de ruído e distúrbios não medidos nos dados utilizados na identificação. Esses métodos permitem que o método dos momentos de Poisson concorra com os métodos de identificação discretos normalmente empregados como por exemplo ARMAX e Box-Jenkins.

nRPN, nRPNSTAT e nRPNDYN são novos índices, propostos nesta dissertação, que quantificam o grau de não-linearidade total, estática e dinâmica, respectivamente, de sistemas dinâmicos não- lineares. Estes novos índices permitem responder se é necessária a utilização de controladores não-lineares ou se simples controladores lineares apresentarão desempenho satisfatório. Além disso, se um controlador não- linear se fizer necessário, quais características este deve apresentar. Nesta dissertação é realizado um estudo acerca das metodologias para determinação do grau de não- linearidade de plant as industriais. Entre as encontradas, a que apresenta as melhores potencialidades é estudada em detalhes, tendo seus atributos e deficiências explorados. As deficiências foram sanadas e um novo conjunto de índices é proposto (nRPN, nRPNSTAT e nRPNDYN), o qual, além de apresentar resultados conclusivos em termos quantitativos, também traz novas informações valiosas para a escolha do controlador adequado, como a parcela estática e dinâmica da não- linearidade. A metodologia proposta é aplicada e testada em um conjunto bastante amplo e complementar de possíveis sistemas dinâmicos não-lineares (tais como: pH, reatores químicos e coluna de destilação de alta pureza) onde o comportamento nãolinear já foi explorado na literatura O desenvolvimento e o estudo do comportamento dinâmico de uma planta laboratorial de cinco tanques acoplados são apresentados. Este sistema possui dinâmica e ganhos variáveis, descontinuidades, inversão no sinal do ganho multivariável e diferentes graus de acoplamento, ilustrando de forma simples vários possíveis comportamentos dinâmicos encontrados em plantas industriais. Devido a esta grande versatilidade a planta foi construída e testada, confirmando os comportamentos previstos nos estudos teóricos.

Neste trabalho, desenvolveram-se modelos matemáticos simplificados para o cálculo de duas variáveis importantes no estudo da extrusão de polímeros: vazão mássica (M&) e pressão na saída da extrusora (Psaída), em função das propriedades dos materiais e das condições de operação do equipamento. Podem-se utilizar esses modelos como ferramentas simples para a definição de critérios de ajustes que se devem fazer em uma extrusora de parafuso único para obter-se o desempenho desejado quando se alimenta o equipamento com um novo material. Para desenvolverem-se os modelos simplificados, utilizaram-se dados experimentais da extrusão de poliestireno (PS) e de polipropileno (PP), bem como resultados preditos por um programa computacional de simulação de extrusão disponível comercialmente. Mediram-se os dados experimentais de vazão mássica e de pressão na saída da extrusora em um equipamento de parafuso único de 45 mm de diâmetro. Realizaram-se esses testes, variando-se a velocidade de rotação do parafuso de 70 a 100 rpm para ambos os polímeros. No primeiro conjunto de simulações, utilizou-se o simulador Flow 2000 (Compuplast Inc.) para ajustarem-se os valores preditos de M& e de Psaída aos dados obtidos experimentalmente através da estimação dos fatores de atrito barril-polímero tanto para o PP quanto para o PS. Posteriormente, realizou-se um planejamento de experimentos, do tipo fatorial fracionado , para obter-se um segundo conjunto de simulações, considerando-se as propriedades dos materiais (reológicas e térmicas) e as condições de operação da extrusora (velocidade de rotação do parafuso e perfil de temperatura nas zonas de aquecimento da extrusora) como fatores de investigação. Com as novas simulações no Flow 2000, ajustaram-se os parâmetros dos modelos simplificados aos valores de vazão mássica e de pressão na saída da extrusora preditos no simulador. Elaboraram-se os modelos simplificados levando-se em conta as interações entre os fatores cujos efeitos consideraram-se significativos nas análises de variância (ANOVA). Obteve-se um modelo linear com 37 termos para o cálculo da vazão mássica e um modelo linear com 41 termos para o cálculo da pressão na saída da extrusora. Posteriormente, aplicou-se uma técnica de regressão multivariável para selecionar apenas os termos importantes dessas 1402IV2− XVI equações, conduzindo a um modelo linear com 10 termos para o cálculo da vazão mássica e a um modelo com 6 termos para o cálculo da pressão na saída da extrusora. Conseguiu-se boa concordância entre os dados experimentais e os valores preditos quando se aplicaram os modelos simplificados.

O uso racional de utilidades (energia elétrica, vapor e água) representa atualmente o grande desefaio para garantir a competitividade e sustentabilidade das indústrias. O custo da energia e os limites ambientais associados à atividade são estratégicos para perpetuar o negócio industrial, ainda mais se depender da vida do setor de refino de petróleo. O trabalho proposto apresenta a minimização do custo da geração de energia em uma empresa petroquímica, que possui produção de energia elétrica e vapor através de um sistema de cogeração. Neste caso, há um equilíbrio a ser buscado entre suprimento de energia de potência e térmica e isto não pode ser rapidamente definido a partir de procedimentos ou malhas de controle localizadas. A aplicação de um modelo de otimização é proposto, baseando-se nas demandas de energia e flexibilidades existentes e rapidamente expondo a configuração que minimiza o custo de produção de vapor e energia elétrica (custo da central termoelétrica). Esta meta foi obtida aplicando duas ferramentas: uma planilha para resgatar, reconciliar, tratar, organizar dados e em seguida exportar informações para a segunda ferramenta, o modelo de otimização no qual o problema da central termoelétrica é construído e solucionado através de algoritmos adequados. Após, o trabalho também envolveu o estudo de 14 casos. O resultado obtido consolida o potencial para otimização: é possível reduzir em cerca de 7% (pelo menos) o custom da unidade termoelétrica, com pequenas intervenções de manutenção e ajustes operacionais. Estes dados são associados com a possibilidade de redução da geração de vapor e demandas de energia elétrica e também com a redução do uso de válvulas redutoras de pressão, na maioria dos casos. / The rational use of utilities (electric energy, steam, and water) represents nowadays the great challenge to assure the competitiveness and sustainability of industries. The energy cost and the associated environmental limits are strategic to perpetuate the industrial business, even more if it depends upon the oil sector life. The proposed work presents the minimization of the cost of steam and power generation in a petrochemical company, which has production of electric energy and steam by cogeneration system. In this case, there is a balance to achieve between work and heat supply and this cannot be readily defined by heuristics or localized control loops. The application of an optimization model is proposed based on needs of energy demands and existing flexibilities, and readily exposing the scenario that minimizes the power and high pressure steam level production. This goal was obtained applying two tools: one flow sheet for rescue, reconciliate, treat, organize data and then export information for the second, the optimization model in which the thermoelectric problem is constructed and then solved by proper algorithms. After, this work also develop and study of 14 cases. The results obtained consolidate the potential for optimization: it is possible to reduce in about 7% (at least) the cost of the thermoelectric unity, with little maintenance and operational adjustments. These figures are associated with the possibility of reducing the steam generating and power needs and also reducing the use of letdown stations, in most cases.

O eixo diferencial é amplamente utilizado em veículos automotores para transmitir a potência do motor para as rodas, sendo esta transmissão realizada por um par de engrenagens conhecidas como pinhão e coroa. O eixo é dito diferencial porque permite que uma roda do veículo gire com rotação diferente da outra, característica necessária nas curvas, por exemplo. Um dos principais objetivos na montagem de um eixo diferencial é proporcionar a coroa e ao pinhão o mesmo posicionamento relativo no qual o par foi usinado, quanto mais próximo da cota nominal melhor será o desempenho do eixo, diminuindo ruído, vibração e aumentando a sua vida útil. A tarefa de montar o pinhão com um posicionamento preciso em linhas de produção é uma tarefa difícil devido a quantidades de peças intermediárias, deslocamento do produto montado (pré-carga), e a tolerância da cota de posicionamento (±0,030mm). Para melhorar este processo e torná-lo viável para produção em linhas de montagem em grande escala a empresa Muri Linhas de Montagem desenvolveu em conjunto com a Dana Traction Tecnologies um sistema de montagem em quatro passos: a) medição dos componentes individuais com e sem pré-carga; b) seleção do anel espaçador para compensar as diferenças dimensionais dos componentes; c) montagem das peças; d) verificação final do conjunto montado. Este sistema já foi desenvolvido e implementado de forma parcialmente satisfatória, pois o processo de seleção do anel espaçador apresenta falhas e desvios de montagem não conhecidos, o que leva a introdução de uma variável de ajuste no processo. Este sistema já foi desenvolvido e implementado de forma parcialmente satisfatória, pois o processo de seleção do anel espaçador apresenta falhas e desvios de montagem não conhecidos, o que leva a introdução de uma variável de ajuste no processo. O problema é que esta variável de ajuste não é constante necessitando ajuste manual e contínuo durante a produção. No decorrer deste trabalho será explicado o funcionamento do eixo diferencial com os seus requisitos de montagem, após será demonstrada a linha de montagem desenvolvida com os processos de medição e cálculo do anel espaçador. Por terceiro será realizada uma avaliação teórica e experimental para buscar as variáveis de processo desconhecidas e por quarto e último será proposto um novo método de medição e montagem mais preciso do que a método atual. / The differential axle is used in the automotive industry to transmit the engine power to the vehicle wheels, this transmission is done by 2 gears called pinion and ring gear. The axle is called differential because permit that the wheels turn with different speed, this is necessary in the curves for example. One of the main objectives to assembly a differential axle is assembly the pinion gear in the same relative position that the pair was machined, as close as possible to active the best performance, reduce noise, vibration and increase the product life. Assembly the axle with precision in assembly lines is difficult because there are internal components positioned between the pinion and ring gear, there is a product displacement (preload) and also a position tolerance (±0,030mm). To improve this process and became viable to produce in assembly lines on large flow volumes the company Muri Linhas de Montagem developed in partnership with Dana Traction Tecnologies one assembly system that works in four steps: a) measure all intermediary components with and without pre-load; b) select the correct spacer to compensate the product dimensional differences; c) assembly the parts; d) check the final assembly. This system is already developed and work partially satisfactory because there is a failure in the process to select the spacer with some unknown deviations, to compensate those deviations is used a variable in the calculations. The problem is this variable is not constant and demand a manual adjust during the assembly process. During this paper work will be explain how works the differential axle with the assembly requirements, after will be show the assembly line developed with the measurement process and spacer calculation procedure. In the next part will be performed an evaluation to find the unknown process variation, in the last part will be proposed a new methodology to measure and assembly the axle with more accuracy comparing with the old methodology.

O controlador PID é o algoritmo de controle mais difundido nas unidades industriais em todo o mundo. Além de estar disponível em praticamente todas as plataformas de controle comerciais, ele é um algoritmo robusto, de fácil entendimento, e capaz de prover performance satisfatória para uma grande variedade de processos industriais. No entanto, seu desempenho para a grande maioria dos casos está muito aquém do esperado. Estimativas mostram que apenas 20% das malhas de controle industriais estejam funcionando de forma adequada, diminuindo a variabilidade do processo. Assim, devido a necessidade de melhorar o desempenho das malhas de controle industriais, foi desenvolvida uma metodologia complexa, porém intuitiva, através da qual é possível sintetizar controladores do tipo PID de qualquer parametrização, baseada em um problema de otimização no domínio da freqüência, que minimiza a diferença entre a resposta do sistema em malha fechada frente à uma variação do tipo degrau unitário, e uma resposta específica desejada, representada pela função desempenho alcançável (FDA). Diversas alternativas foram propostas para automatizar todo o procedimento e facilitar a sua utilização, porém foi deixado um grau de liberdade para o usuário alterar o nível de desempenho desejado, através da variação de um único parâmetro que indica a velocidade da resposta da FDA. Adicionalmente, foi feita uma série de recomendações sobre como variar o nível de desempenho desejado, respeitando as restrições inerentes do processo e os limites de estabilidade para diferentes tipos de casos. Para processos mais complexos, onde o comportamento dinâmico não é representado de forma satisfatória por um único modelo nominal, foi desenvolvida uma metodologia alternativa, adequada para um sistema multi-modelos, caracterizado por um conjunto de modelos lineares representativos de diferentes pontos de operação. Esta metodologia é uma extensão da original, porém, com a resolução de um problema de otimização multiobjetivo com formulação do tipo minmax, onde minimiza-se o máximo desvio entre a resposta em malha fechada de cada modelo considerado e uma única FDA. Neste tipo de abordagem, é estabelecido um compromisso entre o desempenho atingido em todos os pontos de operação, quando controlado por um único controlador PID. / The PID controller is the most widespread used control algorithm in the industrial plants of the whole world. Besides it is available in practically all commercial control system, it is a robust algorithm, easily understood, able to provide satisfactory performance for a number of industrial processes. Nevertheless, its performance for most of the cases is below than the desired. Estimates show that just 20% of the industrial control loops are found to work well, decreasing the processes variability. Thus, due to the need of improving the performance of the industrial control loops, a complex, but intuitive methodology was developed, by which it is possible to synthesize PID controllers of any parameterization, based on a frequency domain optimization problem, that minimizes the difference between the unit step response of the closed loop system and a desired specific response, represented by the attainable performance function (APF). A lot of alternatives were purposed to turn the whole procedure automatic and simple, however, there is a degree of freedom, by which the practitioner can vary the desired performance level, changing only one parameter, which indicates the APF response speed. In addition, a number of recommendations about how to vary the desired performance level were done, respecting the inherent process constraints and the stability boundaries for different kind of cases. For more complex processes, where the dynamic behavior is not well represented by only one nominal model, an alternative methodology was developed, suitable for a multi-model system, that is specified by a set of linear models, each one representing a different operating point. This methodology is the sequence of the original one, but formulated as a min-max multiobjective optimization problem, where it is minimized the maximum deviation between the closed loop response of each considered model and only one APF. In this kind of approach, it is possible to have a trade-off among the performance reached in all operating points, when controlled by only one PID controller.

Metodologias para processamento de sinais têm sido amplamente pesquisadas com o objetivo de se extrair informação útil a partir de dados de processo. Nesta categoria incluem-se os filtros digitais, os quais atenuam os ruídos ou erros aleatórios, componentes de alta freqüência do sinal. Técnicas de atenuação de ruído como aproximações polinomiais, filtros de Fourier e Ondaletas recebem destaque dentre as técnicas de filtragem. Os filtros das Ondaletas, por exemplo, baseiam-se nos componentes de freqüência de um sinal, realizando a transformada deste sinal de modo a eliminar as freqüências acima de um limite especificado. Este sinal suavizado (tendência), de ruído atenuado, também pode ser processado de modo a obter a sua derivada. A diferenciação numérica de sinais digitais tem muitas aplicações em processamento de sinais analíticos, devido às suas propriedades, e em várias técnicas de engenharia. Um sinal suavizado desempenha papel importante na identificação de estados estacionários por exemplo, que encontra grande aplicação em análise de processos, otimização, identificação de modelos e reconciliação de dados. Como estas aplicações requerem dados sob condição estacionária, ou muito próxima dela, um método eficiente de detecção de estacionários se faz necessário. A derivada de um sinal ruidoso fica mascarada pelo próprio ruído e não apresenta de forma clara os pontos de inflexão desejados para a análise do sinal. Faz-se necessário, portanto, tratar o sinal com um filtro de modo a reduzir ou eliminar o ruído antes de processá-lo ou, então, obter a tendência que este sinal apresenta. Neste trabalho busca-se comparar algumas metodologias comumente utilizadas para filtragem de sinais de processos e obtenção de tendências, e analisa-se a aplicação de uma destas metodologias para análise e identificação de estados estacionários. Uma metodologia multi-variável simples e eficiente baseada em PCA é proposta, tomando como referência o método proposto por JIANG et al. (2003), que determina um índice de estado estacionário pontual. Para validação do método, estudaram-se alguns sistemas industriais e os resultados obtidos são apresentados, mostrando-se satisfatórios. / Methodologies for signal processing are being widely used to extract useful information from process data. In this category are included the digital filters, which atenuate noise or randomic errors (i.e., high frequency components of the signal). Techniques to atenuate noise as polinomial approximation, Fourier filters and Wavelets are commonly used. The Wavelet filters, for example, are based on the high frequency components of the signal, performing the wavelet transform on it to eliminate de frequencies above a specified level. This smooth signal, with less noise (trend), can be also processed to obtain its derivative. Numerical diferentiation of digital signals has many applications in analytic signal processing, due to its properties, and in several engineering techniques. A smooth signal plays an important role in steady-state identification, for example, in applications for process analysis, process optimization, model identification and data reconciliation. As these applications require data under steady-state condition, or very close to it, an efficient method for steady-state detection is necessary. The derivative of a noisy signal is hidden by the noise itself and does not show clearly the inflexion points desired for the signal analysis. Thus it is necessary to filter the signal to reduce or even eliminate the noise before processing it. In this work, some commonly used filter techniques are used and compared, and the application of one of them to steady-state detection is analised. An easy and efficient multivariable metodology based on PCA is proposed, according to the method presented by JIANG et al. (2003), which determines a pointwise steady-state index. To validate the method, some industrial systems were studied and the obtained results are presented, showing that they are satisfatory.

Neste trabalho foi desenvolvido um ferramental para prensagem hidráulica de esferóides de alumina baseado no processo de fabricação de tijolos de alumina e placas cerâmicas para revestimento, pois até a realização deste, o processo conhecido no mercado para a produção de esferas de alumina, era através da prensagem conhecida como isostática. Estas esferas e esferóides tem o objetivo de servir como corpos moedores para as indústrias cerâmicas e de mineração para o beneficiamento de matérias primas fazendo a redução das partículas do material sólido, sendo que estas esferas e esferóides são introduzidas no interior de um equipamento chamado de moinho de bolas. Então, foi construído um mecanismo adaptativo que faz com que uma prensa hidráulica, que até o momento fazia apenas a compactação de material com formatos quadrados e retangulares com uma determinada espessura, possa produzir compactados com um formato esférico. O desenvolvimento se deu primeiramente na confecção de um modelo com uma cavidade de compactação obtido através de uma prensa hidráulica manual, sendo possível com esta obter-se esferóides para os testes laboratoriais de densidade aparente e de desgaste, sendo estas duas as principais premissas para o bom desempenho destes corpos moedores nos interiores dos moinhos de desagregação de matéria prima. Com os resultados de densidade aparente e de desgaste que foram alcançados com os esferóides produzidos a partir do modelo de uma cavidade, foi possível construir um protótipo do ferramental, só que agora com oito cavidades de compactação, e foi instalado este protótipo em uma prensa hidráulica automática para verificar se o seu desempenho em real situação de produção industrial, conseguia-se os mesmos resultados obtidos com o modelo de uma cavidade. Com isso, constatou-se a viabilidade da utilização do método proposto do ferramental para a prensagem hidráulica de esferóides de alumina, considerando que todos os ensaios experimentais e industriais realizados, se comportaram com igualdade ou superioridade às esferas produzidas no método tradicional que é através da prensagem isostática. / This paper describes a work where a tool for hydraulic pressing of spheroids made of alumina was developed based on the process of making bricks and making ceramic boards for covering. This work was carried on because so far, the known process in the market for production of spheres made of alumina was through a pressing known as isostatic. These spheres and spheroids have the function to serve as grinder bodies in the ceramic and mining industries for the improvement of raw material resulting in the reduction of particles of the solid material, that is, these spheres and spheroids are put inside an equipment called balls mill. So, an adaptive mechanism was built to make the hydraulic press, which so far performed only the compaction of material of square and rectangle shapes with a determined thickness, producing compacted ones of a spherical shape. At first, the development started with the construction of a model containing one cavity of compaction obtained through a manual hydraulic press, which provided spheroids for the lab tests of ostensible density and abrasion, which are considered the main premises for the good performance of these grinding bodies inside the mills for dissolution of raw material. With the effects as of density apparent and as of wear and tear than it is to have been ranging with the spheroids produced from the template from a crevice, and compared with the values of the balls manufacture traditionally , he went feasible assembly um target system from the tool was built containing not one, but eight cavities of compaction and it was installed inside an automatic hydraulic press to verify if its performance in real situation of industrial production, could reach the same results obtained with the model with one cavity. By doing that, the viability of the use of the tool proposed method for hydraulic pressing of spheroids made by alumina was proved, taking to account that all the experimental and industrial tests performed behaved equally as, or superiorly than the spheres produced in the traditional method, which is through isostatic pressing.

A modelagem termodinâmica apresenta um papel fundamental no desenvolvimento e simulação de processos químicos. Esta modelagem utiliza como ferramenta modelos matemáticos conhecidos como equações de estado. Com o desenvolvimento da mecânica estatística aplicada, novos tipos de equações de estado foram propostas para uso prático em engenharia. Uma delas, PC-SAFT (Perturbed Chain-Statistical Associating Fluid Theory), tem mostrado um grande potencial de aplicação, sendo capaz de vencer a grande maioria das limitações apresentadas por outros modelos e podendo ser aplicada para vários tipos de sistemas contendo fluidos puros e misturas, incluindo casos nos quais estão presentes polímeros. O presente trabalho teve por objetivo o desenvolvimento de ferramentas para utilização da equação de estado PC-SAFT para o cálculo de propriedades termodinâmicas e de equilíbrio de fases. Tais ferramentas deverão ser utilizadas para realização deste tipo de cálculo para polímeros e suas soluções em um simulador de processos que está sendo desenvolvido com colaboração do nosso grupo de pesquisa. As propriedades termodinâmicas para as quais foram implementadas rotinas de cálculo foram a pressão, energia livre de Helmholtz residual, pressão de vapor, ponto crítico de substâncias puras, coeficiente de fugacidade e potencial químico residual. Em relação ao equilíbrio de fases, foram implementadas rotinas para o equilíbrio líquidovapor de misturas binárias e multicomponentes, sendo que para misturas binárias ainda foi implementada uma rotina para misturas contendo polímeros. No caso dos polímeros, foram também implementadas rotinas para o cálculo do equilíbrio líquido-líquido de soluções contendo polímeros polidispersos. Foram propostas metodologias para que uma inicialização adequada de todos os processos considerados que requerem o uso de métodos iterativos e, também, uma mudança de variável para melhorar o desempenho da estimação de parâmetros do modelo. Todas as rotinas foram desenvolvidas no software MatLab e a validação das rotinas desenvolvidas foi feita através da comparação de suas predições com dados experimentais e numéricos da literatura. A grande quantidade de informação gerada nos testes também permitiu identificar os principais problemas numéricos associados a algumas das etapas de cálculo. / The thermodynamic modeling plays a fundamental role in the development and simulation of chemical processes. This modeling is based on the use of mathematical models known as equations of state. With the development of applied statistical mechanics, new types of equations of state were proposed for use in practical engineering. PC-SAFT (Perturbed Chain-Associating Statistical Fluid Theory) is among these models and has shown great potential of application, being able to overcome most of the limitations presented by other models and applicable to various types of systems containing pure fluids and mixtures, including polymeric systems. The main goal of this work was the development of tools for the application of the PC-SAFT equation of state for the calculation of thermodynamic properties and phase equilibrium. These tools are meant to be used for performing this kind of calculation for polymers and their solutions in a process simulator which is being developed with collaboration of our research group. The thermodynamic properties for which calculating routines were implemented were pressure, residual Helmholtz free energy, vapor pressure, critical point of pure substances, fugacity coefficient, and residual chemical potential. Regarding the phase equilibria, the following routines were implemented: vapor-liquid equilibrium of binary and multi-component mixtures of non-polymeric materials, vapor-liquid equilibrium of binary mixtures containing polymers, and liquid-liquid equilibrium for solutions of polydisperse polymers. Routines have also been implemented for the determination the parameters of the PC-SAFT equation of state for pure substances. Methodologies for the definition of appropriate initialization values for the calculation routines requiring the use iterative processes were also developed and a change of variable was proposed to improve the performance of the parameter estimation problem. All routines were implemented in the MatLab 5.0 and their validation was accomplished by comparing our predictions with experimental and numerical data from the literature. The vast amount of information generated in the tests performed also allowed identifying the main numerical problems related to some of the calculation processes.