Research Article  Open Access
Numerical and Experimental Dynamic Analysis of IC Engine Test Beds Equipped with Highly Flexible Couplings
Abstract
Driveline components connected to internal combustion engines can be critically loaded by dynamic forces due to motion irregularity. In particular, flexible couplings used in engine test rig are usually subjected to high levels of torsional oscillations and timevarying torque. This could lead to premature failure of the test rig. In this work an effective methodology for the estimation of the dynamic behavior of highly flexible couplings in real operational conditions is presented in order to prevent unwanted halts. The methodology addresses a combination of numerical models and experimental measurements. In particular, two mathematical models of the engine test rig were developed: a torsional lumpedparameter model for the estimation of the torsional dynamic behavior in operative conditions and a finite element model for the estimation of the natural frequencies of the coupling. The experimental campaign addressed torsional vibration measurements in order to characterize the driveline dynamic behavior as well as validate the models. The measurements were achieved by a coderbased technique using optical sensors and zebra tapes. Eventually, the validated models were used to evaluate the effect of design modifications of the coupling elements in terms of natural frequencies (torsional and bending), torsional vibration amplitude, and power loss in the couplings.
1. Introduction
Flexible couplings enable the transmission of torque from a driver to a driven part of rotating equipment, by accommodating a certain amount of shaft misalignment. This is obtained by reducing the reaction forces due to axial, lateral, angular displacements that are usually present between the coupled shafts. Flexible couplings with torsional compliance (also known as highly flexible couplings) are used to reduce the transmission of shock loads from one shaft to another and/or to alter the elastodynamic characteristics of the driveline by controlling the natural frequencies of the rotating units. The literature is rich of research works regarding highly flexible couplings, with particular reference to installation in test bench for automotive applications, as in the work of Rabeih and Crolla [1] or Reitz et al. [2]. In particular, Gequn et al. [3] performed the dynamic analysis of the driveline—including a multibody model of the coupling—in order to estimate the influence of bench installation condition on the engine main bearing load. The dynamic behavior of driveline including flexible couplings was also challenged from experimental standpoints, trying to evaluate the degree of misalignment of the coupled shafts; the papers of Dewell and Mitchell [4], Cho and Jeong [5], and Patel and Darpe [6] present good descriptions of such a scenario. The same issue was addressed by Xu and Marangoni [7], using an analytical approach based on the Component Mode Synthesis method. CruzPeragón et al. [8] investigated the nonlinear behavior of highly flexible couplings and proposed a methodology based on numerical models and experimental measurements in order to estimate the main coupling parameters.
In this work, the dynamic analysis of the coupling elements in internal combustion (IC) engine test rigs is accounted from both the numerical and experimental standpoints. The goal of the research is to define an effective methodology which aims at foreseeing the dynamic behavior of highly flexible couplings in real operational conditions in order to prevent a number of problems, such as high level of torsional oscillations, whirling of coupling shafts, damage of driven or driver components, and catastrophic failure of couplings or shafts. An industrial application is used as operative framework. The test rig under investigation consists of a few main components: the engine, gears, a transmission shaft, a highly flexible coupling, and a break. Despite the peculiarity of the case study described, this paper aims at proving the complexity of a complete dynamic analysis performed on this kind of mechanical system. In particular, the full understanding of the dynamic behavior requires more than one modeling approach besides an experimental activity.
The paper is organized as follows: Section 2 illustrates the adopted methodology, composed of experimental testing and numerical modeling; Section 3 reports the main results of the study and the discussion; Section 4 is finally devoted to some concluding remarks.
2. Materials and Methods
2.1. Methodology Outline
The proposed methodology is schematically draft in the block diagram of Figure 1 and it accounts for both numerical and experimental activities. Firstly, an experimental campaign was carried out with the aim of quantifying and characterizing the dynamic behavior of the driveline. In particular, torsional vibration measurements were performed by a coderbased technique using highquality optical sensors and equidistantly spaced markers (zebra tape) on the rotating components. The optical sensors were mounted before and after the coupling rubber elements, at the engineside and brakeside, respectively. As a second step, a 3D finite element (FE) model of the driveline was developed in order to estimate all the natural frequencies and mode shapes of the system in the bandwidth of interest and to evaluate which modes can negatively affect the dynamic behavior of the driveline in operational conditions. The numerical model took into account the entire driveline with particular attention to the stiffness and inertia properties of the highly flexible coupling. The FE model was experimentally validated by using data acquired during the experimental campaign mentioned above. It will be shown that the early failure of the elastic coupling is due to a resonance phenomenon of the driveline excited at particular operational conditions. Thus, the numerical model was used in order to foresee the effect of a number of design modifications proposed to reduce the negative dynamic effects. The goal was to move the natural frequencies of the driveline outside the excitation range of the engine harmonics. An iterative process entailing the proposal of new design variants and their numerical modal analyses was thus performed. Eventually, a lumpedparameter (LP) model of the entire test rig (engine, driveline, and brake) was developed in order to estimate the torsional vibration of the system in operational conditions for the different design modifications suggested by the FE analysis. In the torsional model, developed in MatlabSimulink environment, a precise evaluation of the variable inertia properties and torque of the engine as well as the dynamic behavior of the driveline were included. The torsional LP model (of the current unaltered system) was experimentally assessed by comparison with the experimental measurements. Furthermore, the torsional model enables the evaluation of the power losses in the coupling in operating conditions, to be considered as a good feature in order to evaluate the effectiveness of design modifications.
The used methodology and the obtained results have a general meaning from a qualitative point of view. Thus, the adopted approach could be generalized to provide an effective procedure to obtain improvements in the dynamic behavior of IC engine test rig drivelines.
2.2. Case Study
A schematic of the test rig being studied is depicted in Figure 2. The test rig consists of a few main components. The crankshaft of the engine drives a transmission shaft by means of a twostage compound gearbox with a fixed gear ratio of about 1/3. The transmission shaft is connected to the electromagnetic brake through the highly flexible coupling. The coupling is composed of two rubber elements working in series and clumped to three metallic parts, namely, two outer flanges, which are joined with the transmission shaft (engineside) and with the brake, respectively, and a middle flange that joins the rubber elements to the outer flanges (Figure 3).
2.3. Experimental Tests: Setup and Test Protocols
An experimental campaign was carried out, with the aim of characterizing the current system dynamic behavior, determining the response signature, and detecting the source of critical problems, as proposed by Troncossi et al. [9]. In particular, torsional vibration measurements were achieved by a coderbased technique using highquality optical sensors and zebra tapes on the rotating components (Figure 4). The optical sensors were mounted on the two opposite sides of the coupling rubber elements, at the engineside and brakeside, respectively, in order to track the torsional oscillations affecting the rubber elements. Two test campaigns were performed, corresponding to two variants of the highly flexible coupling, equipped with rubber elements having different hardness properties (and thus different stiffness). The purpose was to obtain a wide experimental database useful for the validation of numerical models of the driveline and for the coupling characterization from a dynamical standpoint. For each campaign, a brand new flexible coupling was used, avoiding any damage of a previous use. The effect of ageing of the rubber elements was not taken into account, since the time to failure due to dynamic loads was so fast to exclude any influence of the rubber age.
Torsional vibration measurements on the coupling were carried out by using two optical sensors (Optel Thevon), acquiring TTL signals from zebra tapes with line width of 2 mm. The two sensors were equipped with two different probes (probe Optel Thevon MULTI TBYO 6M HM6X100 SURG, and probe Optel Thevon MULTI SLIT YO 6M HM6X80 SURG), both fixed to a stiff bracket. The zebra tapes were mounted in the two endsections of the coupling, at the brakeside and engineside (in Figure 4 the engineside tape is visible), providing the instantaneous angular speed (IAS) measurements conventionally denoted as and (IAS at the brakeside and engineside, resp.). The high number of lines per revolution (94 and 78, for the brakeside and engineside, resp.) guaranteed a suitable angle resolution in the torsional measurements. Compensation for the data distortions due to the butt joint of the tapes was performed according to the algorithm proposed by Janssens et al. [10]. In addition, the tacho signal from a phonic wheel with 8 teeth fixed to the engine crankshaft was acquired and the corresponding IAS will be referred to as in the following. The counting rate of the optical sensors was 800 MHz and the sampling frequency of the IAS signals was 10240 Hz.
The system was tested in different operational conditions, namely, runup and stationary regimes, and for different engine speeds. Runup tests were conducted in order to find out resonant bandwidths of the system for a continuous change of engine speed. In addition, stationary tests were performed to have more precise information about the natural frequencies for a number of different constant regimes. Two test campaigns were performed, corresponding to rubber elements of the highly flexible coupling having different hardness, namely, 45 Sh and 70 Sh in the Shore scale, respectively. For the sake of data reliability, three different runs were performed for each test condition. The torsional oscillations of the two coupling ends were evaluated in terms of their relative velocity, denoted as and determined as . Data from runup tests were analyzed in the Time and TimeFrequency domains, whereas data measured for stationary tests were analyzed in the Time and Order domains. Velocity signals acquired in stationary tests (each one lasting 20 s) were resampled with reference to the crankshaft angular position [11]. Therefore, signal was taken as reference and the synchronous average [12] of all the data was performed for each thermodynamic cycle (corresponding to two crankshaft revolutions). Statistical parameters of interests (e.g., RMS values and Irregularity Ratio) were then computed on the anglebased data averages. For the frequency analysis, the Fast Fourier Transform (FFT) of the signals was computed by averaging blocks of data corresponding to 20 thermodynamic cycles. In order to associate the frequency content of the acquired signals with the excitations, the order analysis is preferred to the frequency analysis [11, 13]. Since the major excitations of all the driveline components are due to the engine firing, an order tracking was performed by resampling data still referring to crankshaft tacho measure . Timefrequency analyses of the runup test signals, having a duration of 44 s, were performed by computing the Short Time Fourier Transform (STFT) [14, 15], calculating the spectra in time windows of 0.1 s (thus providing a frequency resolution of 10 Hz).
2.4. Finite Element Analysis: Model Development
A FE model was developed in order to estimate the natural frequencies and mode shapes of the driveline in the bandwidth of interest (Figure 5). The FE model accounted the transmission shaft and the highly flexible coupling. The engine and the brake were included as lumped mass elements and corresponding inertia (points () and () in Figure 5): their stiffness characteristics were neglected. Particular attention was devoted to the modeling of the transmission shaft and the rubber elements of the coupling, since they were identified as the most flexible components of the test rig. The flanges of the coupling were considered as lumped mass and corresponding inertia (points (), (), and () in Figure 5). Rigid links (namely, rigid spiders in Figure 5) were used in order to connect the lumped masses (the flanges, engine, and brake) to the 3D mesh of the transmission shaft and rubber elements; Table 1 collects the features of the 3D mesh. Constraints were added in order to represent the real boundary conditions (points (A) and (B) in Figure 5). At the engineside, the three orthogonal displacements were clamped while at the brakeside only the rotational coordinate was included.

