Ce travail de recherche présente une stratégie de modélisation et de dimensionnement pour les assemblages boulonnés en milieu aéronautique. Dans le cas du dimensionnement de grandes structures, les bureaux d’étude se servent habituellement de techniques dites de « zoom structural ». Elles permettent de gérer les ressources numériques à disposition, avec des modèles simples et de grande taille à l’échelle globale, et des modèles plus complexes à l’échelle locale. Cette étude s‘insère dans le cadre d’un processus de zoom structural dans lequel les structures étudiées présentent un très grand nombre de fixations. Un premier travail consiste à proposer une modélisation pertinente des fixations à l’échelle globale. Dans une seconde partie, les dispersions de comportement observées dans les fixations boulonnées sont prises en compte par l’utilisation de la théorie des méconnaissances. Une technique de recalage de la modélisation complète est ensuite proposée afin de tirer partie d’éventuels essais réels ou numériques. Enfin, la stratégie proposée est appliquée sur des cas d’assemblages boulonnés de type aéronautique. Les avantages, mais aussi les limites de l’utilisation d’une telle stratégie y sont discutés. / The objective of this work is to develop a modeling strategy for assemblies involving multiple joints in an aeronautical predesign context. In the process of designing large structures, industrial engineers usually get around computational limitations by using submodeling techniques. In order to take into account phenomena which vary during the global representation of the submodeling process, we propose a modeling approach for the scattering of these phenomena based on the Lack-Of-Knowledge (LOK) theory. The LOK model associated with a relevant deterministic joint model is updated based on experimental or simulated data. The proposed strategy is discussed and illustrated by the case of multiple-joint assemblies typical of those used in aeronautics.

Ce travail de recherche présente une stratégie de modélisation et de dimensionnement pour les assemblages boulonnés en milieu aéronautique. Dans le cas du dimensionnement de grandes structures, les bureaux d'étude se servent habituellement de techniques dites de " zoom structural ". Elles permettent de gérer les ressources numériques à disposition, avec des modèles simples et de grande taille à l'échelle globale, et des modèles plus complexes à l'échelle locale. Cette étude s'insère dans le cadre d'un processus de zoom structural dans lequel les structures étudiées présentent un très grand nombre de fixations. Un premier travail consiste à proposer une modélisation pertinente des fixations à l'échelle globale. Dans une seconde partie, les dispersions de comportement observées dans les fixations boulonnées sont prises en compte par l'utilisation de la théorie des méconnaissances. Une technique de recalage de la modélisation complète est ensuite proposée afin de tirer partie d'éventuels essais réels ou numériques. Enfin, la stratégie proposée est appliquée sur des cas d'assemblages boulonnés de type aéronautique. Les avantages, mais aussi les limites de l'utilisation d'une telle stratégie y sont discutés.

