System Identification Algorithm for Computing the Modal Parameters of Linear Mechanical Systems

The goal of this investigation is to construct a computational procedure for identifying the modal parameters of linear mechanical systems. The methodology employed in the paper is based on the Eigensystem Realization Algorithm implemented in conjunction with the Observer/Kalman Filter Identification method (ERA/OKID). This method represents an effective and efficient system identification numerical procedure based on the time domain. The algorithm developed in this work is tested by means of numerical experiments on a full-car vehicle model. To this end, the modal parameters necessary for the design of active and semi-active suspension systems are obtained for the vehicle system considered as an illustrative example. In order to analyze the performance of the methodology developed in this investigation, the system identification numerical procedure was tested considering two case studies, namely a full state measurement and an incomplete state measurement. As expected, the numerical results found for the identified dynamical model showed a good agreement with the modal parameters of the mechanical system model. Furthermore, numerical results demonstrated that the proposed method has good performance considering a scenario in which the signal-to-noise ratio of the input and output measurements is relatively high. The method developed in this paper can be effectively used for solving important engineering problems such as the design of control systems for road vehicles.


Introduction
In abstract terms, a mathematical model of a dynamical system is a specific pattern that establishes relationships between the physical input and output quantities that can also be observed experimentally.A mathematical model of a dynamical system, therefore, allows for predicting the system dynamic behavior in response to a given set of external excitations [1].The complexity of a dynamical model reflects the model dimension, its general flexibility and the capacity of representing the phenomena of interest in terms of mathematical laws [2][3][4][5].The complexity of a given dynamical model is strongly correlated with the mathematical structure of the model and to the information content present in the observed dataset used in the model development.In this context, the estimation process is an iterative process in which the mathematical structure of a dynamical model is devised and, once a suitable mathematical model is constructed for capturing the physics of the system of interest, one can identify the unknown parameters in the model by using an experimental dataset [6][7][8][9].This process is guided by the intuition of the analyst and is based on the useful information embedded in the dataset available from the experimental campaign [10][11][12][13][14][15][16].In the validation process, on the other hand, the analyst verifies that the selected dynamical model and the identified model parameters are able to reproduce the input-output datasets that are different from those used for constructing the system model by means of the system identification methods.The main advantage of having an analytical model of a mechanical system is represented by the possibility of performing numerical experiments by using dynamic simulations, which are inexpensive, safe and can be repeated for reproducing a large number of different realistic scenarios [17][18][19].However, the accuracy and reliability of the numerical results computed through dynamic simulations completely depend on the quality of the mathematical model of the physical system identified by using the observed dataset [20][21][22].
This investigation is grounded in the mathematical tools of the scientific discipline called system identification.Originally developed in the field of control theory, system identification has become an independent field of scientific research that is focused on the determination of a mathematical model for a physical system by using experimental data [23][24][25].In order to achieve this goal, the analyst should combine the information arising from experimental data with the background knowledge on the functioning of a particular dynamical system.The mathematical techniques used in the field of the system identification serve a threefold objective, namely: they can be employed for identifying a set of unknown parameters of a nonlinear mathematical model by using an experimental dataset; they can be used for the determination of the model order corresponding to a given system; and they allow for obtaining finite-dimensional linear models of a general mechanical system starting from an interrelated sequence of input-output data [26,27].The methods of system identification can be applied to deterministic systems, as well as to stochastic processes [28,29].Furthermore, several effective system identification methods have been developed in the literature such as, for example, adaptive algorithms, least square methods, gradient-based procedures and multi-innovation identification methodologies [30][31][32][33][34][35].In engineering applications, the methods of system identification are principally used for performing the experimental modal analysis, which represents the principal topic of interest of this investigation.To this end, two general approaches can be employed, namely the time-domain analysis and the frequency-domain analysis [36].In the time-domain analysis, the input-output dataset obtained in an experimental campaign is directly provided to the system identification numerical procedures based on the time domain.On the other hand, in the frequency domain analysis, the experimental data are filtered and used by means of a Fourier transform, and subsequently, they are used in a system identification numerical procedure based on the frequency domain.In both cases, the analyst is interested in the calculation of the number of normal modes and natural frequencies necessary for constructing a reliable dynamic model of a given mechanical system without directly deriving the system equations of motion [37][38][39][40].
In mechanical engineering applications, various studies focused on the active and passive control problem of machines and structures based on the mathematical techniques developed in the field of system identification and optimal control can be found [41][42][43][44].The optimal control problem is particularly challenging in the case of rigid and flexible multibody systems, which represent mechanical systems formed by continuum bodies, kinematic joints, force fields and control actuators [45][46][47][48][49][50].Another important application of the system identification scientific discipline is the refinement of a finite element model through dynamic testing for the design of vibration control systems based on open-loop and closed-loop control schemes [51][52][53][54].The methodologies of time domain analysis employed in the general area of applied system identification are of interest for this investigation because the state-space models obtained through these numerical procedures can be readily used for developing effective control actions employing well-established and robust algorithms, such as, for example, the pole placement method, the linear-quadratic regulator technique and the H-infinity method, as well as more advanced nonlinear control approaches [55][56][57][58][59][60][61].In the case of the parametric identification problem, one needs to construct analytical expressions in terms of the desired parameters and correlate them with experimental data in order to identify the parameters of interest [62,63].For this purpose, simple least-square methods, as well as more advanced gradient-based optimization techniques relying on the adjoint method can be used for minimizing the error between the time responses observed in experimental measurements and the time evolutions predicted by the mathematical model based on the unknown parameters [64][65][66][67][68][69].On the other hand, in the case of the model identification problem in the time domain, one needs to find a linear dynamical model directly that embodies the best fit for an input and output dataset measured for a given mechanical system.The model identification problem is based on analytical techniques and numerical linear algebra methods that are of interest for this work and, therefore, are discussed in this paper.
This manuscript is subdivided as follows.In Section 2, the mathematical background related to the system identification numerical procedure developed in this investigation is provided.In Section 3, a full-car mechanical model of a standard road vehicle is analyzed and used as a demonstrative example for the system identification numerical procedure considered in this work.In Section 4, the summary of the paper, the conclusions of these investigations and a discussion about the future directions of research are reported.

