## Journal of the Brazilian Society of Mechanical Sciences

*Print version* ISSN 0100-7386

### J. Braz. Soc. Mech. Sci. vol.24 no.3 Rio de Janeiro July 2002

#### http://dx.doi.org/10.1590/S0100-73862002000300010

**Seismic motion simulation based on Cassino Parallel Manipulator**

**J. C. M. Carvalho ^{I}; M. Ceccarelli^{II}**

^{I}Universidade Federal de Uberlândia Faculdade de Engenharia Mecânica 38400-902 Uberlândia , MG. Brazil. E-mail: jcmendes@mecanica.ufu.br

^{II}Università di Cassino Dipartimento di Meccanica, Strutture, Ambiente e Territorio Via Di Biasio 43, 03043 Cassino (Fr), Italy. E-mail: Ceccarelli@unicas.it

**ABSTRACT**

In this paper we present a study of feasibility by using Cassino Parallel Manipulator (CaPaMan) as an earthquake simulator. We propose a suitable formulation to simulate the frequency, amplitude and acceleration magnitude of seismic motion by means of the movable platform motion by giving a suitable input motion. In this paper we have reported numerical simulations that simulate the three principal earthquake types for a seismic motion: one at the epicenter (having a vertical motion), another far from the epicenter (with the motion on a horizontal plane), and a combined general motion (with a vertical and horizontal motion).

**Keywords:** Parallel manipulator, earthquake, earthquake simulation, seismic motion

**Introduction**

Parallel structures have characteristics such as high stiffness and load capacity and can operate at high velocities and accelerations with high precision. This justifies the great number of theoretical and experimental studies with the aim to use them more than serial manipulators. It has been found applications for flight simulators, milling machines, toys, load cells, packing and assembly.

At the Laboratory of Robotics and Mechatronics in Cassino, Italy, a novel parallel structure with 3 d.o.f., named as CaPaMan – Cassino Parallel Manipulator – has been conceived and suitable formulation for kinematics and dynamics has been developed by theoretical investigations.

Earthquakes are well known from the destruction and lost lives due to the failing of buildings. In general a combined horizontal and vertical ground motion characterizes an earthquake, and it can be conveniently simulated by using a suitable parallel structure. The simulation of a seismic motion has great interest because it helps to understand the earthquake effects and helps to solve problems in predicting the behavior of buildings and civil engineering structures in general.

In this paper we present a study of feasibility for an earthquake simulator by using the CaPaMan structure. The frequency, amplitude and acceleration magnitude of seismic motion are simulated by the movable platform motion excited by a suitable input motion.

In particular, we propose a specific formulation in order to simulate the three principal types of seismic motion and by considering the position related to the epicenter. If the earthquake is considered in the epicenter the seismic motion takes place along the vertical direction. By considering the seismic motion far from the epicenter, the motion ground occurs mainly in a horizontal plane, and otherwise takes place a general motion, i.e. a vertical motion that is combined with a horizontal motion.

**CaPaMan Structure**

CaPaMan (Cassino Parallel Manipulator) structure is characterized by three 4-bar articulated mechanisms AP disposed in an equilateral triangle of the fixed base FP as shown in Fig.1. On the coupler link of each 4-bar mechanism is installed a prismatic joint SJ that is connected with the movable platform MP by a rod CB and a spherical joint BJ. The kinematic variables are the input crank angles a_{k} (k=1,2,3) of each 4-bar mechanism.

In order to describe the CaPaMan kinematic and dynamic behavior, an inertial frame OXYZ has been assumed fixed to FP and a moving frame HX_{p}Y_{p}Z_{p} has been attached to the platform. The kinematic feasibility of CaPaMan was analysed by Ceccarelli (1997) and Ceccarelli and Figliolini (1997). The position and orientation of MP has been expressed using the coordinates y_{k} and *z*_{k} of the articulation points H_{k} that can be easily expressed as a function of the AP input angles a_{k} (k=1,2,3). Since the motor units for robotic applications are usually controlled in position and velocity, the input motion can be given by a cubic function of time t between given initial a_{ki} and final a_{kƒ} angles at initial *t** _{ki}* and final

*t*

_{kƒ}times, respectively. An analytical model for dynamics of CaPaMan by using the Newton-Euler approach has been formulated to compute the input torques which are necessary for a given trajectory of movable platform (Carvalho and Ceccarelli, 1999, 2000). Experimental tests has been confirmed the kinematic feasibility and actuation efficiency of CaPaMan (Ceccarelli et al., 1999a).

**Earthquake Characteristics**

The intensity, shape and duration, that are functions of the terrain properties, can be used to classify an earthquake. Newmark-Rosenblueth (1971) classifies earthquakes in three types:

