Bien que de plus en plus employée au sein des bureaux d'études dans le cadre de la conception mécanique, l'optimisation reste à l'heure actuelle encore relativement peu utilisée dans le cadre des assemblages de structures. La résolution mécanique de ce type de problèmes nécessite la mise en oeuvre de méthodes numériques de résolution capables de prendre en compte des non-linéarités de type frottement et contact. Ces méthodes sont, en raison du temps de calcul, généralement incompatibles avec la recherche d'un optimum global nécessitant un nombre important de résolutions. C'est pour pallier à ce problème que ce travail propose l'emploi d'une approche d'optimisation bi-niveaux de modèles faisant intervenir deux outils: (1) la stratégie multiparamétrique, basée sur la méthode LaTIn, assure d'importantes réductions des temps de calcul associés aux multiples résolutions du problème mécanique et (2) un métamodèle de type cokrigeage construit à partir d'un nombre limité de réponses et gradients calculés par le solveur mécanique. Ce métamodèle est alors capable de fournir à coût extrêmement faible des réponses approchées de la fonction objectif. Une optimisation globale est ensuite réalisée sur ce dernier, assurant l'obtention d'un optimum global. Le métamodèle de cokrigeage est étudié en détail sur des exemples analytiques et mécaniques comportant divers nombres de paramètres. Par ailleurs, une étude complète de l'emploi de la stratégie multiparamétrique est proposée et de nombreux exemples d'assemblages sont considérés, permettant ainsi d'illustrer les performances significatives de la procédure d'optimisation proposée.