It has been known for decades that fatigue crack propagation in elastic-plastic media is very sensitive to load history since the nonlinear behavior of the material can have a great influence on propagation rates. However, the raw computation of millions of fatigue cycles with nonlinear material behavior on tridimensional structures would lead to prohibitive calculation times. In this respect, we propose a global model reduction strategy, mixing both the a posteriori and a priori approaches in order to drastically decrease the computational cost of these types of problems. First, the small scale yielding hypothesis is assumed, and an a posteriori model reduction of the plastic behavior of the cracked structure is performed. This reduced model provides incrementally the plastic state in the vicinity of the crack front, from which the instantaneous crack growth rate is inferred. Then an additional a priori model reduction technique is used to accelerate even more the time to solution of the whole problem. This a priori approach consists in building incrementally and without any previous calculations a reduced basis specific to the considered test-case, by extracting information from the evolving displacement field of the structure. Then the displacement solutions of the updated crack geometries are sought as linear combinations of those few basis vectors. The numerical method chosen for this work is the finite element method. Hence, during the propagation the spatial discretization of the model has to be updated to be consistent with the evolving crack front. For this purpose, a specific mesh morphing technique is used, that enables to discretize the evolving model geometry with meshes of the same topology. This morphing method appears to be a key component of the model reduction strategy. Finally, the whole strategy introduced above is embedded inside an adaptive approach, in order to ensure the quality of the results with respect to a given accuracy. The accuracy and the efficiency of this global strategy have been shown through several examples; either in bidimensional and tridimensional cases for model crack propagation, including the industrial example of a helicopter structure.

Today in the scientific and technological world, physical and artificial processes are often described by mathematical models which can be used for simulation, optimization or control. As the mathematical models get more detailed and different coupling effects are required to include, usually the dimension of these models become very large. Such large-scale systems lead to large memory requirements and computational complexity. To handle these large models efficiently in simulation, control or optimization model order reduction (MOR) is essential. The fundamental idea of model order reduction is to approximate a large-scale model by a reduced model of lower state space dimension that has the same (to the largest possible extent) input-output behavior as the original system. Recently, the system-theoretic method Balanced Truncation (BT) which was believed to be applicable only to moderately sized problems, has been adapted to really large-scale problems. Moreover, it also has been extended to so-called descriptor systems, i.e., systems whose dynamics obey differential-algebraic equations. In this thesis, a BT algorithm is developed for MOR of index-1 descriptor systems based on several papers from the literature. It is then applied to the setting of a piezo-mechanical system. The algorithm is verified by real-world data describing micro-mechanical piezo-actuators. The whole algorithm works for sparse descriptor form of the system. The piezo-mechanical original system is a second order index-1 descriptor system, where mass, damping, stiffness, input and output matrices are highly sparse. Several techniques are introduced to reduce the system into a first order index-1 descriptor system by preserving the sparsity pattern of the original models. Several numerical experiments are used to illustrate the efficiency of the algorithm.

This Thesis is concerned with the approximation of high order controllers or the controller reduction problem. We firstly consider approximating high-order controllers by low order controllers based on the closed-loop system approximation. By approximating the closed-loop system transfer function, we derive a new parametrized double-sided frequency weighted model reduction problem. The formulas for the input and output weights are derived using three closed-loop system configurations: (i) by placing a controller in cascade with the plant, (ii) by placing a controller in the feedback path, and (iii) by using the linear fractional transformation (LFT) representation. One of the weights will be a function of a free parameter which can be varied in the resultant frequency weighted model reduction problem. We show that by using standard frequency weighted model reduction techniques, the approximation error can be easily reduced by varying the free parameter to give more accurate low order controllers. A method for choosing the free parameter to get optimal results is being suggested. A number of practical examples are used to show the effectiveness of the proposed controller reduction method. We have then considered the relationships between the closed-loop system con gurations which can be expressed using a classical control block diagram or a modern control block diagram (LFT). Formulas are derived to convert a closed-loop system represented by a classical control block diagram to a closed-loop system represented by a modern control block diagram and vice versa.

