Orientador: Marco Lúcio Bittencourt / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica / Made available in DSpace on 2018-08-21T10:09:51Z (GMT). No. of bitstreams: 1 ValentedaCosta_GilbertoLuis_M.pdf: 6856772 bytes, checksum: bacab2fefeb1d730bb8f6760c97d0941 (MD5) Previous issue date: 2012 / Resumo: Este trabalho tem como objetivo principal a implementação de uma arquitetura de software para o Método de Elementos Finitos de Alta Ordem (MEF-AO), baseando-se no paradigma de programação orientada a objeto (POO) e no uso de técnicas de otimização de código fonte. O software foi escrito em linguagem C++ e desenvolvido sobre um framework com ferramentas que auxiliaram no desenvolvimento. A modelagem do sistema foi realizada de forma a facilitar e promover o reuso e manutenção do código. Buscou-se, também, a flexibilidade e generalização do MEF-AO ao permitir a variação nos procedimentos da construção das equações e o uso de malhas p não-uniforme. Neste caso, cada elemento pode ser interpolado com uma ordem polinomial diferente, além de permitir o uso de um algoritmo local de solução. Tal característica pode diminuir o número de operações e de armazenamento, pois o número de funções de forma é aumentado apenas onde é necessário o uso de mais pontos para interpolação da malha de solução. No final, o software é avaliado aplicando o problema de projeção para malha de quadrados e hexaedros / Abstract: The main objective of this work is the implementation of a software architecture for the High-Order Finite Element Method (HO-FEM), based on the Object Oriented Paradigm (OOP) and on source-code optimization techniques. The software was written in C++ programming language and developed over a framework which provided tools that assisted the implementation. The system was modeled so to promote code reuse and maintainability. Furthermore, the system modeling also provided flexibility and generalization for the HO-FEM by allowing modifications on the procedures for equation assembling and the use of p-non-uniform meshes. In this case, each element can be interpolated with different polynomial order, and allows the application of an algorithm for local solution. Such features can reduce the number of operations for memory allocation, since the number of shape functions is increased only where a higher density of points is needed by the solution mesh. Finally, the software is assessed by applying the projection problem for meshes of squares and hexahedros / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestre em Engenharia Mecânica

Orientador: Hugo Enrique Hernandez Figueroa / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Elétrica e de Computação / Made available in DSpace on 2018-08-20T11:52:35Z (GMT). No. of bitstreams: 1 Peternella_FellipeGrillo_M.pdf: 1154928 bytes, checksum: 3baba2c7d322af524dd1a6e07ee0e918 (MD5) Previous issue date: 2012 / Resumo: Desde as primeiras formulações em elementos finitos para o eletromagnetismo, os modos espúrios sempre foram um empecilho para a analise modal. Inicialmente, elementos nodais foram utilizados a fim de resolver a equação vetorial de onda no domínio da frequência. Naquele momento, os modos espúrios eram encontrados espalhados pelo espectro, o que dificultava a distinção entre soluções físicas e soluções espúrias. Anos mais tarde, o problema foi aparentemente resolvido pelo uso de elementos de aresta na equação vetorial de onda. Essa nova aproximação não eliminava por completo as soluções não físicas, apenas as confinava em torno do autovalor zero, o que não era problema na maioria das situações na analise modal. Recentemente, novas aplicações de elementos finitos no domínio do tempo resultaram em soluções instáveis, devido a presença dos modos espúrios em baixas frequências. Na década de 1990, Demkowicz[1], propôs uma formulação diferente em elementos finitos usando uma combinação da equação vetorial de onda com a condição do divergente de Gauss, que eliminou por completo os modos espúrios. O objetivo deste trabalho e analisar o impacto da condição de Gauss no FEM e entender como a mesma elimina os modos espúrios / Abstract: Since the earlier finite elements formulation for electromagnetics, spurious solutions have always plagued the modal analysis. Initially, nodal elements were employed to solve the vectorial wave equation in frequency domain. At that moment, the spurious modes could be found spread in frequency spectrum, and it was hard to distinguish between physical solutions and spurious ones. Some years later, this issue was apparently solved by using in the edge elements. This new approach does not completely eliminate non-physical solutions, but only confine them around to zero eigenvalue, which is not a problem in most of the situations in modal analysis. Recently, new applications of finite elements in time domain resulted to instable solution behavior, due to spurious modes at low frequencies. In the nineties, Demkowicz [1] proposed a different FEM formulation using a combination of vectorial wave equation and Gauss divergent equation, which definitely eliminates the spurious modes. The goal of this work is to analyze the impact of the Gauss condition in FEM and understand how it can eliminate spurious modes / Mestrado / Telecomunicações e Telemática / Mestre em Engenharia Elétrica

