Este trabalho apresenta metodologias para a construção e a sintonia de estimadores de estados não-lineares visando aplicações práticas. O funcionamento de um estimador de estados não-linear está calcado em quatro etapas básicas: (a) sintonia; (b) predição; (c) atualização da matriz de covariância de estados; (d) filtragem e suavização dos estados. As principais contribuições deste trabalho para cada uma destas etapas podem ser resumidas como segue: (a) Sintonia. A sintonia adequada da matriz de covariância do ruído de processos é fundamental na aplicação dos estimadores de estado com modelos sujeitos a incertezas paramétricas e estruturais. Sendo assim, foi proposto um novo algoritmo para a sintonia desta matriz que considera dois novos métodos para a determinação da matriz de covariância dos parâmetros. Este algoritmo melhorou significativamente a precisão da estimação dos estados na presença dessas incertezas, com potencialidade para ser usado na atualização de modelos em linha em práticas industriais. (b) Predição. Uma das etapas mais importantes para a aplicação do estimador de estados é a formulação dos modelos usados. Desta forma, foi mostrado como a formulação do modelo a ser usada em um estimador de estados pode impactar na observabilidade do sistema e na sintonia das matrizes de covariância. Também são apresentadas as principais recomendações para formular um bom modelo. (c) Atualização da matriz de covariância dos estados. A robustez numérica das matrizes de covariância dos estados usadas em estimadores de estados sem e com restrições é ilustrada através de dois exemplos da engenharia química que apresentam multiplicidade de soluções. Mostrou-se que a melhor forma de atualizar os estados consiste na resolução de um problema de otimização sujeito a restrições onde as estimativas fisicamente inviáveis dos estados são evitadas. Este também preserva a gaussianidade dos ruídos evitando que estes sejam mal distribuídos. (d) Filtragem e suavização dos estados. Entre as formulações estudadas, observou-se também que a melhor relação entre a acuracidade das estimativas e a viabilidade de aplicação prática é obtida com a formulação do filtro de Kalman estendido sujeita a restrições (denominada Constrained Extended Kalman Filter - CEKF), uma vez que esta demanda menor esforço computacional que a estimação de horizonte móvel, apresentando um desempenho comparável exceto no caso de estimativas ruins da condição inicial dos estados. Como uma solução alternativa eficiente para a estimação de horizonte móvel neste último caso, foi proposto um novo estimador baseado na inclusão de uma estratégia de suavização na formulação do CEKF, referenciado como CEKF & Smoother (CEKF&S). / This work presents approaches to building and tuning nonlinear state estimators aiming practical applications. The implementation of a nonlinear state estimator is supported by four basic steps: (a) tuning; (b) forecast; (c) state covariance matrix update; (d) states filtering and smoothing. The main contributions of this work for each one of these stages can be summarized as follows: (a) Tuning. An appropriate choice of the process-noise covariance matrix is crucial in applying state estimators with models subjected to parametric and structural uncertainties. Thus, a new process-noise covariance matrix tuning algorithm is presented in this work which incorporates two new methods for the parameter covariance matrix computation. The algorithm has improved significantly the state estimation accuracy when the presence of such uncertainties, with potential to be applied in on-line model update in industrial practice. (b) Forecast. One of the most important stages in applying state estimators is the used model formulation. In this way, it has been shown that the model formulation to be used in state estimator can impact on the system observability and noisecovariance matrices tuning. In this work it is also presented the main recommendations to formulate an appropriated model. (c) State covariance matrix update. The numerical robustness of the state covariance matrices used in unconstrained and constrained state estimators is illustrated by two chemical engineering examples tending to multiple solutions. It has been shown that the best technique to update the states consists in solving an optimization problem subjected to constraints, since it prevents from physically unfeasible states. It also preserves the noise gaussianity preventing from bad noise distribution. (d) States filtering and smoothing. Among the studied formulations, it was also noticed that the better relationship between performance and practical application is obtained with an extended Kalman filter formulation subjected to constraints (called Constrained Extended Kalman Filter - CEKF) because it requires small computational effort than MHE with comparable performance, except in case of poor guesses of the initial state. As an efficient solution for moving horizon estimation in the last case, it was proposed a new estimator based on the addition of a smoother strategy into the CEKF formulation, referred as CEKF & Smoother (CEKF&S).