System Identification Algorithm
This section describes the background material necessary for the analytical development and the computer implementation of the system identification numerical procedure employed in this investigation for obtaining the modal parameters of a mechanical system characterized by a linear mathematical structure.

Representation of the Dynamical Model in the Space of States
In this subsection, the key aspects of the mathematical description of a mechanical system considering a representation in the space of states are illustrated.For this purpose, the dynamic equations of a mechanical system having a linear mathematical structure are considered.The equations of motion form a coupled set of n 2 second-order linear ordinary differential equations, where n 2 is the number of the generalized coordinates used in the mathematical model that form an independent set of time-varying parameters.In particular, assume that M, R and K represent the mass, damping and stiffness matrices, respectively, of a mechanical system described by a linear dynamical model.The dynamic equations can be written in matrix form given by: where the variable t is the continuous time, while x(t) is the generalized displacement vector, ẋ(t) is the generalized velocity vector, ẍ(t) is the generalized acceleration vector and F(t) is an external force that is known as an explicit function of time.One can record the time response of the mechanical system of interest considering an output vector denoted with y(t), which is formed by m measurements.Thus, the output vector is a function of time that is described by a set of output equations given by: where C a , C v and C d are output influence matrices associated with the generalized accelerations, velocities and displacements, respectively.The influence matrices associated with the output vector define the relation between the configuration vectors and the output measurements available in practical applications that are contained in the measurement vector.One can denote with z(t) the state vector of the mechanical system defined as follows: The external generalized forces acting on the mechanical system are defined by the r input quantities, which are grouped in the input vector that is denoted with u(t).One can, therefore, use the state vector z(t), as well as the generalized force vector u(t) for rewriting the differential set of dynamic equations and the algebraic measurement equations as follows: where A c denotes the state matrix represented in the continuous time domain, B c denotes the state influence matrix represented in the continuous time domain, C is referred to as the measurement influence matrix and D is called the direct transmission matrix.One can explicitly compute these state-space matrices starting from the system mechanical model represented in the configuration space in the following way: where B 2 is a constant Boolean matrix that characterizes the structure of the external generalized forces as linear functions of the inputs.The Boolean matrix B 2 associated with the input vector u(t) defines the following matrix equation: The set of Equation ( 4) constitutes a state-space model of a mechanical system characterized by a linear dynamics derived considering a continuous-time representation.Moreover, one can construct the discrete-time counterpart of the continuous-time state space model of a dynamical system associated with Equation ( 4), leading to the following set of equations: where the variable k denotes the discrete time, A is the state matrix expressed in the discrete time domain and B is the state influence matrix represented in the discrete time domain.The transformation in the space of states from the continuous-time model of the mechanical system to the discrete-time model is necessary because, in practical applications, the input and output data are acquired in a discrete fashion.Consequently, the discrete-time model of the mechanical system of interest given by Equation ( 7) represents the basic set of equations necessary for the development of a class of system identification numerical procedures based on the time domain that are suitable for handling general mechanical models of linear dynamical systems.In particular, starting from the continuous-time matrices A c and B c , one can construct the corresponding discrete-time matrices A and B by performing a time sampling of the state-space model.By doing so, one can write: where ∆t denotes the constant time step used in the discretization process for creating an equally-spaced time grid.The matrices C and D that appear in the measurement equations, on the other hand, are not affected by the time discretization because the measurement equations are algebraic equations.

