Les problèmes de contrôle optimal se rencontrent dans l'industrie aérospatiale et dans la mécanique. Leur étude conduit à des difficultés mathématiques importantes. Dans cette thèse, nous nous intéressons aux conditions d'optimalité pour certains problèmes de contrôle avec des contraintes exprimées en termes d'inclusions différentielles. Nous considérons aussi des problèmes de contrôle associés aux modèles mathématiques issus de la Mécanique du Contact. Cette thèse est structurée en deux parties et six chapitres. La première partie, contenant les Chapitres 1, 2 et 3, représente un résumé de nos résultats, en Français. Nous y présentons les problèmes étudiés, les hypothèses sur les données, les notations utilisées ainsi que l’énoncé des principaux résultats. Les démonstrations sont omises. La deuxième partie du manuscrit représente la partie principale de la thèse. Elle contient les Chapitres 4, 5 and 6, chacun ayant fait l'objet d'une publication (parue ou soumise) dans une revue internationale avec comité de lecture.Nous y présentons nos principaux résultats, accompagnés des démonstrations et des références bibliographiques. / Optimal control problems arise in aerospace industry and in mechanics. They are challenging and involve important mathematical difficulties. In this thesis, we are interested to derive optimality conditions for optimal control problems with constraints under the form of differential inclusions. We also consider optimal control problems in the study of some boundary value problems arising in Contact Mechanics. The thesis is structured in two parts and six chapters. Part I represents an abstract of the main results, in French. It contains Chapters 1, 2 and 3. Here we present the problems we study together with the assumptions on the data, the notation and the statement of the main results. The proofs of these results are omitted, since them are presented in Part II of the manuscript.Part II represents the main part of the thesis. It contains Chapters 4, 5 and 6. Each of these chapters made the object of a paper published (or submitted) in an international journal. Here we present our main results, together with the corresponding proofs and bibliographical references.

The Generalized Finite Element Method (GFEM) is a numerical method based on the Partition of Unity (PU) concept and inspired on both the Partition of Unity Method (PUM) and the hp-Cloud method. According to the GFEM, the PU is provided by first-degree Lagragian interpolation functions, defined over a mesh of elements similar to the Finite Element Method (FEM) meshes. In fact, the GFEM can be considered an extension of the FEM to which enrichment functions can be applied in specific regions of the problem domain to improve the solution. This technique has been successfully employed to solve problems presenting discontinuities and singularities, like those that arise in Fracture Mechanics. However, most publications on the method are related to linear analyses. The present thesis is a contribution to the few studies of nonlinear analyses of Solid Mechanics by means of the GFEM. One of its main topics is the derivation of a segment-to-segment generalized contact element based on the mortar method. Material and kinematic nonlinear phenomena are also considered in the numerical models. An Object-Oriented design was developed for the implementation of a GFEM nonlinear analyses framework written in Python programming language. The results validated the formulation and demonstrate the gains and possible drawbacks observed for the GFEM nonlinear approach. / O Método dos Elementos Finitos Generalizados (MEFG) é um método numérico baseado no conceito de partição da unidade (PU) e inspirado no Método da Partição da Unidade (MPU) e o método das Nuvens-hp. De acordo com o MEFG, a PU é obtida por meio de funções de interpolação Lagragianas de primeiro grau, definidas sobre uma rede de elementos similar àquela do Método dos Elementos Finitos (MEF). De fato, o MEFG pode ser considerado uma extensão do MEF para a qual se pode aplicar enriquecimentos em regiões específicas do domínio, buscando melhorias na solução. Esta técnica já foi aplicada com sucesso em problemas com descontinuidades e singularidades, como os originários da Mecânica da Fratura. Apesar disso, a maioria das publicações sobre o método está relacionada a análises lineares. A presente tese é uma contribuição aos poucos estudos relacionados a análises não-lineares de Mecânica dos Sólidos por meio do MEFG. Um de seus principais tópicos é o desenvolvimento de um elemento de contato generalizado do tipo segmento a segmento baseado no método mortar. Fenômenos não lineares devidos ao material e à cinemática também são considerados nos modelos numéricos. Um projeto de orientação a objetos para a implementação de uma plataforma de análises não-lineares foi desenvolvido, escrito em linguagem de programação Python. Os resultados validam a formulação e demonstram os ganhos e possíveis desvantagens da abordagem a problemas não lineares por meio do MEFG.

