Next Article in Journal
Strigolactone (GR24) Application Positively Regulates Photosynthetic Attributes, Stress-Related Metabolites and Antioxidant Enzymatic Activities of Ornamental Sunflower (Helianthus annuus cv. Vincent’s Choice) under Salinity Stress
Next Article in Special Issue
Enhancement of Agricultural Materials Separation Efficiency Using a Multi-Purpose Screw Conveyor-Separator
Previous Article in Journal
Data Classification and Demand Prediction Methods Based on Semi-Supervised Agricultural Machinery Spare Parts Data
Previous Article in Special Issue
The Improvement of a Traction Model for Agricultural Tire–Soil Interaction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Vibration Performance Analysis and Multi-Objective Optimization Design of a Tractor Scissor Seat Suspension System

1
College of Vehicle and Traffic Engineering, Henan University of Science and Technology, Luoyang 471003, China
2
School of Vehicle and Traffic Engineering, Henan Institute of Technology, Xinxiang 453000, China
3
Advanced Manufacturing of Mechanical Equipment Henan Province Collaborative Innovation Center, Henan University of Science and Technology, Luoyang 471003, China
4
General R&D Institute of China FAW Group Co., Ltd., Changchun 130022, China
5
State Key Laboratory of Automotive Simulation and Control, Jilin University, Changchun 130022, China
*
Author to whom correspondence should be addressed.
Agriculture 2023, 13(1), 48; https://doi.org/10.3390/agriculture13010048
Submission received: 18 November 2022 / Revised: 19 December 2022 / Accepted: 20 December 2022 / Published: 23 December 2022

Abstract

:
The combination of characteristic parameters is the key and difficult point to improving the vibration attenuation of scissor seat suspension. This paper proposes a multi-objective optimization method based on entropy weight gray correlation to optimize the combination of characteristic parameters with better vibration attenuation. The differential equation of seat suspension motion is derived through mechanical analysis, and a simplified driver seat suspension single degree of freedom model is constructed. The range of spring stiffness and damper damping is calculated theoretically. Through main effect analysis and analysis of contribution, the main influencing factors of seat suspension vibration attenuation are studied, and the influence correlation of the main factors is analyzed. On this basis, the spring stiffness and damper damping are taken as control variables, and the upper plane acceleration, displacement, and transfer rate of the seat suspension are taken as optimization objectives. The Optimal Latin Hypercube Sampling (OLHS) was used to sample the Design of Experiments (DoE), fit the RBF surrogate model, and screen the optimal solution based on the MNSGA-II algorithm and entropy weight gray relation ranking method. The comparative analysis of the performance before and after optimization shows that the vibration reduction performance response indexes of the acceleration, displacement, and transmissibility of the optimized seats are increased by 66.41%, 2.31%, and 8.19%, respectively. The design and optimization method proposed in this study has a significant effect on the vibration reduction of seats, which provides a reference for the optimization of the vibration reduction performance of seat suspension.

1. Introduction

Agricultural mechanization significantly impacts agricultural development, and its level is the main research focus on the agricultural machinery industry [1,2]. Tractors are an important piece of agricultural machinery, and they can be connected to various machines to perform a wide range of agricultural tasks [3]. Due to the strong vibration caused by the complex working conditions of the tractor [4], the driver’s long-term work in this environment will impact his physical and mental health [5]. Tractors are bound to produce certain vibrations during farming operations, and the main reasons for these vibrations are the vibration of the whole vehicle caused by uneven road surfaces [6]; vibrations caused by sudden acceleration, braking, and steering of the vehicle while driving on the road [7]; the vibration caused by coupling when the engine and transmission mechanism is rotating [8]; among them, the uneven road surface is the main cause of vehicle vibration. Therefore, it is meaningful to study the vertical vibration caused by uneven road surfaces [9]. In order to improve tractor ride comfort and driving safety, extensive research has been conducted to reduce tractor vibration from tractor tires, chassis suspension, engine mounting, and seats [10,11]. Improving tire and suspension parameters improves ride comfort while affecting tractor handling stability, braking performance, load-bearing capacity [12], etc. Comparatively speaking, because the seat system has the advantages of simple structure, convenient manufacture, and low cost, using seat suspension to improve tractor rides comfort is the most direct and effective method [13,14]. Tractors face complex operating conditions, and it is particularly important to improve seat suspension vibration attenuation [15].
In terms of seat suspension models, Desai R. et al. [12] validated that seat suspension with spring parallel and diagonal damper arrangement has better ride comfort than other arrangements, which provides guidance for seat suspension model design. Liao X. et al. [16] presented a mathematical model of a seat suspension system based on a negative stiffness structure, discussed the influence of different parameters on the dynamic stiffness of seat suspension, and obtained the ideal configuration parameters range of the suspension system. Du X.-M. et al. [17] introduced a simplified semi-active scissors seat suspension model to reduce the high amplitude vibration transmitted to the driver and proposed a robust state feedback H control based on an MR damper, which significantly reduced the vibration transmitted to the driver. For scissors-type seats with good stability and high reliability [18], they are widely used in agricultural tractors and commercial vehicles [19]. Therefore, the scissors-type seat is used as the research object in this study. In terms of seat suspension damping and stiffness control, Deng L. et al. [20] studied the influence of variable stiffness and damping on seat vibration attenuation, designed non-linear stiffness control and skyhook damping control, which significantly improved the vibration attenuation of seat suspension under different vibration excitations. Zhao Y. et al. [21] proposed a new fast system parameter identification method based on vibration test data for the seat system and estimated seat parameters such as stiffness and damping of the seat suspension system. Maciejewski I. et al. [22], constrained by seat suspension travel, designed an active finite-time sliding mode controller and implemented a multi-objective control strategy, which improved driver comfort with limited active power. He S et al. [23] took nonlinear damping coefficient and suspension stiffness as design variables and frequency-weighted RMS of driver’s seat acceleration as the objective function and proposed an optimization method of ride comfort based on nonlinear damping and an intelligent algorithm, which effectively improves the practicability of seat model. In the tractor seat suspension, the spring stiffness and damper damping play a key role in affecting the vibration attenuation of the seat suspension [20,21,22,23,24]. Research on the combination matching of seat suspension spring stiffness and damper damping has become an outstanding problem to be solved urgently in order to improve seat vibration attenuation [25]. From the above study, stiffness and damping matching are key factors in improving the vibration attenuation of seat suspension.
Tractor seat suspension damping optimization covers both linear and non-linear factors, which makes the optimization of the seat suspension damping process consume a lot of human and material resources [26]. In order to improve the efficiency of optimization, researchers adopt the method of combining the agent model with an optimization algorithm to carry out a multi-objective optimization design [27]. The approximate model uses an efficient mathematical model or empirical formula to replace the actual analysis model without reducing the accuracy, thus greatly reducing the calculation cost of simulation cost and significantly improving the efficiency of optimization designs. At present, the commonly used approximate model methods include Response Surface Methodology (RSM), Radial Basis Function Neural Network (RBFNN), and Kriging [28]. Jiang R. et al. [29] took quality and fatigue life as optimization objectives, carried out a multi-objective optimization design for the control arm based on the Kriging proxy model and NSGA-II algorithm, determined the optimal design of the control arm from Pareto solution by Entropy Weight Gray Relation Analysis (GRA), and proves the reliability and validity of the proposed multi-objective optimization method. Pandey M. et al. [30] used Radial Basis Function (RBF) integrated proxy modeling and NSGA-II algorithm to model the dynamic performance function of a truck with a three-piece bogie. The results show that the proposed optimization solution is significantly superior to the existing one.
Tractor seat vibration attenuation involves multi-disciplinary optimization designs of seat suspension structural dynamics, combination matching of characteristic parameters, structural safety, etc. Due to the interdependence and coupling among various disciplines, it is difficult to balance each performance using the traditional single-objective optimization method and to highlight the comprehensive vibration attenuation of seat suspension. Traditional multi-objective optimization methods mainly include a weighted average method, Benson method, interactive rule method, etc. The multi-objective optimization problem is transformed into a single-objective optimization problem through the weights between the objectives for solving, which is difficult to reflect the actual situation of the problem, and easy to cause local convergence and miss the optimal global solution [31]. Scholars have developed modern intelligent optimization algorithms based on population algorithms, mainly including improved Nondominated Sorting Genetic Algorithm (NSGA-II), Genetic Algorithm (GA), Particle Swarm Optimization (PSO) [32], etc. Modern intelligent optimization algorithm obtains Pareto solution by establishing a reasonable multi-dimensional mapping relationship between design space and target space, which more comprehensively reflects the actual situation of the problem, and improves the convergence accuracy and efficiency [33]. Prasad V. et al. [34] took cargo damage as a design criterion, proposed a GA-based fast convergence optimization algorithm, and compared the proposed algorithm with NSGA-II. This algorithm reduces calculation costs and performs better in cargo safety. Hua Y. et al. [35] aimed to stabilize suspension capability, used a multi-objective optimization method based on the RSM model, and proposed an improved Multiple Objective Particle Swarm Optimization (MOPSO) algorithm is proposed. The optimized design scheme has better vibration attenuation capability.
When considering vehicle use and research costs, the researchers focused on seat suspension in passenger and commercial vehicles [36,37]. Most tractor research enterprises neglect the consideration of driver’s ride comfort, which lays a hidden danger to the driver’s physical and mental health and driving safety [38]. Compared with heavy-duty engineering vehicles, tractors are required to meet more driving conditions due to the requirements of work tasks and carrying capacity, while the driving conditions of engineering vehicles are worse [39]. Therefore, the seat suspension of heavy-duty engineering vehicles needs to be designed for specific operating conditions. Relatively speaking, the design of tractor-seat suspension should meet the requirements of more driving conditions and lower costs. At present, the research on tractor-seat suspension mainly focuses on active and semi-active seat suspension. Vibration is attenuated by controlling stiffness and damping of seat suspension with a magnetorheological damper [40,41]. The purpose of this paper is to optimize the seat suspension system to improve the vibration attenuation according to the spring stiffness of the seat suspension and the damper damping range so that the tractor can have good vibration attenuation under most driving conditions. With the improvement of agricultural mechanization technology, people have higher requirements for the driving comfort of agricultural machinery. Therefore, as an important part of agricultural machinery, the demand for improving the comfort of tractor seats is increasingly significant.
To sum up, the vibration attenuation design of the seat suspension system mainly focuses on stiffness and damping control, and there is little research and application on tractors. Since the research mainly focuses on the active and semi-active control of the stiffness and damping of the seat suspension system, the high cost limits its application in agricultural machinery. The development of seat suspension based on the optimization of stiffness, damping, and structure design is more suitable for the application of tractors. Therefore, this paper combines the approximate model method and the improved multi-objective optimization algorithm to optimize the tractor scissors seat suspension system.
In this paper, a scissors-type seat suspension is taken as a model, and a seat suspension mechanical model and a three-dimensional Adams model are established to analyze its dynamic performance. To study the influence law of seat suspension vibration performance and identify the main influencing factors of seat suspension vibration reduction. The multi-objective optimization method is used to optimize the vibration attenuation of seat suspension. The Pareto solution is sorted by the GRA method of Entropy Weight, and the optimal solution is selected. The vibration attenuation of the optimized before and after scissors seat suspension is compared and analyzed.