Les travaux de thèse s'inscrivent dans le cadre du projet de recherche européen MAAXIMUS (More Affordable Aircraft through eXtended, Integrated and Mature nUmerical Sizing) et portent sur la simulation numérique de problèmes d'assemblages de composants en matériaux composites stratifiés. Ces assemblages sont sources de deux types de non-linéarités. D'une part, l'interface entre les pièce conduit au traitement de non-linéarités fortes (contact, frottement). D'autre part, dans les composants de l'assemblage, le comportement du matériau stratifié est complexe, du fait des nombreux phénomènes de dégradations interagissant entre eux. Ces deux aspects ont une influence forte sur la réponse globale de l'assemblage, ce qui implique la résolution de systèmes de très grandes tailles, nécessitant généralement l'utilisation de moyens de calcul parallèles.Le couplage entre ces deux problématiques nécessite donc l'utilisation d'algorithmes de calcul parallèle dédiés et robustes, à même de traiter de nombreuses non-linéarités très fortes. Pour cela, la méthode LATIN (pour LArge Time INcrement) présente de nombreux avantages, déjà mis en évidence dans le cas de calcul d'assemblages de pièces élastiques lors de travaux précédents. Le but de ces travaux est donc d'élargir le cadre de la méthode au cas des pièces au comportement endommageable et anélastique.Un dernier aspect, qui sera abordé au sein de ces travaux, traite des fortes variabilités des coefficients intervenant dans les lois non-linéaires. Par conséquent, il est nécessaire de pouvoir traiter un très grand nombre de problèmes affectés de valeurs de coefficients différents. Pour cela, la stratégie multiparamétrique, intimement liée à la méthode LATIN, doit être étendue au cas de comportements matériau non-linéaires. Elle sera alors appliquée au travers de plusieurs paramètres variables : coefficients de frottement, précharges des éléments de fixation, seuil d'endommagement des matériaux... / The presented work, within the framework of the European research project MAAXIMUS (More Affordable Aircraft through eXtended, Integrated and Mature nUmerical Sizing), is dedicated to the numerical simulation of assemblies with components made of laminated composites. These assemblies involve two types of high non-linearities. First ones are linked to the interfaces between parts (unilateral contact and friction). Second ones are linked to the constitutive material behaviour, from its initial properties to the complex evolution of degradations. These two non-linearities have a strong influence on the response of the assembly, which involves solving systems with a high number of degrees of freedom and generally requires the use of parallel computing resources.The coupling between the two sources of non-linearities requires dedicated and robust algorithms, able to run on parallel architectures and to deal with many very strong non-linearities. The efficiency of the LATIN method (LArge Time INcrement) has already been highlighted in the case of assemblies with elastic components. A first aim of this work is thus to extend the method to the case of damageable and anelastic components' behaviour.A second aim is to deal with the variability of the coefficients involved in the non-linear laws. Each set of parameters (friction coefficients, preload of fasteners, damage threshold of material laws...) requiring a given calculation, the multiparametric strategy of the LATIN method must be extended to the case of non-linear materials in order to efficiently reduce the computation time.

Dans le contexte économique et compétitif actuel, ce qui permet aux entreprises de se distinguer c'est leur réactivité à fournir des produits de plus en plus complexes sen répondant à des besoins toujours plus diversifiés. Pour faire face à cette exigence de réactivité et malgré ces contraintes, les approches d'ingénierie intégrée sont apparues pour une meilleure gestion du cycle de vie des produits. Dans ce contexte la prise en compte des variabilités et leurs interdépendances a été démontrée comme étant indispensable à l'amélioration des performances (coût, risque, qualité, …). Ce travail de thèse est consacré à l'élaboration d'un plan d'inspection optimal pour : la maîtrise de ces variabilités et leurs effets sur le comportement du produit, l'usage et les fonctionnalités de celui-ci. Pour cela un cadre méthodologique pour la co-conception des activités de contrôles de conformités et celles de suivi de fabrication a été formalisé. Ce cadre permet une prise de décision opérationnelle grâce à l'utilisation des outils opérationnels et assure la satisfaction des objectifs stratégiques (réduction des coûts, amélioration de la qualité, augmentation de la productivité, …). Suite à une synthèse approfondie de la littérature relative à la planification la qualité et de l'amélioration (conjointe) de la performance, les outils AMDEC, KC flowdown, et Activity Based Costing (ABC) ont été retenus. Cependant ces outils présentent certains manques à combler et ont fait donc l'objet de certaines modifications (adaptations) pour convenir au mieux aux besoins du cadre méthodologique proposé. L'intégration d'AMDEC et de KC flowdown a abouti à un nouvel outil nommé ACDE (Analyse de la Causalité, des Défaillances, et leurs Effets). Ainsi l'adaptation d'ABC à l'estimation de coût de la qualité par « Quality-Driven ABC » a été proposée. L'application de ces derniers a été proposée dans le cadre de la planification d'inspection « au juste nécessaire ». / In the current economic and competitive environment, what distinguish the companies is their reactivity to provide increasingly complex products with increasingly diverse needs. To achieve this reactivity despite constraints, integrated engineering approachs have been emerged to better manage and organize the product life cycle. In this context, it is essential to take into account the characteristics' variations and their interdependencies to improve the performance (cost, risk, quality ...). This thesis was devoted to the generation of an optimal inspection plan including process variations monitoring and product variations control. A methodological framework for the co-conception of product control and process monitoring activities has been formalized and validated. This framework allows operational decision-making by the intervention of operational tools and ensures the satisfaction of strategic objectives (cost reduction, improved quality, increased productivity ...). FMEA, KC flowdown, and Activity Based Costing (ABC) were selected following a comprehensive literature review on quality planning and performance improvement. However, these tools present some lack to fill up and therefore are subjected to certain adaptations in order to best fit the needs of the proposed methodological framework. Integration of FMEA and KC flowdown resulted in a new tool called ACDE. As well the adaptation of ABC to estimate the cost of quality by "Quality-Driven ABC" was proposed. The application of these latter has been proposed to limit the inspection planning to “bare essentials activities.

