No presente trabalho são apresentadas soluções analíticas das equações de difusão de nêutrons bidimensionais com dois grupos de energia, a saber, nêutrons rápidos e térmicos em uma placa com propriedades homogêneas. Alem disso, são resolvidos detalhadamente os problemas onde a placa homogênea é substituída por duas e quatro regiões, tornando-os não-homogêneos. A partir da aplicação da transformada de Laplace e da Técnica da Transformada Integral Generalizada (GITT), respectivamente, é resolvida em uma forma analítica o problema de autovalor resultante para o fluxo de nêutrons. No problema heterogêneo são usados filtros para homogenizar as condições de contorno não-homogêneas. Esta é a condição para a aplicação da GITT. Os três problemas mencionados acima são resolvidos aplicando primeiramente a GITT, o qual reduz a dimensão da equação de difusão, seguida da aplicação da transformada de Laplace, o qual reduz a ordem da equação. Deste procedimento, resulta um sistema de equações algébricas dependente das constantes de integração. 0 sistema é resolvido usando a técnica da eliminação de Gauss. Os fluxos transformados pela GITT são recuperados invertendo-se analiticamente a transformada de Laplace usando a expansão de Heaviside, os quais ainda dependem das constantes de integração. A partir da aplicação das condições de contorno e de interface (para os problemas não-homogêneos) obtém-se um sistema de equações algébricas homogêneas, de onde é determinado o fator de multiplicação efetivo Keff pelo método da bissecção. As constantes de integração são determinadas fazendo use da potencia prescrita da placa. Assim, os fluxos de nêutrons transformados pela GITT ficam determinados e os fluxos de nêutrons rápidos e térmicos são recuperados através da formula da inversa da GITT, usando a expansão do potencial. Resultados são comparados com a solução do método de diferenças finitas. / In the present work we present analytical solutions of the bi-dimensional neutron diffusion equation with two energy groups, i.e. fast and thermal neutrons in a sheet with homogeneous properties. Further we solve the detailed problem where the homogeneous sheet is substituted by two and four regions, rendering the problem a non-homogeneous one. Upon application of the Laplace transform and Generalized Integral Transform Tecnique (GITT), respectively, we solve in an analytical fashion the resulting eigenvalue problem for the neutron flux. In the heterogeneous problem, we use filter functions in order to homogenize the non-homogeneous boundary conditions. This is a condition for the application of GITT. We solve the three problems mentioned above applying first GITT, which reduces the dimension of the diffusion equation followed by the Laplace transform, which reduces the order of the equation. This procedure yields a non-homogeneous algebraic system depending on integration constants. The system is solved using the elimination technique by Gauss. The transformed fluxes by GITT are recovered upon inverting analytically the Laplace transform using Heaviside's expansion which depend still on the integration constants. Upon application of the boundary and interface conditions (for the non-homogeneous problem) one obtains a system of homogeneous algebraic equations, where we determine the effective multiplication factor keff by the bisection method. The integration constants are determined making use of the predefined power of the sheet. Thus the neutron fluxes transformed by GITT are determined and the fast and thermal neutron flux are recovered by the inverse formula of GITT, using the potential expansion. Results are compared to the solution by the finite difference method.

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).

