Rocking Response of the Rigid Block Under Rectangular Pulse Excitation: A Comparison between ODE and Optimization-Based Solvers


Rigid block
ODE solver
Mathematical programming
Differential equation of motion

How to Cite

Gagliardo, R. ., Portioli, F. ., Cascini , L. ., & Landolfo, R. . (2022). Rocking Response of the Rigid Block Under Rectangular Pulse Excitation: A Comparison between ODE and Optimization-Based Solvers. Journal of Advances in Applied & Computational Mathematics, 8, 109–116.


In this paper, the response of the rigid block under rectangular pulse excitation is investigated using two different modeling approaches and solvers. The first approach relies on the numerical integration of the differential equation of motion. The second approach is based on the formulation of the dynamic problem in terms of a special class of mathematical programming problem that is the linear complementarity problem. A validation study is carried out comparing the solutions given by the proposed formulation with the ones given by the numerical integration of the differential equation of motion obtained from ODE solvers available in MATLAB®. Potentialities and limitations of the mathematical programming formulation are discussed in terms of energy dissipation and restitution coefficient at impacts and in terms of solution times.


Lagomarsino S. Seismic assessment of rocking masonry structures. Bull Earthq Eng 2015;13:97-128.

Makris N, Roussos Y. Rocking Response and Overturning of Equipment Under Horizontal Pulse-Type Motions. Pacific Earthquakes, Eng Res Cent 1998.

Makris N, Zhang J. Rocking Response and Overturning of Anchored Equipment under Seismic Excitations. Pacific Earthq Eng Res Cent 1999:1-82.

Spanos PD, Koh A-S. Rocking of rigid blocks due to harmonic shaking. J Eng Mech 1984;110:1627-42.

Psycharis IN, Jennings PC. Rocking of slender rigid bodies allowed to uplift. Earthq Eng & Struct Dyn 1983;11:57-76.

Ishiyama Y. Motions of rigid bodies and criteria for overturning by earthquake excitations. Earthq Eng & Struct Dyn 1982;10:635-50.

Aslam M, Fodden WG, Scalise DT. Earthquake rocking response of rigid bodies. J Struct Div ASCE 1980;106:377-92.

Malomo D, Mehrotra A, DeJong MJ. Distinct element modeling of the dynamic response of a rocking podium tested on a shake table. Earthq Eng Struct Dyn 2021;50:1469-75.

Grillanda N, Chiozzi A, Milani G, Tralli A. Tilting plane tests for the ultimate shear capacity evaluation of perforated dry joint masonry panels. Part II: Numerical analyses. Eng Struct 2021;228.

Giresini L, Solarino F, Paganelli O, Oliveira DV, Froli M. ONE-SIDED rocking analysis of corner mechanisms in masonry structures: Influence of geometry, energy dissipation, boundary conditions. Soil Dyn Earthq Eng 2019;123:357-70.

Giouvanidis AI, Dong Y. Seismic loss and resilience assessment of single-column rocking bridges. Bull Earthq Eng 2020;18:4481-513.

Giouvanidis AI, Dimitrakopoulos EG. Rocking amplification and strong-motion duration. Earthq Eng Struct Dyn 2018;47:2094-116.

Galassi S, Ruggieri N, Tempesta G. A Novel Numerical Tool for Seismic Vulnerability Analysis of Ruins in Archaeological Sites. Int J Archit Herit 2020;14:1-22.

Funari MF, Mehrotra A, Lourenço PB. A tool for the rapid seismic assessment of historic masonry structures based on limit analysis optimisation and rocking dynamics. Appl Sci 2021;11:1-22.

De-Felice G, Malena M. Failure pattern prediction in masonry. J Mech Mater Struct 2019;14:663-82.

Casapulla C, Giresini L, Lourenço PB. Rocking and kinematic approaches for rigid block analysis of masonry walls: State of the art and recent developments. Buildings 2017;7.