É conhecida ao longo da história da humanidade a grande quantidade de acidentes vinculados à ação de cargas cíclicas. Há mais de 100 anos, metodologias levam em conta o fenômeno de fadiga em projetos de novos componentes. As metodologias existentes para avaliar a vida em fadiga são empíricas e têm sido aplicadas com sucesso na determinação da vida do componente em fadiga quando o mesmo é submetido a carregamento uniaxial ou a carregamento multiaxial proporcional. Mas a evidencia experimental tem mostrado que quando a solicitação é multiaxial e não proporcional, as leis antes mencionadas deixam de ser adequadas. Dessa forma os critérios clássicos utilizados não preveem corretamente a vida do componente. Neste contexto o presente trabalho avalia algumas das metodologias de fadiga multiaxial não proporcionais disponíveis na bibliografia especializada. Os resultados obtidos com estes métodos são comparados entre sim e com os resultados obtidos utilizando a metodologia de fadiga multiaxial proporcional. Para comparar as metodologias citadas é apresentado um exemplo de solicitação simples, e também a análise do componente de uma máquina agrícola solicitada por um histórico de tensões típicas, obtida a partir do teste da máquina em campo. Finalmente são apresentadas as conclusões sobre as metodologias utilizadas e sobre a importância de levar em conta a não proporcionalidade de um estado de tensões multiaxial oscilantes. / It is known throughout human history the large amount of accidents linked to action of cyclic loading. For over 100 years, methodologies take into account the fatigue phenomenon of new components projects. Existing methodologies for assessing the fatigue life are empirical and have been successfully applied in the determination of fatigue life of the component when it is subjected to uniaxial loading or proportional multiaxial loading. But the experimental evidence has shown that when the request is multiaxial and not proportional, the laws mentioned above are no longer adequate. Thus the classical criteria used not correctly predict component life. In this context, this paper analyzes some of the methodologies of multiaxial fatigue non proportionate available in the relevant literature. The results obtained with these methods are compared with the results so obtained and using the proportional multiaxial fatigue methodology. To compare the methods mentioned is an example of simple request, and also the analysis component of an agricultural machine requested by a history of typical strains, obtained from the machine field testing on a test track. Finally conclusions on the methodologies used and the importance of taking into account the non-proportionality of a state of oscillating multiaxial stresses are presented.

