Lors d'une opération d'usinage, l'intégrité de la surface usinée et l’optimisation du procédé sont conditionnées par les paramètres de coupe (vitesses de coupe et d’avance, géométrie et matériau de l'outil...). Certaines conditions de coupe peuvent induire des effets indésirables tels que des vibrations importantes, des efforts de coupe excessifs et une usure prématurée de l'outil, conduisant à des qualités de surfaces médiocres. Dans l’industrie, l’utilisation d'approches empiriques pour opérer ce choix se révèle couteux et difficilement exploitable. Le développement d’outils de simulation basés sur des modèles prédictifs s’avère nécessaire. Ces modèles permettent de maitriser et de comprendre les phénomènes thermomécaniques aux interfaces outil-copeau et outil-pièce qui conditionnent l'intégrité de la surface usinée ainsi que la durée de vie de l'outil de coupe. L'objectif de la thèse est la modélisation des effets thermomécaniques en usinage avec des approches hybrides 'Analytique-c. Ceci permet d'analyser l'interaction entre les conditions de coupe du procédé d'usinage, le comportement du matériau et les conditions tribologiques aux interfaces outil-copeau et outil-pièce (contact collant-glissant, partage de la source de chaleur due au frottement). Le travail comporte également une validation expérimentale pour la coupe orthogonale à sec. Le modèle proposé est basé sur les développements suivants: (i) mise en place d'une 'Approche 1D par Tranche' pour la prise en compte de l'écoulement de la matière dans la zone primaire de cisaillement du copeau, (ii) modélisation du problème thermique transitoire non linéaire dans le système 'copeau-outil-pièce' en couplant une formulation EF de type Petrov-Galerkin avec la méthode de Newton-Raphson et une intégration implicite dans le temps, (iii) une nouvelle formulation de la distribution de pression le long de la face de coupe de l'outil, (iv) une nouvelle approche pour gérer le partage de la source de chaleur par friction à l'interface outil-copeau. La démarche proposée permet de mettre en place une modélisation thermomécanique de l'interaction outil-matière applicable aux procédés industriels comme le perçage très utilisé dans le domaine aéronautique. Comparée aux simulations basées sur la méthode des Éléments Finis, l'approche développée requiert un temps de calcul de l'ordre de quelques minutes avec une précision comparable / In dry machining, the thermomechanical process of chip formation, the tool wear and the surface integrity depend strongly on the tribological conditions along the tool rake face. Besides, the friction conditions at the tool-chip interface and along the round cutting edge are very complex. It should be noted that to understand the friction effects in machining, we have to analyse the inherent relationship among, the cutting conditions (cutting and feed velocities, tool geometry), the workpiece material behaviour, the thermomechanical characteristics of the tool material, the frictional heat partition in the sliding zone and the friction conditions at the tool-chip and tool-workpiece interfaces. Due to the problem complexity, it appears that despite many works on machining, the understanding of the effect of friction conditions requires further investigations. In the present work, to identify the interaction between the thermomechanical phenomena at the tool-chip interface and the material flow in the primary shear zone; an analytical model has been coupled with a finite element approach. For the tool rake face, a new pressure model was developed. The transient nonlinear thermal problem in the workpiece-tool-chip system has been solved by using a FE model based on the Petrov-Galerkin formulation. The coupling between the primary shear zone (PSZ) (chip formation), the secondary shear zone (SSZ) (sticking zone) and the frictional heat at the sliding zone has been taking into account. The model allows to determine in a fast and simple way several significant machining parameters as: (i) the cutting forces, (ii) the temperature distribution in the tool-chip-workpiece system, (iii) the heat flux from the PSZ to the workpiece, (iv) the tool-chip contact length, (v) the frictional heat partition and (vi) the apparent friction coefficient. The proposed model allows to analyze different industrial machining processes such as drilling and milling