Orientador: Sérgio Tonini Button / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica / Made available in DSpace on 2018-08-20T23:55:31Z (GMT). No. of bitstreams: 1 Judice_GustavoHenrique_M.pdf: 13713684 bytes, checksum: ca6f8c3e7b2f080bdd41cef45bd2a552 (MD5) Previous issue date: 2012 / Resumo: O processo de conformação denominado Laminação Transversal com Cunha, industrialmente conhecido como Cross Wedge Rolling (CWR), tem se destacado na indústria pela sua flexibilidade, produtividade e economia de material, porém pode apresentar defeitos internos que, em alguns casos, inviabiliza a utilização dos produtos. É um defeito de difícil avaliação quanto a sua formação e propagação e ainda não se desenvolveu um mapa de processamento que represente de forma precisa e consistente as condições em que o defeito pode aparecer. Este estudo apresenta resultados dos campos de deformação plástica, tensão, temperatura e do aspecto do desenvolvimento de defeitos internos, acrescentados ao tarugo, através de simulações do processo de laminação transversal com cunha. O estudo se desenvolveu na criação de modelos de elementos finitos pelo software comercial Abaqus® versão 6.9-2, que seguiram condições geométricas de um processo real. Inicialmente, modelos bidimensionais foram simulados apresentando diferentes densidades de malha. Posteriormente, modelos com um vazio central foram simulados e confirmou-se a boa correlação de seus resultados com os encontrados na literatura. Na sequência, um modelo com propriedades de porosidade indicou que a região central do tarugo apresentou a maior densidade de vazios. Também foram simulados modelos com vazios e inclusões de diferentes tamanhos, disposições e formas, com destaque para a região em volta das inclusões, por apresentar a formação de vazios. Por último, foram simulados dois modelos tridimensionais, um apresentando um vazio interno no centro do tarugo, que apresentou os campos de tensão, de deformação plástica e de temperatura parecidos com o do modelo bidimensional correspondente, porém com valores maiores. O outro modelo tridimensional apresentou propriedades de porosidade e indicou que a maior probabilidade de formação de defeitos se concentrou na extensão longitudinal central do tarugo. Os maiores níveis de temperatura, deformação plástica e tensão foram presentes na região central do tarugo, tanto nos modelos bidimensionais quanto nos tridimensionais, indicando que essa região apresentou os maiores níveis de solicitação, o que pode favorecer o surgimento dos defeitos internos / Abstract: Cross Wedge Rolling is a metal forming process which has become very important in the industry for its flexibility, productivity and material savings, although presenting the critical problem of internal defects generation which restricts the use of rolled products, since it is very difficult to asset the defects formation and propagation and it was still not found a processing map to represent CWR with precise and consistent conditions with which the defect should or should not appear. In this study some numerical results are presented in terms of plastic strain, stress, temperature and the development of internal defects within the billet. These results were obtained with different finite element models developed with the commercial software Abaqus Version 6.9-2 by simulating materials with different internal characteristics, representing geometric conditions identified in the actual process. Initially two-dimensional models were simulated with different mesh densities; then, models with a central void were simulated and a good relationship was found with the results found in the literature. Then a model was simulated to characterize the porosity of the material, and was verified that the porosity is concentrated at the center of billet. Other models were also simulated with voids and inclusions of different sizes, distributions and shapes which were added manually. These models showed similar results, and the models with inclusions developed voids within the metallic matrix. Finally, two models (three-dimensional) were simulated, one with an internal void in the center of the billet, which showed the development of stress fields, plastic deformation and temperature similar to that corresponding two-dimensional model, but with higher values. The other three-dimensional model, which presented properties of porosity, showed that the highest probability of defect formation focused in the central and longitudinal extent region of the billet. All models showed that the central region of the billet presents the highest levels of temperature, plastic strain and stress probably favoring the initiation of internal defects / Mestrado / Materiais e Processos de Fabricação / Mestre em Engenharia Mecânica

