Différentiation automatique de codes mécaniques : application à l'analyse de sensibilité des tôles sandwich aux paramètres de modélisation

In engineering, for a better understanding of the mechanical behavior of a structure submitted to some perturbation of the modeling parameters, one often proceed to a sensitivity analysis. This provides quantitative and qualitative information on the behavior of the model under study and gives access to gradients that may be used in identification and optimization methods. In this thesis, we demonstrate that this information may be obtained at a low development effort by applying an Automatic Differentiation (AD) tool to the computer code that implements the model. We adapt the AD techniques to the Asymptotic Numerical Method (ANM), in its Diamant version for sensitivity computations of numerical solutions of nonlinear problems discretized through a finite element method. We discuss in a generic manner both the theoretical aspects and the implementation of several algorithms written in Matlab. Applications are concerned with sandwich beams and sandwich plates in both the static and dynamic (free vibration) cases. Sensitivities are computed with respect to geometric and mechanical parameters, and with respect to elementary stiffness matrix. The generality of our developments allows to take into account several viscoelastic laws with no additional effort. Three kinds of viscoelastic models are studied: constant complex modulus, low damping and higher damping. In comparison with the finite difference approximation often used in mechanics, our approach provides more accurate results for the sensitivity of the structure response to a perturbation of the modeling parameters. It also allows a reduction of the computation effortEn ingénierie, pour mieux comprendre le comportement mécanique d'une structure soumise à une certaine perturbation des paramètres de conception, on procède souvent à une analyse de sensibilité. Celle-ci fournit des informations quantitatives et qualitatives sur le comportement du modèle étudié et offre un accès aux gradients utilisables dans ces méthodes d'identification et d'optimisation. Dans cette thèse, nous démontrons que ces informations peuvent être obtenues à coût de développement faible en appliquant un outil de Différentiation Automatique (DA) au code informatique qui implémente le modèle. Nous adaptons la technique DA à la méthode asymptotique numérique, dans sa version Diamant, pour le calcul de la sensibilité des solutions numériques de problèmes non-linéaires discrétisés par la méthode des éléments finis. Nous discutons de manière générique à la fois les aspects théoriques et l'implémentation de plusieurs algorithmes écrits en Matlab. Les applications concernent des poutres et des plaques sandwich dans les cas statiques et dynamique (vibration libre). Les sensibilités sont calculées par rapport aux paramètres géométriques, mécanique et par rapport à des matrices de rigidité élémentaires. La généralité de nos développements permet de prendre en compte plusieurs lois viscoélastiques sans effort supplémentaire. Trois types de modèles viscoélastiques sont étudiés : module complexe constant, faible amortissement et fort amortissement. Comparée à l'approximation par différences finis souvent utilisée en mécanique, notre approche fournit des résultats plus précis pour la sensibilité de la réponse d'une structure lorsque les paramètres de conception sont perturbés. Elle permet aussi de réduire le temps de calcu