Using parametric model order reduction for inverse analysis of large nonlinear cardiac simulations

Predictive high-fidelity finite element simulations of human cardiac mechanics commonly require a large number of structural degrees of freedom. Additionally, these models are often coupled with lumped-parameter models of hemodynamics. High computational demands, however, slow down model calibration and therefore limit the use of cardiac simulations in clinical practice. As cardiac models rely on several patient-specific parameters, just one solution corresponding to one specific parameter set does not at all meet clinical demands. Moreover, while solving the nonlinear problem, 90% of the computation time is spent solving linear systems of equations. We propose to reduce the structural dimension of a monolithically coupled structure-Windkessel system by projection onto a lower-dimensional subspace. We obtain a good approximation of the displacement field as well as of key scalar cardiac outputs even with very few reduced degrees of freedom, while achieving considerable speedups. For subspace generation, we use proper orthogonal decomposition of displacement snapshots. Following a brief comparison of subspace interpolation methods, we demonstrate how projection-based model order reduction can be easily integrated into a gradient-based optimization. We demonstrate the performance of our method in a real-world multivariate inverse analysis scenario. Using the presented projection-based model order reduction approach can significantly speed up model personalization and could be used for many-query tasks in a clinical setting

A nonsmooth frictional contact formulation for multibody system dynamics

We present a new node-to-face frictional contact element for the simulation of the nonsmooth dynamics of systems composed of rigid and flexible bodies connected by kinematic joints. The equations of motion are integrated using a nonsmooth generalized-α time integration scheme and the frictional contact problem is formulated using a mixed approach, based on an augmented Lagrangian technique and a Coulomb friction law. The numerical results are independent of any user-defined penalty parameter for the normal or tangential component of the forces and, the bilateral and the unilateral constraints are exactly fulfilled both at position and velocity levels. Finally, the robustness and the performance of the proposed algorithm are demonstrated by solving several numerical examples of nonsmooth mechanical systems involving frictional contact.Fil: Galvez, Javier. Université de Liège; BélgicaFil: Cavalieri, Federico José. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Centro de Investigaciones en Métodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones en Métodos Computacionales; ArgentinaFil: Cosimo, Alejandro. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Centro de Investigaciones en Métodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones en Métodos Computacionales; Argentina. Université de Liège; BélgicaFil: Brüls, Olivier. Université de Liège; BélgicaFil: Cardona, Alberto. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Centro de Investigaciones en Métodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones en Métodos Computacionales; Argentin

Hybrid simulation of a structure to tsunami loading

A new hybrid simulation technique has been developed to assess the behavior of a structure under hydrodynamic loading. It integrates the computational fluid dynamics and structural hybrid simulation and couples the fluid loading and structure response at each simulation step. The conventional displacement-based and recently developed force-based hybrid simulation approaches are adopted in the structural analysis. The concept, procedure, and required components of the proposed hybrid simulation are introduced in this paper. The proposed hybrid simulation has been numerically and physically tested in case of a coastal building impacted by a tsunami wave. It is demonstrated that the force error in the displacement-based approach is significantly larger than that in the force-based approach. The force-based approach allows for a more realistic and reliable structural assessment under tsunami loading

Damage Analysis of Reinforced Concrete Structures with Substandard Detailing

The goal of this study is to investigate seismic behaviour of existing R/C buildings designed and constructed in accordance with standards that do not meet current seismic code requirements. In these structures, not only flexure, but also shear and bond-slip deformation mechanisms need to be considered, both separately and in combination. To serve this goal, a finite element model is developed for inelastic seismic analysis of complete planar R/C frames. The proposed finite element is able to capture gradual spread of inelastic flexural and shear deformations as well as their interaction in the end regions of R/C members. Additionally, it is capable of predicting shear failures caused by degradation of shear strength in the plastic hinges of R/C elements, as well as pullout failures caused by inadequate anchorage of the reinforcement in the joint regions. The finite element is fully implemented in the general inelastic finite element code IDARC2D and it is verified against experimental results involving individual column and plane frame specimens with nonductile detailing. It is shown that, in all cases, satisfactory correlation is established between the model predictions and the experimental evidence. Finally, parametric studies are conducted to illustrate the significance of each deformation mechanism on the seismic response of the specimens under investigation. It is concluded, that all deformation mechanisms, as well as their interaction, should be taken into consideration in order to predict reliably seismic damage of R/C structures with substandard detailing

Serviceability performance of steel-concrete composite beams

YesFor composite beams with low degrees of shear connection, additional deflections occur due to slip in the shear connectors, which can be significant for beams with low degrees of shear connection. A design formula is presented for the effective stiffness of composite beams taking account of the stiffness of the shear connectors, which is compared to measured deflections of 6 symmetric beams and an 11m span composite beam of asymmetric profile. It is shown that the comparison is good when using a shear connector stiffness of 70 kN/mm for single shear connectors and 100 kN/mm for pairs of shear connectors per deck rib. Results of push tests on a range of deck profiles confirm these initial elastic stiffnesses. To ensure that the slip at the serviceability limit state does not lead to permanent deformations of the beam, it is proposed that the minimum degree of shear connection should not fall below 30% for un-propped beams and 40% for propped beams of symmetric cross-section.European Commissio

Earthquake-reactivated landslide scenarios in Southern Italy based on spectral-matching input analysis

n/

On the coupling of local 3D solutions and global 2D shell theory in structural mechanics

Most of mechanical systems and complex structures exhibit plate and shell components. Therefore, 2D simulation, based on plate and shell theory, appears as an appealing choice in structural analysis as it allows reducing the computational complexity. Nevertheless, this 2D framework fails for capturing rich physics compromising the usual hypotheses considered when deriving standard plate and shell theories. To circumvent, or at least alleviate this issue, authors proposed in their former works an in-plane-out-of-plane separated representation able to capture rich 3D behaviors while keeping the computational complexity of 2D simulations. However, that procedure it was revealed to be too intrusive for being introduced into existing commercial softwares. Moreover, experience indicated that such enriched descriptions are only compulsory locally, in some regions or structure components. In the present paper we propose an enrichment procedure able to address 3D local behaviors, preserving the direct minimally-invasive coupling with existing plate and shell discretizations. The proposed strategy will be extended to inelastic behaviors and structural dynamics