Orientador : Marco Lucio Bittencourt / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-01T19:19:56Z (GMT). No. of bitstreams: 1 Pereira_CarlosEduardoLeite_M.pdf: 5572640 bytes, checksum: 45c5cb6b60e51ad992bd747fb7f167fe (MD5) Previous issue date: 2002 / Resumo: O presente trabalho tem por objetivo desenvolver a formulação variacional dos modelos de placa de Kirchhoff, incluindo o efeito de membrana, e de Reissner-Mindlin. Ambos os modelos são formulados utilizando-se a mesma metodologia. Primeiramente, definem-se as hipóteses cinemáticas dos modelos, considerando-se pequenos deslocamentos. Em seguida, deduzem-se os campos de deformação e de movimentos rígidos compatíveis com a cinemática. O conceito de trabalho interno é utilizado com o objetivo de encontrar os esforços internos compatíveis com a cinemática dos modelos. Em seguida, o Princípio dos Trabalhos Virtuais (PTV) é utilizado, a fim de caracterizar os esforços externos e os Problemas de Valor de Contorno (PVC) de equilíbrio. Conhecidos esses PVCs, a equação constitutiva de um material elástico, linear, isotrópico e homogêneo é utilizada para obter as relações entre os campos de tensão e de deformação dos modelos de placa, assim como os PVCs ou forma forte em termos da cinemática. As formas fracas das equações diferenciais dos modelos são então aproximadas pelo Método dos Elementos Finitos. Finalmente, considerações sobre os elementos de placa são apresentadas / Abstract: The present work aims to present the variational formulation of the Kirchhoff plate model with membrane effect and the Reissner-Mindlin plate model. Both models are formulated using the same methodology. First, it is defined the kinematic hypotheses of the models for small displacements. After that, the strain and rigid displacement fields compatible with the kinematic hypotheses are obtained. The internal work concept is used to find the internal forces compatible with the kinematic hypotheses of the models. The Virtual Work Principle is used to find the external forces and the equilibrium Boundary Value Problems (BVP). The constitutive equation of the elastic, linear, isotropic and homogeneous material is used to obtain the relations between the stress and strain fields. Based on that, BVPs in terms of the kinematics are obtained. The weak forms of the differential equation are approximated by Finite Element Method. Finally, comments about plate elements are presented. / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestre em Engenharia Mecânica

Orientador : Marco Lucio Bittencourt / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-01T16:50:37Z (GMT). No. of bitstreams: 1 Ferreira_WallaceGusmao_M.pdf: 8465867 bytes, checksum: 43bfa97c360052855e2188e207de0c0d (MD5) Previous issue date: 2002 / Resumo: O presente trabalho tem por objetivo desenvolver a análise estrutural em problemas de fadiga mecânica, usando o Método dos Elementos Finitos (MEF) e recursos de Engenharia de Software. Inicialmente é feita uma revisão dos conceitos básicos de fadiga mecânica, descrevendo as frentes atuais de pesquisa e as metodologias convencionais de projeto. O modelo constitutivo de material para o tratamento de fadiga é definido com base na Mecânica do Dano e a solução numérica do problema de equilíbrio é feita através do MEF. Os resultados obtidos com a implementação do modelo de dano e da análise estrutural de componentes submetidos à fadiga são compatíveis com aqueles apresenta dos na literatura. Os programas computacionais são implementados usando conceitos de Engenharia de Software. Como exemplo da aplicação desses conceitos, é desenvolvido um ambiente gráfico para geração de malhas de elementos finitos, utilizando componentes de software disponíveis comercialmente / Abstract: This work deals with the structural analisys of problems involving mechanical fatigue using the Finite Element Method and Software Engineering concepts. Initially, the basic concepts of mechanical fatigue, the conventional design techniques and the major research areas are reviewed. The material constitutive model to treat the fatigue process is developed based on Damage Mechanics and the numerical solution for the equilibrium equation is obtained by the Finite Element Method. The results obtained with the implementation of the constitutive mo deI and for structural analisys have good agreement with the literature references. All the computational programs are implemented using the concepts of Software Engineering. As example of the application of these concepts, a graphical environment for finite element mesh generation is developed based on third party software components / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestre em Engenharia Mecânica

