Identification of kinematic excitation function by structural analysis of the system’s dynamic

Krzysztof J. Kaliński1 , Michał Mazur2

1, 2Gdańsk University of Technology, Gdańsk, Poland

1Corresponding author

Vibroengineering PROCEDIA, Vol. 3, 2014, p. 267-271.
Accepted 15 September 2014; published 10 October 2014

Copyright © 2014 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


The paper presents a method of the kinematic excitation courses identification in excitation points, based on the car road test acceleration at different measurement points. For the purpose of the laboratory fatigue life investigation of contemporary complex structures (e.g. cars bodies) and components of these structures (i.e. cars roofs), only a few first vibration modes are usually taken into account. During real life tests (i.e. road tests), accelerations at the selected points and along the defined directions are recorded. Subsequently, appropriate information about the measured acceleration allows us to identify such kinematic excitation function on the laboratory stand, whose result is the same as during the real road test. Thanks to the proposed approach, further time-consuming fatigue experiments can be performed on the laboratory stand.

Keywords: structural dynamics, finite element method, fatigue life, modal analysis.

1. Introduction

Fatigue life investigation plays principal role in scope of operation of various transportation means. The reason lies in necessity of considering the fatigue phenomenon caused by a number of load cycles in a range of 109 (so called, the gigacycle fatigue strength). Also this is well known that the experimental research is only one source of trustworthy data for the strength calculations. In fact, this is inconvenient to perform fatigue life tests of the complex structures like car roofs in the terrain. Moreover, the latter is impossible all the more so as in a range of the gigacycles. Such tests are usually performed on special laboratory stands.

There is a need to find out how the examined structure should be kinematically excited on the laboratory stand in order to achieve similar levels of accelerations in the selected measuring points as if the tests were performed on the real road. Such task of open-loop control requires solving the problem of dynamics of the investigated system. It requires integrating dynamics equation of the system of large dimensions with small step of integration; thus time of the computation may be large. Fortunately only a few first vibration modes of the investigated system are usually considered, so that the appropriate modal model can be introduced.

2. Methodology

General concept of the investigation depends on measurement of car body accelerations during operational tests in the terrain conditions within limited period of time. They allow to determine kinematic excitations of the vehicle wheels. Subsequently, the determined functions of kinematic excitation will allow to perform the car body fatigue test on laboratory stand, in a range of gigacycles.

At a current step of investigation, are inaccessible as well car for operational research as suitable laboratory stand for fatigue investigations. However, the authors dispose parametres of the car body modal model, which satisfies ability of the fatigue life computer simulation. Such simulation research is based on one technique of the mechatronic design, which is called virtual prototyping [1]. Following steps are required.

– Desired function of kinematic excitation is obtained from fragmentary road tests.

– Solution of dynamic equation of the car body calculation model.

– Determination of the desired function of vibration acceleration as response to the desired kinematic excitation.

– Identification (reconstruction) of simulated function of the excitation.

– Solution of dynamic equation of the car body calculation model.

– Determination of simulated function of vibration acceleration as response to the simulated function of excitation.

– Comparison of the simulated and the desired functions of vibration acceleration.

Matrix equation of the discrete system dynamics:

M q ¨ + L q ˙ + K q = f ,

may be written in such a way that are separated: subsystem of the unknown motion and known load, subsystem of the known motion and unknown load, and subsystems of the mutual interactions [2]:

M i i M i g M i g T M g g q ¨ i q ¨ g + L i i L i g L i g T L g g q ˙ i q ˙ g + K i i K i g K i g T K g g q i q g = f i f g .

Next we express Eq. (2) as a set of two matrix equations:

M i i q ¨ i + M i g q ¨ g + L i i q ˙ i + L i g q ˙ g + K i i q i + K i g q g = f i ,               M i g T q ¨ g + M g g q ¨ g + L i g T q ˙ i + L g g q ˙ g + K i g T q i + K g g q g = f g ,

or, after rearrangement:

M i i q ¨ i + L i i q ˙ i + K i i q i = f i - M i g q ¨ g - L i g q ˙ g - K i g q g ,               M i g T q ¨ g + M g g q ¨ g + L i g T q ˙ i + L g g q ˙ g + K i g T q i + K g g q g = f g .

First equation of the above subsystem is the matrix equation of the dynamics of the system with unknown motion at the reduced excitation (i.e. both load and kinematic). Second equation is used to compute dynamic reactions. Let us investigate only the first equation of the system of Eq. (4), which will be transformed to modal coordinates:

Ψ i T M i i Ψ i a ¨ i + Ψ i T L i i Ψ i a ˙ i + Ψ i T K i i Ψ i a i = Ψ i T f i - M i g q ¨ g - L i g q ˙ g - K i g q g ,

by the use of transformation:

q i = Ψ i a i ,

where: Ψi=Ψi1Ψi2Ψi,im – matrix of modal displacements’ vectors of the subsystem with unknown motion, ai=ai1ai2ai,im – vector of modal coordinates of the subsystem with unknown motion, im – the number of modal coordinates.

If the matrix of modal displacements’ vectors Ψi is M-orthonormal, then Eq. (5) becomes:

a ¨ i + 2 Z i Ω i a ˙ i + Ω i 2 a i = Ψ i T f i - M i g q ¨ g - L i g q ˙ g - K i g q g ,

where: Ωi=diagωi1ωi2ωi,im – matrix of angular frequencies of the subsystem with unknown motion, Zi=diagζi1ζi2ζi,im – matrix of dimensionless damping coefficients of the subsystem with unknown motion.

