Les travaux de recherche que je mène depuis une dizaine d'années, au Laboratoire d'Ingénierie des Matériaux de Bretagne, portent principalement sur deux domaines : la modélisation des procédés de mise en forme des matériaux métalliques et la simulation du comportement de matériaux polymères. Sachant que les conditions d'emplois et de mise en forme des matériaux sont de plus en plus sévères, l'objectif de ces thématiques est d'améliorer la caractérisation et la modélisation du comportement mécanique des matériaux. Tout d'abord, en cherchant à intégrer des phénomènes mécaniques mal maîtrisés (effets de la vitesse de chargement élevée, sensibilité à la température et dépendance du chemin de chargement) ensuite, en développant des méthodes expérimentales ad hoc permettant de mesurer ces phénomènes, et enfin, en mettant en oeuvre des stratégies numériques originales, en tentant de simuler au mieux ces phénomènes. Ces travaux, s'effectuant dans des domaines de transformations finies, un formalisme en grandes transformations s'appuyant sur la description du mouvement du solide par les coordonnées matérielles entraînées a d'abord été précisé. L'étude de l'intégration temporelle du taux de déformation dans différents repères objectifs a montré l'importance du repère en rotation logarithmique et que la dérivée dans le repère corotationnel n'en était qu'une approximation. Ces restrictions à l'esprit, la modélisation de la prédiction du retour élastique lors de la mise en forme des alliages d'aluminium et de la prise en compte de la température au cours de ces étapes a pu être évaluée à travers le test d'ouverture d'un anneau découpé dans un godet cylindrique. L'influence combinée des modèles de comportement, des paramètres numériques, des modes de découpe numérique et du type de modélisation thermomécanique a été analysée, montrant qu'il reste encore des progrès à réaliser pour la modélisation du retour élastique. Enfin, un modèle de comportement d'Hyperélasto-Visco-Hystérétique, discernant les phénomènes mécaniques observés grâce à plusieurs contributions en contrainte, a été appliqué sur deux matériaux polymères. Une méthode d'identification simple, rapide et nécessitant un seul essai de chargement interrompu par des phases de relaxation a permis de montrer l'existence dans ces matériaux des effets de la viscosité, de l'hyperélasticité et d'une partie hystérétique non-visqueuse. La prise en compte des effets de la température s'est également avérée concluante.

