A finite difference method for the static limit analysis of masonry domes under seismic loads

Abstract

The static limit analysis of axially symmetric masonry domes subject to pseudo-static seismic forces is addressed. The stress state in the dome is represented by the shell stress resultants (normal-force tensor, bending-moment tensor, and shear-force vector) on the dome mid-surface. The classical differential equilibrium equations of shells are resorted to for imposing the equilibrium of the dome. Heyman’s assumptions of infinite compressive and vanishing tensile strength, alongside with cohesive-frictional shear response, are adopted for imposing the admissibility of the stress state. A finite difference method is proposed for the numerical discretization of the problem, based on the use of two staggered rectangular grids in the parameter space generating the dome mid-surface. The resulting discrete static limit analysis problem results to be a second-order cone programming problem, to be effectively solved by available convex optimization softwares. In addition to a convergence analysis, numerical simulations are presented, dealing with the parametric analysis of the collapse capacity under seismic forces of spherical and ogival domes with parameterized geometry. In particular, the influence that the shear response of masonry material and the distribution of horizontal forces along the height of the dome have on the collapse capacity is explored. The obtained results, that are new in the literature, show the computational merit of the proposed method, and quantitatively shed light on the seismic resistance of masonry domes.

Energy-based modelling of in-plane fragility curves for the 2D ultimate capacity of Italian masonry buildings

“…However, from a computational perspective, when a non-associative behaviour is considered, searching for a solution to the boundary value problem for which normal forces vanish in the presence of normal detachments requires more sophisticated computational and modelling strategies. Nonlinear programming methods (Nodargi and Bisegna, 2021b;Kao et al, 2022;Nodargi and Bisegna, 2022) have to be used for such a purpose, or heuristic, linear or convex algorithm procedures, to lower the computational burden, although the convergence is not mathematically proven.…”

Pure compressive solutions for masonry domes under gravity loads

“…Masonry domes, widely used in monumental buildings, are great works of architecture and engineering with a long history dating back to the dawn of civilization. Their static behavior has been widely studied through different methods which can be broadly classified into equilibrium methods [1][2][3] and deformation methods [4][5][6]. The present work aims at comparing equilibrium solutions for masonry domes under gravity loads obtained by applying a membrane approach and a Finite Element (FE) analysis.…”

A solution for the frictional resistance in macro-block limit analysis of non-periodic masonry

“…Therefore, despite the extraordinary computational power available thanks to state of the art CPU processors and advanced software, often structural engineers adopt analytical approaches based on limit analysis (LA) theorems. These have the great advantage of requiring only a few material properties but, inevitably, rely on a very simplified material model [39][40][41][42][43]. In literature, LA has been formulated at both macro and micro scales.…”