A pesquisa foi desenvolvida com a finalidade de se analisar o comportamento elástico e inelástico de estruturas de cabos, sob carregamento estático, pela técnica dos elementos finitos. Para a discretização das estruturas de cabos adotou-se um elemento retilíneo isoparamétrico, e na formulação dos princípios variacionais foram considerados a não linearidade física do material e a não linearidade geométrica, inerente às estruturas de cabos. Empregou-se o método incrementai-iterativo de Newton-Raphson para a resolução do sistema de equações não lineares. Alternativamente a solução foi obtida da minimização direta da função energia potencial total da estrutura de cabos, utilizando-se o método do Gradiente Conjugado. Vários exemplos de estruturas de cabos comumente empregadas em coberturas de grandes vãos foram estudados. Os resultados obtidos segundo os dois procedimentos foram analisados e a eficiência elos métodos foi comprovada. / The research was developed with the purpose of analyzing the elaslic and the inelaslic behaviour of cable slructures, under static load. The finite element method with straight isoparametric element was used. In the formulation of the variational principles, the material and the geometric non-linearity were considered. The Newton-Raphson method was used to resolve the non-linear equations system. Alternatively the solution was obtained of the direct minimization of the total potential energy of the cable structure; in this case using the Conjugated Gradient Method. Several examples of cable structures, commonly employed in roofs of large spaces were studied. The results were obtained according to both processes and the efficiency of the methods was verified.

A pesquisa foi desenvolvida com a finalidade de se analisar o comportamento elástico e inelástico de estruturas de cabos, sob carregamento estático, pela técnica dos elementos finitos. Para a discretização das estruturas de cabos adotou-se um elemento retilíneo isoparamétrico, e na formulação dos princípios variacionais foram considerados a não linearidade física do material e a não linearidade geométrica, inerente às estruturas de cabos. Empregou-se o método incrementai-iterativo de Newton-Raphson para a resolução do sistema de equações não lineares. Alternativamente a solução foi obtida da minimização direta da função energia potencial total da estrutura de cabos, utilizando-se o método do Gradiente Conjugado. Vários exemplos de estruturas de cabos comumente empregadas em coberturas de grandes vãos foram estudados. Os resultados obtidos segundo os dois procedimentos foram analisados e a eficiência elos métodos foi comprovada. / The research was developed with the purpose of analyzing the elaslic and the inelaslic behaviour of cable slructures, under static load. The finite element method with straight isoparametric element was used. In the formulation of the variational principles, the material and the geometric non-linearity were considered. The Newton-Raphson method was used to resolve the non-linear equations system. Alternatively the solution was obtained of the direct minimization of the total potential energy of the cable structure; in this case using the Conjugated Gradient Method. Several examples of cable structures, commonly employed in roofs of large spaces were studied. The results were obtained according to both processes and the efficiency of the methods was verified.

The purpose of this thesis is to develop corotational beam elements for the nonlinear dynamic analyse of flexible beam structures. Whereas corotational beam elements in statics are well documented, the derivation of a corotational dynamic formulation is still an issue. In the first journal paper, an efficient dynamic corotational beam formulation is proposed for 2D analysis. The idea is to adopt the same corotational kinematic description in static and dynamic parts. The main novelty is to use cubic interpolations to derive both inertia terms and internal terms in order to capture correctly all inertia effects. This new formulation is compared with two classic formulations using constant Timoshenko and constant lumped mass matrices. In the second journal paper, several choices of parametrization and several time stepping methods are compared. To do so, four dynamic formulations are investigated. The corotational method is used to develop expressions of the internal terms, while the dynamic terms are formulated into a total Lagrangian context. Theoretical derivations as well as practical implementations are given in detail. Their numerical accuracy and computational efficiency are then compared. Moreover, four predictors and various possibilities to simplify the tangent inertia matrix are tested. In the third journal paper, a new consistent beam formulation is developed for 3D analysis. The novelty of the formulation lies in the use of the corotational framework to derive not only the internal force vector and the tangent stiffness matrix but also the inertia force vector and the tangent dynamic matrix. Cubic interpolations are adopted to formulate both inertia and internal local terms. In the derivation of the dynamic terms, an approximation for the local rotations is introduced and a concise expression for the global inertia force vector is obtained. Four numerical examples are considered to assess the performance of the new formulation against two other ones based on linear interpolations. Finally, in the fourth journal paper, the previous 3D corotational beam element is extended for the nonlinear dynamics of structures with thin-walled cross-section by introducing the warping deformations and the eccentricity of the shear center. This leads to additional terms in the expressions of the inertia force vector and the tangent dynamic matrix. The element has seven degrees of freedom at each node and cubic shape functions are used to interpolate local transversal displacements and axial rotations. The performance of the formulation is assessed through five examples and comparisons with Abaqus 3D-solid analyses. / <p>QC 20131017</p>

