A multiscale method with patches for the propagation of localized uncertainties in structural dynamics

Florent Pled (1) , Mathilde Chevreuil (2) , Anthony Nouy (2)

1 MSME - Laboratoire de Modélisation et Simulation Multi Echelle
2 GeM - Institut de Recherche en Génie Civil et Mécanique

Uncertainty quantification in structural computation has gained an ever increasing interest in the scientific community during the last decades. The propagation of localized uncertainties through stochastic computational models is critical for the design and analysis of mechanical structures and components. The input uncertainties in the computational model may stem from intrinsic variabilities in the material properties or epistemic variabilities resulting from partial or limited information about the geometry, the boundary or initial conditions. Classical monoscale approaches require either local refinement or enrichment techniques to cope with the high complexity of multiscale stochastic problems. On the contrary, multiscale coupling approaches based on substructuring, domain decomposition or multigrid methods have been designed to solve such intractable high-dimensional problems. A multiscale method based on patches has been recently proposed in [1] for the solution of linear stochastic multiscale models (with localized uncertainties) and extended to a broader class of non-linear stochastic multiscale models (with localized uncertainties and non-linearities) in [2]. It relies on a partition of the domain into several subdomains of interest (called patches) containing the different sources of uncertainties and possible non-linearities, and a complementary subdomain. The multiscale solution is then computed using a global-local iterative algorithm that involves the solution of a sequence of linear global problems (with possibly deterministic operators and uncertain right-hand sides) over a deterministic domain, and (non-)linear local problems (with uncertain operators and right-hand sides) over the patches. In the present work, the method is extended to linear structural dynamics problems featuring localized uncertainties. The convergence of the iterative algorithm is analyzed. The proposed multiscale approach offers the possibility to use different global and local computational models with suitable discretizations as well as stand-alone global and local solvers. At the local level, the stochastic problems are solved using sampling-based approaches combined with adaptive sparse approximation methods [3] in order to efficiently derive sparse representations of high-dimensional stochastic local solutions with controlled accuracy. The performances of the multiscale domain decomposition method are illustrated through numerical experiments carried out on a linear vibration stochastic problem with localized random material heterogeneities. References [1] M. Chevreuil, A. Nouy, and E. Safatly. A multiscale method with patch for the solution of stochastic partial differential equations with localized uncertainties. Computer Methods in Applied Mechanics and Engineering, 255(0):255–274, 2013. [2] A. Nouy and F. Pled. A multiscale method for semi-linear elliptic equations with localized uncertainties and non-linearities. ESAIM:M2AN (submitted), 2017. [3] A. Chkifa, A. Cohen, and C. Schwab. Breaking the curse of dimensionality in sparse polynomial approximation of parametric PDEs. Journal de Mathématiques Pures et Appliquées, 103(2):400–428, 2015.

Domaines: Mécanique des structures [physics.class-ph]

Florent Pled, Mathilde Chevreuil, Anthony Nouy. A multiscale method with patches for the propagation of localized uncertainties in structural dynamics. 2nd International Conference on Uncertainty Quantification in Computational Sciences and Engineering (UNCECOMP 2017), Jun 2017, Rhodes, Greece. ⟨hal-01541080⟩