No desenvolvimento de projetos de iluminação de ambientes, um dos objetivos que se destaca é a busca pelo conforto visual, que emprega metodologias de resolução conhecidas, como o Método dos Lumens e o Método Ponto a Ponto. A luz visível está contida no espectro da radiação térmica e, portanto, o fluxo luminoso pode ser relacionado ao fluxo de radiação térmica. Determinar as posições e as competências das fontes de luz necessárias na superfície de projeto ganha importância quando o comportamento, em termos de uniformidade ou de fluxo radiante, é especificado. O presente trabalho visa a estabelecer diferentes valores de fluxo em duas regiões distintas da superfície do projeto. Por meio do posicionamento das fontes de luz, é estabelecido um fluxo maior na região denominada principal e um fluxo menor na região denominada secundária. A modelagem matemática da radiação térmica (Método das Radiosidades) é aplicada ao projeto de iluminação, considerando-se as características da visão humana e o comportamento das fontes de luz. Na modelagem, é considerada uma cavidade retangular tridimensional com superfícies cinza e com condição de parede fria, na qual o poder emissivo das paredes é nulo. As fontes de luz são representadas por unidades de malha no teto. A relação de equações é resolvida por metodologia inversa, usando o algoritmo de Otimização Extrema Generalizada (GEO). Este algoritmo é classificado como um método de otimização estocástica de busca global para a resolução de sistemas considerados inicialmente mal condicionados. A posição e a potência das fontes luminosas são determinadas pela resolução do sistema de equações, de forma a proporcionar um fluxo de radiação duas vezes maior na região principal em relação à região secundária. A função objetivo do processo consiste em minimizar a diferença entre o fluxo desejado e os valores de fluxo de radiação incidente nas duas regiões da superfície de projeto. Em virtude das características de simetria do problema, a relação é estabelecida para apenas um quarto da cavidade. Assim, por exemplo, aplicar a metodologia com 9 fontes de luz a um quarto da região resulta em 36 fontes de luz em toda a cavidade. Os resultados mostram que é possível encontrar um arranjo de fontes de luz preestabelecendo-se duas condições de potência. / In the development of environmental illumination projects, one of the main goals to be achieved is the visual comfort, which is usually done by known methodologies, like the Lumens Method and the Point by Point Method. Since the visible light is contained in the spectrum of thermal radiation, the luminous flux can be related to the thermal radiation flux. The determination of the position and power of the light sources required by the design surface gains an higher importance whenever a behavior is specified, should it be in terms of uniformity or in therms of radiant flux. In this work, we describe a method that allows the establishment of different flux values in two distinct regions of the design surface, which are referred by the names main region and secondary region. Through the spatial arrangement of the light sources, the method sets a more intense flux in the main region and a less intense one in the secondary region. The mathematical model of thermal radiation, known as Radiosity Method, is applied to the illumination design, along with the characteristics of the human vision and the behavior of light sources. In this model, a rectangular three-dimensional cavity is considered. It has gray surfaces and exhibits the conditions of a cold wall, in which the emissivity power of the walls is null. The light sources are represented by a mesh unit in the ceiling. The system of equations is solved by inversemethodology, using the Generalized Extremal Optimization (GEO) algorithm. This algoritm is classified as being a stochastic optimization method of global search to solve systems that are initially considered ill-conditioned. By solving this system, the position and power of light sources can be determined, and this is done in such a way that the flux radiation in the main region is twice more intense then the one in the secondary region. The target function of the whole process is to minimize the difference between the desired flux and the incident flux radiation values for each one of the two design surface regions. We further explore the problem symmetry, solving the equation system for only a quarter of the cavity. This way, if the methodology is applied with nine light sources into a quarter of the region, the entire cavity will behave as if it has 36 light sources. Our results show that, given two prescribed conditions of power, it is possible to find an arrangement of light sources.

