Multi-Objective Optimisation of Tyre and Suspension Parameters during Cornering for Different Road Roughness Proﬁles

: Effective emission control technologies and novel propulsion systems have been developed for road vehicles, decreasing exhaust particle emissions. However, work has to be done on non-exhaust trafﬁc related sources such as tyre–road interaction and tyre wear. Given that both are inevitable in road vehicles, efforts for assessing and minimising tyre wear should be considered. The amount of tyre wear is because of internal (tyre structure, manufacturing, etc.) and external (suspension conﬁguration, speed, road surface, etc.) factors. In this work, the emphasis is on the optimisation of such parameters for minimising tyre wear, but also enhancing occupant’s comfort and improving vehicle handling. In addition to the search for the optimum parameters, the optimisation is also used as a tool to identify and highlight potential trade-offs between the objectives and the various design parameters. Hence, initially, the tyre design (based on some chosen tyre parameters) is optimised with regards to the above-mentioned objectives, for a vehicle while cornering over both Class A and B road roughness proﬁles. Afterwards, an optimal solution is sought between the Pareto alternatives provided by the two road cases, in order for the tyre wear levels to be less affected under different road proﬁles. Therefore, it is required that the tyre parameters are as close possible and that they provide similar tyre wear in both road cases. Then, the identiﬁed tyre design is adopted and the optimum suspension design is sought for the two road cases for both passive and semi-active suspension types. From the results, signiﬁcant conclusions regarding how tyre wear behaves with regards to passenger comfort and vehicle handling are extracted, while the results illustrate where the optimum suspension and tyre parameters have converged trying to compromise among the above objectives under different road types and how suspension types, passive and semi-active, could compromise among all of them more optimally.


