Comportamento tribo-mecânico e desgaste adesivo de materiais em nanoescala: análises por dinâmica molecular e mecânica do contínuo.

As formulações baseadas na mecânica do contínuo, embora precisas até certo ponto, por vezes não podem ser utilizadas, ou não são conceitualmente corretas para o entendimento de fenômenos em escalas reduzidas. Estas limitações podem aparecer no estudo dos fenômenos tribológicos em escala nanométrica, que passam a necessitar de novos métodos experimentais, teóricos e computacionais que permitam explorar estes fenômenos com a resolução necessária. Simulações atomísticas são capazes de descrever fenômenos em pequena escala, porém, o número necessário de átomos modelados e, portanto, o custo computacional - geralmente torna-se bastante elevado. Por outro lado, os métodos de simulação associados à mecânica do contínuo são mais interessantes em relação ao custo computacional, mas não são precisos na escala atômica. A combinação entre essas duas abordagens pode, então, permitir uma compreensão mais realista dos fenômenos da tribologia. Neste trabalho, discutem-se os conceitos básicos e modelos de atrito em escala atômica e apresentam-se estudos, por meio de simulação numérica, para a análise e compreensão dos mecanismos de atrito e desgaste no contato entre materiais. O problema é abordado em diferentes escalas, e propõe-se uma abordagem conjunta entre a Mecânica do Contínuo e a Dinâmica Molecular. Para tanto, foram executadas simulações numéricas, com complexidade crescente, do contato entre superfícies, partindo-se de um primeiro modelo que simula o efeito de defeitos cristalinos no fenômeno de escorregamento puro, considerando a Dinâmica Molecular. Posteriormente, inseriu-se, nos modelos da mecânica do contínuo, considerações sobre o fenômeno de adesão. A validação dos resultados é feita pela comparação entre as duas abordagens e com a literatura.

Aplicação de métodos geoestatísticos no estudo das distribuições espaciais de condutividade hidráulica em áreas do Sistema Aquífero Guarani (SAG) e do Sistema Aquífero Bauru (SAB) no Estado de São Paulo

A condutividade hidráulica (K) é um dos parâmetros controladores da magnitude da velocidade da água subterrânea, e consequentemente, é um dos mais importantes parâmetros que afetam o fluxo subterrâneo e o transporte de solutos, sendo de suma importância o conhecimento da distribuição de K. Esse trabalho visa estimar valores de condutividade hidráulica em duas áreas distintas, uma no Sistema Aquífero Guarani (SAG) e outra no Sistema Aquífero Bauru (SAB) por meio de três técnicas geoestatísticas: krigagem ordinária, cokrigagem e simulação condicional por bandas rotativas. Para aumentar a base de dados de valores de K, há um tratamento estatístico dos dados conhecidos. O método de interpolação matemática (krigagem ordinária) e o estocástico (simulação condicional por bandas rotativas) são aplicados para estimar os valores de K diretamente, enquanto que os métodos de krigagem ordinária combinada com regressão linear e cokrigagem permitem incorporar valores de capacidade específica (Q/s) como variável secundária. Adicionalmente, a cada método geoestatístico foi aplicada a técnica de desagrupamento por célula para comparar a sua capacidade de melhorar a performance dos métodos, o que pode ser avaliado por meio da validação cruzada. Os resultados dessas abordagens geoestatísticas indicam que os métodos de simulação condicional por bandas rotativas com a técnica de desagrupamento e de krigagem ordinária combinada com regressão linear sem a técnica de desagrupamento são os mais adequados para as áreas do SAG (rho=0.55) e do SAB (rho=0.44), respectivamente. O tratamento estatístico e a técnica de desagrupamento usados nesse trabalho revelaram-se úteis ferramentas auxiliares para os métodos geoestatísticos.

Métodos de adequação e diagnóstico em modelos de sobrevivência dinâmicos