A specific procedure was carried out in order to evaluate the input parameters to be included in the FE model. The coupling manufacturer provided the global torsional stiffness of the coupling, but the 3D FE model requires Young’s Modulus of the rubber elements. Therefore, an iterative procedure based on a static FE model of the coupling was performed in order to estimate suitable Young’s Modulus. A unitary torque was applied to the FE model of the coupling and the corresponding rotational displacement was computed. The torsional stiffness was thus calculated by means of a static analysis [15] performed with MSC.Nastran, SOL 101. Several simulations were performed by changing Young’s Modulus in order to reach the torsional stiffness value provided by the manufacturer. After a few iterations, correct Young’s Modulus was determined. The calculation was repeated three times in order to find Young’s Modulus for three different stiffness categories (i.e., 45 Sh, 60 Sh, and 70 Sh).
The 3D FE model requires the complete inertia tensor of the coupling flanges, whereas only the moment of inertia around the rotational axis was available from the component technical manual. In order to estimate the missing parameters, an experimental technique based on frequency response functions (FRFs) measurements was performed for the indirect measurement of the rigid body inertia properties; such a methodology is based on the wellknown Inertia Restrain Method, a technique suitable for a wide range of applications when the mass distribution of components or assemblies is not known (e.g., Mucchi et al. used it in medical [16] and mechanical [17] fields). This method requires that, in the FRFs, the mass line between the highest rigid body mode and the lowest flexible body mode is rather flat.
The FE model has been used in the presented research in order to simulate three design modifications (namely, MOD 1, MOD 2, and MOD 3) conceived for moving the resonances outside the band of the engine main excitation (i.e., the band of interest). The shifting in frequency of the resonances outside the band of interest and the reduction of the number of excited modes were the tasks for improving the system dynamic behavior. At the same time, a number of design constraints had to be respected, as geometrical dimensions and final weight. The used methodology is as follows. The resonances close to the lower threshold of the bandwidth were decreased in frequency, by adding mass or reducing stiffness in specific zones, depending on the mode shape involved. The resonances close to the upper threshold were increased in frequency by reducing mass or increasing stiffness of the transmission shaft and rubber elements. In Section 3.4, the results related to the three design modifications will be outlined and then assessed by the elastodynamic analysis.
2.5. Elastodynamic Analysis: Model Development
A detailed elastodynamic model has been developed in order to capture the local modes of the coupling and to simulate the working torsional behavior, proving the effectiveness of the suggested design corrections. The elastodynamic model of the driveline was focused on the torsional dynamics only through a lumpedparameter torsional model.
The main elements of the driveline are the IC engine, the highly flexible coupling, the brake, and the shafts linking them to each other. Since the focus of the model was the coupling, it was further divided into three main parts: the two halves facing the IC engine and the brake, respectively, and the middle flange connected to the halves by means of the rubber elements. These macroelements of the driveline are indeed the same that were used in the FE model in Section 2.4. A preliminary analysis of the inertia and stiffness of the linking shafts suggested substituting them with a puretorsional spring of the same stiffness, whereas the inertia was divided between the two linked elements.
Usually, a key point of engine modeling is the loss torque due to friction associated with the piston assembly [18]. The determination of this torque is quite complex, but different types of methods have been proposed in the literature (e.g., as the ones reviewed by Richardson [19]). Performance measurement of the engine by means of power indicators was investigated by Heywood [20]. Specific experimental tests with the engine dragged at constant speed and a dynamometer recording the torque or direct friction measurement with accurate but expensive devices were discussed in [21, 22]. Rezeka and Henein [23] introduced a methodology based on the measurement of instantaneous angular velocity in the flywheel, modeling the instantaneous torque losses by means of six weighted coefficients. Later contributions detailed the computation of these coefficients in presence of local deformations [24], specific operating conditions [25], or specific configurations of the engine [26]. The six coefficients proposed by Rezeka and Henein [23] were obtained by a linear regression fitting, making the evaluation for multicylinder engines difficult. The extension to a general ncylinder engine is due to CruzPeragón et al. [27] who proposed a nonlinear identification procedure of the six components of the torque losses. In this paper, an experimental computation of torque losses is out of the scope of the driveline model and it was not taken into account, but only a viscous coefficient will be used to assess the dissipative effects on the torque. The aim of the model is to provide a mathematical tool that includes the effects of parameters such as pressures, inertia, and geometry in order to adapt the model to any new configuration of the driveline.
The analysis was focused on the oscillations around steady working conditions only. In particular, the brake system was not taken into account, since the brakeside shaft is affected by negligible oscillations (see also Section 3.1, Figure 7). Thus the brakeside part of the coupling, the brake itself, and the corresponding linking shaft were considered as the fixed frame (ground). As a consequence, the mechanical system was simplified to a model with three degrees of freedom (DOFs) only, corresponding to the torsional displacements of the IC engine (), of the engineside coupling part (), and of the middle flange of the coupling (). Figure 6 shows a schematic drawing of the 3 DOFs model, where , , and are the inertia of the IC engine, coupling part, and middle flange, respectively. Connections , , and refer to the transmission shaft and the two rubber elements of the coupling, respectively. External torques acting on the single inertia are named as (). It is worth noting that rotational coordinates and correspond to the two measured locations: the crankshaft phonic wheel and engineside flange.
(a)
(b)
Inertia and were easy to compute since they coincide with the inertia of the single mechanical parts with respect to their center of gravity, namely, the half coupling and the middle flange. These values could be assumed from the FE analysis.
The computation of inertia was not trivial. It is the sum of different contributions of the internal components in the IC engine, some of them rotating with different angular speed. In order to reduce the complexity of the model, all main components—in terms of inertia—were reduced to the output axis of the engine, that is, the axis directly connected to the transmission shaft. The main parts of the engines considered within inertia are the alternator, the crankshaft, the gearboxes, the timing system, and the clutch.
Each linking element (, , or ) (Figure 6) was modeled as a spring and a viscous damper element working in parallel. The resulting viscoelastic torque () can be expressed aswhere and are the stiffness and viscous coefficient, respectively, whereas and are the relative rotation and relative speed of the linking element ends.
The characteristics of the transmission shaft (link ) are quite different from the rubber elements of the two coupling halves (links and ).
The dynamic behavior of rubber could be very complex, due to nonlinearity of the material response. Different models have been proposed in the literature. Qi and Boyce [28] proposed a constitutive model capturing the major features of the stressstrain behavior of thermoplastic polyurethanes, including nonlinear hyperelastic behavior, time dependence, hysteresis, and softening. Wei and Kukureka [29] proposed a resonance technique for determining the stiffness and damping properties of a composite or composite structure. Chen et al. [30] suggested a frequencydomain method for estimating the mass, stiffness, and damping matrices of the model of a structure, based on the extraction of normal modes from the complex modes, by means of a transformation matrix. CruzPeragón et al. [8] developed a methodology to identify the coupling characteristics to validate dynamic models of engine assemblies with flexible coupling. The method is based on static and dynamic tests, nonlinear models, and techniques for parameter identification. These methods are powerful tools for the characterization of the material and its implementation in nonlinear models of mechanical systems. In this paper the focus of the lumpedparameter model is on the frequency response of the flexible coupling, rather than a characterization of the rubber material. Since natural frequencies are not defined for nonlinear models, a linear model of the rubber coupling has been considered.
While the viscous damper model is correct for the steel material (in the elastic domain), for the rubber material an equivalent viscous coefficient was computed starting from a hysteretic damping model [31]. The resulting viscous coefficient iswhere is the relative damping value of the rubber and is the angular speed of the hysteretic loop. In this model is the speed of the IC engine cycle.
The only nonzero external torque is , acting on inertia (). is the result of the combustion cycles of the engine in each cylinder. The torque contribution of each cylinder comes from the expansion phase in the IC cycle and the inertia momentum of each part of the crankslide mechanism. With reference to the nomenclature in “Nomenclature,” the contribution of the four cylinders to the external torque is properly combined considering the relative phases:whereIt should be noted that the pressures in the combustion chambers () were estimated by experimental measurements previously performed by the engine manufacturer.
A detailed description of the formula in (3) and (4) can be found in several references about the theory of machines and mechanisms, for example, in the book of Uicker et al. [32]. The quantities in (3) and (4) without subscript index were supposed to be the same for all the cylinders, whereas the other ones differed along the combustion phase () of each cylinder with respect to a certain reference (e.g., the angular displacement of the inertia ):The superscript asterisk recalls that the angular displacement in (5) should be the absolute displacement, not just the oscillation of inertia around the equilibrium configuration.
The equations of torsional motion of the three DOFs system were arranged in matrix form:where and come from (2) and comes from (3). Equations were numerically integrated in Simulink environment. Other details about Simulink implementation can be found in [33].
3. Results
3.1. Experimental RunUp Test Results
Hereafter, the results relative to the current system (with 45 Sh rubber elements of the coupling) will be firstly discussed, with the aim of highlighting the dynamic effects that likely led to the early collapse of the rubber elements. Then, the main data resulting from the substitution of the 45 Sh rubber elements with harder ones (70 Sh) will be shown and compared with the previous ones. Due to confidentiality agreement with the industrial partner, no data related to the shaft velocities can be explicitly reported. Therefore, data relative to runup tests will be shown as normalized to maximum value.
A preliminary comparison among the data acquired in the three different runs confirmed the extreme repeatability of results, being negligible any difference. The following results correspond to the second runs performed (for both the 45 Sh and 70 Sh rubber element cases). In Figure 7 the time series of the acquired IAS are plotted as scaled to the maximum value of the crankshaft velocity . In the case of 45 Sh rubber elements (Figure 7(a)), the oscillation of the velocity of the engineside flange of the coupling is much higher and more irregular than the brakeside one (), which exhibits a very smooth trend (thus meaning that the brakeside shaft is affected by negligible oscillations). In particular, is subjected to a very high increment after the sixteenth second of the run. The analysis of this oscillation is the main tool to understand the dynamic phenomenon underlying the system behavior. To this aim, the analysis was focused on the relative velocity between the two coupling ends, , which was analyzed in the TimeFrequency domain to highlight the presence of resonant bands of the system. Figure 8 presents the STFT in the frequency range between 0 Hz and (not made explicit for confidential reason). In Figure 8(a) it can be noted that the frequency content of is dominated by order 1 of the crankshaft rotation and—to a smaller extent—by orders 0.5, 1.5, and 2, being order 0.5 associated with the engine thermodynamic cycle. Starting from the sixteenth second of the runup, the amplitude of significantly increases with a frequency content firstly associated with the crankshaft order 0.5 (16–18 s) and then dominated by order 1 for a long time interval (18–30 s). Two natural frequencies, whose value is here symbolically reported as and , are excited in these two phases, being widespread in a large bandwidth. It could be noted that the second one is excited also by the crankshaft order 1.5, with lower energy, in the interval 8–12 s. Other possible higher resonances, and , are slightly excited by the crankshaft orders 1.5 and 2 at about 18–20 s, but with low energy. The RMS value of computed for the entire duration of the runup is about 156 rpm.
(a)
(b)
The further test campaign carried out with the coupling carrying on the 70 Sh rubber elements led to significantly different results (Figure 7(b)). As expected, the stiffening effect of the harder rubber increased the system natural frequencies and brought about benefits on the dynamic behavior of the driveline: the oscillations of speed of the engineside flange were indeed significantly smaller than in the previous campaign, whereas was basically the same (Figure 7). As a consequence, relative velocity was significantly lower, as it can be seen in Figure 8(b), where the full scales of time plot and colormap are the same as in Figure 8(a) in order to highlight the important reduction. The RMS value of computed for the entire duration of the runup was 74 rpm, that is, less than the half of the previous campaign value. Some dynamic effects were still present, though resulting in speed oscillations smaller than in the case of the 45 Sh rubber elements. In particular, the signal was dominated by crankshaft order 1, which excited natural frequency at the time interval 10–14 s, corresponding to a velocity that was about 50% of the final runup velocity. Moreover, at the very beginning of the runup (for a velocity of about 30% of the maximum value), order 0.5 seemed to excite resonance . The asterisk symbol is introduced in order to have explicit reference to the natural frequency values for the 70 Sh coupling, which are obviously different with respect to the 45 Sh case.
3.2. Experimental Stationary Test Results
Stationary tests at different velocities were performed. The results of a limited, significant selection are presented and discussed in the paper and the corresponding runs are conventionally referred to as Regime A, Regime B,…, Regime F (for confidentiality reasons), where Regime A and Regime F are about 50% and 94%, respectively, of the maximum speed achieved in the runup tests and that will be kept as reference in the paper. After resampling and synchronously averaging the velocity signals based on the rotation of the crankshaft, the angle and order analyses as well as the time statistics are available. In particular, a more accurate estimation of the system natural frequencies with respect to the runup results was achieved.
Starting again from the analysis of the 45 Sh results, Figures 9(a) and 9(c) report the crankshaftanglebased trend of over five thermodynamic cycles and its orderbased spectral analysis for Regime C (limited to the first two crankshaft orders), which is about 66% of the runup maximum velocity (i.e., the velocity achieved at the twentieth second of the runup in Figure 7). The consistency with the timefrequency analysis of the runup tests appears to be evident: signal is indeed dominated by crankshaft orders 1 and 2, partly exciting the resonant frequency and , and at a lower extent by orders 0.5 and 1.5, which excite resonant frequencies and , respectively. A minor contribution of the transmission shaft orders is also somehow appreciable (corresponding to multiples of about onethird of the crankshaft order 1). It is worth noticing the presence of a low frequency resonance corresponding at about 0.1 crankshaft order. This natural frequency was excited in each stationary regime, but at a negligible level since the excitation coming from the engine does not contain the corresponding harmonics content. Therefore, the excitation of such a resonance was likely due to low frequency mechanical noise present in each test.
(a)
(b)
(c)
(d)
Table 2 reports the most significant statistical values of computed for the stationary tests, that is, the RMS value and the “Irregularity Ratio” (IR) between the peaktopeak amplitude and the mean value of the transmission shaft velocity. Since the results obtained from the three runs acquired for each condition proved very repeatable, data are reported referring to one run only (the second one acquired). The maximum and the mean absolute values of the differences among the runs in the computation of IR are +0.8% (Regime B) and 0.29%, respectively. Consistently with the timefrequency analysis, Regime C represents the most critical regime in terms of oscillations that the rubber elements are subjected to (in the case of the 45 Sh rubber elements).

