Les bruits automobiles induits par le frottement sont à l'origine de nombreuses plaintes clients et occasionnent des coûts de garantie considérables pour les constructeurs automobiles. Les objectifs de la thèse consistent à comprendre la physique à l'origine de ces bruits et proposer des méthodologies de simulation afin de les éradiquer. Un système générique est tout d'abord étudié. Ce système discret met en jeu un contact entre deux masses et une loi de frottement de Coulomb présentant une discontinuité à vitesse relative nulle. Des calculs de valeurs propres complexes de ce système linéarisé autour de sa position d'équilibre glissant sont menés et montrent la présence d'instabilités par flottement voire par divergence. Les simulations temporelles montrent quant à elles que les non-linéarités de contact permettent de stabiliser les niveaux vibratoires en cas d'instabilité selon quatre régimes distincts. De plus, malgré ses trois degrés de liberté, ce système est capable de reproduire les mécanismes de stick-slip, sprag-slip et couplage modal ainsi que les bruits de crissement, grincement et craquement rencontrés sur les systèmes automobiles. Des études paramétriques sont également présentées et mettent en avant des bifurcations de Hopf ainsi que l'effet déstabilisant potentiellement induit par l'amortissement. Des méthodologies permettant de catégoriser les réponses en termes de bruit et de mécanisme sont par la suite proposées. Les occurrences et risques de ces derniers sont alors analysés et des tendances sont dégagées. Enfin, la relation entre les bruits et les mécanismes est établie. L'attention est ensuite portée sur un système automobile particulier. Afin d'étudier son comportement crissant, les analyses de stabilité et les simulations temporelles sont désormais menées sur des modèles éléments-finis. Les simulations temporelles permettent d'observer l'établissement de vibrations auto-entretenues et d'identifier, parmi tous les modes instables prédits lors des analyses de stabilité, celui qui est réellement à l'origine de l'instabilité. L'effet du coefficient de frottement sur les motifs de coalescence et les cycles limites est également investigué. Le risque de crissement est ensuite évalué pour des conditions d'utilisation variées du système. La méthodologie, basée sur des analyses de stabilité, permet de retrouver les principaux constats expérimentaux obtenus sur banc d'essai. Le rôle des géométries et des matériaux constituant le système est également discuté. Enfin, une solution permettant de réduire de façon significative le risque de crissement est proposée.

[pt] Uma prospecção bem sucedida de petróleo e gás requer muitos esforços para se sobrepor os desafios encontrados, tais como vibrações axiais, laterais e torcionais. Estes fenômenos podem causar a falha prematura de componentes do sistema de perfuração, disfunção nos equipamentos de medição e aumento no tempo e custo no processo de perfuração. Em particular, vibrações torcionais estão presentes em grande parte dos processos de perfuração e podem alcançar um estado crítico: stick-slip. Um melhor entendimento sobre este fenômeno proporciona ferramentas para evitar o aumento do tempo e do custo da prospecção, assegurando o investimento e sucesso do processo de perfuração. Neste trabalho, é descrito um procedimento experimental com um atrito não linear objetivando induzir stick-slip e é feito uma modelagem analítica simples do problema. O modelo de atrito é baseado em um atrito seco imposto por um dispositivo de freio desenvolvido. O comportamento não linear da bancada experimental é analisada e o modelo numérico é validado comparando diagramas de bifurcações numérica e experimentais. / [en] A successful oil and gas prospecting requires many efforts to overcome the encountered challenges, some of these challenges include drill string axial, lateral and torsional vibrations. These phenomena may cause premature component failures of the drilling system, dysfunction of measurement equipments, and increase time and costs of the prospecting process. Torsional vibrations are present in most drilling processes and may reach a severe state: stick-slip. An improved understanding about the stick-slip phenomenon provides tools to avoid the increase of prospecting time and costs, assuring the investment and success of the drilling process. Firstly, a numerical analysis of the drill string is performed with different friction models. These models are proposed in order to get familiar with the drill string dynamics. Also, it is described the experimental procedure with a nonlinear friction aiming to induce stick-slip and is performed a simple analytical modeling of the problem. The friction model is based on dry friction imposed by a break device. The nonlinear behavior of the experimental apparatus is analyzed and the numerical model is validated comparing experimental and numerical bifurcation diagrams.