Definition of the Markov Parameter Set
In this subsection, the Markov parameter set associated with a mechanical system is introduced.To this end, the set of Equation ( 7) that forms the state-space representation of a linear dynamical system based on a discrete-time description is considered.For simplicity, assume that the set of initial conditions is homogeneous, namely z(0) = 0. Thus, by recursively solving the state-space Equation (7) and considering the entire time history of the input vector u(k), one obtains: In order to compute the response of the mechanical system to an impulse excitation corresponding to a given entry of the input vector of the input variables, one can consider separately the influence on the time response of an impulsive excitation applied to each entry of the input vector and exploit the principle of superposition of effects that characterize each linear dynamical system.By doing so, one can assemble the analytical results into a set of impulse-response matrices Y k having dimensions m by r known as system Markov parameters and given by: Since the structural matrices that define the discrete-time dynamical model in the state-space as reported in Equation ( 7) are incorporated in the set of Markov parameters, these parameters can be used for identifying a mathematical model of a linear dynamical system obtained from a set of input-output experimental data.Employing the set of Markov parameters, one can rewrite the measurement equations as follows: The set of Equation ( 11) represents an input-output mathematical relationship associated with the state-space model of the mechanical system based on the discrete time domain in which the contribution to the output vector y(k) at time step k is obtained from the time history of the input vector u(k) weighted by the sequence of the system Markov parameters Y k given by Equation (10).

Representation of the Observer Model in the Space of States
In this subsection, the development of a state-space model associated with a state estimator used for predicting the dynamic behavior of a mechanical system is discussed.The state-space description of the dynamical system is based on the definition of the state vector z(k), which allows for obtaining the relation between the input and output vectors u(k) and y(k) associated with the mechanical system.In practical applications, the state vector is not directly accessible from experimental measurements, and therefore, the introduction of a state observer denoted with the rectangular matrix G becomes necessary.By introducing the state observer in the state-space model of a mechanical system based on the discrete time domain, one obtains: where ẑ(k) is an estimated or observed discrete-time state vector and ŷ(k) is an estimated or observed discrete-time measurement vector.The resulting state matrix and influence matrix Ā and B are modified by the introduction of the state estimator represented by the observer matrix G as follows: The set of Equation ( 12) defined in the discrete time domain forms an observer-based state-space model of the linear mechanical system of interest.In Equation (12), the observer matrix G should be designed in order to change the eigenvalues of the state matrix A and lead to an asymptotically stable dynamical system characterized by the observer state matrix Ā.For a sufficiently large time step k, if the estimated output vector ŷ(k) tends to the actual measurement vector y(k) obtained from experimental measurements, then the estimated state ẑ(k) should approach the real state z(k), which is not directly measurable.Furthermore, the discrete-time observer state-space model given by Equation ( 12) is also useful to cope with the uncertainties in the mathematical model that describes the mechanical system, for handling the noise disturbances which affect the input and output measurements and in the case of unknown initial conditions.

Definition of the Markov Parameter Set Associated with the Sate Observer
In this subsection, the Markov parameter set associated with an observer-based state-space model of a mechanical system is introduced.The set of the observer Markov parameters can be calculated in the discrete time domain by recursively solving the observer state-space model given by Equation ( 12) leading to the following set of equations: where the constant matrices Ȳk that appear in this sequential set are referred to as the Markov parameters of the observed system.In general, the set of the observer Markov parameters Ȳk can be rewritten as: The recursive solution of the discrete-time observer state-space model given by Equation ( 12) based on a zero set of initial conditions leads to the following equations: Furthermore, by using the definition of the observer Markov parameter given by Equation ( 14), one can rewrite the input-output sequence description associated with a linear dynamical system as follows: In this equation, the output vector ŷ(k) estimated by using the observed model tends to the actual measured output y(k) for a large time step k provided that the estimation error closely approaches zero.Therefore, one can assume that, for a sufficiently large number denoted with p, the observer Markov parameters Ȳk can be truncated for a sufficiently large time step k such that k > p.Therefore, the sequence of observer Markov parameters Ȳk represents a finite and compact set of rectangular matrices, which can be effectively used for computing the sequence of system Markov parameters Y k starting from input-output measurements of the vectors u(k) and y(k).