A análise de dados de sobrevivência tem sido tradicionalmente baseada no modelo de regressão de Cox (COX, 1972). No entanto, a suposição de taxas de falha proporcionais assumida para esse modelo pode não ser atendida em diversas situações práticas. Essa restrição do modelo de Cox tem gerado interesse em abordagens alternativas, dentre elas os modelos dinâmicos que permitem efeito das covariáveis variando no tempo. Neste trabalho, foram revisados os principais modelos de sobrevivência dinâmicos com estrutura aditiva e multiplicativa nos contextos não paramétrico e semiparamétrico. Métodos gráficos baseados em resíduos foram apresentados com a finalidade de avaliar a qualidade de ajuste desses modelos. Uma versão tempo-dependente da área sob a curva ROC, denotada por AUC(t), foi proposta com a finalidade de avaliar e comparar a qualidade de predição entre modelos de sobrevivência com estruturas aditiva e multiplicativa. O desempenho da AUC(t) foi avaliado por meio de um estudo de simulação. Dados de três estudos descritos na literatura foram também analisados para ilustrar ou complementar os cenários que foram considerados no estudo de simulação. De modo geral, os resultados obtidos indicaram que os métodos gráficos apresentados para avaliar a adequação dos modelos em conjunto com a AUC(t) se constituem em um conjunto de ferramentas estatísticas úteis para o próposito de avaliar modelos de sobrevivência dinâmicos nos contextos não paramétrico e semiparamétrico. Além disso, a aplicação desse conjunto de ferramentas em alguns conjuntos de dados evidenciou que se, por um lado, os modelos dinâmicos são atrativos por permitirem covariáveis tempo-dependentes, por outro lado podem não ser apropriados para todos os conjuntos de dados, tendo em vista que estimação pode apresentar restrições para alguns deles.

Estudo fenomenológico e numérico do escoamento estratificado óleo-água ondulado e com mistura na interface

Escoamentos bifásicos estão presentes em diversos processos naturais e industriais, como na indústria de petróleo. Podem apresentar-se em diferentes configurações topológicas, ou, padrões de escoamento, entre eles o escoamento estratificado ondulado e o estratificado com mistura na interface. Os escoamentos bifásicos estratificados óleo-água têm sido utilizados como uma forma conveniente de evitar a formação de emulsões de água em óleo em oleodutos e possuem uma ocorrência comum em poços de petróleo direcionais. Quando a onda interfacial ultrapassa determinado limite geométrico e cinemático, surge o fenômeno do entranhamento de gotas, representado por misturas entre as fases junto à interface que promovem um aumento na queda de pressão. Modelos têm sido apresentados pela literatura na tentativa de descrever o fenômeno do entranhamento de gotas. Neste trabalho é apresentada uma nova proposta de modelagem matemática unidimensional para o entranhamento de gotas com o objetivo de melhorar a previsão dos parâmetros envolvidos, em especial, da fração volumétrica de óleo e da queda de pressão bifásica. Também foi utilizada simulação numérica computacional, CFD (Computational Fluid Dynamics), com o uso de software comercial para obtenção dos valores dos parâmetros do escoamento estratificado ondulado óleo-água (fração volumétrica de óleo, queda de pressão, amplitude e comprimento da onda interfacial). Os resultados da modelagem fenomenológica para entranhamento e os de CFD foram comparados com bancos de dados experimentais. Os resultados em CFD mostram concordância com os resultados experimentais, tanto na análise qualitativa das propriedades geométricas das ondas interfaciais, quanto na comparação direta com os dados para fração volumétrica e queda de pressão. Os resultados numéricos da modelagem fenomenológica para fatores de entranhamento apresentam boa concordância com dados da literatura.

Simulação e modelagem computacional de dados de espalhamento à baixos ângulos - enfoque em estruturas de alta simetria

