Lors de l'opération d’usinage à grande vitesse, la résistance à l'usure des outils de coupe est améliorée par l’utilisation des revêtements mono ou multicouches sur les faces actives de l’outil. Cependant, le chargement thermomécanique généré à l'interface outil-pièce affecte considérablement les zones de contact. Par cet effet, plusieurs modes d'usure tels que la fissuration, l’abrasion, l’adhésion et le délaminage du revêtement peuvent se manifester. L'étude du comportement des revêtements et de leurs différents modes de dégradation permet de mieux comprendre leur impact sur la durée de vie de l'outil et ainsi optimiser le procédé d'usinage. Dans ce travail de thèse, une approche numérique multi-étapes a été proposée pour prédire l'usure des outils de coupe revêtus. Cette approche est composée par trois principales étapes. La première consiste à effectuer une simulation éléments finis de l’usinage pour une courte durée (jusqu’à la stabilisation du chargement à l’interface outil/pièce). La deuxième étape consiste à récupérer ce chargement et de l’utiliser comme une entrée du modèle XFEM/Level-set. Ce dernier permet d’analyser le comportement des couches de revêtement sans recours à un maillage conforme aux interfaces. Par conséquence, la distorsion du maillage est évitée lorsque le profil d'outil usé est mis à jour, ainsi que le temps de calcul CPU est drastiquement réduit. La dernière étape de cette approche consiste à calculer le taux d’usure et ainsi prédire le déplacement des nœuds de l’outil de coupe affectés par l’usure. Les essais expérimentaux ont permis d’une part d’identifier les paramètres de contact outil/pièce, et d’autre part de valider l’approche proposée / In high speed machining, wear resistance of the cutting tools is improved by depositing single or multilayered coatings on their surface. However, the thermomechanical loading generated at the tool-workpiece interface greatly affects the contact zones. For this purpose, several wear modes such as cracking, abrasion, adhesion and delamination of the coating can be occurred. The study of the coatings behavior and their different degradation modes lead to better understanding of their impact on the tool life and machining process under optimal conditions. In this PhD thesis work, a multi-step numerical approach has been proposed to predict wear of the coated cutting tools. This approach involves three main steps. The first is to perform a finite element simulation of the orthogonal cutting for a short time (until the loading stabilization at the tool/workpiece interface). The second step is to recover this loading and use it as an input for the XFEM/Level-set model. The latter allow to take into account the coating layers presence without any need of mesh conforming to the interfaces. As a result, the mesh distortion is avoided when the worn tool profile is updated, as well as the CPU calculation time is drastically reduced. The final step of this approach is to convert the wear rate equation into a nodal displacement, thus representing the cutting tool wear. Based on the experimental tests, a procedure for identifying tool/workpiece contact parameters, and for calibrating the wear equation for each coating layer has been proposed. Experimental trials have been also used to validate the proposed approach