Introduction
Due to increased environmental issues, the automotive industry has increased its focus on energy efficient driving, without neglecting the aspects of comfort and vehicle stability. In this direction, effective emission control technologies and novel propulsion systems have been developed, decreasing exhaust particle emissions. However, work also has to be done on decreasing non-exhaust traffic related sources (i.e., the tyre-road interaction and the tyre wear, which both are inevitable in road vehicles).
The non-exhaust traffic-related sources have a great impact on the pavement condition (i.e., degradation and permeability) [1] and environmental pollution. Regarding the latter, in 2013, it was estimated that the tyre wear, in a couple of European countries (Germany, Netherlands, Sweden, Italy, United Kingdom, Denmark and Norway) was around 300,000 tonnes in total with about 40 % coming from passenger vehicles. At the same time, a similar amount of wear is disposed per year in the environment from vehicles in India, where the population size is 5.5 times larger [2]. According to Grigoratos et al. [3], a significant percentage of these fall in the PM 10 fraction, which means that they have a diameter larger than 10 µm and they finally end up in air, water, soils, etc. [4]. Hence, the need to develop more environmentally friendly vehicle systems that can decrease tyre wear has risen. This is especially the case now that the more environmentally friendly, but much heavier, electric vehicles are expected to have increased particle pollution from tyre wear compared to conventional vehicles, a fact that could potentially cancel the benefits of removing the exhaust emissions [5,6].
Tyre wear occurs because of friction during the sliding between the tyre tread surface and the road [7]. With the tread being the component responsible for the vehicle-road interaction, the aspect of minimising wear is a crucial criterion during the tyre design. According to Huang et al. [8], one of the existing categorisations of wear is between normal and abnormal. The first leads to uniform wear along the tyre circumference and over its width, while the latter is defined by uneven and irregular wear. Uneven wear mostly describes the non-uniform wear distribution over the tyre width, whereas the irregular wear mainly considers the circumferential wear. The total amount of wear is related to internal (tyre design, manufacturing, etc.) and external (vehicle, road, driving condition, environmental circumstances, etc.) factors ( Figure 1) [9]. According to Maitre et al. [10], the driving condition is the most dominant in terms of its impact to wear while the tyre design, the environmental circumstances, the vehicle and the road follow. This work covers the majority of these factors, as highlighted in Figure 1. Specifically, the optimisation of tyre and suspension design of a passenger vehicle is investigated to minimise wear under multiple road roughness surfaces. Various wear models have been presented in the literature, to investigate how different factors influence wear and to study the wear behaviour in detail. Bin Ma et al. [11] simplified the kinetic sliding friction coefficient taking into consideration the road roughness. Afterwards, they coupled the contact model with a 9DOF vehicle model to study tyre marks using wear quantity during the vehicle's pro and post-crash phases. Da Silva et al. [12] developed a qualitative formula for tyre wear evaluation and conducted a sensitivity analysis of a few tyre and vehicle parameters during cornering manoeuvres using a simplified single-track model. Huang et al. [8] proposed a theoretical tyre model for predicting the 3D tyre wear with regards to the roughness of the road and vehicle dynamic characteristics. While most models focus on the lateral and longitudinal direction, Sueoka et al. [13] presented a computationally efficient analytical model to study tyre wear due to vertical excitations and considered the suspension systems as well. Later, Li et al. [14] incorporated in this model a formula of tyre wear considering the temperature effect and the dynamic characteristics of the vehicle, which allowed them to analyse the effects of speed, ambient temperatures, tyre pressure and sprung mass on tyre wear. Thereby, this model is able to evaluate more accurately the wear performance while also considering the suspension effect, which is not widely discussed in the literature as far as the wear performance is concerned. This model is chosen in this work.
The estimation of tyre wear has been extensively experimentally studied in order to identify effective methods and to capture the effects of different factors on it. For example, Stalnaker et al. [15] developed a methodology for estimating tyre wear indoors, in order to establish consistent test results. Similarly, Knuth et al. [16] described a simulation method of indoor testing which could accurately capture the tyre-vehicle-driving interaction. In the same direction, Lupker et al. [17] provided a tool which could estimate numerically global tyre wear as well as qualitatively determine the wear distribution. This wear model was later further validated and investigated in a sensitivity analysis [18]. Only recently, Farroni et al. [19] developed a physical model of tyre wear to analyse the impact of thermal and frictional effects on vehicle performance. Similarly, Emami et al. [20] designed and developed a new portable test setup to study friction and wear, while Lepine et al. [21] presented a novel empirical tyre wear model for heavy vehicles that can be used to predict the wear for multi-axle vehicles based on route data and a vehicle model. In addition to the wear models and testing methods, Yamazaki et al. [22] investigated experimentally the impact of alignments such as camber angle and toe angle to the wear performance. Case studies were considered for various alignment configurations simulated real-life configurations observed in road tests. Tandy et al. [23] studied how increased shoulder wear is finally affecting the driving behavior, illustrating that tyre lateral force and overturning moment capacities increase significantly with the usage of the vehicle and hence its wear. However, in order to reduce the time-consuming experimental procedures, Tamada [24] considered the prediction of uneven tyre wear conducting progress simulation in a wearing out finite element (FE) tyre model.
Even if the modelling and estimation of tyre wear has been extensively investigated, very few works have considered the optimisation towards tyre wear minimisation. Up to now, most of these works used FEA models to assess the wear performance. For example, Koishi et al. [25] investigated the trade-off between uneven wear and wear life using multi-objective optimisation, where the objectives were evaluated using response surface methodology in order to save computational time as wear simulation was significant high, even for a super computer. Similarly, Serafinska et al. [26] suggested a multiobjective optimisation approach for uniform wear, by minimising the ratio of the contact pressure in the tyre shoulder and the contact pressure in the tyre footprint central part. However, the FEA models focus mainly on the detailed tyre modelling and they require high computational power, but they do not consider the rest of the vehicle subsystems and their interaction with the tyre. This has led to an unclear understanding regarding the tradeoff between important vehicle performance aspects and wear. For instance, the conflict with regards to comfort and vehicle stability is widely studied [27,28]; however, only Anderson et al. [29] investigated the trade-off between wear performance and handling. To the authors knowledge, there is not extensive literature on the optimisation of the vehicle and tyre parameters with regards to comfort, vehicle stability and wear performance where this work focuses. This is a considerably critical subject, as, during the conceptual design, the effort is placed upon efficient simulation and optimisation of both tyre and suspension parameters in order to decrease the development costs for physical testing.
Considering the above, in this work, the trade-off among comfort, vehicle stability and wear performance is investigated using a vehicle model combined with a wear model, which considers the vehicle dynamic characteristics, temperature effects, tyre dimensions, vehicle velocity and tyre slip angles. The emphasis is on the optimisation of both tyre and suspension parameters for minimising wear on a passenger vehicle, which is equipped with either passive or semi-active suspensions. The aim is to seek for a tyre design that is not significantly affected by different road profiles and to investigate how different control algorithms in the semi-active suspension influence the tyre wear. Hence, initially, tyre parameters (e.g., pressure, tyre width, outside radius, crown thickness and chordwise radius) are optimised (Scenario 1) for a vehicle being equipped with a passive suspension and driven over an S-Path, on which road roughness of Class A and B are assigned. Afterwards, a common optimal solution is sought among the alternatives provided by the two optimisation cases, requiring to have close design variable values and provide similar wear in both cases. Then, the identified tyre design is adopted and the optimum suspension design is sought (Scenario 2) for the two cases, but also for different suspension types. For each suspension type, a common solution with regards to their design variables is identified among the optimal alternatives provided by the two optimisation cases, requiring to have close values in their design variables. Significant conclusions regarding how tyre wear behaves with regards to passenger comfort and vehicle stability are extracted, while the results illustrate where the optimum suspension and tyre parameters have converged trying to compromise among the above objectives under different road types, and which suspension could compromise among all of them more optimally.
This paper is organised as follows. Firstly, all the models (vehicle, suspensions, tread and wear) are described and the road path and profiles used as excitations are displayed. Secondly, the validation of the models with IPG/CarMaker is illustrated. Thirdly, the formulation of the multi-objective optimisation is displayed. Then, the results are outlined. Finally, conclusions are extracted.

