Research Article  Open Access
LowCost Monitoring System of Sensors for Evaluating Dynamic Solicitations of Semitrailer Structure
Abstract
Analysis of the fatigue life of a semitrailer structure necessitates identification of the loads and dynamic solicitations in the structure. These forces can be introduced in computer simulation software (multibody + finite element) for analysing the response of different design solutions to them. These numerical models must be validated and some parameters need to be measured directly in a field test with real vehicles under various driving conditions. In this study, a lowcost monitoring system is developed for application to a real fleet of semitrailers. According to the definition of the numerical model, the guidance of a virtual vehicle is defined by the threedimensional kinematics of the kingpin. For characterisation of these movements, a monitoring system having a lowcost inertial measurement unit (IMU) and global positioning system (GPS) antennas is developed with different configurations to enable analysis of the best costbenefit (result accuracy) solution, and an extended Kalman filter (EKF) that characterises the kinematic guidance of the kingpin is proposed. A semitrailer was equipped with the experimental lowcost monitoring system and highprecision sensors (IMU, GPS) in order to validate the results obtained by the experimental lowcost monitoring system and the inertialextended Kalman filter developed. The validated system has applicability in the lowcost monitoring of a fleet of real vehicles.
1. Introduction
Vehicle weight reduction is one of the effective approaches for reducing fuel consumption and emissions [1]. A simulation has revealed that, for every 10% weight reduction from the weights of an average new car and light truck, their fuel consumption would reduce by 6.9% and 7.6%, respectively. In the case of a commercial vehicle, a reduction in its kerb weight corresponds to an increased payload capacity. Therefore, weight reduction has a greater impact for this type of vehicle than for other types of vehicles.
In semitrailers, the body frame represents more than 75% of the total weight; therefore, this component has the most room for improvement. Weight reduction of the body frame should be achieved while keeping the safety level of the vehicle and the ability of the vehicle and of its components to maintain their functions unchanged or while limiting the vehicle’s degradation with aging to acceptable levels.
Today, for heavyduty trucks, an adequate endurance specification is assumed when 10% of the population of vehicles produced exceed the 800,000 km without any failure [2]. Phenomena that can influence the endurance of a vehicle can be classified into the following categories: fatigue, wear, corrosion, and shocks and collisions.
Weight reduction of the body frame of a vehicle affects its fatigue life; therefore, to ensure that the aging and endurance remain unaffected by this reduction, it is necessary to establish the vehicle design on the basis of a detailed study of the fatigue life. Fatigue life evaluation necessitates an accurate identification of the loads that the body frame will be subjected to during operation and performance of optimization by the use of simulation techniques based on the Finite Element Method. The methodology for fatigue life calculation is established by means of a virtual test, as shown in Figure 1.
As shown in this figure, a real vehicle is employed for performing discretisation of the body frame by means of finite element software (such as ANSYS or ABAQUS). This chassis information is then introduced in multibody dynamics simulation software (such as ADAMS of MSC Software Co.) to develop a virtual experimental program. This program includes different driving conditions (speed, manoeuvres, obstacles, roads, etc.) with the aim of defining the load cases to be introduced in fatigue life prediction algorithms.
To determine the loads during vehicle operation, a methodology based on the simulation of virtual models using data acquired by lowcost instrumentation has been developed. The effectiveness of this methodology—which combines acquired real data and a virtual model to optimise designs—for application to other types of vehicles has been proved [3]. The main advantage of combining a virtual model and real data is the reduction in the complexity of the data acquisition system and the required number of sensors. In a vehicle such as that analysed, the loads necessary for calculation of the fatigue life are due to the actions of tractor, suspension, load, and so forth [4, 5]. This implies the need for complex instrumentation, thereby making this methodology unsuitable for application to a vehicle fleet.
The present paper presents the development and implementation of a lowcost instrumentation system that is nonintrusive and robust; this system was developed with the aim of applying the abovedescribed methodology to trailertype vehicles. A semitrailer is a vehicle that can be coupled to different tractor vehicles. Therefore, in order for the virtual model to be able to calculate the loads regardless of which tractor vehicle is connected to the semitrailer, the inputs to the virtual model are proposed to be the actions of the tractor vehicle at the point of connection.
The multibody model inputs developed for use with the methodology can be forces, acceleration, velocities, and/or positions. Given the need to develop a lowcost system to provide inputs to the model, the kinematic guidance of the connection point (i.e., kingpin) between the tractor and the semitrailer was selected as the input in the present study.
2. Materials and Methods
2.1. Related Work
The most appropriate device for measuring the velocity of the kingpin is a Global Navigation Satellite System (GNSS) receiver, because two of the outputs are the absolute speed and the course angle (as of April 2013, only the United States NAVSTAR global positioning system (GPS) and the Russian GLONASS were global operational GNSSs). However, the velocity acquired with a lowcost GNSS receiver causes two problems—a high level of noise and low precision, which worsen under the influence of external factors such as multipath and low satellite visibility [8]. Given the environment in which these vehicles operate, these problems occur frequently. Noise in the GNSS signal implies variations in the measured velocity, which in turn implies significant changes in acceleration (Figure 2). During simulation, these variations in the input of the virtual model translate into loads much higher than actual loads. Fatigue analysis based on thus obtained loads would be unreliable, and filtering of the input signal would not solve this problem. Furthermore, speed measurement of the kingpin using a GNSS receiver cannot be performed for all types of vehicles, since the kingpin is usually hidden under the box or the load, and this approach is applicable only in the case of a semitrailer or container.
In consideration of both of these problems, the lowcost system developed in the present study is equipped with an inertial filter—specifically, an inertialextended Kalman filter (iEKF)—for the estimation of velocity of the kingpin using a GNSS receiver, which can be placed at a different location. The proposed lowcost system is based on the concept of fusion of an inertial measurement unit (IMU) and GNSS receivers with lowcost sensors. Sensor fusion refers to the combination of data acquired from multiple sensors with related information, which facilitates achievement of more specific inferences than those possible by using a single independent sensor [9]. The following are the advantages of using a multisensor approach over a singlesensor approach [10]:(i)Robust operational performance: any of the sensors has the potential to provide information during unavailability malfunctioning or failure of the other sensors.(ii)Extended temporal coverage: one sensor is able to detect or measure an event at times when the other sensors are unable to.(iii)Extended spatial coverage: one sensor can cover what other sensors cannot.(iv)Increased confidence: multiple independent measurements provide information on the same event.(v)Improved detection performance: multiple independent measurements of the same event are integrated effectively.(vi)Reduced vulnerability to denial: this occurs as a result of increased dimensionality of the measurement space.
The following are the advantages of sensorfused data over singlesensor data [9]:(i)Improved estimate: if several identical sensors are used, combination of observations from these sensors may lead to an improved estimate. Furthermore, a statistical advantage is gained by the addition of independent observations, under the assumption that the data are combined in an optimal manner.(ii)Improved observation process: use of the relative placement or motion of multiple sensors may result in an improved observation process. For example, two sensors that measure the angular directions to an object can be coordinated to determine the position of the object via triangulation.(iii)Improved observability: broadening of the baseline of physical observables can result in significant improvements. If the observations are associated correctly, the combination of sensors provides a better estimation than what could be obtained by any of the sensors independently. This results in a smaller error region.
The fusion of an IMU and GNSS receiver is frequently used for accurate determination of attitude. In [11], with the aim of determining the full attitude even for a nonaccelerating vehicle, a navigation system that incorporates inertial and twoantenna carrier phase differential GPS measurements was developed. Aided inertial navigation is a wellstudied problem [12–16]. Several GPS pseudorange and Doppleraided inertial systems have been documented in the literature; see, for example, [17].
A lowcost system that fuses inertial sensing and a twoantenna GPS with an unscented Kalman filter (UKF) has been developed in [18]. Many previous works on lowcost attitude estimation [19–28] have focused on the fusion of inexpensive sensors, including microelectromechanical system (MEMS) based gyroscopes and accelerometers, magnetometers, and the GPS.
In lowcost land vehicle navigation applications, in addition to accurate position and velocity information being required, an accurate attitude needs to be provided by the navigation system [28–31]. Use of the GPS signal for attitude determination has several advantages such as the small size, low cost, lack of cumulative errors, and high accuracy of the GPS. However, attitude determination using the GPS also has drawbacks such as its susceptibility to the external environment and instability in dynamic applications [31]. The conventional method of attitude determination is to use an inertial navigation system (INS). However, the INS has a drawback in that errors accumulate over time, especially in the case of a lowcost MEMSINS. If the GPS and MEMSINS are combined effectively, the reliability and accuracy of attitude determination can be improved [32–34].
The lowcost system proposed in this work differs from conventional systems in that it entails estimation of the velocities (to be used as inputs in computer simulation models) at a point of the vehicle and not the position or attitude. In any case, given the advantages of the fusion of these sensors, estimation of other states that are inputs to safety systems is a routine methodology for a vehicle.
Several works have proposed estimation methods of several key vehicle states—sideslip angle, longitudinal velocity, roll, and grade—for example, by combining automotivegrade inertial sensors with a GPS receiver [35]. Kinematic estimators that are independent of uncertain vehicle parameters integrate inertial sensors with GPS antennas to provide high update estimates of the vehicle states and sensor biases. Use of a twoantenna GPS system would enable quantification of the effects of pitch and roll on the measurements, which have been demonstrated to be quite significant in sideslip angle estimation.
The sideslip angle is an important state of vehicle lateral dynamics, and its estimation has been studied by many researchers [36–45]. For example, Bevly et al. proposed an estimation method that utilises a singleantenna GPS along with an IMU [43, 44]. This method essentially integrates the biascorrected yaw rate to obtain the vehicle heading angle, which can provide the sideslip angle along with the vehicle velocity angle measured by the singleantenna GPS. In [46], a new method was proposed that employs two lowcost singleantenna GPS receivers to estimate the vehicle sideslip angle. In the present study, four kinematics equations that relate the velocities of two different GPS receivers to the heading angle and to the longitudinal/lateral velocities of a vehicle are presented. These equations are then analysed in order to determine the theoretical limit of the sideslip/heading angle estimation. Fusion of a GPS with an IMU through an EKF is performed.
2.2. Mathematical Model
A mathematical model of a semitrailer is established in this section. This model is then used in a computer simulation to implement an EKF with a lowcost IMU and GPS receivers. According to the definition of the mathematical model, the longitudinal and lateral acceleration values (vehicle frame) and the yaw rate provided by the IMU are used as inputs. The GPS speed (reference frame) is considered as an output or measurement vector.
2.2.1. Reference and Vehicle Frames
For the definition of the vehicle dynamics model and computer simulation model, the kinematics of the kingpin should be defined as an input by characterising its speed. According to the usual approach of using a multibody system, it is necessary to define the kinematics in local coordinates in accord with the motion of the vehicle. Given that the employed sensing system includes various types of sensors, two reference systems are defined.
Reference Frame. The following is the notation of coordinates used for the reference frame (sometimes termed the local vertical, local horizontal (LVLH) frame): east, north, up (ENU). The local ENU coordinates are formed from a plane tangential to the Earth’s surface, fixed at a specific location; therefore, this plane is sometimes known as a “local tangent” or “local geodetic” plane.
Vehicle Frame. In this case, the sensitive axes of the accelerometer sensor are made to coincide with the axes of the semitrailer. The origin coincides with the location of the kingpin of the vehicle.
This frame is in accordance with the ISO standards, specifically ISO 8855. In this coordinate system (Figure 3), the forward movement of the vehicle is described along the positive axis; the lateral movement of the vehicle is described along the axis; it is positive when the vehicle is oriented to the left (from the driver’s position); and the vertical movement of the vehicle is described along the axis.
2.2.2. Extended Kalman Filter
In order to apply the proposed methodology, the equations of the virtual model must be expressed in spacestate form. The statespace representation of a physical system is composed of a set of inputs, outputs, and state variables related via firstorder differential equations, which are combined in a firstorder matrix differential equation. The state variables or states are expressed as vectors and the algebraic equations are written in matrix form, provided that it is a linear system. A general way to represent a nonlinear system with inputs, outputs, and states is as follows:Here, the first equation is the state equation and the second one is the output equation, where is the “state vector,” ; is the “input vector,” ; is the “output vector,” ; and and are random variables that represent the process noise and measurement noise, respectively.
Applying the concept of Taylor series expansion, the estimation around the current estimate can be linearised using the partial derivatives of the process and the measurement functions in order to compute estimates even in the face of nonlinear relationships. It can be assumed that the system has a state vector and that it is governed by the nonlinear stochastic difference equation. In discrete form, the equation can be expressed aswhere the measurement equation is
In this case, the nonlinear function “” is the difference equation that relates the state at the previous time step () to that at the current time step (k). The state vector is defined as follows:
Further, , which represents the inputs from the IMU (acceleration and yaw rate), is expressed as
, which represents the process and measurement noises, is given as
Finally, represents the “output vector” and . Its value is obtained via different study cases and it depends on the adopted instrumentation system and layout, as described later.
In practice, the individual values of the noises and at each time step are unknown. However, the state and measurement vectors without these values can be, respectively, approximated aswhere is some a posteriori estimate of the state (from a previous time step ).
2.2.3. Development of the Mathematical Model of Semitrailer Using iEKF
It is considered that the vehicle (semitrailer) performs a flat movement and that it is guided by the tractor unit (truck), as described in Figure 3. The implemented vehicular dynamics model takes into account the inertial effects, but it does not consider dynamic roll, suspension performance, inclination of the road, and so forth. For the formulation of the system dynamics model in the statespace form, the acceleration in the mobile system of the kingpin is considered; accordingly, the model can be expressed in terms of linear and angular velocities as follows: where and are longitudinal and lateral acceleration, respectively, at the kingpin (vehicle frame) and are the heading angle and yaw rate, respectively, and are the longitudinal and lateral velocities, respectively (vehicle frame).
Then, the velocities in the vehicle frame are expressed asIn discrete form, they can be expressed as However, in the simplified model, these inaccuracies are included in noise system. The model can be expressed in matrix form aswhere , and are white noise in the longitudinal/lateral velocities and in the heading angle, respectively.
2.2.4. Measurement Equations for the Implemented Filter (iEKF)
In order to define the best costbenefit layout of the proposed measurement system, different sensor configurations are considered in this study (Figure 4). Three lowcost GPS receivers, denoted as GPS_1, GPS_2, and GPS_KP, are considered, as shown in Figure 4.
The readings acquired from each of the receivers are the speed and path angle of each of the points at which the receivers are located. GPS_1: , , coordinates in the vehicle frame GPS_2: , , coordinates in the vehicle frame GPS_KP: , , coordinates in the vehicle frame
The measured speed for each GPS receiver is considered to represent the horizontal velocity (E, N) of each point.
The speed of each of the points in the reference frame (E, N) can be expressed in terms of the GPS measurement (Figure 5). For instance, the speed of point 1 can be expressed asFurthermore, the speed of a point of the semitrailer can be expressed in terms of the speed of another point (e.g., the kingpin). The velocity of point 1 in the vehicle frame is expressed asThe velocity of point 1 in the reference frame can be expressed asThe GPS planar velocities (east/north velocities of any GPS antenna) can be expressed as functions of the longitudinal/lateral velocities, yaw rate, and heading angle as follows:
2.2.5. Implementation of iEFK
For estimation of a system/process with nonlinear difference and measurement relationships, the linearisation and estimation equations are rewritten aswhere and are the actual state and measurement vectors, respectively; and are the approximate state and measurement vectors, respectively; is an a posteriori estimate of the state at step ; is the Jacobian matrix of partial derivatives of with respect to ; that is, is the Jacobian matrix of partial derivatives of with respect to ; that is, is the Jacobian matrix of partial derivatives of with respect to ; that is, is the Jacobian matrix of partial derivatives of with respect to ; that is,Note that, for simplicity of notation, the timestep subscript is not attached to the Jacobians , , , and , despite the fact that they are different at each time step.
iEKF Time Update Equations. The complete set of EKF equations is shown below. Note that we have substituted for to ensure consistency with the earlier “super minus” a priori notation and that we now attach the timestep subscript to the Jacobians , , , and to reinforce the notion that they are different at (and therefore must be recomputed at) each time step. The time update equations project the state and covariance estimates from the previous time step to the current time step . and are the process Jacobians at step , and is the process noise covariance at step .
iEKF Measurement Update Equations. The measurement update equations correct the state and covariance estimates using the measurement . and are the measurement Jacobians at step , and is the measurement noise covariance at step . (Note that we now attach the subscript to the Jacobians, thereby permitting their change with each measurement.) can be expressed in matrix form as follows: The matrices and are, respectively, given asTable 1 presents the study cases corresponding to the different sensor configurations analysed in this study.

