As equações da cinétiica pontual de um reator nuclear térmico são integradas numericamente, utilizando um método matricial de continuação analitica. Essas equações são essencialmente não-negativas e possuem um autovalor dominante vinculado à reatividade do sistema. Também, descrevem-se os métodos de Hansen e Porsching.

[pt] Neste trabalho é apresentada a formulação do modelo matemático de um elemento com número reduzido de graus-de-liberdade para a análise linear tridimensional de vigas. O elemento possui quatro graus-de-liberdade por ponto nodal, sendo três deslocamentos e uma rotação local da seção reta do elemento em relação ao eixo longitudinal da viga. Os efeitos de rotação local da seção reta do elemento em relação ao eixo longitudinal da viga. Os efeitos de rotação reta devido à flexão são considerados na formulação através da condição de não haver deformações transversais por cisalhamento. As condições de contorno de engastamento e de compatiblidade entre elementos contíguos são garantidas através de um procedimento de penalidade, introduzindo modificações no funiconal da energia potencial total da estrutura, assegurando a compatibilidade entre os elementos. Análises de algumas estruturas simples são apresentadas para demonstrar a eficiência da formulação, e suas soluções comparadas com resultados analíticos ou numéricos. A formulação é proposta de forma a ser facilmente estendida a análises mais gerais que incluam efeitos dinâmicos e de não linearidade do material. / [en] Throughout this thesis the formulation of the mathematical modelo f a general tri-dimensional beam element with reduced degrees-of-freedom is presented. The element is a four node displacement based finite element with three linear and one locally referrer rotational displacements at each node. Bending effects are taken into account through the imposition of non-transversal shear condition. Coupling between adjoining elements anda t the element clamped restraint is asssured by a penalty procedure imposed at the common node location. The analysis of some simple structures are presented and the solutions compared with other analytical or numerical results. The formulation hás been proposed so that it can be generally extended to material non-linear and dynamic analysis.

[pt] Este trabalho tem como objetivo apresentar uma alternativa para se analisar e avaliar opções sobre DI Futuro. Para tanto, faremos uso da teoria clássica sobre derivativos, e em particular, do modelo sugerido por Black [2] para a avaliação de opções sobre futuros de commodities. O contrato em questão, não possui solução analítica devido ao comportamento não linear do seu pay- off. A teoria define que a equação diferencial que descreve o comportamento do preço do ativo é função do ativo objeto. Neste trabalho, algumas simplificações foram assumidas, face a não adoção de um modelo estocástico que determine o comportamento futuro da taxa livre de risco, neste caso definida como um parâmetro determinístico do modelo. É fato de que tal simplificação não invalida os resultados, pelo contrário, McConnell e Schwartz [17] mostram que a relação custo benefício em se adotar modelos mais sofisticados não compensa frente aos resultados obtidos quando praticidade e ganhos são comparados. De posse da equação diferencial que governa o comportamento do preço do derivativo, se faz presente a necessidade de se usar um procedimento numérico - Método de Diferenças Finitas Explícito (MDFE). / [en] The main objective of this paper is to describe an alternative model to value Brazilian DI Future option. And so, we will make use of the classical derivatives theory, in particular, to the model introduced by Black for options on commodities future contracts. For such instrument, the analytical solution is not possible to be obtained due to the non-linear formulation of the pay-off (Risk Neutral Valuation). The theory defines the differential equation that describes the asset price behavior, in this case the financial operation agreed, as function of the underlying variables that govern its behavior. In the present work some simplifications had been carried through, regarding the non-adoption of a stochastic model to represent the future behavior of the risk-free rate, being defined as a deterministic parameter in the model. One must bear in mind that such simplification does not invalidate the results; on the contrary, McConnell e Schwartz [17] shows that the trade-off between the practicability and the profit in term of the results makes questionable the use of the more sophisticated model. Having the differential equation that governs the behavior of the derivative contract price, a numerical procedure is carried out - Explicit Finite Differences Method (EFDM).