Vehicle Model
The vehicle model is a one-dimensional model with three degrees of freedom which only includes vertical motions ( Figure 2). The simulation model considers four basic subsystems of the vehicle: the chassis, the suspension systems, the unsprung mass and the tyres. The chassis is considered as a rigid body of mass (m s ) and is connected with the unsprung mass (m u ) through the suspension system, which is modelled as a spring and a damper (K u and C u ). In this work, the damper is considered both passive and semi-active. Regarding the tyres, a more advanced modelling is applied compared to the common quarter car model. Firstly, the tyre sidewall is modelled as a spring and a damper (K T and C T ) and connects the unsprung mass with the tyre tread (m t ), whose mass is varying dynamically according to analytical equations presented later in details. Finally, the tread element consists of linear springs and dampers (K t and C t ), which receive the unevenness of the road profile (z R ) as an excitation. Sprung mass : m szs + C u (ż s −ż u ) + K u (z s − z u ) = 0 (1) while the parameters used are selected according to IPG/CarMaker 8.0 to generate a digital twin of the software demo vehicle. The parameter values are presented in Table 1. Table 1. Vehicle model parameters.

Suspension Designs
Semi-active suspensions do not consume significant amount of energy and can provide a more energy efficient solution compared to fully active ones. At the same time, the reliability of passive suspensions can be maintained. Therefore, different control strategies have been developed exploiting the benefits of semi-active suspensions. In this work, the suspension system of the vehicle model operates as passive (PS) and semi-active. Several different control algorithms are used in the semi-active suspension case: SH-2 [30], ADD [31], SH-ADD-2 [32], PDD [33], SH-PDD [34] and GH-2 [35]. The operational conditions of the semi-active control algorithms are illustrated in Table 2. Table 2. Semi-active control algorithms applied in the vehicle suspension systems ( Figure 2).

Case
Description The parameter d ST refers to the vehicle suspension travel, i.e., z s − z u , whileḋ ST is the vehicle suspension velocity, i.e.,ż s −ż u . b The parameter α f is the crossover frequency between SH-2 and ADD ∈ [10-60] rad sec .

Tread Model
In order to evaluate the wear quantity, the energy loss has to be calculated, hence the contact between the tyre and the road should be modelled with regards to that and not only with spring and dampers as usual. Therefore, this model considers that the tread. which is in contact with the road surface (z R ), is a cuboid element with a, b and h being its length, width and height, respectively ( Figure 3a). Regarding the first two, they are evaluated according to empirical equations (Equation (4)) [9,36], which consider tyre dimensions, i.e., the outside radius (d) and the crown thickness (b 0 ), and its deformation (δ): where ψ (=115) and s (=0.67) are empirical coefficients. The tyre deformation (δ) is evaluated according to Equation (5): where α 1 is 1/(2πdr w ), with r w being the tyre chordwise curvature radius; α 2 (=3.70 * 10 −8 ) is a constant for car tyre; and F z is the wheel load. According to δ (Equation (5)), the tread deformation is related with the pressure (P) (Figure 3b), the tyre radius (d) and the chordwise radius (r w ). The force-deflection curves that occur through Equation (5) are plotted in Figure 3b for different tyre pressures (P). Finally, the mass of the tread element can be derived by Equation (6), where ρ (=1156 kg/m 3 ) is the rubber density of the tread.

Tyre Wear Quantity
Tyre wear is defined according to the frictional or wear energy, E w [Nm/m 2 ], which is the energy due to slippage per unit area and can be evaluated as the product of the shear stress ( − → τ ) and slip ( − → S ): where b is the contact length; τ x and τ y are, respectively, shear stresses in xand y-directions; and S x and S y are, respectively, slip distances in xand y-directions. Wear energy during cornering: In order to calculate the wear energy (E w ) during cornering, a solid tyre is considered using a brush model (Figure 4a), while, for the wear evaluation during cornering, a wear model is used (Figure 4b).
In the considered wear model, parabolic pressure distribution (q z ) and small slip angles are assumed. The contact pressure is defined as: where p m is the maximum pressure that occurs from the vertical load F z and is equal to 3F z /2ab if the contact patch is assumed as cuboid (Figure 3a). According to the wear model (Figure 4b), the lateral slip angle is α , while the lateral shear stress at a distance x measured from the the contact patch (0-b) is τ y and is expressed by Equation (9): where C y [N/m] is the lateral shear elastic tread coefficient per unit area. The shear force τ z (=µ s q z b h ) applied to the contact patch per unit area is equal to the lateral shear stress applied to the sliding point b h (Equation (9)), i.e., the maximum frictional force per unit area (τ y max ). Therefore, combining Equations (8) and (9): where µ s is the static friction coefficient, while ζ = C y btanα 4µ s p m is introduced as a new coefficient, to help in the simplification of the equations. Based on ζ, the sliding point (b h ) is defined as: where the first part is when the tyre is in the adhesion region, while the second part considers that sliding occurs across the whole contact area when ζ h is greater than or equal to unity. Then, considering the force equilibrium in the sliding region (b h ≤ x ≤ b), the sliding distance is defined according to Equation (12), where µ d is the kinematic friction coefficient: Considering Equation (7), the wear energy per unit area is defined as: where Equations (8) and (11) can be substituted, and the equation of wear energy can be derived at the two conditions explained in Equation (11). In addition, in this work, µ d and µ s are assumed to be equal to µ, while F z is considered from the vehicle model: where F Z 0 is the static load applied. Mass Loss: The wear energy per contact area described above can be related with the mass loss according to Equation (15), which was derived by Lupker et al. [17] from experiments: where f 1 and f 2 are constants at a given temperature. The mass loss (∆m) evaluated at each time step is used to calculate the m t dynamically. As an extension of this expression, Li et al. [14] redefined f 1 by incorporating in it the effect of temperature: where f 2 is 1.5106 and f 0 (=2.010 −10 ) is the constant f 1 at T 0 (=60 • ). T t is defined as the steady state temperature of the tyre tread: where V[m/s] is the vehicle velocity; r is rim radius; c is the tyre width; γ (=0.12) is the tyre hysteresis coefficient; β (=1.40) is the correction coefficient; and T ∞ is the ambient temperature.
Wear quantity: The total wear quantity that is caused by the mass loss (Equation (16)) is defined as the height change (∆h): Finally, the total tread wear is defined: where H(t − t 0 ) is the total tread wear occurred during the first tyre revolution taking place as long as t ≤ t 0 , with t 0 (=2πr r /V) being the rotating period and r r the rolling radius.

Ride Comfort
Ride comfort (RC) is assessed via the root mean square of the sprung mass' vertical acceleration (z s ), when measurements are not applicable to the seat or the occupant's body [37]. More specifically, RC is calculated as follows: where T ([s]) is the duration andz w z (m/s 2 ) is the weighted acceleration as a function of time. In this work, no weighting is applied in the accelerations measured at the sprung mass. Therefore,z w z is considered equal toz s .