Computational Algorithm for the Markov Parameter Set
In this subsection, a computational algorithm for the calculation of the set of Markov parameters Y k associated with the mechanical system of interest is presented.For this purpose, consider the input-output sequence description given by Equation (11).Assuming a sequential set of l different discrete time steps k from k = 0 to k = l − 1, one can rewrite the input-output sequence description given by Equation (11) using a matrix notation as follows: where: In Equation ( 18), the matrix denoted with ỹ is an m × l matrix formed by the measured output vectors, the matrix denoted with Y is a matrix of dimensions m × rl, which contains the sequence of unknown Markov parameters, and the matrix denoted with U is a rl × l matrix composed of known input measurements.Therefore, Equation (18) represents the correlation between the input and output sequences of time histories described by Equation (11).Since for a sufficiently large p, it is possible to replace the input-output description given by Equation (11) with the observer input-output relationship given by Equation ( 17) with a good approximation, one can approximate the input-output description given by Equation ( 18) employing the following matrix equation: where: Therefore, a good approximation for the numerical solution of the set of Markov parameters associated with the observed mechanical system that appear in Equation ( 20) can be found in the following way: where the matrix Û+ indicates the Moore-Penrose pseudoinverse matrix of the matrix Û.Subsequently, the Markov parameters Y k can be calculated by using the following recursive equations: This set of recursive equations can be effectively used for calculating the set of Markov parameters Y k from the set of Markov parameters related to the observed system Ȳk identified employing input and output experimental measurements.

Definition of the Markov Parameter Set Associated with the Observer Gain
In this subsection, the definition of an additional set of Markov parameters called observer gain Markov parameters is provided.This additional set of Markov parameters Y 0 k is defined in a recursive fashion employing the following equations: The set of observer gain Markov parameters denoted with Y 0 k can be obtained from the identified set of Markov parameters associated with the observed state-space model Ȳk by means of the following equations: One can consider a combination of the system Markov parameters with the observer gain Markov parameters given by the following matrix: The combined set of Markov parameters encapsulated in the matrix Γ k can be effectively employed for identifying the discrete-time state-space model characterized by the structural matrices A, B, C and D associated with a given mechanical system.This approach offers several advantages.For instance, an optimal rectangular matrix associated with the state estimator G is obtained directly by using the system identification numerical procedure.Furthermore, the use of the state estimator allows for reducing the number of Markov parameters necessary for the subsequent computer implementation of this system identification numerical procedure.

State-Space System Identification Numerical Procedure
In this subsection, the system identification numerical procedure developed in the paper is discussed.The proposed numerical procedure is based on the use of the Eigensystem Realization Algorithm (ERA) in conjunction with the Observer/Kalman Filter Identification (OKID) method [70].The ERA/OKID procedure is a system identification method based on the time domain that is capable of computing the discrete-time set of the system and the observer gain matrices employing the identified set of Markov parameters included in the matrix Γ k .This algorithm can be effectively employed for performing the experimental modal analysis of linear mechanical systems such as the full-car vehicle model of interest for this investigation.For this purpose, the first step of the proposed system identification numerical procedure is the determination of the combination of the system and observer gain Markov parameters discussed in the previous subsection of the paper.Subsequently, the next step of the identification procedure is focused on the construction of the generalized Hankel matrix H(k − 1) defined as: For instance, in the case of the first time step, one has k = 1.Thus, it is possible to write: where α and β are integer numbers.By decomposing the generalized Hankel matrix denoted with H(0) employing the Singular Value Decomposition (SVD), one obtains the following matrix factorization: where the matrices R and S are formed by orthonormal column vectors, while Σ is a non-zero rectangular matrix that is given by: where: where σ i , 1 = 1, 2, . . ., n represent the identified singular values of the generalized matrix indicated with H(0).Furthermore, denoting with R n and S n the rectangular matrices composed of the first n columns of the matrices R and S, one can write the matrix H(0) and the pseudoinverse matrix H+ (0) as: By examining the magnitudes of the singular value σ i , one can estimate the order n of the model of the mechanical system under examination.To this end, one can keep in the identified dynamical model the normal modes associated with relatively high magnitudes and discard the singular values corresponding to relatively low magnitudes.On the other hand, assuming a given model order n of the linear mechanical system under examination, one needs to assemble another generalized Hankel matrix denoted with H(1) necessary for the calculation of the discrete-time structural matrices that identifies the dynamical model of the mechanical system.Therefore, one can define the generalized matrix H(1) given by: At this final stage, one can mathematically formulate in the discrete time domain a mechanical model represented in the state-space obtained by using the proposed system identification numerical procedure that is mathematically expressed by the following equations: where Â, B, and Ĉ respectively represent the identified structural matrices represented in the space of states that are associated with the identified model in the discrete time domain, D denotes the identified direct transmission matrix that is directly obtained in the previous step of the system identification numerical procedure where the computation of the Markov parameter set is carried out and Ĝ represents the identified observer gain matrix.The matrices E m and E r+m indicate constant Boolean matrices that are defined as follows: Finally, the modal factorization of the identified state matrix Â can be readily computed as follows: where Λ indicates the identified eigenvalue matrix and Ψ denotes the identified eigenvector matrix both represented in the discrete time domain.By converting the identified discrete-time state-space model to the identified continuous-time state-space model, the modal decomposition of the identified system state matrix Âc leads to the identification of the diagonal matrix Λc of dimension n by n containing the system eigenvalues, which can be used for obtaining the system natural frequencies, as well as the damping ratios of the mechanical system.