Este trabalho apresenta um estudo experimental do efeito da razão de bloqueio em um canal aerodinâmico no Número de Strouhal e o fenômeno da biestabilidade. É estudado o escoamento turbulento sobre um cilindro e também sobre dois cilindros dispostos lado a lado, onde fez-se uso de cilindros de vários diâmetros a fim de variar a razão de bloqueio no canal. Dois tipos de espaçamento entre os tubos são estudados. A técnica experimental consiste na medição de flutuações de velocidades em um canal aerodinâmico utilizando a técnica de anemometria de fio quente. Os dados obtidos da medição no canal aerodinâmico são tratados com o uso de ferramentas estatísticas, espectrais e de ondaletas. Para um cilindro, os resultados confirmam o aumento do número de Strouhal com o aumento da razão de bloqueio, na faixa do número de Reynolds utilizada. Razões de bloqueio mais baixas podem acarretar diminuição do número de Strouhal, também em função do número de Reynolds. Os resultados confirmam a existência do fenômeno da biestabilidade no escoamento sobre dois cilindros dispostos lado a lado, assim como a influência da razão de bloqueio sobre o fenômeno, sendo mais predominante em um dos espaçamentos estudados. Para razões de bloqueio altas, o fenômeno da biestabilidade pode não ocorrer ou pouco ocorre. Para razões de bloqueio mais baixas podemos ter também a não ocorrência do fenômeno em função das dimensões dos tubos serem mais reduzidas, fazendo com que a captação do fenômeno também seja reduzida. / This paper presents an experimental study of the effect of blockage ratio in a aerodynamic channel on Strouhal number and the bistability phenomenon. It studied the turbulent flow around a single cylinder and also on two cylinders arranged side-by-side, where use has been made of cylinders of various diameters to vary the blockage ratio of the channel. Two types of pitch to diameter ratio are studied. The experimental technique consists of measuring velocity fluctuations in aerodynamic channel using the technique of hot-wire anemometry. The data obtained from measuring the aerodynamic channel are treated with the use of statistical tools, spectral and wavelet analysis. For a single cylinder, the results confirm the increase of the Strouhal number with increasing blockage ratio in the range of Reynolds number used. Lower blockage ratio can cause a decrease of the Strouhal number, also a function of Reynolds number. The results confirm the existence of the bistability phenomenon in the flow on two cylinders arranged side-by-side, as well as the influence of blockage ratio on the phenomenon, being more predominant in one of the row spacings. For higher blockage ratios, the bistability phenomenon may not occur or occurs shortly. For lower blockage ratio, due to the smaller dimensions of the tubes the phenomenon may not occur on the capture of the phenomenon is also reduced with the present experimental technique.

Neste trabalho, é estudado o controle da transmissão do som numa placa retangular e fina. Para tanto, é encontrada a resposta dinâmica da placa, excitada por forças harmônicas pontuais e piezomomentos, obtida usando uma base não-clássica e uma análise modal. A radiação sonora emitida pela vibração da placa é encontrada. A potência sonora radiada pode ser calculada aplicando controle ativo diretamente na estrutura, na forma de uma entrada vibratória, uma vez conhecida a resposta na superfície da placa, obtendo-se uma signicativa redução analitica. Os piezocerâmicos, modelados como quatro momentos pontuais, são unidos a superfície da placa como atuadores. A potência sonora transmitida antes e depois do controle é comparada, usando diferentes número de atuadores. Uma estratégia clássica de controle linear quadrático (LQR) e empregada no contexto de um procedimento de otimização da posição dos atuadores do sistema.

O movimento da água em solos não saturados (zona vadosa) pode ser descrito pela equação de Richards, obtida pela combinação da equação de conservação de massa e da lei de Darcy-Buckingham. Usualmente essa equação é representada sob três formas básicas: h - baseada na carga de pressão, θ - baseada na umidade volumétrica e mista - baseada na combinação das duas variáveis. Essas equações são matematicamente equivalentes; porém, em geral, suas aproximações numéricas podem levar a resultados bastante diferentes. Neste trabalho foram validados e comparados modelos numéricos, implementados em um código computacional na linguagem C, para a aproximação da equação de Richards, nas suas três formas, em uma dimensão. Para isso, empregou-se o método de diferenças finitas centrado no espaço e atrasado no tempo. Visando melhorar a conservação de massa dos modelos utilizados e o tempo de execução do programa, principalmente em problemas que apresentam mudanças bruscas nas condições de contorno e/ou camadas estratificadas, foram também implementadas três estratégias para a escolha do passo de tempo, assim como técnicas para aproximar as não-linearidades nos pontos intermediários da malha computacional e naqueles situados nas interfaces das camadas. A validação dos modelos acima se deu através da comparação das soluções aproximadas com a solução semi-analítica de Philip e a totalmente analítica de Srivastava e Yeh para colunas de solos homogêneas e estratificadas. A partir de dados disponíveis na literatura avaliou-se também a precisão e a eficiência das metodologias estudadas nas: i) aproximações da condutividade e da capacidade hidráulica e ii) na aplicação das estratégias de passo de tempo, quanto à precisão dos resultados simulados, à conservação de massa e ao tempo de execução do programa. De acordo com o desempenho das metodologias analisadas foram apontadas as melhores opções para a solução de cada problema, indicando-se assim alternativas na implementação de modelos numéricos para o escoamento da água na zona vadosa.