L’objectif de cette thèse est de proposer des éléments finis poutres corotationnels 2D et 3D pour l’analyse du comportement dynamique non-linéaire des structures à barres. La contribution majeure de cette thèse est l’utilisation de fonctions d’interpolations cubiques à la fois pour la détermination de l’expression des efforts internes mais aussi celle des termes d’inertie. En négligeant le carré du déplacement transversal dans le repère local, une expression analytique des termes dynamiques en 2D est obtenue. Sur base d’une étude comparative approfondie sur la paramétrisation des rotations et les algorithmes d’intégration temporelle et d’une approximation des rotations locales, nous proposons deux éléments finis poutre 3D précis et robustes. Contrairement au premier, le second élément 3D prend en compte les déformations de gauchissement et l'excentricité du centre de cisaillement. Les diverses comparaisons réalisées démontrent la pertinence des formulations proposées. / The purpose of this thesis is to propose several corotational beam formulations for both 2D and 3D nonlinear dynamic analyse of flexible structures. The main novelty of these formulations is that the cubic interpolation functions are used to derive not only the internal force vector and the tangent stiffness matrix but also the inertial force vector and the dynamic matrix. By neglecting the quadratic terms of the local transversal displacements, closed-form expressions for the inertial terms are obtained for 2D problems. Based on an extensive comparative study of the parameterizations of the finite rotations and the time stepping method, and by adopting an approximation of the local rotations, two consistent and effective beam formulations for 3D dynamics are developed. In contrast with the first formulation, the second one takes into account the warping deformations and the shear center eccentricity. The accuracy of these formulations is demonstrated through several numerical examples.

Flexible beam structures are popular in civil and mechanical engineering. Many of these structures undergo large displacements and finite rotations, but with small deformations. Their dynamic behaviors are usually investigated using finite beam elements. A well known method to derive such beam elements is the corotational approach. This method has been extensively used in nonlinear static analysis. However, its application in nonlinear dynamics is rather limited. The purpose of this thesis is to investigate the nonlinear dynamic behavior of flexible beam structures using the corotational method. For the 2D case, a new dynamic corotational beam formulation is presented. The idea is to adopt the same corotational kinetic description in static and dynamic parts. The main novelty is to use cubic interpolations to derive both inertia terms and internal terms in order to capture correctly all inertia effects. This new formulation is compared with two classic formulations using constant Timoshenko and constant lumped mass matrices. This work is presented in the first appended journal paper. For the 3D case, update procedures of finite rotations, which are central issues in development of nonlinear beam elements in dynamic analysis, are discussed. Three classic and one new formulations of beam elements based on the three different parameterizations of the finite rotations are presented. In these formulations, the corotational method is used to develop expressions of the internal forces and the tangent stiffness matrices, while the dynamic terms are formulated into a total Lagrangian context. Many aspects of the four formulations are investigated. First, theoretical derivations as well as practical implementations are given in details. The similarities and differences between the formulations are pointed out. Second, numerical accuracy and computational efficiency of these four formulations are compared. Regarding efficiency, the choice of the predictor at each time step and the possibility to simplify the tangent inertia matrix are carefully investigated. This work is presented in the second appended journal paper. To make this thesis self-contained, two chapters concerning the parametrization of the finite rotations and the derivation of the 3D corotational beam element in statics are added. / QC 20120521

