Parallel computing features of a multiscale strategy with space and time homogenization
Résumé
The behavior of large structures often results from complex phenomena taking place at a fine scale, with small length of variations both in space and time. For instance local cracking or local buckling have a global effect on large structures. In order to predict the behavior of such structures, sophisticated models have been developed to describe the material at a very fine scale compared to the scale of the structure. Corresponding numerical simulations require the use of fine discretizations over both space and time domains, which often leads to systems with a very large number of degrees of freedom whose calculation cost is generally prohibitive. Recently, a new multiscale computational strategy has been proposed for the analysis of such problems [1,2,3]. This strategy is based on the LArge Time INcrement method (LATIN method [4]), which enables one to work globally over the time-space domain. It can be seen as a parallel mixed domain decomposition strategy in space, including automatic space and time homogenization, and for which no periodicity condition is needed, in opposition with standard homogenization techniques. The first idea consists in splitting the structure into substructures and interfaces, and dividing the studied time interval into subintervals. Interface fields are separated into a macro part and a micro complement. Beacuse of the Saint-Venant's principle, an adapted choice of macro quantities, combined with the resolution of a global macro problem, ensures that micro quantities have only local effects. This choice provides to the method the numerical scalability with respect to space variables. A technique to extend this property to time aspects is in progress. The second idea consists in using an iterative technique involving the resolution of a set of independant problems on the refined "micro" scale and the resolution of a number of problems on the homogenized "macro" scale. A new approximation technique for the resolution of micro problems is presented as an improvement over the version introduced in [1]. It leads to drastic saving in terms of computation cost. A three-dimensional numerical example, with (visco)plastic material and unilateral contact with friction is provided to show the performance of the strategy combined with the new approximation technique, and to highlight the parallel properties.