Submitted by maria angelica Varella (angelica@sibi.ufrj.br) on 2018-02-06T13:29:23Z No. of bitstreams: 1 149225.pdf: 1936793 bytes, checksum: 7973f12e8760eecebb44bac866edcedd (MD5) / Made available in DSpace on 2018-02-06T13:29:23Z (GMT). No. of bitstreams: 1 149225.pdf: 1936793 bytes, checksum: 7973f12e8760eecebb44bac866edcedd (MD5) Previous issue date: 1978-11 / Este trabalho tem como objetivo o desenvolvimento de um algoritmo computacional, na área de Mecânica dos Fluidos, para projeto de instalações e sistemas hidráulicos. Consiste, principalmente, no desenvolvimento de um método de projeto de sistemas hidráulicos, compatível com a utilização de métodos computacionais em relação às variáveis de projeto, características do movimento de fluidos incompressíveis e Newtonianos, em regime estacionário. Com tal objetivo, desenvolveu-se um algoritmo numérico, utilizando-se equações de Darcy-Weisbach, Hagen-Poiseuille e HazenWilliams, para o cálculo da perda de carga, e de Colebrook, para o cálculo do fator de atrito. A resolução da equação de Colebrook foi pelo método iterativo de Newton-Raphson. As características para a determinação dos coeficientes de resistência em elementos do sistema hidráulico, tais como válvulas, reduções, derivações, etc., compõe um segundo algoritmo para as perdas de carga secundá- rias. Este algoritmo numérico, correspondente a um sistema iterativo convergente, atende aos vários anéis e aos padrões de projeto dos sistemas. O programa global, composto das rotinas de cálculo, utiliza também um esquema iterativo de correção dos incrementas das vazões, em cada anel do sistema, ao mesmo tempo em que considera as alterações necessárias resultantes da utilização dos elementos secundários. Os elementos secundários são calculados segundo as equações convencionais da literatura existente, na forma de energia cinética, de modo a facilitar a introdução de novos parâmetros quando estes se fizerem necessários. Além disso, o programa atende especificações restritivas em pressão de ejeçao nos vários pontos de intersucção dos anéis, possibilitando uma rápida análise de performance do comportamento do sistema, corrigindo automaticamente os diâmetros dos tubos necessários. / This work aims at the development of a algorithm within the Fluids Mechanical area to the computational installations and hydraulic systems design. It mainly consists of the development of a method of hydraulic systems design compatible to the use of computational methods connected with the design variation and with the incompressible fluids' move characteristics and Newtonian one, in a stationary mood. Held to such purpose, a numerical algorithm was developed through the use of equations from Darcy – Weisbach Hazen Poiseuille and Hazen-Williams, to figure out the loss of charge, and from Colebrook, to figure out the factor of friction. Colebrook's equation solution was gotten by the use of iterative method. The characteristics to establish Newton-Raphson's coefficients of resistence within the hydraulic systems elements, such as valves, reductions, derivations and so one, compose a second algorithm to the secondary lasses of charge. This numerical algorithm corresponds to a convergent iterativa system and fits on several loops and on patterns of the systems design. The whole schedule, composed of usual calculation, also makes use of a small iterative scheme for correcting the increments of drainages in each system loop, while, at the sarne time, it controls the necessary changes which come from the use of the secondary elements. The secondary elements are calculated according to the conventional equations from the existing literature, in a cinetic energy mood, so that to facilitate the inclusion of further details whenever they are necessary. Besides, the schedule is okay for restrictive specifications in pressure of ejection on different points of the loops' intersuction and for this reason the schedule makes possible a quick study on the system's behavior acting by putting in order, automatically, the diameter of tubes.