Dans cette thèse, nous nous sommes intéressés au comportement vibratoire d’un véhicule soumis à une excitation moteur dans deux plages de régimes différentes : basses fréquences (0 – 50 Hz) et moyennes fréquences (200 – 800 Hz). Le but était de fournir des méthodologies numériques permettant de prendre en compte les phénomènes de couplage vibratoires existant entre les différents sous-systèmes constitutifs d’une caisse automobile.En basses fréquences, nous avons adopté une approche globale où chaque sous-système a été caractérisé séparément. Tout d’abord, le comportement de la caisse a été caractérisé expérimentalement et numériquement par une méthode jusqu’alors réservée au domaine aéronautique dite d’appropriation modale. Les résultats numériques ont alors été confrontés aux résultats expérimentaux. Par ailleurs, le comportement non-linéaire en amplitude et en fréquence des pièces de filtration moteur (SMO) a été déterminé sur banc de mesure. Un fort comportement non-linéaire a pu être constaté et ces caractérisations ont été exploitées en construisant des nappes raideur-fréquence-amplitude. Dans un second temps, des méthodes numériques permettant de réaliser l’assemblage non - linéaire de la caisse et du groupe moto – propulseur (GMP) via les pièces de filtration non linéaires ont été mises en place. Pour ce faire, nous avons développé une méthode dite de Balance Harmonique (HBM) qui permet de prédire la dynamique non-linéaire de systèmes complexes. Afin d’appliquer cette méthode à une structure industrielle, nous avons utilisé une méthode de condensation sur les degrés de liberté non-linéaires, technique bien adaptée aux cas de structures linéaires reliées localement par des éléments de liaison non-linéaires. Cette méthode a tout d’abord été validée sur un périmètre restreint comprenant un GMP relié à un bâti rigide par ses pièces de filtration. A cette occasion, des phénomènes non-linéaires importants ont été constatés expérimentalement. Un modèle numérique de GMP a été construit et l’utilisation de la méthode HBM a permis de retrouver ces constats. Enfin, après avoir réalisé l’assemblage non-linéaire des trois sous-systèmes GMP - SMO - Caisse, la structure a été excitée de plusieurs manières différentes : appropriation numérique non-linéaire et excitation réelle d’un GMP. En moyennes fréquences, nous avons présenté dans ce mémoire une étude importante pour le groupe Renault concernant la caractérisation des pièces de filtration en moyennes fréquences. Au cours de la thèse, une méthodologie numérique basée sur la méthode FBS permettant de déconfiner (ou découpler) une suspension moteur initialement reliée à un banc de mesure a été proposée. La faisabilité numérique du déconfinement a ainsi été démontrée. Cette méthode permet donc, en dépit de phénomènes de couplage avec le banc de mesure, d’obtenir le comportement vibratoire intrinsèque de la pièce. / In this thesis we studied the vibratory behaviour of a whole vehicle under engine excitation at low frequencies (0 – 50 Hz) and medium frequencies (200 – 800 Hz). The aim of the thesis was to provide numerical methodologies to take into account coupling effects between all the sub-systems constituting a whole car. In low frequencies, we used a global approach where each subsystem was characterized separately before coupling. First the car body was characterised both experimentally and numerically using a modal appropriation method that is commonly used in the aeronautic field. Numerical shapes of the modes were correlated to experimental shapes. In addition, the amplitude and frequency non linear behaviour of the engine mounts was measured on a test bench. A strong non linear behaviour was observed and stiffness – frequency – amplitude layers were constructed based on those data.Secondly, numerical methods were developed in order to calculate the coupled non linear response between the engine, the engine mounts and the car body. We used a harmonic balance method that allows calculating the non linear dynamics of complex mechanical systems. In order to apply this method to large industrial finite element models, a condensation method on non linear degrees of freedom was developed. This technique is well adapted to problems of linear structures linked together with localnon linear joints. This method was validated on the isolated engine linked to a bench by the engine mounts. Strong non linear phenomena on the rigid body modes of the engine were observed experimentally.A numerical model of the engine was developed and the HBM method allowed reproducing these non linear phenomena. Eventually, the non linear model of the whole vehicle was coupled and excited by different efforts. First we calculated the response of the assembly using the appropriation method. Then, the structure was excited by a real four – cylinder engine excitation.In medium frequencies, we presented an important study for the group Renault concerning the stiffness measurement of the engine mounts. A numerical methodology based on the FRF Based Substructuring(FBS) method was developed. This method was applied to uncouple an engine mount initially coupled to a test bench. The numerical feasibility of the method was proved and allowed to get the own vibratory behaviour of the engine mount despite coupling phenomena with the test bench.

In this diploma thesis is created a cranktrain model in the MBS system ADAMS Engine. The thesis defines basic parameters and proportions of the torsional damper for in-line 5-cylinder engine. This thesis also puts mind to calculation of torsional vibration. In the MBS system ADAMS Engine is demonstrated affect of torsional damper on the torsional vibration of the crankshaft.