Taludes escavados em filitos muitas vezes estão sujeitos à instabilidade, principalmente quando há outras famílias de descontinuidades presentes na estrutura. Este problema está presente na mina de Candiota, onde foi identificado, em um dos taludes, potencialidade para ruptura por tombamento. Este trabalho aborda uma discussão sobre métodos analíticos e numéricos aplicado a tombamento primário, com objetivo de avaliar o fator de segurança do talude em estudo por meio de método numérico por elementos finitos, utilizando o Phase2 da Rocscience, e métodos analíticos. Também são apresentadas doze simulações para avaliar se o movimento de tombamento em modelagens numéricas obedece a um efeito de escala. Com uso de métodos numéricos dois modelos foram construídos: um contínuo equivalente e outro contínuo (mas com as descontinuidades inseridas no modelo). No modelo com as descontinuidades a rocha foi tratada como elástica por meio do critério de Hoek-Brown e as descontinuidades foram tratadas como elasto-plásticas por meio do critério de Mohr-Coulomb. Os resultados das modelagens da mina de Candiota foram semelhantes e mostraram fatores de segurança em níveis de estabilidade, sendo que no modelo contínuo (mas com descontinuidades inseridas) se observou de forma mais clara as tensões cisalhantes induzidas nas descontinuidades quando ocorre o movimento de tombamento. O uso dos métodos analíticos mostrou-se ineficientes para o problema apresentado. Os resultados das doze modelagens sugerem que o aumento do espaçamento das descontinuidades influencia nos fatores de resistência dos taludes e na forma do movimento de tombamento. / Slopes in phyllites are often subject to instability, especially when there are other families of discontinuities in the structure. This problem is present in Candiota mine, which was identified in one of the slopes potential to toppling. This work show a discussion of analytical and numerical methods applied to toppling, to evaluate the slope safety factor being studied by numerical method finite element using the Phase2 of Rocscience, and analytical methods. Also shown are twelve simulations to evaluate if the toppling in numerical modeling follows a scaling effect. With use of numerical methods two models were built: an equivalent continuous and another solid (but with discontinuities inserted in the model). In the model with the discontinuities the rock was treated as elastic by the Hoek-Brown criterion and discontinuities were treated as elastic-plastic by the Mohr-Coulomb criterion. The results of Candiota mine were similar and showed safety factors in stability levels, and in the continuous model (with discontinuities) was observed more clearly the shear stress induced in discontinuities when the movement of overturning occurs. The use of analytical methods proved inefficient for the problem presented. The results of the twelve modeling suggests increasing the spacing of the discontinuities of slope influences the resistance factor and the form of the toppling movement.

Esta dissertação apresenta um estudo da caracterização numérica do comportamento mecânico de um músculo esquelético estriado sujeito à combinação de diferentes contrações. A partir dos mecânismos fisiológicos da contração do tecido muscular esquelético, são apresentadas as principais características da produção de força do músculo em diferentes tipos de contrações (isométrica, concêntrica, excêntrica) e suas combinações. A representação numérica da resposta mecânica do músculo esquelético foi investigada por diferentes autores. Dos diversos modelos encontrados na literatura, alguns foram escolhidos para a continuidade deste estudo, segundo sua capacidade representativa e facilidade de implementação em métodos numéricos de solução para grandes deformações. Os modelos disponíveis na literatura foram implementados e avaliados para diferentes sequências de contrações, sendo também realizado um ajuste de parâmetros para dados experimentais de contrações isométricas e isométricas-excêntricas-isométricas. Nenhuma das propostas avaliadas apresentou resultados satisfatórios. A principal deficiência dos modelos foi a incapacidade de reproduzir a dependência do histórico de carregamentos, ou seja, da combinação de contrações. Desta maneira, uma nova proposta de modelo constitutivo foi desenvolvida e implementada para testes uniaxiais e em um código de elementos finitos para testes tridimensionais. A partir destas implementações, foram executados testes numéricos para diferentes sequências de contrações, avaliando a capacidade representativa do modelo proposto neste trabalho. O modelo proposto apresentou bons resultados para contrações isométricas e a combinação de contrações com diferentes níveis de alongamento, ou encurtamento, a uma mesma velocidade. A principal contribuição deste modelo é a capacidade de representar a resposta associada à fadiga muscular e o ganho, ou perda de força, observados experimentalmente. / This work present a striated skeletal muscle numerical characterization subjected to different contractions. From the physiological contraction mechanism of the skeletal muscle tissue the force produced in different kinds of contraction (isometric, concentric and eccentric), and their combination, are presented. The numerical representation of these response was investigated by different authors. From the diverse models found in literature, some were selected to be studied according to their representation capability and the implementation ease in numerical methods for large strains. The models available in literature were implemented and evaluated for different contractions combinations and a parameter identification for experimental results of two isometric contraction and a isometric-eccentric-isometric contraction. None of the proposed models presented satisfactory results. The main deficiency of these models were the incapability of reproducing the loading history dependence, in other words, the contraction combination. In this way, a new material model was proposed and implemented to uniaxial and tridimensional finite element method tests. From these, different contractions sequences evaluated the proposed model representation capabilities. The proposed model present good results to isometric contraction, as well as contraction combination with different stretch, or shortening, level in the same speed. The main contribution of this model is the capability of represent the the response associated to the muscle fatigue and the force gain or loss, experimentally observed.