Vehicle Handling
The vehicle handling is normally assessed through cornering, braking and traction abilities, which all are enhanced if the tyre loads/deflections are maintained at a low level. Therefore, the root mean square of the tyre deflection (z u − z Road ) is considered as an index (Equation (21)) for vehicle stability (Figure 2), taking advantage of the available outputs from the vehicle model: At the same time, the suspension travel depicts the ability of the suspension system to support and hold the vehicle load. This is achieved when rattle space requirements of the vehicle are kept small. Therefore, the root mean square of the suspension travel (z s − z u ) is considered as an index (Equation (22)) for vehicle stability as well:

Road Profiles and Path
The combination of the vehicle model with the above analysis on the wear quantity allows the assessment of the wear under different conditions (i.e., road surfaces, paths and vehicle velocities), tyres and vehicle types. More specifically, different road surfaces can be used for exciting the vehicle model and different paths can be considered by applying in the wear model the vehicle's corresponding slip angles and velocities while driving on any path. The vehicle's corresponding slip angles and velocities have to be obtained from offline simulations prior to the simulation of the vehicle model.
In this work, road profiles of Class A and B (Figure 5a) are used to excite the vehicle model and are designed according to ISO-8608 [38]. Then, the S-path illustrated in Figure 5b, which consists of two turns of different radius (70 and 60 m), is used to extract the slip angles and the vehicle velocity while driving on it. In this work, both are obtained from simulating the demo vehicle of IPG/CarMaker 8.0, IPG AUTOMOTIVE, Braunschweig, Germany ( Figure 6) while driving over the S-path, having assigned on it the two road profiles.

Validation of the Model with IPG CarMaker 8.0
In order to secure that the optimisation results in this work are realistic, the model is validated using IPG/CarMaker 8.0 for the path illustrated in Figure 5b and the road profile of Class B assigned to it (Figure 5a).
As mentioned above, the vehicle parameters selected for the model are extracted from the software demo vehicle to build its digital twin. The suspension and tyre parameters of the models are adjusted according to Tables 1 and 3. More specifically, the suspension parameters (Table 1) of the front wheels are adjusted into the section Car/VehicleDataSet/Suspensions of the software. Regarding the suspensions, for both the stiffness and the damping coefficient, a characteristic value option is selected, while the values of the buffer are set zero in the software, as the buffer is not included in the model. Then, regarding the tyres, the parameters related with them and considered in the vehicle and wear model (P, d, c, r, K T and C T ), are adjusted according to section Tires/TireDataSet/GeneralandModelParameters of the software. The only parameters that are not provided are the crown thickness (b 0 ), the chordwise curvature (r w ) and the lateral shear elastic tread coefficient (C y ). Regarding the first two, its consideration is not possible, but C y is extracted after evaluating it using the cornering stiffness obtained from the software simulation results. More specifically, the side force applied in the wheel for a rectangular contact area, as in this case, can be evaluated for small slip angles: Thus, according to Equation (23), the cornering stiffness C Fa is equal to C y ab 2 2 α. Therefore, after fitting the simulation results (lateral slip angles and side forces) to a linear equation, as shown in Figure 6b, the value of C y is extracted. Having created the digital twin of the demo vehicle ( Figure 6), the wear energy (E y [Nm/m 2 ]) is selected for comparison in order to validate the simulation results. Therefore, the obtained measurements are used as follows for calculating the wear energy (E w y carmaker ): • The sliding distance (s d ) is calculated according to Equation (12). • The parameter (φ) is calculated as follows: • The sliding force is obtained by integrating the maximum possible frictional force distribution (µq z ) over the sliding region b h ≤ x < b: • Finally, the wear energy per contact area according to Equation (26) is evaluated by multiplying the sliding force with the slipping distance (S y ) and dividing it by the contact area: The comparison of E w y carmaker with E w y is illustrated in Figure 7, which proves a good convergence in the two simulation results. More specifically, Figure 7 illustrates high similarity in terms of the trend of the two curves. Slight differences occur in the magnitude with the model used in this work underestimating a few of the peaks. However, this underestimation is not an issue as the trend is accurately followed, and, given the simplicity of the model in comparison with IPG/CarMaker 8.0, the results should be considered valid.

