Measuring the wheel-rail forces of a roller coaster

. During a roller coaster ride, the vehicles and their passengers are exposed to multiple times the acceleration of earth’s gravity. The corresponding forces must be transmitted by the wheels to the track and the support structure. To validate the load assumptions of static and fatigue analysis, it is necessary to measure these forces. Currently, it is deemed sufﬁcient to measure only the accelerations, providing a limited view of the whole. This article presents a method to measure the actual forces and moments that act on a wheel. Since the diameters of roller coaster wheels are comparatively small (approx. 200 mm) and the loads high (approx. 30 kN), industrially available multi-axis transducers cannot be used. Therefore, a transducer design based on strain gauges was developed and successfully implemented in a roller coaster vehicle. Due to its scalability, its application is not limited to just roller coaster wheels: it can be used for all types of non-driven wheels as well.


Introduction
Roller coasters are gravity-driven amusement rides that are guided along a spatial track.This track is designed to expose the passengers and the vehicle to multiple times the acceleration of earth's gravity, which is the sensation of such a ride.The movement on this spatial curve and its design in order to produce the designated accelerations is well researched.The reader is referred to the following exemplary literature for further studies: Pfeiffer (2005), Pombo and Ambrosio (2007), Tändl (2009), andBraccesi et al. (2015).
The validation of these models is mostly based on acceleration measurements, which can be easily conducted according to ASTM F2137-16 (2016).Due to their nature, however, they can only provide a limited view of the system.The local wheel-rail forces, especially the tangential forces and the moments acting at the contact point, cannot be derived from these.Therefore, the design of a force and torque transducer to measure these components is presented in this article.

Roller coaster vehicle and wheel bogies
A roller coaster vehicle is usually composed of a vehicle body and frame, which holds the passenger's seat and the restraint system.The wheel bogies are connected to the frame via a central axle or king pins.The concept with an axle is shown in Fig. 1.
In contrast to railway vehicles, the wheels are cylindrical or concave and the track is usually made out of tubular steel pipes (Schwarzkopf, 1972).Each wheel bogie usually holds at least three different types of wheels.The running wheels carry the main load of the vehicle, whereas the side wheels guide the vehicle on the track.The upstop wheels keep the vehicle on the track during negative vertical accelerations.Variations of this similar set-up have been in use since a first vehicle concept like this was patented in the 1920s (Miller, 1921).

Measurement of wheel forces
Since all forces and moments resulting from the vehicle's motion are transferred from the wheel to the track, the measurement of these is of special interest.On the one hand, Published by Copernicus Publications on behalf of the AMA Association for Sensor Technology.
they can be used to assess the vehicle dynamics; on the other hand, the loads for the vehicle and the track can be determined.Therefore, different measurement methods have been developed and are currently in use for cars and commercial and rail vehicles.Most of these methods are based on two fundamental force measurement principles, namely the use of strain gauges and a deformation element or piezoelectric force transducers (Weiler, 1993;Weiblen and Hofmann, 1998).Representative applications for both principles are briefly presented here.
A commonly used measurement device for the wheel forces of cars is a dynamometer consisting of a rectangular array comprising four triaxial piezoelectric force transducers (Evers et al., 2002) -Fig.2.
This arrangement provides 12 force signals which can be used to calculate the normal force and the two tangential forces at the contact point as well as the three moments acting on the wheel.Piezoelectric transducers are characterised by their high linearity and ratio between range and response threshold as well as low cross-talk and hysteresis.Despite these advantages, the cost of the sensors is high and the sensor array needs a designated space, which makes it difficult to fit into small wheels (Bonfig, 1995).
A different type of piezoelectric measurement wheel is presented in Bastiaan (2018).Here piezoelectric strain sensors are incorporated into the rubber tire.Neural networks are employed to determine the relationship between the strain measurement and the wheel forces and moments.However, the estimation is not very accurate.
Rail vehicles usually also employ strain gauges fitted onto the front and backside of the wheel discs of a wheelset as described in Schwabe and Berg (2007), Magel et al. (2008) and Joch et al. (2012).The normal and the tangential force cause strain at the measurement points which depends on their magnitude, the location of the contact point and the wheel radius.A weighted summation of the strain signals is then used to determine the corresponding forces.The waviness of the signal due to the wheelset rotation has to be corrected by low-pass filtering.
Apart from this solution, there exist stationary measuring devices as described in Stephanides et al. (2008) as well.
Here, a deformation element equipped with strain gauges is installed at certain locations of the railway track that measures the wheel-rail forces, when the train passes over it.