Les préoccupations environnementales actuelles ont conduit les constructeurs automobiles à proposer de nouvelles technologies dans le but de réduire les émissions de CO2. Parmi ces technologies, le downsizing appliqué aux moteurs turbocompressés à allumage commandé est une des solutions privilégiées, car permettant d'atteindre des points de fonctionnement fortement chargés, avec un meilleur rendement thermique. Cependant, les fortes charges favorisent l'apparition de cliquetis, un phénomène potentiellement dommageable pour le moteur et qui l'empêche d'exploiter tout son potentiel. Du fait des variabilités cycliques de combustion dans le moteur, le cliquetis, qui dépend des conditions locales dans la chambre de combustion, peut apparaître uniquement sur quelques cycles, à différents endroits et instants. Dans cette thèse, une approche par Simulation aux Grandes Échelles (SGE) a été choisi, dans le but d'étudier et d'améliorer notre compréhension du cliquetis. L'étude se base sur la SGE d'un moteur industriel, le RENAULT 1.2 TCe 115. Un premier ensemble de 30 cycles a été simulé sur un point de fonctionnement de référence, correspondant à un point cliquetant dans la base de données banc d'essais fournie par RENAULT. Les résultats de simulation ont été comparés aux résultats expérimentaux, aussi bien en termes de variabilités cycliques de combustion que de cliquetis. A la suite, un balayage d'avance à l'allumage a été simulé pour étendre la base de données LES à des points plus faiblement et plus fortement cliquetants. La base de données résultante se compose de 150 cycles de combustion, utilisés pour développer des méthodologies et outils, dans le but de mieux caractériser le cliquetis et d'approfondir sa compréhension. L'accès numérique à toute grandeur dans la chambre de combustion, combiné à la description séparée dans cette simulation entre la flamme de pré-mélange initiée par la bougie et l'auto-inflammation dans les gaz frais, ont permis de caractériser le cliquetis en se focalisant sur son origine : l'auto-inflammation. A la suite, les méthodologies et outils développés ont soutenu une analyse détaillée des mécanismes qui contrôlent l'apparition du cliquetis. En particulier, le lien entre le cliquetis et les variabilités cycliques de combustion a été exploré. Les résultats mettent notamment en évidence l'impact des variabilités cycliques, aussi bien de la vitesse de propagation que de la forme de la flamme de pré-mélange, sur le cliquetis. / The rising concerns about the environment have led car manufacturers to come up with new engine technologies, in order to reduce the impact of internal combustion engines on CO2 emissions. In this context, downsizing of turbocharged spark-ignition engines has become a commonly used technology, the advantage of which is to operate the engine under thermally more efficient high loads. However, these high loads favour the appearance of potentially damaging knock phenomena, which prevent the engine to fully exploit its potential. Because of cyclic combustion variability (CCV) in the engine, knock, which depends on the local conditions inside the combustion chamber, can appear at different locations and timings and not in all engine cycles. In this thesis, a Large-Eddy Simulation (LES) approach was selected to investigate and further improve our understanding of the appearance of knock. The study is based on the LES of a production engine, the RENAULT 1.2 TCe 115. For this engine, a set of 30 cycles was initially simulated at a single operating point, corresponding to a knocking point in the test bench database from RENAULT. The results were compared to experimental findings, both in terms of CCV and knock. Subsequently, a spark-timing sweep was simulated in order to enlarge the LES database to also include weaker and stronger knock levels. The resulting LES, which consists of 150 combustion cycles, was used to develop methodologies and tools with the objective to better characterize and understand knock. The computational access to any quantity inside the combustion chamber, together with the separate description with the present LES approach between the spark-triggered premixed flame propagation and auto-ignition, were exploited to characterize knock focusing on its source: autoignition in the fresh gases. Then, the developed methodologies and tools supported a detailed analysis of the mechanisms that control the knock onset. In particular, its link with CCV was explored. The results point out the impact of the cyclic variability in the premixed flame propagation speed and shape on knock.