Représentations séparées pour la simulation multi-échelle du comportementmécanique et de l’endommagement des matériaux composites.Résumé: Le développement de méthodes numériques performantes pour simuler les structurescomposites est un défi en raison de la nature multi-échelle et de la complexité des mécanismed’endommagement de ce type de matériaux. Les techniques classiques de discrétisationvolumique conduisent à des coûts de calcul importants et sont restreintes en pratique à despetites structures.Dans cette thèse, un nouvelle stratégie basée sur une représentation séparée de la solution estexplorée. L’objectif est de proposer un cadre numérique efficace et fiable pour analyser les endommagementsdans les composites stratifiés sous chargements statiques et dynamiques. Ladécomposition propre généralisée (PGD) est utilisée pour construire la solution.Pour traiter l’endommagement, et plus particulière le délaminage, un modèle de zone cohésivea été implémenté dans la PGD. Une approches multi-échelle innovante est également proposéepour simuler le comportement mécanique des composites à microstructure périodique. L’idéeprincipale est de séparer deux échelles : l’échelle du motif périodique (microstructure) et l’échellemacroscopique. Les résultats de la PGD sont très proches des résultats obtenus par la méthodeéléments fini classique. Finalement, la PGD permet de réduire significativement la complexitédes modèles tout en gardant une précision satisfaisante. / Separated representations for the multiscale simulation of the mechanicalbehavior and damages of composite materials.Abstract: The development of efficient simulations for composite structures is very challengingdue to the multiscale nature and the complex damage process of this materials. When usingstandard 3D discretization techniques with advanced models for large structures, the computationalcosts are generally prohibitive.In this thesis, a new strategy based on a separated represenation of the solution is explored todevelop a computationally efficient and reliable numerical framework for the analysis of damagesin laminated composites subjected to quasi-static and dynamic loading. The PGD (Proper GeneralizedDecomposition) is used to build the solution.To treat damage, and especially delamination, a cohesive zone model has been implemented inthe PGD solver. A novel multiscale approach is also proposed to compute the mechanical behaviorof composites with periodic microstructure. The idea is to separate two scales: the scaleof periodic pattern and the macroscopic scale. The PGD results have been compared with theresults obtained with the classcial finite element method. A close agreement is found between thetwo approach and the PGD has significantly reduced the model complexity.

Dans le cadre d’une Cifre avec Airbus Helicopters, le projet a pour but le développement d’un modèle pour le pré-dimensionnement de pièces lamifiées élastomère-métal dont le rôle est critique en termes de conception et de sécurité pour les architectures de rotors d’hélicoptères. Pour cela, un premier volet de la thèse a consisté à caractériser le comportement élasto-dissipatif du matériau d’étude (un butadiène chargé de noir de carbone) via une campagne d’essais statiques et dynamiques, sous différents modes de sollicitations (uniaxiales et biaxiales) et sur un spectre assez large de fréquences, d’amplitudes et de températures. A partir de ces observations, un modèle phénoménologique de comportement hyper-viscoélastique est proposé. Sur la base d’un modèle de Maxwell généralisé, il permet de traduire les phénomènes observés sur la gamme de sollicitations visées. Un accent particulier a été porté sur la prise en compte de l’effet Payne en adoptant une approche originale. Les paramètres du modèle sont identifiés par une méthode robuste et rapide. Le modèle est ensuite développé à la fois dans un code commercial de calcul par éléments finis et dans un outil de calcul basé sur une méthode de réduction de modèles. Enfin, une étude du comportement en fatigue est réalisée à travers une campagne d’essais originaux servant de point de départ à la proposition d’une loi d’endommagement continu. / In association with Airbus Helicopters, the aim of the project is to develop a model to pre-size elastomer-metal laminated devices whose role is critical in terms of design and safety for helicopters rotor architectures. To do so, the first part of this thesis consisted in characterizing the elasto-dissipatice behavior of the studied material (a carbon black filled butadiene rubber) thanks to static and dynamic tests, with several kind of loading (uni-axial and bi-axial) and over a wide range of frequences, amplitudes and temperatures. From these observations, a phenomenological hyper-viscoelastic model is proposed. Based on a generalized Maxwell model, it is able to describe the phenomena over the loading range of concern. A particular focus is made to take in account the Payne effect thanks to an original approach. The model parameters are identified with a fast and robust method. The model is then implemented in a commercial finite element code and in a tool based on a model reduction method. Last, a study of the behaviour in fatigue is performed with an original characterization campaign from which a continuous damage law is proposed.