A pesquisa voltada para a retroanálise de taludes de mineração está se tornando cada vez mais importante. As dimensões das atuais cavas de mineração visam sempre obter a melhor relação entre custo, benefício e segurança. Isso exige que por vezes o trabalho ocorra próximo destes limites. Neste trabalho o estudo de retroanálise consiste em um talude misto composto por aluvião e rocha em uma mineração de carvão localizada no município de Arroio dos Ratos, RS. O objetivo foi realizar um estudo de retroanálise onde se utilizou dois métodos de análise método dos Elementos Finitos (MEF) e o método do Equilíbrio Limite (EL). A metodologia adotada foi dividida em três fases: uma etapa de campo, para coleta do material e dados geométricos do talude, uma etapa de laboratório para caracterizar todo o material e a última etapa computacional onde foram gerados os modelos e simulações numéricas através dos softwares comerciais da empresa RocScience, sendo eles o SLIDE e o RS2. Os resultados encontrados mostraram que a porção do talude composta por aluvião apresentava seis camadas intercaladas de areias bem graduadas e argilas de média a alta plasticidade. Em relação as simulações, observou uma diferença entre o MEF e EL devido ao grande número de camadas presentes no talude e a interação entre elas. Foram criados quatro diferentes cenários para observar o comportamento do talude, com presença e ausência de água. O MEF mostrou uma superfície de ruptura próxima a encontrada porem com a presença de água o Fator de Segurança (FS) sofreu uma grande depreciação. O EL considerou uma superfície de ruptura satisfatória, nas simulações com presença de água FS sofreu uma depreciação mínima mesmo com máxima saturação de água no talude. Apesar das diferenças ambos os métodos são confiáveis mas devem ser usados de forma adequado ao tipo de problema existente. Situações mais simples e com taludes mais homogêneos o EL trona-se mais indicado, para problemas de maior complexidades e interações entre camadas o MEF se mostrou satisfatório e com boa resposta. O resultados através do MEF se aproximaram mais da situação real encontrada em campo. / A research focus on retro analysis on mine slopes is becoming increasingly important. The dimensions of current mining caves aim to always obtain the best cost benefit ratio as well as the highest safety level. This demands that the work sometimes is done within these limits. In this paper the study of retro analysis is made up of a mixed slope composed of alluvial and rocks, in a coal mine in Arroio dos Ratos, RS. The goal was to carry out a retro analysis study using two analysis methods: The Finite Element Method (FEM) and Limit Equilibrium analysis (LE). The methodology adopted was split into three phases: a field stage for the sample collection and slopes geometric , a laboratory stage to classify all the material and a final computacional stage were model and numerical simulations were created by way of commercial software, Slide and RS2, by RocScience. The results showed that the portion of the slope composed of alluvium had six layers interspersed with well graded sands and high plasticity clays. The simulations observed a difference between the FEM and LE due to the number of layers present in the slope and the interaction among them. Four different scenarios were created in order to observe the behavior of the slope. The FEM presents surface failure similar to one found in the field. The LE always takes into consideration the overall critical surface which is not observed in the field, it also show high security factor levels above one 1.3 that does not reflect a real situation since the slope is collapsed. In spike of the differences, both methods are reliable but they must be used according to the existing problems. For simpler situations with homogeneous slopes LE is recommended while move complex problems and interactions among lawyers calls for FEM.