Calibration of multi-axis force and torque transducers
The calibration provides the relationship between the p signals of the measurement points (e.g. the bridge voltage of strain gauges, or the electric charge of a piezoelectric element), arranged in the vector u, and the n acting forces and moments on the transducer, arranged in the vector k.It is clear that p needs to be greater than or equal to n in order to determine all components.The relationship can be expressed according to Weiler (1993) as follows: or The measurement matrix Q = Q ij , Q ikl comprises the matrix Q ij ∈ R n×p with the coefficients for the linear elements u j and Q ikl ∈ R n×p 2 with the coefficient for the quadratic u 2 j and mixed elements u k • u l (k = l).For an ncomponent transducer with p measuring points, the number of coefficients can be reduced by utilising the commutativity of the multiplication and is shown in Table 1.(Giesecke, 2007).
It is obvious that the number of coefficients increases when n and p increase.The required calibration tests are at least the given values in Table 1 divided by p.It is not possible to determine all coefficients uniquely if p > n.If the deformation element is stiff, the quadratic and mixed coefficients are relatively small, and for standard industrial applications they can be neglected (Giesecke, 2007).Equation (2) then simplifies to Eq. (3) and u p,1 only contains the linear elements of  Evers et al. (2002) and Bonfig (1995).the measuring signals.
The highest possible rank of Q n,p is n; therefore, the transducer needs to be calibrated with at least m = n linear independent load situations, which yields the following equation: (4) Matrix K n,m contains the applied forces and moments measured by a reference transducer; U p,m contains the signals of the measuring points.The unknown measurement matrix Q n,p can be determined by a right-hand multiplication with the Moore-Penrose inverse of U p,m .If p > n, the solution is determined by a least square fit (Knabner and Barth, 2013); for p = n it is unique.
Generally, it is desirable that most of the elements of Q n,p are equal to zero, which can be achieved by a mechanical isolation of the measuring points.This means that a force along one axis only produces a signal at one of the measuring points.The global deviations due to the summation of the signals are less, when only one signal is required to calculate the corresponding force (Giesecke, 2007).

Specification
The main focus of the transducer lies on the measurement of the normal force F z and the tangential force F y , independent of the contact point location.The transducer needs to fit into the available space of a roller coaster wheel and should easily replace an existing wheel without the need to modify the adjacent parts.The measurement of the moments M x and M z is desirable but not essential.It is especially important that the manufacturing costs are low, since a standard roller coaster car can consist of up to four wheel bogies with each holding up to six wheels.Thus at least six transducers are required to measure all forces acting on one wheel bogie.The major transducer specifications, sorted into required (R) and optional (O), are summed up in Table 2.

Concept
Out of the aforementioned requirements, the most critical constraints for the transducer are the available space in a roller coaster wheel (diameter 200 mm) in conjunction with the maximum normal force of 30 kN and the relatively low manufacturing costs.The space constraint makes it nearly impossible to fit a telemetry system and a power supply for the transducer rotating with the wheel itself.However, all roller coaster wheels have one advantage in comparison to most other vehicle wheels: they spin freely and do not need to transmit longitudinal forces onto the rails.Therefore, the transducer can be fixed to the wheel suspension and the wheel bandage rotates around it, as shown in Fig. 3.The desired manufacturing costs practically prohibit the use of triaxial piezoelectric transducers, since the cost of an array is about EUR 20 000.Therefore a solution using strain gauges is developed.
The ring-shaped space for the deformation element is determined by the bearing with an inner diameter of 140 mm.Since the contact point can shift, especially in the y direction, it makes sense to base the transducer on the measurement of shear stresses.The shear force in a beam is constant if no other force is applied on it and it is independent of its lever arm (Giesecke, 2007).In order to create shear forces in the deformation element when a force F y or F z is present at the contact point (see Fig. 3), the measurement beams for the shear strain must be aligned in the x direction.To increase the eigenfrequency of the element, a design with four beams is chosen, which is similar to a dynamometer (Fig. 2). Figure 4 shows this configuration, including all acting forces and moments at the contact point.
If the four beams were pin-jointed, no moments would be acting along the measurement beams B 1 -B 4 .Then the force and moment equilibriums are (5) y : ↑ z : • l e , (8) By measuring the eight shear forces B y and B z , it is obviously possible to calculate the forces F y and F z as well as the moments M x and M z because they are independent of the non-measured normal forces B x .F x and M y would require the measurement of these normal forces in the beam.However, the wheel spins freely, so these are small and for this application not of interest.In order to withstand the maximum normal force F z of 30 kN and to minimise hysteresis, an integrated design for the deformation element is chosen.To maximise the shear strains in the four beams, their sections are designed similarly to an I-beam according to the recommendations in Giesecke (2007).

