On s'intéresse ici au calcul des déformations post-sismiques, dues principalement au phénomène de relaxation dans l'asthénosphère visqueuse. Ces simulations font intervenir un comportement viscoélastique et des éléments de zone cohésive pour modéliser le glissement sur la faille. La prise en compte de la géométrie particulière de cette zone, ainsi que la nécessité de couvrir une zone géographique allant de la Chine à l'Australie conduit à des calculs par éléments finis de très grande taille, que l'on aborde par décomposition de domaines sur calculateur parallèle.