The subject of the thesis is the investigation of asymptotic properties of solutions of the Cauchy problem for structurally damped sigma-evolution operators with time dependent, monotonous, dissipation term. An appropriate energy for solutions of the sigma-evolution equations is defined and some estimates for energies of higher order are proved. In the scale invariant case the optimality of these estimates is shown. Further, the influence of properties of the time dependent dissipation on L^p-L^q estimates for the energy with p and q bigger or equal to 2 and from the conjugate line is clarified. Also smoothing properties of the operators under consideration are investigated. The connection between the regularity of the data and the regularity of the solution in terms of L^2 based Gevrey spaces is considered. Finally, L^1-L^1-estimates in the special case delta = sigma/2 and decreasing dissipative coefficient. / Thema der vorliegenden Dissertation ist die Untersuchung asymptotischer Eigenschaften von Lösungen des Cauchy Problems für strukturell gedämpfte sigma-Evolutions-Operatoren mit zeitabhängigem, monotonen Dissipationskoeffizienten. Es wird eine geeignete Energie definiert und für diese Abschätzungen, auf für entsprechende Energien höherer Ordnung gezeigt. Darüber hinaus wird der Einfluss des Dissipationskoeffizienten auf L^p-L^q Abschätzungen auf und entfernt von der konjugierten Linie untersucht. Im skaleninvarianten Fall wird die Schärfe der Abschätzungen bewiesen. Weiterhin wird der Zusammenhang zwischen der Regularität der Daten und der der Lösung in Termen von L^2-basierten Gevrey-Räumen untersucht. Schließlich werden L^1-L^1-Abschätzungen für den Spezialfall delta = sigma/2 und monoton fallenden Dissipationskoeffizienten gezeigt.

Ces travaux de thèse sont composés principalement de deux grandes parties de volumes inégaux dont voici le résumé des contenus : Partie I : Elle concerne les problèmes de stabilité et de bifurcation rencontrés dans les matériaux indépendants du temps physique (i.e., élastiques ou élasto-plastiques). Une analyse approfondie et détaillée des différentes études dans ce domaine nous permet de faire une présentation originale et synthétique de la théorie de la stabilité et de la bifurcation. Une illustration de cette théorie est ensuite donnée à travers une collaboration industrielle avec la SNECMA. On y effectue des simulations numériques du formage des aubes de réacteurs d'avions accompagnées d'analyses des états critiques de flambage. Les conclusions de ces simulations montrent clairement la nécessité de tenir compte du comportement visqueux du matériau et motivent ainsi la seconde partie de la thèse. Partie II : La deuxième partie des travaux de thèse est la plus importante quant au volume et au fond. Elle concerne les problèmes de stabilité des matériaux dépendant du temps physique (i.e., visco-élastiques ou visco-plastiques). L'objectif de cette étude est double : D'une part, une modélisation théorique du problème de stabilité des solides visqueux est proposée. Cette modélisation a pour but de poser le problème étudié comme un problème de stabilité de l'évolution quasi statique du solide visqueux considéré. En effet, l'absence d'équilibres pour de telles structures visqueuses nous amène naturellement à étudier la stabilité de leurs évolutions quasi statiques. Ces dernières dépendent bien entendu de la vitesse du chargement, ce qui constitue une source de difficultés dans une telle analyse. Cette façon de poser le problème en termes de stabilité de trajectoires est, à notre connaissance, originale et généralise la notion plus classique de la stabilité d'un équilibre. Elle consiste, en quelque sorte, à mesurer la sensibilité des mouvements aux perturbations extérieures. D'autre part, la recherche de critères assurant la stabilité de l'évolution visqueuse est menée en étroite liaison avec les résultats mathématiques relatifs aux équations différentielles non autonomes. Il est mis en évidence, notamment, que la principale difficulté dans cette analyse vient du caractère non autonome des équations différentielles régissant l'évolution quasi statique visqueuse. Les conditions suffisantes de stabilité auxquelles on aboutit sont obtenues principalement par deux méthodes. La première consiste à construire des fonctionnelles de Lyapunov adéquates à la physique du problème ou à utiliser la méthode de linéarisation autour de la réponse quasi statique. L'autre démarche est plus directe et permet d'établir le critère de stabilité de l'évolution par des considérations de continuité de la réponse par rapport aux données initiales. Cette deuxième démarche a l'avantage de s'appliquer aussi bien aux évolutions visco-plastiques qu'aux évolutions élasto-plastiques. Enfin, l'application de ces différents résultats aux modèles de tiges de Shanley et pour le cas des poutres permet d'illustrer notre méthodologie sur des exemples concrets et de comprendre la portée des critères de stabilité ainsi obtenus.