Increasing environmental awareness has driven a lot of research to look into various environmentally friendly lubricants which can replace more conventional mineral oil based lubricants. This is true in particular for marine and hydropower applications where the risk of lubricant leakage can be damaging to the local environment. Glycerol is an organic compound produced as a byproduct when producing Biodiesel. It is environmentally friendly and has been used by the pharmaceutical and food industries for a long time. Recent research has shown that glycerol has very good tribological properties and can be used as a lubricant in some applications. However, its high viscosity and high static friction are disadvantageous and needs improvement. The purpose of this study is to investigate possible surface active environmentally adapted additives in glycerol and water solutions. 14 additives have been selected for testing and were tested in a pin-on-disc start-stop friction test to see their effect on the static and dynamic friction coefficients. From these 14 additives two were selected due to their superior performance and were further studied. The effect of concentration of water and additive in glycerol was also investigated using the same test conditions. A phosphor based additive has been found to be the best performing.

Ces travaux portent sur la commande et l’analyse de la stabilité d’un système électropneumatique constitué d’un axe linéaire commandé par deux servodistributeurs régulant le débit massique entrant dans chaque chambre de l’actionneur. La problématique générale est motivée par l’apparition d’un phénomène de redécollage sur ce système électropneumatique difficilement pris en compte par les études actuelles en automatique. Ce problème, rencontré depuis de nombreuses années, concerne toutes les commandes linéaires et non linéaires mono et multidimensionnelles étudiées au laboratoire. Il se traduit par des mouvements saccadés du vérin au voisinage de l’équilibre. Ce phénomène est dû à la présence de frottements secs et aux dynamiques des pressions dans les chambres pneumatiques de l’actionneur, qui continuent à évoluer (intégrer le débit massique entrant délivré par les servodistributeurs), même après l’équilibre mécanique. La première partie de ce mémoire propose une commande non linéaire commutée afin d’éviter le phénomène de redécollage de l’actionneur électropneumatique notamment vis-à-vis des variations de frottements secs qui peuvent à tout moment causer ce phénomène. Cette technique est finalement mise en œuvre et son efficacité est constatée. La plus grande partie de ce mémoire traite l’analyse de l’actionneur électropneumatique avec sa loi de commande commutée. La présence de frottements secs et l’application d’une loi de commande commutée nous a amené à concilier une démarche d’analyse de stabilité, en considérant une classe de systèmes commutés appelée systèmes affines par morceaux. La principale difficulté de cette démarche réside dans l’obtention de fonctions de Lyapunov adéquates, qui se transforme en un problème d’optimisation sous contraintes LMI (Linear Matrix Inequality) en utilisant la S-procédure. Afin d’analyser la stabilité d’un système PWA (PieceWise Affine), la première démarche proposée permet le calcul d’une fonction de Lyapunov quadratique par morceaux sous la forme d’un problème d’optimisation sous contraintes LMI, en imposant des conditions suffisantes de stabilité. Ces dernières permettent, contrairement aux méthodes classiques, d’assurer la convergence de trajectoires d’état non pas vers un point d’équilibre, mais vers un ensemble des points d’équilibre d’un système PWA. L’approche proposée permet aussi l’étude de la robustesse vis-à-vis des variations paramétriques dans le système. Nous proposons aussi une deuxième approche pour la construction d’un type de fonctions de Lyapunov dites polynomiales par morceaux, via l’utilisation des "sum of square" et de la "power transformation", afin d’analyser la stabilité d’un ensemble de points d’équilibre d’un système PWA, en présence de phénomènes de glissement et de variations paramétriques. Cette approche propose des conditions suffisantes moins conservatives que celles imposées par les fonctions de Lyapunov quadratique par morceaux. En effet, sur des exemples de systèmes PWA présentant de dynamiques discontinues sur les frontières entre les cellules, pouvant générer à tout moment des phénomènes de glissement, ces dernières s’avèrent inefficaces et ne permettent pas d’assurer la stabilité des systèmes PWA en présence de ces phénomènes. Par conséquent, les résultats sur la fonction de Lyapunov quadratique par morceaux sont étendus pour pouvoir calculer des fonctions de Lyapunov polynomiales par morceaux d’ordre supérieur, en résolvant un problème d’optimisation sous contraintes LMI. Ces dernières permettent de garantir des conditions plus générales et moins conservatives par rapport à celles développées dans la littérature. Ces deux approches ont été appliquées afin d’analyser la stabilité de l’ensemble des points d’équilibre du système électropneumatique, en considérant à la fois un modèle de frottements sous la forme d’une saturation et un autre sous la forme d’un relais présentant une dynamique discontinue. [...] / This work focuses on the control and stability analysis of an electro-pneumatic system, i.e. a linear pneumatic cylinder controlled by two servo valves regulating the mass flow entering each chamber of the actuator. The general problem is motivated by the appearance of stick-slip on the electro-pneumatic system, hardly taken into account by the current studies in automatic control. This problem, encountered throughout the years, concerns all mono- and multidimensional linear and non-linear controls systems studied at the laboratory. In pneumatic cylinders, the phenomenon consists in a displacement of the rod a while after it has come to a rest ; this is due to the fact that the force acting on the rod initially becomes smaller that the threshold which is necessary for a motion, and then this threshold is overcome later on. In this case, stick-slip is caused by the presence of dry friction and by the pressure dynamics in the chambers, which continue to evolve (integrating the net incoming mass flow from the servovalves) even after the rod has stopped. The first part of this thesis proposes a nonlinear switching control law in order to avoid stick-slip on pneumatic cylinder, taking into account with the variations of dry friction that may occur at any time causing this phenomenon. This technique is implemented and its effectiveness is recognized. The greatest part of this thesis deals with the stability analysis of the pneumatic cylinder with its switched control law. The presence of dry friction and the application of a switched control law requires an appropriate method for approaching the stability analysis ; this method is based on considering the closed-loop system as belonging to a class of switched systems called piecewise affine systems (PWA). The main difficulty in this approach lies in obtaining adequate Lyapunov functions for proving stability, which turns into an optimization problem under LMI constraints (Linear Matrix Inequality) using the S-procedure. In order to analyze the stability of a PWA system, a first method is proposed allowing the computation of a piecewise quadratic Lyapunov function through an optimization problem under LMI constraints. The methods takes into account, in contrast to conventional methods, that the states might converge not to a single point but to a set of equilibrium points. The proposed approach allows also the study of robustness with respect to parametric variations in the system. A second method is also proposed for the construction of a type of Lyapunov functions called piecewise polynomial, using the “sum of squares” and “power transformation” techniques. This approach proposes less conservative sufficient conditions than those imposed by the piecewise quadratic Lyapunov functions, yielding a more succesfull stability test when for PWA systems featuring sliding modes and parametric variations. In fact, on PWA systems with discontinuous dynamics (which can generate sliding phenomena), piecewise quadratic Lyapunov functions might prove ineffective to prove the stability. Therefore, the results on piecewise quadratic Lyapunov functions are extended in order to compute piecewise polynomial Lyapunov functions of higher order, by solving an optimization problem under LMI constraints. These functions are more general and allow less conservative conditions compared to those formerly developed in the literature. Both of these methods have been applied to the stability analysis of the set of equilibrium points of the pneumatic cylinder, considering first a friction model in saturation form and then a model in relay form with a discontinuous dynamics. The application of the methods is successful, i.e. the robust stability is proven under dry friction threshold variations, with possibility of sliding modes.