Made available in DSpace on 2016-06-02T19:55:23Z (GMT). No. of bitstreams: 1 2572.pdf: 799090 bytes, checksum: 4abf7acc4d34d2e0c1659207352d35aa (MD5) Previous issue date: 2009-08-04 / Universidade Federal de Sao Carlos / Clavulanic acid (CA) is a β-lactam antibiotic produced industrially by submerged cultures of the filamentous bacterium Streptomyces clavuligerus. Due to its importance, several studies on the purification of clavulanic acid have been conducted at the Department of Chemical Engineering of the Federal University of São Carlos DEQ-UFSCar. As a natural sequence of these works, a need has been identified for greater knowledge about the CA adsorption process by mathematically modeling the incorporation of the phenomenon of CA degradation in the sites of resin through a first order reaction. Because amino acids stand out among the principal molecules that are found in fermented broth, it has become imperative for the development of the process to evaluate purification in terms of the amino acids present during the purification process. Some of these amino acids, e.g., tyrosines, present physicochemical properties similar to those of clavulanic acid. The process of batch absorption can be simulated by a mathematical model, providing kinetic and transport parameters involved in absorption. A mathematical model of the clavulanic acid purification process in a fixed-bed absorption column was developed to determine the rupture curve. The modeling allowed for the determination of several important process parameters, such as transport parameters and the intrinsic kinetic constants of absorption. A sensitivity analysis demonstrated that the mass transfer coefficient strongly influences the results of the model. The simulation of the process led to the conclusion that the process is limited both by the degradation stage during absorption and by the mass transfer, indicating that both factors must be taken into account in the analysis of this type of process in a fixed-bed column. It was also found that the degradation of CA in the resin sites is a factor that does not allow the rupture curve to reach the value of C/Co equal to 1, which explains the low yield of the CA process during absorption. This work also aimed to apply the sequencing method to solve the partial differential equations that represent the mathematical model which describes the adsorption of CA in Amberlite IRA-400 resin and to compare it with the finite difference method which is more commonly employed. The application of the proposed methodology in the clavulanic acid purification process in a packed column showed a good performance, with a rupture curve very similar to that obtained by the finite difference method, thus demonstrating the efficiency of the proposed method. / O ácido clavulânico (AC) é um antibiótico beta-lactâmico produzido industrialmente por culturas submersas da bactéria filamentosa Streptomyces clavuligerus. Devido a esta importância, alguns trabalhos sobre purificação do ácido clavulânico vêm sendo desenvolvidos no DEQ-UFSCar. Como seqüência natural desses trabalhos, surgiu a necessidade de um conhecimento maior do processo de adsorção do AC por meio da modelagem matemática a incorporação do fenômeno de degradação do AC nos sítios da resina, através de uma reação de primeira ordem. Como dentre as principais moléculas que se encontram no caldo fermentado destacam-se os aminoácidos, tornou-se imprescindível para o desenvolvimento do processo avaliar a purificação em termos dos aminoácidos presentes durante o processo de purificação. Alguns destes aminoácidos, como por exemplo, as tirosinas apresentam propriedades físico-químicas próximas à do ácido clavulânico. O processo de adsorção em batelada pode ser simulado através de um modelo matemático, obtendo-se parâmetros cinéticos e de transporte envolvidos na adsorção. Uma modelagem matemática do processo de purificação do ácido clavulânico em coluna de leito fixo foi desenvolvida para obter a curva de ruptura. A modelagem permitiu a determinação de vários parâmetros importantes do processo tais como parâmetros de transporte e as constantes cinéticas intrínsecas de adsorção. Análise de sensibilidade demonstrou que a coeficiente de transferência de massa exerce forte influência nos resultados do modelo. Pela simulação do processo foi possível concluir que o processo é limitado tanto pela etapa de degradação durante a adsorção quanto pela transferência de massa, indicando que ambos devem ser levados em conta na análise do processo deste tipo de processo em coluna de leito fixo. E, também, verificou-se a degradação do AC nos sítios da resina é um fator que não permite que a curva de ruptura não atinja o valor de C/Co igual a 1 o que explica o baixo rendimento do processo do AC durante o processo de adsorção. Este trabalho teve também como finalidade aplicar o método da seqüência na resolução das equações diferenciais parciais que representam o modelo matemático que descreve a adsorção de AC na resina Amberlite IRA- 400 e compará-lo com o método de diferenças finitas que é mais comumente utilizado. A aplicação da metodologia proposta ao processo de purificação do ácido clavulânico em coluna recheada mostrou seu bom desempenho, cuja curva de ruptura é bem próxima da curva de ruptura obtida pelo método de diferenças finitas, mostrando, assim, a eficiência da metodologia proposta.

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.