Type 1 - Type 1 earthquake is characterized by a single shock. Generally occurs no long of the epicenter and in a compact soil as those that has been occurred in Agadir in 1960, Skopje in 1963 and Port Hueneme in 1957 (Gavarini, 1984; Sarà, 1985). An example is shown in Fig.2.

– Type 2 – It is an earthquake characterized by a long duration, it is very irregular, and it has a large frequency spectrum. Generally occurs at a medium distance of the epicenter and in a compact soil as the El Centro earthquake (California/USA) that has been occurred in 1940 as shown in Fig.3.

– Type 3 – This earthquake type is characterized by a predominant frequency as the Mexico City earthquake occurred in 1964, shown in Fig.4. The terrain acts as a filter for the reflection waves of the seismic motion.

There is a fourth earthquake type that is characterized by an important deformation on the ground (with a magnitude of a meter) which we do not consider in this paper.

A typical acceleration seismogram is characterized by an initial phase that corresponds to the beginning of the seismic motion; an intermediate phase during which occurs the maximum acceleration peaks and displacements, and a final phase that represents the end of the earthquake as shown in Figs. 2, 3 and 4. Characteristics of the terrain and distance of the epicenter can give to seismic motion a preferential direction: along a vertical line, on a horizontal plane, or both horizontal and vertical directions.

The feasibility of use CaPaMan as an earthquake simulator is shown through results of numerical simulations of three earthquake types, and by considering a vertical, or a horizontal or a vertical and horizontal motion.

**Earthquake Simulation by CaPaMan**

An earthquake simulation by using CaPaMan can be performed by properly giving an input motion to links b_{k} (k=1,2,3). Crank angles a_{k} (k=1,2,3) give the amplitude of seismic motion, and acceleration peaks for the displacement of MP. The random characteristics of an earthquake give the possibility to define limits for maximum and minimum intervals for input angles a_{k} and time cycle.

In order to obtain a general earthquake shape, the initial and final phases of an accelerogram has been divided in two parts, as shown in Figs. 5 and 6.

Figure 5 represents a typical shape of an earthquake defined by limits of input angles. The beginning of the seismic motion is limited by cycle i_{2}; and the final phase from cycle i_{3} to cycle i_{f} . Between cycles i_{2} and i_{3} occurs the acceleration peaks of earthquake.

In order to simulate a single shock earthquake, a maximum and minimum shock input angle a _{shock _ max/min }has been defined between cycles i_{2} and i_{3} as shown in Fig.6.

Acceleration peaks for earthquakes of types 1 and 2 have been obtained from cycle time as shown in Fig.7. Small cycle times associated with large input angles at interval i_{2}i_{3} correspond to large accelerations peaks. Type 3 earthquake has a constant frequency.

Referring to Figs. 5 and 6 each input angle a* _{k}* (k=1,2,3) has been given according to the following expressions:

For earthquakes of types 2 and 3 the expression of each input angle in the interval i_{2} and i_{3} takes the form:

and for single shock earthquake it holds as

Referring to Fig.7, input time for each cycle can be obtained as

Between cycles i2 and i3 the input time for single shock earthquake is given as

for earthquake of type 2 is given as

and, for earthquake of type 3 is given by the earthquake frequency *F _{req}* as

*R* in Eqs. (1) to (12) is a random function in the interval (0,1) which has been successfully used to simulate seismic motions.

Kinematics characteristics of CaPaMan have been analyzed in several papers (Ceccarelli, 1997; Ceccarelli & Figliolini, 1997; Carvalho & Ceccarelli, 1999; Ceccarelli et al., 1999a and 1999b). And they can be easily used to simulate a seismic motion along a vertical direction, a horizontal plane and a general direction, i.e. vertical and horizontal motion by simulating the distance of the epicenter.

**Numerical Simulations**

In this section, we reported results for three numerical simulations: a single shock earthquake in a vertical direction Fig.8; an earthquake of type 2 in a horizontal direction Fig.9, and an earthquake of type 3 in a general direction Fig.10.

Reported simulations refer to the specific dimensions of the built prototype in Cassino, Fig.1: c_{k}=a_{k}=200mm; b_{k}=d_{k}=80mm; h_{k}=116mm; r_{f}=r_{p}=109,5mm and m=2,912kg, k=(1,2,3).

Results of the numerical simulation of a single shock earthquake in a vertical direction are shown in Fig.8 and they show accelerations of the same magnitude as those presented in Fig.2. The input parameters defining the seismic wave have been given as: a* _{0}*=90deg; a

*a*

_{i1max}=*a*

_{i2max}=*16,5deg; a*

_{i4max}=*a*

_{i1min}=*a*

_{i2min}=*-16,5deg; a*

_{i4min}=*20deg; a*

_{shock–max}=*-20deg; t*

_{shock–min}=_{max}=0,4s; t

_{min}=0,3s; t

_{ifmax}=t