Cette thèse est consacrée à l’étude de méthodes de domaines fictifs pour les éléments finis. Ces méthodes, initialement conçues pour l’approximation de problèmes d’interactions fluide/structure, consistent à prolonger un domaine réel par un domaine de géométrie simple appelé domaine fictif. On applique ces méthodes à un problème de contact unilatéral sans frottement en petite déformation entre deux corps élastiques séparés par une distance initiale non nulle et possédant par ailleurs des conditions aux bords de type Dirichlet et Neumann. Les deux premiers chapitres sont consacrés à l’introduction des méthodes de domaines fictifs et du problème unilatéral de contact de deux corps élastiques. Le chapitre 3 est consacré à l’analyse a priori et à l’étude numérique de ce problème de contact en domaine fictif avec les conditions aux bords de Dirichlet et de contact qui sont prises en compte à l’aide d’une méthode de type Nitsche. Des résultats théoriques de consistance de la méthode discrète, d’existence et d’unicité sont présentés. Afin d’obtenir une estimation d’erreur a priori optimale, une stabilisation de la méthode de domaine fictif est nécessaire. Ces résultats sont validés numériquement sur des cas tests en dimensions deux et trois. Le chapitre 4 est consacré à l’étude d’un estimateur d’erreur de type résidu d’un problème de contact sans domaine fictif entre un corps élastique et un corps rigide. Les résultats théoriques sont également validés sur deux cas tests numériques : un domaine rectangulaire avec seulement une partie de la zone de contact en contact effectif ainsi qu’un contact de type Hertz en dimensions deux et trois. Le chapitre 5 est une généralisation du chapitre 4 à l’approche domaine fictif et au cas de deux corps élastiques. / This thesis is dedicated to the study of the fictitious domain methods for the finite element methods. These methods, initially designed for the fluid-structure interaction, consist in immersing the real domain in a simply-shaped and a geometrically bigger domain called the fictitious domain. We apply these methods to a unilateral frictionless contact problem in small deformation of two deformable elastics bodies separated by an initial gap and satisfying boundary Dirichlet and Neumann conditions. The first two chapters are devoted to the introduction of these methods and to the unilateral contact problem. The chapter 3 is dedicated to a theoretical study for Dirichlet and contact boundary conditions taken into account with a Nitsche type method. Some theoretical results are presented: the consistency of the discrete method, existence and uniqueness results. To obtain an optimal a priori error estimate, a stabilized fictitious domain method is necessary. These results are numerically validated using Hertz contact in two and three dimensions. The chapter 4 is devoted to the study of a residual-based a posteriori error estimator, without the fictitious domain approach, between an elastic body and rigid obstacle. The numerical study of two tests cases will be performed: a rectangular domain with only a part of the potential zone of contact in effective contact as well as a Hertz contact in two and three dimensions. The chapter 5 is a generalization of the chapter 4 to the fictitious domain approach and the care of to two elastics bodies.