The thesis is concerned with the contact mechanics behaviour of non-conforming solids. The geometry of the solids considered gives rise to various contact configurations, from concentrated contacts with circular and elliptical configuration to those of finite line nature, as well as those of less concentrated form such as circular flat punches. The radii of curvature of mating bodies in contact or impact give rise to these various nonconforming contact configurations and affect their contact characteristics, from those considered as semi-infinite solids in accord with the classical Hertzian theory to those that deviate from it. Furthermore, layered solids have been considered, some with higher elastic modulus than that of the substrate material (such as hard protective coatings) and some with low elastic moduli, often employed as tribological coatings (such as solid lubricants). Other bonded layered solids behave in viscoelastic manner, with creep relaxation behaviour under load, and are often used to dampen structural vibration upon impact. Analytic models have been developed for all these solids to predict their contact and impact behaviour and obtain pressure distribution, footprint shape and deformation under both elastostatic and transient dynamic conditions. Only few solutions for thin bonded layered elastic solids have been reported for elastostatic analysis. The analytical model developed in this thesis is in accord with those reported in the literature and is extended to the case of impact of balls, and employed for a number of practical applications. The elastostatic impact of a roller against a semi-infinite elastic half-space is also treated by analytic means, which has not been reported in literature. Two and three-dimensional finite element models have been developed and compared with all the derived analytic methods, and good agreement found in all cases. The finite element approach used has been made into a generic tool for all the contact configurations, elastic and viscoelastic. The physics of the contact mechanical problems is fully explained by analytic, numerical and supporting experimentation and agreement found between all these approaches to a high level of conformance. This level of agreement, the development of various analytical impact models for layered solids and finite line configuration, and the development of a multi-layered viscoelastic transducer with agreed numerical predictions account for the main contributions to knowledge. There are a significant number of findings within the thesis, but the major findings relate to the protective nature of hard coatings and high modulus bonded layered solids, and the verified viscoelastic behaviour of low elastic modulus compressible thin bonded layers. Most importantly, the thesis has created a rational framework for contact/impact of solids of low contact contiguity.

During life of railway vehicles, shape irregularities develop on wheels and rails because of wear. The shape irregularities then affect forces in wheel-rail contact and cause further damage of contact surfaces, vibrations and noise and increase risk of derailment. A numerical simulation of railway vehicle motion with more details on contact surfaces geometry was created to investigate dynamic contact loads in wheel-rail contact. A variety of methods can be used to evaluate forces in rolling contact, the method chosen for this study was algorithm CONTACT based on boundary element method. Four studies are presented in this papers: contact loads from a wheel with a flat and with a wavy tread pattern, loads on wavy rail and load in a curve. The first three studies investigated effects of existing wear patterns, the last one looked for cause of common wear pattern developing on rails. Results of the studies with worn components used showed that the worst kind of shape irregularities is a flat present on wheel. This type of shape cause loss of contact and following impacts. The study of ride in curve showed that cause of high wear in curves, especially those with small radii, is caused by vibration of wheelset. This vibration is then caused by different length of inner and outer rail and wheels travelling along a different path.

