The present thesis is focused on the development of a thorough mathematical modelling and computational solution framework aimed at the numerical simulation of journal and sliding bearing systems operating under a wide range of lubrication regimes (mixed, elastohydrodynamic and full film lubrication regimes) and working conditions (static, quasi-static and transient conditions). The fluid flow effects have been considered in terms of the Isothermal Generalized Equation of the Mechanics of the Viscous Thin Films (Reynolds equation), along with the massconserving p-Ø Elrod-Adams cavitation model that accordingly ensures the so-called JFO complementary boundary conditions for fluid film rupture. The variation of the lubricant rheological properties due to the viscous-pressure (Barus and Roelands equations), viscous-shear-thinning (Eyring and Carreau-Yasuda equations) and density-pressure (Dowson-Higginson equation) relationships have also been taken into account in the overall modelling. Generic models have been derived for the aforementioned bearing components in order to enable their applications in general multibody dynamic systems (MDS), and by including the effects of angular misalignments, superficial geometric defects (form/waviness deviations, EHL deformations, etc.) and axial motion. The bearing exibility (conformal EHL) has been incorporated by means of FEM model reduction (or condensation) techniques. The macroscopic in fluence of the mixedlubrication phenomena have been included into the modelling by the stochastic Patir and Cheng average ow model and the Greenwood-Williamson/Greenwood-Tripp formulations for rough contacts. Furthermore, a deterministic mixed-lubrication model with inter-asperity cavitation has also been proposed for full-scale simulations in the microscopic (roughness) level. According to the extensive mathematical modelling background established, three significant contributions have been accomplished. Firstly, a general numerical solution for the Reynolds lubrication equation with the mass-conserving p - Ø cavitation model has been developed based on the hybridtype Element-Based Finite Volume Method (EbFVM). This new solution scheme allows solving lubrication problems with complex geometries to be discretized by unstructured grids. The numerical method was validated in agreement with several example cases from the literature, and further used in numerical experiments to explore its exibility in coping with irregular meshes for reducing the number of nodes required in the solution of textured sliding bearings. Secondly, novel robust partitioned techniques, namely: Fixed Point Gauss-Seidel Method (PGMF), Point Gauss-Seidel Method with Aitken Acceleration (PGMA) and Interface Quasi-Newton Method with Inverse Jacobian from Least-Squares approximation (IQN-ILS), commonly adopted for solving uid-structure interaction problems have been introduced in the context of tribological simulations, particularly for the coupled calculation of dynamic conformal EHL contacts. The performance of such partitioned methods was evaluated according to simulations of dynamically loaded connecting-rod big-end bearings of both heavy-duty and high-speed engines. Finally, the proposed deterministic mixed-lubrication modelling was applied to investigate the in fluence of the cylinder liner wear after a 100h dynamometer engine test on the hydrodynamic pressure generation and friction of Twin-Land Oil Control Rings.

O presente trabalho é focado no desenvolvimento de uma modelagem matemática abrangente e um conjunto de ferramentas computacionais destinadas à simulação numérica de mancais hidrodinâmicos sujeitos a diferentes regimes de lubrificação (misto, elastohidrodin^amico e completamente hidrodinâmico) e condições operacionais (estático, quasi-estático e transiente). O escoamento do fluido lubrificante é descrito matematicamente através da Equação Generalizada da Mecânica dos Filmes Finos Viscosos (equação de Reynolds), juntamente com o modelo de cavitação p - Ø proposto por Elrod-Adams que incorpora automaticamente as condições de contorno complementares de JFO para a completa conservação dos fluxos mássicos na presença de ruptura de filme lubricante (cavitação). As variações das propriedades reológicas do lubrificante devido aos efeitos de piezo-viscosidade (equação de Barus e Roelands), \shear-thinning" (equações de Eyring e Carreau-Yasuda) e piezo-densidade (equação de Dowson-Higginson) são também contempladas na modelagem. Modelos genéricos de mancais de deslizamento s~ao propostos de forma a possibilitar suas aplicações em problemas de dinâmica multicorpos, assim como a consideração dos efeitos de desalinhamento angular, defeitos de superfícies (forma/ondulação, deformações EHL, etc.) e movimento axial. A flexibilidade dos mancais em situações envolvendo elastohidrodinâmica é incorporada por meio de técnicas de condensação de sistemas empregadas no MEF. A in fluência na escala macroscópica dos fenômenos de lubrificação mista é incluída na modelagem através das formulações estatísticas de Patir & Cheng para a perturbação do escoamento devido à rugosidade das superfícies, e Greenwood-Williamson/Greenwood-Tripp para o contato rugoso. Adicionalmente, um modelo determinístico é também proposto no trabalho de forma a possibilitar simulações do regime misto de lubrificação diretamente na escala (microscópica) de rugosidade. A partir da extensa modelagem matemática desenvolvida, três principais contribuições são propostas. Inicialmente, uma formulação genérica para a solução da equação de Reynolds com o modelo de cavitação p - Ø foi desenvolvida de acordo com o Método dos Volumes Finitos Baseado em Elementos (MVFbE). Este novo método permite a solução de problemas de lubrificação com geometrias complexas e discretizados por malhas não estruturadas.O método foi validado a partir da comparação com diversos exemplos da literatura, e posteriormente utilizado para a simulação de mancais texturizados. A segunda contribuição da pesquisa diz respeito à proposição de técnicas mais robustas e estáveis para a solução particionada do problema fluido-estrutura estabelecido em contatos EHL conformes. As técnicas de acoplamento apresentadas são: Método de Gauss-Seidel com Ponto Fixo (PGMF), Método de Gauss-Seidel com Aceleração de Aitken (PGMA) e Método de Quasi-Newton de Interface (IQN-ILS). O desempenho de tais métodos foi avaliado por meio de simulações dinâmicas de dois mancais de biela, o primeiro utilizado em motores de combustão interna (MCI) de alta velocidade (high-speed engine), e o segundo em motores sujeitos a elevado carregamento (heavy-duty diesel engine). Finalmente, a modelagem determinística proposta para o regime misto de lubri- ficação foi aplicada na investigação da influência do desgaste superficial de cilindros de motores de combustão interna, obtidos após 100h de teste em dinamômetro, na geração de sustentação hidrodinâmica e atrito em anéis de controle de óleo (TLOCRs).