Acary V, Brogliato B. Numerical methods for non-smooth dynamical systems. Berlin: Springer; 2008.

Jean M. The non-smooth contact dynamics method. Comput Methods Appl Mech Eng 1999;177:235-57.

Rafiee A, Vinches M, Bohatier C. Modelling and analysis of the Nîmes arena and the Arles aqueduct subjected to a seismic loading, using the Non-Smooth Contact Dynamics method. Eng Struct 2008;30:3457-67.

Lancioni G, Lenci S, Piattoni Q, Quagliarini E. Dynamics and failure mechanisms of ancient masonry churches subjected to seismic actions by using the NSCD method: The case of the medieval church of S. Maria in Portuno. Eng Struct 2013;56:1527-46.

Andersen ED, Roos C, Terlaky T. On implementing a primal-dual interior-point method for conic quadratic optimization. Math Program Ser B 2003;95:249-77.

Tasora A, Anitescu M. A matrix-free cone complementarity approach for solving large-scale, nonsmooth, rigid body dynamics. Comput Methods Appl Mech Eng 2011;200:439-53.

Stewart DE, Trinkle JC. An implicit time-stepping scheme for rigid body dynamics with inelastic collisions and coulomb friction. Int J Numer Methods Eng 1996;39:2673-91.<2673::AID-NME972>3.0.CO;2-I.

Anitescu M, Potra FA. A time-stepping method for stiff multibody dynamics with contact and friction. Int J Numer Methods Eng 2002;55:753-84.

Livesley RK. Limit analysis of structures formed from rigid blocks. Int J Numer Methods Eng 1978;12:1853-71.

Baggio C, Trovalusci P. Limit analysis for no-tension and frictional three-dimensional discrete systems. Mech Struct Mach 1998;26:287-304.

Baggio C, Trovalusci P. Collapse behavior of three-dimensional brick-block systems using non-linear programming. Struct Eng Mech 2000;10:181-95.

Cascini L, Gagliardo R, Portioli F. LiABlock_3D: A Software Tool for Collapse Mechanism Analysis of Historic Masonry Structures. Int J Archit Herit 2020;14:75-94.

Orduña A, Lourenço PB. Three-dimensional limit analysis of rigid blocks assemblages. Part I: Torsion failure on frictional interfaces and limit analysis formulation. Int J Solids Struct 2005;42:5140-60.

Orduña A, Lourenço PB. Three-dimensional limit analysis of rigid blocks assemblages. Part II: Load-path following solution procedure and validation. Int J Solids Struct 2005;42:5161-80.

Portioli F, Casapulla C, Cascini L. An efficient solution procedure for crushing failure in 3D limit analysis of masonry block structures with non-associative frictional joints. Int J Solids Struct 2015;69-70:252-66.

Portioli FPA. Rigid block modeling of historic masonry structures using mathematical programming: a unified formulation for non-linear time history, static pushover, and limit equilibrium analysis. Bull Earthq Eng 2020;18:211-39.

Fishwick RJ. Limit analysis of rigid block structures—University of Portsmouth, 1996.

Ferris MC, Tin-Loi F. Limit analysis of frictional block assemblies as a mathematical program with complementarity constraints. Int J Mech Sci 2001;43:209-24.

Krabbenhoft K, Huang J, Da Silva MV, Lyamin AV. Granular contact dynamics with particle elasticity. Granul Matter 2012;14:607-19.

Lim K-W, Krabbenhoft K, Andrade JE. A contact dynamics approach to the Granular Element Method. Comput Methods Appl Mech Eng 2014;268:557-73.

Portioli F, Cascini L. Contact Dynamics of Masonry Block Structures Using Mathematical Programming. J Earthq Eng 2018;22:94-125.

Zhang J, Makris N. Rocking response of free-standing blocks under cycloidal pulses. J Eng Mech 2001;127:473-83.

Creative Commons License

This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

Copyright (c) 2021 R. Gagliardo, F.P.A. Portioli, L. Cascini , R. Landolfo