Cette thèse propose une approche pour générer un mouvement corps complet avec contacts non coplanaires, permettant à un robot de se déplacer dans un environnement, de manipuler des objets complexes ou de collaborer avec différents agents. Les méthodes développées dans cette thèse tentent de prendre en compte une grande variété de robots, de l'humanoïde au manipulateur à base fixe en passant par les objets sous actionnés. En premier lieu, nous abordons le problème du choix des positions des points de contacts qu'un robot sous-actionné doit prendre pour se déplacer dans l'environnement. Nous calculons, en un seul problème d'optimisation non-linéaire, une séquence de postures qui satisfait une séquence de contacts donnés. Cette formulation permet de trouver la position des contacts optimale, car le choix de la position d'un contact d'une posture va prendre en compte les postures précédentes et suivantes. Elle permet aussi d'effectuer des tâches pour certaines postures qui prendront en compte l'aspect prioritaire du déplacement. Nous introduisons ensuite une méthode de génération de mouvement qui, en se basant sur la programmation quadratique, permet de résoudre le problème de géométrie inverse et de la dynamique inverse pour un robot à base fixe ou mobile, tout en satisfaisant des contraintes d'égalités et d'inégalités.Cette génération de mouvement est assez rapide pour fonctionner à la vitesse de la boucle de contrôle des robotsHRP2-10 et HRP4, et peut donc être utilisé en temps réel. À l'aide d'une machine à état, nous transformons la séquence de postures calculée à priori en une série de tâches à effectuer par le générateur de mouvement, ce qui permet à notre robot de se déplacer dans un environnement complexe. Nous étendons alors notre méthode de génération de mouvement pour calculer la commande d'un nombre arbitraire de robots. Cette extension nous permet de gérer des tâches de manipulation d'objets complexes, de collaboration entre plusieurs agents et de mouvement dans un environnement dynamique. Nous pouvons aussi spécifier directement les tâches dans le repère de l'objet manipulé pour faciliter l'élaboration de notre consigne. Dans l'optique de valider cette méthode sur un robot réel, nous formulons le problème d'estimation des paramètres inertiels d'un objet manipulé grâce à l'algèbre vectorielle spatiale. Finalement, nous validons nos travaux sur les robots HRP2-10 et HRP4. Sur le premier robot, nous validons la génération de posture et la génération de mouvement mono-robot sur le scénario demonté d'une échelle verticale aux normes industrielles. La manipulation d'objets et l'estimation des paramètres inertiels sont validées par la suite sur le robot HRP4. / This PhD proposes a whole body motion generation approach with non coplanar contacts that allowsa robot to move in an environment, manipulate complex objects or collaborate with differentagents.Methods developed in this PhD try to manage many kinds of robots, from the humanoid to thefixed base manipulator and also handling underactuated objects.Firstly, we address the problem contacts positioning that an underactuated robot should taketo move in its environment.We compute in one non-linear optimization problem a sequence of postures that fulfill aninputed contact list. This formulation allows to find the optimal contact placement regardingprevious and next stances. It also allows to execute a task for some posture while taking into accountthe priority of the motion.Next, we introduce a motion generation method that uses quadratic programming to solveinverse kinematics and dynamics problems for a fixed or mobile base robot under equality andinequality constraints.This motion generation is fast enough to fit the HRP2-10 and HRP4 control loop andcan be used in real-time.With a finite state machine we turn the posture sequence into a list of tasks that should beexecuted by the motion generation to allow a robot to move in a complex environment.We extend this motion generation scheme to compute the motion of an arbitrary number of robots.This extension allows us to manage complex object manipulation tasks, multi-agent collaboration andmotion in a dynamic environment. We can also specify a task in the manipulated object frameto ease motion design.To validate this method on a real robot, we formulate inertial parametersestimation of manipulated objects with spatial vector algebra.Finally, we validate our works on the HRP2-10 and HRP4 robot. On the first one,we validate the posture and mono-robot motion generation on a scenario where the robot climbs anindustry standard vertical ladder.On the second one, we validate object manipulation and inertial parameters estimation.

Utilizing survey data collected through the City Planning and Indigeneity on the Prairies (CPIP) Project in the summer of 2014 this thesis examines what, if any, affect contact between Aboriginal and non-Aboriginal people has on non-Aboriginal people’s racial attitudes. Two hypotheses are tested: contact with Aboriginal people reduces new and old-fashioned racism attitudes amongst non-Aboriginal people; and contact has greater effect on the racial attitudes of those with lower education levels than those with higher education levels. Key findings indicate that both hypotheses are in fact correct but each with its own caveats.

Dans le domaine ferroviaire de nombreux incidents ont montré le problème de l'échauffement de la caténaire au droit du contact avec le pantographe durant le captage de courant à l'arrêt. L'échauffement à l'interface pantographe/caténaire peut entraîner la rupture du fil de contact de la caténaire. Le travail présenté dans cette étude, issu de la collaboration entre la SNCF et l'institut FEMTO-ST, vise à améliorer la compréhension des phénomènes physiques mis en jeu. L'objectif à terme est d'anticiper une maintenance coûteuse comme le remplacement du fil de contact. Un modèle électrothermique 2D transitoire du fil de contact a été développé et couplé à un modèle thermique 1D transitoire pour obtenir une modélisation quasi 3D. La modélisation, alimentée par des données expérimentales, permet de déterminer la répartition du courant électrique et donc la production de chaleur interne. Une modélisation électrothermique de la bande vient compléter le système. / In the railroad and trains domain, many incidents show the main problem of overheating of the catenary at the contact with the pantograph when the train was stopped whereas all the electrical systems of the train should nevertheless be fed. Analysis of these incidents has shown that the overheating of the interface catenary-pantograph during the ream conditioning was sometimes at the origin of the break of the contact wire. In order to forecast such very expensive problems for the company, the French National Railway Company (SNCF) and the research institute FEMTO-ST carried out theoretical and experimental studies to better understand this phenomenon. First a quasi 3D transient electrothermal modeling tool has been developed for the contact wire. It has also permitted to estimate the distribution of current in the wire in order to obtain the internal heat power generation. An electrothermal modelisation of the strip complete the system. Finally the heat transfer equation in the wire with particular boundary conditions has been solved in all the finite differences network thanks to the Euler's implicit method.