From the analysis of all the data retrieved from both runup and stationary tests performed on the current test rig (with 45 Sh rubber elements), it can be concluded that five natural frequencies were likely present in the bandwidth 0–. In particular,(i) was never significantly excited;(ii) was excited by crankshaft order 0.5 at a velocity corresponding about to Regime B;(iii) was excited by both orders 1 and 1.5 depending on the crankshaft speed;(iv)the mediumhigh natural frequencies and were excited with low energy by crankshaft orders 1.5 and 2.The most important contribution to the coupling relative oscillations was provided by natural frequency , which is spread in quite a large band. Efforts to improve the dynamic response of the highly flexible coupling (and the entire driveline as a consequence) should be thus focused at moving these resonant bands away from the bandwidth that can be excited by crankshaft orders 0.5 to 2 (while not introducing, at the same time, other resonances in this bandwidth).
Table 2 also reports the statistical parameters computed for the stationary tests (second run) conducted with the 70 Sh rubber elements. The maximum and the mean absolute values of the differences among the runs in the computation of IR are −0.6% (Regime A) and 0.21%, respectively. The global decrement of the oscillation amplitudes appears evident for Regimes B to F, for which the current coupling presented its critical response. In spite of this great improvement, the new coupling response was slightly worse for Regime A. Figures 9(b) and 9(d) report the crankshaftanglebased trend of and its orderbased spectral analysis for Regime C, for a direct comparison with the analogous analysis presented for the 45 Sh case. Figure 10 reports the same quantities relative to Regime A at 70 Sh. In these figures the resonant frequencies are appreciable as quite widespread hills among the narrowband peaks corresponding to the crankshaft and transmission shaft orders. Since the harder rubber entails a higher stiffness of the coupling, it is reasonable to conclude that natural frequencies and of the 45 Sh rubber coupling (corresponding, resp., about to 0.10 and 0.45 crankshaft orders, for Regime C, Figure 9(c)) moved to the higher values and for the 70 Sh rubber coupling (corresponding about to 0.20 and 0.77 orders for Regime C, Figure 9(d), and to 0.25 and 0.96 for Regime A, Figure 10(b)). No other natural frequencies seemed to be excited for the 70 Sh rubber coupling for any regime achieved in the runup and in the stationary tests.
(a)
(b)
In order to analyze if the stiffening effect of the harder rubber induced secondary effects in other parts of the driveline, the IAS of the crankshaft is compared both for the runup and stationary tests. The colormaps in Figure 11 (sharing the same full scale) report the comparison between the STFTs of , whereas Table 3 reports the RMS values of the crankshaft velocity oscillations, , and the IR values computed for the seven stationary tests. The data analysis reveals that the oscillations of were slightly higher at lowmedium regimes in the case of the 70 Sh rubber, but to a negligible extent so that no problematic operations were induced on the whole system.