_{ifmin}=0,25s; number_of_cycles=50. The vertical linear motion of MP is obtained when the input crank angles a

_{k}(k=1,2,3) are equal as shown in Fig.8a. Figure 8b presents time history of vertical acceleration. Fig.8c gives the values of required input torques and Fig.8d gives the consequent displacement of point H of MP.

The earthquake of type 2 has been simulated for a horizontal motion and the results are shown in Fig.9. The parameters for the seismic wave are: a * _{0}*=90deg; a

*a*

_{i1max}=*a*

_{i2max}=*16,5deg; a*

_{i4max}=*a*

_{i1min}=*a*

_{i2min}=*-16,5deg; t*

_{i4min}=_{max}=0,4s; t

_{min}=0,3s; t

_{ifmax}=t

_{ifmin}=0,2s; number_of_cycles=70. Figure 9 shows results for a horizontal motion of MP along the axis x when the input angle a

_{1}is assumed as constant, a

_{2}is given by Eqs. (1) to (5) and a

_{3}= p - a

_{2}is ensured.

The results of the numerical simulation are shown in Fig.9 and a small acceleration component in the vertical direction has been computed in agreement with what occurs in a real horizontal earthquake, Figs.9g and 9h.

Figure 10 shows results of a general seismic motion for the type 3 earthquake. The values of the parameters for the seismic wave are assumed as: F_{req}=2,4cycle/sec; a* _{0}*=70deg; a

*a*

_{i1max}=*a*

_{i2max}=*35deg; a*

_{i4max}=*a*

_{i1min}=*a*

_{i2min}=*-35deg; number_of_cycles=310.*

_{i4min}=

**Conclusions**

In this paper we have presented a study of feasibility by using CaPaMan (Cassino Parallel Manipulator) as earthquake simulator. A suitable formulation has been presented as based on seismogram of the three principal types of earthquakes that can de identified as a seismic motion along a vertical direction, a horizontal plane and a combined general motion direction.

Several numerical tests have been carried out and results prove the feasibility of CaPaMan to simulate any type of seismic motion "in" or "out" of the epicenter since a suitable actuation has been computed in the form of input torques for the three motors of CaPaMan.

Practical experiments, reported in (Ceccarelli et al., 1999b), have been shown the feasible seismic simulated motion of the mobile platform by using a simple programming of the earthquake operation.

**Acknowledgments**

The first author is thankful to the Brazilian Ministry of Science and Technology – CNPq for the grant, which permitted him to spend a sabbatical year in 1998 at the Laboratory of Robotics and Mechatronics in Cassino, Italy.

**References**

Braga, F. (Ed.), 1997, "Elementi di Ingegneria Antisismica", International Centre for Mechanical Sciences, Udine, Italy. [ Links ]

Carvalho, J.C.M. & Ceccarelli, M., 1999, "A Dynamic Analysis for Cassino Parallel Manipulator", Proceedings of the 10^{th} World Congress on the Theory of Machines and Mechanisms, Oulu, Finland, Vol.3, pp. 1202-1207. [ Links ]

Carvalho, J.C.M. & Ceccarelli, M., 2000, "A Closed Form Formulation for the Inverse Dynamics of Cassino Parallel Manipulator", Multibody System Dynamics, Accepted for publication. [ Links ]

Ceccarelli, M., 1997, "A New 3 dof Spatial Parallel Mechanism", Mechanism and Machine Theory, Vol. 32, No. 8, pp. 895-902. [ Links ]

Ceccarelli, M. & Figliolini, G., 1997, "Mechanical Characteristics of CaPaMan (Cassino Parallel Manipulator)", Proceedings of the 3^{rd} Asian Conference on Robotics and its Application, Tokyo, Japan, pp. 301-308. [ Links ]

Ceccarelli, M., Pugliese, F., Carvalho, J.C.M., 1999a, "An Experimental System for Measuring CaPaMan Characteristics", Proceedings of the 8^{th} International Workshop on Robotics in Alpe-Adria-Danube Region, Munich, Germany, pp. 31-36. [ Links ]

Ceccarelli, M., Carvalho, J.C.M., Pugliese, F., Lanni, C., 1999b, "CaPaMan (Cassino Parallel Manipulator) as Sensored Earthquake Simulator", Proceedings of the 1999 IEEE/RSJ International Conference on Intelligent Robots and Systems, Kyongju, Korea, Vol.3, pp. 1501-1506. [ Links ]

Gavarini, C., 1984, "Dinamica Delle Strutture", Ed. Scientifiche Associate, Roma, Italy. [ Links ]

Newmark, N.M., Rosenbluet, E.R., 1971, "Fundamentals of Earthquake Engineering", Prentice-Hall, New York, USA. [ Links ]

Sarà, G. (Ed.), 1985, "Ingegneria Antisismica", vol 1, Liquore Editore, Napoli, Italy. [ Links ]