During life of railway vehicles, shape irregularities develop on wheels and rails because of wear. The shape irregularities then affect forces in wheel-rail contact and cause further damage of contact surfaces, vibrations and noise and increase risk of derailment. A numerical simulation of railway vehicle motion with more details on contact surfaces geometry was created to investigate dynamic contact loads in wheel-rail contact. A variety of methods can be used to evaluate forces in rolling contact, the method chosen for this study was algorithm CONTACT based on boundary element method. Four studies are presented in this papers: contact loads from a wheel with a flat and with a wavy tread pattern, loads on wavy rail and load in a curve. The first three studies investigated effects of existing wear patterns, the last one looked for cause of common wear pattern developing on rails. Results of the studies with worn components used showed that the worst kind of shape irregularities is a flat present on wheel. This type of shape cause loss of contact and following impacts. The study of ride in curve showed that cause of high wear in curves, especially those with small radii, is caused by vibration of wheelset. This vibration is then caused by different length of inner and outer rail and wheels travelling along a different path.

No description available.

No description available.

In this work, we develop an interface-preserving level-set method in the finite-element framework for interfacial flows with moving contact lines. In our method, the contact line is advected naturally by the flow field. Contact angle hysteresis can be easily implemented without explicit calculation of the contact angle or the contact line velocity, and meshindependent results can be obtained following a simple computational strategy. We have implemented the method in three dimensions and provide numerical studies that compare well with analytical solutions to verify our algorithm. We first develop a high-order numerical method for interface-preserving level-set reinitialization. Within the interface cells, the gradient of the level set function is determined by a weighted local projection scheme and the missing additive constant is determined such that the position of the zero level set is preserved. For the non-interface cells, we compute the gradient of the level set function by solving a Hamilton-Jacobi equation as a conservation law system using the discontinuous Galerkin method. This follows the work by Hu and Shu [SIAM J. Sci. Comput. 21 (1999) 660-690]. The missing constant for these cells is recovered using the continuity of the level set function while taking into account the characteristics. To treat highly distorted initial conditions, we develop a hybrid numerical flux that combines the Lax-Friedrichs flux and a penalty flux. Our method is accurate for non-trivial test cases and handles singularities away from the interface very well. When derivative singularities are present on the interface, a second-derivative limiter is designed to suppress the oscillations. At least (N + 1)th order accuracy in the interface cells and Nth order accuracy in the whole domain are observed for smooth solutions when Nth degree polynomials are used. Two dimensional test cases are presented to demonstrate superior properties such as accuracy, long-term stability, interface-preserving capability, and easy treatment of contact lines. We then develop a level-set method in the finite-element framework. The contact line singularity is removed by the slip boundary condition proposed by Ren and E [Phys. Fluids, vol. 19, p. 022101, 2007], which has two friction coefficients: βN that controls the slip between the bulk fluids and the solid wall and βCL that controls the deviation of the microscopic dynamic contact angle from the static one. The predicted contact line dynamics from our method matches the Cox theory very well. We further find that the same slip length in the Cox theory can be reproduced by different combinations of (βN; βCL). This combination leads to a computational strategy for mesh-independent results that can match the experiments. There is no need to impose the contact angle condition geometrically, and the dynamic contact angle automatically emerges as part of the numerical solution. With a little modification, our method can also be used to compute contact angle hysteresis, where the tendency of contact line motion is readily available from the level-set function. Different test cases, including code validation and mesh-convergence study, are provided to demonstrate the efficiency and capability of our method. Lastly, we extend our method to three-dimensional simulations, where an extension equation is solved on the wall boundary to obtain the boundary condition for level-set reinitializaiton with contact lines. Reinitialization of ellipsoidal interfaces is presented to show the accuracy and stability of our method. In addition, simulations of a drop on an inclined wall are presented that are in agreement with theoretical results. / Doctor of Philosophy / When a liquid droplet is sliding along a solid surface, a moving contact line is formed at the intersection of the three phases: liquid, air and solid. This work develops a numerical method to study problems with moving contact lines. The partial differential equations describing the problem are solved by finite element methods. Our numerical method is validated against experiments and theories. Furthermore, we have implemented our method in three-dimensional problems.