2. Methods

2.1. Main Effect Analysis

The main effect of a factor on the response is the average of the response of a factor to all tests at a certain level, so changing the level of a single factor, using the average of the effects of all possible combinations of each level and other factors on the results, provides the main effect. Typically, experimental designs calculate the main effect of a factor of response by constructing a multiple quadratic regression model based on the results of input factor and output response samples.
y ¯ ( x ) = y ¯ ( x 1 , x 2 , , x m ) = β 0 + i = 1 m β i x i + i = 1 m β i i x i 2 + i = 1 m 1 i < j = 2 m β i j x i x j
As an example, the polynomial composition of the model with two input variables is as follows:
y = c 0 + c 1 x 1 + c 2 x 2 + c 3 x 1 2 + c 4 x 2 2 + c 5 x 1 x 2
Derivation of the above second order polynomial yields:
d y = c 1 d x 1 + c 2 d x 2 + 2 c 3 x 1 d x 1 + 2 c 4 x 2 d x 2 + c 5 d x 1 x 2
Then, x 1 , x 2 the Linear Main Effect of the linear term is:
M x 1 = c 1 d x 1 , M x 2 = c 2 d x 2
x 1 , x 2 the Linear Main Effect of the second order term is:
M x 1 2 = 2 c 3 x 1 d x 1 , M x 2 2 = 2 c 4 x 2 d x 2
Among them, M x 1 2 = 2 c 3 x 1 d x 1 , M x 2 2 = 2 c 4 x 2 d x 2 .
x 1 , x 2 the interaction effect is.
M x 1 x 2 = c 5 d x 1 x 2 = c 5 ( x 1 d x 2 + x 2 d x 1 )
Among them,
d ( x 1 x 2 ) = [ M a x ( x 1 ) M i n ( x 2 ) + M i n ( x 1 ) M a x ( x 2 ) ] [ M a x ( x 1 ) M a x ( x 2 ) + M i n ( x 1 ) M i n ( x 2 ) ]
The main effect has the following laws: (1) numerical value (absolute value): the larger the value of the main effect, the greater the degree of influence of the factor on the response. (2) Direction: a positive main effect indicates that the response increases with the increase of the factor; conversely, a negative main effect indicates that the response becomes smaller with the increase of the factor. (3) Order: the presence of a second-order main effect means that the effect of the factor of the response is not linear. (4) Interaction effect: a large interaction effect means that two parameters varying simultaneously will greatly affect the response.

2.2. Analysis of Contribution

Contribution analysis mainly uses the Regression of DOE to calculate the contribution. According to the ranking of the contribution of design variables to performance objectives, design variables in the high discrete or high non-linear analysis are screened to reduce calculation costs and improve the efficiency of optimization design. The analysis and calculation steps are as follows.
Step 1: Normalization process
The DOE method is used to obtain experimental samples between design variables and performance objectives. The design variables have different design spaces, and the contribution values vary from the design space, requiring normalization of the sample data inputs using Equation (7).
x i = x i x ¯ σ = ( x i 1 N i = 1 N x i ) × [ 1 N i = 1 N ( x i x ¯ ) 2 ] 1
where: x ¯ is the mean of the sample data, σ is the standard deviation, N is the total number of sample data, x i denotes the original input, and x i denotes the normalized input.
Step 2: Contribution Analysis
If there are k design variables ( x 1 , x 2 , , x k ), then any performance objective can be formulated in a multiple regression model as:
P ( x 1 , x 2 , , x k ) = μ + i = 1 k Q i ( x i ) + + i = 2 k j = 1 k 1 R i j ( x i , x j ) + ε
where: P ( x 1 , x 2 , , x k ) is any performance target, i = 1 k Q i ( x i ) is the main effect of the design variable, i = 2 k j = 1 k 1 R i j ( x i , x j ) is the cross effect of any two design variables, μ is the total mean, μ is a constant term, and ε is the error.
The main effects of the design variables can be expressed in Equation (9).
i = 1 k Q i ( x i ) = i = 1 k β ^ i x i
Therefore, the contribution values of the design variables can be defined by Equation (10).
N x i = 100 β ^ i i | β ^ i | i = 1 , 2 , , k
where: β ^ is the main effect coefficient x i of the design variable calculated by the least squares method; N x i is the contribution of the design variable x i .

2.3. Gray Relation Analysis

Gray Relation Analysis (GRA) is a method to measure the degree of approximation between experimental sequences and ideal sequences using Gray Relation Grade (GRG), which is widely used in multi-objective and multi-decision optimization problems and has comprehensive advantages in solving complex decision problems. The specific steps of GRA are as follows.
Step 1: Data pre-processing
Due to the different orders for magnitude of the experimental data, GRA may not be able to obtain a reliable optimized solution, and the experimental data need to be transformed into a set of dimensionless data x onto (0 < x < 1) for further quantitative analysis. Depending on the characteristics of the response characteristics, different data preprocessing techniques can be used.
If the target has the characteristic of “bigger is better”. The normalization method can be expressed as follows:
x i * ( k ) = x i ( k ) min k x i ( k ) max k x i ( k ) min k x i ( k )
If the objective is of a “lower is better” nature. The normalization method can be expressed as follows:
x i * ( k ) = max k x i ( k ) x i ( k ) max k x i ( k ) min k x i ( k )
If a specific value is ideal for the target. The normalization method can be expressed as:
x i * ( k ) = 1 | x i ( k ) T | max [ max k x i ( k ) T , T min k x i ( k ) ]
where: x i * ( k ) is the kth response characteristic value of the ith experiment after normalization; x i ( k ) is the initial design value of the response characteristic; min k x i ( k ) and max k x i ( k ) are the minimum and maximum values of all response characteristics x i ( k ) , k = 1, 2, … n, i = 1, 2, ... m; m is the number of experiments; n is the number of response characteristics; and T is the specific value.
Step 2: Calculate Gray Relative Correlates (GRC)
The gray correlation coefficients are obtained. The gray correlation coefficient of the kth response characteristic of the ith experiment is expressed as:
γ ( x 0 * ( k ) , x i * ( k ) ) = min + ζ max 0 i ( k ) + ζ max
where: x 0 * ( k ) is the reference experimental sequence; x i * ( k ) is the initial experimental sequence; 0 i ( k ) = | x 0 * ( k ) x i * ( k ) | is the absolute difference between x 0 * ( k ) and x i * ( k ) ; max = max i max k 0 i ( k ) and min = min i min k 0 i ( k ) are the maximum and minimum values of 0 i ( k ) , respectively; ζ is the identification coefficient, ζ [0, 1], generally defined as 0.5.
Step 3: Gray Correlation Gauge (GRG)
The gray correlation coefficient (GRC) was averaged to calculate the gray correlation degree (GRG), expressed as:
Γ ( x 0 * , x i * ) = 1 n k = 1 n γ ( x 0 * ( k ) , x i * ( k ) )
where: Γ is the gray correlation degree, and n is the number of response characteristics.
The relative significance of each response characteristic may be different, and the simple averaging method of Equation (15) may lead to inaccurate gray correlations (GRG). Therefore, the weighted gray correlation measure (GRG) is often calculated by assigning different weights to the response characteristics by Equation (16).
Γ ( x 0 * , x i * ) = k = 1 n ω k γ ( x 0 * ( k ) , x i * ( k ) )
where, ω k is the weight of the kth response characteristic, k = 1 n ω k = 1 .

