Cette thèse traite de la réduction d’ordre de modèle du phénomène communément rencontré dans la modélisation de microsystèmes, à savoir, dans la littérature anglaise, le « squeeze-film damping ». Dans un premier chapitre sont présentées les différentes méthodes de réduction d’ordre de modèle. Dans le cas des systèmes linéaires, elles ont un cadre théorique bien établi. Ces méthodes peuvent être adaptées pour les systèmes non-linéaires. La validité des modèles réduits résultants sera alors réduite à un certain espace des phases, leur établissement faisant intervenir certaines trajectoires particulières servant d’apprentissage. On présente finalement la méthode des modes normaux non-linéaires dont les modèles résultants ne dépendent pas d’une trajectoire d’apprentissage. Au chapitre 2, on s’intéresse plus particulièrement au phénomène de « squeeze-film damping » régi par l’équation de Reynolds. Après avoir détaillé son établissement à partir de certaines hypothèses, on décrit les différentes méthodes de résolution de l’équation linéaire puis non-linéaire de la littérature. On compare ensuite les résultats d’un modèle de l’équation de Reynolds à des simulations éléments finis de l’équation de Navier-Stokes afin de valider les hypothèses faites pour la dérivation de l’équation de Reynolds. On propose ensuite une résolution originale par changement de variable. On étudie aussi plusieurs autres résolutions possibles ainsi que plusieurs bases de projection parmi celles décrites dans le premier chapitre. Le chapitre 3 est consacré à la modélisation du problème couplé que constitue le micro-interrupteur MEMS qui est un candidat au remplacement des interrupteurs à base de transistors dans les communications RF. Sa modélisation fait intervenir trois domaines, la mécanique, l’électrostatique, et la fluidique à travers l’équation de Reynolds. Après voir décrit les différents modèles de la littérature, on propose un modèle réduit couplé dont le modèle fluidique est basé sur le modèle établi au chapitre 2. Ce modèle est validé par rapport à des modèles différences finies et à des résultats expérimentaux de la littérature.Enfin le quatrième chapitre traite de la réduction du coût d’évaluation du modèle réduit couplé de micro-interrupteur du chapitre 3. La première méthode proposée consiste à trouver une fonction d’approximation de la projection de la force fluidique sur le premier mode mécanique, fonction des coordonnées modales mécaniques position et vitesse. Cette méthode ne se révèle valable que dans le cas incompressible. Dans le cas compressible, la résolution de l’équation de Reynolds restant obligatoire, on utilise la méthode de Rewienski et al. qui consiste à linéariser par morceaux les fonctions régissant la dynamique. Une autre méthode de linéarisation par morceaux, tirant parti d’une particularité du modèle du chapitre 2 et permettant de s’affranchir d’une trajectoire d’apprentissage, est également proposée. / This thesis deals with reduced-order modelling of squeeze-film damping, a fluidic phenomenon that is commonly encountered in MEMS. In the first chapter, reduced-order modelling methods are presented. For linear systems, well-established theories exist. They can be adapted to nonlinear systems. However, the resulting reduced-order models are valid in a certain region of the state-space only, depending on the training trajectory. The method of nonlinear normal modes, which does not depend on a training trajectory is also introduced. Chapter two is focused on the squeeze-film damping phenomenon governed by the Reynolds equation. We first establish the equation from appropriate hypotheses, and then present the different resolutions of its linear and nonlinear form found in literature. The results from a model based on the Reynolds equation are then compared to results from a finite element Navier-Stokes model, in order to validate the various hypotheses made. An original method of resolution based on a change of variable is then proposed. Several other method of resolution are studied as well as different projection bases amongst those presented in chapter one.Chapter three is dedicated to the modelling a micro-switch, a candidate to the replacement of switches based on transistors in RF communications. Its modelling implies the coupling of three domains: mechanics, electrostatics, and fluidics with Reynolds equation. Following a description of the models from literature, a coupled model is proposed, the fluidic model being the one presented in chapter two. This model is validated compared to finite difference models as well as experimental data from the literature.Finally, the fourth chapter aims at reducing the evaluation cost of the coupled micro-switch model established in chapter three. The first method consists in finding an approximation function of the projection of the fluidic force on the first linear mechanical mode as a function of the mechanical modal coordinates, position and speed. This method is applicable in the incompressible case only. In the compressible case, the Reynolds equation has to be solved. The method of Rewienski and al., which consists in piecewise-linearizing the functions governing the dynamics, is used. Another method based on a piecewise-linear approach, taking advantage of the particular structure of the model presented in chapter two, thus not depending on a training trajectory, is proposed.