The measurement vector and the related matrix for the study cases described in Table 1 are presented in Table 2.

The complete set of iEKF equations is shown in Figure 6.
2.2.6. Sensors and Covariance Matrices
In practice, the matrices of the process noise covariance and measurement noise covariance may change with each time step or measurement. However, they can be assumed as constant. The matrix can be expressed as
If the errors are considered as being independent, the covariance can be neglected; that is,
For characterisation of the vehicle’s dynamics, it is instrumented using a lowcost IMU at the point of its coupling with the truck, which is the kingpin (indicated as “KP” in Figure 3). The following readings acquired from the IMU (which includes three accelerometers and three gyroscopes) are considered: : longitudinal acceleration, vehicle frame : lateral acceleration, vehicle frame : yaw rateIf the effect of gravity is neglected and the movement plane is considered (while neglecting the effects of pitch and roll), the following expressions allow relating the dynamic behaviour of the vehicle with the IMU measurement at the KP:where , and are the biases of the longitudinal/lateral accelerometer and rate gyro, respectively, and , and are white noise in the longitudinal/lateral accelerometer and that in the rate gyro, respectively.
The three bias terms (, and ) of the IMU measurements are assumed to be constant throughout the measurement windows [46]. To evaluate this case, the following can be considered:where the matrix can also be expressed as follows:The obtained dispersion of measurements and the statistical characterisation of the lowcost IMU and GPS receivers after completion of the virtual experimental program are as presented in Table 3.