Este trabalho descreve um modelo matemático para simulação da laringe humana durante a fonação. O objetivo foi produzir uma técnica computacional de grande escala de processamento para capturar os fenômenos fisiológicos que ocorrem na laringe durante a vocalização e servir de base para estudos mais aprofundados sobre esta importante estrutura do corpo humano. Usando o método dos elementos finitos como base para discretizar as equações dos tecidos musculares da laringe e das equações de Navier-Stokes, e um modelo de descrição da colisão entre as pregas vocais, o sinal glotal foi obtido a partir de diferentes geometrias de laringes com diferentes propriedades viscoelásticas. Os resultados confirmaram a teoria mioelástica-aerodinâmica que descreve a dinâmica da fonação, reproduzindo inclusive fenômenos fisiológicos que os modelos existentes são incapazes de simular. Estudos adicionais foram feitos para verificar a viabilidade do modelo para simular algumas doenças que danificam a laringe. / This work describes a mathematical model to simulate the human larynx during a phonation. The objective was to produce a large-scale computational technique to capture several physiological phenomena that take place on the larynx during the vocalization and to assist further studies about this important structure of the human body. Using the finite element methods as the way to discretize the muscle tissue equations of the larynx and the Navier-Stokes equations and a model to describe the collision between both vocal folds, the glottal signal for different larynx geometries with different viscoelastic properties was obtained. The results confirmed the myoelastic-aerodynamic theory which describes the dynamic of the phonation, also reproducing physiologic phenomena that current models are unable to simulate. Additional studies were conducted to confirm the feasibility of the model to simulate some diseases that affect the larynx.

This bachelor’s degree project aimed to improve the wheel profile of a rail mounted trolley and determine the cause of wheel failure. The proceedings of this project where modelled after an approach for solving wear problems with an emphasis on designing for sustainability. A case study and root cause analysis (RCA) was performed and the flanged wheels were deemed insufficient for the given heavy-haul system. Possible areas of wheel profile improvement were identified and further researched with multiple literature reviews. Throughout the projects duration several limitations were introduced that reduced the concept testing to exclusively theoretical prediction models. Archard’s model was implemented to predict wear and operating time for the proposed material and wheel tread profile concepts. The wheel flange dimensions were chosen based on recommendations from wheel and rail interference handbooks among other sources. The final wheel and rail profile suggestion improved operating time by approximately 300% and wear resistance by 50% compared to its predecessor. This result was achieved by applying the same theoretical prediction model to both current and suggested profiles. The findings of this project are meant to aid SCA among others in similar cases and additionally highlight the value of product improvement from a technological, sociological, and environmental perspective.