Recentemente, surgiu na literatura uma solução analítica das equações de cinética pontual que considera a reatividade como função do tempo, utilizando o método da decomposição. O presente trabalho dá um passo a frente, considerando as equações de cinética pontual em conjunto com efeitos de realimentação de temperatura. Mas, primeiro, uma breve introdução do modelo de cinética pontual e dos aspectos relevantes desta abordagem são apresentados. O trabalho prossegue acrescentando a realimentação de temperatura como uma perturbação de primeira ordem na equação reatividade, e através de manipulações algébricas, o conjunto de equações de cinética pontual passa a ser expresso como uma equação diferencial não linear de segunda ordem. Esta equação é, então, resolvida pelo método de decomposição, ou seja, expandindo as variáveis dependentes como séries infinitas, construindo-se então um sistema recursivo que permite calcular cada um dos termo destas séries. A não linearidade é tratada utilizando os polinômios Adomian. Os resultados aqui obtidos são comparados com a literatura, apresentando variações percentuais máximas da ordem de [0,1%]. Faz-se uma breve análise da convergência e da estabilidade da solução, usando um método baseado no critério de Lyapunov. / An analytical solution of the point kinetics equations to calculate reactivity as a function of time by the decomposition method has recently appeared in the literature. The present work goes one step forward, by considering the neutron point kinetics equations together with temperature feedback effects. But first, a brief introduction of the point kinetics model and the relevant aspects of this approach are presented. The work proceeds adding the temperature feedback as a first order perturbation in the reactivity equation, and using algebraic manipulation, the set of point kinetics equations becomes expressed as a non linear second order differential equation. This equation is then solved by the decomposition method, that is, expanding the dependent variables as infinite series, building a recursive system that allows the evaluation of each term of these series. The non linearity is treated using the Adomian polynomials. The results obtained are compared with literature, with maximum percentage changes of about [0,1%]. A brief analysis of the convergence and stability of the solution is made, using a method based on the Lyapunov criterion.

SOUSA, R. C. Análise e comparação de sete métodos numéricos utilizados na determinação dos parâmetros da curva de Weibull aplicados aos dados de velocidade do vento coletados na cidade de Paracuru-CE e Camocim-CE. 2011. 83 f. Dissertação (Mestrado em Engenharia Mecânica) – Centro de Tecnologia, Universidade Federal do Ceará, Fortaleza, 2011. / Submitted by Marlene Sousa (mmarlene@ufc.br) on 2016-03-31T16:59:17Z No. of bitstreams: 1 2011_dis_rcsousa.pdf: 6262897 bytes, checksum: 65414f8d08a48b7497a547e1d037cc3a (MD5) / Approved for entry into archive by Marlene Sousa(mmarlene@ufc.br) on 2016-03-31T19:18:15Z (GMT) No. of bitstreams: 1 2011_dis_rcsousa.pdf: 6262897 bytes, checksum: 65414f8d08a48b7497a547e1d037cc3a (MD5) / Made available in DSpace on 2016-03-31T19:18:15Z (GMT). No. of bitstreams: 1 2011_dis_rcsousa.pdf: 6262897 bytes, checksum: 65414f8d08a48b7497a547e1d037cc3a (MD5) Previous issue date: 2011-02-23 / This dissertation deals with the analysis and comparison of 7 (seven) numerical methods for the assessment of effectiveness in determining the parameters for the Weibull distribution, using wind speed data collected in Camocim and Paracuru cities, in the period from August 2004 to April 2006, obtained by the Department of Infrastructure of the Ceará state. Using the method of analysis of variance, RMSE test and the chi-square and analysis of power for test to compare the proposed methods, this study aims to determine which method was most effective in determining the shape parameters of Weibull distribution for the available data, aiming to establish criteria to a better utilization of wind power in the of Ceará state, which is in national prominence in the use of renewable sources for electricity generation. / Esta Dissertação de Mestrado trata da análise e comparação de 7 (sete) métodos numéricos para avaliação da eficiência na determinação dos parâmetros para a distribuição de Weibull, através de dados de velocidade de vento coletados nas cidades de Paracuru-Ce e Camocim-Ce, no período compreendido entre agosto de 2004 a abril de 2006, obtidos junto à Secretaria de Infraestrutura do Ceará – SEINFRA. Utilizaram-se as metodologias de análise de variância, A raiz quadrada dos erros quadráticos médios, do teste do qui-quadrado e da análise da potência, para comparação dos métodos propostos. Este trabalho visou à determinação de qual método fora o mais eficiente na determinação dos parâmetros da curva de distribuição de Weibull para os dados disponíveis, objetivando se estabelecer critérios para melhor aproveitamento do potencial eólico do estado do Ceará, que tem ocupado posição de destaque nacional na utilização de fontes renováveis para geração de energia elétrica.