La mise en évidence de la variabilité pharmacocinétique et/ou pharmacodynamique permet l'optimisation de l'utilisation des cytotoxiques. L'association des thérapies ciblées à la chimiothérapie conventionnelle peut apporter des avantages supplémentaires en termes de bénéfice thérapeutique, mais aussi peut provoquer des interactions médicamenteuses et augmenter les variabilités interindividuelles. Les interactions médicamenteuses sont maintenant connues pour avoir un impact significatif sur l'élimination des médicaments.Le premier volet de ces travaux consiste à mesurer les concentrations sanguines de deux molécules (le lapatinib et la vinorelbine) chez les patientes et à les utiliser pour la modélisation pharmacocinétique. Cette modélisation a clairement montré l'influence du poids et du nombre des plaquettes sur la pharmacocinétique de la vinorelbine permettant ainsi de diminuer les variabilités pharmacocinétiques. Elle a également mis en évidence l'influence du lapatinib sur la pharmacocinétique de la vinorelbine. Cependant, compte tenu de l'absence de groupe témoin, nous n'avons pas réussi à obtenir une significativité statistique pour ces résultats. En parallèle, aucune influence de la vinorelbine sur le comportement pharmacocinétique du lapatinib n'a été mise en évidence.Le second volet concerne la modélisation de la réponse et de la tolérance des patientes pour cette association médicamenteuse (modélisation pharmacodynamique). La neutropénie est la toxicité dose-limitante et comme le lapatinib n'est pas connu pour être myélotoxique, nous avons modélisé cette toxicité par rapport au comportement pharmacocinétique de la vinorelbine. Là encore, nous avons observé une très forte tendance quant à l'influence du lapatinib sur la myélotoxicité de la vinorelbine. Concernant l'efficacité de cette association, la modélisation de l'évolution de la croissance tumorale a mis en évidence une synergie d'action entre ces deux molécules. A notre connaissance c'est la première fois qu'une telle modélisation tant du point de vue de la tolérance que de celui de la réponse, est réalisée lors d'une association de chimiothérapie classique et d'une thérapeutique ciblée.La modélisation PK-PD de population des médicaments anticancéreux peut apporter une aide précieuse aux cliniciens dans ce domaine. Elle peut également être essentielle dans la prise de décision clinique précoce.