Submitted by ANA KARLA PEREIRA RODRIGUES (anakarla_@hotmail.com) on 2017-08-23T16:00:02Z No. of bitstreams: 1 arquivototal.pdf: 9090532 bytes, checksum: d4fefb1d97e9c6d585d5d18a33abf752 (MD5) / Made available in DSpace on 2017-08-23T16:00:02Z (GMT). No. of bitstreams: 1 arquivototal.pdf: 9090532 bytes, checksum: d4fefb1d97e9c6d585d5d18a33abf752 (MD5) Previous issue date: 2016-11-30 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES / In this thesis we study controllability results, asymptotic behavior and inverse problem related to some problems of the theory of partial di erential equations. Two particular systems are the focus of the study: the Mindin-Timoshenko system, describing the vibrational motion of a plate or a beam, and the phase eld system describing the temperature and phase of a medium having two distinct physical states. The rst chapter is devoted to the study of the 1-D Mindlin-Timoshenko system with discontinuous coe cient. A Carleman inequality is obtained under the assumption of monotonicity on the beam speed. Subsequently, two applications are provided: the controllability of the control system acting on the boundary and Lipschitzian stability of the inverse problem of recovering a potential from a single measurement of the solution. In the second chapter we consider a contact problem characterized by the behavior of a two-dimensional plate whose board makes contact with a rigid obstacle. The formulation of this problem is presented by the 2-D Mindlin-Timoshenko system with boundary conditions and suitable damping terms. Concerning such system, is proved via penalty techniques, the existence of solution and that the system energy has exponential decay when the time approaches in nity. In the third chapter, the study is aimed at a nonlinear phase- eld system de ned in a real open interval. Here we present some controllability results when a single control acts, by means of Dirichlet conditions, on the temperature equation of the system on one of the endpoints of the interval. To prove the results is used the method of moments, plus a spectral study of operators associated to the system and xed point theory to deal with the nonlinearity. / Nesta tese estudamos resultados de controlabilidade, comportamento assintótico e problema inverso relacionados a alguns problemas da teoria de equações diferenciais parciais. Dois sistemas particulares são foco do estudo: o sistema de Mindin-Timoshenko, que descreve o movimento vibratório de uma placa ou viga, e o sistema de campo de fases que descreve a temperatura e a fase de um meio onde ocorrem dois estados físicos distintos. O primeiro capítulo é dedicado ao estudo do sistema de Mindlin-Timoshenko 1-D com coe ciente descontínuos. Uma desigualdade de Carleman é obtida sob a hipótese de monotonicidade sobre velocidade da viga. Posteriormente, são fornecidas duas aplicações: a controlabilidade do sistema com controles agindo na fronteira e a estabilidade Lipschitziana do problema inverso de recuperar um potencial através de uma única informação obtida sobre a solução. No segundo capítulo consideramos um problema de contato caracterizado pelo comportamento de uma placa bidimensional cujo bordo faz contato com um obstáculo rígido. A formulação deste problema é apresentada pelo sistema de Mindlin-Timoshenko 2-D com condi ções de fronteira e termos de amortecimento (damping) adequados. Sobre tal sistema, é provada, através de técnicas de penalização, a existência de solução e, posteriormente, que sua energia possui decaimento exponencial quando o tempo tende ao in nito. No terceiro capítulo o estudo é voltado a um sistema de campo de fases não-linear de nido em um intervalo aberto real. Neste espaço apresentamos alguns resultados de controlabilidade quando um único controle age, sob condições de Dirichlet, na equação da temperatura em um dos bordos do intervalo. Para provar os resultados é utilizado o método dos momentos, além de uma estudo espectral de operadores associados ao sistema e teoria de ponto xo para lidar com a não-linearidade.

Este trabalho descreve um modelo matemático para simulação da laringe humana durante a fonação. O objetivo foi produzir uma técnica computacional de grande escala de processamento para capturar os fenômenos fisiológicos que ocorrem na laringe durante a vocalização e servir de base para estudos mais aprofundados sobre esta importante estrutura do corpo humano. Usando o método dos elementos finitos como base para discretizar as equações dos tecidos musculares da laringe e das equações de Navier-Stokes, e um modelo de descrição da colisão entre as pregas vocais, o sinal glotal foi obtido a partir de diferentes geometrias de laringes com diferentes propriedades viscoelásticas. Os resultados confirmaram a teoria mioelástica-aerodinâmica que descreve a dinâmica da fonação, reproduzindo inclusive fenômenos fisiológicos que os modelos existentes são incapazes de simular. Estudos adicionais foram feitos para verificar a viabilidade do modelo para simular algumas doenças que danificam a laringe. / This work describes a mathematical model to simulate the human larynx during a phonation. The objective was to produce a large-scale computational technique to capture several physiological phenomena that take place on the larynx during the vocalization and to assist further studies about this important structure of the human body. Using the finite element methods as the way to discretize the muscle tissue equations of the larynx and the Navier-Stokes equations and a model to describe the collision between both vocal folds, the glottal signal for different larynx geometries with different viscoelastic properties was obtained. The results confirmed the myoelastic-aerodynamic theory which describes the dynamic of the phonation, also reproducing physiologic phenomena that current models are unable to simulate. Additional studies were conducted to confirm the feasibility of the model to simulate some diseases that affect the larynx.

