L'identification des paramètres influants pour le comportementdes liaisons par contact (coefficients de frottement, jeux, ...) n'est pas chose aisée. La mesure de tels paramètres est biensouvent fortement entachée d'erreurs très importantes.Les phénomènes mis en jeux (contact et frottement)sont fortement non-linéaires. La construction par simulationsnumériques de réponses moyennes de l'assemblage ou bien desurfaces de réponse nécessite un grand nombre de calcul trèscoûteux. Des techniques efficaces pour l'obtention rapide de telles informations sont proposées dans cet article.Pour décrire les assemblages nous utilisons une décomposition detype sous-structures et interfaces. L'algorithme de résolution utilisé est basé sur la méthode LATIN.Pour construire à moindre frais la surface de réponse,nous utilisons un résultat de calcul obtenu préalablementpour un jeu de paramètres comme initialisation pour un calculavec un autre jeu de paramètres.Dans un second temps, nous proposons une méthode de constructiondes caractéristiques stochastiques de la réponse de l'assemblage,par projection de la solution dans une base de type chaospolynomial.