La mise en évidence de la variabilité pharmacocinétique et/ou pharmacodynamique permet l’optimisation de l’utilisation des cytotoxiques. L’association des thérapies ciblées à la chimiothérapie conventionnelle peut apporter des avantages supplémentaires en termes de bénéfice thérapeutique, mais aussi peut provoquer des interactions médicamenteuses et augmenter les variabilités interindividuelles. Les interactions médicamenteuses sont maintenant connues pour avoir un impact significatif sur l'élimination des médicaments.Le premier volet de ces travaux consiste à mesurer les concentrations sanguines de deux molécules (le lapatinib et la vinorelbine) chez les patientes et à les utiliser pour la modélisation pharmacocinétique. Cette modélisation a clairement montré l’influence du poids et du nombre des plaquettes sur la pharmacocinétique de la vinorelbine permettant ainsi de diminuer les variabilités pharmacocinétiques. Elle a également mis en évidence l’influence du lapatinib sur la pharmacocinétique de la vinorelbine. Cependant, compte tenu de l’absence de groupe témoin, nous n’avons pas réussi à obtenir une significativité statistique pour ces résultats. En parallèle, aucune influence de la vinorelbine sur le comportement pharmacocinétique du lapatinib n’a été mise en évidence.Le second volet concerne la modélisation de la réponse et de la tolérance des patientes pour cette association médicamenteuse (modélisation pharmacodynamique). La neutropénie est la toxicité dose-limitante et comme le lapatinib n’est pas connu pour être myélotoxique, nous avons modélisé cette toxicité par rapport au comportement pharmacocinétique de la vinorelbine. Là encore, nous avons observé une très forte tendance quant à l’influence du lapatinib sur la myélotoxicité de la vinorelbine. Concernant l’efficacité de cette association, la modélisation de l’évolution de la croissance tumorale a mis en évidence une synergie d’action entre ces deux molécules. A notre connaissance c’est la première fois qu’une telle modélisation tant du point de vue de la tolérance que de celui de la réponse, est réalisée lors d’une association de chimiothérapie classique et d’une thérapeutique ciblée.La modélisation PK-PD de population des médicaments anticancéreux peut apporter une aide précieuse aux cliniciens dans ce domaine. Elle peut également être essentielle dans la prise de décision clinique précoce. / The identification of sources of pharmacokinetic and / or pharmacodynamic variabilities, optimizes the use of cytotoxic agents. The combination of targeted therapies with conventional drugs may provide additional benefits in terms of therapeutic benefit, but also can cause drug interactions and increased variability. Drug interactions are known to have a significant impact on drug disposition.The first part of this work is to measure blood concentrations of two molecules (lapatinib and vinorelbine) in patients and to perform pharmacokinetic modeling. This modeling clearly showed the influence of weight and the platelet number on the pharmacokinetics of vinorelbine reducing thereby the pharmacokinetic variabilities. It also highlighted the influence of lapatinib on the pharmacokinetics of vinorelbine. However, due to the lack of control group, we failed to obtain statistical significance for these results. In parallel, no effect of vinorelbine on the pharmacokinetics of lapatinib has been detected.The second part concerns the modeling of the response and tolerance of patients for this drug combination (pharmacodynamic modeling). Neutropenia was detected as the dose-limiting toxicity of the combination and the fact that lapatinib is not known to be myelotoxic, we performed toxicity modeling based on the pharmacokinetics of vinorelbine. Again, we observed a strong tendency on the influence of lapatinib on the myelotoxicity of vinorelbine. Regarding the efficacy of this combination, tumor growth modeling showed a synergistic action between the two molecules. To our knowledge this is the first time that such a model, both in terms of tolerance and response, is implemented in a combination of conventional chemotherapy and targeted therapy.Population PK-PD modeling of cancer drugs can provide valuable assistance to clinicians in this field. It can also be essential in the early clinical decision making.

Les travaux de thèse s'inscrivent dans le cadre du projet de recherche européen MAAXIMUS (More Affordable Aircraft through eXtended, Integrated and Mature nUmerical Sizing) et portent sur la simulation numérique de problèmes d'assemblages de composants en matériaux composites stratifiés. Ces assemblages sont sources de deux types de non-linéarités. D'une part, l'interface entre les pièce conduit au traitement de non-linéarités fortes (contact, frottement). D'autre part, dans les composants de l'assemblage, le comportement du matériau stratifié est complexe, du fait des nombreux phénomènes de dégradations interagissant entre eux. Ces deux aspects ont une influence forte sur la réponse globale de l'assemblage, ce qui implique la résolution de systèmes de très grandes tailles, nécessitant généralement l'utilisation de moyens de calcul parallèles.Le couplage entre ces deux problématiques nécessite donc l'utilisation d'algorithmes de calcul parallèle dédiés et robustes, à même de traiter de nombreuses non-linéarités très fortes. Pour cela, la méthode LATIN (pour LArge Time INcrement) présente de nombreux avantages, déjà mis en évidence dans le cas de calcul d'assemblages de pièces élastiques lors de travaux précédents. Le but de ces travaux est donc d'élargir le cadre de la méthode au cas des pièces au comportement endommageable et anélastique.Un dernier aspect, qui sera abordé au sein de ces travaux, traite des fortes variabilités des coefficients intervenant dans les lois non-linéaires. Par conséquent, il est nécessaire de pouvoir traiter un très grand nombre de problèmes affectés de valeurs de coefficients différents. Pour cela, la stratégie multiparamétrique, intimement liée à la méthode LATIN, doit être étendue au cas de comportements matériau non-linéaires. Elle sera alors appliquée au travers de plusieurs paramètres variables : coefficients de frottement, précharges des éléments de fixation, seuil d'endommagement des matériaux...