Road vehicle tyres, railway wheels and ball bearings all generate rolling contact forces which are transferred within a finite area of contact between the rolling element and the substrate. Either it is visible or not for the human eye, a certain degree of roughness is always present on the contacting surfaces and it influences the generation of both vertical and lateral contactforces. The purpose of this investigation is to enhance the understanding and modelling of the influence from small-scale surface roughness on the generation of rolling contact forces. To this end, a computationally efficient method to include roughness-induced contact nonlinearities in the dynamic modelling of rolling contacts is proposed. The method is implemented in a time domain model for vertical wheel–track interaction to model rolling-induced rail vibrations, showing good agreement with measurements. Furthermore, a test rig is developed and used for the investigation of tyre–road rolling contact forces. Detailed studies are performed on the influence of substrate roughness on the resulting contact forces for a tyre tread block which is rolling at different operating conditions. The choice of substrate as well as the rolling velocity and the slip ratio is observed to have significant influence on the resulting friction coefficient. For high slip ratios, stick–slip oscillations appear, exhibiting frequency content which is largely dependent on the choice of substrate. The outcomes of this study can potentially be used to improve future tyre–road contacts with respect to wear, traction and noise generation. / <p>QC 20161013</p> / Centre for Eco2 Vehicle Design

[pt] Os custos de uma campanha exploratória, incluindo a perfuração de poços de petróleo estão cada vez mais elevados. As vibrações nas colunas de perfuração estão presentes a todo momento nas operações de perfuração, o que podem ocasionar a quebra destas colunas e do BHA (Bottom Hole Assembly), danos em brocas e demais equipamentos. Além do custo do possível dano e/ou perda destes equipamentos, existem os custos da sonda parada e a despesas operacionais das companhias contratadas, o chamado custo de Stand by, e o mais indesejável, a perda do poço. Fenômenos como Bit Bounce, Stick-Slip e Forward/Backward Whril, que advém das vibrações axiais, torcionais e laterais respectivamente, são motivos de vários estudos. Esta dissertação pretende fazer a validação do modelo matemático de dinâmica torcional utilizado na dissertação de mestrado de CAYRES (2013) para gerar mapas de estabilidade para vibrações torcionais. Nesta dissertação parâmetros reais de uma coluna de perfuração foram inseridos no modelo com objetivo de obter mapas de estabilidade de vibrações torcionais para a perfuração de um poço na Bacia de Solimões. Tais mapas de estabilidade foram gerados a partir deste modelo matemático, que na tese de CAYRES (2013), é baseado em procedimento experimental, com um atrito não linear induzindo o fenômeno stick-slip que está associado às vibrações torcionais. Este fenômeno tem relação entre torque e velocidade angular na broca provocando vibrações indesejáveis na coluna de perfuração. A geração destes mapas permitiu que a análise de sensibilidade e a averiguação da ocorrência ou não de vibrações torcionais fossem feitas, bem como a inserção de valores de WOB (peso sobre a broca) e RPM (velocidade de rotação) de um caso real de perfuração de um poço na Bacia do Solimões. Nesta análise de sensibilidade foram considerados, além dos aspectos acima citados, o tipo de broca e formação geológica, permitindo assim a validação do modelo de CAYRES (2013) e contribuindo para a otimização dos processos operacionais da perfuração de poços de petróleo. / [en] The costs of an exploration campaign, including the drilling of oil wells are increasingly high. The vibrations in the drill strings are present at all times in drilling operations that can cause damage of these drill pipe and BHA (Bottom Hole Assembly), as drill bits and other equipment. Besides the cost of a possible damage or loss of such equipment, there are the costs when the rig is not operating and the expenses of all companies involved in the operations, called the cost of stand by, and the worst case scenario of the undesirable loss of the well. The Phenomena like Bit Bounce, Stick-Slip and Forward/Backward Whril, that arises from axial, torsional and lateral vibrations respectively are the motivation of several studies. This dissertation aims to show the validation of the mathematical modeling of the torsional dynamics used in CAYRES (2013) that generate stability maps for torsional vibrations. In this dissertation real parameters from a drill string were inserted in this mathematical modeling, aiming to obtain stability maps for torsional vibrations during a drilling operation of a well in Solimões Basin. Such stability maps were generated from a dynamic model that in CAYRES (2013) thesis is based on an experimental procedure with a nonlinear friction inducing stick-slip phenomenon associated to torsional vibrations. This phenomenon is related to torque and angular velocity on the bit generating undesirable vibrations. Build stability maps and checking when torsional vibrations happen, and then comparing with a real case of WOB (weight on bit) and RPM (rotary speed) in Solimões Basin. In this sensitivity analysis were considered, besides the aspects mentioned above, the kind of bit and geological formation in order to contribute to the optimization of the oil well drilling operational process.