Orientador: Prof. Dr. Wesley Góis / Dissertação (mestrado) - Universidade Federal do ABC, Programa de Pós-Graduação em Engenharia Mecânica, 2015. / Na área de otimização de estruturas, um método muito explorado é o Método de Otimização Topológica (MOT), que combina o Método dos Elementos Finitos (MEF) com um algoritmo de otimização para encontrar a distribuição ótima de material no interior do domínio de projeto da estrutura. Apesar de o MOT ser um método altamente eficiente para a determinação do layout ótimo da estrutura, o mesmo sofre com problemas de instabilidade numérica, sendo caracterizados por soluções dependentes da malha aplicada (dependência de malha), problemas de mínimo local e pelo surgimento de regiões contendo o padrão de tabuleiro de xadrez (ou checkerboard pattern, em inglês) na topologia final da estrutura. O problema de instabilidade de tabuleiro de xadrez é descrito por regiões no domínio da estrutura contendo elementos com presença e ausência de material, distribuídos de forma intercalada e cujo formato final se assemelha ao de um tabuleiro de xadrez. Esta distribuição produz, numericamente, uma rigidez maior do que uma configuração de mesmo volume de material distribuída uniformemente (o aumento na rigidez é comumente atribuído a uma "pseudo-rigidez" incumbida a elementos quadrilaterais (e triangulares) de baixa ordem (DIAZ e SIGMUND, 1995; JOG e HABER, 1996)). Este fenômeno é causado pelo mau condicionamento do conjunto de soluções das equações de equilíbrio obtidas via MEF em sua formulação clássica, e é convencionalmente contornado por meio da aplicação de filtros de sensibilidade, solução esta de fácil implementação computacional e que garante o controle eficaz da formação do tabuleiro de xadrez. Nesse contexto, este trabalho apresenta uma variante do método SIMP (Solid Isotropic Material with Penalization), desenvolvida por meio da aplicação do Método dos Elementos Finitos Generalizados (MEFG) ¿ metodologia alternativa, que mais recentemente vem sendo explorada, como solução para as limitações da formulação clássica do Método dos Elementos Finitos (MEF) ¿ de modo a estudar o impacto da mesma numa possível redução dos problemas de instabilidade numérica previamente citada. O algoritmo desenvolvido ao longo do trabalho é baseado no algoritmo apresentado em SIGMUND (2001), assim como alguns exemplos utilizados como referência. Por fim, exemplos numéricos apresentados ao longo do trabalho comprovam a redução parcial, ou em certos casos, completa do problema de instabilidade de tabuleiro de xadrez, demonstrando a viabilidade da utilização do MEFG como uma alternativa para obtenção da distribuição ótima de material no interior do domínio de projeto da estrutura. / Within the structures optimization study area, one of the extensively explored methods is the Topology Optimization Method (TOM), which combines the Finite Element Method (FEM) with an optimization algorithm to find the optimal distribution of material inside the structure design domain. Although the TOM is a highly efficient method for determining the optimal structural layout, it suffers from numerical instability problems, characterized by the mesh dependency, local minima and the checkerboard pattern. The latter, characterized by the alternating density values obtained from adjacent elements, while the patch of element maintains the connectivity resembling that of a checkerboard; which produces a numerically greater stiffness than a uniform distribution configuration (the increase in stiffness is commonly attributed to a "pseudo-stiffness" entrusted to low-order quadrilateral (and triangular) elements (DIAZ and Sigmund 1995; JOG and HABER, 1996)). This phenomenon occurs due to the finite element formulation used in the optimization process, and can be conventionally avoided through the use of sensitivity filter, a solution of easy implementation which ensures an effective control of the "checkerboard pattern". Thus, the main objective of this study is to present a variant of the SIMP method (Solid Isotropic Material with Penalization), developed by the application of the Generalized Finite Element Method (GFEM) - an alternative methodology of the Finite Element Method (FEM), which most recently has been explored as a solution to the limitations of the FEM classic formulation - on the Topology Optimization Method (TOM), in order to study his impact on the reduction of the numerical instability problems previously cited. The algorithm developed throughout this work is based on the algorithm shown in Sigmund (2001), and has some of his examples used as reference. At last, the numerical examples provided throughout this work demonstrate a partial reduction, or in some cases, a complete removal of checkerboard pattern, showing the use of GFEM as a possible alternative for obtaining the optimal material distribution.

No description available.

Cette thèse contient deux matières différentes. Dans la première partie, deux cas sont considérés. L'un est le modèle plus lisse de la plaque mince et l'autre est les équations des limites elliptiques avec des données limites incertaines. Dans cette partie, les convergences stochastiques des méthodes des éléments finis sont prouvées pour chaque problème.Dans la deuxième partie, nous fournissons une analyse mathématique du problème inverse linéarisé dans la tomographie d'impédance électrique multifréquence. Nous présentons un cadre mathématique et numérique pour une procédure d'imagerie du tenseur de conductivité électrique anisotrope en utilisant une nouvelle technique appelée Tentomètre de diffusion Magnéto-acoustographie et proposons une approche de contrôle optimale pour reconstruire le facteur de propriété intrinsèque reliant le tenseur de diffusion au tenseur de conductivité électrique anisotrope. Nous démontrons la convergence et la stabilité du type Lipschitz de l'algorithme et présente des exemples numériques pour illustrer sa précision. Le modèle cellulaire pour Electropermécanisme est démontré. Nous étudions les paramètres efficaces dans un modèle d'homogénéisation. Nous démontrons numériquement la sensibilité de ces paramètres efficaces aux paramètres microscopiques critiques régissant l'électropermécanisme. / This thesis contains two different subjects. In first part, two cases are considered. One is the thin plate spline smoother model and the other one is the elliptic boundary equations with uncertain boundary data. In this part, stochastic convergences of the finite element methods are proved for each problem.In second part, we provide a mathematical analysis of the linearized inverse problem in multifrequency electrical impedance tomography. We present a mathematical and numerical framework for a procedure of imaging anisotropic electrical conductivity tensor using a novel technique called Diffusion Tensor Magneto-acoustography and propose an optimal control approach for reconstructing the cross-property factor relating the diffusion tensor to the anisotropic electrical conductivity tensor. We prove convergence and Lipschitz type stability of the algorithm and present numerical examples to illustrate its accuracy. The cell model for Electropermeabilization is demonstrated. We study effective parameters in a homogenization model. We demonstrate numerically the sensitivity of these effective parameters to critical microscopic parameters governing electropermeabilization..