Le sous-courant équatorial (EUC), en alimentant l'upwelling équatorial, peut avoir une forte influence sur la température de surface (SST) du Pacifique équatorial Est et donc sur la variabilité associée à El Niño. L'EUC et les cellules de circulation méridienne (shallow subtropical/tropical overturning cells, STCs/TCs) l'alimentant sont étudiés en combinant données in situ et modélisation. Les trajectoires de masses d'eau sont calculées dans des simulations réalistes (le modèle OPA forcé par les vents des réanalyses NCEP sur 1948-1999 ou des satellites ERS sur 1992-1999). Leur analyse met en évidence des cheminements des masses d'eaux propres aux évènements El Niño-La Niña de 1997-1998, avec des recharges et décharges de la bande équatoriale complexes et asymétriques. Cette analyse montre aussi l'apport d'eaux froides par les STCs et l'EUC lors de la brusque transition vers La Niña en mai 1998. Les données de courant et de température des mouillages TAO/TRITON le long de l'équateur à 170°W, 140°W et 110°W sont méthodiquement bouchées sur 1980-2002. On montre que des séries continues du débit, de la température, de la profondeur et de l'énergie cinétique de l'EUC sur toute son extension méridienne peuvent alors être construites. Leur analyse révèle que la forte variabilité interannuelle du débit de l'EUC est une réponse linéaire et quasi-stationnaire à la tension de vent zonale équatoriale intégrée zonalement dans le Pacifique Ouest et central. La température de l'EUC, indispensable pour l'estimation du transport de chaleur, varie elle linéairement avec la différence des profondeurs de la thermocline et de l'EUC dans le Pacifique central. Le modèle numérique, validé entre autre à l'aide des séries de l'EUC, est utilisé pour étudier sur 1951-1999 la circulation équatoriale associée à l'EUC: la convergence dans la pycnocline, l'upwelling équatorial et la divergence en surface à 5°N et 5°S. Leurs variabilités en débit sont quasi-égales à celle de l'EUC, qui est donc un bon indicateur de la force des STCs. Ces variabilités sont principalement causées par la tension de vent zonale intégrée zonalement sur tout le bassin, en accord avec des théories linéaires. Des déphasages avec la SST équatoriale, notamment l'avance de 5 mois de l'upwelling et de la divergence sur la SST, révèlent des relations de cause à effet très intéressantes, confirmées par les bilans de chaleur. La différence de température entre la divergence et la convergence a des variations interannuelles et à plus long-terme égales à celles de la SST équatoriale. Les conséquences sur les bilans et échanges de masse et de chaleur dans la bande équatoriale sont ensuite quantifiées. La variabilité du transport de chaleur méridien associé à la convergence/divergence est due aussi bien aux variations de débit que de température de la convergence et de la divergence. Ainsi, pendant un événement El Niño, la baisse des débits aura tendance à réchauffer la bande équatoriale (recharge), alors que l'augmentation de la différence entre les températures de la divergence et de la convergence aura l'effet contraire (décharge). Pour le Pacifique Est, les variations interannuelles du débit dominent celles de la température de l'EUC pour le transport de chaleur de l'EUC. Les liens avec les théories d'El Niño et sa variabilité décennale sont discutés.