The primary focus of this thesis is to examine the detailed seismic structure of the northern Cascadia margin, including the Cascadia basin, the deformation front and the continental shelf. The results of this study are contributing towards understanding sediment deformation and tectonics on this margin. They also have important implications for exploration of hydrocarbons (oil and gas) and natural hazards (submarine landslides, earthquakes, tsunamis, and climate change). The first part of this thesis focuses on the role of gas hydrate in slope failure observed from multibeam bathymetry data on a frontal ridge near the deformation front off Vancouver Island margin using active-source ocean bottom seismometer (OBS) data collected in 2010. Volume estimates (∼ 0.33 km^3) of the slides observed on this margin indicate that these are capable of generating large (∼ 1 − 2 m) tsunamis. Velocity models from travel time inversion of wide angle reflections and refractions recorded on OBSs and vertical incidence single channel seismic (SCS) data were used to estimate gas hydrate concentrations using effective medium modeling. Results indicate a shallow high velocity hydrate layer with a velocity of 2.0 − 2.1 km/s that corresponds to a hydrate concentration of 40% at a depth of 100 m, and a bottom simulating reflector (BSR) at a depth of 265 − 275 m beneath the seafloor (mbsf). These are comparable to drilling results on an adjacent frontal ridge. Margin perpendicular normal faults that extend down to BSR depth were also observed on SCS and bathymetric data, two of which coincide with the sidewalls of the slump indicating that the lateral extent of the slump is controlled by these faults. Analysis of bathymetric data indicates, for the first time, that the glide plane occurs at the same depth as the shallow high velocity layer (100±10 mbsf). In contrast, the glide plane coincides with the depth of the BSR on an adjacent frontal ridge. In either case, our results suggest that the contrast in sediments strengthened by hydrates and overlying or underlying sediments where there is no hydrate is what causing the slope failure on this margin. The second part of this dissertation focuses on obtaining the detailed structure of the Cascadia basin and frontal ridge region using mirror imaging of few widely spaced OBS data. Using only a small airgun source (120 cu. in.), our results indicate structures that were previously not observed on the northern Cascadia margin. Specifically, OBS migration results show dual-vergence structure, which could be related to horizontal compression associated with subduction and low basal shear stress resulting from over-pressure. Understanding the physical and mechanical properties of the basal layer has important implications for understanding earthquakes on this margin. The OBS migrated image also clearly shows the continuity of reflectors which enabled the identification of thrust faults, and also shows the top of the igneous oceanic crust at 5−6 km beneath the seafloor, which were not possible to identify in single-channel and low-fold multi-channel seismic (MCS) data. The last part of this thesis focuses on obtaining detailed seismic structure of the Vancouver Island continental shelf from MCS data using frequency domain viscoacoustic full waveform inversion, which is first of its kind on this margin. Anelastic velocity and attenuation models, derived in this study to subseafloor depths of ∼ 2 km, are useful in understanding the deformation within the Tofino basin sediments, the nature of basement structures and their relationship with underlying accreted terranes such as the Crescent and the Pacific Rim terranes. Specifically, our results indicate a low-velocity zone (LVZ) with a contrast of 200 m/s within the Tofino basin sediment section at a depth 600 − 1000 mbsf over a lateral distance of 10 km. This LVZ is associated with high attenuation values (0.015 − 0.02) and could be a result of over pressured sediments or lithology changes associated with a high porosity layer in this potential hydrocarbon environment. Shallow high velocities of 4 − 5 km/s are observed in the mid-shelf region at depths > 1.5 km, which is interpreted as the shallowest occurrence of the Eocene volcanic Crescent terrane. The sediment velocities sharply increase about 10 km west of Vancouver Island, which probably corresponds to the underlying transition to the Mesozoic marine sedimentary Pacific Rim terrane. High attenuation values of 0.03 − 0.06 are observed at depths > 1 km, which probably corresponds to increased clay content and the presence of mineralized fluids. / Graduate / 0373 / 0372 / 0605 / subbarao@uvic.ca