Optimisation Configuration
In this work, the emphasis is on the optimisation of tyre and suspension parameters for minimising tyre wear. However, as comfort and vehicle stability are still important criteria in the design process, they are also included as objectives by configuring a multiobjective optimisation problem. In order to decrease the complexity of the optimisation, the tyre and suspension parameters are optimised separately considering two optimisation scenarios. Initially, the optimal tyre design solutions (Scenario 1) are sought for the vehicle equipped with a passive suspension, while driving over a road class A (Case 1) and B (Case 2) S-path ( Figure 5). Then, after selecting one of the optimal tyre design solutions, which can provide similar tyre wear for both road profiles, the optimal suspension designs (Scenario 2) are sought for different suspension types, both passive and semi-active (PS, SH-2, ADD, SH-ADD-2, PDD, SH-PDD and GH-2), for the two road cases.

Objectives
The multi-objective optimisation problem is formulated in order to minimise tyre wear, enhance ride comfort and increase vehicle handling at the same time. In both scenarios, wear and comfort are represented from the same objective function (F 1 and F 2 ), whereas vehicle handling (F 3 ) is represented with the most suitable function according to the scenario. In Scenario 1, where the tyre is optimised, TD, defined in Equation (21), is selected as the objective function, while, in Scenario 2, where the focus is on the suspension, ST is selected, as defined in Equation (22).
In both scenarios, the optimisation is conducted using the GAMULTIOBJ toolbox from MATLAB 2017b, where the PopoulationSize is selected as 250 and the EliteCount as 5. The rest of the options are selected according to the default suggestions from the software.

Design Variables: Scenario 1
In this scenario, the following design variables are selected: P, c, d, b 0 andr w (29) while the upper and lower bounds set are illustrated in Table 4. In addition, K t and K T might be considered indirect design variables as they are evaluated according to the problem design variables. More specifically, during the optimisation, the force-deflection curve (Figure 3b) is evaluated according to Equation (5), which is a relation of P, d and r w , for F z varying from 500 to 7000 N for each potential optimal solution. Then, having obtained the force-deflection curve, the K t is evaluated as the curve slope at the linear region, while the sidewall stiffness (K T ) is evaluated by 0.8 * K t as the sidewall is always less stiff than the tread.

Design Variables: Scenario 2
Regarding the second scenario, where the optimum designs are sought for the different suspension types, the design variables shown in Equation (??) are selected: For PS: For SH-2, ADD, GH-2, PDD and SH-PDD: For SH-ADD2: while the upper and lower bounds set in the optimisation are illustrated in Table 4.

Results
As described above, in this work, two optimisation cases, with regards to different excitations, are studied in two optimisation scenarios where the focus is on tyre and suspension design, respectively. In this section, the results of the two multi-objective optimisation scenarios are displayed. More specifically: • Scenario 1 investigates the tyre optimisation in a vehicle being employed with a passive suspension and driven over an S-Path, in which road surfaces of Class A (Case 1) and B (Case 2) are assigned. The optimisation results are illustrated in Figures 8 and 9, which display the Pareto fronts with regards to the objectives (Figure 8a,b) and the design variables (Figure 9a-e). Having obtained the optimal solution alternatives for the two cases, common solutions are sought between them, requiring to provide similar wear in the two cases (<8% difference) and have close design variable values (<9% difference). The common solutions identified with the above characteristics are illustrated in Figures 8 and 9 alongside the alternatives, while their values are listed in Table 5. In addition, in Figure 10, the final circumferential tread profiles for the optimum designs are compared with the initial tyre design, after driving over the Class B S-Path twenty times back and forth. • Scenario 2 explores the suspension optimisation for two road cases and various suspension types (PS, SH-2, ADD, SH-ADD-2, GH-2, PDD and SH-PDD). The optimisation results are illustrated in Figure 11, which presents the Pareto fronts of the optimum solutions for the two cases. For each suspension type, a common solution is identified among the optimal alternatives provided by the two optimisation cases, requiring them to have close design variables values (<3-15% difference according to the cases). Finally, the values of the identified common solutions are presented in Table 6 alongside with the threshold of difference that was allowed.
The goal of these scenarios is for conclusions to be extracted regarding: (a) how tyre wear behaves with regards to passenger comfort and vehicle handling; (b) where optimum tyre and suspension parameters have converged trying to compromise among the above objectives under different road types; and (c) which suspension could compromise among all of them optimally.

