L'identification des paramètres influants pour le comportement des liaisons par contact (coefficients de frottement, jeux, ...) n'est pas chose aisée. La mesure de tels paramètres est bien souvent 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 simulations numériques de réponses moyennes de l'assemblage ou bien de surfaces de réponse nécessite un grand nombre de calcul très coû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 de type 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éalablement pour un jeu de paramètres comme initialisation pour un calcul avec un autre jeu de paramètres. Dans un second temps, nous proposons une méthode de construction des caractéristiques stochastiques de la réponse de l'assemblage, par projection de la solution dans une base de type chaos polynomial.