[pt] Diversos acidentes envolvendo dutos submarinos já aconteceram no Brasil como o do Gasoduto PCR-01/Lubnor e o da plataforma P-55 de 2,3 km no campo de Roncador, na Bacia de Campos. Esses acidentes ocorrem em muitos casos devido à flambagem, à fadiga, ao intenso processo corrosivo, à deficiência de projeto, a falta de manutenção dos mesmos. Os dutos submarinos representam uma alta porcentagem no custo total de uma unidade de extração de petróleo. Sua importância se torna vital para a extração de petróleo em águas territoriais brasileiras. Seu projeto estrutural é uma atividade que pode ser dividida em várias etapas que devem ser seguidas à risca para garantir a conformidade do mesmo com as normas nacionais e internacionais sua operação em um dos mais inóspitos cenários da natureza: o fundo do mar. Nesta situação, o duto fica sujeito a carregamentos extremos, tanto estáticos quanto dinâmicos. Por este motivo, ele é constituído de uma estrutura de camadas bastante complexa. O objetivo deste trabalho é a avaliação das condições de operação de dutos submarinos sujeitos à ação de ondas e correntes marinhas. Um estudo paramétrico em elementos finitos é realizado para compreender a influência das condições de contorno, diferentes comprimentos de vão livre do duto submarino e variações do tipo de solo na resposta dinâmica do duto submarino. È desenvolvida também uma solução analítica para o problema. As soluções numérica, analítica e as recomendações da norma são confrontadas e tiradas conclusões sobre os principais fatores a serem considerados pelos engenheiros envolvidos na construção de estruturas deste tipo. / [en] Several accidents involving subsea pipelines have happened in Brazil as the pipeline PCR-01 / Lubnor and the P-55 platform 2.3 km in the Roncador field in the Campos Basin. These accidents occur in many cases due to buckling, fatigue, the intense corrosive process, design deficiency, lack of maintenance. The submarine pipelines represent a high percentage in the total cost of an oil extraction unit. Its importance becomes vital for the extraction of oil in Brazilian waters. Its structural design is an activity that can be divided into several steps that must be followed strictly to ensure its conformity with national and international standards operation in one of the harshest scenes of nature: the seabed. In this situation, the product is subject to extreme loads, both static and dynamic. For this reason, it consists of a complex layer structure. The objective of this study is the evaluation of the operating conditions of submarine pipelines subject to the action of waves and currents. A parametric study of finite element is performed to understand the influence of the boundary conditions, different span lengths of submarine pipeline and variations of soil type on the dynamic response of submarine pipeline. It is also developed an analytical solution to the problem. The numerical solutions, analytical and standard recommendations are compared and conclusions drawn on the main factors to be considered by the engineers involved in the construction of such structures.

[pt] O desenvolvimento consistente do método convencional dos elementos de contorno (CBEM), com a adição de conceitos da versão simplificada do método híbrido dos elementos de contorno (HBEM), proveniente do potencial variacional de Hellinger-Reissner, conduz-se a um processo computacionalmente mais econômico, sem a necessidade de ter sua precisão numérica reduzida para problemas de grande escala, podendo ser bidimensional ou tridimensional, de potencial ou elasticidade. Conseguiu-se mostrar que as matrizes de potencial duplo e simples do CBEM, H e G, respectivamente, cuja avaliação numérica requer a manipulação de integrais singulares e impróprias, podem ser obtidas de maneira expedita, eliminando-se quase toda a integração numérica, com exceção de algumas integrais regulares. Uma importante característica da formulação proposta, que advém da base variacional do HBEM, é a facilidade da obtenção de resultados em pontos internos, de maneira direta e sem a utilização de qualquer integral de contorno, já que a solução fundamental é a própria solução do problema. O presente trabalho pertence a um projeto cujo resultado final deve ser um código computacional para problemas de grande escala (milhões de graus de liberdade). Nesta fase, alguns exemplos numéricos foram testados para avaliar a aplicabilidade do método expedito, o seu esforço computacional e a convergência do resultado para as variáveis envolvidas no método. Para isso, foram implementados algoritmos para problemas bidimensionais de potencial e elasticidade - usando elementos lineares, quadráticos e cúbicos - e tridimensionais - usando elementos triangulares e quadrilaterais, lineares e quadráticos nos dois casos. Os códigos computacionais foram implementados focando na solução de problemas de grande escala. Espera-se que numa etapa final o projeto possa ser bem mais eficaz, com a incorporação de procedimentos do método fast multipole. / [en] The consistent development of the conventional boundary elements method (CBEM) by adding the concepts of the hybrid boundary element simplified method (HBEM) , from the Hellinger-Reissner variational potential leads to a computationally less intensive procedure, although not necessarily less accurate for large scale, two-dimensional or three-dimensional problems of potential and elasticity. It was shown that both single-layer and double-layer potential matrices, G and H, respectively, are obtained in an expeditious way that vanish almost any numerical integration, except for a few regular integrals, even G and H evaluation requires the handling of singular and improper integrals. The proposed formulation comes from the HBEM variational base and its evaluation at internal points is straightforward without the application of any boundary integral, since the fundamental solution is the analytical one. This work belongs to a project that aims a computer code for large-scale problems (millions of degrees of freedom). At this stage, some numerical examples were analyzed to evaluate the applicability of the method expeditious its computational effort and convergence of the results for the variables involved in the method. It was developed by the algorithms implementation for potential and elasticity problems. In the case of two-dimensional were employed linear, quadratic and cubic elements and to the three-dimensional case were employed triangular, quadrilateral, linear and quadratic elements in both cases. The computational codes were always implemented focused on solving largescale problems. It is expected that in a final stage of the project with the incorporation procedure of the method fast multipole, it can be more efficiently.