Optimum Tyre Design Solutions Optimisation Results
Regarding Scenario 1, Figure 8a,b illustrates the relation between the optimisation objectives and reveals the trade-offs among them. More specifically, Figure 8a reveals the conflicting relation between wear and comfort levels, where the increase of the RC index, which yields the increase of discomfort, leads to the decrease of tyre wear. At the same time, the conflicting relation of comfort and handling, captured by the optimisation in Figure 8b, serves as an indication of the validity of the results. In addition, the conflicting relation of the two objectives, handling and wear, with comfort reveal the linear relation between them, implying that increased levels of tyre deflections lead to increasing wear.
As far as the design variables are concerned, according to Figure 9a-e, the tyre pressure (P), the outside radius (d) and the chordwise curve (r w ) illustrate a conflicting relation with wear. For these, the algorithm has provided optimal alternatives ranging from the lower bound to the upper bound, which allowed capturing the conflicting relation between these variables and wear. On the other hand, the Pareto alternatives of tyre width (c) and crown thickness (b 0 ), as shown in Figure 9b,d, respectively, have not captured their relation with wear, as the algorithm converged close to one value. In the case of the tyre width, this value is different from the upper or lower bound, whereas the crown thickness alternatives have converged close to the lower bound trying to decrease the length of the tread (Figure 3a). However, both values are realistic according to the literature.
Regarding the different road cases, Case 2 alternatives are providing larger values in all the objectives, as the more intense road roughness profile has increased discomfort, vehicle instability and tyre wear. Regarding the latter, according to Figure 8a, Case 2 Pareto alternatives have increased wear around ∼25% compared to Case 1, which illustrates how road profiles increase wear significantly if different tyre designs are adopted. In terms of convergence, Case 2 has provided more scattered fronts in both the objective and the design variables fronts, illustrating that the less intense excitation (Case 1 with Class A road roughness) allowed the algorithm to converge to a more solid Pareto front. However, in both cases, the alternatives mostly have captured the same conflicting relations, while the design variables are within the same range of values.
Finally, as mentioned above, common solutions are sought requiring for them to provide similar wear in the two cases (<8% difference) and have close design variables values (<9% difference). This search led to two common solutions between the optimal alternatives of the two cases, as pointed out in Figures 8 and 9 and Table 5. Based on these, in Case 1 where road Class A is considered, the common solutions are closer to the upper edge of the Pareto front, being more comfort oriented in this case, while in Case 2 they are closer to the lower edge, being more stability oriented. Regarding their design variable values (Table 4), they are all the same with only the pressure (P t ) being significantly different. This outcome of the optimisation implies that the optimum tyre design has two pressures (P) with which it can drive optimally over both a road Class A and B S-path with the same levels of wear. The similarity of these two solutions in terms of the overall wear is illustrated in Figure 10, where the circumferential tread profile that occurs with these solutions is compared along with the one that would occur with the initial tyre design, as shown in Table 3. In addition, Figure 10 displays how the optimisation has improved the tread circumferential profile by not only minimising the overall wear, but also smoothening the variations compared to the initial solution. (a) (b) Figure 10. The tread profile induced to the tyre for the two common optimal tyre designs, after the vehicle is driven over the Class B S-path 20 times back and forth.