Dimensionnée de manière semi-empirique, la voie ferrée est un système mécanique dont le comportement dynamique reste difficile à appréhender et à quantifier. Un outil numérique peut alors être à la fois une aide à la conception, en évaluant la performance de nouvelles structures, et un élément de diagnostic sur les voies existantes, en complément de mesures terrain adaptées.L'outil développé dans ce travail s'appuie sur les techniques de résolution des équations de propagation dans les guides d'ondes dans le domaine des fréquences/nombres d'onde. Sa particularité est l'introduction d'un modèle réduit basé sur une sous-structuration périodique de la structure, en ne conservant, pour la résolution des équations de dispersion, que quelques nombres d'onde judicieusement choisis. En s'appuyant sur cette technique de réduction de modèle et son extension à des modèles temporels avec contact mobile, les coûts de calcul et de stockage sont largement diminués. Cela en fait un outil performant et utilisable dans des études d'ingénierie portant sur la voie ferrée.Des étapes de vérifications numériques sur les hypothèses sous-jacentes à la réduction ont été réalisées en construisant un modèle fréquentiel complet par transformée de Floquet. En parallèle, des campagnes d'essais dynamiques (mesures de réceptance et d'accélération sur traverse au passage des trains) ont été réalisées sur une zone de transition entre voie ballastée et voie sur dalle sur ligne à grande vitesse. Les résultats obtenus expérimentalement permettent de valider les simulations dans les domaines temporel et fréquentiel sur les différentes zones de la transition. Le modèle permet de mieux comprendre l'effet de la sous-structure sur le comportement dynamique. Une des grandes avancées de ce travail est également l'introduction d'une stratégie permettant d'utiliser la méthode de réduction pour coupler des zones présentant des tranches différentes. Des outils de post-traitement ont été développés pour mettre en lumière les effets dynamiques générés par la transition. / Railway tracks have evolved over years based on empirical results and their dynamic behavior still can be difficult to assess or to quantify. A numerical model can thus provide assistance in track design by assessing the mechanical performance of new structures, or allowing the diagnostic of existing track parts, as a complement to relevant in-situ measurements.The numerical model developed in this work combines a reduction strategy with numerical techniques used to solve the propagation equations in waveguides in the frequency / wave number domains. The peculiarity of the exposed methodology is the introduction of a model based on a periodic sub structuring of the track, keeping only a few wavelengths carefully chosen to solve dispersion equations. Based on this model reduction technique and its extension to time domain models with moving contact, the computational time and storage capacity required are greatly reduced. Thus, this model is efficient and useful for engineering purposes in railway tracks studies.Numerical validation of the reduction is carried out by building a complete reference model in the frequency domain. In parallel, measurement campaigns (receptance measurement and sleeper acceleration under passing trains) were performed on a transition zone between ballasted and slab tracks on a High Speed Line. These tests are compared to simulation results in both time and frequency domains on different areas of the transition. After successful validation, the model is used to improve understanding of the role of the substructure on the dynamic behavior.The final major development of this work is the introduction of a strategy to extend the reduction to piecewise periodic structures and the development of post-processing tools to highlight the dynamic effects generated by the transition zone.

Des simulations numériques rapides et précises du contact élastohydrodynamique (EHD) sont recherchées pour aider au développement de produits. L'objectif de cette thèse est de proposer un modèle compact pour le problème du contact EHD en appliquant des méthodes de réduction de modèle. Dans ce but l'équation de Reynolds (non-linéaire), l'équation d'élasticité (linéaire) et l'équilibre de la charge, sont résolus dans un système d'équations unique par la méthode de Newton. La réduction s'effectue par projection sur un sous espace de faible dimension, qui repose sur des solutions du système complet. De plus, une approximation du système est effectuée, dans laquelle les matrices du système réduit sont approximées. Pour le problème du contact EHD stationnaire, un algorithme de génération automatique des modèles compacts est présenté. L'algorithme fournit des modèles réduits stables et rapides sur une région de paramètres définies. La méthode de Newton réduite est également étendue aux fluides non-newtoniens. Les résultats du modèle réduit sont en très bon accord avec ceux du système complet, malgré un temps de calcul clairement plus petit. Par ailleurs, une nouvelle formulation pour le problème de contact EHD transitoire est introduite, dans laquelle la région de calcul est adaptée à la taille du contact. Ceci permet d'obtenir des modèles réduits efficaces, en particulier pour des excitations à grandes amplitudes. Alternativement, la méthode "Trajectory-Piecewise-Linear" (TPWL) est appliquée au problème du contact EHD transitoire. Cette méthode permet une accélération du calcul conséquente. / In today's product development process, fast and exact simulational models of complex physical problems gain in significance. The same holds for the elastohydrodynamic (EHD) contact problem. Thus, the objective of this work is to generate a compact model for the EHD contact problem by the application of model order reduction. Thereto, the EHD contact problem, consisting of the nonlinear Reynolds equation, the linear elasticity equation and the load balance, is solved as a monolithic system of equations using Newton's method. The reduction takes place by projection onto a low-dimensional subspace, which is based on solutions of the full system. Moreover, a so-called system approximation is executed at which the reduced system matrices are substituted by less complex surrogates. For the stationary EHD contact problem, an algorithm for the automated generation of the compact model is presented. This algorithm provides fast and numerically stable reduced systems on a given parameter range. Additionally, the reduced Newton method is extended to the consideration of Non-Newtonian fluids whereat highly accurate results are obtained requiring a very low computational time. Furthermore, a new formulation for the transient EHD contact problem is introduced, at which the computational area is adapted to the current contact size. This kind of morphing enables efficient reduced models in particular for excitations of large amplitude. Beside of the reduced Newton-method with system approximation, the method Trajectory Piecewise Linear (TPWL) is applied to the transient EHD contact problem. Here, further speed-up potential arises. Despite a distinctly lower computational time, the reduced model is in very good accordance with the full system.

Ce travail de recherche est une contribution au développement de la méthode Décomposition Propre Généralisée (PGD) à la résolution de problèmes de diffusion-réaction raides dédiés à la mécanique des matériaux. Ce type d’équations est notamment rencontré lors de l’oxydation des matériaux polymères et il est donc nécessaire de mettre en place un outil pour simuler ce phénomène afin de prédire numériquement le vieillissement de certains matériaux composites à matrice organique utilisés dans l’aéronautique. La méthode PGD a été choisie dans cette thèse car elle permet un gain en temps de calcul notable par rapport à la méthode des éléments finis. Néanmoins cette famille d’équations n’a jamais été traitée avec cette méthode. Cette dernière se résume à la recherche de solutions d’Équations aux Dérivées Partielles sous forme séparée. Dans le cas d’un problème 1D transitoire, cela revient à chercher la solution sous la forme d’une représentation séparée espace-temps. Dans le cadre de cette thèse, un outil numérique a été mis en place permettant une flexibilité telle que différents algorithmes peuvent être testés. La diffusion Fickienne 1D est tout d’abord évaluée avec en particulier une discussion sur l’utilisation d’un schéma de type Euler ou Runge-Kutta à pas adaptatif pour la détermination des fonctions temporelles. Le schéma de Runge-Kutta permet de réduire notablement le temps de calcul des simulations.Ensuite, la mise en place de l’outil pour les systèmes d’équation de type diffusion-réaction nécessite des algorithmes de résolution de systèmes non linéaires, couplés et raides. Pour cela, différents algorithmes ont été implémentés et discutés.Dans le cas d’un système non linéaire, l’utilisation de la méthode de Newton-Raphson dans les itérations pour la recherche du nouveau mode permet de réduire le temps de calcul en limitant le nombre de modes à considérer pour une erreur donnée. En ce qui concerne les couplages, deux stratégies de résolution ont été évaluées. Le couplage fort mène aux mêmes conclusions que dans le cas non linéaire. Les systèmes raides mais linéaires ont ensuite été traités en implémentant l’algorithme de Rosenbrock pour la détermination des fonctions temporelles. Cet algorithme permet contrairement à Euler et à Runge-Kutta de construire une solution avec un temps de calcul raisonnable liée à l’adaptation du maillage temporel sous-jacent à l’utilisation de cette méthode. La résolution d’un système d’équations de diffusion-réaction raides non linéaires utilisée pour la prédiction de l’oxydation d’un composite issu de la littérature a été testée en utilisant les différents algorithmes mis en place. Néanmoins, les non linéarités et la raideur du système génèrent des équations différentielles intermédiaires à coefficients variables pour lesquelles la méthode de Rosenbrock montre ses limites. Il sera donc nécessaire de tester ou développer d’autres algorithmes pour lever ce verrou.Mots / This work presents the development of the Proper Generalized Decomposition (PGD) method for solving stiff reaction-diffusion equations in the framework of mechanics of materials. These equations are particularly encountered in the oxidation of polymers and it is therefore necessary to develop a tool to simulate this phenomenon for example for the ageing of organic matrix composites in aircraft application. The PGD method has been chosen in this work since it allows a large time saving compared to the finite element method. However this family of equations has never been dealt with this method. The PGD method consists in approximating a solution of a Partial Differential Equation with a separated representation. The solution is sought under a space-time separated representation for a 1D transient equation.In this work, a numerical tool has been developed allowing a flexibility to test different algorithms. The 1D Fickian diffusion is first evaluated and two numerical schemes, Euler and Runge-Kutta adaptive methods, are discussed for the determination of the time modes. The Runge-Kutta method allows a large time saving. The implementation of the numerical tool for reaction-diffusion equations requires the use of specific algorithms dedicated to nonlinearity, couplingand stiffness. For this reason, different algorithms have been implemented and discussed. For nonlinear systems, the use of the Newton-Raphson algorithm at the level of the iterations to compute the new mode allows time saving by decreasing the number of modes required for a given precision. Concerning the couplings, two strategies have been evaluated. The strong coupling leads to the same conclusions as the nonlinear case. The linear stiff systems are then studied by considering a dedicated method, the Rosenbrock method, for the determination of the time modes. This algorithm allows time saving compared to the Runge-Kutta method. The solution of a realistic nonlinear stiff reaction-diffusionsystem used for the prediction of the oxidation of a composite obtained from the literature has been tested by using the various implemented algorithms. However, the nonlinearities and the stiffness of the system generate differential equations with variable coefficients for which the Rosenbrock method is limited. It will be necessary to test or develop other algorithms to overcome this barrier.

Dans le domaine de l’ingénierie (par exemple l’aéronautique, l’automobile, la biologie, les circuits), les systèmes dynamiques sont le cadre de base utilisé pour modéliser, contrôler et analyser une grande variété de systèmes et de phénomènes. En raison de l’utilisation croissante de logiciels dédiés de modélisation par ordinateur, la simulation numérique devient de plus en plus utilisée pour simuler un système ou un phénomène complexe et raccourcir le temps de développement et le coût. Cependant, le besoin d’une précision de modèle améliorée conduit inévitablement à un nombre croissant de variables et de ressources à gérer au prix d’un coût numérique élevé. Cette contrepartie justifie la réduction du modèle. Pour les systèmes linéaires invariant dans le temps, plusieurs approches de réduction de modèle ont été effectivement développées depuis les années 60. Parmi celles-ci, les méthodes basées sur l’interpolation se distinguent par leur souplesse et leur faible coût de calcul, ce qui en fait un candidat prédestiné à la réduction de systèmes véritablement à grande échelle. Les progrès récents démontrent des façons de trouver des paramètres de réduction qui minimisent localement la norme H2 de l’erreur d’incompatibilité. En général, une approximation d’ordre réduit est considérée comme un modèle de dimension finie. Cette représentation est assez générale et une large gamme de systèmes dynamiques linéaires peut être convertie sous cette forme, du moins en principe. Cependant, dans certains cas, il peut être plus pertinent de trouver des modèles à ordre réduit ayant des structures plus complexes. A titre d’exemple, certains systèmes de phénomènes de transport ont leurs valeurs singulières Hankel qui se décomposent très lentement et ne sont pas facilement approchées par un modèle de dimension finie. En outre, pour certaines applications, il est intéressant de disposer d’un modèle structuré d’ordre réduit qui reproduit les comportements physiques. C’est pourquoi, dans cette thèse, les modèles à ordre réduit ayant des structures de retard ont été plus précisément considérés. Ce travail a consisté, d’une part, à développer de nouvelles techniques de réduction de modèle pour des modèles à ordre réduit avec des structures de retard et, d’autre part, à trouver de nouvelles applications d’approximation de modèle. La contribution majeure de cette thèse couvre les sujets d’approximation et inclut plusieurs contributions au domaine de la réduction de modèle. Une attention particulière a été accordée au problème de l’approximation du modèle optimale pour les modèles structurés retardés. À cette fin, de nouveaux résultats théoriques et méthodologiques ont été obtenus et appliqués avec succès aux repères académiques et industriels. De plus, la dernière partie de ce manuscrit est consacrée à l’analyse de la stabilité des systèmes retardés par des méthodes interpolatoires. Certaines déclarations théoriques ainsi qu’une heuristique sont développées permettant d’estimer de manière rapide et précise les diagrammes de stabilité de ces systèmes. / In the engineering area (e.g. aerospace, automotive, biology, circuits), dynamical systems are the basic framework used for modeling, controlling and analyzing a large variety of systems and phenomena. Due to the increasing use of dedicated computer-based modeling design software, numerical simulation turns to be more and more used to simulate a complex system or phenomenon and shorten both development time and cost. However, the need of an enhanced model accuracy inevitably leads to an increasing number of variables and resources to manage at the price of a high numerical cost. This counterpart is the justification for model reduction. For linear time-invariant systems, several model reduction approaches have been effectively developed since the 60’s. Among these, interpolation-based methods stand out due to their flexibility and low computational cost, making them a predestined candidate in the reduction of truly large-scale systems. Recent advances demonstrate ways to find reduction parameters that locally minimize the H2 norm of the mismatch error. In general, a reduced-order approximation is considered to be a finite dimensional model. This representation is quite general and a wide range of linear dynamical systems can be converted in this form, at least in principle. However, in some cases, it may be more relevant to find reduced-order models having some more complex structures. As an example, some transport phenomena systems have their Hankel singular values which decay very slowly and are not easily approximated by a finite dimensional model. In addition, for some applications, it is valuable to have a structured reduced-order model which reproduces the physical behaviors. That is why, in this thesis, reduced-order models having delay structures have been more specifically considered. This work has focused, on the one hand, in developing new model reduction techniques for reduced order models having delay structures, and, on the other hand, in finding new applications of model approximation. The major contribution of this thesis covers approximation topics and includes several contributions to the area of model reduction. A special attention was given to the H2 optimal model approximation problem for delayed structured models. For this purpose, some new theoretical and methodological results were derived and successfully applied to both academic and industrial benchmarks. In addition, the last part of this manuscript is dedicated to the analysis of time-delayed systems stability using interpolatory methods. Some theoretical statements as well as an heuristic are developed enabling to estimate in a fast and accurate way the stability charts of those systems.

La solution numérique des systèmes dynamiques est un moyen efficace pour étudier des phénomènes physiques complexes. Cependant, dans un cadre à grande échelle, la dimension du système rend les calculs infaisables en raison des limites de mémoire et de temps, ainsi que le mauvais conditionnement. La solution de ce problème est la réduction de modèles. Cette thèse porte sur les méthodes de projection pour construire efficacement des modèles d'ordre inférieur à partir des systèmes linéaires dynamiques de grande taille. En particulier, nous nous intéressons à la projection sur la réunion de plusieurs sous-espaces de Krylov standard qui conduit à une classe de modèles d'ordre réduit. Cette méthode est connue par l'interpolation rationnelle. En se basant sur ce cadre théorique qui relie la projection de Krylov à l'interpolation rationnelle, quatre algorithmes de type Lanczos rationnel pour la réduction de modèles sont proposés. Dans un premier temps, nous avons introduit une méthode adaptative de type Lanczos rationnel par block pour réduire l'ordre des systèmes linéaires dynamiques de grande taille, cette méthode est basée sur l'algorithme de Lanczos rationnel par block et une méthode adaptative pour choisir les points d'interpolation. Une généralisation de ce premier algorithme est également donnée, où différentes multiplicités sont considérées pour chaque point d'interpolation. Ensuite, nous avons proposé une autre extension de la méthode du sous-espace de Krylov standard pour les systèmes à plusieurs-entrées plusieurs-sorties, qui est le sous-espace de Krylov global. Nous avons obtenu des équations qui décrivent cette procédure. Finalement, nous avons proposé une méthode de Lanczos étendu par block et nous avons établi de nouvelles propriétés algébriques pour cet algorithme. L'efficacité et la précision de tous les algorithmes proposés, appliqués sur des problèmes de réduction de modèles, sont testées dans plusieurs exemples numériques. / Numerical solution of dynamical systems have been a successful means for studying complex physical phenomena. However, in large-scale setting, the system dimension makes the computations infeasible due to memory and time limitations, and ill-conditioning. The remedy of this problem is model reductions. This dissertations focuses on projection methods to efficiently construct reduced order models for large linear dynamical systems. Especially, we are interesting by projection onto unions of Krylov subspaces which lead to a class of reduced order models known as rational interpolation. Based on this theoretical framework that relate Krylov projection to rational interpolation, four rational Lanczos-type algorithms for model reduction are proposed. At first, an adaptative rational block Lanczos-type method for reducing the order of large scale dynamical systems is introduced, based on a rational block Lanczos algorithm and an adaptive approach for choosing the interpolation points. A generalization of the first algorithm is also given where different multiplicities are consider for each interpolation point. Next, we proposed another extension of the standard Krylov subspace method for Multiple-Input Multiple-Output (MIMO) systems, which is the global Krylov subspace, and we obtained also some equations that describe this process. Finally, an extended block Lanczos method is introduced and new algebraic properties for this algorithm are also given. The accuracy and the efficiency of all proposed algorithms when applied to model order reduction problem are tested by means of different numerical experiments that use a collection of well known benchmark examples.