O objetivo desta tese é introduzir uma nova metodologia para a cinética bidimensional multi- grupo de difusão de nêutrons em reatores nucleares. A metodologia apresentada usa uma aproximação polinomial em um domínio homogêneo retangular com condições de contornos não homogêneas. Como ela consiste em uma série de Taylor truncada, sua estimativa de erro varia de acordo com o tamanho do retângulo. Os coeficientes são obtidos principalmente pelas suas relações com o termo independente, que _e determinado pela equação diferencial. Estas relações são obtidas apenas pelas condições de contorno, e é demonstrado serem linearmente independentes. Um esquema numérico é feito para assegurar uma rápida convergência. Estes procedimentos feitos para um retângulo homogêneo são feitos para construir soluções para problemas de autovalor e dependentes do tempo de geometria ortogonal global com parâmetros seccionalmente constantes pelo método iterativo SOR. O autovalor dominante e sua autofunção são obtidos pelo método da potência no problema de autovalor. A solução para casos dependentes do tempo usam o método de Euler modificado na variável tempo. Quatro casos-teste clássicos são considerados para ilustração. / The objective of this thesis is to introduce a new methodology for two{dimensional multi{ group neutron diffusion kinetics in a reactor core. The presented methodology uses a polyno- mial approximation in a rectangular homogeneous domain with non{homogeneous boundary conditions. As it consists on a truncated Taylor series, its error estimates varies with the size of the rectangle. The coefficients are obtained mainly by their relations with the independent term, which is determined by the differential equation. These relations are obtained by the boundary conditions only, and these relations are proven linear independent. A numerical scheme is made to assure faster convergence. The procedures done for one homogeneous rectangle are used to construct the solution of global orthogonal geometry with step{wise constant parameters steady state and time dependent problems by the iterative SOR algo- rithm. The dominant eigenvalue and its eigenfunction are obtained by the power method in the eigenvalue problem. The solution for the time dependent cases uses the modi ed Euler method in the time variable. Four classic test cases are considered for illustration.

O semieixo automotivo é um componente usualmente metálico pelo qual o torque é transmitido do diferencial para as rodas do veículo. Diferentemente dos semieixos convencionais, os fabricados em materiais compósitos possuem algumas variáveis de projeto que devem ser levadas em consideração, como a orientação das fibras, a sequência de empilhamento, a espessura e o número de camadas. Neste trabalho, cilindros compósitos para aplicação em semieixos de transmissão foram processados por enrolamento filamentar com o objetivo de substituir peças em aço convencional por compósito laminado em fibra de carbono com matriz epóxi. Para tal, os parâmetros da aplicação do semieixo em compósito, como frequência natural, torque critico de flambagem e capacidade de transmissão de torque foram considerados. Foi realizado um ensaio preliminar de torção e foram selecionados três laminados: [±22/±45], [±89/±45] e [±45]2. Cilindros de carbono/epóxi foram fabricados pelo processo de enrolamento filamentar e ensaios de compressão radial e longitudinal dos cilindros foram realizados. Análises mecânicas foram realizadas por meio de abordagens analíticas e numéricas, sendo que a última foi baseada no método dos elementos finitos. Os resultados analíticos foram comparados com os numéricos e experimentais, apresentando uma boa correlação entre eles. O cilindro com o laminado [±45]2 apresentou o melhor desempenho de resistência sob torção. Para torque critico de flambagem, assim como em compressão radial e longitudinal, o laminado [±45]2 apresentou um desempenho satisfatório para a aplicação, apesar de não ser o melhor desempenho dentre os laminados. Conclui-se que o laminado [±45]2 é o que apresenta melhor desempenho, em geral e que a redução de massa do semieixo desenvolvido em compósito utilizando este laminado, em relação ao de aço é de 47%. / The automotive half shaft is usually a metallic-based component in which the torque is transmitted from the differential to the wheels. Unlike conventional half shafts, in composite half shafts there are several parameters to be improved, namely fiber orientation, stacking sequence, thickness and number of layers. In this study, composite cylinders for application in half shafts have been processed by filament winding in order to replace conventional steel parts by laminated composites with carbon fiber reinforced epoxy composites. For that, it must satisfy several parameters of the composite half shaft, such as critical speed, critical buckling torque and load carrying ability. A preliminary torsional test has been performed and three laminates were selected: [±22/±45], [±89/±45] and [±45]2. Carbon/epoxy cylinders were manufactured by filament winding process and radial and longitudinal compression tests were performed. Mechanical analysis has been accomplished by using analytical and numerical approaches, since the last one is based on the finite element method. Analytical results were compared with the numerical and experimental ones, presenting a good correlation between them. The cylinder [±45]2 has the best performance under torsional loading. For critical buckling torque, as well as radial and longitudinal compression, the laminate [±45]2 presented a satisfactory performance for the application, but not the best performance within the laminates. The main conclusion is that the laminate [± 45]2 has the best performance and the weight reduction of the composite half shaft developed using this laminate is of 47%, when compared to a steel-based shaft.

A presente discussão e uma extensão do modelo de cinética pontual de nêutrons, onde a reatividade e decomposta em termos de contribuição de curtas e longas escalas de tempo. A primeira representa o controle operacional do reator, enquanto a segunda e devido a alteração da composição química do combustível nuclear, como consequência do burn-up. E um primeiro passo em uma nova direção, uma vez que considera os efeitos dos principais venenos na cinética de nêutrons, ou seja, Xenônio-135 e Sam ario-149. O modelo proposto consiste em um sistema de equações não-lineares acoplado para a densidade de nêutrons, para os precursores de nêutrons atrasados e para as cadeias de decaimento dos venenos produtos de fissão. O sistema de equações e resolvido através de um método de decomposição, que expande os termos não-lineares em uma série infinita, obtendo um sistema recursivo, onde a inicialização da recursão e uma equação linear homogênea e os passos de recursão subsequentes consideram contribuições não-lineares como termo fonte constru dos em passos de recursão anteriores. A construção hierárquica do modelo também e realizada, onde graus espaciais de liberdade são considerados. São apresentados casos de estudos com várias estruturas temporais afim de mostrar a robustez da abordagem atual para este tipo de problema. / The present discussion is an extension to Neutron point kinetics models, where the reactivity is decomposed in a short and a long term contribution. The rst one represents operational reactor control, whereas the second one is due to the change of the chemical composition of the nuclear fuel as a consequence of burn-up. This is a rst step into a new direction where we consider only the e ects of the principal neutron poisons on neutron kinetics, i.e, Xenon-135 and Samarium-149. The proposed model consists in a system of coupled nonlinear equations for the neutron density, the delayed neutron precursors and the neutron poison decay chains. The equation system is solved using a decomposition method, which expands the non-linear terms in an in nit series, obtaining a recursive system, where the recursion initialization is a homogeneous linear equation and the subsequent recursion steps consider the non-linear contributions as source terms constructed from previous recursion steps. A hierarchical construction of the model is also performed, where spatial degrees of freedom are considered. We present case studies with severe time structure in order to show the robustness of the present approach for this kind of problems.