(a)
(b)
It is worth recalling that the performed experimental analysis did not permit determining the vibration modes associated with the mentioned natural frequencies. In other words, it is not possible to state that only torsional modes were excited, since it cannot be excluded that a flexural mode of the transmission shaft could induce coupled oscillations in the IAS.
3.3. Finite Element Analysis: Model Validation
Table 4 collects the elastic properties of the rubber element estimated by the FE static analysis described in Section 2.4, the inertia properties estimated by the experimental procedure, and the moment of inertia of the flanges around the rotational axis provided by the manufacturer. These data have been used as input data in the numerical modal analysis (Sol 103 in MSC.Nastran) presented in Section 2.4.
 
from catalogue of manufacturer. from experiments. 
The simulation results obtained through a numerical modal analysis (Sol 103 in MSC.Nastran) regarding 45 Sh and 70 Sh rubber configurations were compared with measurements (Sections 3.1 and 3.2). The experimental STFTs with the two rubber configurations clearly show a few resonance regions in the frequency range 0– that are collected in Tables 5(a) and 6(a). Tables 5(b) and 6(b) collect the natural frequencies estimated by the FE analysis. Regarding the 45 Sh rubber configuration, four numerical modes are very close to the resonance regions detected by the experimental STFT. The first matching concerns the first torsional mode of the driveline, where the elements in the zone between the engineside and the middle flange move out of phase with respect to the brake (3% versus 2.9% of ). The second matching regards the second torsional mode, which is characterized by a high amplitude displacement of the middle flange (14%–16% versus 14.9% of ). Finally, the last two experimental frequencies (32%–38% and 68% of ) correspond, respectively, to the 5th (36.8% of ) and the 8th (69.7% of ) numerical modes. In particular, the 5th numerical mode is a double mode (two roots at same frequency due to symmetry) and involved the rotation of the middle flange around axes and . This is the mode mainly excited in operational conditions. The remaining modes detected in the FE analysis are not present in the experimental map, since they are not excited by the engine harmonics in operational conditions due to their particular shape. Similar considerations can be done for the 70 Sh rubber configuration (Table 6). The comparison between the experimental and the numerical results shows that the mode which determines the highest peaks in the experimental maps regards the third local mode of the coupling. Thus, attention should be paid in order to move this mode far away from the excitation harmonics.
(a) 45 Sh: experimental frequencies  
 