The Generalized Finite Element Method (GFEM) is a numerical method based on the Partition of Unity (PU) concept and inspired on both the Partition of Unity Method (PUM) and the hp-Cloud method. According to the GFEM, the PU is provided by first-degree Lagragian interpolation functions, defined over a mesh of elements similar to the Finite Element Method (FEM) meshes. In fact, the GFEM can be considered an extension of the FEM to which enrichment functions can be applied in specific regions of the problem domain to improve the solution. This technique has been successfully employed to solve problems presenting discontinuities and singularities, like those that arise in Fracture Mechanics. However, most publications on the method are related to linear analyses. The present thesis is a contribution to the few studies of nonlinear analyses of Solid Mechanics by means of the GFEM. One of its main topics is the derivation of a segment-to-segment generalized contact element based on the mortar method. Material and kinematic nonlinear phenomena are also considered in the numerical models. An Object-Oriented design was developed for the implementation of a GFEM nonlinear analyses framework written in Python programming language. The results validated the formulation and demonstrate the gains and possible drawbacks observed for the GFEM nonlinear approach. / O Método dos Elementos Finitos Generalizados (MEFG) é um método numérico baseado no conceito de partição da unidade (PU) e inspirado no Método da Partição da Unidade (MPU) e o método das Nuvens-hp. De acordo com o MEFG, a PU é obtida por meio de funções de interpolação Lagragianas de primeiro grau, definidas sobre uma rede de elementos similar àquela do Método dos Elementos Finitos (MEF). De fato, o MEFG pode ser considerado uma extensão do MEF para a qual se pode aplicar enriquecimentos em regiões específicas do domínio, buscando melhorias na solução. Esta técnica já foi aplicada com sucesso em problemas com descontinuidades e singularidades, como os originários da Mecânica da Fratura. Apesar disso, a maioria das publicações sobre o método está relacionada a análises lineares. A presente tese é uma contribuição aos poucos estudos relacionados a análises não-lineares de Mecânica dos Sólidos por meio do MEFG. Um de seus principais tópicos é o desenvolvimento de um elemento de contato generalizado do tipo segmento a segmento baseado no método mortar. Fenômenos não lineares devidos ao material e à cinemática também são considerados nos modelos numéricos. Um projeto de orientação a objetos para a implementação de uma plataforma de análises não-lineares foi desenvolvido, escrito em linguagem de programação Python. Os resultados validam a formulação e demonstram os ganhos e possíveis desvantagens da abordagem a problemas não lineares por meio do MEFG.

This thesis is focused on stress in brain aneurysms. It consists of three parts. First part is aimed for gaining information about the topic from scientific resources. Next part consists of analyses of geometry of cerebral aneurysms on the computed wall stress. Analyses are performed on four basic geometrical models and results are being discussed. The risky areas in each investigated shape have been identified as well as the comparisons of stress between those shapes have been performed and the most dangerous shape among investigated shapes has been determined. Third part investigates the influence of surrounding tissue on the brain aneurysm. Conclusion of this thesis is that brain gray tissue has positive yet negligible effect on the computed wall stress.

Title: Approximation, numerical realization and qualitative analysis of contact problems with friction Author: Tomáš Ligurský Department: Department of Numerical Mathematics Supervisor: prof. RNDr. Jaroslav Haslinger, DrSc., Department of Numerical Mathe- matics Abstract: This thesis deals with theoretical analysis and numerical realization of dis- cretized contact problems with Coulomb friction. First, discretized 3D static contact prob- lems with isotropic and orthotropic Coulomb friction and solution-dependent coefficients of friction are analyzed by means of the fixed-point approach. Existence of at least one solution is established for coefficients of friction represented by positive, bounded and con- tinuous functions. If these functions are in addition Lipschitz continuous and upper bounds of their values together with their Lipschitz moduli are sufficiently small, uniqueness of the solution is guaranteed. Second, properties of solutions parametrized by the coefficient of friction or the load vector are studied in the case of discrete 2D static contact problems with isotropic Coulomb friction and coefficient independent of the solution. Conditions under which there exists a local Lipschitz continuous branch of solutions around a given reference point are established due to two variants of the...