Strojniški vestnik - Journal of Mechanical Engineering 59(2013)1, 14-24 © 2013 Journal of Mechanical Engineering. All rights reserved. D0l:10.5545/sv-jme.2012.755 Received for publication: 2012-08-23 Received revised form: 2012-10-06 Accepted for publication: 2012-11-16 Model Development and Dynamic Load-Sharing Analysis of Longitudinal-Connected Air Suspensions Yikai Chen15* - Jie He2 - Mark King3 - Wuwei Chen4 - Changjun Wang5 - Weihua Zhang1 1 Hefei University of Technology, School of Transportation Engineering, China 2 Southeast University, School of Transportation, China 3 Queensland University of Technology, Centre for Accident Research and Road Safety- Queensland, Australia 4 Hefei University of Technology, School of Mechanical and Automotive Engineering, China 5 Traffic Management Research Institute of the Ministry of Public Security, China The objective of this research was to investigate the effect of suspension parameters on dynamic load-sharing of longitudinal-connected air suspensions of a tri-axle semi-trailer. A novel nonlinear model of a multi-axle semi-trailer with longitudinal-connected air suspension was formulated based on fluid mechanics and thermodynamics and was validated through test results. The effects of suspension parameters on dynamic load-sharing and road-friendliness of the semi-trailer were analyzed. Simulation results indicate that the road-friendliness metric DLC (Dynamic Load Coefficient), is generally in accordance with the load-sharing metric - DLSC (Dynamic Load Sharing Coefficient). When the static height or static pressure increases, the DLSC optimization ratio declines monotonically. The effect of employing larger air lines and connectors on the DLSC optimization ratio gives varying results as road roughness increases and as driving speed increases. The results also indicate that if the air line diameter is always assumed to be larger than the connector diameter, the influence of air line diameter on load-sharing is more significant than that of the connector. Keywords: modeling, driving condition, dynamic load-sharing, longitudinal-connected, air suspension, heavy truck 0 INTRODUCTION Extensive studies of "road-friendly" heavy vehicles have been performed during the last few decades to reduce road damage and increase the rated load of vehicles. However, the load-sharing ability of multi-axle heavy vehicles, which has a strong correlation with road-friendliness, has been far from adequately investigated. Load-sharing is defined as the equalization of the axle group load across all wheels/ axles [1]. When a multi-axle heavy vehicle with leaf suspensions travels on a rough road or hits a bump/ pothole, such as a bridge-head, or speed control humps, unequally distributed loads among the axles of an axle group tend to appear due to the ineffectiveness of the load-sharing mechanism (centrally pivoted walking beam, trunnion shaft, etc.) and the high stiffness of leaf springs [2]. This phenomenon causes overloading of a single axle of the axle group, which has at least two disadvantages: (a) it increases the possibility of a tire bursting as well as reducing the maneuverability and stability of the vehicle; (b) it accelerates the rutting and fatigue that contributes to pavement damage [3]. As a consequence, the improvement of load-sharing within axle groups has attracted much attention among both vehicle manufacturers and road management departments. Load-sharing performances of axle groups are specified in regulations for road-friendly vehicles in many countries. The DIVINE (Dynamic Interaction between Vehicle and Infrastructure Experiment) project undertaken by OECD (Organization for Economic Cooperation and Development) suggests that in order to qualify as a road-friendly tandem suspension, the average load variation per unit of relative vertical suspension displacement must be less than 0.3 kN/mm [4]. The Australian specification for road-friendly suspensions nominates that road-friendly suspensions must have static load-sharing, i.e., load-sharing when the vehicle is stationary, to a defined value, between axles in an axle group or tires in an axle group. However, the formal methodology to determine the static load-sharing value on a heavy vehicle is not defined [5]. In Europe, an air suspension needs to have fully-functioning hydraulic shock absorbers to pass a static road-friendliness test [6], and heavy vehicles with road-friendly suspensions are allowed higher static axle loads. A common problem with these regulations is that only the static load-sharing of vehicles is specified and that there is no requirement for suspensions to retain their dynamic load-sharing performance, i.e., load-sharing when the vehicle is driving. Many other load-sharing metrics have also been proposed by researchers. LSC (Load-Sharing Coefficient) [7] and DLSC (Dynamic Load-Sharing Coefficient) [8] have been used to evaluate static and dynamic load-sharing, respectively. Noting that perfect load equalization would give a LSC of 1.0 [9], LSC values for steel suspensions were documented in the range 0.791 to 0.957 [7]. Air suspensions with conventional-size longitudinal air lines were placed in the middle of this range with LSCs of 0.904 to 0.925 [4]. More recent studies commissioned by the National Road Transport Commission of Australia found that installation of larger air lines on multi-axle air suspensions increased longitudinal air flow between air springs on adjacent axles [10] and [11] Follow-up tests funded by the Queensland Department of Main Roads discovered that an improvement in DLSC of 4 to 30% for a tri-axle coach and 37 to 77% for a tri-axle semi-trailer were obtained by alternating the conventional-size longitudinal air connection (three 6.5 mm inside diameter connectors connecting 6.5 mm inside diameter air lines) with a larger air connection (three 20 mm inside diameter connectors connecting 50 mm inside diameter air lines) [12] and [13]. However, due to the limitations of laboratory equipment, only vehicle speed and a limited number of air connections were considered in most tests, so that the effect of some other factors like the static absolute air pressure of the air spring and static height of the air spring on load-sharing have not yet been reported before. Limitations of laboratory and on-road tests can be addressed by developing realistic models of longitudinal-connected air suspensions. Potter et al developed a simplified tandem bogie model, and by changing the damping coefficient and torsional stiffness of the leveling beam of the model [2], it can represent load-leveling steel suspensions, independent steel suspensions, longitudinal-connected air suspensions, and independent air suspensions, respectively. Davis proposed a model of a tri-axle semitrailer with longitudinal-connected air suspensions [12], and used a variable "load-sharing fraction" to represent the load-sharing ability of the suspension. However, the physical meaning of the variable was unclear. A more realistic model of a similar tri-axle semi-trailer was developed by Roebuck et al based on aerodynamics and thermodynamics [14]. In the model, the volumetric flow rate [m3/s] between two air springs was assumed to be simply proportional to the difference in air pressure with a constant coefficient Cflow [m3/(kPa s)]; in addition, the volumes and effective areas of the air springs were simplified as constants while the vehicle was travelling. Unfortunately, these simplifications of nonlinearities reduced the precision of the proposed models. A more realistic model of longitudinal-connected multi-axle air suspensions is urgently needed for precise analysis and optimization of load-sharing in multi-axle semi-trailers. The rest of this paper is organized in the following order. In Section 1, a novel nonlinear model of longitudinal-connected tri-axle air suspensions is derived based on fluid mechanics and thermodynamics. The accuracy of the model is validated and load-sharing criteria are chosen in Section 2. Based on the model, the effects of air suspension parameters (static height and static absolute air pressure of air spring, inside diameters of air line and connector) on dynamic load-sharing are analyzed in Section 3. Finally, Section 4 presents a summary of the results and draws conclusions. 1 INTEGRATED MODEL OF VEHICLE AND ROAD EXCITATION 1.1 Mathematic Model of the Tri-Axle Semi-Trailer A basic half model representing a typical tri-axle semitrailer with longitudinal-connected air suspensions in most western countries was employed, as shown in Fig. 1. This model includes five degrees of freedom (DOF), which are vertical displacements of sprung mass and three unsprung masses, z, x1, x2, x3, as well as the pitch angle of the sprung mass, Front Rear Fig. 1. Schematic of the tri-axle semi-trailer with longitudinal-connected air suspensions The equations of motion of the semi-trailer are given by: mt1 Xj = (q1 - x1 )kt1 + c1 (z - i1 - (>l) — 1 (1) - (Pi - Po)Ai + 3 mg, mt2= (q2 - x2)kt2 + c2(z - x2) - 1 (2) - (P 2 - Po) As 2 + 3 mg, mt3x3 = (q3 - x3)kt3 + c3(z - xx3 +l) - 1 (3) - (Ps3 - Po)As3 + 3mg, J< = (Ps3 - P0) As3 - C3( Z - X3 + < - 1 mg) P - Po) Ai - Ci(z- X - <) - 3 mg) mi = (Psl - P0) Ai + (P - PoA + + (Ps3 - Po)As3 - Cl(Z - X1 -(l) - - c2(z - x2) -c3(z - x1 -(>l) - mg, l- (4) (5) where mt1, mt2, mt3 and q1, q2, q3 are the unsprung masses and road excitation of the three axles, respectively. m is the sprung mass of the semi-trailer. J is the moment of inertia of the gross sprung mass around the lateral axis. Psl, Ps2, Ps3 and Asl, As2, As3 are the dynamic absolute air pressure and dynamic effective area of the three air springs, respectively. P0 is atmospheric pressure. l is the wheelbase, c1, c2, c3 are the damping coefficients of three dampers, and ktl, kt2, kt3 are the stiffness of the three tires. 1.2 Road Roughness Excitation As the angular frequency, a>, is related to f as m = 2 nf Eq. (7) is rewritten as: Gq (®) = 4n2 rs1, Vsims10 (14) where Vsj is a function of the dynamic height of the front air spring, and is given by: Vsl = (z - x -p) Asl + Vsl (15) where Asj can also be approximated as a function of dynamic height based on experimental data. msj depends on the air flow rate inside the front connector, Gj [kg/s], is given by: -j'Gd. (16) Since only small variations of temperature, air pressure and air spring volume exist when the semitrailer is travelling, the air flow inside the front connector is considered to be an incompressible steady flow, which satisfies the following formula, according to Bernoulli's equation [23]: Pi 1 2 Pf 1 + 1 2 ôVsi =-V/1' 2 p 2 (17) where Pfj is the dynamic absolute air pressure inside the front connector; vsj, Vfj are the air flow speed [m/s] inside the front air spring and the front connector, respectively; and the air inside all the air springs, connectors and the air line is assumed to have a same constant density, p, when the semi-trailer is travelling; Afj is the effective area of the front connector, which is equal to the actual area multiplied by a contraction coefficient, 0.7 [24]. because into Eq. A Noting that As1 is related to Af1 ,1'v,1 = A lfl-vfl inserting vs1 = A V / A,1 (17) yields: V i (18) where v^ is modified with a coefficient, cd (0.8), to reflect the friction of the connector. Therefore, the actual air flow speed inside the front connector, vfj, and Gj are given by: Vi = cjp, -J i - (4) 'Ai (19) Gi = sgKPfi -Psl)cdAfi x \2p\ Pi - Pf i| A 1 - f (20) Substituting Eqs. (15), (16) and (20) into Eq. (14) yields: Pi = Pio C [(* - - P) Ai + Vsl0 ] m-1 x X [mi0 + f0' Sgn(Pf 1 - Psi )cdAf 1 X x^2p\Psi -Pf\j( 1 -(A^))dt]14. (21) The dynamic absolute air pressures inside the three connectors and the longitudinal air line are assumed to have the same value, Pfj, during travel. Thus, similar expressions for the air flow rate inside the middle and rear connectors (G2 and G3), as well as the dynamic absolute air pressure inside the middle and rear air springs (Ps2 and Ps3) are derived as follows: G2 = SgnCP/1 - Pl)CdAf 2 x x, 2p P G3 = Sgn( Pf 1 - P)CdAf 3 X (22) X, 2P|PS3 - Pf, A2 Af3 1 - (f) (23) P2 = P2OV14 [(Z - x2) As2 + Vs20 ] ' W-J X X [«20 + Sgn(Pf j - Ps 2 )cdAf 2 X A X I2p|Ps2 - Pf M (1 - (f))dt ] (24) 2 P3 = P3oOZ - X +V) A3 + Vs30 ] X [m3o + Jo' sgn( Pf ; - Ps3)cdAf 3 X -14 -14 „ m x Table 1. Parameters of the Tri-axle Semi-trailer Model A f3 xj 2p|Ps3 - Pf i|/ (1 - (Af))d'] (25) where Ps20, Vs20, ms20 are the static absolute air pressure, static volume and static air mass of the middle air spring, Ps30, Vs30, ms30 are the corresponding parameters for the rear air spring; and Af2, Af are the areas of the middle and rear connectors, respectively. The three air springs have the same static absolute air pressure, static volume and static air mass, and the three connectors have the same inside diameter. The air line is made of steel, so the volume of air line is constant. Based on Eq. (14), the dynamic absolute air pressure inside the air line, P/1, is expressed as a function of the static air mass of the air line, ml0, the static absolute air pressure inside the air line, Ps10, and the gross air flow rate inside the air line, Gh shown as follows: m,„ + G,dt p = ( ' 0 Jo ' y.4p rf 1 ^ > slO' Gl = -G1 - G2 - G3 (26) (27) Substituting Eqs. (20), (22), (23) and (27) into Eq. (26) yields: P! = m,o-1-4po (m,o +|o' {- sgn(pfl - pJ x x cdAfp\P„ -Pfy A2 j 1 - ( f - - sgn(Pf! - P2)CdAf2 x^ ¡2p\P.2 - p„i "l - ( _ As2 - - Sgn( Pfl - P, 3) CdAf 3 ji p p 3 - p y 1 -( _ As3 _ i dt)lA. (28) Based on the equations in Section 1.1, Section 1.2 and Eqs. (21), (24), (25) and (28), an integrated model of road excitation and a fully-loaded tri-axle semitrailer with longitudinal-connected air springs was developed with Matlab/ Simulink. Parts of the key parameters are tabulated in Table 1. The expression of the effective area of each air spring as a function of the dynamic height of corresponding air spring will be obtained based on test results in Section 2.2. Parameter Value Dimension Description 0.0783 m2 Static effective area of each air spring Vs10 0.0125 m3 Static volume of each air spring hso 0.16 m Static height of each air spring df 0.0065 m Inside diameter of each connector ds 0.0065 m Inside diameter of the longitudinal air line Static absolute air pressure inside Ps 10 464288 Pa each air spring, each connector and the air line Po 101325 Pa Atmosphere pressure m 8700 kg Gross sprung mass of the semitrailer ma 336 kg Unsprung mass of each air spring Moment of inertia of the gross J 5684 kg-m2 sprung mass around the lateral axis kn 1960000 N/m Stiffness of dual tires on each hub P 6.5417 kg/m3 Density of air inside air springs, air connectors and the air line Damping coefficient of each crebound 288600 N-s/m damper when dynamic height of respective suspension is increasing Damping coefficient of each cbump 184500 N-s/m damper when dynamic height of respective suspension is decreasing 2 LOAD-SHARING CRITERIA AND MODEL VALIDATION 2.1 Load-Sharing Criteria Criteria need to be chosen to evaluate the load-sharing of the semi-trailer. A metric often used to characterize the magnitude of dynamic forces of a wheel in an axle group is the DLC (Dynamic Load Coefficient) [7], defined as: DLC (i) = o. F (i)' mean V / (29) where denotes the standard deviation of wheel-force i, and Fmean(i) denotes the mean wheel-force of wheel i. Although DLC is usually referred to as a road-friendliness criterion and has been criticized for its mutually exclusivity with another load-sharing criterion, LSC [1], it still has been widely used as one measure to differentiate suspension types from each other (e.g., steel vs. air) [12], [13] and [25]. De Pont points out that LSC does not address dynamic load-sharing [8]. The DLSC was proposed as an alternative to LSC, to account for the dynamic nature of wheel-forces and instantaneous load-sharing during travel, and is defined as [8]. DLSC = j (DLSi(j) -1 fjDLSi( j))2 j=i k j=i k (30) The dynamic load-sharing of wheel i, DLSi(j), is: nFi ( j) DLSi ( j ) = - X F ( j) (31) where n is the number of wheels on one side of an axle group; k is the number of terms in the dataset; and Fi(j) is the instantaneous force at wheel i. In this study, the average DLSC of tires on the same side of the semi-trailer axle group was employed as a metric of load-sharing. The average DLC of tires on the same side of the semi-trailer axle group was used to evaluate road-friendliness as well as to analyze the relationship between load-sharing and road-friendliness. 2.2 Model Validation The prototype of the tri-axle semi-trailer was tested on various road sections for verification of the integrated model of vehicle and road excitation, as shown in Fig. 2. The tests were part of a joint project titled "Heavy vehicle suspensions - testing and analysis" between the Queensland University of Technology (QUT) and the Department of Transport and Main Roads, Queensland (TMR) [26]. The setups of the tests are shown in Fig. 2. Two types of longitudinal connections were used to connect the passive air suspensions on the same side: conventional (three 6.5 mm inside diameter connectors connecting a 6.5 mm inside diameter air line) and large (three 20 mm inside diameter connectors connecting a 50 mm inside diameter air line). Strain gauges (one per hub) were mounted on the neutral axis of each axle between the spring and the hub to record the shear force on the hubs, i.e., air spring force, and accelerometers were mounted as closely as possible to each hub and to the corresponding upper positions at the chassis to derive the height of each air spring. In addition, six air pressure transducers were employed to obtain the pressures inside the air springs, and a TRAMANCO P/L on-board CHEK-WAY telemetry system was used to record all the data. The dynamic force of each tire was derived based on the shear force on the respective hub and the acceleration on the respective axle. The effective area of each air spring was obtained by dividing the respective shear force by the respective pressure inside the air spring, and the volume of each air spring was derived by multiplying the respective effective area by the respective spring height. The tests comprised of driving the semi-trailer over three typical urban road sections at speeds ranging from 60 to 80 km/h; the sections of road varied from smooth with long undulations to rough with short undulations. The IRI (International Roughness Index) values of each road section were provided by TMR, and the IRI is related to Gq(n0) in Eq. (6) as IRI = 0.78x10^Gq(n0) [27]. Ten seconds of dynamic signal data were recorded per road section, and this was done for both experimental cases (i.e., conventional longitudinal connection vs. large longitudinal connection) for the fully loaded condition. Thus, the dynamic effective area of each air spring is approximated as a function of the dynamic height of corresponding air spring, y, based on the experimental results, shown as: A = -7.670500y3 + 2.866880y2 - -0.354226y + 0.093002. (32) The effective area multiplied by the dynamic spring height yields: Vs1 = -7.670500y4 + 2.866880y3 - -0.354226y + 0.093002y. (33) The comparisons between the test and simulation results in terms of load-sharing performance are listed in Table 2. As shown in Table 2, a reasonable agreement exists between test and simulation results for both types of connection, under various road roughness conditions and vehicle speeds. The absolute error ratio of the DLC between each test and corresponding simulation is less than 10%; except for test/simulation 4, the absolute error ratios of the DLSC are less than 20%. It is also noted that all the simulation results are smaller than the corresponding test results, which mainly dues to wear of the suspensions of the test vehicle after a period of use and some simplifications of the model of longitudinal-connected air suspensions. =i Fig. 2. Setups of the tests; a) prime mover, b) test semi-trailer with test load, c) large longitudinal air line, d) strain gauge, e) accelerometer mounted on the axle, f) air pressure transducer, g) CHEK-WAY telemetry system underneath the semi-trailer [26] Table 2. Comparison of Load-Sharing Performances between Tests and Simulations Test/ simulation number Type of longitudinal connection IRI Velocity [km/h] Load-sharing criteria Test results Simulation results Error ratio (compared with the test results) [%] 1 conventional 6.213 60 DLC 0.0791 0.0733 -7.3 DLSC 0.0505 0.0431 -14.6 2 large 6.213 60 DLC 0.0699 0.0637 -8.9 DLSC 0.0440 0.0357 -18.8 3 conventional 7.602 70 DLC 0.1034 0.1001 -3.2 DLSC 0.0851 0.0721 -15.3 4 large 7.602 70 DLC 0.0983 0.0926 -5.8 DLSC 0.0819 0.0645 -21.2 5 conventional R ftftn 80 DLC 0.1773 0.1679 -5.3 8.880 DLSC 0.1506 0.1256 -16.6 6 large ft ftftn 80 DLC 0.1775 0.1626 -8.4 8.880 DLSC 0.1474 0.1256 -14.8 It can be concluded from Table 2 that the of vehicle and road excitation in this study can be simulation results correlated well with the employed for further analysis. measurements. Therefore, the integrated model 3 EFFECT OF SUSPENSION PARAMETERS ON LOAD-SHARING 3.1 Effect of Static Height and Static Pressure The effect of the static height and static pressure of air springs for the fully-loaded semi-trailer are discussed in this section, for a constant vehicle speed of 20 m/s and a standard "B" class road profile. Note that the effective area of each air spring at a given height is nearly constant with various static absolute air pressures, i.e., Eq. (32) always holds. The new static pressure of each air spring under a new static height is related to the corresponding static spring height as follows: Pi0(new) = 2900 X 9.8 X (-7.670500^_) + + 2. S0(new) - 0.354226hs0(new) + + 0.093 0 02) 1 +1.013 25 *105. (34) 0.t —•—Type 1 --□■-Type 2 ' a) 0.15 0.16 0.17 0.18 0.19 0.2 0.21 0.22 Static height (m) 0.09 0.085 0.08 0.075^____ 0.07 0.065 0.15 0.16 0.17 0.18 0.19 0.2 0.21 0.22 b) Static height (m) 18% 16% £ 14% "Ü3 12% 2 ' is 10% .EÜ I 8% o 6% 4% 2 0.15 0.16 0.17 0.18 0.19 0.2 0.21 0.22 c) Static height (m) Fig. 3. Effect of static height on load-sharing; a) DLSC, b) DLC, c) DLSC optimization ratio and DLC The influence of the static height (varying from 0.15 to 0.22 m) and corresponding static pressure (varying from 463431 to 495204 Pa) on load-sharing are shown in Fig. 3. It can be seen in Figs. 3a and b that both the DLSC and DLC reduce as the static height increases. In fact, when the static height increases, the absolute value of dAs10 / dy increases and dVs10 / dy decreases based on Eqs. (32) and (33). These result in the reduction of dynamic stiffness of each air spring and accordingly the decline of DLSC and DLC. Another finding is that compared with the DLC of the semi-trailer with connection "2", which decreases at a relatively constant rate, the DLSC of the semitrailer with connection "2" decreases more slowly and becomes constant when the static height exceeds 0.19 m. This indicates that when a large connection is employed, the load-sharing will not change much as the static height increases and the air pressure of air springs decreases, but the dynamic tire forces will continue to reduce. The optimization ratios of both DLSC and DLSC decline as the static height increases, as shown in Fig. 3c. As the static height increases, Vs10 increases, and the volume of air line / Vs10 decreases, and the effect of employing large air line and connectors becomes less prominent. 3.2 Effect Of Inside Diameter of Air Line and Connector The effects of size of air line and connector (varying from 10 to 100 mm) on load-sharing are plotted in Fig. 4, with a constant vehicle speed of 20 m/s for the fully-loaded semi-trailer and a standard "B" class road profile. The diameters of the connectors are always less than or equal to those of the air lines in all the simulations. It can be seen in Figs. 4a and 4b that with a fixed diameter of air line, both DLSC and DLC reduce quickly as the air line diameter increases from 10 to 30 mm. For example, with a 100 mm diameter air line, reductions up to 11.4 and 8.4% are observed for DLSC and DLC, respectively. When the diameter of air line increases beyond 30 mm, both DLSC and DLC decrease more slowly and finally become constants. However, when the diameter of the connector is fixed, the change of DLSC with air line diameter is different from that of the DLC. With a 10 mm diameter connector, as the air line diameter increases from 10 to 100 mm, the DLSC only decreases 1.0%, while the DLC decreases 6.8%. Thus, although the load-sharing of the semi-trailer improves very slowly by increasing the size of the connector, the dynamic tire force and accordingly the road-friendliness of the semi-trailer are effectively improved. 0.052' 0.05- O M 0.048-o 0.046- a) 100 Diameter of connector (mm) Diameter of air line (mm) 0.085 0.08 ^ 0.075 Q 0.07 40^ 60 0.08 0.078 0.076 0.074 0.072 0.07 an 40 20 100 -|00 80 60 40 b) Diameter of air line (mm) Diameter of connector (mm) Fig. 4. Effect of size of air line and connector on load-sharing optimization ratios; a) DLSC, b) DLC The change in dynamic tire force can be revealed more clearly with spectral analysis. Fig. 5 shows the PSD for dynamic tire force for three types of longitudinal connections among air suspensions, i.e., type "1", type "2" (both as specified above), and type "3" (three 100 mm inside diameter connectors connecting a 100 mm inside diameter air line). It is evident that the dynamic tire forces are effectively isolated by using large air lines and connectors, especially at frequencies lower than 0.3 Hz and frequencies from 0.8 to 5.0 Hz. Compared with the semi-trailer with air connection "1", the peak values around the body bounce frequency (about 1.9 Hz) of the front dual tires decreased by 18.6 and 60.6% through employing air connections "2" and "3", and the corresponding optimization ratios of the middle, rear dual tires are 22.3, 50.5, 36.6 and 58.1%, respectively. 4 CONCLUSIONS In this study, the effects of suspension parameters on dynamic load-sharing of longitudinal-connected air suspensions of a tri-axle semi-trailer were investigated comprehensively. A novel nonlinear model of longitudinal-connected tri-axle air suspensions a) 10 10 Frequency (Hz) b) 10 10 Frequency (Hz) 10' 10 N X 10 - Q œ CL 10 - 10 10" cj Frequency (Hz) Fig. 5. Effect of size of air line and connector on optimization ratios of load-sharing; a) front axle, b) middle axle, c) rear axle was formulated based on fluid mechanics and thermodynamics and validated through test results. The effects of road surface conditions, driving speeds, air line diameters and connector diameters on the dynamic load-sharing capability of the semi-trailer were analyzed in terms of the DLSC and DLC, and the following conclusions can be drawn: 1. The road-friendliness metric DLC, is generally in accordance with the load-sharing metric DLSC. 0.05 0.048 0.046 0.044 0.065 2. When the static height or the static pressure is increased, the DLSC optimization ratio declines. The reason for this phenomenon is that the static height increases with more static pressure, Vs10 increases, and the volume of air line divided by Vs10 decreases. Thus, the effect of employing a large air line and connectors becomes less prominent. 3. Assuming that the diameter of the air line is always larger than that of the connector, the influence of the diameter of air line is more significant than that of the connector. When the semi-trailer is driving at 20 m/s on a standard "B" class road and the connector diameter is fixed at 10 mm, the DLSC only decreases 1.0% as the air line diameter increases from 10 to 100 mm. However, when the air line diameter is fixed at 100 mm, the reduction reaches 11.4% as the connector diameter increases from 10 to 30 mm; as the connector diameter continues to increase, the DLSC declines at a very slow rate and becomes constant. Based on the proposed model, investigation of the control methods of the tri-axle air suspension system with longitudinal air line and their influence on load-sharing will be undertaken in the future. 5 ACKNOWLEDGEMENT This work was supported by the National Natural Science Foundation of China (Grant Nos. 51078087, 51178158, 51075112), the Natural Science Foundation of Anhui Province (Grant No. 11040606Q39), and the Fundamental Research Funds for the Central Universities (Grant No. 2012HGQC0015 and 2011HGBZ0945). The assistance of Dr Lloyd Davis from the Queensland Department of Transport and Main Roads is also greatly acknowledged. 6 REFERENCES [1] Davis, L.E., Bunker, M.B. (2008). Load-Sharing in Heavy Vehicle Suspensions - New Metrics for Old. Queensland University of Technology, QUT Digital Repository, Brisbane. [2] Potter, T.E.C., Cebon, D., Cole, D.J., Collop, A,C. (1996). Road damage due to dynamic tyre forces, measured on a public road. International Journal of Heavy Vehicle Systems, vol. 3, no. 1-4, p. 346-62. [3] Cebon, D. (1987). Assessment of The Dynamic Wheel Forces Generated by Heavy Road Vehicles. Symposium on Heavy Vehicle Suspension Characteristics, p. 199212. [4] Cantieni, R., Krebs, W., Heywood, R. (1998). Dynamic Interaction between Vehicles and Infrastructure Experiment (DIVINE), Technical Report No. DSTI/ DOT/RTR/IR6( 98)1 /FINAL. Organisation for Economic Co-operation and Development (OECD), Paris. [5] Department of Transport and Regional Services -Australia. (2004). Vehicle Standards Bulletin VSB 11 Certification of Road-Friendly Suspension Systems. Australian Department of Infrastructure, Transport, Regional Development and Local Government, Canberra. [6] Costanzi, M., Cebon, D. (2007). An investigation of the effects of lorry suspension performance on road maintenance costs. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, vol. 221, no. 11, p. 1265-1277, D0I:10.1243/09544062JMES639. [7] Sweatman, P.F. (1983). A Study of Dynamic Wheel Forces in Axle Group Suspensions of Heavy Vehicles. Australian Road Research Board, Special Report. Report No. SR27. Vermont South, Victoria. [8] de Pont, J.J. (1997). Assessing Heavy Vehicle Suspensions for Road Wear. Research report No 95. Transfund New Zealand, Wellington. [9] Potter, T.E.C., Cebon, D., Collop, A.C., Cole, D.J. (1996). Road-damaging potential of measured dynamic tyre forces in mixed traffic. Proceedings of the Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering, vol. 210, no. 3, p. 215-225, D0I:10.1243/PIME_PR0C_1996_210_265_02. [10] Estill & Associates Pty Ltd. (2000). Operational Stability and Performance of Air Suspension on Various Vehicle Configurations. Australia: Department of Transport and Works, South Perth. [11] Roaduser Systems Pty Ltd. (2005). Stability and On-Road Performance of Multi-Combination Vehicles with Air Suspension Systems. Australia: National Road Transport Commission, Canberra. [12] Davis, L.E. (2010). Heavy Vehicle Suspensions- Testing and Analysis. Ph.D. thessis, Queensland University of Technology, Brisbane. [13] Davis, L.E., Bunker, J.M. (2011). Altering heavy vehicle air suspension dynamic forces by modifying air lines. International Journal of Heavy Vehicle Systems. vol. 18, no. 1, p. 1-17, D0I:10.1504/ IJHVS.2011.037957. [14] Roebuck, R.L., Cebon, D., Dale, S.G. (2006). Optimal control of a semi-active tri-axle lorry suspension. Vehicle System Dynamics, vol. 44, no. supl., p. 892903, DOI:10.1080/00423110600907493. [15] Papagiannakis, A.T., Zelelew, H.M, Muhunthan, B. (2007). Wavelet analysis of energy content in pavement roughness and truck dynamic axle loads. Transportation Research Record: Journal of the Transportation Research Board, vol. 2005/2007, p. 153-159, DOI:10.3141/2005-16. [16] Taskin, Y., Hacioglu, Y., Yagiz, N. (2007). The use of fuzzy-logic control to improve the ride comfort of vehicles. Strojniski vestnik - Journal of Mechanical Engineering, vol. 50, no. 10, p. 462-468. [17] Guclu, R. (2004). The fuzzy-logic control of active suspensions without suspension-gap degeneration. Strojniski vestnik - Journal of Mechanical Engineering, vol. 53, no. 4, p. 233-240. [18] Shi, X.M., Cai, C.S. (2009). Simulation of dynamic effects of vehicles on pavement using a 3D interaction model. Journal of Transportation Engineering - ASCE, vol. 135, no. 10, p. 736-744, D0I:10.1061/(ASCE) TE.1943-5436.0000045. [19] ISO 8068:1995 (E) (1995). Mechanical Vibration-Road Surface Profiles-Reporting of Measured Data. International Organization for Standardization, Geneva. [20] Chen, J.P., Chen, W. W., Zhu, H., Zhu, M.F. (2010). Modeling and simulation on stochastic road surface irregularity based on Matlab/Simulink. Transactions of the Chinese Society of Agricultural Machinery, vol. 41, no. 3, p. 11-15. [21] Prabakar, R.S., Sujatha, C., Narayanan, S. (2009). Optimal semi-active preview control response of a half car vehicle model with magnetorheological damper. Journal of Sound and Vibration, vol. 326, no. 3-5, p. 400-420, D0I:10.1016/j.jsv.2009.05.032. [22] Wang, J.S., Zhu ,S.H. (2009). Linearized model for dynamic stiffness of air spring with auxiliary chamber. Journal of Vibration and Shock, vol. 28, no. 2, p. 72-76. [23] White, F.M. (2011). Fluid Mechanics. McGraw-Hill, Columbus. [24] Chen, S.M. (2011). Hydraulic and Pneumatic Transmission. China Machine Press, Beijing. [25] Davis, L.E., Bunker, J.M. (2009). Suspension Testing of 3 Heavy Vehicles - Dynamic Wheel Force Analysis - Report. Australia: Department of Main Roads-Queensland Government, Brisbane. [26] Davis, L.E., Bunker, J.M. (2009). Heavy Vehicle Suspension Testing and Analysis-Dynamic Load Sharing. Australia: Department of Main Roads -Queensland Government, Brisbane. [27] Chen, H.X., He, Z.Y. (2008). A study on simulation of road roughness based on international roughness index. Highway, vol. 11, p. 155-160.