Immediate occupancy and functionality of critical facilities including hospitals, emergency operations centers, communications centers, and police and fire stations is of utmost importance immediately after a damaging earthquake, as they must continue to provide fundamental health, emergency, and security services in the aftermath of an extreme event. Although recent earthquakes have proven the acceptable performance of the structural system in such buildings, when designed according to recent seismic design codes, in many cases damage to the nonstructural components and systems was the main cause of disruption in their functionality. Seismic isolation is proven to be an effective technique to protect building structures from damaging earthquakes. It has been the method of choice for critical facilities, including hospitals in Japan and the United States in recent years. Seismic isolation appears to be an ideal solution for protecting the nonstructural components as well. While this claim was made three decades ago, the supporting research for freestanding (unanchored) equipment and contents (EC) is fairly new. With the focus on freestanding EC, this study investigates the seismic performance of sliding and wheel/caster-supported EC in fixed-base and base-isolated buildings. The study adopts a comparative approach to provide a better understanding of the advantages and disadvantages of using each structural system. The seismic response of sliding EC is investigated analytically in the first part of the thesis, while the response of EC supported on wheels/casters is examined through shake table experiments on two pieces of hospital equipment. The study finds base isolation to be generally effective in reducing seismic demands on freestanding EC, but it also exposes certain situations where isolation in fact increases demands on EC. Increasing the frictional resistance for sliding EC or locking the wheel/casters in the case of wheel/caster-supported EC is highly recommended for EC in base-isolated buildings to prevent excessive displacement demands. Furthermore, the study suggests several design probability functions that can be used by practicing engineers to estimate the peak seismic demands on sliding and wheel/caster-supported EC in fixed-base and base-isolated buildings. / Dissertation / Doctor of Philosophy (PhD)