Demonstrative Example
In this section, a full-car mechanical model of a road vehicle is considered as a demonstrative example for the implementation of the system identification numerical procedure developed in this investigation.The computational method elaborated in this work is applied to the system identification problem of the vertical dynamics of a standard road vehicle.By doing so, the proposed method represents also a useful mathematical tool for identifying the modal parameters necessary for the optimal design of active and semi-active suspension systems.In order to demonstrate this fact, consider the mechanical system showed in Figure 1.
In the mechanical model of the full-car system represented in Figure 1, the following numerical data illustrated in the same figure are assumed: The full-car vehicle model of the mechanical system considered as a case study has n 2 = 7 degrees of freedom.In particular, three degrees of freedom are related to the motion of the chassis, and four degrees of freedom describe the configuration of the wheels.In order to represent the vertical dynamics of a standard road vehicle, one can assume the following vector of generalized coordinates: where x c (t) denotes the generalized coordinate vector associated with the three degrees of freedom of the vehicle chassis, x w (t) denotes the generalized coordinate vector associated with the four degrees of freedom of the vehicle wheels, z G (t) is the vertical displacement of the vehicle centroid, φ(t) represents the vehicle roll angular displacement, ϑ(t) denotes the vehicle pitch angular displacement, while z R 1 (t), z R 2 (t), z R 3 (t) and z R 4 (t) identify the vertical displacement of the four wheels of the vehicle.
Using the analytical methods of Lagrangian mechanics, the system equations of motion can be written in a compact matrix form given by: where M, R and K are respectively the system generalized mass, damping and stiffness matrices.
For the full-car vehicle system, the system structural matrices are defined as follows: where: The entries of the mass, damping and stiffness matrices M, R and K of the full-car vehicle model are respectively defined as follows: Moreover, the matrices A c , B c , C and D, which define the dynamical model of the mechanical system represented in the space of the states can be obtained employing Equation (5).By using the numerical data mentioned before, one can obtain the following set of eigenvalues of the matrix A c of the full-car vehicle model referred to the continuous-time state-space model: where i = √ −1 and the eigenvalue vector λ c has n c = 14 complex conjugate entries since the full-car vehicle system is a damped mechanical system described by a vector of n 2 = 7 generalized coordinates; therefore, the dimension of the system state vector is n = 14.