Orientador : Loir Afonso Moreira / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-01T20:15:40Z (GMT). No. of bitstreams: 1 Baptistella_PedroHenrique_M.pdf: 4895786 bytes, checksum: 524a48114b7bde3e42723a076929150c (MD5) Previous issue date: 2001 / Resumo: Em 1970, Ahmad et aI. apresentaram uma formulação, para a análise de placas e cascas, baseada na degeneração de um elemento sólido tridimensional, através da redução de sua dimensão na direção da espessura. O campo de deslocamento era, então, interpolado a partir das funções de forma da família Serendipity e, para a obtenção das características do elemento, utilizava-se integração numérica consistente. Os resultados obtidos por esse elemento foram excelentes para situações de placas e cascas moderadamente grossas, atendendo, portando, à teoria de Mindlin- Reissner ~ contudo, com a redução da espessura, o elemento tomava-se excessivamente rígido e os resultados não tendiam, como era de se esperar, àqueles da teoria clássica de Kirchho./f para placas e cascas finas. O objetivo deste trabalho é o estudo comparativo, na análise de estabilidade de placas e cascas, entre o desempenho, utilizando-se integração numérica consistente, do elemento finito isoparamétrico, quadrilateral, quadrático da família Serendipity e do elemento finito isoparamétrico, quadrilateral, quadrático da família Lagrangiana, obtidos, ambos, a partir da formulação de Ahmad. Pretende-se estabelecer, dentre os dois elementos, qual é o mais adequado para iniciar o desenvolvimento de um elemento finito subparamétrico do tipo hierárquico, baseado no conceito da aproximação p, visando retirar a característica de rigidez excessiva do elemento de Ahmad, para situações de placas e cascas finas / Abstract: In 1970, Ahmad et aI. presented a formulation for the analysis of plates and shells, based on the degeneration of a three-dimensional solid element, by reducing its dimention in the thickness direction. The displacement field was then interpolated through the shape function of Serendipity fami1y and, in order to achieve the element characteristics, consistent integration was used. The obtained results were excellent for moderately thick plates and shells aplications, being in concordance with the Mind/in-Reissner's theory; however, with the reduction of the thickness, the element became excessively rigid and the results did not tend to those ofthe Kirchhoff's theory for thin plates and shells. The purpose of this work is the comparative study inthestabilityanalysis of plates and shells concerning the performance, considering consistent numerical integration, of the isoparametric, quadrilateral, quadratic finite element of Serendipity fami1y and of the isoparametric, quadrilateral, quadratic finite element of Lagrangian fami1y, achieved, both, ftom the Ahmad' s formulation, one intends to estabilish, between both elements, which is the most adequa te to start the development of a subparametric hierarchical element, based on p aproximation concept, objectifying to remove the characteristic of excessive rigidy of the Ahmad' s elements for the situation of thin plates and shells / Mestrado / Mestre em Educação Física