The measurement noise covariance is a diagonal matrix because all measurements are independent of each other and their size is dependent on the analysed measures, according to the study cases presented in Table 1.
3. Results of the Implementation in Real Vehicle
For application of the proposed instrumentation system to a real vehicle and for analysing the best costbenefit layout, an experimental test program was developed. To this end, a semitrailer was equipped with sensors required to apply the proposed methodology and to validate its results. Details of this test program are as follows:(i)The developed measurement system included a lowcost IMU (Figure 7(a)) and three GPS receivers (Figure 7(b)). The chosen lowcost IMU was the UM7LT orientation sensor of CH Robotics [6], and the GPS receiver was VBOX Sport by Racelogic [7].(ii)The validation system was composed of a highaccuracy IMU (ADMAGEntryLevel, GeneSys, [47]), as shown in Figure 7(c).
(a)
(b)
(c)
In order to perform the validation, the relative errors of the velocities (longitudinal and lateral) estimated in the vehicle frame were calculated using the following expressions:where and are the relative errors of the longitudinal and lateral velocities, respectively, estimated in the vehicle frame in step . The manoeuvre shown in Figure 8 was performed using the instrumented vehicle. Then, following the completion of the experiment, the errors for the different sensor configurations considered were calculated. In Figures 9–13, different plots of measured variables are shown.
The plots of relative errors (longitudinal and lateral) are shown in Figures 14 and 15 and maximum values of errors for the different sensor configurations are presented in Table 4.