Polymer electrolyte membrane (PEM) fuel cell (FC) stacks require elastomeric gaskets for each cell to keep the reactant gases within their respective regions[1]. If any gasket degrades or fails, the reactant gases can leak or mix with each other directly during operation or standby, affecting the overall operation and performance of the FC. The elastomeric gaskets used as FC seals are exposed to a range of environmental conditions, and concurrently, subjected to mechanical compression between the bipolar plates forming the cell. The combination of mechanical stress and environmental exposure may result in degradation of the seal material[2] over a period of time. In order to address the durability and make reliability predictions, the long-term stability of the gaskets in FC assemblies is critical. The aim of this study is to investigate the performance of elastomeric seals in a simulated FC environment in the presence of mechanical stresses. The overall scope of the study includes mechanical and viscoelastic properties characterization, and lifetime durability predictions based on an accelerated characterization approach. With the help of finite element analysis software, ABAQUS, a fixture was designed to perform strain-based accelerated characterization of seal material in air, deionized (DI) water, 50v/50v ethylene glycol/water solution, and 0.1M sulfuric acid solution. Dogbone samples were strained to different levels in the custom fixture and submerged in liquid solutions at 90°C and in air at 90°C and 120°C. It was observed that mechanical properties such as tensile strength, strain to break, 100% modulus, crosslink density, and tensile set degrade due to aging and the extent of change (increase or decrease) depends significantly on the strain level on the specimen. Trouser tear tests were conducted on reinforced specimens in air and deionized water (DI) to evaluate the tear resistance of an elastomeric seal material intended for proton exchange membrane fuel cells. Plots relating the crack growth rate with tearing energy were obtained at various temperatures and provided significant insight into the rate and temperature dependence of the tearing strength of the seal material. Stick-slip crack propagation was observed at all temperatures and loading rates, although the behavior was suppressed significantly at low loading rates and high temperatures. Crack growth rate versus tearing energy data at different temperatures was shifted to construct a master curve and an estimate on the threshold value of tear energy was obtained which may be helpful in designing components where material tear is of concern. Strain energy release rate (SERR) value, calculated using the J-integral approach for a pre-existing crack in ABAQUS, was used to estimate the crack growth rate in a given seal cross-section to predict lifetime. In order to assess the viscoelastic behavior and to investigate the long term stress relaxation behavior of the seal material, compression stress relaxation (CSR) tests were performed on molded seals, called as SMORS, over a range of environmental conditions using a custom-designed fixture. The effect of temperature and environment was evident on material property changes and presented in terms of momentary properties and stress relaxation behavior. Various mechanisms involved in material degradation, chain scission and crosslinking, were suggested and insights were gained into how cure state and level of antidegradants in a material dictate the material behavior during the first phase of environmental exposure leading to change in material properties. Ring samples made of silicone were also tested using the fixture to obtain insight additional into material degradation due to aging. Results presented from testing on SMORS showed a lot more variation in data as compared to neat silicone rings due to the complexity involved in making SMORS. For understanding the deformation behavior of an elastomeric seal and its sealing performance, finite element characterization of seal cross-section was carried out on O-ring and SMORS cross-section. The effect of a seal's layout on distribution and magnitude of contact stresses and contact width was investigated for the O-ring and the information obtained thereby helped to analyze a complex assembly such as SMORS, where several interfaces and boundary conditions are involved. Stress/strain profiles were generated to visualize their concentration and distribution in the seal cross-section. Frictionless and rough interfacial conditions between seal material and platens were assumed and it was found that its effect on contact width and peak contact pressure was insignificant. Results obtained from FEA on SMORS were validated through comparison with contact mechanics approach and experimental data and it was found that Lindley's equation correlates well with experimental data whereas ABAQUS overestimates the load values at a given compression. Lindley's approach may be used to develop contact pressure profiles that may help estimate peak contact pressure at a given time so leaking can be avoided. / Ph. D.