(b) 45 Sh: numerical results  

(a) 70 Sh: experimental frequencies  
 
(b) 70 Sh: numerical results  

3.4. Finite Element Analysis: Dynamic Behavior Improvement
Experimental results clearly showed that, within the frequency band of interest, that is, 18%–60% of —where the main excitation harmonics due to the engine lie—resonances occur. Moreover, FE simulations permitted defining the corresponding mode shapes of such resonances (Section 3.3). This section presents the results of the three design modifications MOD 1, MOD 2, and MOD 3 proposed to move the resonances outside the band of interest.
The first design modification (MOD 1) regarded the 45 Sh rubber element configuration; MOD 1 addressed an increased weight of the middle flange that reduced at 17.8% the frequency of the third local mode of the coupling. Furthermore, the steel transmission shaft was replaced by a stiffer titanium shaft keeping the third torsional mode outside the band of interest. Eventually, a flywheel was introduced on the engine shaft in order to keep the second torsional frequency at low frequency. The comparison between Tables 5 and 7(a) shows that targets are successfully reached, but MOD 1 leads to a rather heavy design, which could entail problems at highspeed conditions. It is worth noting that Mode 6 remains in the band of interest after the modification, but due to its particular shape it is not excited in operational conditions.
(a) Simulation results 45 Shore MOD 1  
 