2.4. Entropy Method

The Entropy Method is an objective assignment method for obtaining weight coefficients and is widely used in various fields [42]. “Entropy” was originally a physical thermodynamic concept but was later developed by C E Shannon into the entropy theory of information theory, and the “entropy” defined by C E Shannon is called “information entropy.” Entropy is a measure of uncertain information. The smaller the amount of information, the greater the uncertainty, the larger the entropy accordingly, and the smaller the weight in calculating the comprehensive evaluation value, while the larger the amount of information, the smaller the uncertainty, the smaller the entropy accordingly, and the larger the weight in calculating the comprehensive evaluation value.
The specific steps of the entropy method are as follows.
Step 1: dimensionless processing
In order to eliminate the influence of the difference in magnitude and order of magnitude between the initial performance indicators on the results, the initial performance indicators need to be dimensionless, consistent with the normalization method in the previous section. The evaluation matrix P is:
p = ( p i j ) m × n = ( p 11 p 21 p m 1 p 12 p 22 p m 2 p 1 n p 2 n p m n )
where, p i j is the initial data after the normalization process.
Step 2: Calculation of index information entropy
According to the definition of information entropy, the indicator information entropy of E K is:
E k = H Σ i = 1 m f k i ln ( f k i )
where, H = 1 / ln ( m ) , f k i = p k i / Σ i = 1 m p k i , when f k i = 0 , let f k i ln ( f k i ) = 0 . m is the number of evaluation objects, f k i is the kth performance index under the ith evaluation object accounted for the proportion of the performance index.
Step 3: Calculation of the entropy weight of the index
The information entropy of each indicator is calculated from the indicator information entropy formula E 1 , E 2 ,…, E n , then the entropy weight of the kth performance indicator is:
W k = 1 E k n Σ k = 1 n E k
In the above equation, n is the number of evaluation indicators, 0 W k 1 and, Σ k = 1 n W k = 1 .

2.5. RBF Approximation Model