Esta tese apresenta uma abordagem para a criação rápida de modelos em diferentes geometrias (complexas ou de alta simetria) com objetivo de calcular a correspondente intensidade espalhada, podendo esta ser utilizada na descrição de experimentos de es- palhamento à baixos ângulos. A modelagem pode ser realizada com mais de 100 geome- trias catalogadas em um Banco de Dados, além da possibilidade de construir estruturas a partir de posições aleatórias distribuídas na superfície de uma esfera. Em todos os casos os modelos são gerados por meio do método de elementos finitos compondo uma única geometria, ou ainda, compondo diferentes geometrias, combinadas entre si a partir de um número baixo de parâmetros. Para realizar essa tarefa foi desenvolvido um programa em Fortran, chamado de Polygen, que permite modelar geometrias convexas em diferentes formas, como sólidos, cascas, ou ainda com esferas ou estruturas do tipo DNA nas arestas, além de usar esses modelos para simular a curva de intensidade espalhada para sistemas orientados e aleatoriamente orientados. A curva de intensidade de espalhamento é calculada por meio da equação de Debye e os parâmetros que compõe cada um dos modelos, podem ser otimizados pelo ajuste contra dados experimentais, por meio de métodos de minimização baseados em simulated annealing, Levenberg-Marquardt e algorítmicos genéticos. A minimização permite ajustar os parâmetros do modelo (ou composição de modelos) como tamanho, densidade eletrônica, raio das subunidades, entre outros, contribuindo para fornecer uma nova ferramenta para modelagem e análise de dados de espalhamento. Em outra etapa desta tese, é apresentado o design de modelos atomísticos e a sua respectiva simulação por Dinâmica Molecular. A geometria de dois sistemas auto-organizado de DNA na forma de octaedro truncado, um com linkers de 7 Adeninas e outro com linkers de ATATATA, foram escolhidas para realizar a modelagem atomística e a simulação por Dinâmica Molecular. Para este sistema são apresentados os resultados de Root Mean Square Deviations (RMSD), Root Mean Square Fluctuations (RMSF), raio de giro, torção das hélices duplas de DNA além da avaliação das ligações de Hidrogênio, todos obtidos por meio da análise de uma trajetória de 50 ns.

Caracterização limnológica do rio do Peixe (microrregião geográfica de Botucatu - SP), em duas épocas do ano (períodos de seca e chuva)

O rio do Peixe, um dos principais afluentes do rio Tietê (Reservatório de Barra Bonita/Hidrovia Tietê Paraná), juntamente com seus formadores (microrregião geográfica de Botucatu - SP) constituem os principais mananciais das cidades que se localizam em suas bacias. A maior parte da bacia hidrográfica do rio do Peixe desenvolve-se sob rochas sedimentares arenosas, favoráveis às atividades de extração de areia, fato que aliado ao manejo incorreto do solo favorece a ocorrência de processos erosivos. Foram realizadas amostragens de água durante dez dias consecutivos e de sedimento durante três dias alternados, em dois períodos hidrológicos (seca e chuva), para a determinação de variáveis físicas, químicas e biológicas,com o objetivo principal de caracterizar do ponto de vista limnológico esse rio e sua foz. As análises das diferentes variáveis seguiram os métodos que são utilizados rotineiramente no Laboratório de Limnologia do CRHEA/EESC/USP. As águas do rio do Peixe apresentaram, principalmente altas concentrações de sólidos em suspensão. Os maiores valores de turbidez, condutividade, alcalinidade, HCO3-, CO3-, nitrito, nitrogênio amoniacal, nitrogênio Kjedhal, fosfato inorgânico, carbono orgânico e sólidos em suspensão foram obtidos durante o período de chuvas, influenciados principalmente pelo escoamento superficial que promove o carreamento de materiais da bacia hidrográfica. Foi possível ainda através das características limnológicas, agrupar as estações de amostragem ao longo do contínuo do rio em três regiões (alto, médio e baixo curso), e diferenciar do ponto de vista trófico as estações da foz com o rio Tietê. Um aumento da biomassa fitoplanctônica, no período de seca, no baixo curso do rio do Peixe foi favorecido pela menor vazão na barragem de Barra Bonita.

Considerações sobre os métodos de dimensionamento estrutural de pavimentos com camadas recicladas a frio com espuma de asfalto.

A reciclagem profunda com espuma de asfalto tem sido uma alternativa de sucesso para a restauração de pavimentos degradados. Em relação às soluções tradicionais de reabilitação, como os recapeamentos, tem a vantagem de proporcionar a correção de defeitos em camadas inferiores, com a manutenção ou pequena elevação do greide da pista, além de ganhos ambientais, como um menor consumo de materiais virgens da natureza e redução do volume de material descartado. Entretanto, no Brasil não há método para dimensionamento estrutural para esta tecnologia, o que dificulta seu emprego. Para o desenvolvimento de um procedimento de dimensionamento que contemple este tipo de solução, foram estudados métodos presentes na bibliografia internacional: guia da AASHTO de 1993 e Caltrans, dos EUA, TRL386 e TRL611, da Inglaterra, as duas versões do guia sul-africano TG2 e os métodos oriundos do Austroads, tanto o procedimento interino de 2011 como adaptações de órgãos da Austrália e Nova Zelândia. Observou-se divergência de opiniões quanto ao comportamento do material reciclado com espuma de asfalto. Alguns órgãos e autores consideram o comportamento do mesmo mais próximo às misturas asfálticas, sendo o mecanismo de falha o trincamento, e outros o definem como semelhante a um material granular modificado com alta coesão e ruptura devido às deformações permanentes. Correlaciona-se tal associação ao teor de espuma usualmente utilizado nas obras rodoviárias. Outros aspectos que se destacam para este tipo de base são o ganho de resistência ao longo do tempo devido à cura, mesmo com início da operação da rodovia e a importância da infraestrutura remanescente no dimensionamento. Tais fatos foram corroborados pelos estudos de caso e resultados do trecho experimental construído na Rodovia Ayrton Senna - SP 070, monitorado por meio de ensaios deflectométricos com FWD durante um ano. Como resultado do trabalho, foi proposto um procedimento para o dimensionamento estrutural de pavimentos com camadas recicladas a frio com espuma de asfalto utilizando dados deflectométricos que atende o método do Manual de Pavimentação do Departamento Nacional de Infraestrutura de Transportes (DNIT) e incorpora diferentes aspectos na análise mecanicista. Outras conclusões são a viabilidade técnica a longo prazo da solução mencionada e a importância do controle tecnológico, com ênfase para o monitoramento deflectométrico nos primeiros meses de operação do pavimento para averiguar a evolução da cura do material.

Avaliação do comportamento mecânico de corpos de prova de misturas asfálticas a quente resultantes de diferentes métodos de compactação.

Os procedimentos de dosagem Marshall e Superpave definem os teores de ligante de projeto baseados em parâmetros volumétricos. Nessa situação, sistemáticas de dosagem com tipos de compactação diferentes podem conduzir a teores de ligante de projeto distintos que definirão a vida útil dos revestimentos asfálticos. O objetivo principal desse trabalho é avaliar o comportamento mecânico de misturas asfálticas moldadas por diferentes métodos de compactação de laboratório e analisar a relação com os resultados de amostras obtidas a partir de misturas compactadas por rolagem pneumática na mesa compactadora francesa. A fase experimental consistiu na dosagem de misturas pelos métodos Marshall e Superpave (este último com dois tamanhos de moldes), além da compactação na Prensa de Cisalhamento Giratória (PCG) e da moldagem de placas na mesa compactadora. Avaliou-se o efeito do tipo de compactação, do tamanho do molde e do número de giros do Compactador Giratório Superpave (CGS) no teor de projeto, nos parâmetros volumétricos, no comportamento mecânico e no desempenho quanto à fadiga e à resistência ao afundamento em trilha de roda. Adicionalmente, foi avaliada a eficiência do método Bailey de composição granulométrica quanto à resistência à deformação permanente em função do tipo de agregado. Constatou-se que o método Bailey, por si só, não garante resistência à deformação permanente, sendo essa dependente do tipo de agregado incluindo seus parâmetros de forma. O principal produto da pesquisa, com efeitos práticos no projeto de misturas asfálticas, traduz-se na recomendação do método Superpave com molde de 100 mm (para TMN <= 12,5 mm) para volume de tráfego médio a alto em detrimento ao método Superpave com 150 mm, tendo em vista que o primeiro apresenta densificação mais semelhante às amostras preparadas na compactação por rolagem (similar ao que ocorre em pista) o que resulta em comportamento mecânico também mais próximo da realidade de campo. A utilização dos moldes de 150 mm de diâmetro no CGS pode ser viabilizada desde que se adote um número de giros menor do aquele proposto para projeto pelo Asphalt Institute (2001). Por fim, é fundamental que os ensaios e os cálculos para obtenção dos parâmetros volumétricos e escolha do teor de projeto sigam ao normatizado pela ASTM, pelo Asphalt Institute (2001) e pela ABNT.