Les travaux de recherche présentés dans cette HDR sont répartis en trois thématiques complémentaires. La première concerne la modélisation des instabilités locales dans les matériaux métalliques (localisation, striction ...). Il s'agit de développer des outils théoriques et numériques de prédiction de ces phénomènes et leur validation au travers de courbes limites de formage pour des aciers ferritiques ou dual-phase. Deux approches complémentaires ont été adoptées : la première est micromécanique se basant sur la plasticité cristalline couplée à des techniques de transition d'échelles, tandis que la seconde est phénoménologique utilisant des modèles de comportement avancés que nous avons couplé à de l'endommagement. Le critère de bifurcation de Rice, qui est relié à la perte d'ellipticité des équations gouvernant le problème aux limites, a été particulièrement analysé au travers de ces deux approches. L'influence sur l'apparition d'une localisation de plusieurs paramètres matériau ou liés à la microstructure ainsi que l'effet de certains mécanismes déstabilisants ont été analysés. Il est mis en évidence, notamment, le rôle tout à fait essentiel de la formation de points de vertex aux points courants de la surface de plasticité dans la détection d'une localisation de la déformation plastique en bande de cisaillement. Le deuxième axe de recherche est relatif aux instabilités de type structure (flambage, plissement ...). Ces phénomènes se manifestent le plus souvent en présence de structures minces ou relativement élancées. Dans ce cadre, nous proposons une approche basée sur l'étude de stabilité des évolutions quasi statiques. Pour des matériaux sensibles à la vitesse de déformation, l'absence de bifurcation ainsi que d'états d'équilibre nous amènent naturellement à poser le problème comme celui de la stabilité de trajectoires quasi statiques. Pour les matériaux élasto-plastiques, cette approche est justifiée par le fait que bien souvent nous sommes en présence d'une évolution quasi statique, pour un trajet de chargement donné, même si chaque point de cette évolution représente un état d'équilibre. Le critère de stabilité unifié que nous proposons est donné par la positivité de la seconde variation de l'énergie totale et est valable pour des solides visco-élastiques, visco-plastiques ou élasto-plastiques. Plus récemment, nous avons étendu ce critère à des modèles à gradients. Enfin, la comparaison de ce critère avec les résultats existants relatifs au flambage plastique nous a permis de rediscuter le critère de non-bifurcation de Hill en relation avec le choix du modèle de plasticité approprié. Il est mis en évidence, de nouveau, le rôle important joué par le développement de points de vertex sur la surface de plasticité lors de la prédiction du flambage avec la théorie du module tangent élasto-plastique. Le troisième axe de recherche concerne le développement d'éléments finis de type coques volumiques. L'idée de ce nouveau concept est de combiner les avantages des éléments coques et 3D pour formuler des éléments particulièrement adaptés à la simulation des structures minces. Ces nouveaux éléments de coques volumiques ont l'avantage de prendre en compte la flexion suivant une direction choisie tout en conservant la formulation classique des éléments volumiques. Ils permettent ainsi de modéliser des structures comportant des parties minces et des zones épaisses, sans les difficultés habituelles de raccord de maillages contenant des éléments coques et 3D. Leur utilisation en emboutissage est très prometteuse, en particulier pour des problèmes où les effets dans l'épaisseur sont d'importance majeure. Étant sous-intégrés, pour améliorer leurs performances, les modes de hourglass générés sont alors efficacement stabilisés par des techniques récentes. Les modes de verrouillages (membrane et cisaillement transverse) sont éliminés par des techniques de projection pouvant se mettre sous le formalisme " Assumed Strain Method ".