(b) Simulation results 70 Shore MOD 2  
 
(c) Simulation results 60 Shore MOD 3  

The second design modification (MOD 2) concerned the 70 Sh rubber element configuration, where the middle flange was lightened and a flywheel on the engine shaft and a stiffer transmission shaft were placed. Table 7(b) collects the resulting natural frequencies. The 2nd torsional mode (2nd TG) still remains in the band of interest: in spite of this drawback, the global benefits of this modification will be illustrated in Section 3.6. Moreover, as stated in Section 3.2, the second local mode of the coupling (2nd LC) is not excited in operational conditions.
The last design modification (MOD 3) took into consideration the use of rubber elements with 60 Sh hardness. The standard middle flange was lightened by milling some parts and by replacing steel screws with titanium ones. The transmission shaft was included in titanium and a flywheel was mounted on the engine shaft. Table 7(c) collects the natural frequencies. The third local mode of the coupling is above the band of interest threshold, whereas the second torsional mode is still inside the band of interest. However, the 60 Sh rubber has the highest relative damping value; thus the vibration amplitude at this resonance is expected to be reduced with respect to 70 Sh. Moreover, as stated in Section 3.2, the second local mode of the coupling (2nd LC) is not excited in operational conditions. The quality of these improvements can be appreciated also in Section 3.6.
3.5. Elastodynamic Analysis: Model Validation
It has to be highlighted that since in the frequency range of interest only puretorsional modes occur (as well as notexcited axial and local modes, see Section 3.3), the elastodynamic model was focused on the torsional dynamics only.
The validation of the elastodynamic model described in Section 2.5 was based on the experimental analysis results of Sections 3.1 and 3.2. It should be noted that experimental data collected in the 45 Sh rubber element configuration could not be used for the validation of such a model. In fact, those data are mainly affected by the presence of the local mode of the coupling, which cannot be predicted by the puretorsional model described above. The model validation was thus performed on data acquired with the 70 Sh rubber coupling.
The natural frequencies computed by the LP model are in good accordance with those from FE model and experimental activity. The comparison is shown in Table 8. The output parameter for both experimental data and simulation results is the Irregularity Ratio (IR), between the peaktopeak amplitude and the mean value of the transmission shaft velocity, expressed as percentage. Results in correspondence of the rotational coordinates at different speed of the engine are shown in Table 9.
(a) Experimental resonances  
 
(b) FE results  
 
(c) LP results  


Simulation results are in good accordance with the experimental data, both clearly identifying a higher oscillation at Regime A that decreases as the speed increases from A to E, and then it increases again at Regime F. The advantage of the elastodynamic model is that the causes of the increased oscillation could be easily investigated through the simulation. In particular, the analysis of the resonance frequencies of the system shows that mode shape 2 is particularly burdensome for coordinate (Figure 12). The corresponding natural frequency is quite close to the 1st crankshaft order of Regime A, thus justifying both simulation results and experimental analysis. Moreover, the same natural frequency is still close to the half of the crankshaft order of Regime F.
3.6. Elastodynamic Analysis: Simulation Results and Discussion
The elastodynamic model of the driveline was used to simulate the dynamic behavior of the system in three different design changes. The same designation of the modifications discussed in Section 3.4 is used. Results are given as Irregularity Ratio: Tables 10, 11, and 12 collect the IR for MOD 1, MOD 2, and MOD 3.