Resistência à fadiga de tubo API 5L X65 cladeado e soldado circunferencialmente com eletrodos de Inconel® 625

As recentes descobertas de petróleo e gás na camada do Pré-sal representam um enorme potencial exploratório no Brasil, entretanto, os desafios tecnológicos para a exploração desses recursos minerais são imensos e, consequentemente, têm motivado o desenvolvimento de estudos voltados a métodos e materiais eficientes para suas produções. Os tubos condutores de petróleo e gás são denominados de elevadores catenários ou do inglês \"risers\", e são elementos que necessariamente são soldados e possuem fundamental importância nessa cadeia produtiva, pois transportam petróleo e gás natural do fundo do mar à plataforma, estando sujeitos a carregamentos dinâmicos (fadiga) durante sua operação. Adicionalmente, um dos problemas centrais à produção de óleo e gás das reservas do Pré-Sal está diretamente associado a meios altamente corrosivos, tais como H2S e CO2. Uma forma mais barata de proteção dos tubos é a aplicação de uma camada de um material metálico resistente à corrosão na parte interna desses tubos (clad). Assim, a união entre esses tubos para formação dos \"risers\" deve ser realizada pelo emprego de soldas circunferenciais de ligas igualmente resistentes à corrosão. Nesse contexto, como os elementos soldados são considerados possuir defeitos do tipo trinca, para a garantia de sua integridade estrutural quando submetidos a carregamentos cíclicos, é necessário o conhecimento das taxas de propagação de trinca por fadiga da solda circunferencial. Assim, neste trabalho, foram realizados ensaios de propagação de trinca por fadiga na região da solda circunferencial de Inconel® 625 realizada em tubo de aço API 5L X65 cladeado, utilizando corpos de prova do tipo SEN(B) (Single Edge Notch Bending) com relações entre espessura e largura (B/W) iguais a 0,5, 1 e 2. O propósito central deste trabalho foi de obter a curva da taxa de propagação de trinca por fadiga (da/dN) versus a variação do fator de intensidade de tensão (&#916;K) para o metal de solda por meio de ensaios normatizados, utilizando diferentes técnicas de acompanhamento e medição da trinca. A monitoração de crescimento da trinca foi feita por três técnicas: variação da flexibilidade elástica (VFE), queda de potencial elétrico (QPE) e análise de imagem (Ai). Os resultados mostraram que as diferentes relações B/W utilizadas no estudo não alteraram significantemente as taxas de propagação de trinca por fadiga, respeitado que a propagação aconteceu em condições de escoamento em pequena escala na frente da trinca. Os resultados de propagação de trinca por fadiga permitiram a obtenção das regiões I e II da curva da/dN versus &#916;K para o metal de solda. O valor de &#916;Klim obtido para o mesmo foi em torno de 11,8 MPa.m1/2 e os valores encontrados das constantes experimentais C e m da equação de Paris-Erdogan foram respectivamente iguais a 1,55 x10-10 [(mm/ciclo)/(MPa.m1/2)m] e 4,15. A propagação de trinca no metal de solda deu-se por deformação plástica, com a formação de estrias de fadiga.

Investigação de camadas aqüíferas por métodos geoelétricos no Estado de São Paulo