Modelling
The CAD model of the deformation element with the locations of the strain gauges, each holding two 45 • patterns, is shown in Fig. 5.
Each of the corresponding strain gauges indexed with F (front)/B (back) and O (outside)/I (inside) are connected to a full bridge to measure eight shear strains.Disturbance forces and moments as well as the temperature strains are compensated by the connection of the bridges (Giesecke, 2007).To analyse the design of the element, a finite element analysis with ANSYS is performed.Figure 6 shows the results when a force F y or a force F z , acting along the depicted coordinate axis, is centrally applied at the contact point of the wheel.It can be observed that the measuring points are mechanically isolated.They show only a shear stress if the corresponding force is applied; i.e. if the force F y is acting, only the gauges Y show a shear strain in the corresponding plane.If F z is acting, only the gauges Z measure a shear strain.However, if the contact point is laterally shifted, a force F z will also cause shear stresses for the gauges Y .This is obvious, because a moment M x is present, which creates a shear force according to Eq. ( 8).A shift of F y in the x direction produces a moment M z and therefore only shear strains at Y .Overall it can be observed that the force F z can be calculated by using only the signals of the gauges Z (similar to Eq. 7).The force F y as well as the moments M x and M z require a different combination of the signals of the gauges Y for their calculation (Eqs.6, 8 and 10).Therefore, even with a shift of the contact point, the force F z is mechanically isolated transmitted to the measuring points from all other acting forces and moments.By changing the dimensions of the web of the I-beam profile, the resulting strains for the force F z can be easily scaled to an optimum value of 1000 µm m −1 (according to Hottinger Baldwin Messtechnik GmbH, 2008) with a small impact on the static strength of the element.Since the expected magnitude of the force F y is only one-tenth of the force F z , the corresponding strain is by its nature considerably lower and only 300 µm m −1 for gauges Y 3 and Y 4 and 160 µm m −1 for Y 1 and Y 2 , respectively.The strain could be increased by reducing the thickness of the I-beam flange.However, since this directly influences the static strength considered for the case when the force F z is applied, the size reduction is limited.
Furthermore, the FE analysis shows that the first eigenfrequency of the deformation element with no load is above 1200 Hz.
The full design of a measurement wheel with a diameter of 360 mm is shown in Fig. 7.The deformation element with the strain gauges in Fig. 7 is oriented in the same relation to the coordinate system as shown in Fig. 5.
The outer wheel consists of an aluminium rim with a vulcanised polyurethane bandage.It is fixed with a ring onto the middle wheel that holds the bearings, so that the outer wheel can be easily replaced.The bearings in O-arrangement are contactless sealed with a self-made labyrinth seal to minimise friction.The inner ring of the bearing is fitted onto the deformation element with the strain gauges.To prevent its rotation, a feather key is present on the upper, load-free side of the wheel hub.The same deformation element can be used used for a 200 mm diameter wheel with a thin section bearing.

Calibration
In order to determine the measurement matrix according to Sect.2.3, two basic methods are available.The first one is a theoretical calibration by using a finite element analysis and www.j-sens-sens-syst.net/7/469/2018/ the second one is an experimental calibration with an actual test rig and reference force transducers.