O objetivo precípuo deste estudo é a implementação computacional de modelos constitutivos elásticos e elastoplásticos para materiais com gradação funcional em regime de grandes deslocamentos e elevadas deformações. Para simular numericamente um problema estrutural, são empregados aqui elementos finitos sólidos (tetraédrico e hexaédrico) com ordem de aproximação polinomial qualquer. Grandezas da Mecânica Não Linear do Contínuo, como deformação e tensão, são utilizadas na formulação deste estudo. Para reproduzir os grandes deslocamentos, é empregada a análise não linear geométrica. A descrição adotada aqui é a Lagrangiana total, e o equilíbrio da estrutura é expresso pelo Princípio da Mínima Energia Potencial Total. Com relação à resposta elástica do material, são usadas leis constitutivas hiperelásticas, nas quais a relação tensão-deformação é obtida a partir de um potencial escalar. O comportamento elastoplástico do material é definido pela decomposição da deformação nas parcelas elástica e plástica, pelo critério de plastificação de von-Mises, pela lei de fluxo associativa, pelas condições de consistência e de complementaridade, pelo parâmetro de encruamento isotrópico e pelo tensor das tensões inversas, relacionado ao encruamento cinemático. Duas formulações elastoplásticas são utilizadas aqui: a de Green-Naghdi, na qual a deformação é decomposta de forma aditiva; e a hiperelastoplástica, em que o gradiente é decomposto de forma multiplicativa. É empregado também o conceito de material com gradação funcional (GF), a qual é definida como a variação gradual (contínua e suave) das propriedades constitutivas do material. A solução numérica do equilíbrio de forças é feita via método iterativo de Newton-Raphson. Para satisfazer o critério de plastificação, são utilizadas as estratégias de previsão elástica, e de correção plástica via algoritmos de retorno. Basicamente foram desenvolvidos cinco programas computacionais: o gerador automático das funções de forma; o gerador de malhas de elementos finitos sólidos; o código para análise de materiais em regime elástico; o código para análise de materiais em regime elastoplástico; e o programa de pós-processamento. Além desses, o aluno teve contato com os programas EPIM3D e DD3IMP ao longo de seu estágio de doutorado na Universidade de Coimbra (Portugal). Os programas EPIM3D e DD3IMP são empregados para analisar, respectivamente, materiais em regime elastoplástico, e processos de conformação de metais. Para o problema da barra sob tração uniaxial uniforme, são descritas equações e soluções analíticas para materiais homogêneos e com GF em regime elastoplástico. Para reduzir o tempo de simulação, foi empregada a programação em paralelo. De acordo com os resultados das simulações numéricas, as principais conclusões são: o refinamento da malha de elementos finitos melhora a precisão dos resultados para materiais em regimes elástico e elastoplástico; as formulações elastoplásticas de Green-Naghdi e hiperelastoplástica parecem ser equivalentes para pequenas deformações; a formulação hiperelastoplástica é equivalente ao modelo mecânico dos programas EPIM3D e DD3IMP para materiais em regime de pequenas deformações elásticas; foram constatados ganhos significativos, em termos de tempo de simulação, com a paralelização dos códigos computacionais de análise estrutural; e os programas desenvolvidos são capazes de simular - com precisão - problemas complexos, como a membrana de Cook e o cilindro fino transversalmente tracionado. / The main objective of this study is the computational implementation of elastic and elastoplastic constitutive models for functionally graded materials in large deformation regime. In order to numerically simulate a structural problem, the finite elements used are solids (tetrahedric and hexahedric) of any order of approximation. Entities from Nonlinear Continnum Mechanics, as strain and stress, are used in the present formulation. To reproduce the finite displacements, the geometrically nonlinear analysis is employed. The description adopted here is the total Lagrangian, and the structural equilibrium is expressed by means of the Principal of Minimum Total Potential Energy. Regarding the elastic material response, hyperelastic constitutive laws are used, in which the stress-strain relation is obtained from a scalar potential. The elastoplastic material behavior is defined by the strain decomposition in the elastic and plastic parts, by the von-Mises yield criterion, by the associative flow law, by the consistency and complementarity conditions, by the isotropic hardening parameter, and by the backstress tensor, related to the kinematic hardening. Two elastoplastic formulations are used here: the Green-Naghdi one, in which the strain is additively decomposed; and the hyperelastoplasticiy, in which the gradient is multiplicatively decomposed. The concept of functionally graded (FG) material, in which the constitutive properties vary gradually (continuous and smoothly), is also used. The numerical solution of the forces equilibrium is obtained via Newton-Raphson iterative procedure. In order to satisfy the yield criterion, the strategies of elastic prediction and plastic correction (via return algorithms) are used. Basically, five computer codes have been developed: the automatic shape functions generator; the solid mesh generator; the code for analysis of materials in the elastic regime; the code for analysis of materials in the elastoplastic regime; and the post-processor. Besides these, the student had contact with the programs EPIM3D and DD3IMP during his doctoral stage in the University of Coimbra (Portugal). The programs EPIM3D and DD3IMP are employed to analyze, respectively, materials in the elastoplastic regime, and sheet-metal forming processes. For the problem of the bar under uniform uniaxial tension, equations and analytical solutions are described for homogeneous and FG materials. To reduce the simulation time, the parallel programming has been employed. According to the numerical simulation results, the main conclusions are: the results accuracy is improved with mesh refinement for materials in the elastic and elastoplastic regimes; the Green-Naghdi elastoplastic formulation and the hyperelastoplasticity appear to be equivalent for small strains; the hyperelastoplastic formulation is equivalent to the mechanical model of the programs EPIM3D and DD3IMP for materials the small elastic strains regime; simulation time reduction has been obtained with the parallelization of the computer codes for structural analysis; the developed programs are capable of simulating, precisely, complex problems, such as the Cook\'s membrane and the pulled thin cylinder.