The power loss was used as a further indicator. According to standard DIN 740 the relative damping is the ratio between the power loss of one vibration cycle and elastic deformation energy. The elastic deformation energy depends on the main frequency in the oscillation spectrum of the coupling and can be easily computed, while the relative damping is usually given in the manufacturer’s catalog (it is the parameter in (2)). The power losses of one vibration cycle for different hardness of the rubber and different speeds are collected in Table 13.

Modification MOD 1 keeps torsional vibrations at low level, with a stable and limited oscillation around the reference speed. Consequently, the power loss has low values compared to the other modifications. As a drawback, the FE analysis shows that all the resonance frequencies are shifted to lower values; for example, five resonances lay in the band (0%–18% of ). Even if they are outside the frequency band of interest (18%–60% of ), it is not excluded that the test rig will be used at lower speed regimes in the future, with the consequent need of further design modifications to avoid resonance problems again. Moreover, the MOD 1 requires a sensible increase of the coupling’s mass and then a further structural load on the supports of the test rig.
Modifications MOD 2 and MOD 3 are in the opposite direction of MOD 1. They shifted the resonance frequencies at higher values out of the selected frequency bandwidth. Torsional vibrations are still acceptable: the speed oscillation is less than 5% with respect to regime. In MOD 2 the FE model shows that the 3rd local mode frequency is outside the limit of 60% of but still close (60.9%), while MOD 3 takes a little bit higher safety factor (62.2%). Comparing the damping power in Table 13, MOD 2 shows an increased value at Regimes A and B—probably due to a close resonance—which is reduced in MOD 3. These considerations led to choose MOD 3 as the optimal design improvement.
4. Conclusions
The paper presents an effective methodology based on numerical models and experimental measurements, which aims at foreseeing the dynamic behavior of highly flexible couplings of IC engine test rigs in real operational conditions. The main goal is to prevent high level of torsional oscillations, whirling of coupling shafts, and early failures of any components. A procedure to optimize the dynamic behavior of the driveline is thus proposed and applied to an industrial application, useful to show that fully understanding of the dynamic behavior of a real mechanical system requires more than one modeling approach besides an experimental activity. In the test bench being studied, the output shaft of the engine is connected to an electromechanical brake through a transmission shaft, which hosts a highly flexible coupling with rubber elements.
The results of the reported case study led to draw the following conclusions:(i)the experimental activity showed the presence of a resonance close to the working condition of the coupling, which led to an early breaking of the rubber elements of the joint;(ii)the FE model enabled the characterization of the shape of the mode determining such a resonance and the suggestion of three different design modifications to avoid the resonance in working conditions;(iii)the LP model allowed choosing the design modification less burdensome in terms of torsional vibration of the driveline.The presented methodology could be a very useful tool in prototype design and optimization, as well as identifying the origin of unwanted dynamic effects in test rig. Although these quantitative results concern a particular test rig, the used methodology and the drawn conclusions have a general meaning from the qualitative point of view. They can be thus applied to a large variety of mechanical systems and applications offering useful guidelines in order to foresee the influence of operational conditions and design modifications on vibration generation.
Nomenclature
:  Pressure in th cylinder 
:  Bore of the piston 
:  Length of the crank 
:  Mass of the piston 
COG:  Center of gravity 
:  Mass of the piston rod (p.r.) 
:  Length of the p.r. 
:  Inertia of the p.r. 
:  Distance COG and head of p.r. 
:  Distance COG and foot of p.r. 
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work has been developed within the Advanced Mechanics Laboratory (MechLav) of Ferrara Technopole, realized through the contribution of Regione EmiliaRomagnaAssessorato Attività Produttive, Sviluppo Economico, Piano telematicoPORFESR 2007–2013, Activity I.1.1.
References
 E. M. A. Rabeih and D. A. Crolla, “Coupling of driveline and body vibrations in trucks,” in Proceedings of the 1996 SAE International Truck and Bus Meeting and Exposition, vol. 1203, SAE Conf. Trans., Detroit, 1996. View at: Publisher Site  Google Scholar
 A. Reitz, J. W. Biermann, and P. Kelly, “Special test benches to investigates driveline related NVH phenomena,” 8th Aachen Colloquium, Automobile and engine Technology, 1999. View at: Google Scholar
 S. Gequn, L. Min, and W. Haiqiao, “Research on the influence of bench installation conditions on simulation of engine main bearing load,” SAE InternatIonal Journal of Engines, pp. 1885–1890, 2009. View at: Publisher Site  Google Scholar
 D. L. Dewell and L. D. Mitchell, “Detection of a misaligned disk coupling using spectrum analysis,” Journal of Vibration, Acoustics, Stress, and Reliability in Design, vol. 106, no. 1, pp. 9–16, 1984. View at: Publisher Site  Google Scholar
 H.W. Cho and M. K. Jeong, “Enhanced prediction of misalignment conditions from spectral data using feature selection and filtering,” Expert Systems with Applications, vol. 35, no. 12, pp. 451–458, 2008. View at: Publisher Site  Google Scholar
 T. H. Patel and A. K. Darpe, “Experimental investigations on vibration response of misaligned rotors,” Mechanical Systems and Signal Processing, vol. 23, no. 7, pp. 2236–2252, 2009. View at: Publisher Site  Google Scholar
 M. Xu and R. D. Marangoni, “Vibration analysis of a motorflexible couplingrotor system subject to misalignment and unbalance, part I: theoretical model and analysis,” Journal of Sound and Vibration, vol. 176, no. 5, pp. 663–679, 1994. View at: Publisher Site  Google Scholar
 F. CruzPeragón, J. M. Palomar, F. A. Díaz, and F. J. JiménezEspadafor, “Practical identification of nonlinear characteristics of elastomeric couplings in engine assemblies,” Mechanical Systems and Signal Processing, vol. 23, no. 3, pp. 922–930, 2009. View at: Publisher Site  Google Scholar
 M. Troncossi, E. Mucchi, and A. Rivola, “Torsional Vibration Analysis of a Test Rig Driveline Equipped with a Flexible Coupling,” in Proceedings of 5th European Conference of Mechanical Engineering (ECME '14, pp. 198–206, Florence (Italy, 2014. View at: Google Scholar
 K. Janssens, P. Van Vlierberghe, W. Claes, B. Peeters, T. Martens, and P. D'Hondt, “Zebra tape butt joint detection and correction algorithm for rotating shafts with torsional vibrations,” in Proceedings of the ISMA2010, Leuven, Belgium, September 2022, 2010. View at: Google Scholar
 K. R. Fyfe and E. D. S. Munck, “Analysis of computed order tracking,” Mechanical Systems and Signal Processing, vol. 11, no. 2, pp. 187–202, 1997. View at: Publisher Site  Google Scholar
 S. Gade, H. Herlufsen, H. KonstantinHansen, and N. J. Wismer, “Order Tracking Analysis,” Technical Review No. 2, Brüel & Kjær, 1995. View at: Google Scholar
 S. Braun, “The synchronous (time domain) average revisited,” Mechanical Systems and Signal Processing, vol. 25, no. 4, pp. 1087–1102, 2011. View at: Publisher Site  Google Scholar
 A. V. Oppenheim, R. W. Schafer, and J. R. Buck, DiscreteTime Signal Processing, PrenticeHall, 1999. View at: Publisher Site
 R. X. Gao and R. Yan, “From fourier transform to wavelet transform: a historical perspective,” in Wavelets  Theory and Applications for Manufacturing (Chapter 2), Springer, BerlinHeidelberg, 2011. View at: Publisher Site  Google Scholar
 E. Mucchi, G. Bottoni, and R. Di Gregorio, “Indirect measurement of the inertia properties of a knee prosthesis through a simple frequencydomain technique,” Journal of Medical Devices, vol. 3, Article ID 0445015, 2009. View at: Publisher Site  Google Scholar
 E. Mucchi, S. Fiorati, R. Di Gregorio, and G. Dalpiaz, “Determining the rigidbody inertia properties of cumbersome systems: Comparison of techniques in time and frequency domain,” Experimental Techniques, vol. 35, no. 3, pp. 36–43, 2011. View at: Publisher Site  Google Scholar
 E. Ciulli, “A review of internal combustion engine losses. part 2: studies for global evaluations,” ARCHIVE: Proceedings of the Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering 19891996 (vols 203210), vol. 207, no. 3, pp. 229–240, 1993. View at: Publisher Site  Google Scholar
 D. E. Richardson, “Review of power cylinder friction for diesel engines,” Journal of Engineering for Gas Turbines and Power, vol. 122, pp. 608–618, 2000. View at: Publisher Site  Google Scholar
 J. B. Heywood, Internal Combustion Engine Fundamentals, vol. 26, McGrawHill, New York, USA, 1988.
 K. Liu, X. J. Liu, and C. L. Gui, “Scuffing failure analysis and experimental simulation of piston ringcylinder liner,” Tribology Letters, vol. 5, no. 4, pp. 309–312, 1998. View at: Publisher Site  Google Scholar
 H. M. Uras and D. J. Patterson, “Measurement of piston and ring assembly friction instantaneous IMEP method,” SAE Paper, Article ID 830416, 1983. View at: Publisher Site  Google Scholar
 S. F. Rezeka and N. A. Henein, “A new approach to evaluate instantaneous friction and its components in internal combustion engines,” SAE Technical Papers, Article ID 840179, 1985. View at: Publisher Site  Google Scholar
 N. G. Chalhoub, H. Nehme, N. A. Henein, and W. Bryzik, “Effects of structural deformations of the crankslider mechanism on the estimation of the instantaneous engine friction torque,” Journal of Sound and Vibration, vol. 224, no. 3, pp. 489–503, 1999. View at: Publisher Site  Google Scholar
 Y. H. Zweiri, J. F. Whidborne, and L. D. Seneviratne, “Instantaneous friction components model for transient engine operation,” Proceedings of the Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering, vol. 214, no. 7, pp. 809–824, 2000. View at: Publisher Site  Google Scholar
 Y. H. Zweiri, J. F. Whidborne, and L. D. Seneviratne, “Detailed analytical model of a singlecylinder diesel engine in the crank angle domain,” Proceedings of the Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering, vol. 215, no. 11, pp. 1197–1216, 2001. View at: Publisher Site  Google Scholar
 F. CruzPeragón, J. M. Palomar, F. A. Díaz, and F. J. JiménezEspadafor, “Fast online identification of instantaneous mechanical losses in internal combustion engines,” Mechanical Systems and Signal Processing, vol. 24, no. 1, pp. 267–280, 2010. View at: Publisher Site  Google Scholar
 H. J. Qi and M. C. Boyce, “Stressstrain behavior of thermoplastic polyurethanes,” Mechanics of Materials, vol. 37, no. 8, pp. 817–839, 2005. View at: Publisher Site  Google Scholar
 C. Y. Wei and S. N. Kukureka, “Evaluation of damping and elastic properties of composites and composite structures by the resonance technique,” Journal of Materials Science, vol. 35, no. 15, pp. 3785–3792, 2000. View at: Publisher Site  Google Scholar
 S. Y. Chen, M. S. Ju, and Y. G. Tsuei, “Estimation of mass, stiffness and damping matrices from frequency response functions,” Journal of Vibration and Acoustics, vol. 118, no. 1, pp. 78–82, 1996. View at: Publisher Site  Google Scholar
 S. S. Rao, Mechanical Vibrations, Prentice Hall, New York, 2003.
 J. Uicker, G. Pennock, and J. Shigley, Theory of Machines and Mechanisms, Oxford University Press, New York, USA, 2010.
 M. Cocconcelli, A. Agazzi, E. Mucchi, G. Dalpiaz, and R. Rubini, “Dynamic analysis of coupling elements in IC engine test rigs,” in Proceedings of the ISMA2014  USD2014 Conference, pp. 1005–1018, Leuven, Belgium, September 2014. View at: Google Scholar
Copyright
Copyright © 2017 M. Cocconcelli 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.