Theoretical calibration
For the theoretical calibration, the geometry of the deformation element without the other parts of the wheel is analysed with a finite element analysis using different load cases to determine the average strains acting on the area of the strain gauges.The main advantage of this method is that no actual manufactured parts as well as a calibration test rig are required.Also, the geometry and the load application are ideal.Thus this method describes the best possible conditions.
As described in Sect. 4 there are p = 8 measuring points for n = 4 force and moment components.Therefore, the measurement matrix is rectangular and at least eight different load cases would be required to determine all coefficients.However, there are only n = 6 force and moment components available which can act on the transducer.Thus, the resulting system of equations is under-determined and the solutions for the measurement matrix are infinite.Furthermore, a calibration with the disturbance forces F x and M y is required.According to Sect.2.3 there are only p = n = 4 measuring points required to measure n = 4 components F y , F z , M x and M z .By employing the relationships of Eqs. ( 6), ( 7), ( 8) and ( 10), the signals of the measuring points can be reduced to This formula uses the symmetries of the deformation element and for the calibration only the measurable components F y , F z , M x and M z must be used.It must be noted that the proposed signal reduction in Eq. ( 11) cannot be realised completely electrically -only the summation of the signals Z is possible by a parallel or antiparallel connection of the strain gauge's bridge and supply voltages.
The measurement matrix Q 4,4 , given in Eq. ( 12), is calculated out of four different load cases with the finite element analysis.The rows are labelled with the component that is calculated with the sum of the corresponding column elements multiplied by the signals of u 4,1 , which are indicated at the column locations.
The absolute values of only six coefficients are significantly larger than zero.Again, the mechanical isolation of the measuring points Z regarding F z can be observed, since the last column shows only one coefficient substantially larger than zero.As expected by Eqs. ( 6), ( 8) and ( 10), the measuring points Y are not mechanically isolated regarding the measurement of the force F y .They show a signal when the force F y or a moment M z caused by F y or a moment M x caused by a y shift of the force application point of F z is present.By using the matrix, the deformation element can be tested with different load cases, load application points and disturbance forces.In all load cases, the forces at the contact point can be measured with deviations less than 0.0001 %, the moments with 0.05 % compared to the nominal value.This means that even under ideal conditions, there are deviations.The reasons for these are the distortion of the deformation element and the resulting change in the lever arms as well as a nonsymmetric finite element grid.

Experimental calibration
The mechanical properties of the deformation element, its manufacturing, and the positioning of the strain gauges are always subject to certain tolerances.Therefore, an experimental calibration is necessary.This is carried out with a specially designed calibration rig as shown in Fig. 8 that allows a static calibration.It basically consists of a sliding carriage on linear guides and a rail pipe console that can be shifted in the y direction.The measurement wheel is installed on the sliding carriage and the deformation element is oriented in the same relation to the coordinate system as shown in Fig. 7.The force is applied with two hydraulic cylinders and is measured with reference transducers.The moment M x is applied by shifting the rail pipe console, the moment M z by shifting the application point of the force F y .It must be noted that the application of F z and M x is close to the actual wheel-rail contact.For F y and M z a static calibration cannot provide this.These components occur in the contact area of the cylindrical wheel with the cylindrical rail due to slip, which requires a rolling or sliding contact (Knothe and Stichel, 2003).The force F y is therefore applied with a thrust pin and an eccentric pad over a larger contact area at the side of the wheel so as not to damage the polyurethane bandage.
For each component, apart from the force F z , several load cases with the maximum value in two directions (positive and negative) are carried out.For F z , only the positive maximum force is applied, since the wheel has a unilateral contact with the rail.Also, the load is applied under different angular positions, since the turning of the wheel can cause inner stresses on the deformation element.In total 121 load cases are used for the calibration and the measurement matrix is determined with the right multiplication of the Moore-Penrose inverse U in Eq. ( 4).The experimentally identified measurement matrix is given in Eq. ( 13).
The overall agreement between the finite element calibration and the experiment is good.As expected, the mechanical and electrical isolation of the measuring points is less distinct.The corresponding coefficients are still close to zero but several times larger than in Eq. ( 12).The largest deviation to the finite element analysis can be observed at the coefficients for the moments; however, the load application in the calibration rig is based on a distance measurement, which is less accurate.Especially the moment M z is introduced by a force F y where the load application point is not very distinct due to the area of the pad.