Modal parameters Ψi, Ωi, Zi of the subsystem with unknown motion can be obtained with the use of the experimental modal analysis techniques [3, 4].

Let us consider measurements signals’ vector acquired by the accelerometers mounted on vibrating subsystem with the unknown motion. These signals can be expressed as:

y i = y 1 y 2 y i c T ,

where: ic – number of the signals.

Vector of measured signals is related to the vector of generalized accelerations and the vector of modal accelerations of the subsystem with the unknown motion by following output equation [5]:

y i = C i Ψ i a ¨ i = C Ψ i a ¨ i ,

where: Ci – matrix of outputs of the subsystems with unknown motion, CΨi=CiΨi – modal matrix of outputs of the subsystem with unknown motion.

By using following transformation:

C Ψ i T y i = C Ψ i T C Ψ i a ¨ i ,

we obtain:

a ¨ i = C Ψ i T C Ψ i - 1 C Ψ i T y i .

Calculation of the vector of modal accelerations a¨i is possible only if a number of its components is not greater than the number of components of the vector of measurement signals y. Otherwise matrix expression inverted in Eq. (10) will be singular.

Let us consider Eq. (7) considering a few simplifications, i.e.:

– Excitation by the generalized forces vector is neglected fi0,

– Matrices of damping and inertia of the interactions subsystem are zeroes. (Mig0, Lig0).

Thus Eq. (7), after simplifications, has a form:

a ¨ i + 2 Z i Ω i a ˙ i + Ω i 2 a i = - Ψ i T K i g q g .

And after transformations:

K i g T Ψ i a ¨ i + 2 Z i Ω i a ˙ i + Ω i 2 a i = - K i g T Ψ i Ψ i T K i g q g ,

we obtain vector of generalized coordinates of the subsystem with known motion:

q g = - K i g T Ψ i Ψ i T K i g - 1 K i q T Ψ i a ¨ i + 2 Z i Ω i a ˙ i + Ω i 2 a i .

It is, at the same time, vector of sought kinematic excitations, required as the input for the measured output responses vector y of the subsystem with unknown motion. It is only possible to compute vector of generalized coordinates qg if its dimension is not greater than dimension of the modal coordinates’ vector ai. Otherwise the matrix expression inverted in Eq. (14) will be singular.

In this simplified approach only the knowledge about the stiffness properties (components of stiffness matrix Kig) of the subsystem of mutual interaction is required for the computation of the kinematic excitation functions qgt. Elements of the stiffness matrix Kig can be usually easily obtained as they are results of the stiffness properties of the elements being kinematically excited.

Fig. 1. Model of the structure. Red x – excitations points, blue dots – measurement points

Fig. 2. Exemplary time displacements of the structure for node 56, green line – simulated plot, blue dotted line – measured plot

3. Computer simulation

Model of the structure is presented in Fig. 1. It is mainly based on beam elements, but there are also applied mass, spring and damping elements. The model was kinematically excited along vertical direction in four nodes 56, 262, 267 and 273. The responses were registered along three directions in ten nodes 1, 11, 17, 34, 72, 86, 88, 89, 288 and 296. Time modal history was synthesized with the use of the twelve first modes for the period of time 2 s and step of integration 0.001 s.

All twelve modal displacements were identified with the use of Eq. (11). Then the generalized coordinates of the subsystem with the known motion qgt were computed with the use of Eq. (14). Results for the selected node 56 are presented in Fig. 2. The obtained signals are nearly identical with the given one.

4. Conclusions

The results of computer simulation (virtual prototyping) showed very good compliance of the given signals of kinematic excitation with the computed ones. The latter also concerns conformity of the desired and simulated vibrations of the car body. It is possible to effectively compute kinematic excitations’ functions, even if the simulated and given signals are not the same. Numerical computations are effective due to the use of modal model of the car body. In case of real structures, modal model can be easily obtained with the use of experimental modal analysis techniques. Thus, the modal model parameters obtained during real experiments, can be applied for a purpose of virtual prototyping.

It should be noted that despite ability of good adjustment of the kinematic excitation function to real working conditions of the structure, virtual prototyping would not give an answer to the question what the fatigue endurance is. In order to determine fatigue endurance, suitable programme of experiments should be performed on the laboratory stand (number of cycles in a range of 109), at kinematic excitation determined on a stage of virtual prototyping.

The economic criterion is fulfilled in the end. Laboratory investigations are cheaper than testing of ready transportation means during their routine operation.


The authors kindly acknowledge firms vMACH Engineering GmbH of Markt Indersdorf, Germany, and Macross GmbH of Munich, Germany, for providing the research with some experimental data concerned modal model of the car body, as well as – real kinematic excitations of chassis vehicles.


  1. Kaliński K. J., Buchholz C. Mechatronic design of strongly nonlinear systems on a basis of three wheeled mobile platform. Mechanical Systems and Signal Processing, 2014, [CrossRef]
  2. Kruszewski J., Gawroński W., Wittbrodt E., Najbar F., Grabowski S. The Rigid Finite Element Method. Arkady, Warszawa, 1975, (in Polish). [CrossRef]
  3. Maia N. M. M., Silva J. M. M. Theoretical and Experimental Modal Analysis. Research Studies Press, Taunton, Somerset, England, 1997. [CrossRef]
  4. Heylen W., Lammens S., Sas P. Modal Analysis Theory and Testing. KU Leuven, Leuven, 2007. [CrossRef]
  5. Kaliński K. J. A Surveillance of Dynamic Processes in Mechanical Systems. Publishing House of Gdańsk University of Technology, Gdańsk, 2012, (in Polish). [CrossRef]