Optimum Suspension Design
Regarding Scenario 2 and the suspension optimisation, the Pareto fronts of the different suspension types with regards to the three optimisation objectives are illustrated in Figure 11, while the common solutions identified are pointed out in Figure 11a.  Based on the results in Figure 11b,e, the relation between comfort and handling is consistent with what was expected, illustrating the conflict between the two objectives. At the same time, the optimisation algorithm has converged to solutions with similar wear characteristics for most of the suspension types as the wear objective has small differences (Figure 11a,c,d), especially for Case 1. On the one hand, the above outcomes imply that for a suspension system, after having decided the appropriate suspension type and the optimal tyre design with regards to wear, the engineers can focus on the rest of the objectives to tune the suspension optimally. On the other hand, different control algorithms could make a significant difference in terms of the tyre wear performance. According to the Pareto fronts, the increase of the wear from Case 1 (Figure 11a, Class A) to Case 2 (Figure 11d, Class B) is ∼13% for most of the semi-active suspensions, illustrating their lack of adaptability in different road profiles. However, SH-2 and SH-ADD-2 (Figure 11c) provide the minimum levels of wear, with insignificant differences between them, and maintained the same tyre wear levels for the two cases.
Regarding the compromise between comfort and vehicle handling, the Pareto fronts of all the suspension types in both cases have similar trends, and they validate the theoretical basis of each control algorithm. More specifically, SH-ADD2 is trying to balance the characteristics of SH-2 and ADD, with its Pareto alternatives being between the SH-2 and ADD Pareto fronts. In this direction, the SH-ADD2 Pareto front connects the upper part of the SH-2 front (the stability oriented) and the lower part of the ADD front (the comfort oriented). Similarly, SH-PDD is trying to balance the characteristics of SH-2 and PDD with its Pareto front alternatives being between those of SH-2 and PDD, while it also connects the upper part of the PDD front (the comfort oriented) and the lower part of the SH-2 front (the stability oriented). At the same time, GH-2 is providing solutions with more vehicle stability but less comfort. Finally, the PS Pareto alternatives in both cases include solutions that could replace few semi-active algorithms (SH-ADD2, PDD and ADD) but with slightly either less comfort or less vehicle handling.
As far as the common solutions are concerned, in most of the suspension types, the common identified solution between the two road cases is one that enhances comfort more, being in the top part of the Pareto font (Figure 11a). Only in PDD, the common solution is located at the bottom part of the Pareto front, as shown in Figure 11a, which is the one that emphasises vehicle handling. The design variables are are illustrated in Table 6. Based on Table 6, for a few of the suspension types (PS, GH-2, PDD and SH-ADD2), the design variables have converged to a soft spring, while the rest (SH-2, ADD and SH-PDD) to a much stiffer one. Regarding, SH-2 and SH-ADD-2, which are the ones that provide similar tyre wear levels for both road profiles, they have converged to different suspension tuning and also the common solution was identified with a high allowed threshold, which could be the reason these two algorithms provide minimum tyre wear for the two road cases. The first has converged to a stiff spring and damper coefficient for both C min and C max , however the comfort levels are maintained to good levels for both road cases. On the other hand, SH-ADD-2 has converged to a soft suspension design, with spring (K 3 ) and damping coefficients (C min ) close to the lower bounds allowed by the optimisation. However, the vehicle handling is maintained to similar levels with the SH-2 solution.

Conclusions
To sum up, this work considers a vehicle model combined with a tyre wear model in order to optimise tyre and suspension parameters to minimise tyre wear, enhance comfort and improve stability. In addition to the search for the optimum tyre and suspension design, the optimisation is also used as the tool to identify any conflicting or not relations between the design variables and the optimisation objectives.
In summary, the following conclusions are extracted.
• Comfort illustrates the same conflicting relation with wear as with vehicle stability. This means that the increase of the suspension travel, hence degradation of handling, leads to an increase in wear but to an improvement in comfort. • Wear could increase up to 21% in different road profiles, while the appropriate tyre design could provide only 2% increase if the road roughness changes from Class A to B. The appropriate tyre design could be extracted using the optimisation method proposed in the current work. • For the same tyre design, two pressures could optimally combine comfort, wear and stability in two road roughness cases but also maintain wear at the same levels in these two cases. The recursive feasibility of this outcome should be tested, but it is an interesting remark for the current case study. • Regarding suspension types, according to the results, the type and their control algorithm should be selected with regards to tyre wear damage. In the current case study, SH-2 and SH-ADD-2 seemed to be able to provide the best wear performance. After the selection of the control algorithm, the tuning of the suspension parameters, i.e., stiffness and damping coefficient, should take place mainly with regards to comfort and vehicle handling. This is because according to the results it seems that different configurations do not affect wear to a great extent.
Further work is in progress to develop active suspension systems using novel control algorithms to decrease wear. At the same time, tyre and suspension design should be investigated in depth with regards to the tyre wear using more detailed wear and vehicle models. Finally, their experimental validation and investigation is an aspect that should be considered further.