Linearity
To assess the linearity deviation, the actual values of the force and moment components, calculated by multiplying the measurement matrix by signals of the strain gauges, are compared to the target value of the reference transducers.If these are plotted together, they should ideally show a straight line through the origin with a slope of 1.The deviations from the ideal measurement are presented in Fig. 9.
The maximum deviation of the actual to the target values divided by the nominal value is defined as the linearity deviation (Dubbel et al., 1997).Depending on the manufacturing quality of the wheel, the linearity deviations for F y are less than 2 %, and for F z less than 0.6 %.These deviations are regarded as good for a non-industrial transducer with an arbitrary force application point.The deviations for the moments are higher, and range between 6 % and 18 % for M x and between 10 % and 57 % for M z .The reasons for these deviations are the inner stresses acting on the deformation element when turning the outer wheel.They cause a shift of the zero point for all strain gauge signals.This shift is only compensated for the force F y and F z , because they are calculated by an unweighted sum (refer to Eqs. 6 and 7).For the moments the sum must be weighted (refer to Eqs. 8 and 10), and so the zero shift cannot be compensated.This disturbance can be reduced if the manufacturing of the press fit for the bearing has smaller tolerances.Besides, a heat treatment of the deformation element can reduce stresses due to machining.Since the disturbance is generated by turning the wheel, a low-pass filtering can also reduce this.However, the usable frequency range of the transducer decreases, respectively.

Normally distributed deviations
To get a more global view of the deviations and to incorporate effects of the cross-talking, the deviations can also be treated as normally distributed (Parthier, 2008).To prove this, the difference between the actual and target values are calculated for all load cases and sorted according to their value and relative quantity.This creates a histogram as shown in Fig. 10.For all deviations the mean value and the standard deviation are calculated and the corresponding normal distribution is also shown in Fig. 10.
It can be observed that they are approximately normally distributed; thus, an interval with the positive and negative limits twice the standard deviation holds 95.5 % of all deviations (Parthier, 2008).These limits in relation to the nominal value are herein defined as normally distributed deviations.For F y they are less than 2 % and for F z less than 0.5 %, depending on the manufacturing quality.These values are similar to the linearity deviations.The measurement of the moments again shows higher deviations which are between 3 % and 11 % for M x and between 5 % and 36 % for M z , www.j-sens-sens-syst.net/7/469/2018/being slightly lower than the linearity deviation.Overall, it can be stated that the forces can be measured with low and acceptable deviations.The measurement of the moments has an orientating character if they are not filtered.

Reproducibility of the measurement results
The four manufactured measurement wheels are installed in a wheel bogie of an existing roller coaster train as running and side wheels (Fig. 1).Furthermore, the position of the train is captured by counting the rotations of the wheel.Several runs are carried out under similar conditions.However, it must be noted that the velocity of the train varies on each run since it is influenced by a lot of parameters, such as the temperature of all bearings and the wheel bandage as well as the current wind and weather situation.A variation of the velocity v significantly changes the load on the wheels.Its influence is quadratic since e.g. the centrifugal acceleration a of a body following a circular path with radius r is calculated according to Dankert and Dankert (2011): To determine the random deviation, the concepts of the DIN 1319-1 (1995) are employed; i.e. a fictive mean measurement is created out of five measurements by calculating the mean value of each measurement on each synchronised waypoint.The deviations to this mean measurement are calculated for all five measurements on each waypoint.Lastly, out of all deviations the standard deviation in relation to the nominal values (given in Sect.3) is determined as shown in Fig. 11.
The standard deviations are similar to the previously discussed normally distributed deviations (Sect.6.2), so only the forces can be measured with an acceptable deviation.As previously stated, the deviations of the moments are a result of inner stresses that shift the zero point if the wheel is turned.Thus, they can be reduced if the data are low-pass filtered, which is achieved here with a moving average filter with a width of 200 ms.The random deviations decrease significantly, as shown by their standard deviation in relation to the nominal value -Fig.12. Thus, by sacrificing the dynamic range of the transducer, even the moments can be measured with an acceptable random deviation.In addition, the remaining deviations are actually not random: they are caused by the change in the velocity of the train.