The Radial Basis Function (RBF) is a typical pre-feedback control algorithm, which mainly contains an input layer, an implicit layer, and an output layer, and its structure is shown in Figure 1. Among them, the input layer is mainly used to introduce variable information, and the number of variables determines the number of dimensions of the input layer. The implicit layer is the core component of the RBFNN algorithm, which maps the input information on the output layer according to certain mathematical relationships to a built-in nonlinear algorithm. The number of neurons in the implicit layer is directly related to the learning ability, fault tolerance, and approximation accuracy of RBF. The output layers output the corresponding structural response to a linear mapping function as the implied layer, and the output layers complete the optimization calculation process.
The most distinctive feature of the RBF method is that its implicit layer nodes use the Euclidean distance function as computational input, converting a multidimensional optimization problem into a one-dimensional optimization problem with distance as a variable. The RBF method uses the radial basis function as the transfer function and maps the radial distance r from the input sample x to the center of the implicit layer as a variable with the following specific mathematical expression.
{ f ( x ) = i = 1 n w m k v k ( x ) + b m v k ( x ) = x x k
where: f ( x ) is an approximation of the true response value f ( x ) ; x = [ x 1 , x 2 x n ] is a n-dimensional inputs vector; n is the number of basis functions; m is the number of output responses; w m k is the connection weight between the kth implied layer node and the mth output layers node; b m is the deviation of the mth output response; v k ( x ) is a radial basis function indicating the distance between the sample x and the ith sample x i in the design space.

2.6. Optimize the Design Process

Figure 2 shows the main process and content of this study.

3. Scissor Seat Suspension Dynamics Model

3.1. Seat Suspension Structure Form

The scissor seat suspension application tractor is shown in Figure 3a; the physical object is shown in Figure 3b–d. The 3D model is shown in Figure 3e, and the structural sketch is shown in Figure 3f. The seat suspension consists of upper and lower frames, scissor bars, springs, shock absorbers, and other components. Mutual articulation point E connects the seat suspension scissor rods. The scissor rod 3 is articulated with the seat’s lower frame of point B, and its upper-end D point can slide into the upper frame’s linear slots. The shear rod 4 articulates with the upper frame of the seat on pointing A, and its lower end C can slide into the linear slots of the seat’s lower frame. Spring 2 is a linear coil spring with stiffness k. Spring 2 is placed horizontally in the upper frame. Spring 2 is connected to point G at the upper end of the scissor rod 4, and the other end is connected to the right side of the upper frame of the seat suspension. Damper 6 is located between the upper and lower frames of the seat suspension and has a damping factor of c. The shear bar has the dimensions B D = l 1 + l 2 , A C = l 3 + l 4 , E A = l 4 , E B = l 1 = E C = l 3 , E D = l 2 . The angle between the shear bar and the bottom is ϕ .

3.2. Seat Suspension Mechanics Model

The following assumptions are made to simplify the equation’s calculation: ① Regardless of the influence of the seat cushion between the human body and the upper plate of the seat suspension, 75% of the human body mass and the upper plate quality of the seat suspension is called the bearing mass m. ② In simplifying the equation, assume that the scissor rods have no mass. ③ The seat under excitation produces vibrations in a static equilibrium position.
The scissor rod motion is shown in Figure 4 under excitation y 1 , with the seat base plate playing the role of displacement excitation y 1 = Q sin ω t , where Q is the excitation amplitude, and ω is the excitation circle frequency. The upper panel of the seat’s response are y 2 , The relationship between spring deformation Δ x , y 1 and y 2 is then as shown in Equation (21):
Δ x = l 0 { cos [ θ ( ϕ + Δ ϕ ) ] cos ( θ ϕ ) } ( l 3 + l 4 ) [ sin ( ϕ + Δ ϕ ) sin ϕ ] ( y 2 y 1 )
where ϕ is the angle in static equilibrium between the shear bar and the horizontal direction of the seat floor, and Δ ϕ is the change in angle under excitation.
Figure 5 depicts the overall structural force analysis of the seat suspension:
The inertial force on the upper plate of the seat suspension during the vibration is m y ¨ , and the equilibrium equation can be established using D’Alembert’s principle:
{ F B x + F C x F d cos ( γ + Δ γ ) = 0 F B y + F C y m y ¨ F d sin ( γ + Δ γ ) = 0 2 F C y l 1 cos ( ϕ + Δ ϕ ) m y ¨ 2 l 1 cos ( ϕ + Δ ϕ ) = 0 F d = c ( y ˙ 2 y ˙ 1 ) sin ( ϕ + Δ ϕ ) F k = k Δ x F C x = f d F C y
where F d is the damper’s damping force (N), γ is the angle between the axial direction of the damper and the horizontal direction of the upper plate of the suspension in static equilibrium (°). Δ γ is the change in damper tilt due to the excitation effect. f d is the sliding friction factor. F k is the spring force (N).
Figure 6 depicts the shear rod force analysis, and the equilibrium equations of the shear rods AC and BD can be established.
Force analysis of shear bar AC (see Figure 6a):
{ F C x + F E x F A x k Δ x = 0 F C y + F E y F A y = 0 F A x l 4 sin ( ϕ + Δ ϕ ) + F A y l 4 cos ( ϕ + Δ ϕ ) + F C y l 3 cos ( ϕ + Δ ϕ ) + F C x l 3 sin ( ϕ + Δ ϕ ) + k Δ x [ l 4 sin ( ϕ + Δ ϕ ) + l 0 sin ( θ ϕ Δ ϕ ) ] = 0
Force analysis of shear bar BD (see Figure 6b):
{ F B x + F D x F E x 1 = 0 F B y F E y 1 F D y = 0 F B x l 1 sin ( ϕ + Δ ϕ ) F B y l 1 cos ( ϕ + Δ ϕ ) F D x l 2 sin ( ϕ + Δ ϕ ) F D y l 2 cos ( ϕ + Δ ϕ ) = 0
When F E x = F E x 1 , F E y = F E y 1 , the Equations (23) and (24) is solved by:
F E x = l 1 , 2 F B x 2 l 2 l 1 , 2 cot ( ϕ + Δ ϕ ) F B y 2 l 2 l 3 , 4 F C x 2 l 4 l 3 , 4 cot ( ϕ + Δ ϕ ) F C y 2 l 4 k Δ x l 0 sin ( θ ϕ Δ ϕ ) 2 l 4 sin ( ϕ + Δ ϕ )
F E y = l 1 , 2 F B y 2 l 2 l 1 , 2 tan ( ϕ + Δ ϕ ) F B x 2 l 2 l 3 , 4 tan ( ϕ + Δ ϕ ) F C x 2 l 4 l 3 , 4 F C y 2 l 4 k Δ x l 0 sin ( θ ϕ Δ ϕ ) 2 l 4 cos ( ϕ + Δ ϕ )
From F D x = F D y f d and Equation (23) we have:
F B x F E x + ( F B y F E y ) f d = 0
According to assumption ③, Δ ϕ and Δ γ approaches zero, then ϕ + Δ ϕ ϕ , γ + Δ γ γ . Substitute F E x , F E y into Equation (27) to combine the above equations, and organize and simplify to obtain:
x 1 F B x + x 2 F B y + x 3 F C x + x 4 F C y + k Δ x l 0 sin ( θ ϕ ) 2 l 4 sin ϕ + k Δ x l 0 sin ( θ ϕ ) 2 l 4 cos ϕ f d = 0
In the formula, x 1 = ( 1 l 1 , 2 2 l 1 + l 1 , 2 2 l 2 tan ϕ f d ) , x 2 = ( l 1 , 2 2 l 1 cot ϕ l 1 , 2 2 l 2 f d + f d ) , x 3 = ( l 3 , 4 2 l 4 + l 3 , 4 2 l 4 tan ϕ f d ) , x 4 = ( l 3 , 4 2 l 4 cot ϕ + l 3 , 4 2 l 4 f d )
Simplifying the collation yields:
m y ¨ 2 + A 1 y ˙ 2 + A 2 y 2 = A 1 y ˙ 1 + A 2 y 1
In the formula, A 1 = c [ ( l 2 l 1 ) ( sin 2 γ + 2 f d sin 2 γ ) tan ϕ + l 1 , 2 ( 2 sin 2 γ + f d sin 2 γ tan 2 ϕ ) ] [ 2 l 1 , 2 ( f d tan ϕ + 1 ) ] , A 2 = k l 0 2 l 3 , 4 2 ( sin θ cos θ tan ϕ ) 2 .
The coefficients of Equation (28) are all constants derived from Equation (29), and the differential equation of motion of the seat are a linear constant coefficient differential equation. A 1 and A 2 represent the equivalent damping and equivalent stiffness of the seat suspension system, respectively, and their magnitudes are related to the seat’s structural parameters.

4. Seat Suspension Vibration Performance Impact Analysis

4.1. Main Performance Parameters of the Seat Suspension

According to the literature, the range of load-bearing masses acting on the seat suspension system is m (45 < m < 85) kg [43]. In order to study the effect of different masses on the vibration transfer characteristics of the seat suspension, the vibration transfer function study is carried out with a minimum value of 45 kg, a middle value of 70 kg and a maximum value of 85 kg as the representatives of the load-bearing mass. A certain type of scissor seat structure produced by a company is used as the research object, and the static equilibrium position of the seat when 70 kg acts on the seat suspension system is used as the initial position of the seat suspension system.

4.1.1. Stiffness Coefficient k

The inherent frequency of the seat vibration system f 0 depends on the spring stiffness k , the load mass m , the angle between the shear bar and the horizontal direction of the seat base of static equilibrium ϕ , and the geometric parameters of the shear bar, θ , β , l 0 and l 3 , 4 .
From Equation (29), the seat vibration damping inherent circular frequency ω 0 and inherent frequency f 0 can be obtained as:
ω 0 = A 2 m = l 0 l 3 , 4 ( sin θ cos θ tan ϕ ) k m
f 0 = ω 0 2 π = l 0 2 π l 3 , 4 ( sin θ cos θ tan ϕ ) k m
The inherent frequency f 0 of the seat suspension system should be (1 < f 0 < 2 Hz) [44], and the value of the seat suspension springs stiffness coefficient k is (68.36 < k < 273.44) N/mm when the bearing mass of the seat suspension is 70 kg, as obtained from Equation (31).

4.1.2. Damping Ratio ζ

The damping ratio ζ to seat vibration depends on the damping coefficient of the damper c , the stiffness of the horizontal linear spring k , the load mass m , the inclination angle γ of the damper at the static equilibrium position and the compression angle ϕ of the suspension, the sliding friction factor f d and the geometric parameters l 0 , l 1 , l 2 , l 1 , 2 and θ of the seat suspension.
From Equation (29), the damping ratio when the seat is vibrating is:
ζ = A 1 2 A 2 m = c [ ( l 2 l 1 ) ( sin 2 γ + 2 f d sin 2 γ ) tan ϕ + l 1 , 2 ( 2 sin 2 γ + f d sin 2 γ tan 2 ϕ ) ] [ 4 l 0 ( f d tan ϕ + 1 ) ( sin θ cos θ tan ϕ ) m k ]
The damper damping ratio ζ of the seat suspension system should be in the range of (0.18 < ζ < 0.35) [44]. From Equation (32), when the seat suspension load mass is 70 kg, the seat suspension damper damping coefficient c takes the value range of (1.4 < c < 5.5) N⋅s/mm.

4.1.3. Displacement s

Maximum displacement of the seat response surface is one of the performance evaluation parameters of seat suspension. According to the requirements of human comfort, the smaller the maximum displacement of the response surface, the better the vibration absorption performance. Let the action on the system excitation unit amplitude y 1 = e i ω t , and the response y 2 = H ( ω ) e i ω t , the two equations will be derived and brought into the Formula (29) to obtain, eliminating e i ω t , the complex frequency response function is:
H ( ω ) = ω 0 2 + i 2 ζ ω 0 ω ω 0 2 ω 2 + i 2 ζ ω 0 ω
From Equation (33), the modal value of the complex frequency response function is:
| H ( ω ) | = 1 + ( 2 ζ λ ) 2 ( 1 λ 2 ) 2 + ( 2 ζ λ ) 2
By:
y 1 = Q sin ω t
The displacement function of the plane vibration on the seat is obtained as:
y 2 ( t ) = Q 1 + ( 2 ζ λ ) 2 ( 1 λ 2 ) 2 + ( 2 ζ λ ) 2 sin ( ω t φ )
where λ is the frequency ratio, λ = ω / ω 0 ; φ is the phase difference angle between excitation and response. From Equations (32) and (35), the spring stiffness, k, and the damper damping, c, have a large effect on the magnitude of the performance evaluation parameter displacement, s.

4.1.4. Acceleration a

Performance evaluation parameters include the maximum acceleration of the plane response on the seat, and the smaller the maximum acceleration of the plane response on the seat, the better the vibration absorption performance. The acceleration calculation formula is shown in Equation (36).
The seat upper planes vibration response displacement from Equation (35) is y 2 ( t ) , and the seat upper planes response acceleration a = y ¨ 2 ( t ) :
a = Q ω 2 1 + ( 2 ζ λ ) 2 ( 1 λ 2 ) 2 + ( 2 ζ λ ) 2 sin ( ω t φ )
From Equations (32) and (36), it can be seen that in the seat suspension system, the spring stiffness, k, and the damper damping, c, have a large effect on the magnitude of the performance evaluation parameter, acceleration, a.

4.1.5. Transmission Rate η

Performance evaluation parameters include the seat vibration transmission rate; the smaller the seat vibration transmission rate, the better the vibration attenuation. From Equations (34) and (35), the transmission efficiency calculation formula Equation (37) is obtained.
The vibration transmittance η of the seat is the ratio of the response amplitude Q 1 + ( 2 ζ λ ) 2 ( 1 λ 2 ) 2 + ( 2 ζ λ ) 2 to the excitation amplitude Q , i.e.,:
η = 1 + ( 2 ζ λ ) 2 ( 1 λ 2 ) 2 + ( 2 ζ λ ) 2
From Equations (32) and (37), it can be seen that in the seat suspension system, the spring stiffness, k, and the damper damping, c, have a large influence on the magnitude of the performance evaluation parameters η.

4.2. Load Quality Impact Analysis

The parameters of the seat of the initial position are shown in Table 1. The sliding friction factor f d is taken as 0.3. The initial position of the seat is ensured as a static equilibrium position in different suspension loads masses by the preload of the coil springs of the seat suspension.
In order to make f 0 , ζ within a reasonable range of values, the values of seat suspension spring stiffness coefficient k and damper damping coefficient c at the equilibrium position are 172.28 N/mm and 3.37 N⋅s/mm, respectively. When the excitation amplitude is 2 mm, the input frequency is f (0 < f < 10) Hz, and the seat suspension carrying mass is 45 kg, 70 kg, and 85 kg. Substituting into Equations (31), (32), and (37), the seat vibration transmission characteristics under different seat suspension masses m can be obtained, as shown in Figure 7.
The seat suspension vibration performance parameters for different suspension masses are shown in Table 2.
Under the rule of satisfying the suspension inherent frequency f 0 (1 < f 0 < 2) Hz. The damping ratio ζ (0.18 < ζ < 0.35) constraints. Figure 7 and Table 2, show that the greater the load-bearing mass of three different seat suspensions, the greater the vibration amplitude in the vibration amplification area and the more obvious the vibration attenuation of the seat in the vibration isolation area.

4.3. Spring Stiffness Coefficient Influence Analysis

In the vibration isolation zone, when the input sinusoidal excitation amplitude is 2 mm and the input frequency is the seat suspension inherent frequency corresponding to different stiffnesses with fixed damping (c = 3.37 Ns/mm), respectively. The vibration performance parameters with different stiffness factors are shown in Table 3. The seat response acceleration and displacement curves (m = 70 kg) are shown in Figure 8a,b.
When the load-bearing mass is 70 kg, with a reasonable range of values for the inherent frequency f 0 and damper damping ratio ζ , the spring stiffness coefficient k is obtained from Equation (31) in the range of (68.36 < k < 273.44) N/mm. From Figure 8, it can be seen that in the range of k, with the increase of spring stiffness coefficient, the acceleration and displacement of the plane on the seat increase, and the vibration attenuation of the seat suspension decreases.

4.4. Damper Damping Influence Analysis

In the vibration isolation zone, when the input sinusoidal excitation amplitude is 2 mm and the input frequency is the seat suspension inherent frequency corresponding to different damping and fixed stiffness (k = 172.28 N/mm), respectively. The vibration performance parameters with different damping coefficients are shown in Table 4. The seat response acceleration and displacement curves (m = 70 kg) are shown in Figure 9a,b.
When the load-bearing mass is 70 kg, with a reasonable range of values for the inherent frequency f 0 and damper damping ratio ζ , the damper damping coefficient c is obtained from Equation (32) in the range of (1.4 < c < 5.5) N⋅s/mm. From Figure 9, it can be seen that in the range of value c, with the increase of damping coefficient, the acceleration and displacement of the plane on the seat decrease, and the vibration attenuation of the seat suspension is improved.

4.5. Analysis of Factors Influencing Different Performance Indicators

(1)
Main Effect Analysis
According to the main effect analysis formula for Section 2.1, Table 5 shows the table of the main factor effect of the control factor when the single optimization target is in the upper plane of the seat, mainly listing the analysis results in the optimization target response acceleration a, response displacement s and transfer rate η. Figure 10 shows the main effect plots of a, s, and η.
From Figure 10a,d, it can be seen that the main effect of acceleration a is that acceleration a increases with the increase of spring stiffness k, and damper damping c is non-linear to the fluctuation a. k plays a decisive role in a, followed by c. From Figure 10b,e, it can be seen that the main effect analysis regarding the displacement s yields that k becomes linearly and positively correlated with s. The effect of the damper c on s becomes nonlinear, with k having a decisive effect on s and c being the second most important. From Figure 10c,f, it can be seen that the effects of both k and c on the transfer rate η are volatile and become nonlinear. Therefore, the seat suspension system is a multi-factor interaction system. In response to the fluctuation of the relationship between both influencing factors k and c on η, it is not clear that k or c has a decisive influence on η. Contribution analysis is needed to verify further the results of the main effect analysis in Figure 10.
(2)
Contribution analysis
According to the contribution degree analysis formula for Section 2.2, 80 sample points were designed and calculated by the Latin hypercube method, and the contribution degree analysis method was selected to calculate the acceleration a, displacement s, and transmission rate η of the upper plane of the seat suspension as response indicators. The results of some sample points were listed as shown in Table 6.
From Figure 10, it can be seen that k has a decisive effect on both a and s; c has the second most important effect on a and s, consistent with the contribution analysis in Text 4.5(2). From Figure 11, c has a decisive effect on the transfer rate η, and k has the second most important effect on η.
In summary, corresponding to the seat suspension of a certain bearing mass, the suspension system stiffness coefficient k and damping c determine the seat vibration performance and also affect the response indexes such as seat vibration acceleration a, vibration displacement s, and transmission rate η. The stiffness coefficient k has a greater influence on the acceleration a and vibration displacement s, and the damping c has a greater influence on the transmission rate η. Response indices are interrelated and influenced, and one of the optimum may degrade other performance indices, resulting in a reduction in the overall seat vibration performance. Therefore, multi-objective optimization of seat vibration attenuation is required in order to obtain a compromise solution with excellent comprehensive performance.

5. Seat Suspension Vibration Attenuation Multi-Objective Optimization

5.1. Approximate Model Approach

The Optimal Latin Hypercube Sampling (OLHS) was used for the Design of Experiments (DoE) sampling to obtain experimental samples of different variable parameters (k, c). A total of 80 sample points were selected to fit the Kriging and RBF approximation models for each performance index.
In order to ensure the accuracy of multi-objective optimization, the accuracy of the approximation models needs to be verified. The higher the precision, the higher the confidence in the optimal solution. The coefficient R2 is often used as an evaluation index in engineering. The closer the R2 value is to 1, the higher the overall forecast accuracy of the approximate model is. Its mathematical expression is as follows.
R 2 = i = 1 n ( y ^ i y ¯ ) 2 i = 1 n ( y i y ¯ ) 2
where: n is the number of test sample points; y ^ i and y i are the predicted and actual response values of the approximate model corresponding to the ith sample point, respectively; y ¯ is the average of the actual responses to all sample points.
Within the range of values of design variables, another 20 sample points were selected using the central composite design. The accuracy of the Kriging and RBF approximation models were verified. It used the cross-validation method. The accuracy verification results of the approximation model are shown in Figure 12a–d.
As can be seen from Table 7, the value of the coefficient of determination R2 of the RBF proxy model is larger than that of the Kriging proxy model. Therefore, the RBF prediction model has better accuracy, and the RBF approximation model is used as the proxy model for the initial samples in this paper.

5.2. Multi-Objective Optimization Methods

5.2.1. Improved NSGA-II Algorithm

The NSGA-II was proposed by K Deb after introducing an elite retention strategy, fast non-dominated sorting method, and crowded distance comparison method based on the NSGA algorithm [45]. NSGA-II algorithm avoids the shortcomings of the NSGA algorithm due to the lack of diversity in the population due to shared fitness and improves the optimization efficiency and computational accuracy.
Although the crowded distance comparison method can maintain population diversity better, it is still inadequate for optimization problems with more than two objective functions. In order to further improve the population diversity of multi-objective optimization problems, a fixed threshold ε elimination strategy is used in MSGA-II to replace the crowding distance comparison method, which can solve the multi-objective optimization problems more reasonably. The principle of the MNSGA-II algorithm is shown in Figure 13.
The MNSGA-II algorithm operates as follows:
From Figure 14, we can see the operation steps of MNSGA-II algorithm.

5.2.2. Multi-Objective Optimal Design

Based on the constructed RBF proxy model on the multi-objective optimization platform Isight. Multi-objective optimization is carried out by using an improved NSGA-II optimization algorithm. The optimized population size is 20, evolutionary algebra is 200, and crossover probability is 0.9. The Pareto solution is obtained by optimization through 4001 iterations. As shown in Figure 15a shows.
The green dots in Figure 15 show the compromise solutions recommended by the optimization platform after the multi-objective optimization. According to Figure 15a–d, it can be seen that the acceleration a, displacement s, and transfer rate η, the three objective values, are interrelated and affect each other; one optimal often brings another worse performance, so the three selections are a result of compromise. Therefore, it is necessary to find a method to balance the performance of the three and to select the solution with the better comprehensive performance of the Pareto solution set, and the green solution recommended by the Isight platform does not meet the value range of the seat suspension damping value ζ, so it needs to be reselected. In this paper, the entropy-weighted gray correlation analysis ranking method is used to rank the comprehensive performance of Pareto solutions.

5.3. Entropy-Weighted Gray Correlation Ranking

The three objective functions correspond to the 300 Pareto frontier solutions of Figure 15a: a, s, and η, all of which are required to be as small as possible. The objective functions are first normalized to derive the gray correlation coefficient, and then the gray correlation degree of each scheme is obtained based on the objective weight values. The gray correlation coefficient and gray correlation degree values of each Pareto solution are calculated according to the formulas in Section 2.3 and Section 2.4. The results are shown in Table 8, which ranks the Pareto frontier solutions and shows the design scheme with the best overall performance. Figure 16 shows the gray correlation values of the Pareto frontier solutions. The gray correlation for the 203rd design solution has a maximum value of 0.928. The recommended solution corresponding to the maximum value of gray correlation is shown in Figure 15a.
The gray correlation degree value is based on Equation (16), and the GRG of each design solution can be obtained. Based on the magnitude of the GRG, each design solution can be ranked to obtain the best solution in terms of overall performance. The design solution to the highest GRG value is calculated to represent the solution to the best overall performance.
According to the gray correlation sorting, the sorting first is 203 group, the gray correlation value is 0.928, the corresponding k and c are substituted into the Formula (31) test, and the recommended solution is discarded because it does not meet the ζ taking value range. The second-ranked group is 165, and the gray correlation value is 0.882. At this time, k = 68.36 N/mm and c = 2.77 N⋅s/mm meet the requirements of the ζ-taking range and are selected as the preferred solution for the optimization of this paper.

5.4. Comparison of Vibration Vibration Attenuation before and after Optimization

To verify the effectiveness of the multi-objective optimization, the combination of the Pareto preferred solution k = 68.36 N/mm and c = 2.77 N⋅s/mm was substituted into Adams with the initial k = 172.28 N/mm and c = 3.37 Ns/mm combination of the seat suspension for simulation analysis and comparative verification. The comparison of seat vibration attenuation parameters before and after optimization is shown in Table 9, and the results of response indexes (a, s, η) are shown in Figure 17.
When the sinusoidal displacement excitation amplitude is 2 mm and the input frequency is the seat suspension natural frequency corresponding to different stiffness k and damping c, respectively. The results are shown in Figure 17a–c, where the acceleration a steady-state response value is significantly reduced, the displacement s steady-state response value is reduced, and the transmittance η is improved in general. The response indexes a, s, and η vibration attenuation are improved by 66.41%, 2.31%, and 8.19%, respectively, from Table 9.

6. Conclusions

This paper presents a multi-objective optimization design for seat suspension vibration attenuation based on an improved NSGA-II algorithm. A seat suspension dynamics model is established. The accuracy of the RBF approximation model is verified by the cross-validation method. Under different spring stiffness coefficients and damper damping coefficients, Pareto solutions are ordered by means of gray correlation analysis with entropic weight. The response index of plane vibration on the seat is analyzed, and the following conclusions are obtained:
(1)
The equivalent stiffness and equivalent damping of the seat suspension are derived. The equivalent stiffness is mainly related to the spring stiffness, and the equivalent damping is mainly related to the damper damping. The scissor seat suspension under micro-amplitude vibration conditions can be approximated as a linear vibration system. In the range of values of seat suspension inherent frequency f 0 and damping ratio ζ , the larger the seat suspension load mass, the larger the vibration amplitude of the seat on the vibration amplification zone and the greater the vibration attenuation in the vibration isolation zone. When the seat suspension load mass increased by 55.56% and 88.89%, the seat natural frequency decreased by 20.20% and 27.27%, the damping ratio decreased by 22.86% and 28.57%, and the seat vibration peak increased by 5.37% and 13.42%.
(2)
Within the constraint range of seat suspension inherent frequency f 0 and damping ratio ζ : under the condition of certain value of damper damping, as the value of spring stiffness increases, the value of acceleration and displacement of seat upper plane response increases, and the vibration attenuation ability of seat decreases. Under the condition of a certain value of spring stiffness, as the value of damper damping increases, the value of acceleration and displacement of seat upper plane response decreases, and the vibration attenuation of seat suspension improves.
(3)
Through main effect analysis and contribution analysis, the relationship between the control variable and response index is obtained. The effect of k on the response to a and s is decisive, and the effect of c on the response from a and s is second. c on the response to η is decisive, and the effect of k on the response of η is second.
(4)
RBF proxy model is constructed by combining experimental design sampling. Based on MNSGA-II multi-objective optimization algorithm, multi-objective optimization of seat suspension vibration performance is conducted, and Pareto frontier disaggregation is obtained. The entropy-weighted gray correlation ranking method is used to obtain the preferred solution that satisfies the range of variables. The effectiveness of the multi-objective optimization in this paper is verified by simulation analysis, and the response indexes a, s, and η vibration attenuation are improved by 66.41%, 2.31%, and 8.19%, respectively, which effectively improves the seat suspension vibration attenuation. The method used in this paper provides a reference for the study of vibration attenuation control of seat suspension.

Author Contributions

Conceptualization, S.Z. and W.W.; methodology, S.Z. and W.W.; software, W.W.; validation, L.X., S.Z. and X.C.; formal analysis, Y.C.; investigation, X.C.; resources, S.Z.; data curation, W.W.; writing—original draft preparation, W.W.; writing—review and editing, W.W.; visualization, W.W.; supervision, L.X.; project administration, L.X.; funding acquisition, L.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Major Science and Technology Project of Henan Province, grant number 221100240400, and by the Key Scientific Research Project of Colleges and Universities in Henan Province, grant number 20A460013.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rahman, A.; Ali, R.; Kabir, S.N.; Rahman, M.; Al Mamun, R.; Hossen, A. Agricultural Mechanization in Bangladesh: Status and Challenges towards Achieving the Sustainable Development Goals (SDGs). Ama-Agric. Mech. Asia Afr. Lat. Am. 2020, 51, 106–120. [Google Scholar]
  2. Chen, Y.C.; Chen, L.W.; Chang, M.Y. A Design of an Unmanned Electric Tractor Platform. Agriculture 2022, 12, 112. [Google Scholar] [CrossRef]
  3. Mano, Y.; Takahashi, K.; Otsuka, K. Mechanization in land preparation and agricultural intensification: The case of ricefrmng in the Cote d’Ivoire. Agric. Econ. 2020, 51, 899–908. [Google Scholar] [CrossRef]
  4. Brunetti, J.; D’Ambrogio, W.; Fregolent, A. Analysis of the Vibrations of Operators’ Seats in Agricultural Machinery Using Dynamic Substructuring. Appl. Sci. 2021, 11, 4749. [Google Scholar] [CrossRef]
  5. Han, H.-W.; Dae, K.K.; Hong, I.W.; Man, K.J.; Je, C.S.; Young-Jun, P. Analysis on Whole-Body Vibration of 100 kW Class Agricultural Tractor Operator for Evaluating Ride Comfort. J. Agric. Life Sci. 2022, 56, 135–142. [Google Scholar] [CrossRef]
  6. Li, Z.; Li, X.; Liu, B.; Wang, J. Influence of vehicle body vibration induced by road excitation on the performance of a vehicle-mounted piezoelectric-electromagnetic hybrid energy harvester. Smart Mater. Struct. 2021, 30, 55019. [Google Scholar] [CrossRef]
  7. Jung, H.; Park, G.; Kim, J.K. Piezoelectric-based dither control for automobile brake squeal suppression under various braking conditions. J. Vib. Control. 2021, 27, 2192–2204. [Google Scholar] [CrossRef]
  8. Huang, Y.; Tong, S.; Tong, Z.; Cong, F. Signal Identification of Gear Vibration in Engine-Gearbox Systems Based on Auto-Regression and Optimized Resonance-Based Signal Sparse Decomposition. Sensors 2021, 21, 1868. [Google Scholar] [CrossRef]
  9. Adam, S.A.; Jalil, N.A.; Rezali, K.M.; Ng, Y.G.; Sound and Vibration Research Group. The effect of posture and vibration magnitude on the vertical vibration transmissibility of tractor suspension system. Int. J. Ind. Ergon. 2020, 80, 103014. [Google Scholar] [CrossRef]
  10. Yoo, H.; Oh, J.; Chung, W.-J.; Han, H.-W.; Kim, J.-T.; Park, Y.-J.; Park, Y. Measurement of stiffness and damping coefficient of rubber tractor tires using dynamic cleat test based on point contact model. Int. J. Agric. Biol. Eng. 2021, 14, 157–164. [Google Scholar] [CrossRef]
  11. Oh, J.; Chung, W.J.; Han, H.W.; Kim, J.T.; Son, G.H.; Park, Y.J. Evaluation of Tractor Ride Vivrations by Cab Suspension System. Trans. Asabe 2020, 63, 1465–1476. [Google Scholar] [CrossRef]
  12. Desai, R.; Guha, A.; Seshu, P. A comparison of different models of passive seat suspensions. Proc. Inst. Mech. Eng. Part D-J. Automob. Eng. 2021, 235, 2585–2604. [Google Scholar] [CrossRef]
  13. Gao, P.; Liu, H.; Xiang, C.; Yan, P. Vibration reduction performance of an innovative vehicle seat with a vibration absorber and variable damping cushion. Proc. Inst. Mech. Eng. Part D-J. Automob. Eng. 2022, 236, 689–708. [Google Scholar] [CrossRef]
  14. Cvetanovic, B.; Cvetković, D.; Praščević, M.; Cvetković, M.; Pavlović, M. An analysis of the impact of agricultural tractor seat cushion materials to the level of exposure to vibration. J. Low Freq. Noise Vib. Act. Control. 2017, 36, 116–123. [Google Scholar] [CrossRef] [Green Version]
  15. Zhu, Z.; Yang, Y.; Wang, D.; Cai, Y.; Lai, L. Energy Saving Performance of Agricultural Tractor Equipped with Mechanic-Electronic-Hydraulic Powertrain System. Agriculture 2022, 12, 436. [Google Scholar] [CrossRef]
  16. Liao, X.; Du, X.; Li, S. Design of cab seat suspension system for construction machinery based on negative stiffness structure. Adv. Mech. Eng. 2021, 13, 16878140211044931. [Google Scholar] [CrossRef]
  17. Du, X.M.; Yu, M.; Fu, J.; Peng, Y.X.; Shi, H.F.; Zhang, H. H control for a semi-active scissors linkage seat suspension with magnetorheological damper. J. Intell. Mater. Syst. Struct. 2019, 30, 708–721. [Google Scholar] [CrossRef]
  18. Deng, L.; Sun, S.; Christie, M.; Li, W. Numerical Study of Rotary Magnetorheological Seat Suspension on the Impact Protection. Adv. Appl. Nonlinear Dyn. Vib. Control. 2021, 2022, 1003–1017. [Google Scholar]
  19. Zhang, X.; Wang, C.; Guo, K.; Sun, M.; Yao, Q. Dynamics modeling and characteristics analysis of scissor seat suspension. J. Vib. Control. 2017, 23, 2819–2829. [Google Scholar] [CrossRef]
  20. Deng, L.; Sun, S.; Christie, M.; Ning, D.; Jin, S.; Du, H.; Zhang, S.; Li, W. Investigation of a seat suspension installed with compact variable stiffness and damping rotary magnetorheological dampers. Mech. Syst. Signal Process. 2022, 171, 108802. [Google Scholar] [CrossRef]
  21. Zhao, Y.; Alashmori, M.; Bi, F.; Wang, X. Parameter identification and robust vibration control of a truck driver’s seat system using multi-objective optimization and genetic algorithm. Appl. Acoust. 2021, 173, 107697. [Google Scholar] [CrossRef]
  22. Maciejewski, I.; Blazejewski, A.; Pecolt, S.; Krzyzynski, T. A sliding mode control strategy for active horizontal seat suspension under realistic input vibration. J. Vib. Control. 2022, 0, 10775463221082716. [Google Scholar] [CrossRef]
  23. He, S.; Chen, K.; Xu, E.; Wang, W.; Jiang, Z. Commercial Vehicle Ride Comfort Optimization Based on Intelligent Algorithms and Nonlinear Damping. Shock. Vib. 2019, 2019, 2973190. [Google Scholar] [CrossRef]
  24. Deng, H.; Deng, J.; Yue, R.; Han, G.; Zhang, J.; Ma, M.; Zhong, X. Design and verification of a seat suspension with variable stiffness and damping. Smart Mater. Struct. 2019, 28, 65015. [Google Scholar] [CrossRef]
  25. Ning, D.; Du, H.; Sun, S.; Zheng, M.; Li, W.; Zhang, N.; Jia, Z. An Electromagnetic Variable Stiffness Device for Semiactive Seat Suspension Vibration Control. IEEE Trans. Ind. Electron. 2020, 67, 6773–6784. [Google Scholar] [CrossRef]
  26. Desai, R.; Guha, A.; Seshu, P. Modelling and simulation of an integrated human-vehicle system with non-linear cushion contact force. Simul. Model. Pract. Theory 2021, 106, 102206. [Google Scholar] [CrossRef]
  27. Jiang, R.; Jin, Z.; Ci, S.; Liu, D.; Sun, H. Reliability-based multilevel optimization of carbon fiber reinforced plastic control arm. Struct. Multidiscip. Optim. 2022, 65, 314. [Google Scholar] [CrossRef]
  28. Jiang, R.; Jin, Z.; Liu, D.; Wang, D. Multi-Objective Lightweight Optimization of Parameterized Suspension Components Based on NSGA-II Algorithm Coupling with Surrogate Model. Machines 2021, 9, 107. [Google Scholar] [CrossRef]
  29. Jiang, R.; Sun, T.; Liu, D.; Pan, Z.; Wang, D. Multi-Objective Reliability-Based Optimization of Control Arm Using MCS and NSGA-II Coupled with Entropy Weighted GRA. Appl. Sci. 2021, 11, 5825. [Google Scholar] [CrossRef]
  30. Pandey, M.; Regis, R.G.; Datta, R.; Bhattacharya, B. Surrogate-assisted multi-objective optimization of the dynamic response of a freight wagon fitted with three-piece bogies. Int. J. Rail Transp. 2021, 9, 290–309. [Google Scholar] [CrossRef]
  31. Yang, L.; Wang, R.; Sun, Z.; Meng, X.; Zhu, Z. Multi-objective optimization design of hydropneumatic suspension with gas-oil emulsion for ride comfort and handling stability of an articulated dumper truck. Eng. Optim. 2021, 24, 1–20. [Google Scholar] [CrossRef]
  32. Grotti, E.; Mizushima, D.M.; Backes, A.D.; Awruch, M.D.D.F.; Gomes, H.M. A novel multi-objective quantum particle swarm algorithm for suspension optimization. Comput. Appl. Math. 2020, 39, 105. [Google Scholar] [CrossRef]
  33. Zhao, Y.; Xu, H.; Deng, Y.; Wang, Q. Multi-objective optimization for ride comfort of hydro-pneumatic suspension vehicles with mechanical elastic wheel. Proc. Inst. Mech. Eng. Part D-J. Automob. Eng. 2019, 233, 2714–2728. [Google Scholar] [CrossRef]
  34. Prasad, V.; Pawaskar, D.N.; Seshu, P. Controller design and multi-objective optimization of heavy goods vehicle suspension system by geometry-inspired GA. Struct. Multidiscip. Optim. 2021, 64, 89–111. [Google Scholar] [CrossRef]
  35. Hua, Y.; Zhu, H.; Xu, Y. Multi-Objective Optimization Design of Bearingless Permanent Magnet Synchronous Generator. IEEE Trans. Appl. Supercond. 2020, 30, 1–5. [Google Scholar] [CrossRef]
  36. Varela, M.; Gyi, D.; Mansfield, N.; Picton, R.; Hirao, A.; Furuya, T. Engineering movement into automotive seating: Does the driver feel more comfortable and refreshed? Appl. Ergon. 2019, 74, 214–220. [Google Scholar] [CrossRef] [Green Version]
  37. Nagarkar, M.P.; Patil, G.J.V.; Patil, R.N.Z. Optimization of nonlinear quarter car suspension-seat-driver model. J. Adv. Res. 2016, 7, 991–1007. [Google Scholar] [CrossRef] [Green Version]
  38. Kim, J.H.; Dennerlein, J.T.; Johnson, P.W. The effect of a multi-axis suspension on whole body vibration exposures and physical stress in the neck and low back in agricultural tractor applications. Appl. Ergon. 2018, 68, 80–89. [Google Scholar] [CrossRef]
  39. Kim, J.H.; Marin, L.S.; Dennerlein, J.T. Evaluation of commercially available seat suspensions to reduce whole body vibration exposures in mining heavy equipment vehicle operators. Appl. Ergon. 2018, 71, 78–86. [Google Scholar] [CrossRef]
  40. Liu, P.; Ning, D.; Luo, L.; Zhang, N.; Du, H. An Electromagnetic Variable Inertance and Damping Seat Suspension with Controllable Circuits. IEEE Trans. Ind. Electron. 2022, 69, 2811–2821. [Google Scholar] [CrossRef]
  41. Zheng, E.; Fan, Y.; Zhu, R.; Zhu, Y.; Xian, J. Prediction of the vibration characteristics for wheeled tractor with suspended driver seat including air spring and MR damper. J. Mech. Sci. Technol. 2016, 30, 4143–4156. [Google Scholar] [CrossRef]
  42. Zhang, S.; Liu, K.; Zhang, S.; Xu, L. Study on Aerodynamic Performance and Lightweight Multiobjective Optimization Design of Wheel With Entropy Weighted Grey Relational Analysis Methods. IEEE Access 2022, 10, 93421–93438. [Google Scholar] [CrossRef]
  43. Mehta, C.R.; Tewari, V.K. Biomechanical model to predict loads on lumbar vertebra of a tractor operator. Int. J. Ind. Ergon. 2015, 47, 104–116. [Google Scholar] [CrossRef]
  44. Yi-jie, S.H.U.I.; Rakheja, S.; Shangguan, W.B. Calculation and analysis for equivalent stiffness and damping of a scissor type seat-seat suspension-human body system. J. Vib. Shock. 2016, 35, 38–44. [Google Scholar]
  45. Verma, S.; Pant, M.; Snasel, V. A Comprehensive Review on NSGA-II for Multi-Objective Combinatorial Optimization Problems. IEEE Access 2021, 9, 57757–57791. [Google Scholar] [CrossRef]
Figure 1. Radial basis neural network structure diagram.
Figure 1. Radial basis neural network structure diagram.
Agriculture 13 00048 g001
Figure 2. Flow chart of seat suspension vibration attenuation optimization design.
Figure 2. Flow chart of seat suspension vibration attenuation optimization design.
Agriculture 13 00048 g002
Figure 3. Structure diagram of a type of scissor seat suspension. (a) Seat suspension for tractors. (b) Seat suspension shall be applied to the seat. (c) Picture of seat suspension (front view). (d) Picture of seat suspension (left view). (e) 3D model diagram of seat suspension. (f) Structure diagram of seat suspension.
Figure 3. Structure diagram of a type of scissor seat suspension. (a) Seat suspension for tractors. (b) Seat suspension shall be applied to the seat. (c) Picture of seat suspension (front view). (d) Picture of seat suspension (left view). (e) 3D model diagram of seat suspension. (f) Structure diagram of seat suspension.
Agriculture 13 00048 g003
Figure 4. Shear y 1 bar motion under excitation.
Figure 4. Shear y 1 bar motion under excitation.
Agriculture 13 00048 g004
Figure 5. Overall force analysis of the seat.
Figure 5. Overall force analysis of the seat.
Agriculture 13 00048 g005
Figure 6. Shear bar force analysis diagram. (a) AC rod force analysis. (b) BD rod force analysis.
Figure 6. Shear bar force analysis diagram. (a) AC rod force analysis. (b) BD rod force analysis.
Agriculture 13 00048 g006
Figure 7. Vibration transmission characteristics under different seat suspension load mass.
Figure 7. Vibration transmission characteristics under different seat suspension load mass.
Agriculture 13 00048 g007
Figure 8. Seat response acceleration curve with different spring stiffness values. (a) Acceleration Response. (b) Displacement Response.
Figure 8. Seat response acceleration curve with different spring stiffness values. (a) Acceleration Response. (b) Displacement Response.
Agriculture 13 00048 g008
Figure 9. Acceleration and displacement curves of seat response with different damper damping values. (a) Acceleration Response. (b) Displacement Response.
Figure 9. Acceleration and displacement curves of seat response with different damper damping values. (a) Acceleration Response. (b) Displacement Response.
Agriculture 13 00048 g009
Figure 10. Main effect analysis diagram. (a) Influence of spring stiffness k on acceleration a. (b) Influence of spring stiffness k on displacement s. (c) Spring stiffness k to transmissibility η Impact analysis of. (d) Influence of damper damping c on acceleration a. (e) Influence of damper damping c on displacement s. (f) Damper damping c pair transmissibility η Impact of.
Figure 10. Main effect analysis diagram. (a) Influence of spring stiffness k on acceleration a. (b) Influence of spring stiffness k on displacement s. (c) Spring stiffness k to transmissibility η Impact analysis of. (d) Influence of damper damping c on acceleration a. (e) Influence of damper damping c on displacement s. (f) Damper damping c pair transmissibility η Impact of.
Agriculture 13 00048 g010aAgriculture 13 00048 g010b
Figure 11. Analysis of the contribution of each indicator.
Figure 11. Analysis of the contribution of each indicator.
Agriculture 13 00048 g011
Figure 12. RBF and Kriging fitting model for vibration attenuation index. (a) Accuracy verification of acceleration a approximate model. (b) Accuracy verification of displacement s approximate model. (c) Transmissibility η Accuracy verification of approximate model. (d) Accuracy verification of approximate model of natural frequency f0.
Figure 12. RBF and Kriging fitting model for vibration attenuation index. (a) Accuracy verification of acceleration a approximate model. (b) Accuracy verification of displacement s approximate model. (c) Transmissibility η Accuracy verification of approximate model. (d) Accuracy verification of approximate model of natural frequency f0.
Agriculture 13 00048 g012
Figure 13. Principle of MNSGA-II algorithm.
Figure 13. Principle of MNSGA-II algorithm.
Agriculture 13 00048 g013
Figure 14. MNSGA-II algorithm step flow chart.
Figure 14. MNSGA-II algorithm step flow chart.
Agriculture 13 00048 g014
Figure 15. RBF optimization of each vibration attenuation index. (a) Preferred solutions for maximum acceleration, maximum displacement, and transfer rate. (b) Preferred solution for maximum acceleration and maximum displacement. (c) Preferred solution for maximum acceleration and transfer rate. (d) Preferred solution for maximum displacement and transfer rate.
Figure 15. RBF optimization of each vibration attenuation index. (a) Preferred solutions for maximum acceleration, maximum displacement, and transfer rate. (b) Preferred solution for maximum acceleration and maximum displacement. (c) Preferred solution for maximum acceleration and transfer rate. (d) Preferred solution for maximum displacement and transfer rate.
Agriculture 13 00048 g015aAgriculture 13 00048 g015b
Figure 16. Gray correlation of the Pareto front for vibration attenuation.
Figure 16. Gray correlation of the Pareto front for vibration attenuation.
Agriculture 13 00048 g016
Figure 17. Comparison of each response index before and after optimization. (a) Comparison of acceleration before and after optimization. (b) Comparison of displacement before and after optimization. (c) Comparison of delivery rate before and after optimization. (d) Comparison of acceleration before and after optimization on frequency.
Figure 17. Comparison of each response index before and after optimization. (a) Comparison of acceleration before and after optimization. (b) Comparison of displacement before and after optimization. (c) Comparison of delivery rate before and after optimization. (d) Comparison of acceleration before and after optimization on frequency.
Agriculture 13 00048 g017
Table 1. Seat suspension static balance position structure parameters.
Table 1. Seat suspension static balance position structure parameters.
l 0 /mm l 1 /mm l 2 /mm l 3 /mm l 4 /mm θ /(°) γ /(°) ϕ /(°)
5013210713212011520.611.6
Table 2. Vibration performance parameters under different seat suspension load quality.
Table 2. Vibration performance parameters under different seat suspension load quality.
Load Capacity of Different Suspensions/kgVibration Transfer Function Peak Frequency/HzVibration Transfer Function Peak Amplitude/mmInherent Frequency/HzDamping RatioVibration Transmission Rate at 3 Hz
m = 451.782.981.980.350.87
m = 701.283.141.580.270.51
m = 851.003.381.440.250.41
Table 3. Vibration performance parameters at different spring stiffness factors.
Table 3. Vibration performance parameters at different spring stiffness factors.
Different Stiffnesses k/(N/mm)Maximum Response Acceleration/(mm/s2)Maximum Response Displacement/mmVibration Transmission RateSeat Suspension Inherent Frequency/Hz
k = 7087.02114.481.7081.01
k = 170255.12115.611.3361.58
k = 272469.13116.391.0831.99
Table 4. Vibration performance parameters under different damping coefficients of dampers.
Table 4. Vibration performance parameters under different damping coefficients of dampers.
Different Damping c/(N⋅s/mm)Maximum Response Acceleration/(mm/s2)Maximum Response Displacement/mmVibration Transmission RateSeat Suspension Inherent Frequency/Hz
c = 1.9352.05116.581.2281.59
c = 3.7248.36115.521.0431.59
c = 5.5223.19115.221.2721.59
Table 5. Response results of each index under different variables.
Table 5. Response results of each index under different variables.
Main Effect Evaluation IndicatorsAcceleration
a/(mm/s2)
Displacement
s/mm
Transmission
Rate/
kckckc
Level1128.964391.331114.821116.6501.8533.651
2182.668235.396115.262115.4511.9022.294
3252.884302.738115.684115.6742.2512.062
4335.443215.925116.062115.3222.4201.671
5464.419218.989116.561115.3102.8031.541
The main
Effect value
335.455175.4061.7401.3400.9502.110
Sorting121212
Interaction effect value132.6720.5800.705
The overall
average
272.88115.682.24
Table 6. Results of experimental design of contribution.
Table 6. Results of experimental design of contribution.
No.Design VariablesPerformance Response
k(N/mm)c(N⋅s/mm)a(mm/s2)s(mm)η
195.0554.394119.380114.9301.449
2261.1203.403442.620116.3102.457
3264.0083.403442.620116.3102.457
7973.9904.694104.120114.8101.323
80258.4151.648712.660118.2004.718
Table 7. Proxy model accuracy evaluation index.
Table 7. Proxy model accuracy evaluation index.
Performance
Response
Kriging Surrogate ModelRBF Surrogate Model
R2R2
a/(mm/s2)0.98510.9921
s/mm0.96240.9772
η/0.97730.9928
f 0 /0.98350.9932
Table 8. Gray correlation coefficients and gray correlation values.
Table 8. Gray correlation coefficients and gray correlation values.
No.Number of Gray Correlations for Each Response IndexGray CorrelationSort by
a(mm/s2)s(mm)η
10.5160.6670.6320.589179
20.5240.5600.6920.568186
1650.9260.8870.7630.8822
2030.9670.8960.8760.9281
3000.7630.8530.6340.76482
Table 9. Comparison of seat vibration attenuation parameters before and after optimization.
Table 9. Comparison of seat vibration attenuation parameters before and after optimization.
Vibration Performance Parametersk/(N/mm), c/(N⋅s/mm) Before and After OptimizationPerformance Improvement Rate/%
k = 172.28, c = 3.37k = 68.36, c = 2.77
Maximum response acceleration/(mm/s2)2598766.41
Maximum response displacement/(mm)115.64112.972.31
Rate of vibration transmission/1.2571.3608.19
Seat suspension inherent frequency/Hz1.561.0035.90
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, S.; Wei, W.; Chen, X.; Xu, L.; Cao, Y. Vibration Performance Analysis and Multi-Objective Optimization Design of a Tractor Scissor Seat Suspension System. Agriculture 2023, 13, 48. https://doi.org/10.3390/agriculture13010048

AMA Style

Zhang S, Wei W, Chen X, Xu L, Cao Y. Vibration Performance Analysis and Multi-Objective Optimization Design of a Tractor Scissor Seat Suspension System. Agriculture. 2023; 13(1):48. https://doi.org/10.3390/agriculture13010048

Chicago/Turabian Style

Zhang, Shuai, Weizhen Wei, Xiaoliang Chen, Liyou Xu, and Yuntao Cao. 2023. "Vibration Performance Analysis and Multi-Objective Optimization Design of a Tractor Scissor Seat Suspension System" Agriculture 13, no. 1: 48. https://doi.org/10.3390/agriculture13010048

APA Style

Zhang, S., Wei, W., Chen, X., Xu, L., & Cao, Y. (2023). Vibration Performance Analysis and Multi-Objective Optimization Design of a Tractor Scissor Seat Suspension System. Agriculture, 13(1), 48. https://doi.org/10.3390/agriculture13010048

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop