La vitesse critique de flottement est un facteur essentiel à la conception aéronautique car elle caractérise le régime de vol au-delà duquel l’aéronef risque de subir un mécanisme de ruine. L’objectif de cette thèse est d’étudier l’impact des incertitudes d’origines aléatoires et épistémiques sur la limite de stabilité linéaire pour des configurations aéroélastiques idéalisées. Dans un premier temps, un problème de propagation directe d’incertitudes aléatoires relatives à des paramètres de fabrication d’une aile en forme de plaque en matériau composite stratifié a été considéré. La représentation du matériau par la méthode polaire lève la contrainte de grande dimensionnalité du problème stochastique initial et permet l’utilisation du Chaos Polynômial. Cependant, la corrélation introduite par cette paramétrisation nécessite une adaptation de la base polynômiale. Enfin, un algorithme d’apprentissage automatique a été employé pour traiter des discontinuités dans le comportement modal des instabilités aéroélastiques. Le second volet de la thèse concerne la quantification d’incertitudes de modélisation de caractère épistémique qui sont introduites au niveau de l’opérateur aérodynamique. Ces travaux, menés à partir d’un formalisme Bayésien, permettent non seulement d’établir des probabilités de modèle, mais aussi de calibrer les coefficients des modèles dans un contexte stochastique afin d’obtenir des prédictions robustes pour la vitesse critique. Enfin, une étude combinée des deux types d’incertitude permet d’améliorer le processus de calibration. / The critical flutter velocity is an essential factor in aeronautic design because it caracterises the flight envelope outside which the aircraft risks to be destroyed. The goal of this thesis is the study of the impact of uncertainties of aleatory and epistemic origin on the linear stability limit of idealised aeroelastic configurations. First, a direct propagation problem of aleatory uncertainties related to manufacturing parameters of a rectangular plate wing made of a laminated composite material was considered. The representation of the material through the polar method alleviates the constraint of the high number of dimensions of the initial stochastic problem, which allows the use of polynomial chaos. However, the correlation which is introduced by this parametrisation requires an adaption of the polynomial basis. Finally, a machine learning algorithm is employed for the treatment of discontinuities in the modal behaviour of the aeroelastic instabilities. The second part of the thesis is about the quantification of modelling uncertainties of epistemic nature which are introduced in the aerodynamic operator. This work, which is conducted based on a Bayesian formalism, allows not only to establish model probabilities, but also to calibrate the model coefficients in a stochastic context in order to obtain robust predictions for the critical velocity. Finally, a combined study of the two types of uncertainty allows to improve the calibration process.

Dans ce mémoire, nous étudions l'exposition aux ondes électromagnétiques l'échelle d'une population. Cette étude est réalisée pour une exposition à une onde plane orientée frontalement sur des modèles anatomiques du corps humains (fantômes) avec une puissance incidente de 1W/m2 et une fréquence de 2.1 GHz. Pour ce faire, l'idée est de construire un modèle du WBSAR (Whole Body averaged Specific Absorption Rate) en fonction de lamorphologie. Les facteurs morphologiques influençant le WBSAR ont été identifiés à l'aide des fantômes existant dont le nombre se limite à 18 et des modèles de régressions. Cette analyse préliminaire montre que les facteurs morphologiques externes (taille, poids ...) seuls ne suffisent pas pour construire un tel modèle, mais l'introduction des facteurs morphologiques internes (muscles, graisse ...) est nécessaire. L'absence de données statistiques sur les facteurs morphologiques internes de populations nous a conduits à intégrer des connaissances apriori sur ces facteursafin d'étudier l'exposition d'une population donnée. Des lois paramétriques usuelles et des mélanges de gaussiennes sont utilisés pour modéliser ces facteurs internes afin d'étudier leur influence sur le quantile du WBSAR à 95 %.L'utilisation des fantômes homogènes où le fantôme homogène est obtenu en remplaçant tous les tissus internes par un unique tissu équivalent. Ces fantômes homogènes permettent de s'affranchir de l'influence de la morphologie interne dans un modèle du WBSAR et facilitent également leur déformation par une technique de morphing. Ainsi nous avons pu enrichir la base de fantômes existant. Pour donner une estimation du quantile du WBSAR à 95 %, nous avons mis en place un plan d'expériences séquentiel qui repose sur un modèle paramétrique du WBSAR et l'inférence bayésienne et qui permet de raffiner la région à 95 %. Cette approche nous a permis de trouver le quantile à 95 % pour une population et un exemple de fantôme correspondant. Cependant, cette approche ne permet pas d'obtenir une estimation de toute la distribution du WBSAR. Afin d'obtenir une estimation de toute cette distribution, nous avons construit une surface de réponse en utilisant les polynômes de chaos. Dans l'objectif d'obtenirdes résultats cohérents nous avons effectué un changement de variables permettant de traduire les connaissances physiques dans cette surface de réponse.

Dans ce mémoire, nous étudions l’exposition aux ondes électromagnétiques l’échelle d’une population. Cette étude est réalisée pour une exposition à une onde plane orientée frontalement sur des modèles anatomiques du corps humains (fantômes) avec une puissance incidente de 1W/m2 et une fréquence de 2.1 GHz. Pour ce faire, l’idée est de construire un modèle du WBSAR (Whole Body averaged Specific Absorption Rate) en fonction de lamorphologie. Les facteurs morphologiques influençant le WBSAR ont été identifiés à l’aide des fantômes existant dont le nombre se limite à 18 et des modèles de régressions. Cette analyse préliminaire montre que les facteurs morphologiques externes (taille, poids ...) seuls ne suffisent pas pour construire un tel modèle, mais l’introduction des facteurs morphologiques internes (muscles, graisse ...) est nécessaire. L’absence de données statistiques sur les facteurs morphologiques internes de populations nous a conduits à intégrer des connaissances apriori sur ces facteursafin d’étudier l’exposition d’une population donnée. Des lois paramétriques usuelles et des mélanges de gaussiennes sont utilisés pour modéliser ces facteurs internes afin d’étudier leur influence sur le quantile du WBSAR à 95 %.L’utilisation des fantômes homogènes où le fantôme homogène est obtenu en remplaçant tous les tissus internes par un unique tissu équivalent. Ces fantômes homogènes permettent de s’affranchir de l’influence de la morphologie interne dans un modèle du WBSAR et facilitent également leur déformation par une technique de morphing. Ainsi nous avons pu enrichir la base de fantômes existant. Pour donner une estimation du quantile du WBSAR à 95 %, nous avons mis en place un plan d’expériences séquentiel qui repose sur un modèle paramétrique du WBSAR et l’inférence bayésienne et qui permet de raffiner la région à 95 %. Cette approche nous a permis de trouver le quantile à 95 % pour une population et un exemple de fantôme correspondant. Cependant, cette approche ne permet pas d’obtenir une estimation de toute la distribution du WBSAR. Afin d’obtenir une estimation de toute cette distribution, nous avons construit une surface de réponse en utilisant les polynômes de chaos. Dans l’objectif d’obtenirdes résultats cohérents nous avons effectué un changement de variables permettant de traduire les connaissances physiques dans cette surface de réponse. / In this report, a study of the exposure to electromagnetic waves at the level of population is proposed. This study is conducted in the case of a plane wave oriented frontally on anatomical models (phantoms) with an incident power of 1W/m2 and at the frequency of 2.1 GHz. For this purpose, the idea is to build a surrogate model of the WBSAR (Whole Body averaged Specific Absorption Rate) as a function of morphology. An identificationof morphological factor is performed using the limited number of phantoms (18phantoms)and different regressions. This analysis shows that to use only external morphological factors (height, weight ...) is not sufficient to build such surrogate model and that internal morphological factors (muscles, fat ...) are also important. The absence of statistical data on internal morphology conducts us to introduce some prior knowledge on these internal factors to study the exposure for a given population. Some parametric laws and Gaussian mixture are used to study their influence on the quantile of the WBSAR at 95 The lack of data on internal morphology led us to consider homogeneous phantoms instead of the heterogeneous one. In homogeneous phantoms, all the internal tissues are substituted by an equivalent tissue. These homogeneous phantoms allow us to overcome the influence of internal morphology for a WBSAR model and to use a morphing technique to increase the database of phantoms. To give an estimation of the WBSAR quantile at 95 %, a sequential experiment design is used. This approach is based on a parametric model and Bayesian inference. This experiment design allows us to estimate the quantile at 95 for a given population and an example of one corresponding phantom. Nevertheless, this approach does not allow one to estimate the whole distribution of the WBSAR for a given population. To give an estimate of this distribution a surface of response is established using polynomial chaos. This approach gives incoherent results with the physical phenomenon. To obtain coherent results, physical knowledge is introduced by variables changing in this responsesurface.

Le problème de la dimension stochastique élevée est récurrent dans les analyses probabilistes des structures. Il correspond à l'augmentation exponentielle du nombre d'évaluations du modèle mécanique lorsque le nombre de paramètres incertains est élevé. Afin de pallier cette difficulté, nous avons proposé dans cette thèse, une approche à deux étapes. La première consiste à déterminer la dimension stochastique efficace, en se basant sur une hiérarchisation des paramètres incertains en utilisant les méthodes de criblage. Une fois les paramètres prépondérants sur la variabilité de la réponse du modèle identifiés, ils sont modélisés par des variables aléatoires et le reste des paramètres est fixé à leurs valeurs moyennes respectives, dans le calcul stochastique proprement dit. Cette tâche fut la deuxième étape de l'approche proposée, dans laquelle la méthode de décomposition de la dimension est utilisée pour caractériser l'aléa de la réponse du modèle, par l'estimation des moments statistiques et la construction de la densité de probabilité. Cette approche permet d'économiser jusqu'à 90% du temps de calcul demandé par les méthodes de calcul stochastique classiques. Elle est ensuite utilisée dans l'évaluation de l'intégrité d'une toiture à ossature bois d'une habitation individuelle installée sur un site d'aléa sismique fort. Dans ce contexte, l'analyse du comportement de la structure est basée sur un modèle éléments finis, dans lequel les assemblages en bois sont modélisés par une loi anisotrope avec hystérésis et l'action sismique est représentée par huit accélérogrammes naturels fournis par le BRGM. Ces accélérogrammes permettent de représenter différents types de sols selon en se référant à la classification de l'Eurocode 8. La défaillance de la toiture est définie par l'atteinte de l'endommagement, enregistré dans les assemblages situés sur les éléments de contreventement et les éléments d'anti-flambement, d'un niveau critique fixé à l'aide des résultats des essais. Des analyses déterministes du modèle éléments finis ont montré que la toiture résiste à l'aléa sismique de la ville du Moule en Guadeloupe. Les analyses probabilistes ont montré que parmi les 134 variables aléatoires représentant l'aléa dans le comportement non linéaire des assemblages, 15 seulement contribuent effectivement à la variabilité de la réponse mécanique ce qui a permis de réduire la dimension stochastique dans le calcul des moments statistiques. En s'appuyant sur les estimations de la moyenne et de l'écart-type on a montré que la variabilité de l'endommagement dans les assemblages situés dans les éléments de contreventement est plus importante que celle de l'endommagement sur les assemblages situés sur les éléments d'anti-flambement. De plus, elle est plus significative pour les signaux les plus nocifs sur la structure.

Le problème de la dimension stochastique élevée est récurrent dans les analyses probabilistes des structures. Il correspond à l’augmentation exponentielle du nombre d’évaluations du modèle mécanique lorsque le nombre de paramètres incertains est élevé. Afin de pallier cette difficulté, nous avons proposé dans cette thèse, une approche à deux étapes. La première consiste à déterminer la dimension stochastique efficace, en se basant sur une hiérarchisation des paramètres incertains en utilisant les méthodes de criblage. Une fois les paramètres prépondérants sur la variabilité de la réponse du modèle identifiés, ils sont modélisés par des variables aléatoires et le reste des paramètres est fixé à leurs valeurs moyennes respectives, dans le calcul stochastique proprement dit. Cette tâche fut la deuxième étape de l’approche proposée, dans laquelle la méthode de décomposition de la dimension est utilisée pour caractériser l’aléa de la réponse du modèle, par l’estimation des moments statistiques et la construction de la densité de probabilité. Cette approche permet d’économiser jusqu’à 90% du temps de calcul demandé par les méthodes de calcul stochastique classiques. Elle est ensuite utilisée dans l’évaluation de l’intégrité d’une toiture à ossature bois d’une habitation individuelle installée sur un site d’aléa sismique fort. Dans ce contexte, l’analyse du comportement de la structure est basée sur un modèle éléments finis, dans lequel les assemblages en bois sont modélisés par une loi anisotrope avec hystérésis et l’action sismique est représentée par huit accélérogrammes naturels fournis par le BRGM. Ces accélérogrammes permettent de représenter différents types de sols selon en se référant à la classification de l’Eurocode 8. La défaillance de la toiture est définie par l’atteinte de l’endommagement, enregistré dans les assemblages situés sur les éléments de contreventement et les éléments d’anti-flambement, d’un niveau critique fixé à l’aide des résultats des essais. Des analyses déterministes du modèle éléments finis ont montré que la toiture résiste à l’aléa sismique de la ville du Moule en Guadeloupe. Les analyses probabilistes ont montré que parmi les 134 variables aléatoires représentant l’aléa dans le comportement non linéaire des assemblages, 15 seulement contribuent effectivement à la variabilité de la réponse mécanique ce qui a permis de réduire la dimension stochastique dans le calcul des moments statistiques. En s’appuyant sur les estimations de la moyenne et de l’écart-type on a montré que la variabilité de l’endommagement dans les assemblages situés dans les éléments de contreventement est plus importante que celle de l’endommagement sur les assemblages situés sur les éléments d’anti-flambement. De plus, elle est plus significative pour les signaux les plus nocifs sur la structure. / The problem of the curse of dimensionality is frequently encountered in practical applications. It can be defined as the significant increase of the number of mechanical model calls with the number of uncertain parameters. To overcome this difficulty, a two-steps stochastic approach has been developed in this work. The first step of this approach consists in calculating the stochastic effective dimension by the means of Morris screening method. Once the most significant uncertain parameters on the variability of the mechanical responses are identified, they are modeled as random variables and the remaining parameters are fixed to their respective mean values. This allows us to reduce significantly the stochastic dimension of the problem in the second step of the approach where the decomposition method is used to estimate the statistical characteristics of the mechanical responses. The efficiency and the accuracy of this approach are evaluated through an academic problem dealing with the assessment of the integrity of a three-span five-story frame structure subjected to horizontal loads. We have demonstrate that we can reduce about 90% of the computation time required by the classical stochastic methods. Then, the proposed approach is used to the analysis of the integrity of timber roofs under seismic loading. The behaviour of this structure is described through a finite element model where the timber joints are modeled by anisotropic hysteresis law, and the seismic action is represented by eight real earthquake ground motion records. These accelerograms provided by the French institution involved in geosciences BRGM allow us to take into account different soil types according to the classification provided by the europeen design code dealing with seismic events Eurocode 8. The failure of timber roofs is reached when the damage levels in the timbers joints localized on the buckling and bracing members reach the critical value. It is shown, through a deterministic analysis, that the structure resists the seismic hazard representing the city of Le Moule in Guadeloupe. The stochastic analysis has shown that, among the 134 random variables representing the uncertainty in the nonlinear behaviour of the timber joints, only 15 have a significant effect on the variability of the structural response, which allow us to reduce the stochastic dimension in the computation of the statistical moments. According to the estimates of the mean and the standard deviation, we have shown that the variability of bracing members damage is greater than the variability of buckling members damage. Moreover, the variability of the bracing members damage is more significant for the earthquake ground motion records having the lowest collapse PGA.