For repeated transient solutions of geometrically nonlinear structures the numerical effort often poses a major obstacle. Thus, the introduction of a reduced order model, which takes the nonlinear effects into account and accelerates the calculations considerably, is often necessary.This work yields a method that allows for rapid, accurate and parameterisable solutions by means of a reduced model of the original structure. The structure is discretised and its dynamic equilibrium described by a matrix equation. The projection on a reduced basis is introduced to obtain the reduced model. A comprehensive numerical study on several common reduced bases shows that the simple introduction of a constant basis is not sufficient to account for the nonlinear behaviour. Three requirements for an rapid, accurate and parameterisable solution are derived. The solution algorithm has to take into account the nonlinear evolution of the solution, the solution has to be independent of the nonlinear finite element terms and the basis has to be adapted to external parameters.Three approaches are provided, each responding to one requirement. These approaches are assembled to the integrated method. The approaches are the update and augmentation of the basis, the polynomial formulation of the nonlinear terms and the interpolation of the basis. A Newmark-type time-marching algorithm provides the frame of the integrated method. The application of the integrated method on test-cases with geometrically nonlinear finite elements confirms that this method leads to the initial aim of a rapid, accurate and parameterisable transient solution.

Pour les solutions transitoires répétées des structures géométriquement nonlinéaires l’effort numérique présente souvent une contrainte importante. Ainsi, l’introduction d’un modèle d’ordre réduit, qui prend en compte les effets nonlinéaires et qui accélère considérablement les calculs, s’avère souvent nécessaire.Ce travail aboutit à une méthode qui permet des solutions transitoires accélérées, fidèles et paramétrables, à travers d’un modèle réduit de la structure initiale. La structure est discrétisée et son équilibre dynamique décrit par une équation matricielle. La projection sur une base réduite est introduite afin d’obtenir un modèle réduit. Une étude numérique complète sur plusieurs bases communes démontre que la simple introduction d’une base constante ne suffit pas pour prendre en compte le comportement nonlinéaire. Trois exigences sont déduites pour une solution transitoire accélérée, fidèle et paramétrable. L’algorithme de solution doit permettre un suivi de l’évolution nonlinéaire de la solution transitoire, la solution doit être autonome des termes nonlinéaires en éléments finis et la base doit être adaptée à des paramètres externes.Trois approches sont mises en place, chacune répondant à une exigence. Ces approches sont assemblées dans la méthode intégrée. Les approches sont la mise-à-jour et augmentation de la base , la formulation polynomiale des termes nonlinéaires et l’interpolation de la base. Un algorithme de type Newmark forme le cadre de la méthode intégrée. L’application de la méthode intégrée sur des cas test en élément finis géométriquement nonlinéaires confirme qu’elle répond au but initial d’obtenir des solutions transitoires accélérées, fidèles et paramétrables. / For repeated transient solutions of geometrically nonlinear structures the numerical effort often poses a major obstacle. Thus, the introduction of a reduced order model, which takes the nonlinear effects into account and accelerates the calculations considerably, is often necessary.This work yields a method that allows for rapid, accurate and parameterisable solutions by means of a reduced model of the original structure. The structure is discretised and its dynamic equilibrium described by a matrix equation. The projection on a reduced basis is introduced to obtain the reduced model. A comprehensive numerical study on several common reduced bases shows that the simple introduction of a constant basis is not sufficient to account for the nonlinear behaviour. Three requirements for an rapid, accurate and parameterisable solution are derived. The solution algorithm has to take into account the nonlinear evolution of the solution, the solution has to be independent of the nonlinear finite element terms and the basis has to be adapted to external parameters.Three approaches are provided, each responding to one requirement. These approaches are assembled to the integrated method. The approaches are the update and augmentation of the basis, the polynomial formulation of the nonlinear terms and the interpolation of the basis. A Newmark-type time-marching algorithm provides the frame of the integrated method. The application of the integrated method on test-cases with geometrically nonlinear finite elements confirms that this method leads to the initial aim of a rapid, accurate and parameterisable transient solution.