The purpose of this thesis is to propose several corotational beam formulations for both 2D and 3D nonlinear dynamic analyse of flexible structures. The main novelty of these formulations is that the cubic interpolation functions are used to derive not only the internal force vector and the tangent stiffness matrix but also the inertial force vector and the dynamic matrix. By neglecting the quadratic terms of the local transversal displacements, closed-form expressions for the inertial terms are obtained for 2D problems. Based on an extensive comparative study of the parameterizations of the finite rotations and the time stepping method, and by adopting an approximation of the local rotations, two consistent and effective beam formulations for 3D dynamics are developed. In contrast with the first formulation, the second one takes into account the warping deformations and the shear center eccentricity. The accuracy of these formulations is demonstrated through several numerical examples.

O objetivo precípuo deste estudo é a implementação computacional de modelos constitutivos elásticos e elastoplásticos para materiais com gradação funcional em regime de grandes deslocamentos e elevadas deformações. Para simular numericamente um problema estrutural, são empregados aqui elementos finitos sólidos (tetraédrico e hexaédrico) com ordem de aproximação polinomial qualquer. Grandezas da Mecânica Não Linear do Contínuo, como deformação e tensão, são utilizadas na formulação deste estudo. Para reproduzir os grandes deslocamentos, é empregada a análise não linear geométrica. A descrição adotada aqui é a Lagrangiana total, e o equilíbrio da estrutura é expresso pelo Princípio da Mínima Energia Potencial Total. Com relação à resposta elástica do material, são usadas leis constitutivas hiperelásticas, nas quais a relação tensão-deformação é obtida a partir de um potencial escalar. O comportamento elastoplástico do material é definido pela decomposição da deformação nas parcelas elástica e plástica, pelo critério de plastificação de von-Mises, pela lei de fluxo associativa, pelas condições de consistência e de complementaridade, pelo parâmetro de encruamento isotrópico e pelo tensor das tensões inversas, relacionado ao encruamento cinemático. Duas formulações elastoplásticas são utilizadas aqui: a de Green-Naghdi, na qual a deformação é decomposta de forma aditiva; e a hiperelastoplástica, em que o gradiente é decomposto de forma multiplicativa. É empregado também o conceito de material com gradação funcional (GF), a qual é definida como a variação gradual (contínua e suave) das propriedades constitutivas do material. A solução numérica do equilíbrio de forças é feita via método iterativo de Newton-Raphson. Para satisfazer o critério de plastificação, são utilizadas as estratégias de previsão elástica, e de correção plástica via algoritmos de retorno. Basicamente foram desenvolvidos cinco programas computacionais: o gerador automático das funções de forma; o gerador de malhas de elementos finitos sólidos; o código para análise de materiais em regime elástico; o código para análise de materiais em regime elastoplástico; e o programa de pós-processamento. Além desses, o aluno teve contato com os programas EPIM3D e DD3IMP ao longo de seu estágio de doutorado na Universidade de Coimbra (Portugal). Os programas EPIM3D e DD3IMP são empregados para analisar, respectivamente, materiais em regime elastoplástico, e processos de conformação de metais. Para o problema da barra sob tração uniaxial uniforme, são descritas equações e soluções analíticas para materiais homogêneos e com GF em regime elastoplástico. Para reduzir o tempo de simulação, foi empregada a programação em paralelo. De acordo com os resultados das simulações numéricas, as principais conclusões são: o refinamento da malha de elementos finitos melhora a precisão dos resultados para materiais em regimes elástico e elastoplástico; as formulações elastoplásticas de Green-Naghdi e hiperelastoplástica parecem ser equivalentes para pequenas deformações; a formulação hiperelastoplástica é equivalente ao modelo mecânico dos programas EPIM3D e DD3IMP para materiais em regime de pequenas deformações elásticas; foram constatados ganhos significativos, em termos de tempo de simulação, com a paralelização dos códigos computacionais de análise estrutural; e os programas desenvolvidos são capazes de simular - com precisão - problemas complexos, como a membrana de Cook e o cilindro fino transversalmente tracionado. / The main objective of this study is the computational implementation of elastic and elastoplastic constitutive models for functionally graded materials in large deformation regime. In order to numerically simulate a structural problem, the finite elements used are solids (tetrahedric and hexahedric) of any order of approximation. Entities from Nonlinear Continnum Mechanics, as strain and stress, are used in the present formulation. To reproduce the finite displacements, the geometrically nonlinear analysis is employed. The description adopted here is the total Lagrangian, and the structural equilibrium is expressed by means of the Principal of Minimum Total Potential Energy. Regarding the elastic material response, hyperelastic constitutive laws are used, in which the stress-strain relation is obtained from a scalar potential. The elastoplastic material behavior is defined by the strain decomposition in the elastic and plastic parts, by the von-Mises yield criterion, by the associative flow law, by the consistency and complementarity conditions, by the isotropic hardening parameter, and by the backstress tensor, related to the kinematic hardening. Two elastoplastic formulations are used here: the Green-Naghdi one, in which the strain is additively decomposed; and the hyperelastoplasticiy, in which the gradient is multiplicatively decomposed. The concept of functionally graded (FG) material, in which the constitutive properties vary gradually (continuous and smoothly), is also used. The numerical solution of the forces equilibrium is obtained via Newton-Raphson iterative procedure. In order to satisfy the yield criterion, the strategies of elastic prediction and plastic correction (via return algorithms) are used. Basically, five computer codes have been developed: the automatic shape functions generator; the solid mesh generator; the code for analysis of materials in the elastic regime; the code for analysis of materials in the elastoplastic regime; and the post-processor. Besides these, the student had contact with the programs EPIM3D and DD3IMP during his doctoral stage in the University of Coimbra (Portugal). The programs EPIM3D and DD3IMP are employed to analyze, respectively, materials in the elastoplastic regime, and sheet-metal forming processes. For the problem of the bar under uniform uniaxial tension, equations and analytical solutions are described for homogeneous and FG materials. To reduce the simulation time, the parallel programming has been employed. According to the numerical simulation results, the main conclusions are: the results accuracy is improved with mesh refinement for materials in the elastic and elastoplastic regimes; the Green-Naghdi elastoplastic formulation and the hyperelastoplasticity appear to be equivalent for small strains; the hyperelastoplastic formulation is equivalent to the mechanical model of the programs EPIM3D and DD3IMP for materials the small elastic strains regime; simulation time reduction has been obtained with the parallelization of the computer codes for structural analysis; the developed programs are capable of simulating, precisely, complex problems, such as the Cook\'s membrane and the pulled thin cylinder.

