1. Introduction
The exploitation environment of equipment designed for use in special conditions (e.g., forestry industry, power industry, mining industry, or the military) is changing rapidly due to, among other things, the availability of new technologies, the need to ensure higher comfort of work, and the need to meet new emission limitations. Currently-produced engineering objects are energy efficient, comfortable, relatively easy to use and eco-friendly. However, implementation of new solutions to such objects often results in higher purchase costs. Unit maintenance cost can be balanced by sufficiently efficient and failure-free operation of such objects, which is the effect of appropriate training of the operator and the suitability of the technical parameters of the object to the local conditions of use. Another way to improve efficiency ratings is to modernize existing equipment. New and modernized equipment is expected to be of consistently high quality and sufficiently long service life. This is evident in many industries, including special high-mobility wheeled vehicles designed for off-road use. With regard to vehicles manufactured in the 1970s and 1980s, their operation potential was planned for 30 years and the mileage for 200–250 thousand kilometers. However, at present, these expectations are changing in the direction of increasing the mileage to 500–600 thousand kilometers and decreasing the service life to 20 years [
1]. The verification of assumptions about the maintenance profile demands appropriate changes in the field of vehicle production and in its equipment, with newer systems supporting the work of crews, e.g., [
2,
3,
4], as well as in the testing and verification of operating parameters [
5]. Invariably, great difficulties appear when evaluating the durability of such vehicles with respect to the target mileage [
6,
7].
Analysis and verification of the durability of a new vehicle is carried out in several stages, most often in accordance with the adopted design model, e.g., shown in
Figure 1 [
8]. Three basic stages of the verification of design assumptions can be distinguished: analysis of computational models and virtual tests carried out with the use of dedicated computer software (Final Element Method—FEM, Multibody Dynamic Simulation—MBDS) [
9,
10]; stand tests of subassemblies and, if possible, of all vehicles [
11,
12]; and prototype tests of qualification mileage in normal and contract research conditions on proving grounds and test road sections [
13].
The dynamic development of computers and software has caused vehicle manufacturers to increase the share of virtual testing in the process of preparing vehicle designs for serial manufacturing, allowing for the abandonment (or limiting) of physical testing of the prototype [
14]. Significant progress has been attained in the use of test rig examination of subassemblies and vehicle systems over the last 20 years, and even of whole vehicles within the last 10 years. Examples of test rigs for carrying out vehicle durability tests by generating vertical excitations from the road profile (tire-coupled road simulator) and placing the vehicle in a complex load state (spindle-coupled road simulator) are shown in
Figure 2.
These tests usually concern the engine and drivetrain, steering, frame, cab, seats, suspension, hydraulics, and controls. However, the performance of such tests requires the establishment of test stand control signals (obtained by the compilation of the recorded road load data acquisition—RLDA) [
5,
15,
16,
17] and the satisfaction of the test stand requirements, which limit the possibility of testing vehicles with 6 × 6 or larger chassis and over 6000 kg GVW [
11]. Moreover, such stands are very expensive; for this reason, they are used by the largest vehicle manufacturers [
18]. In the case of small-series vehicle production (i.e., special vehicles) by manufacturers who often do not have extensive testing facilities at their disposal, the estimation of vehicle durability is carried out by means of experimental tests performed during test runs over a distance sufficient to observe changes and to predict the target durability of vehicle components in the context of its mileage. Here, however, a problem arises indirectly due to the high quality of the vehicle subassemblies. To observe a visible change in the durability of subassemblies as well as in the entire vehicle where bench testing of components or entire vehicles is not possible, tests under normal conditions must be performed over increasingly longer distances. Long-distance road tests of a vehicle are time-consuming and costly if only a single vehicle is available. Therefore, ways to speed up tests performed under contract research conditions, e.g., by increasing the driving speed during the tests, choosing test road sections that generate higher loads than those generated under the conditions of normal vehicle operation, and not carrying out road tests where the level of loads does not significantly limit the durability of the vehicle component being tested [
15,
19]. Therefore, it becomes particularly important to establish reference test road sections and running speeds for mileage tests under contract research conditions in order to establish the correlation between vehicle driving conditions (load, speed, road type) and the Remaining Useful Life (RUL) of the selected key vehicle components [
20].
This paper presents a method of selection for off-road test sections in the context of generated loads at selected points of a special wheeled vehicle with 8 × 8 running gear, as well as the possibility of accelerating the tests. The method allows for determination of the model speed distributions on selected test road sections (here, off-road) and indicates which chosen off-road test sections generate the most significant loads on the key vehicle components. It is possible to estimate the possible acceleration rate of the tests in relation to the loads recorded on the reference road.
2. Vehicle Life Cycle Environmental Profile
Vehicle test conditions are specified by analyzing the stated intended vehicle operating conditions, e.g., in life cycle environmental profile (LCEP). This usually contains basic data on the operating conditions, assumed vehicle speed, cargo to be carried, and target mileage as well as indicating the types of test roads (main roads, secondary roads, off-road, cross-country) [
21,
22,
23,
24]. An example of assumptions for testing high-mobility wheeled vehicles is shown in
Table 1.
A more accurate analysis of predicted speeds and road types is possible based on an analysis of vehicle use conditions under operational conditions, e.g., [
25,
26,
27]. On the basis of the data from such analyses, it is possible to arrange the test program and to select speeds and road test sections collected in a reference measurement loop for measuring loads at selected vehicle measurement points. The results of the analysis of the recorded load runs can be compared with the loads recorded on other roads of the same class to indicate the test acceleration factor, that is, the reduction in test time (h) or mileage (km).
2.1. Speed Modeling
Vehicle speed has a significant effect on the generation of loads; hence, it should be controlled in test runs. The beta statistical distribution can be used for this purpose, which allows determination of the proportions of the speed spectrum components according to the operating profile and evaluation of the speed distribution obtained from the test runs.
The statistical distribution of beta is defined as follows:
where
f(
x) is the percentage of the adopted speed range throughout the entire vehicle speed range,
x is the value of the speed range, and
A,
B are the parameters of the distribution shape. Given the vehicle movement input data, for example on secondary roads (
Table 1) where the assumed speed is V
avg = 50 km/h, V
max = 80 km/h and the share of the test road section length per 1000 km of mileage is 25% (250 km), one can ascertain data to carry out the tests. This is shown in
Table 2 and in
Figure 3.
Small errors in determined driving time (%), i.e., 99 instead of 100, and in sectional miles (%), i.e., 102 instead of 102, are acceptable; they can be reduced by determining the A and B values more accurately. The determined speed distribution can be illustrated graphically, as shown in
Figure 3.
The received speed distribution can be treated as a reference, and the determined speed distributions obtained for selected test road sections can be compared to it. Comparison of the distributions provides hints for the driver as to, e.g., whether speeds should be increased or decreased during test drives. If the given speeds cannot be reached during the measurement runs, it may be a reason to reevaluate the road section category and reclassify it (e.g., from off-road to cross-country).
2.2. Vehicle Load
According to the applicable standards and test procedures [
23,
24], high-mobility wheeled vehicles should be loaded with a load equal to their maximum payload during testing. This verifies the durability of the vehicle in a very conservative way. In practice, a vehicle only rarely carries a full load, especially in off-road and cross-country conditions. This has a great impact on the estimated durability, as research shows that the durability of a vehicle loaded with the maximum load compared to a vehicle without load is approximately ten times lower [
28]. At the stage of selection of the test road sections, vehicle tests are carried out with the vehicle fully loaded with cargo, maintaining its correct location in the cargo space. Moreover, considering that special high-mobility wheeled vehicles are equipped with a Central Tire Inflation System (CTIS), it is important to establish and consequently maintain a constant tire pressure for the load and type of road [
29]. In addition, the driving style of the driver should be monitored, e.g., by comparing the loads on the left and right sides of the vehicle to observe whether the driver is avoiding bumps under the left wheels.
2.3. Selection of Road Test Sections
To determine whether the selected test conditions lead to the generation of greater than operating loads and to determine the possible effects of accelerated testing, an appropriate indicator is needed, with values that can be determined at the stage of qualification tests conducted by a designated testing institution when a limited set of operational data concerning the tested vehicle is available. The factor that meets these limitations may be the parameter
dp, called pseudo-damage [
30]. This is derived from the transformation of the general Basquin formula into the number of permissible load cycles
N of monotonic loading with a range
S of an element made of a material, the properties of which are represented by the coefficients
α (strength of the material to the load equal to half a cycle) and
β (slope coefficient of the strength curve of the material in the log-log diagram in the range of the number of load cycles
N from 10
3 to 10
6) [
19,
31], as shown in
Figure 4.
When generalizing Relation (2), the damage value can be included in the Palmgren–Miner linear damage accumulation formula as a quantity,
D, which should be
D ≤ 1:
Equation (3) can be expressed in a simplified form as:
where
. Then, when the purpose of the tests is to compare the effect of selected operating conditions (different road sections) on the damage rate of a vehicle component, the parameters
α and
β can be assumed to remain constant during the tests, and Equation (4) is simplified to
where:
is pseudo-damage.
It can be concluded from Relation (5) that for the evaluation of selected road sections (which have an influence on damage to the vehicle component), the load process can be recorded and then the cycles of this load can be counted using various cycle counting methods [
19,
32]. Taking into account the length of the road sections and assuming that the test conditions are fixed (i.e., the load, tire pressure, and speed distributions are constant), a normalized value can be determined for each road section
as follows:
where
is the value calculated from Relation (4) and
is the length of the test section.
By performing vehicle mileage tests on various test sections of roads of the same class (e.g., off-road), it is possible to determine which road generates more damage to a component by comparing the normalized values of pseudo damage
and
under invariable load and comparable speed distributions in the tests:
In our calculations, a preliminary value of
β = 5 was assumed. This corresponds to a high cycle strength, typical for vehicle components without finishing operations such as grinding or polishing, and without notches, typical for welded joints [
31]. Pictures that exemplify the analyzed test road sections are provided in
Figure 5. The analyzed roads were of the off-road class only.
Equation (5) applies stress–strain history to determine the fatigue life of a component. For preliminary selection of test road sections, this method allows the use of other signals containing data about load cycles, for example, an acceleration signal. However, in order to obtain a direct correlation with the fatigue life of a component, the measurement data should relate directly to the stress–strain history. Then, the choice of α and β values relates strictly to the fatigue properties of the material. In other cases, the obtained value of the pseudo-damage factor dp has no reference to the fatigue strength of the material, and only indicates the fatigue nature of the forcing source (path).
3. Results
The procedure presented for the initial verification of the test road sections was carried out using a heavy-duty high-mobility wheeled vehicle. Due to the classified nature of the tests, their description and vehicle characteristics are limited. The general appearance of the vehicle is shown in
Figure 6.
The tested vehicle was planned for an in-service system characterized as in
Table 1. Technical specifications of the vehicle are presented in
Table 3.
One of the first stages of the study was to select test road sections to perform simulated in-service loads corresponding to off-road operation, which represented 55% of the target mileage. Initial research focused on finding suitable test road sections for testing suspension components. The test sections were selected for the following reasons: ability to drive at assumed speeds, load symmetry for the left and right wheel tracks, and ability to accelerate testing. The test vehicle was equipped with a measuring system that recorded the speed on each test road section and the acceleration at selected measurement points. An example presentation of speed changes on a selected test road section is provided in
Figure 7.
Driving speeds on selected test roads were recorded in a similar manner; speed distributions were derived from these measurements, as summarized in
Figure 8.
A summary of the driving characteristics of the selected road test sections is presented in
Table 4.
An example of the distribution of loads acting on the left (aL) and right (aR) sides of the vehicle is shown in
Figure 9, and a summary is shown in
Figure 10.
The accelerations of the cab to the left are represented by the value aL. Respectively, the accelerations of the cab to the right side are represented by the value aR. If the acceleration distribution is balanced equally for the left and right sides, then the magnitude represented in
Figure 10 as Load Asymmetry has a value of 1. If the value is greater than 1, then the accelerations operating to the left side are dominant, which may be due to larger bumps under the wheels on the right side of the vehicle.
Figure 11 shows a summary of the value of the Im indicator calculated according to Equation (7), representing a multiplication of normalized loads corresponding to individual road test sections for Reference Section 1. The loads corresponded to cyclical changes in the vertical acceleration of the front axle of the vehicle which represent simultaneous deflection of the spring elements of the suspension for the left and right sides of the vehicle. This affected, among other things, the effect of the inclination of the whole vehicle around the transverse axis (i.e., pitching).
Figure 12 shows a compilation of the values of the Im indicator calculated according to Equation (7), representing a multiple of the normalized loads corresponding to the individual road test sections for Reference Section 1. The loads corresponded to the cycles of change in vertical acceleration of the sprung mass on the front axle of the vehicle.
Figure 13 presents a summary of the Im indicator values calculated according to Equation (7), representing a multiple of the normalized loads corresponding to the individual road test sections with respect to Reference Section 1. The loads corresponded to cycles of changes in the angular acceleration of the front axle of the vehicle relative to the longitudinal axis of the vehicle, causing the effect of body tilting (rolling).
4. Discussion
The data presented in
Figure 7 and
Figure 8 show that vehicle speed distributions during the test run were different from the assumptions included in the LCEP; the maximum and average speeds in the road test sections were between 50% and 70% with respect to the reference speed. Both the achieved maximum and average speeds were lower compared to the required values (see
Table 1). The achieved speeds were closest to the LCEP on the Test Section 14, as indicated in
Table 4. The other road sections did not allow the assumed speed limits to be reached. On the basis of this, a preliminary conclusion can be made that the road roughness profile places the roads between off-road and cross-country. Due to the type of road surfaces in question, the roughness profile is not constant and may change seasonally throughout the year. This is a variable factor that should be taken into account when planning vehicle durability tests.
From the chart in
Figure 9 and
Figure 10, it can be seen that the loads were asymmetric during the test drives. The load asymmetry index shows a consistent pattern wherein cab accelerations to the left have values greater than accelerations to the right. In
Figure 9, the results of the test of the symmetry of the loads acting on the left and right wheels of the vehicle are presented for the Test Sections 14 and 16 and the Reference Section 1. The data presented in this case provide an assessment of the load symmetry, which affects the equal load on the left and right side of the vehicle. It can be seen from this figure that there is load asymmetry in each of the assumed ranges in the sections shown. A summary of the determined asymmetry is shown in
Figure 10. This may be due to different bump profiles on the left and right sides of the vehicle, to the driver’s driving style, or to a combination of those factors. An effort can be made to eliminate the effect of load asymmetry by alternating runs in opposite directions on each section of the test road. If this does not correct the observed load asymmetry, it may be improved by discussion with the driver and persuading him or her not to avoid bumps acting on the left wheels of the vehicle.
From the data presented in
Figure 11, it can be seen that the highest vertical loads inducing simultaneous deflection of the left and right suspension appear on Test Sections 3, 5, 6, 12, and 14. In comparison to the loads occurring on Reference Section 1, the test acceleration on these sections is from two to three times higher. At the same time, the acceleration recorded on these test road sections do not result in an increase in the loads acting on the elements that make up the sprung mass of the front of the vehicle compared to Reference Section 1. This is advantageous because, as noted, the subject of the tests was the vehicle’s suspension system. By analyzing the data presented in
Figure 13, it can be seen that Test Sections 3, 5, 6, 8, 14, and 16 caused loads that led to angular motion of the vehicle axles. At the same time, it can be observed that Test Section 12 is better-suited for testing springs than for testing stabilizers. It can additionally be pointed out that Test Sections 10 and 11 did not generate significant loads, leading to the conclusion that they can be excluded from the testing programme.
5. Conclusions
The implementation of the method presented here should begin by determining the assumed test speeds (Vmax and Vavg) from the LCEP data. Then, control spreads should be performed using the beta distribution for the preselected road section types. In the next step, the control rides should be performed and the symmetry of the loads for the left and right paths should be checked. After acceptance of the obtained results, the recording of the loads at the chosen points can be carried out and the possible Im factor can be evaluated. The further procedure depends on the specific objectives of the tests and the decisions to be made. The presented procedure for the initial selection of road test sections allows identification of the most suitable test road sections for the planned tests. The use of the Beta statistical distribution provides a reference distribution of the test speed, which is the basis for evaluating the speeds achieved on the test road sections. The use of acceleration as a load indicator is easy to use and simple to interpret. As can be seen from the presented results, when testing the springs, the highest test acceleration was achieved Test Sections 7, 8, 10, and 11 due to their high deflection values. However, in order to take into account the loads acting on the stabilizers, Test Sections 7, 8, and 14 can be eliminated from the testing programme. Considering the speed distributions on individual test road sections, it can be seen that only Test Section 14 provides speeds close to those in the LCEP; over the other sections, speeds were lower than assumed. This leads to the general conclusion that for advanced accelerated testing specially prepared test road sections should be used, as these can cause particular forms of loading and permit the possibility of significantly accelerated testing.