4. Conclusions and Future Work
Study of the fatigue life of freight vehicles in general and of semitrailers in particular requires knowledge of the loads to which they are subjected. In the present study, dynamics simulations were developed for this purpose, in which the inputs of these simulations were the kinematic variables of the kingpin. In order to determine the states of loads to which a fleet of vehicles is subjected, a vehiclemountable lowcost instrumentation system was developed, which can be used in real movement conditions with acceptable accuracy.
This flexible instrumentation system includes a combination of sensors, that is, an IMU and GPS antennas, as well as a softwarebased extended Kalman filter for estimations in consideration of system nonlinearities.
Various sensor configurations were proposed to enhance the estimation accuracy of the proposed system, and the best sensor configuration was found to be the one corresponding to Case in Table 4.
Future work includes implementation of the developed instrumentation system in a fleet of vehicles in collaboration with a semitrailer manufacturer, with the aim of characterising the loads in service and estimating the vehicle fatigue life as the basis for process redesign and structural optimisation.
Competing Interests
The authors declare no competing interests.
Authors’ Contributions
Pablo Luque, Daniel A. Mántaras, and Aida Rodríguez conceived and designed the main localisation scheme and formulated the mathematical model. Javier G. Jalón and José L. López performed analysis of the IMU errors and collaborated for establishment of the mathematical model. Finally, Ángel Martín, Luis Castejón, Hugo Malón, and Marco Carrera designed and performed the experiments and validated the experimental results.
Acknowledgments
This work was funded by the Spanish Ministry of Economy and Competitiveness under Project no. TRA201238826 (“Diseño Óptimo de Autobuses y Semirremolques Aligerados Basado en Predicción de Vida Frente a Fatiga, Mediante Técnicas de Ensayos Virtuales y Datos Obtenidos en Tiempo Real”). The authors would like to thank LECITRAILER S.A. for providing the vehicle (semitrailer) and truck as well as the necessary facilities for performing the experiments. The authors pay tribute to their friend and colleague Marco Carrera Alegre, Principal Investigator of different projects concerning this topic, who unexpectedly died last May 20 (may he rest in peace).
References
 A. Bandivadekar, K. Bodek, L. Cheah et al., “On the road in 2035 reducing transportation's petroleum consumption and GHG emissions,” Tech. Rep. LFEE 200805 RP, MIT Laboratory for Energy and the Environment, Cambridge, Mass, USA, 2008. View at: Google Scholar
 L. Morello, L. Rossini, G. Pia, and A. Tonoli, The Automotive Body: Volume II: System Design, Springer, London, UK, 1st edition, 2011.
 D. A. Mántaras and P. Luque, “Virtual test rig to improve the design and optimisation process of the vehicle steering and suspension systems,” Vehicle System Dynamics, vol. 50, no. 10, pp. 1563–1584, 2012. View at: Publisher Site  Google Scholar
 P. Luque and D. Mántaras, “Pneumatic suspensions in semitrailers: part I. General considerations and simplified models,” International Journal of Heavy Vehicle Systems, vol. 10, no. 4, pp. 295–308, 2003. View at: Google Scholar
 P. Luque and D. A. Mántaras, “Pneumatic suspensions in semitrailers: part II. Computer simulation,” Heavy Vehicle Systems, vol. 10, no. 4, pp. 309–320, 2003. View at: Publisher Site  Google Scholar
 CH Robotics, http://www.chrobotics.com.
 VBOX Sport, http://www.vboxmotorsport.co.uk/index.php/en/products/performancemeters/vboxsport.
 L. Cong, E. Li, H. Qin, K. V. Ling, and R. Xue, “A performance improvement method for lowcost land vehicle GPS/MEMSINS attitude determination,” Sensors, vol. 15, no. 3, pp. 5722–5746, 2015. View at: Publisher Site  Google Scholar
 D. L. Hall and J. Llinas, Handbook of Multisensor Data Fusion, CRC Press, London, UK, 2nd edition, 2008.
 E. Waltz and J. Llinas, Multisensor Data Fusion, Artech House, London, UK, 1st edition, 1990.
 Y. Yang and J. A. Farrell, “Two antennas GPSaided INS for attitude determination,” IEEE Transactions on Control Systems Technology, vol. 11, no. 6, pp. 905–918, 2003. View at: Publisher Site  Google Scholar
 K. R. Britting, Inertial Navigation Systems Analysis, WileyInterscience, New York, NY, USA, 1st edition, 1971.
 J. A. Farrell and M. Barth, The Global Positioning System and Inertial Navigation, McGrawHill, New York, NY, USA, 1st edition, 1999.
 D. H. Titterton and J. L. Weston, Strapdown Inertial Navigation Technology, On Behalf of the Institution of Electrical Engineers, Peter Peregrinus, 2nd edition, 2004.
 A. Noureldin, T. B. Karamat, and J. Georgy, Fundamentals of Inertial Navigation SatelliteBased Positioning and Their Integration, Springer, London, UK, 1st edition, 2012.
 M. S. Grewal, A. P. Andrews, and C. G. Bartone, Global Navigation Satellite Systems, Inertial Navigation, and Integration, John Wiley & Sons, New Jersey, NJ, USA, 3rd edition, 2013.
 B. Parkinson and J. Spilker, Global Positioning System: Theory and Applications Vol II, American Institute of Aeronautics and Astronautics, Washington, DC, USA, 1st edition, 1996.
 Z. Wu, M. Yao, H. Ma, W. Jia, and F. Tian, “Lowcost antenna attitude estimation by fusing inertial sensing and twoantenna GPS for vehiclemounted satcomonthemove,” IEEE Transactions on Vehicular Technology, vol. 62, no. 3, pp. 1084–1096, 2013. View at: Publisher Site  Google Scholar
 A. M. Sabatini, “Quaternionbased extended Kalman filter for determining orientation by inertial and magnetic sensing,” IEEE Transactions on Biomedical Engineering, vol. 53, no. 7, pp. 1346–1356, 2006. View at: Publisher Site  Google Scholar
 D. GebreEgziabher, R. C. Hayward, and J. D. Powell, “Design of multisensor attitude determination systems,” IEEE Transactions on Aerospace and Electronic Systems, vol. 40, no. 2, pp. 627–649, 2004. View at: Publisher Site  Google Scholar
 X. Yun and E. R. Bachmann, “Design, implementation, and experimental results of a quaternionbased Kalman filter for human body motion tracking,” IEEE Transactions on Robotics, vol. 22, no. 6, pp. 1216–1227, 2006. View at: Publisher Site  Google Scholar
 Y. S. Suh, “Orientation estimation using a quaternionbased indirect Kalman filter with adaptive estimation of external acceleration,” IEEE Transactions on Instrumentation and Measurement, vol. 59, no. 12, pp. 3296–3305, 2010. View at: Publisher Site  Google Scholar
 D. Jurman, M. Jankovec, R. Kamnik, and M. Topič, “Calibration and data fusion solution for the miniature attitude and heading reference system,” Sensors and Actuators A: Physical, vol. 138, no. 2, pp. 411–420, 2007. View at: Publisher Site  Google Scholar
 H. Rehbinder and X. Hu, “Driftfree attitude estimation for accelerated rigid bodies,” Automatica, vol. 40, no. 4, pp. 653–659, 2004. View at: Publisher Site  Google Scholar
 A.L. Lee and J.H. Kim, “3dimensional pose sensor algorithm for humanoid robot,” Control Engineering Practice, vol. 18, no. 10, pp. 1173–1182, 2010. View at: Publisher Site  Google Scholar
 Y.C. Lai and S.S. Jan, “Attitude estimation based on fusion of gyroscopes and single antenna GPS for small UAVs under the influence of vibration,” GPS Solutions, vol. 15, no. 1, pp. 67–77, 2011. View at: Publisher Site  Google Scholar
 S. Hong, M. H. Lee, H.H. Chun, S.H. Kwon, and J. L. Speyer, “Observability of error states in GPS/INS integration,” IEEE Transactions on Vehicular Technology, vol. 54, no. 2, pp. 731–743, 2005. View at: Publisher Site  Google Scholar
 Z. Wu, M. Yao, H. Ma, and W. Jia, “Improving accuracy of the vehicle attitude estimation for lowcost INS/GPS integration aided by the GPSmeasured course angle,” IEEE Transactions on Intelligent Transportation Systems, vol. 14, no. 2, pp. 553–564, 2013. View at: Publisher Site  Google Scholar
 Y. Li, M. Efatmaneshnik, and A. G. Dempster, “Attitude determination by integration of MEMS inertial sensors and GPS for autonomous agriculture applications,” GPS Solutions, vol. 16, no. 1, pp. 41–52, 2012. View at: Publisher Site  Google Scholar
 U. Maeder and M. Morari, “Attitude estimation for vehicles with partial inertial measurement,” IEEE Transactions on Vehicular Technology, vol. 60, no. 4, pp. 1496–1504, 2011. View at: Publisher Site  Google Scholar
 B. Wang, L. Miao, S. Wang, and J. Shen, “An integer ambiguity resolution method for the global positioning system (GPS)based land vehicle attitude determination,” Measurement Science and Technology, vol. 20, no. 7, Article ID 075108, 2009. View at: Publisher Site  Google Scholar
 D. M. Bevly, J. Ryu, and J. C. Gerdes, “Integrating INS sensors with GPS measurements for continuous estimation of vehicle sideslip, roll, and tire cornering stiffness,” IEEE Transactions on Intelligent Transportation Systems, vol. 7, no. 4, pp. 483–493, 2006. View at: Publisher Site  Google Scholar
 D.H. Hwang, S. H. Oh, S. J. Lee, C. Park, and C. Rizos, “Design of a lowcost attitude determination GPS/INS integrated navigation system,” GPS Solutions, vol. 9, no. 4, pp. 294–311, 2005. View at: Publisher Site  Google Scholar
 Z. Wu, M. Yao, H. Ma, and W. Jia, “Lowcost attitude estimation with MIMU and twoantenna GPS for Satcomonthemove,” GPS Solutions, vol. 17, no. 1, pp. 75–87, 2013. View at: Publisher Site  Google Scholar
 M. Wada, K. S. Yoon, and H. Hashimoto, “High accuracy road vehicle state estimation using extended Kalman filter,” in Proceedings of the IEEE Intelligent Transportation Systems, pp. 282–287, Dearborn, Miss, USA, October 2000. View at: Publisher Site  Google Scholar
 A. Nishio, K. Tozu, H. Yamaguchi, K. Asano, and Y. Amano, “Development of vehicle stability control system based on vehicle sideslip angle estimation,” Tech. Rep. 2001010137, SAE, New York, NY, USA, 2001. View at: Google Scholar
 L. Li, F.Y. Wang, and Q. Zhou, “A robust observer designed for vehicle lateral motion estimation,” in Proceedings of the IEEE Intelligent Vehicles Symposium, pp. 417–422, Las Vegas, Nev, USA, June 2005. View at: Publisher Site  Google Scholar
 X. Gao and Z. Yu, “Vehicle sideslip angle estimation by using high gain observer,” in Proceedings of the 9th International Symposium on Advanced Vehicle Control, Kobe, Japan, 2008. View at: Google Scholar
 U. Kiencke and A. Daiß, “Observation of lateral vehicle dynamics,” Control Engineering Practice, vol. 5, no. 8, pp. 1145–1150, 1997. View at: Publisher Site  Google Scholar
 D. Piyabongkarn, R. Rajamani, J. A. Grogg, and J. Y. Lew, “Development and experimental evaluation of a slip angle estimator for vehicle stability control,” in Proceedings of the American Control Conference, pp. 5366–5371, Minneapolis, Minn, USA, June 2006. View at: Google Scholar
 P. Yih, J. Ryu, and J. C. Gerdes, “Vehicle state estimation using steering torque,” in Proceedings of the American Control Conference, Boston, Mass, USA, 2004. View at: Google Scholar
 R. Anderson and D. M. Bevly, “Estimation of slip angles using a model based estimator and GPS,” in Proceedings of the American Control Conference (AAC '04), pp. 2122–2127, Boston, Mass, USA, July 2004. View at: Google Scholar
 D. M. Bevly, “Global Positioning System (GPS): a lowcost velocity sensor for correcting inertial sensor errors on ground vehicles,” Journal of Dynamic Systems, Measurement and Control, vol. 126, no. 2, pp. 255–264, 2004. View at: Publisher Site  Google Scholar
 D. M. Bevly, J. C. Gerdes, and C. Wilson, “The use of GPS based velocity measurements for measurement of sideslip and wheel slip,” Vehicle System Dynamics, vol. 38, no. 2, pp. 127–147, 2002. View at: Publisher Site  Google Scholar
 J. Ryu, E. J. Rossetter, and J. C. Gerdes, “Vehicle sideslip and roll parameter estimation using GPS,” in Proceedings of the AVEC 6th International Symposium on Advanced Vehicle Control, Hiroshima, Japan, 2002. View at: Google Scholar
 J.H. Yoon and H. Peng, “Vehicle sideslip angle estimation using two singleantenna GPS receivers,” in Proceedings of the ASME Dynamic Systems and Control Conference (DSCC '10), pp. 863–870, September 2010. View at: Publisher Site  Google Scholar
 GeneSys, http://www.genesysadma.de/.
Copyright
Copyright © 2016 Pablo Luque et al. 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.