Orientador : Hugo Enrique Hernandez Figueroa / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Eletrica e de Computação / Made available in DSpace on 2018-08-01T22:09:13Z (GMT). No. of bitstreams: 1 RubioMercedes_CosmeEustaquio_D.pdf: 5866050 bytes, checksum: 68669f333d337e8c63e289efcad5b166 (MD5) Previous issue date: 2002 / Resumo: Foram desenvolvidos códigos numéricos, baseados em elementos finitos 2D no domínio da freqüência, para a simulação de junçalimentadas por guias fotônicos, supostamente muito longos. Sobre as fronteiras da janela computacional, nas regiões de radiação não guiada, foram utilizadas PMLs (Perfectly Matched Layers). Sobre as portas ou regiões de acesso, cinco esquemas foram desenvolvidos, utilizando-se cinco condições de contorno diferentes: a primeira baseada na expansão modal dos campos; a segunda considerando-se uma variação de onda TEM (transversal eletromagnética); a terceira baseada na aproximação paraxial dos campos, e as duas ultimas baseadas em aproximações de Padé (Padé(l,l) e Padé(2,2)) ou de ângulo largo. Estas três últimas, até onde sabemos, estão sendo propostas pela primeira vez, e representam alternativas eficientes (relativamente à primeira) e confiáveis (relativamente à segunda), para a simulação de uma extensa gama de junções planares. O desempenho dos esquemas propostos, foi avaliado através da análise de um número considerável de ef) exemplos. Uma estratégia para analisar circuitos fotônicos, combinando os esquemas apresentados e lC) o BPM (Beam Propagation Method), é discutida e fartamente ilustrada. Estruturas de interesse atual, como ressoadores altamente compactos, utilizados como elementos add/drop em sistemas DWTIM) (Dense WavelengthDivision Multiplexing), foram analisados em detalhe / Abstract: Novel 2D finite element formulations in the frequency domain, for the simulation of arbitrary planar photonic junctions, were developed and thoroughly described. Such junctions are composed by multiple access ports, fed by supposedly very long photonics guides. On the boundaries of the computational window, placed over not guided radiation regions, PMLs (Perfectly Matched Layers) wereadopted. On the ports or access wave guides,five different boundary conditions were implemented: The first, basedon a modal expansionof the fields ; the second, considering apure modal wave variation the third, based on the paraxial approach of the fields,and the last two, based on the Padé (Padé(l,l) and Padé (2,2)) or wide angle approximations. These three last boundary conditions, to the best of our knowledge,were proposed here for the first time, and represent efficient alternatives (relatively to the first) and more reliable (relatively to the second) for the simulation of a wide range of planar junctions. The performanceof the proposed schemes,was evaluated through the analysis of a considerable number of examples. A strategy to analyze photonic circuits, combining the presented schemes with the BPM (Beam Propagation Method), is also discussed and widely illustrated. Structures of current interest , such as highly compact resonators, used as add/drop devices in DWDM (Dense Wavelength Division Multiplexing) systems were analyzed in detail / Doutorado / Telecomunicações e Telemática / Doutor em Engenharia Elétrica

Orientador : Loir Afonso Moreira / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-01T21:31:25Z (GMT). No. of bitstreams: 1 Paschoalini_AmarildoTabone_D.pdf: 8353999 bytes, checksum: bfeebd1be75d31e7d86303a81f912b57 (MD5) Previous issue date: 2001 / Resumo: Este trabalho apresenta um elemento subparamétrico hierárquico baseado na versão p do método dos elementos finitos para a análise de placas e cascas. O primeiro nível de aproximação da solução é obtido através do elemento finito quadrilateral quadrático de nove nós da família Lagrangeana, baseado na degeneração de um elemento sólido tridimensional e na formulação de Reissner-Mindlin, com integração numérica consistente. Para outros níveis de aproximação, sucessivos refinamentos hierárquicos são usados, objetivando remover a caracteristica de rigidez excessiva do elemento na análise de placas e cascas finas. Tal formulação é aplicada na análise estática e dinâmica de placas e cascas (finas e moderadamente grossas). Exemplos numéricos são apresentados para mostrar a precisão, eficiência e vantagens da presente formulação. Os resultados numéricos obtidos nos exemplos de aplicação são comparados com soluções analíticas e outras técnicas numéricas, disponíveis na literatura / Abstract: This research presents a subparametric hierarchical finite element based on the p-version concept for the analysis of plates and shells. The first level of approximation for the solution is obtained through the isoparametric quadrilateral quadratic nine-no de Lagrangean shell finite element, based on the degeneration of three-dimensional solid element and the Reissner-Mindlin's formulation, with consistent numerical integration. For other approximation levels, successive hierarchical refinements are used, aiming to remove the characteristic of excessive rigidity of the isoparametric element in the analysis ofthin plates and shells. Such formulation is applied to static and dynamic analysis of plates and shells (thin and moderately thick). Numerical examples are presented to show the accuracy, efficiency and advantages of present formulation. Numerical results obtained for the application examples are compared with analytical solutions and other numerical techniques, available in the open literature / Doutorado / Mecanica dos Sólidos e Projeto Mecanico / Doutor em Engenharia Mecânica