Submitted by Renata Lopes (renatasil82@gmail.com) on 2017-03-07T17:42:36Z No. of bitstreams: 1 denispiresrodriguesalves.pdf: 7315782 bytes, checksum: 6961087b46d541aeae9c4509e2d9a289 (MD5) / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2017-03-10T12:19:54Z (GMT) No. of bitstreams: 1 denispiresrodriguesalves.pdf: 7315782 bytes, checksum: 6961087b46d541aeae9c4509e2d9a289 (MD5) / Made available in DSpace on 2017-03-10T12:19:54Z (GMT). No. of bitstreams: 1 denispiresrodriguesalves.pdf: 7315782 bytes, checksum: 6961087b46d541aeae9c4509e2d9a289 (MD5) Previous issue date: 2016-09-02 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / As treliças espaciais são estruturas compostas de barras usualmente metálicas (aço ou alumínio) que são utilizadas para diversas finalidades, mas principalmente para sustentar a cobertura de grandes vãos presentes em galpões e pavilhões. O presente trabalho tem como objetivo a modelagem computacional do comportamento termo-elastoplástico de treliças espaciais em regime de grandes deslocamentos e deformações, utilizando o modelo da equação de transferência de calor e um modelo constitutivo elastoplástico modificado para incluir a influência térmica. Simulações computacionais do modelo resultante podem ser usadas para o desenvolvimento de projetos de estruturas submetidas a grandes variações de temperatura, como as que ocorrem em um incêndio. O método dos elementos finitos (MEF) foi utilizado para determinar o campo de temperaturas na seção transversal das barras, enquanto que para encontrar os deslocamentos e as deformações nas barras em função da variação térmica e do carregamento foi utilizado o método da rigidez direta. Para resolver as equações de equilíbrio não-lineares resultantes do modelo constitutivo termomecânico foi utilizado o método de Newton-Raphson. O código desenvolvido foi inicialmente validado através de simulações computacionais em estruturas simples onde variações de temperatura alteram o módulo de elasticidade, o módulo plástico e a tensão de escoamento do material e podem causar a plastificação e até a ruptura das barras. Posteriormente são apresentados e discutidos os resultados obtidos a partir de treliças mais complexas, com geometria similar às usualmente utilizadas em aplicações de engenharia, submetido a uma situação simplificada de incêndio-padrão. A alta temperatura causa a diminuição da resistência e da rigidez das barras e informações importantes como o tempo de incêndio suportado pela estrutura e o número de barras plastificadas são extraídas das simulações e podem servir como uma medida de segurança para evitar danos maiores em locais com grandes aglomerados de pessoas. / Space trusses are structures usually composed of metalic rods (steel or aluminum) that are used for several purposes, but mainly to sustain the roof of large spans present in sheds and pavilions. The present work has the aim of computationally model the coupled thermo-elastoplastic behavior of space trusses under large displacements and large strains, using the heat transfer equation model and an elastoplastic constitutive model modified to include the thermal influence. Computer simulations of the resulting mathematical model can be used for the development of structural projects under large variations of temperature, as occurs in fire situations. The finite element method (FEM) was used to determine the temperature field in transversal section of rods. In order to find the displacements and strains due to thermal variation and loadings, it was used the direct stiffness method. The Newton-Raphson method was used to solve the resulting non-linear equilibrium equations of the thermomecanic constitutive model. The developed code was initially validated through computational simulations of simple structures where thermal variations affect the Young modulus, the plastic modulus and the yield stress of the material. The results of more complex trusses, with a geometry similar to the ones usually adopted in engineering applications, under a simplified standard fire situation are also presented. The high temperature causes a decrease in the rods' resistance and stiffness and important informations such as the fire time supported by the structure and the number of plastified rods are achieved from the simulations and can be used as a security measure to avoid greater damage in places with large crowds of people.