Results and comparison to simulations
To evaluate the measurements, they are compared to the results of a multibody simulation (MBS) of the roller coaster train on the same track.The MBS was carried out with SIM-PACK employing the FASTIM algorithm (Kalker, 1982) for the calculation of the tangential forces and material parameters for viscoelastic wheels according to Knothe and Miedler (1995).It must be noted that in the MBS an ideal track geometry without any disturbances is simulated.To dampen disturbance from the track position, the measurements are filtered with the above-mentioned moving average filter and compared to the simulated values.Representatively, Figs. 13 and 14 show the results of the simulation and the measurement for the tangential force F y on a side wheel and normal forces F z on a running wheel.
The overall agreement is good and it is especially remarkable that the FASTSIM algorithm predicts the tangential forces due to slip very well even for viscoelastic wheels that are treated herein as purely elastic.However, it must be noted that according to Braat (1993), viscoelastic effects in rolling contact are less signficant for the relevant high velocities in roller coaster applications (3 m s −1 < v < 35 m s −1 ).Furthermore, the force mostly depends on the track geometry and the resulting slip angle of the wheels, which is underlined by the high reproducibility as described in Sect.7.1.
The normal forces on the wheel are also mostly generated by the track geometry itself.They represent the main loads on the wheel.As can be observed in Fig. 14 the agreement is also good.From positions 0 to 550 m the wheel loads of the measurement are higher than the simulated ones.However, the simulated velocity is also higher than the measured ones, so part of the deviations can be explained by that fact (see also Eq. 14).For the rest of the ride the simulated velocity matches the measured one closely (deviations less than 0.2 m s −1 ).Therefore, the simulated wheel load is closer to the measured one.Also, it is easy to identify the track areas where the upstop wheel is in contact with the track pipe.These are the waypoints 15, 175, 675 and 915 m, because the measured and simulated wheel forces are close to 0 kN (i.e.< 0.2 kN).It must be noted that even the unfiltered measurements show these areas.This emphasises the fact that even if the wheel is turning at high speed the zero shift due to inner stresses is compensated by the arrangement and interconnection of the strain gauges.No force is measured if no force is present, even if the wheel spins fast.
By comparing the normal force F z of the running wheel (Fig. 14) with the tangential force F y of a side wheel (Fig. 13), it can be observed that F y is negative when F z is lower than the resulting force due to the dead weight of the vehicle (approx.5 kN).This can be easily explained, since the side wheel is perpendicularly mounted to the running wheel (Fig. 1).If the normal force F z of the running wheel is lower than the force due to the vehicle's dead weight, the wheel suspension of the running wheel rebounds and the wheel bogie with the side wheel moves in the z direction (Fig. 1).This movement creates a slip velocity along the local y axis of the wheel (Fig. 7) which results in a tangential force F y on the wheel.

Conclusions
The article proposes a simple and low-priced design for a force and torque transducer incorporated into an non-driven wheel with an outer diameter as small as 200 mm while still being able to carry and measure a maximum normal force of 30 kN.Also, it is possible to measure tangential forces in the axial direction, which are only one-tenth of the normal force.Furthermore, the application point of the forces does not need to be specified exactly.The design is easily scalable to measure different forces and can be used for all kinds of non-driven wheels.
By performing a theoretical calibration, the best possible results of the developed transducer under ideal conditions were shown.The experimental calibration showed higher deviation to the reference, but for measuring the forces these were still comparatively low (less than 2 % of the nominal value).Due to manufacturing tolerances of the transducer, the measurement of the moments acting on the wheel has only an orientating character.It can be dramatically improved by low-pass filtering the data; however, this reduces the frequency range of the transducer.Therefore, further studies should focus on improving the manufacturing quality and eliminating the zero shift that results from turning the wheel.Especially improving the circularity or changing the fit of the bearing seat could reduce these deviations significantly.One should investigate whether it is possible to measure the moments with acceptable deviations.
By employing the measurement wheel in a roller coaster vehicle, the real-life application of the transducer design was successfully demonstrated.The comparison to a multibody simulation under similar conditions showed a good agreement.Therefore, the proposed design can be used to assess the loads acting on a roller coaster wheel and thus also acting on the rail, its support structure and the adjacent vehicle frame.Furthermore, it could be possible to use the designed wheel not only for the purpose of validating the load as-sumptions, but also to integrate it permanently into the roller coaster vehicle.Then it would be possible to monitor the condition of the vehicle and track during operation.The wheel could communicate with the control system and the system could react as described in Hauer (2015).

Figure 3 .
Figure 3. Measurement wheel with all acting forces and moments at the contact point.

Figure 4 .
Figure 4. Scheme of the measurement wheel.

Figure 11 .
Figure 11.Standard deviation for five unfiltered measurements from a fictive mean value measurement: F y , F z , M x , M z .

Figure 12 .
Figure 12.Standard deviation for five filtered measurements from a fictive mean value measurement: F y , F z , M x , M z , moving average filter with 200 ms width.

Figure 13 .Figure 14 .
Figure 13.Tangential force F y on a side wheel: measurement and simulation.

Table 1 .
Number of coefficients for an n-component transducer with p measuring points