Orientador: Philippe Remy Bernard Devloo / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Civil / Made available in DSpace on 2018-08-02T02:48:48Z (GMT). No. of bitstreams: 1 Longhin_GustavoCamargo_M.pdf: 6501443 bytes, checksum: c3b310f91b0c026c56f80bf6bc4c6fba (MD5) Previous issue date: 2001 / Resumo: Para obtenção da solução do sistema de equações, geralmente elabora-se um código que armazena a matriz dos coeficientes na memória RAM. Em seguida inicia-se o processo de decomposição desta matriz. A matriz é formada pela contribuição de cada elemento do domínio aos graus de liberdade do problema. Quanto maior o problema, maior o tamanho desta matriz e consequentemente mais memória RAM será necessária para seu armazenamento. Portanto, um procedimento de montagem anterior à decomposição pode ser inviabilizado devido ao tamanho da matriz. Com esta motivação, Bruce Irons desenvolve no início da década de 70, um método de resolução que não exige uma montagem inicial da matriz de rigidez global de equações. Neste método é definida uma estrutura onde uma equação totalmente adicionada é imediatamente decomposta e armazenada num dispositivo independente. A matriz que recebe as equações totalmente adicionadas é denominada matriz frontal e com isso o método também é denominado método frontal. Sobre a estrutura frontal são aplicadas técnicas de optimização por paralelismo. São utilizados equipamentos com memória compartilhada e portanto, utilizão-se as bibliotecas oriundas da especificação posix (pthread no ambiente GNU & Linux) para desenvolvimento multi-threading. São apresentados resultados comparando o método frontal com outros métodos bem como as comparações entre os métodos seriais e os paralelos. Aplica-se um tratamento orientado a objetos para desenvolvimento dos solvers. Nota-se com a orientação a objetos um excelente grau de modularidade, documentação, extendibilidade e manutenção no código elaborado. A utilização de UML (Unified Modeling Language) é também de grande valia no desenvolvimento/planejamento do projeto / Abstract: Obtaining the solution of a system of linear equations, generally results in a code elaboration which stores the matrix coeficients in the RAM memory and afterwards, some decomposition processes starts. The matrix is assembled summing up the contributions of each element from the domain to the problem's degrees of freedom. The bigger the problem, the larger the assembled matrix, therefore a higher requirement regarding the RAM memory capacity. From this, a procedure which does not assemble the matrix of coeficients prior to its decomposition would be more interesting. With that motivation, Bruce Irons developed in the beggining of the seventies a procedure which does not require an initial assembly of the global stiffness matrix. In this method a structure is defined where a totally added equation is immediatelly decomposed and the decomposition results are stored in an independent storage device. The matrix which receives the equations contribution was called frontal matrix and so was the method. On that frontal structure parallel optimization techniques are applied. Shared memory equipments are the hardware basis for the implementation and accordingly, public domain multithreading libraries based on the posix specification are used (pthread under GNU & Linux) for the multi-threading development. Results are shown comparing standart methods against the frontal solver as well as serial codes against parallel ones. Object oriented techniques are applied for the solvers development and planning. As a result, excelent degrees of modularity, extendibiIity, documentation and management are observed. The Unified Modelling Language (UML) utilization as a helpping tool for object oriented development was also very important / Mestrado / Estruturas / Mestre em Engenharia Civil

Orientador : Hugo Enrique Hernandez Figueroa / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Eletrica e de Computação / Made available in DSpace on 2018-08-03T14:16:13Z (GMT). No. of bitstreams: 1 Silva_JosePatrocinioda_D.pdf: 1389226 bytes, checksum: 1aa7239a889f3d8f840db6896db2c349 (MD5) Previous issue date: 2003 / Doutorado