No Estado de São Paulo, a maioria dos poços profundos perfurados localiza-se em terrenos sedimentares. Assim sendo foi realizado este trabalho a fim de verificar a aplicabilidade de métodos geofísicos à procura de água subterrânea no Estado. A água da chuva atingindo a superfície do terreno toma caminhos a saber: evaporação, escoamento e infiltração. Em face da porosidade, permeabilidade e força da gravidade, a água de infiltração toma o caminho descendente até: - atingir uma zona permeável (saturada ou não) e seu fluxo. - atingir uma rocha sedimentar impermeável ou rocha ígnea sã. - saturar fissuras existentes em rochas ígneas sãs. Basicamente que um poço seja produtor de água, é necessária a extração de água acumulada nas seguintes condições: 1- água acumulada em rochas permeáveis saturadas. 2- água acumulada em fissuras de rochas ígneas sãs. 3- água acumulada em depressões existentes ou na superfície de rochas ígneas ou do embasamento cristalino, coberto por rochas porosas e permeáveis. Baseado nestas formas de ocorrência de água subterrânea, as pesquisas de campo foram orientadas no sentido de: - identificar pela interpretação das curvas ou mesmo uma curva de sondagem elétrica, aquela camada aqüífera, saturada de água e já conhecida \"a priori\" como aqüífero, bem como estabelecer a sua profundidade, espessura e comportamento geológico. - verificar a profundidade do nível freático. - mapear a superfície do embasamento cristalino, procurando achar nela, depressões causa das por erosões pretéritas. Procurou-se verificar também: até que ponto, em um dado local uma estratigrafia geológica conhecida, corresponde a uma \"estratigrafia geofísica\" obtida por medidas na superfície. - determinar a resistividade, ou resistividades, das diferentes formações litológicas. - procurar correlacionar variações laterais de resistividade de uma certa camada geológica com as variações na sua granulometria e litologia ou sais ) dissolvidos na solução que a satura. De um modo geral a eletrorresistividade fornece bons resultados quando aplicada à procura de água subterrânea no Estado de São Paulo. Nos sedimentos do Grupo Tubarão (Permo-Carbonífero) verifica-se uma correspondência entre altos valores de resistividade e as altas vazões dos poços (vazão específica). A \"estratigrafia geofísica\" obtida corresponde muito bem à estratigrafia geológica. Espessuras dos derrames de basalto (Eocretáceo), intrusões de diabásio, arenito da Série Bauru (Neocretáceo), arenito da Formação Botucatu (Eocretáceo), sedimentos do Grupo Tubarão (Permo-Carbonífero) e camadas de sedimentos recentes são facilmente determinadas. Na região da Praia Grande existem lençóis de água doce sobrejazendo à água salobra e salgada.

Contribuição ao conhecimento de caracteristicas tecnologicas de materiais rochosos

Neste trabalho foi desenvolvida uma investigação experimental com o objetivo de ampliar o conhecimento sobre caracterização tecnológica de materiais de construção, com ênfase no aspecto da alternabilidade. Foram estudadas rochas de 15 localidades representando os seguintes 11 tipos litológicos: basalto compacto, basalto vesicular, basalto microvesicular, diabásio, quartzito, calcário, granito, charnockito, piroxenito, migmatitoe arenito. Os ensaios de laboratório foram conduzidos em corpos de prova cilídricos, cúbicos e britados obtidos a partir das amostras enviadas, para determinação dos seguintes índices: massa específica aparente, porosidade aparente, absorção, abrasão Los Angeles, resistência ao impacto em corpos de prova regulares, resistência ao impacto Treton, resistência ao esmagamento, resistência à compressão simples axial, índice escleromético, resistência ao desgaste a úmido e durabilidade. Foram executados quase 3000 ensaios para obtenção de índices fisicomecânicos e de durabilidade, por um período de aproximadamente 2 anos. Todos os ensaios foram efetuados em amostras inicialmente lavadas, secas em estufa a 100-110°C por aproximadamente 10 horas e posteriormente, esfriadas à temperatura ambiente durante pelo menos 2-3 horas, antes de cada ensaio. Os seguintes ensaios novos foram introduzidos: índice esclerométrico, impacto em corpos de prova moldados e desgaste úmido, com discussões a respeito de sua aplicabilidade. Os resultados de cada grupo de corpos de prova em diferentes rochas, para estes ensaios novos, foram submetidos à análise estatística para obtenção de valores médios e coeficientes de variação para determinar a precisão e o número de ensaios necessários para cada uma destas propriedades. Os resultados médios foram submetidos à correlação estatística e análises de regressão pelo método dos mínimos quadrados, para os quais foram determinados os coeficientes de correlação ) simples, erros padrão de estimativa e relações funcionais para cada par selecionado de variáveis. Estas relações foram apresentadas graficamente para comparação visual com os correspondentes resultados de ensaios. A partir dos resultados obtidos quanto à precisão e reprodutibilidade dos ensaios e verificada sua aplicabilidade, foram elaboradas propostas para métodos de ensaios. Os ensaios de durabilidade, que constituem a segunda parte deste trabalho, foram executados nas rochas acima mencionadas, cujas amostras foram submetidas a ciclos sucessivos de, saturação em solução de sulfato e secagem em estufa, saturação em etilenoglicol e secagem em estufa, saturação em água e secagem em estufa, lixiviação contínua no exterior Soxhlet e exposição natural, em 2 séries de ensaios. A primeira série consistiu na comparação de processos de alteração e a segunda série, na investigação do comportamento mecânico dos materiais rochosos face a diferentes condições fisicoquímicas de alteração artificial e natural. P parâmetro mecânico utilizado neste caso, foi a resistência ao impacto Treton. Os resultados obtidos na primeira série, foram submetidos à análise de regressão, obtendo-se gráficos de comparação de linhas teóricas de velocidades de alteração pelos diferentes processos e, gráficos mostrando a influência do fator de escala no processo de alteração. Os resultados da segunda série de ensaios, permitiram estabelecer 2 parâmetros numéricos: queda percentual de resistência mecânica e coeficiente de alterabilidade específica, que refletem a suscetibilidade das rochas a estes diferentes processos de alteração. A correlação estatística, entre estes parâmetros numéricos da alteração e a relação entre as tensões máximas de ruptura, no estado saturado e no estado seco, do ensaio de compressão simples axial, indicou ser esta relação, um índice promissor da alterabilidade de rochas. Finalmente, discutida a validade ) destes índices, tendo em vista os resultados encontrados e a ausência na literatura de procedimentos adequados, que permitam classificar as rochas quanto ao comportamento mecânico face à alteração natural ou artificial, foram elaboradas propostas de métodos de ensaios, que poderão preencher esta lacuna.