Case Study A
In this subsection, the first case study, referred to as Case Study A, is discussed.In this case study, it is assumed that the whole set of the system degrees of freedom is perturbed by an external source of energy, and the corresponding time responses are recorded.In order to consider a realistic scenario, it is assumed that only the measurements of the system generalized accelerations are available.Thus, suppose that one is capable of measuring the m = 7 accelerations corresponding to all the system degrees of freedom.Furthermore, assume that the r = 7 forces applied on the mechanical system can be recorded, as well.In this case, the matrices C a and B 2 are respectively given by: 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 In this numerical example, a time interval of interest equal to T = 30 (s) is considered, and it is assumed that the input perturbations can be modeled by white noises.The force acting on the vertical degree of freedom of the full-car system is shown in Figure 2 as an example.By means of a dynamical simulation carried out in the time domain, one can obtain the generalized displacements, velocities and accelerations of the full-car system corresponding to the external excitations.The vertical acceleration of the chassis of the vehicle is shown in Figure 3 as an example.In this numerical example, noise-free data are considered.The use of the system identification procedure applied to the output and input measurements allows for identifying the set of Markov parameters associated with the state-space observed system given by Ȳk .Subsequently, one can determine the Markov parameters associated with the actual state-  In this case, a sharp fall in singular values take places in correspondence with the identified model order of the system state-matrix, as shown in Figure 4.By using the identified singular values, a system state-space model can be determined following the identification system procedure discussed in the paper.The identified state-matrix Âc corresponding to the first case study has the following set of eigenvalues: By comparing Equations ( 52) and ( 55), one can observe that, for noise-free data, the time-domain system identification algorithm developed in this work can exactly capture all the eigenvalues of the system state-matrix.

Case Study B
In this subsection, the second case study, referred to as Case Study B, is discussed.In this second case study, only four system degrees of freedom corresponding to the vertical motion of the wheels are excited by a set of external perturbations, and the corresponding time responses are recorded.Thus, suppose that one is only capable of measuring the m = 4 accelerations of the four wheels.Furthermore, assume that only the r = 4 forces exerted on the vehicle wheels are available as experimental input measurements.In this case, the matrices C a and B 2 are given by: Consider a time span for the dynamic simulations equal to T = 30 (s), and the external input forces are assumed to be white noises as in the preceding case.The external force applied on the first wheel of the full-car system is shown in the Figure 5 as an example.The acceleration of the first wheel corresponding to the external excitations of the full-car system is shown in Figure 6 as an example.In this case, a small amount of artificial noise is used for altering the measured inputs and outputs so that the proposed system identification numerical procedure can be tested in the presence of uncontrollable disturbances.The Hankel matrix singular values obtained in Case Study B are shown in Figure 7.In this case, the transition from the relatively large singular values associated with the normal modes of the full-car vehicle system to the relatively small singular values representing spurious modes induced by the artificial noise is smoother when compared to the previous case study, as shown in Figure 7.The identified state-matrix Âc corresponding to the second case study has the following set of eigenvalues: By comparing Equations ( 52) and ( 58), it can be noticed that, for a dataset affected by artificial noise, the system identification procedure leads to small errors in the first few eigenvalues of the system state-matrix.

Summary, Conclusions and Future Directions of Research
The areas of interest for the research of the authors are the dynamics of rigid and flexible multibody systems, the numerical methods of applied system identification and the theory of nonlinear control [71][72][73][74][75][76][77][78][79][80].Consequently, the research work of the authors is finalized to the development of new methods for performing accurate analytic modeling, effective numerical parameter identification procedures employing experimental data and nonlinear control optimization algorithms for dynamic models of multibody mechanical systems [81][82][83][84][85][86][87][88][89][90].In this work, on the other hand, a system identification numerical procedure for identifying the parameters of linear mechanical systems was developed.The system identification numerical procedure is based on a computational approach called the Eigensystem Realization Algorithm (ERA) combined with a numerical algorithm named the Observer/Kalman Filter Identification (OKID) method.In order to illustrate the effectiveness of the numerical procedure developed in this work, the proposed method was used for obtaining the modal parameters of a full-car vehicle model, which represents a standard road vehicle model used for the design of active and semi-active suspension systems.Therefore, the analysis proposed in this investigation can be considered as the first step in the design of a control system for road vehicles.The numerical results found by means of numerical experiments showed a good agreement with the actual modal parameters of the full-car vehicle model even when the input-output measurements are affected by an external source of disturbances.Future investigations will be devoted to performing a systematic comparison of the system identification algorithm developed in this work with other numerical procedures available in the literature.

Figure 2 .
Figure 2. Vertical input force on the chassis: Case Study A.

2 )Figure 3 .
Figure 3. Vertical output acceleration of the chassis: Case Study A.
space mechanical system indicated with Y k and the Markov parameters associated with the observer gain denoted with Y 0 k .Employing the set of identified Markov parameters, one can assemble the generalized Hankel matrices H(0) and H(1) required for identifying the state-space model.The generalized matrix H(0) obtained in Case Study A features a set of singular values shown in Figure 4.

Figure 5 .
Figure 5. Vertical input force on the first wheel: Case Study B.

2 )Figure 6 .
Figure 6.Vertical output acceleration of the first wheel: Case Study B.