[pt] O processo de fraturamento de materiais quase-frágeis requer atenção especial para a predição da direção de propagação de fraturas. A simulação do fraturamento com o método dos elementos finitos (MEF) tem como desvantagem a dependência da trajetória da fratura com respeito à malha adotada. Além disso, há certa dificuldade para os modelos numéricos representarem a fratura em modo misto por conta dos parametros envolvidos. O Método dos Elementos Finitos Estendido (XFEM) é uma técnica que combina o MEF com funções de enriquecimento para representar descontinuidades no campo de deslocamentos. Neste contexto, discutem-se nesta dissertação os critérios para a nucleação e propagação de fraturas e sua implementação no contexto do XFEM. As implementações foram feitas no framework GeMA, um software desenvolvido no Tecgraf / PUC-Rio. Os critérios de propagação de fraturas implementados baseiamse na abordagem das tensões e permitem controlar diferentes geometrias e tamanhos da área de avaliação na ponta da trinca. Um estudo paramétrico é apresentado para modelar uma viga de concreto sob carregamento não proporcional com fratura em modo misto. Foram consideradas diferentes questões como: discretização da malha, zona de avaliação, iniciação e propagação de fraturas e técnicas de controle de solução. Além disso, outros modelos com diferentes condições de contorno foram analisados para validar os critérios em situações complexas. As constatações paramétricas obtidas através do estudo da viga se monstraram válidas para os demais modelos avaliados. As implementações dos critérios de propagação de fraturas no XFEM, demonstraram excelentes concordâncias nas simulações das trajetórias de fraturamento, comparado com os dados experimentais. / [en] The fracture process of quasi-brittle materials requires special attention for the prediction of the direction of fracture propagation. The fracture simulation with the finite element method (FEM) has as its disadvantage the dependence of the fracture trajectory with respect to the mesh adopted. Besides, there is some difficulty for numerical models to represent the fracture in mixed mode because of the parameters involved. The Extended Finite Element Method (XFEM) is a technique which combines the FEM with enrichment functions to represent discontinuities in the displacement field. In this context, this dissertation discusses the criteria for nucleation and propagation of fractures and their implementation in the context of XFEM. The implementations were made in the GeMA framework, a software developed at Tecgraf / PUC-Rio. The implemented crack growth criteria is based on the stress approach and allows to control different geometries and sizes of the evaluation area in the crack tip. A parametric study is presented for modeling a concrete beam under nonproportional loading with mixed-mode fracture. Different situations were taken into account such as mesh refinement, geometry and size of the evaluation region, crack initiation and propagation and solution control techniques. Also, several models with different loading and boundary conditions were made to validate the criteria under complex situations. The parametric findings obtained through the study of the beam proved to be valid for the other models. The implementations of the fracture propagation criteria in the XFEM demonstrated excellent agreement in the simulations of the fracture trajectories compared to the experimental data.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-04-11T12:25:14Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InfluenciaVariacaoRazao.pdf: 2017712 bytes, checksum: c97eb6c1192134f2ea8090bac4e210d0 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-08-29T17:23:12Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InfluenciaVariacaoRazao.pdf: 2017712 bytes, checksum: c97eb6c1192134f2ea8090bac4e210d0 (MD5) / Made available in DSpace on 2014-08-29T17:23:12Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_InfluenciaVariacaoRazao.pdf: 2017712 bytes, checksum: c97eb6c1192134f2ea8090bac4e210d0 (MD5) Previous issue date: 2008 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / Neste trabalho foi feito uma análise da combustão em um combustor ciclônico através de medidas experimentais da temperatura e concentração de gases na parede interna da câmara de combustão. Com o objetivo de encontrar parâmetros operacionais adequados para o projeto proposto, a alimentação de ar e serragem do combustor foi variada em razões de equivalência pobres (com excesso de ar) enquanto os dados experimentais eram computados. Os perfis encontrados foram confrontados com a teoria da combustão de sólidos e com os campos de temperatura e concentração de gases encontrados numericamente por Cunha (2005) através do software Fluent V.6.0. Nesta comparação foi possível encontrar boas concordâncias qualitativas entre as temperaturas medida e calculada, porém houve diferenças no quesito concentração de gases. Foi possível também identificar a razão de equivalência mínima para que o combustor ciclônico tivesse em seus produtos gasosos baixos teores de poluentes, como CO. Além disso, ao final deste trabalho foi proposta uma metodologia para o dimensionamento de combustores ciclônicos de acordo com a faixa de consumo de particulado que se deseja incinerar. / In this work it was made an analysis of the combustion in a cyclonic combustor through experimental measures of the temperature and gases concentration in the internal wall of the combustion chamber. With the objective of finding an appropriate operational parameters for the proposed design, the feeding of air and sawdust of combustor was varied in a poor equivalence ratio (with excess of air) while the experimental data were computed. The data profiles obtained were confronted with the theory of the combustion of solids and with the temperature fields and gases concentration founded by Cunha (2005) numerically in the code Fluent V.6.0. In this comparison it was possible to find good qualitative value among the temperatures measured and the one from calculations, however there were differences in the topic concentration of gases. It was possible also to identify the minimum equivalence rate that the cyclonic combustor had in their gaseous products low pollutant tenors, like CO. Besides, at the end of this work a methodology was proposed to find the appropriate size of a cyclonic combustor in agreement with the strip of sawdust consumption that it wants to incinerate.