Estudo numérico do efeito de blocagem hidrodinâmica em um túnel de cavitação sobre o escoamento ao redor de um submarino.

O projeto de submarinos é um tema que vem se desenvolvendo há muitos anos. No Brasil, vivemos um momento de pleno avanço na área, sendo encabeçado pela Marinha do Brasil com o desenvolvimento do primeiro submarino de propulsão nuclear brasileiro. No decorrer desse projeto, diversos estudos devem ser desenvolvidos, dentre eles, ensaios com modelos em escala, que serão realizados em um túnel de cavitação. Durante estes experimentos, espera-se que surjam diversas fontes de incerteza, que podem comprometer os resultados obtidos, destacando-se o efeito de blocagem hidrodinâmica. Nesse contexto, esta dissertação apresenta uma investigação do fenômeno de blocagem no modelo de submarino DARPA SUBOFF na seção de testes de um túnel de cavitação. Para isso, simulações tridimensionais de média de Reynolds (RANS) foram realizadas para o modelo de submarino posicionado no túnel (sob blocagem) e sob condição de águas profundas (referência), para diferentes ângulos de ataque e número de Reynolds (Re) constante. Coeficientes de força e momento hidrodinâmicos foram obtidos e os métodos de correção de blocagem clássico e WICS foram empregados para obter fatores que pudessem corrigir esses valores. A eficiência de cada método foi avaliada comparando os valores corrigidos com os de referência. Os resultados mostraram que ambos métodos de correção são eficientes na recuperação de coeficientes de força e momento hidrodinâmicos para razões de blocagem inferiores a 10%, mas ressalvas são feitas quanto à utilização de cada um desses métodos, sendo o método clássico indicado para o estudo de geometrias simples e o método WICS para modelos mais complexos. Paralelamente, a topologia do escoamento ao redor do submarino foi estudada através do método dos pontos críticos, para avaliar se a blocagem hidrodinâmica pode causar alterações nas estruturas presentes no escoamento. A análise topológica permitiu avaliar a evolução nas estruturas do escoamento para os ângulos de ataque e Re simulados. Foram capturadas as linhas de separação no corpo e avaliada sua relação com a tensão de cisalhamento ao longo do casco. A blocagem hidrodinâmica não causou alterações drásticas na topologia do escoamento.

On the development of advanced techniques for mixed-elastohydrodynamic lubrification modelling of journal and sliding bearing systems.

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.