La conception des turbomachines à haute densité énergétique nécessite de plus en plus la maîtrise d'un plus grand nombre de paramètres fonctionnels. La moindre défaillance d'un composant conduit quasi immédiatement la machine à la rupture. C'est en particulier le cas pour le comportement des composants à films minces.L'appellation "film mince" correspond à tout espace de très faible épaisseur situé entre le rotor et le stator de la turbomachine. Leur but est soit de limiter les fuites de manière à optimiser les performances intrinsèques de la machine, soit alors à supporter et stabiliser le rotor. Ces derniers cas sont plus appelés communément "joints lisses ou annulaires" et "paliers fluides".Lorsqu'un fluide circule dans un espace de très faible épaisseur, typiquement quelques centièmes de millimètres sur une distance très longue, son champ de vitesses, donc de pression, dépend fortement des phénomènes visqueux aux parois dont l'une est mise en rotation et l'autre est immobile. Les efforts fluides sur ces parois peuvent être alors importants et doivent être pris en compte dans le dimensionnement de la machine.La connaissance précise de ces écoulements très complexes est indispensable pour déterminer les efforts statiques et dynamiques appliqués au rotor de manière à pouvoir dimensionner un fonctionnement calme. / The design of high performance aerospace turbo pumps requires more control of an increasing number of functional parameters. Any component failure led almost immediately to a machine failure. This is particularly the case for the behavior of thin film lubricated components.The term "thin film" means any thin space between the rotor and the stator of the engine. Their goal is either to limit leakage to maximize the machine intrinsic performance, or to support and stabilize the rotor. These cases are more commonly called "smooth or annular seals" and "fluid film bearings".When a fluid flows in a space of very small thickness, typically a few hundredths of a millimeter, the velocity field, hence the pressure, are highly dependent on the walls viscous forces. Fluid forces on the walls (which one is rotated and the other is stationary) can then be important and should be taken into account in the design of the machine.The precise knowledge of these complex flows is essential to determine the static and dynamic forces applied to the rotor to ensure a quite functioning of the turbo pump.