Mechanical vibrations are inevitable in drilling operations. Torsional stick-slip vibrations are vibrations that occur in drilling columns, which are produced by periodic variations of torque and characterized by large fluctuations in the speed of the drill bit. These vibrations are dangerous, primarily by the cyclical characteristic of the phenomenon that by the amplitude of the same, which can cause fatigue of the pipe, failures in the components of the drill string, deformations in the walls of the well, excessive wear of the drill, low rate of penetration, and collapse of the drilling process. The frequency of these unwanted oscillations can be reduced by the application of automatic control techniques. The objective of this study is to evaluate through numerical simulations, the application of conventional control techniques, such as proportional-integral control (PI), and nonlinear, as the sliding mode control (SMC) and the input-output linearization control (IOLC), to eliminate the presence of stick-slip oscillation in drilling columns. The controllers are designed primarily to maintain a constant speed of rotation system, by manipulating engine torque, thereby inferentially control the speed of the drill, thus providing optimum operation conditions, beyond preserving system stability. Results of simulations using drill string torsional models of two degrees of freedom (2-DOF) and four degrees of freedom (4-DOF) show the performance of the proposed control systems, which are analyzed and qualitatively compared. / Vibrações mecânicas são inevitáveis nas operações de perfuração. Vibrações torcionais stick-slip são vibrações que ocorrem em colunas de perfuração, as quais são produzidas pelas variações periódicas de torque e caracterizadas por grandes oscilações da velocidade da broca. Estas vibrações são prejudiciais, mais pela característica cíclica do fenômeno que pela amplitude da mesma, podendo originar fadiga da tubulação, falhas nos componentes da coluna de perfuração, deformações nas paredes do poço, desgaste excessivo da broca, baixa taxa de penetração e, inclusive, colapso do processo de perfuração. A frequência destas oscilações indesejadas pode ser reduzida pela aplicação de técnicas de controle automático. O objetivo deste trabalho é avaliar, mediante simulações numéricas, a aplicação de técnicas de controle convencional, como o controle proporcional-integral (PI), e não linear, como o controle por modos deslizantes (SMC) e o controle por linearização entrada-saída (IOLC) para eliminar a presença de oscilações stick-slip em colunas de perfuração. Os controladores são desenvolvidos principalmente para manter constante a velocidade do sistema de rotação, mediante a manipulação do torque do motor, para assim controlar inferencialmente a velocidade da broca, fornecendo desta maneira condições ótimas de operação, além de preservar a estabilidade do sistema. Resultados das simulações, usando modelos torcionais de uma coluna de perfuração de dois graus de liberdade (2-DOF) e de quatro graus de liberdade (4-DOF), mostram o desempenho dos sistemas de controle propostos, os quais são analisados e comparados qualitativamente.