[pt] Nesta dissertação estuda-se em detalhe o modelo de coluna submetido a carregamentos não conservativos dependentes do deslocamento. Analisa-se a influência da base elástica de Winkler e de apoios elásticos nos dois problemas clássicos de Beck e Leipholz. Apresenta-se uma metodologia para o método dinâmico de determinação da condição e carga crítica através de uma formulação matricial. A modelagem é desenvolvida para um elemento de viga esbelta através do método de Rayleigh-Ritz com o campo de deslocamentos descrito pela combinação de polinômios cúbicos e funções trigonométricas. Desenvolveu-se um programa em MATLAB de acordo com a metodologia apresentada. Os problemas clássicos foram utilizados para avaliar numericamente a eficácia desta abordagem. Apresentam-se as análises de influência de base elástica e apoios elásticos avaliados tanto separadamente como em conjunto. Isto permite mostrar que as cargas críticas e respectivos modos são afetados de forma distinta para cada tipo de apoio elástico e forma de carregamento. Finalmente, as análises são comparadas com resultados presentes na bibliografia, de modelos mais completos, considerando movimento de fluido ao longo do elemento, com o objetivo de avaliar a faixa de aplicabilidade do modelo de Beck para o caso de tubulações. / [en] In this dissertation, the model of column subjected to displacementdependent non-conservative loads is studied in detail. The influence of Winkler s elastic foundation and elastic supports is carried out for some classical problems such as Beck s and Leipholz s. A methodology based on a matrix formulation is presented for determining the critical condition and load in a dynamic approach. The model developed is based on Rayleigh-Ritz method, using a combination of cubic polynomials and trigonometric functions to describe the displacement field of slender beam elements. This methodology was implemented in a computer program developed in MATLAB. The classical problems were used as numerical tests to evaluate the accuracy of this approach. Later, the analysis of elastic foundation and elastic supports influence are presented separately and together. It was possible to show that both critical load and its respective modes are affected differently for each type of elastic support and load distribution. Finally, the influence analysis of Beck s problem is compared to results presented in in bibliography of a model which includes fluid movement within the element. This comparison has the purpose to estimate a range of applicability of the classical Beck s model to a pipe conveying fluid.

Neste trabalho formulações chamadas de "Sistemas de Parâmetros Concentrados': são usadas na obtenção de uma solução em forma fechada de um problema bidimensional de solidificação de ligas. Os sistemas concentrados são obtidos a partir do sistema diferencial original que descreve a distribuição de temperatura, resultando em modelos matemáticos mais simples que relacionam a temperatura no contorno do meio com uma nova temperatura média gerada por um processo de integração . A determinação de equações que relacionam a temperatura média e do contorno geram diferentes abordagens sendo a principal delas relacionada ao uso de fórmulas de integração numérica de Hermite, que propiciam a introdução de informações do contorno no modelo simplificado. Aqui o modelo bidimensional é abordado em meio simples e composto sendo que as equações unidimensionais simplificadas obtidas pela integração são tratadas pela Técnica das Transformadas Integrais Generalizadas (GITT). Resultados numéricos, obtidos com o software matemático Maple são apresentados. / In this work the so-called "lumped analysis'' is used to obta.in a closedform solution to a multidimensional solidification problem. The lumped systems provide simpler model, than the original one by using an integration scheme that results in average variables. Different approaches are basead on the choice of auxiliary equations to relate the average and the original variable after the integration process. Here the bidimensional problem is solved in homogeneous end in a composite medium and the resulting one-dimensional equation is solved by the Generalized Integral Tranform Technique. Numerical results are obta.ined by the symbolic software Maple.

O presente trabalho trata da formulação, algoritmização e implementação numérica de um problema não-linear de controle de fronteira livre sujeito a restrições também não lineares definido em [Fri 94] e relativo ao modelo de funcionamento de um conversor catalítico monolítico cerâmico. São apresentados resultados de algumas simulações numéricas, usando um programa em FORTRAN77, no ambiente de estações de trabalho SUN e DEC alfa 3000.