The thesis deals with frequency domain finite elements and time domain finite elements. Derivation of the wave equation, generation of the discretization mesh and the composition of matrices for solving wave equation are presented. Exploitation of COMSOL Multiphysics for the analysis of the electromagnetic wave propagation in a wave guide is explained. Finally, in this thesis is described the m-file of MATLAB for analyzing electric field intensity of the waveguide in longitudinal direction to both domain. In last part is concisely exposed perfect match layers.

In this paper a 2D coupled thermal-stress finite element model is established and used to predict thermal phenomena at the disc-pad interface of a disc brake system. The importance of certain critical settings and parameters for the 2D FE model has been identified (such as, a limited degree of freedom for a brake pad in place of accepted practice that considers uniform contact), here a non-uniform pressure distribution resulting from friction bending moment effects due to the introduction of a pivot point. These parameters affect the distributions of both interface temperature and pressure. The simulation results show that when the interface conductance h is 10^6 W/m^2K or higher, the interface temperature distribution is no longer sensitive to friction bending moment effects. However, when h is 30000 W/m^2K or lower, the interface temperature distribution and heat partition ratio are significantly affected by the setting used for the rotational degree of freedom of the pad. The simulation results provide a useful reference for a better design of a disc brake system for different applications.

The Explicit Finite Element Method is a powerful tool in nonlinear dynamic finite element analysis. Recent major developments in computational devices, in particular, General Purpose Graphical Processing Units (GPGPU's) now make it possible to increase the performance of the explicit FEM. This dissertation investigates existing explicit finite element method algorithms which are then redesigned for GPU's and implemented. The performance of these algorithms is assessed and a new asynchronous variational integrator spatial decomposition (AVISD) algorithm is developed which is flexible and encompasses all other methods and can be tuned based for a user-defined problem and the performance of the user's computer. The mesh-aware performance of the proposed explicit finite element algorithm is studied and verified by implementation. The current research also introduces the use of a Particle Swarm Optimization method to tune the performance of the proposed algorithm automatically given a finite element mesh and the performance characteristics of a user's computer. For this purpose, a time performance model is developed which depends on the finite element mesh and the machine performance. This time performance model is then used as an objective function to minimize the run-time cost. Also, based on the performance model provided in this research and predictions about the changes in GPU's in the near future, the performance of the AVISD method is predicted for future machines. Finally, suggestions and insights based on these results are proposed to help facilitate future explicit FEM development.

Doctor of Philosophy / Department of Mechanical and Nuclear Engineering / Liang-Wu Cai / X. J. Xin / Thin walled structures such as plates and shells have application in many fields of engineering because these structures are light weight and can support large loads when designed suitably. In real world, loads may cause these structures to vibrate which can be undesirable causing fatigue and failure of the structure. Such undesirable vibrations need to be reduced or eliminated. In this work, analytical studies of flexural wave propagation for idealized geometries are conducted and finite element method (FEM) is used to explore the effects of composite panel designs of finite size for the reduction of vibration transmission. In the analytical studies, the influence of the material properties on the reflection and transmission characteristics are explored for an infinite bi-material plate, and infinite plate with a strip inhomogeneity. In the analytical study of an infinite thin plate with a solid circular inclusion, the far and near field scattering characteristics are explored for different frequencies and material properties. All the analytical studies presented here and reported in the literature consider infinite plates to characterize the flexural wave propagation. Obtaining closed form solutions to characterize the flexural wave propagation in a finite plate with inclusions is mathematically difficult process. So, FEM is used to explore the composite panel designs. The understanding gained about the material properties influence on the flexural wave propagation from analytical studies helped with the choice of materials for FEM simulations. The concept of phononic crystals is applied to define the design variations that are effective in suppressing vibration transmission. Various design configurations are explored to study the effects of various parameters like scatterer’s material properties, geometry and spatial pattern. Based on the knowledge gained through a systematic parametric study, a final design of the composite sandwich panel is proposed with an optimum set of parameters to achieve the best vibration reduction. This is the first study focused on reducing vibration and wave transmission in composite rotorcraft fuselage panels incorporating the concept of phononic crystals. The optimum sandwich panel design achieved 98% vibration transmission reduction at the frequency of interest of 3000 Hz.

Call number: LD2668 .T4 1978 L44 / Master of Science

Aerosol particle deposition rates on surfaces inside a clean room are predicted by a model developed to account for particle convection, diffusion and sedimentation. External forces acting on the particle also influence the rate of deposition. Both electrical charge build up on product surfaces and temperature gradients in the air near the product surface are known to effect the rate of deposition. A description of an electrostatic and thermophoretic force on the particle is thus included in the model. The equations governing the particle deposition process and the approach used in obtaining a solution to these equations are both described. A finite element numerical solution is detailed, followed by a description of the electrostatic force models. Finally, predictions of the model are presented with a comparison to data experimentally obtained by other researchers.