[pt] Um dos objetivos desta tese é analisar o custo computacional do método de Monte Carlo aplicado a um problema modelo de dinâmica, considerando incertezas na força de atrito. O sistema mecânico a ser estudado é composto por um oscilador de um grau de liberdade que se desloca sobre uma esteira em movimento. Considera-se a existência de atrito seco entre a massa do oscilador e a esteira. Devido a uma descontinuidade na força de atrito, a dinâmica resultante pode ser dividida em duas fases que se alternam, chamadas de stick e slip. Neste estudo, um parâmetro da força de atrito dinâmica é modelado como uma variável aleatória. A propagação de incerteza é estudada por meio da aplicação do método de Monte Carlo, considerando três abordagens diferentes para calcular aproximações da resposta dos problemas de valor inicial que modelam a dinâmica do problema: NV) aproximações numéricas calculadas usando método de Runge-Kutta de quarta e quinta ordens com passo de integração variável; NF) aproximações numéricas calculadas usando método de Runge-Kutta de quarta ordem com passo de integração fixo; AN) aproximação analítica obtida com o método de múltiplas escalas. Nas abordagens NV e NF, para cada valor de parâmetro, uma aproximação numérica foi calculada. Já para a AN, apenas uma aproximação analítica foi calculada e avaliada para os diferentes valores usados. Entre as variáveis aleatórias de interesse associadas ao custo computacional do método de Monte Carlo, encontram-se o tempo de execução e o espaço em disco consumido. Devido à propagação de incertezas, a resposta do sistema é um processo estocástico com uma sequência aleatória de fases de stick e slip. Essa sequência pode ser caracterizada pelas seguintes variáveis aleatórias: instantes de transição entre as fases de stick e slip, suas durações e o número de fases. Para estudar as variáveis associadas ao custo computacional e ao processo estocástico foram construídos modelos estatísticos, histogramas normalizados e gráficos de dispersão. O objetivo é estudar a dependência entre as variáveis do processo estocástico e o custo computacional. Porém, a construção destas análises não é simples devido à dimensão do problema e à impossibilidade de visualização das distribuições conjuntas de vetores aleatórios de três ou mais dimensões. / [en] One of the objectives of this thesis is to analyze the computational cost of the Monte Carlo method applied to a toy problem concerning the dynamics of a mechanical system with uncertainties in the friction force. The system is composed by an oscillator placed over a moving belt. The existence of dry friction between the two elements in contact is considered. Due to a discontinuity in the frictional force, the resulting dynamics can be divided into two alternating phases, called stick and slip. In this study, a parameter of the dynamic friction force is modeled as a random variable. Uncertainty propagation is analyzed by applying the Monte Carlo method, considering three different strategies to compute approximations to the initial value problems that model the system s dynamics: NV) numerical approximations computed with the Runge-Kutta method of 4th and 5th orders, with variable integration time-step; NF) numerical approximations computed with the Runge-Kutta method of 4th order, with a fixed integration time-step; AN) analytical approximation obtained with the multiple scale method. In the NV and NF strategies, for each parameter value, a numerical approximation was calculated, whereas for the AN strategy, only one analytical approximation was calculated and evaluated for the different values of parameters considered. The run-time and the storage are among the random variables of interest associated with the computational cost of the Monte Carlo method. Due to uncertainty propagation, the system response is a stochastic process given by a random sequence of stick and slip phases. This sequence can be characterized by the following random variables: the transition instants between the stick and slip phases, their durations and the number of phases. To study the random processes and the variables related to the computational costs, statistical models, normalized histograms and scatterplots were built. Afterwards, a joint analysis was performed to study the dependece between the variables of the random process and the computational cost. However, the construction of these analyses is not a simple task due to the impossibility of viewing the distributionto of joint distributions of random vectors of three or more.