Concurrent multi-domain simulations in linear structural dynamics using multiple grid multiple time-scale (MGMT) method
Work presented in this paper describes a general algorithm and its finite element implementation for performing concurrent multiple sub-domain simulations in linear structural dynamics. Using this approach one can solve problems in which the domain under analysis can be selectively discretized spatially and temporally, hence allowing the user to obtain a desired level of accuracy in critical regions while improving computational efficiency globally. The mathematical background for this approach is largely based upon the fundamental principles of domain decomposition methods (DDM) and Lagrange multipliers, used to obtain coupled equations of motion for distinct regions of a continuous domain. These methods when combined together systematically yield constraint forces that not only ensure conservation of energy but also enforce continuity of velocities across sub-domain interfaces. Multiple grid (MG) connections between non-conforming sub-domains are modeled using mortar elements whereas coupled multiple time-scale (MT) equations are derived for the classical Newmark integration algorithm (and its constituents). Fully discretized equations of motion for component sub-domains, augmented with an interface condition are then solved using block elimination method and Crout factorization. A proof of stability is provided using Energy method and overall efficiency, accuracy and stability of multiple sub-domain coupling is evaluated using a series of numerical examples. Primary observations are made for the evolution and distribution of kinematic quantities and structural wave propagation across connecting sub-domain. Numerical stability is verified by ensuring energy balance at global as well as component sub-domain level. Discussed examples highlight the greatest advantage of MGMT method; which is high simulation speedups (at the cost of reasonably small errors).