Modeling of the 3D rocking problem
Highlights
► Two models, the CSM and WM, are developed for the 3D Rocking problem. ► The CSM uses concentrated vertical springs and dampers, and the WM distributed. ► Euler angles are used to describe rotations and Lagrange's equations are used to derive the equations of motion. ► Sliding, uplift and geometric non-linearities are included in the behavior of the models. ► Examples illustrate the ability of the models to capture 3D rocking.
Introduction
The rocking motion of a body is a highly non-linear phenomenon, that has attracted a significant number of researchers, since it was first formulated by Housner [1] in 1963. The majority of the studies on rocking have only considered the problem as two dimensional, i.e. the planar rocking of a body, and the associated complex properties. However, rocking is in general a three dimensional phenomenon and thus planar rocking is only a limited aspect of it.
In fact, even in the simplest case of a symmetric body, one cannot expect the rocking motion to be planar, unless the ground excitation is in parallel with one of the axes of symmetry of the body. The latter, however, is rarely true and thus one should expect the motion of the body to include 3D displacements and rotations. Despite this fact, the literature for the full three dimensional formulation of the problem is rather limited in comparison with the work done on the planar aspect of the phenomenon.
Housner proposed the first model, the inverted pendulum model [1], to study the planar rocking motion of a body. Expansions of this model in assemblies of two stacked rigid bodies were presented in [2], [3]. While the inverted pendulum model neglected sliding, uplift and the deformability of the support medium, a series of papers gradually incorporated these effects into the 2D rocking models. Sliding and uplift, for a rigid support medium, were included with some additional assumptions, in the work of [4], [5]. Psycharis and Jennings [6] proposed a model, considering the flexibility of the ground, with either concentrated springs and dampers at the corners, or a Winkler type of foundation, which however was geometrically linear and ignored sliding and uplift. Koh et al.'s model [7] also included a Winkler type of foundation model with geometrically non-linear equations, however still the effects of friction and uplift were not incorporated. The latter model was used to investigate the response to random input in [8]. Finally, fully geometrically non-linear CSM and WM models that included the effects of uplift and sliding were proposed in [9].
On the other hand, the literature on the modeling of the 3D rocking problem is rather limited. Konstantinids and Makris in [10] proposed a 3D rocking model for rigid symmetric prisms on a rigid foundation, however the analysis was stopped right before an impact occurred, while sliding and uplift was also neglected. 3D rocking models for rigid cylindrical objects on a rigid base were presented in [11], [12], with the assumption of no loss of energy during impacts, while the motion of a cylinder on a Winkler foundation was examined under the assumptions of no slipping or twisting in [13]. Apart from these models, it has been suggested in [14] that the 3D rocking problem can be examined with 2D approximations, an approach which as the authors state would require the existence of a constant rocking plane, i.e. the plane perpendicular to the vector of rotation. In addition, researchers [15], [16], [17] have used models based on the DEM theory as developed by Cundall [18], or general FEM models such as those proposed in [19], [20]. These models, however, were developed for 3D contact problems in general and are thus not optimized for the rocking problem.
In the current paper, two 3D rocking models are developed by expanding the benefits of the CSM and WM, proposed in [9], into the 3D rocking case. The suggested models examine rigid body prisms on a deformable base and take geometric non-linearities, sliding and uplift fully into account. The models are later used in suitable examples, where the complex 3D character of rocking is highlighted.
Section snippets
Body properties
The body in Fig. 1 is a rectangular rigid prism with centroid C and mass m. The Cartesian system S is defined by the axes X, Y, Z parallel to the unit vectors , , respectively and intersecting at O, the point of reference. Note that vectors will be denoted in bold from here on. Point O is a point fixed on the ground surface, which is chosen to match the YZ plane. Let us also refer to the Cartesian coordinates of any point i, , as . In the particular case of C the
Defining the generalized coordinates
In order to define the current position and orientation of the body, six degrees of freedom are required. The first three are simply the coordinates of the center of mass C, with respect to the Cartesian base S, i.e. xC, yC and zC. The remaining three DOF are needed to describe the rotation of the body. To this end, the sequential Euler angles , , are used as shown in Fig. 2. The initial Cartesian basis is rotated with respect to the vector by an angle , resulting in the basis
Models
In this section the two models developed to simulate the vertical reactions of the supporting medium are described and they are (i) the Concentrated Springs Model and (ii) the Winkler Model. Both models simulate the ground as vertical springs and dampers. These models produce only vertical forces creating the terms and in the system of equations (28). Note that the horizontal frictional forces at the interface between the ground and the block will be described later in Section 5.
Frictional forces
A Coulomb Model is used to simulate friction. The model considers the frictional forces to act on the bottom corners of the box. A frictional springs formulation is used, that produces only horizontal frictional forces creating the terms , and hence in the system of equations (28).
Example 1: earthquake response simulation
A rectangular concrete rigid prism of dimensions h=0.35 m, b=0.09 m, w=0.1 m is placed on a support medium that can be described by the WM with the values: , . These values would correspond to the experimentally estimated values for an elastic rubber mat foundation of modulus of elasticity 7 MPa and thickness 6 cm. The coefficient of friction has the value . The following equations are used to determine the parameters of the CSM:
These
Conclusions
Two models for examining the 3D rocking motion of a rigid rectangular prism on a flexible support medium, the 3D and the 3D , have been introduced in this study. The models describe the motion of the body by tracking the displacements of the center of mass and the rotation of the body, as expressed by the three sequential Euler rotations . The geometric non-linearities, sliding and uplift are taken fully into account and no terms are omitted in the dynamic equations of
References (31)
- et al.
Dynamic analysis of stacked rigid blocks
Soil Dynamics and Earthquake Engineering
(2001) - et al.
Analysis of block random rocking
Soil Dynamics and Earthquake Engineering
(1986) The behavior of inverted pendulum structures during earthquakes
Bulletin of the Seismological Society of America
(1963)Dynamic behaviour of rocking two-block assemblies
Earthquake Engineering and Structural Dynamics
(1990)Motion of rigid bodies and criteria for overturning by earthquake excitations
Journal of Earthquake Engineering and Structural Dynamics
(1982)- et al.
Base excitation of rigid bodies. 1: formulation
Journal of Engineering Mechanics
(1991) - et al.
Rocking of slender rigid bodies allowed to uplift
Earthquake Engineering and Structural Dynamics
(1983) - et al.
Harmonic rocking of rigid block on flexible foundation
Journal of Engineering Mechanics
(1986) - et al.
Robust modeling of the rocking problem
Journal of Engineering Mechanics
(2012) - et al.
The dynamics of a rocking block in three dimensions