Next Article in Journal
MixerNet-SAGA A Novel Deep Learning Architecture for Superior Road Extraction in High-Resolution Remote Sensing Imagery
Next Article in Special Issue
Physiological Responses Related to Sitting Comfort Due to Changes in Seat Parameters
Previous Article in Journal
Power Drive Architectures for Industrial Hydraulic Axes: Energy-Efficiency-Based Comparative Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of a 7-DOF Biodynamic Model for a Seated Human and a Hybrid Optimization Method for Estimating Human-Seat Interaction Parameters

1
Department of Mechanical Engineering, Kyungpook National University, Daegu 41566, Republic of Korea
2
Korea Aerospace Research Institute, KSLV-II R&D Directorate, Daejeon 34133, Republic of Korea
3
School of Mechanical Engineering, Kyungpook National University, Daegu 41566, Republic of Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(18), 10065; https://doi.org/10.3390/app131810065
Submission received: 31 July 2023 / Revised: 28 August 2023 / Accepted: 4 September 2023 / Published: 6 September 2023
(This article belongs to the Special Issue Seating Comfort and Biomechanical Application)

Abstract

:
Existing biodynamic models adopt apparent mass and seat-to-head transmissibility to predict the response of seated humans to whole-body vibration, limiting their ability to capture the actual response of distinct body segments in different excitation conditions. This study systematically develops a 7-DOF seated human model, a vibration experiment, and a novel hybrid optimization to estimate unknown mechanical parameters and predict the response of different human body segments to vertical vibrations. Experimental results showed that the upper trunk and head were most susceptible to transmitted vibrations. Combining the 7-DOF model and HOM resulted in accelerated optimization, improved numerical stability, and significant minimization of the objective function value compared to conventional algorithms. Notably, the estimated parameters, particularly stiffness, remained consistent regardless of increasing excitation magnitude or change in the body segment data used. Additionally, the model captured the non-linearity in human biodynamics through stiffness softening. These findings are applicable in seating systems optimization for comfort and safety.

1. Introduction

The study of the effects of whole-body vibration (WBV) on human health is a topic of significant interest, particularly due to the increased exposure of individuals to whole-body vibration resulting from vehicular transport. Prolonged exposure to WBV can lead to discomfort and adversely affect the health of passengers, particularly in a constrained sitting posture [1,2]. In this regard, it is essential to investigate how vibration is transmitted to the human body, especially in seated posture.
The region of the resonant frequency of vibration that affects passengers is of critical importance. The effect of WBV on the seated human body is more dominant under low-frequency excitation, specifically below 30 Hz [3]. Several studies have reported the principal resonant frequencies within less than 13 Hz, with the dominant peak in the frequency range of 4–6 Hz under vertical vibrations [4,5,6]. Therefore, studying the effects of vertical vibrations in regions within 0–15 Hz on different body segments is important to design systems that ensure that the human natural frequencies do not coincide with excitation frequency within this range.
One effective technique for investigating the dynamic response of humans to WBV is the development and use of biodynamic seated human models with different postural and excitation conditions. Biodynamic models are used to augment experimental measurements, which are not only time-consuming but limited by safety considerations. Existing biodynamic models are primarily categorized into finite element (FE), multibody (MB), and lumped parameter (LP) models. While complex models like finite elements and multibody systems may provide a more detailed representation of human body dynamics, a lumped-parameter model can be sufficient for investigating the dynamic response of humans to vibration. The simplicity and computational efficiency of lumped-parameter models make them ideal for exploring different scenarios and investigating the effects of different parameters on the system’s response [6]. Furthermore, lumped-parameter models have been extensively validated against experimental data and can be customized to include more detailed features of the human body, making them a viable option for investigating the dynamic response of humans to vibration [7,8,9,10].
The estimation of appropriate human parameters using biodynamics models is another widely reported issue in WBV studies. It is highly difficult to measure human mechanical parameters such as stiffness and damping coefficient from live tissues. Likewise, seats are made from viscoelastic materials whose mechanical properties are also difficult to measure [9,11,12]. Therefore, researchers estimate these viscoelastic properties through the minimization of errors between the biodynamic model responses and a target response function derived from the measured responses from experiments. Most existing studies consider measured vertical apparent mass (AM) as a target response function due to the rapid convergence of the AM optimization problem. A few studies also adopted seat-to-head transmissibility (STHT) in the parameter’s optimization problem [6,8]. However, model parameter identification with AM does not often describe the contributions of low inertia upper body segments to the total response [13]. Studies by Cho and Yoon [8] and Matsumo and Griffin [14] established that transmissibility is more appropriate for predicting transmitted vibration to selected body segments. Therefore, it is necessary to directly measure the response of each body segment to whole-body vibration in order to access the distinct response of the body segments to vibration, and so far, there are limited studies in this area.
Existing studies adopt either gradient-based or evolutionary-based algorithms to estimate the parameters of biodynamic models. However, there is still an exhaustive search for the appropriate optimization method that produces reliable estimates of these human parameters. For instance, Bae and Kang, Shao et al., and Huang et al. (2020) [10,15,16] optimized the parameters of 2-, 3-, 4-, and 5-DOF models using the gradient-based MATLAB function, called “fmincon”. Zhang et al. [17] introduced a reference vector-guided evolutionary algorithm to obtain the parameters of a 7-DOF human model, whereas Zhao et al. [18] explored multi-objective and genetic algorithms (GAs) for parameter identification of a 5-DOF model. From these studies, it was found that the separate or hybrid application of these optimization techniques is usually associated with trade-offs that include an exhaustive search for initial values, randomly generated optimal solutions, high computation time, and generation of estimated parameters that do not reflect the dynamic response of the human body to WBV [6,19]. Therefore, there is still a gap in research to develop model parameter estimation algorithms that are both reliable and computationally efficient.
In this study, we developed a lumped parameter seven-degree-of-freedom (7-DOF) biodynamic model for a seated human under vertical vibrations, including rotational degrees-of-freedom (DOF) and seat (and backrest) modeled as massless torsional and translational spring and damper. We designed a vibration experiment to obtain the transmissibility of different body segments of a seated subject under random base excitation in the vertical direction. A new hybrid optimization method (HOM) that combines gradient-based and genetic algorithms in some systematically defined steps was developed to estimate unknown stiffness and damping parameters of the model. We quantitatively evaluate the model fitting performances of the optimization methods using goodness-of-fit criteria, objective function value at the final solution, and computation speed. The sensitivities of the estimated parameters to the appearance and magnitude of the resonance frequencies were studied based on the relationship between the increment of the excitation magnitudes (0.5 to 2.0 m/s2 rms) and using the transmissibility of different body segments in the optimization problem. Important findings on the key parameters affecting the region of resonance and how vibration is transmitted to each body segment were discussed. Ultimately, the development of the 7-DOF biodynamic model and HOM advances research on the optimization of comfort and safety for individuals when riding in vehicles.

2. Development of the 7-DOF Human Model

2.1. Model Description

The proposed 7-DOF biodynamic model consists of lumped parameters, including mass, spring, and damper elements. The model considers distinct body segments, comprising the upper body parts (including the head, neck, trunk, and upper extremities) and the lower body parts (including the pelvis, thigh, knee, leg, ankle, and foot). Each segment is treated as nonoverlapping, with the center of mass identified close to the proximal joints. The transmission of vibration from the seat to other segments is influenced by the viscoelastic properties of the human-seat coupling. Viscoelastic seats exhibit higher resonance frequencies compared to rigid seats [12,20]. To account for this, we represent the human-seat coupling using effective translational spring-damper combinations attached to the trunk and thigh, mimicking the viscoelastic properties of human skin and seat foam. The coupling between the shoe, foot skin, and foot segment is also modeled as an effective translational spring damper. These spring-damper combinations are attached perpendicularly to their respective body segments. Additionally, the joints at the hip, neck, knee, and ankle are modeled as torsional spring-damper elements, allowing rotation in the x-z plane.
Figure 1 illustrates the 7-DOF seated human model, consisting of five rigid bodies connected by dynamic joints and supported by a backrest. The model is subjected to vertical excitation, with m 1 , m 2 , m 3 , m 4 , m 5 representing the mass and J 1 , J 2 , J 3 , J 4 , J 5 representing the moment of inertia of the head, trunk, thigh, leg, and foot, respectively. The lengths of the segments are denoted as L 1 , L 2 , L 3 , L 4 , L 5 , and the trunk, thigh, and foot have thicknesses T 2 , T 3 , T 5 . The translational spring-damper elements attached to the trunk, thigh, and foot are characterized by stiffness coefficients k 21 , k 22 , k 31 , k 32 , k 51 and damping coefficients c 21 , c 22 , c 31 , c 32 , c 51 . The distances between the joints and the spring-damper combinations on trunk, thigh, and foot are represented by d 21 , d 22 , d 31 , d 32 , and d 51 . The torsional spring-damper elements at the joints have stiffness coefficients k 1 T , k 2 T , k 3 T , k 4 T and damping coefficients c 1 T , c 2 T , c 3 T , c 4 T and z b is the seat base excitation.

2.2. Derivation of Model Equations of Motion

The equations governing motions in the 7-DOF human model were derived through the Lagrangian dynamics approach in the form:
d d t K E q ˙ k K E q k + P E q k + D q ˙ k = 0 k = 1 ,   2 , , 7
where q k = x h , z h , θ 1 , θ 2 , θ 3 , θ 4 , θ 5 are the generalized coordinates corresponding to the horizontal and vertical displacement of the hip joint and rotational displacement of the center of gravity of each body segment of the 7-DOF model.
The kinetic energy (KE) of the system is given as:
K E = 1 2 i = 1 5 m i x ˙ i 2 + z ˙ i 2 + 1 2 i = 1 5 J i θ ˙ i 2 ,
where x i and z i are the position vectors of the center of gravity of the body segments in the x–z plane derived as
x 1 = x h + L 1 2 cos θ 1 + L 2 cos θ 2 ,   x 2 = x h + L 2 2 cos θ 2 ,
x 3 = x h + L 3 2 cos θ 3 ,   x 4 = x h + L 3 cos θ 3 + L 4 2 cos θ 4
x 5 = x h + L 3 cos θ 3 + L 4 cos θ 4 + L 5 2 cos θ 5 ,
z 1 = z h + L 1 2 sin θ 1 + L 2 sin θ 2 ,   z 2 = z h + L 2 2 sin θ 2 ,
z 3 = z h + L 3 2 sin θ 3 ,   z 4 = z h + L 3 sin θ 3 + L 4 2 sin θ 4 ,
z 5 = z h + L 3 sin θ 3 + L 4 sin θ 4 + L 5 2 sin θ 5
Equations of the potential energy (PE) and Rayleigh’s dissipation function (D) of the system are given as:
P E = 1 2 i = 2 , j = 1 5 k i j δ i j 2 + 1 2 i = 1 4 k i T ϕ i 2 + i = 1 5 m i g z i ,
D = 1 2 i = 2 , j = 1 5 c i j δ ˙ i j 2 + 1 2 i = 1 5 c i T ϕ ˙ i 2 ,
where δ i j and ϕ i are the deflections of the translational and torsional springs; g represents gravity; and z i is the change in the value of z-axis position vector of the center of gravity of body segments.
The deformation of the 7-DOF model under the influence of the vertical excitation, z b is presented in Figure 2. The deflection terms, δ i j and ϕ i are taken as the displacement of the position of the translational and torsional spring-damper as the model vibrates from the undeformed to the deformed configurations. Similarly, the position of the center of gravity, P i with coordinates, x i and z i during vibrations is considered. The deflection of the torsional spring damper, ϕ i is derived as the angular deflection during vibrations in the form:
ϕ 1 = θ 1 θ 1 * θ 2 θ 2 * ,   ϕ 2 = θ 2 θ 2 * θ 3 θ 3 * ,
ϕ 3 = θ 4 θ 4 * θ 3 θ 3 * ,   ϕ 4 = θ 5 θ 5 * θ 4 θ 4 * ,
where * denotes the initial value of the angle before vibration. On the other hand, the deflection of the translational spring-damper, δ i j is taken as the parallel displacement of the ends of the translational spring-damper with respect to the undeformed configuration. The coordinates of the ends of the spring damper are derived by setting the hip joint as the reference point.
For instance, the derivation of deflection of the spring-damper at the thigh segment is described in Figure 3a. The position of the hip joint before and after deformation is denoted as C 31 * x h * , z h * and C 31 x h , z h ; B 31 * and B 32 * are the position of the end of the first and second spring-damper in the undeformed configuration and B 31 and B 32 in the deformed configuration attached to the thigh; A 31 * and A 32 * , are the fixed ends of the spring-damper; 31 * and 32 * , are the unstretched length of the first and second spring-damper. The asterisk * represent the undeformed state for each of these points. Considering the first spring damper in Figure 3a, δ 31 is the parallel spring extension from B 31 * A 31 * to B 31 A 31 * with respect to the undeformed configuration, where x–z coordinates for points, A 31 * x 31 * , z 31 * and B 31 x 31 , z 31 are derived as
x 31 * = x h * + T 3 2 sin θ 3 * + d 31 cos θ 3 * + 31 * sin θ 3 * ,
z 31 * = z h * T 3 2 cos θ 3 * + d 31 sin θ 3 * 31 * cos θ 3 * ,
x 31 = x h + T 3 2 sin θ 3 + d 31 cos θ 3 ,
z 31 = z h T 3 2 cos θ 3 + d 31 sin θ 3 + z b
and by rotating the thigh through θ 3 * , we can see that δ 31 is the vertical axis (z-axis) displacement of the spring. Clockwise rotation transformation is applied to B 31 A 31 * in the form
x 31 x 31 * z 31 z 31 * = cos θ 3 * sin θ 3 * sin θ 3 * cos θ 3 * x 31 x 31 * z 31 z 31 *
where prime represents transformed terms of B 31 A 31 * . The transformed vertical axis displacement is thus derived as
δ 31 = x 31 x 31 * sin θ 3 * + z 31 z 31 * cos θ 3 * Δ 31 * = z h z h * cos θ 3 * x h x h * sin θ 3 * T 3 2 cos θ 3 θ 3 * + d 31 sin θ 3 θ 3 * + T 3 2 + z b cos θ 3 *
The transformation matrix technique was applied to derive an expression for δ 32 such that all terms in Equation (13) were retained except d 31 that changed to d 32 . A transformation rotation through θ 5 * was applied to the foot region, as shown in Figure 3b. The expressions for coordinates of the points, A 51 * x 51 * , z 51 * and B 51 x 51 , z 51 applied in the derivation of δ 51 are as follows:
x 51 * = x h * + T 3 2 sin θ 3 * + L 3 cos θ 3 * + L 4 cos θ 4 * + T 5 2 sin θ 5 * + d 51 cos θ 5 * + 51 * sin θ 5 * ,
z 51 * = z h * T 3 2 cos θ 3 * + L 3 sin θ 3 * L 4 sin θ 4 * T 5 2 cos θ 5 * + d 51 sin θ 5 * 51 * cos θ 5 * ,
x 51 = x h + T 3 2 sin θ 3 + L 3 cos θ 3 + L 4 cos θ 4 + T 5 2 sin θ 5 + d 51 cos θ 5 ,
z 51 = z h T 3 2 cos θ 3 + L 3 sin θ 3 L 4 sin θ 4 T 5 2 cos θ 5 + d 51 sin θ 5 + z b t
In contrast, we applied counterclockwise rotation to the trunk region through π θ 2 * , as shown in Figure 3c. The expressions for coordinates of the points, A 21 * x 21 * , z 21 * and B 21 x 21 , z 21 applied in the derivation of δ 21 are as follows:
x 21 * = x h * T 2 2 sin θ 2 * d 21 cos θ 2 * 21 * sin θ 2 * ,
z 21 * = z h * + T 2 2 cos θ 2 * + d 21 sin θ 2 * + 21 * cos θ 2 * ,
x 21 = x h T 2 2 sin θ 2 + d 21 cos θ 2 ,
z 21 = z h + T 2 2 cos θ 2 + d 21 sin θ 2 + z b t ,
The coordinates of the points, A 22 * x 22 * , z 22 * and B 22 x 22 , z 22 were obtained such that all terms in Equations (18)–(21) were retained except for the length terms d 21 and 21 * , which changed to d 22 and 22 * , respectively. Substituting the expressions in Equations (3), (4), and (7) and expressions for δ 11 , δ 12 , δ 21 , δ 22 and δ 51 , into the energy equations in (2), (5), and (6), we solve Lagrange’s equation in (1) and obtained non-linear equations of motion with respect to each generalized coordinate, q k . The expressions for the deflection terms δ 21 , δ 22 , δ 32 , and δ 51 are provided in Appendix A.

2.3. Linearization of Equations of Motion

Non-linear analysis of a vibrating system is not only complex, difficult, and time-consuming, but it also challenges access to familiar modal properties, such as natural frequencies and mode shapes. Therefore, we employed linearization techniques on the 7-DOF model’s non-linear equations of motion by expanding the Taylor’s series at the initial configurations using the following approximations:
θ = θ * + θ ~ , cos θ ~ = 1 , sin θ ~ = θ ~ ,
where θ ~ is a small perturbation around the initial value of absolute angular deflection, θ * .
The linearized equations of motion, expressed as a matrix, are written in the form:
M x ¨ t + C x ˙ t + K x t = f c z ˙ b t + f k z b t ,
where M , C , and K are mass, damping, and stiffness matrices; x t is the generalized coordinates vector; f k and f c are generalized stiffness and damping coefficients vectors; z b t is displacement due to input excitation. Taking the Fourier transform of Equation (23), the linear equations of motion in the frequency domain are expressed as:
X j ω = ω 2 M + j ω C + K 1 j ω f c + f k Z b j ω
The transmissibility from the seat to each segment of the 7-DOF model can be computed from Equation (24) and will be used in the estimation of the unknown parameters. The expressions for the matrices and force vectors are provided in Appendix B.

3. Human Vibration Experiment and Identification of Human Parameters

3.1. Experimental Set-Up and Procedures

We conducted a human vibration experiment at the Jeonbuk Institute of Automotive Technology (JIAT) to measure the dynamic response of a seated human to vertical vibrations using an MTS 248.05 uniaxial hydraulic exciter (Figure 4a). The subject sat on a viscoelastic vehicle seat (Figure 4b), fixed tightly to the vibration table using an aluminum test jig. Random vibration within the excitation frequency range of 2–80 Hz was applied to the subject in the vertical direction. The acceleration of the body segments of the subject under vertical vibrations was recorded using a NexGen monitoring system (Figure 5a), which comprises a wireless receiver, docking station, TK Motion Manager software (eZ-Analyst V5.1.140), and SXT movement monitors (inertial measurement unit, IMU sensors). During the vibration measurement, the IMU sensors were tightly fixed to the skin of the subject at the head, thigh, leg, chest (under the cloth), and at the base of the seat (Figure 5b). The IMU sensors had a sampling frequency of 128 Hz, and vertical excitations of 0.5 and 2.0 m/s2 rms were applied for 5 min per session.
Afterward, the measured acceleration signals were processed using MATLAB signal processing techniques. We applied a Hanning window with 50% overlap and 19 averages and computed the transmissibility from the base to each body segment using the ‘modalfrf’ function in MATLAB. The technical specifications of the uniaxial hydraulic exciter are listed in Table 1.

3.2. Identification of Human Parameters

For physical modeling, simulation, or parameter estimation, physical properties such as mass, mass moment of inertia, length, and thickness of body segments are required [21]. We extracted the mass, length, and thickness of human body segments from the 8th Korean Human Body Dimension Survey Report of Size Korea [22], which provided data for Korean adult males and females aged 18–59 years, with an average height of 171 cm and weight of 70 kg, within the range of body size of the subject used in our experiment. The corresponding mass moment of inertia data was reported by Ma et al. [23], who used a 3D scan to calculate the moment of inertia of various human segments. Also, we obtained the range of initial angles of the body segments in a seated posture from Kim et al. [24] and adapted the appropriate angles to the 7-DOF model. The identified human and geometric parameters used in the 7-DOF model equations of motion are listed in Table 2.

4. Optimization Methods

4.1. Objective Function Formulation

We developed an objective function to estimate the unknown stiffness and damping parameters by minimizing the sum of squared errors between measured and calculated transmissibility:
M i n i m i z e f 1 f 2 T E x p . f T M o d e l f 2 d f f 1 f 2 T E x p . f 2 d f
S u b j e c t   t o 20 k i j 200   k N / m ,   a n d   1   k i T 4   k N m / r a d
10 c i j 4000   N s / m ,   a n d   5 c i T 500   N m s / r a d
where T E x p . and T M o d e l are transmissibility obtained from experiments and those calculated using Equation (24) for each body segment, while f 1,2 is the frequency range of vibration measurements. For our analysis, T E x p . measured at the head, trunk, thigh, and leg were used as feeding data in the objective function to estimate the unknown stiffness and damping parameters. The parameters were subject to upper and lower bounds obtained from reviewing the published range of values of stiffness and damping in the literature [15,16,25,26,27]. The specified bounds limit the estimated stiffness and damping coefficients to only positive values and speed up the convergence of the optimization algorithms.

4.2. Gradient-Based Algorithm (GBA)

Gradient-based algorithms (GBAs) are commonly used in human model parameter estimation due to their speed and ease of implementation. However, they have some limitations, such as being highly sensitive to the initial values of the unknown parameters. This limitation becomes even more significant in a study like ours where there are no established data on the values of stiffness and damping coefficients that can accurately model the viscoelastic properties of the skin and seat subject to dynamic excitations.
In this study, we execute GBA using the ‘fmincon’ MATLAB optimization function to minimize the sum of squared errors between measured and calculated transmissibility in the 7-DOF model. We initially set computation parameters, including a maximum iteration of 50, a step tolerance of 1 × 10−15, and a function tolerance of 1 × 10−12, under the bounds specified in Equation (25). However, this resulted in a poor fitting of the transmissibility and high objective function value at final solutions. The objective function value at final solutions (FVAL) is a real number returned when the algorithm converges. Typically, in a minimization problem, the lower the FVAL, the higher the chances that the obtained solution is close to the optimum solution of the desired problem [28].
To improve the fitting results by GBA, we continuously adjusted the computation parameters using trial-and-error techniques. Finally, we achieved much better fits, with an FVAL of 0.0513 obtained at the final solutions of GBA using an average of 100 iterations, a maximum function evaluation of 2700, and within the computation time of 1300 s. The difficulties in finding the appropriate initial values demonstrate the challenges associated with GBA in estimating human parameters and highlight the importance of carefully selecting and adjusting the computation parameters to obtain accurate and reliable solutions.

4.3. Genetic Algorithm (GA)

Unlike gradient-based algorithms, genetic algorithm (GA) is a gradient-free optimization technique that is used to find global or near-global optima without specifying initial values. The algorithm is based on a random population evolution technique using some set of operators. In this study, we employed GA to estimate stiffness and damping values while fitting the measured and calculated transmissibility. GA was executed using the ‘ga’ function in MATLAB with a population size of 90, 9 elites, crossover probability of 0.8, function tolerance of 1 × 10−12, 100 maximum generations, adaptive feasible mutation function, and defaults values for the creation, fitness scaling, selection, and crossover functions. At this time, the algorithm obtained poor fitting of the transmissibility plots with a mean FVAL of 0.1078 returned after 880 s.
According to Mathworks [28], a larger population size improves the capability of GA to comprehensively search for optimal values close to the global optima, albeit at the expense of the computation speed. The initial population size was increased to 180 and then 270, while other computation parameters were retained, resulting in much-improved fits of the transmissibility plots as presented in Section 5.
However, the optimized parameters obtained by GA were characterized by high instability, making it difficult to perform parameter sensitivity analyses on the optimized stiffness and damping coefficients’ dynamic response to the excitation.

4.4. Hybrid Optimization Method (HOM)

The hybrid optimization method is developed to overcome the limitations of GBA and GA in the parameter estimation process (Figure 6). By combining the strengths of both algorithms, the HOM aims to improve the model fitting and reduce the randomness and exhaustive search for initial values involved in the optimization process. The first major step of HOM eliminates the exhaustive search encountered in GBA by using a typical random initial population-based GA. The computation parameters include an initial population size of 180, 18 elites, a crossover probability of 0.85, and the predefined parameters in the previous GA process. A mean FVAL of 0.0381 was returned after 1100 s in the GA process. In the second major step, GBA is used to explore the local optima of the near-global optimal solutions obtained by GA in the first step. We set the maximum iteration to 50, the maximum function evaluation to 2700, the step tolerance to 1 × 10−15, and the function tolerance to 1 × 10−12, with the optimized parameters from the first GA step as initial values.
For GBA. Interestingly, this step yielded a reduced FVAL of 0.0293 and a better fit of the transmissibility plots in the GBA process. Furthermore, we noticed that the values of some parameters obtained at different body segments varied only slightly while some remained unchanged.
In the third major step, we executed a modified GA (denoted as mGA) to address the randomness in optimized values exhibited by the typical GA, such as the one used in the first step. We defined an initial population matrix in the mGA and assigned values of the optimized parameters obtained by GBA as the elements of the population matrix. During this step, if FVAL by mGA is less than that by GBA in the previous step, a parameter sensitivity score is applied to update the bounds, and the mGA process is reinitiated. The final optimized parameters obtained in the mGA step are adapted as the HOM parameters. The HOM’s advantages over GBA and GA include improved model fitting, reduced randomness, and lower values of the objective function, making it a promising optimization method for parameter estimation in various applications.

5. Results

5.1. Experimental Results

Figure 7 shows the transmissibility obtained from processing the measured acceleration data described in the experimental section. The data is presented for four body segments: head, trunk, thigh, and leg. The results showed that the upper trunk and head segments have the highest magnitude of transmitted vibrations, both at 0.5 and 2.0 m/s2 rms. The thigh and calf on the other hand are less susceptible to the transmitted vibrations.
Furthermore, at a forcing amplitude of 0.5 m/s2 rms, all body segments exhibit a first principal resonant frequency in the range of 4–6 Hz, while the second principal reso-nance frequency is prominent in the head, thigh, and calf regions, within the range of 8–15 Hz. Our experimental approach confirms that the that the peak resonance fre-quencies of human response to vibration in the vertical direction are in the low-frequency region and the resonant frequencies obtained for each of the head, trunk, thigh and calf are in the region reported in literature [29,30,31].
Additionally, a non-linear response of the human body to an increase in excitation magnitude was observed for all the body segments. The ‘softening effect’ was observed as the excitation magnitude increased from 0.5 to 2.0 m/s2 rms, indicating the strong non-linearity of human body biodynamics. Also, the peak magnitude of transmissibility at each body segment reduced, and the resonance frequency became less observable and diminished in magnitude as the vibration exposure time increased.

5.2. Model Fitting Results Obtained by the Optimization Methods

We analyzed the performance of the optimization methods in fitting the model by comparing results obtained by GBA, GA, and HOM. The model exhibited a strong fit to experimental data for all the optimization methods, with HOM showing the overall best fit at both 0.5 and 2.0 m/s2 rms (Figure 8 and Figure 9). Furthermore, we established goodness of fit criteria, ε as suggested by Wong [32], to quantitatively evaluate the model-fitting results:
ε = 1 T E x p . T M o d e l 2 N 2 / T E x p . N ,
where N is the number of experimental data points used in the goodness of fit comparison. A value of ε close to 1 indicates a good fit, with values over 80% considered an acceptable tolerance of error for analysis. As shown in Table 3 and Table 4, an average of 94% and 95% was obtained at forcing amplitude of 0.5 and 2.0 m/s2 rms for all the optimization methods. The HOM achieved the highest ε value, signifying the best fit among the three optimization methods.
The performance of the optimization methods was also evaluated based on the objective function value at the final solutions, FVAL. Figure 10 demonstrates that when comparing the individual use of GBA or GA to the proposed method, HOM consistently achieved the lowest objective function value for all body segments. Moreover, by incorporating the GBA-optimized values in the MGA step, the computation process was expedited. The initial objective function value of approximately 0.400 rapidly decreased to a minimum of 0.0145 within the first 10 generations, as shown in Figure 11. Typically, the objective function values improve rapidly in initial generations and then level off as the optimal value is approached [28]. Detailed FVAL results for GBA, GA, and HOM for each body segment are presented in Table 5 and Table 6.
Furthermore, the computation time of HOM is significantly reduced compared to using GBA and GA separately. Specifically, the computation time of GBA and GA steps in the HOM process is over five times reduced for GBA and about 30% for GA, as shown in Table 7 and Table 8. However, the combined time of HOM is still nearly the same as that of GA and a bit more than GBA.
Overall, the optimization methods performed well in fitting the model, with HOM being the most effective. The use of HOM also produced lower objective function values and faster computation time, which significantly improved the optimization process.

5.3. Estimated Parameters and Parameter Sensitivity

The stiffness and damping coefficients estimated by GBA, GA, and HOM for 10 translational parameters: k 21 , k 22 , k 31 , k 32 and k 51 and c 21 , c 22 , c 31 , c 32 and c 51 and 8 torsional parameters: k 1 T , k 2 T , k 3 T , k 4 T and c 1 T , c 2 T , c 3 T , c 4 T at 0.5 and 2.0 m/s2 rms forcing magnitudes are presented in Figure 12 and Figure 13. The estimates were obtained by minimizing the squared error between the calculated and measured transmissibility for the head, trunk, thigh, and leg segments of the 7-DOF model. We analyzed the estimated parameters and compared the performance of each optimization method to gain insights into the key parameters that affect the appearance and magnitude of resonance frequencies and other modal characteristics of the model.
While HOM produced nearly the same values for the estimated translational and torsional parameters at either 0.5 or 2.0 m/s2 rms for all body segments, GBA and GA showed significant variances for all estimated parameters. We attribute the high randomness in estimated parameters by GA to the nature of its random initial population assignment and that of GBA to the limitations of guessing the appropriate initial values for the estimated parameters.
To improve the results of estimated parameters, we performed a parameter sensitivity analysis by establishing a parameter sensitivity matrix, P with elements [26]:
P i j = R i / R i ψ j / ψ j ,
where ψ = k 21 , k 22 , k 31 , k 32 , k 51 , c 21 , c 22 , c 31 , c 32 , c 51 , k 1 T , k 2 T , k 3 T , k 4 T , c 1 T , c 2 T , c 3 T , c 4 T ;   ψ j and ψ j represent the j th ( j = 1,2 , , 18 ) parameter at its nominal value and its value after a small perturbation; R i represent the responses to be examined, in our case, the first resonance frequency ( i = 1 ), the second resonance frequency ( i = 2 ), the peak magnitude of resonance frequency ( i = 3 ), computed with the nominal value of each term in ψ , extracted from the final optimization results that allow us to fit the model to experiment as presented in Figure 8 and Figure 9; R i represent the differences in the value of the above responses after perturbation. We set the perturbation factor, α to 30% such that ψ j = α ψ j and the new values of each term in ψ j , is ψ j n e w = ψ j ± ψ j with other parameters kept at their nominal values. Then, new values of each examined response, R i n e w are computed with ψ j n e w such that R i = R i n e w R i .
At ±30%, the resultant parameter sensitivity matrix, normalized based on the maximum value in each row, is as follows:
P + 30 % = 0.02 0.05 0.00 0.00 0.00 0.00 1.00 1.00 0.09 0.01 0.01 0.00 0.41 0.31 0.02 0.21 0.05 0.91 0.01 0.00 0.00 0.02 0.01 1.00 0.00 0.00 0.03 0.09 0.00 0.32 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.02 0.01 0.05 0.05 0.01 0.04 0.00 0.00 0.86 0.06 0.03 0.01 T   and   P 30 % = 0.05 0.08 0.00 0.00 0.00 0.00 1.00 1.00 0.05 0.01 0.00 0.00 0.52 0.49 0.01 0.28 0.06 0.85 0.02 0.00 0.00 0.02 0.01 1.00 0.00 0.00 0.04 0.05 0.00 0.38 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.01 0.01 0.03 0.02 0.01 0.02 0.00 0.00 0.78 0.08 0.05 0.01 T
For each examined response, the sensitivity of each parameter was analyzed based on their score in the parameter sensitivity matrix. The greater the value, the more sensitive the parameter is to the considered response. The values from the matrix show that k 31 , k 51 , c 21 , c 31 , c 51 and c 3 T are the most significant parameters to the appearance and magnitude of the resonance frequency. Table 9 summarizes the order of sensitivity of these parameters for the examined responses.
We utilized the sensitivity matrix score to update the bounds in the mGA repeatedly and extracted the final optimized parameters. This process allowed us to obtain numerically stable values for each of k 21 , k 22 , k 31 , k 32 and k 51 , whose estimated values remained unchanged at 0.5 m/s2 rms regardless of the body segment transmissibility data used in the objective function (Figure 12). Similarly, at 2.0 m/s2 rms, all the k i j terms remained consistent except k 31 , whose value reduced slightly for all body segments, confirming the ‘softening effect’ observed in the experimental results. The reduction in k 31 as the excitation magnitude increased, suggested the softening in the stiffness of the tissue beneath the ischial tuberosities, as reported in the literature [27,33]. As observed with the k i j terms, values of c 22 , c 32 and c 51 were consistent, while the variances in values of c 21 and c 31 for some of the body segments can be associated with their sensitivity to the magnitude of resonance frequency.
Also, as presented in Figure 13, the torsional parameters for HOM showed nearly the same values across all body segments and forcing magnitudes. Particularly, the torsional stiffness coefficients, k i T were unchanged. The changes in the torsional damping, especially c 3 T is associated with its high significance in determining the magnitude of the resonance frequencies.

6. Discussions

6.1. Proposed Model

The proposed 7-DOF seated human model was successfully validated for its intended application in predicting the response of distinct body segments to vertical vibration in the low-frequency range of 0–15 Hz. We were able to access the response of each body segment with increasing excitation magnitude, with the upper trunk and head showing the most susceptibility to transmitted vibration. Also, the model accurately captured the resonant frequencies of each body segment observed in the experimental results. Furthermore, the model demonstrated the ability to exhibit non-linearity in human biodynamics despite its simplicity by replicating the “softening effect.”
The high goodness-of-fit scores obtained during the model fitting process suggest that the proposed model is effective with different optimization algorithms. The validated model can be used in future research to optimize the comfort and safety of individuals in various seating arrangements. However, it is important to acknowledge the limitations of the proposed model, including its applicability in the high-frequency range and the limited number of experimental subjects. Further studies are necessary to address these limitations and enhance the model’s applicability in diverse scenarios.

6.2. Estimation of Unknown Parameters Using Proposed HOM

The lack of established data on human parameters, such as stiffness and damping coefficients, highlights the significance of developing an effective parameter optimization algorithm. In this study, the proposed HOM was able to combine the strengths of typical GBA and GA to estimate the unknown stiffness and damping coefficients. Not only did HOM outperform GBA and GA in terms of goodness of fit, but it also achieved the lowest objective function values at faster computation speeds.
Furthermore, the combination of the 7-DOF model and HOM revealed that the unknown human parameters, particularly stiffness, could be estimated using data measured from any of the body segments, with minor variations in optimized values across different excitation magnitudes. Also, Figure 12 demonstrates that the spring-damper at the foot, distal fragment of the thigh, and the lower trunk are stiffer than those at the upper trunk and proximal fragment of the thigh. This suggests that the vibration transmitted through the seat is greater than that through the foot. Also, the upper trunk, being less stiff, experiences more displacement, leading to more displacement of the head segment and the potential for motion sickness, as suggested by Bovenzi [2]. Similarly, the high displacement (lower stiffness) exhibited by the spring damper at the buttock region can account for the reduction in k31 at increased excitation magnitude, implying the softening of the tissue beneath the ischial tuberosities.

7. Conclusions

In this study, we successfully developed a seven-degree-of-freedom (7-DOF) model under vertical vibration and estimated 18 effective stiffness and damping coefficients using a new hybrid optimization method. The goal was to understand how vibration is transmitted to distinct human body segments, especially in the low-frequency region, by analyzing the estimated parameters by GBA, GA, and HOM using transmissibility responses from experiments and those calculated using the 7-DOF model. Results from experiments showed that the upper trunk and head are most susceptible to transmitted vibrations, with the thigh and leg showing lesser magnitudes. Also, the peak magnitude of transmissibility at each body segment reduced, and the resonance frequency became less observable and diminished in magnitude as the vibration exposure time increased.
Our findings also showed that combining the 7-DOF model with the hybrid optimization method significantly accelerated the optimization process, enhanced numerical stability, and resulted in a significant reduction of the objective function value compared to using GBA or GA as standalone algorithms. Additionally, HOM was more suitable for the estimation of the unknown stiffness and damping parameters. Notably, we found that unknown human parameters, particularly stiffness, could be estimated using data measured from any of the body segments, with minor variations in optimized values across different excitation magnitudes. Finally, the ‘softening effect’ was observed as the excitation magnitude increased from 0.5 to 2.0 m/s2 rms., indicating the strong non-linearity of human body biodynamics.
These results have practical implications for optimizing the comfort and safety of individuals exposed to whole-body vibration when riding in vehicles, as ergonomists and design engineers can use the developed model and HOM to identify the body segments most susceptible to vibration. Future work will focus on studying the correlation of transmissibility at different points on each body segment, the effects of subject variabilities, and the addition of lateral excitation.

Author Contributions

Conceptualization, A.O.A. and N.K.; methodology, A.O.A., B.-G.S., and N.K.; software, A.O.A. and B.-G.S.; validation, A.O.A. and J.-J.B.; formal analysis, A.O.A.; investigation, A.O.A.; resources, N.K.; data curation, A.O.A.; writing—original draft preparation, A.O.A.; writing—review and editing, N.K.; visualization, A.O.A.; supervision, N.K.; project administration, N.K.; funding acquisition, N.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean Government (MIST) (No. 2020R1A2C2013256).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board of Kyungpook National University (protocol code KNU-2020-0085 on 14 August 2020).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Deflection Terms and Non-linear Equations of Motion of the 7-DOF Human Body Model

The derived expression for deflection terms, δ i j are as follows:
δ 32 = z h z h * cos θ 3 * x h x h * sin θ 3 * T 3 2 cos θ 3 θ 3 * + d 32 sin θ 3 θ 3 * + T 3 2 + z b cos θ 3 *
δ 21 = x h x h * sin θ 2 * z h z h * cos θ 2 * T 2 2 cos θ 2 θ 2 * d 21 sin θ 2 θ 2 * + T 2 2 z b cos θ 2 *
δ 22 = x h x h * sin θ 2 * z h z h * cos θ 2 * T 2 2 cos θ 2 θ 2 * d 22 sin θ 2 θ 2 * + T 2 2 z b cos θ 2 *
δ 51 = z h z h * cos θ 5 * x h x h * sin θ 5 * T 5 2 cos θ 5 θ 5 * + d 51 sin θ 5 θ 5 * + T 5 2 + z b cos θ 5 * + T 3 2 cos θ 3 * θ 5 * cos θ 3 θ 5 * L 3 sin θ 3 * θ 5 * sin θ 3 θ 5 * + L 4 sin θ 4 * θ 5 * sin θ 4 θ 5 *

Appendix B. Definition of Elements of the Mass, Damping, and Stiffness Matrices and the Force Vector

The mass matrix, M is composed of:
M = M 11 0 M 13 M 14 M 15 M 16 M 17 M 22 M 23 M 24 M 25 M 26 M 27 M 33 M 34 0 0 0 M 44 M 45 0 0 S y m . M 55 M 56 M 57 M 66 M 67 M 77
and the elements of M are as follows:
M 11 = m 1 + m 2 + m 3 + m 4 + m 5 ,   M 13 = 1 2 L 1 m 1 sin θ 1 * ,
M 14 = 1 2 L 2 2 m 1 + m 2 sin θ 2 * ,   M 15 = 1 2 L 3 m 3 + 2 m 4 + 2 m 5 sin θ 3 * ,
M 16 = 1 2 L 4 m 4 + 2 m 5 sin θ 4 * ,   M 17 = 1 2 L 5 m 5 sin θ 5 *
M 22 = m 1 + m 2 + m 3 + m 4 + m 5 ,   M 23 = 1 2 L 1 m 1 cos θ 1 *
M 24 = 1 2 L 2 m 2 + 2 m 3 cos θ 2 * ,   M 25 = 1 2 L 3 m 3 + 2 m 4 + 2 m 5 cos θ 3 *
M 26 = 1 2 L 4 m 4 + 2 m 5 cos θ 4 * ,   M 27 = 1 2 L 5 m 5 cos θ 5 * ,
M 33 = J 1 + 1 4 L 1 2 m 1 ,   M 34 = 1 2 L 1 L 2 m 1 cos θ 1 * θ 2 * ,
M 44 = J 2 + 1 4 L 2 2 m 2 + 4 m 3 ,   M 45 = 1 2 L 2 L 3 m 3 cos θ 2 * θ 3 * ,
M 55 = J 3 + 1 4 L 3 2 m 3 + 4 m 4 + 4 m 5 ,   M 56 = 1 2 L 3 L 4 m 4 + 2 m 5 cos θ 3 * θ 4 * ,
M 57 = 1 2 L 3 L 5 m 5 cos θ 3 * θ 5 * ,   M 66 = J 4 + 1 4 L 4 2 m 4 + 4 m 5 ,
M 67 = 1 2 L 4 L 5 m 5 cos θ 4 * θ 5 * , M 77 = J 5 + 1 4 L 5 2 m 5
The damping matrix, C is composed of:
C = C 11 C 12 C 13 C 14 0 C 16 C 17 C 22 C 23 C 24 0 C 26 C 27 C 33 C 34 0 0 0 C 44 C 45 0 0 S y m . C 55 C 56 C 57 C 66 C 67 C 77
and the elements of C are as follows:
C 11 = c 31 + c 32 s i n 2 θ 3 * + c 21 + c 22 s i n 2 θ 2 * + c 51 s i n 2 θ 5 * ,
C 12 = c 31 + c 32 cos θ 3 * sin θ 3 * c 21 + c 22 cos θ 2 * sin θ 2 * c 51 cos θ 5 * sin θ 5 * ,
C 13 = c 31 d 31 + c 32 l 32 sin θ 3 * + T 3 2 c 51 sin θ 3 * θ 5 * + L 3 c 51 cos θ 3 * θ 5 * ,
C 15 = c 21 d 21 + c 22 d 22 sin θ 2 * ,     C 16 = c 51 L 4 sin θ 5 * cos θ 4 * + θ 5 * ,
C 17 = c 51 d 51 sin θ 5 * ,   C 22 = c 31 + c 32 c o s 2 θ 3 * + c 21 + c 22 c o s 2 θ 2 * + c 51 c o s 2 θ 5 * ,
C 23 = c 31 d 31 + c 32 d 32 cos θ 3 * + c 51 L 3 cos θ 3 * θ 5 * + T 3 2 sin θ 3 * θ 5 * cos θ 5 * ,
C 24 = c 21 d 21 + c 22 d 22 cos θ 2 * ,   C 26 = c 51 L 4 cos θ 5 * cos θ 4 * + θ 5 * ,   C 27 = c 51 d 51 cos θ 5 * ,
C 33 = c 1 T ,   C 34 = c 1 T ,   C 44 = c 21 d 21 2 + c 22 d 22 2 + c 1 T + c 2 T ,
C 45 = c 1 T ,   C 55 = c 31 d 31 2 + c 32 d 32 2 + c 51 L 3 cos θ 3 * θ 5 * + T 3 2 sin θ 3 * θ 5 * 2 + c 2 T + c 3 T ,
C 56 = c 51 L 3 L 4 cos θ 4 * θ 5 * cos θ 3 * θ 5 * + T 3 2 sin θ 3 * θ 5 * c 3 T ,
C 57 = c 51 d 51 L 3 cos θ 3 * θ 5 * + T 3 2 sin θ 3 * θ 5 * , C 66 = c 51 L 4 2 c o s 2 θ 4 * θ 5 * + c 3 T + c 4 T ,
C 67 = c 51 d 51 L 4 cos θ 4 * θ 5 * c 4 T ,   C 77 = c 51 d 51 2 + c 4 T
The stiffness matrix, K is composed of:
K = K 11 K 12 K 13 K 14 0 K 16 K 17 K 22 K 23 K 24 0 K 26 K 27 K 33 K 34 0 0 0 K 44 K 45 0 0 S y m . K 55 K 56 K 57 K 66 K 67 K 77
and the elements of K are as follows:
K 11 = k 31 + k 32 s i n 2 θ 3 * + k 21 + k 22 s i n 2 θ 2 * + K 51 s i n 2 θ 5 * ,
K 12 = k 31 + k 32 cos θ 3 * sin θ 3 * k 21 + k 22 cos θ 2 * sin θ 2 * k 51 cos θ 5 * sin θ 5 * ,
K 13 = k 31 d 31 + k 32 d 32 sin θ 3 * + T 3 2 k 51 sin θ 3 * θ 5 * + L 3 k 51 cos θ 3 * θ 5 * ,
K 15 = k 21 d 21 + k 22 d 22 sin θ 2 * ,   K 16 = k 51 L 4 sin θ 5 * cos θ 4 * + θ 5 * ,
K 17 = k 51 d 51 sin θ 5 * ,   K 22 = k 31 + k 32 c o s 2 θ 3 * + k 21 + k 22 c o s 2 θ 2 * + k 51 c o s 2 θ 5 * ,
K 23 = k 31 d 31 + k 32 d 32 cos θ 3 * + k 51 L 3 cos θ 3 * θ 5 * + T 3 2 sin θ 3 * θ 5 * cos θ 5 * ,
K 24 = k 21 d 21 + k 22 d 22 cos θ 2 * ,   K 26 = k 51 L 4 cos θ 5 * cos θ 4 * + θ 5 * ,   K 27 = k 51 d 51 cos θ 5 * ,
K 33 = k 1 T ,   K 34 = k 1 T ,   K 44 = k 21 d 21 2 + k 22 d 22 2 + k 1 T + k 2 T ,
K 45 = k 1 T ,   K 55 = k 31 d 31 2 + k 32 d 32 2 + k 51 L 3 cos θ 3 * θ 5 * + T 3 2 sin θ 3 * θ 5 * 2 + k 2 T + k 3 T ,
K 56 = k 51 L 3 L 4 cos θ 4 * θ 5 * cos θ 3 * θ 5 * + T 3 2 sin θ 3 * θ 5 * k 3 T ,
K 57 = k 51 d 51 L 3 cos θ 3 * θ 5 * + T 3 2 sin θ 3 * θ 5 * , K 66 = k 51 L 4 2 c o s 2 θ 4 * θ 5 * + k 3 T + k 4 T ,
K 67 = k 51 d 51 L 4 cos θ 4 * θ 5 * k 4 T ,   K 77 = k 51 d 51 2 + k 4 T
The force vector, F e = f k + f c is composed of:
F e = F e 1 F e 2 F e 3 F e 4 F e 5 F e 6 F e 7 T
and the elements of the force vector, F e are as follows:
F e 1 = c 31 + c 32 cos θ 3 * sin θ 3 * + c 21 + c 22 cos θ 2 * sin θ 2 * + c 51 cos θ 5 * sin θ 5 * z ˙ b   + k 31 + k 32 cos θ 3 * sin θ 3 * + k 21 + k 22 cos θ 2 * sin θ 2 * + k 51 cos θ 5 * sin θ 5 * z b
F e 2 = c 31 + c 32 cos 2 θ 3 * + c 21 + c 22 cos 2 θ 2 * + c 51 cos 2 θ 5 * z ˙ b   k 31 + k 32 cos 2 θ 3 * + k 21 + k 22 cos 2 θ 2 * + k 51 cos 2 θ 5 * z b M g
F e 3 = 1 2 L 1 m 1 g cos θ 1 *
F e 4 = c 21 d 21 + c 22 d 22 cos θ 2 * z ˙ b k 21 d 21 + k 22 d 22 cos θ 2 * z b   1 2 L 2 m 3 + 2 m 1 g cos θ 2 *
F e 5 = c 31 d 31 + c 32 d 32 cos θ 3 * + c 51 cos θ 5 * L 3 cos θ 3 * θ 5 * + T 3 2 sin θ 3 * θ 5 * z ˙ b   k 31 d 31 + k 32 d 32 cos θ 3 * + k 51 cos θ 5 * L 3 cos θ 3 * θ 5 * + T 3 2 sin θ 3 * θ 5 * z b   1 2 L 3 m 3 + 2 m 4 + 2 m 5 g cos θ 3 *
F e 6 = c 51 L 4 cos θ 5 * cos θ 4 * θ 5 * z ˙ b k 51 L 4 cos θ 5 * cos θ 4 * θ 5 * z b   1 2 L 4 m 4 + 2 m 5 g cos θ 4 *
F e 7 = c 51 d 51 cos θ 5 * z ˙ b k 51 d 51 cos θ 5 * z b 1 2 L 5 m 5 g cos θ 5 *

References

  1. Bovenzi, M.; Betta, A. Low-back disorders in agricultural tractor drivers exposed to whole-body vibration and postural stress. Appl. Ergon. 1994, 25, 231–241. [Google Scholar] [CrossRef] [PubMed]
  2. Bovenzi, M.; Hulshof, C.T.J. An Updated Review of Epidemiologic Studies on the Relationship Between Exposure to Whole-Body Vibration and Low Back Pain (1986–1997). Int. Arch. Occup. Environ. Health 1999, 72, 351–365. [Google Scholar] [CrossRef] [PubMed]
  3. International Organization for Standardization (ISO). Mechanical Vibration and Shock: Range of Idealized Values to Characterize Seated-Body Biodynamic Response under Vertical Vibration; International Organization for Standardization: Geneva, Switzerland, 2002. [Google Scholar]
  4. Zhou, Z.; Griffin, M.J. Response of the Seated Human Body to Whole-Body Vertical Vibration: Discomfort Caused by Sinusoidal Vibration. Ergonomics 2014, 57, 714–732. [Google Scholar] [CrossRef]
  5. Kim, E.; Fard, M.; Kato, K. Characterisation of the Human-seat Coupling in Response to Vibration. Ergonomics 2017, 60, 1085–1100. [Google Scholar] [CrossRef] [PubMed]
  6. Rakheja, S.; Dewangan, K.N.; Dong, R.G.; Marcotte, P.; Pranesh, A. Whole-Body Vibration Biodynamics—A Critical Review: II. Biodynamic Modelling. Int. J. Veh. Perform. 2020, 6, 52–84. [Google Scholar] [CrossRef]
  7. Matsumoto, Y.; Griffin, M.J. Non-linear Characteristics in the Dynamic Responses of Seated Subjects Exposed to Vertical Whole-Body Vibration. J. Biomech. Eng. 2002, 124, 527–532. [Google Scholar] [CrossRef]
  8. Cho, Y.; Yoon, Y.S. Biomechanical Model of Human on Seat With Backrest for Evaluating Ride Quality. Int. J. Ind. Ergon. 2001, 27, 331–345. [Google Scholar] [CrossRef]
  9. Kim, T.H.; Kim, Y.T.; Yoon, Y.S. Development of a biomechanical model of the human body in a sitting posture with vibration transmissibility in the vertical direction. Int. J. Ind. Ergon. 2005, 35, 817–829. [Google Scholar] [CrossRef]
  10. Bae, J.J.; Kang, N. Development of a Five-Degree-of-Freedom Seated Human Model and Parametric Studies for its Vibrational Characteristics. Shock Vib. 2018, 2018, 1649180. [Google Scholar] [CrossRef]
  11. Tufano, S.; Griffin, M.J. Nonlinearity in the vertical transmissibility of seating: The role of the human body apparent mass and seat dynamic stiffness. Veh. Syst. Dyn. 2013, 51, 122–138. [Google Scholar] [CrossRef]
  12. Dewangan, K.N.; Rakheja, S.; Marcotte, P.; Shahmir, A. Effects of elastic seats on seated body apparent mass responses to vertical whole-body vibration. Ergonomics 2015, 58, 1175–1190. [Google Scholar] [CrossRef] [PubMed]
  13. Wang, W.; Rakheja, S.; Boileau, P.É. Effect of back support condition on seat to head transmissibilities of seated occupants under vertical vibration. J. Low Freq. Noise Vib. Act. Control 2006, 25, 239–259. [Google Scholar] [CrossRef]
  14. Matsumoto, Y.; Griffin, M.J. Modelling the dynamic mechanisms associated with the principal resonance of the seated human body. Clin. Biomech. 2001, 16, S31–S44. [Google Scholar] [CrossRef] [PubMed]
  15. Shao, X.; Xu, N.; Liu, X. Investigation and Application on the Vertical Vibration Models of the Seated Human Body. Automot. Innov. 2018, 1, 263–271. [Google Scholar] [CrossRef]
  16. Huang, Y.; Zhang, P.; Liang, S. Apparent mass of the seated human body during vertical vibration in the frequency range 2–100 Hz. Ergonomics 2020, 63, 1150–1163. [Google Scholar] [CrossRef]
  17. Zhang, S.; Shi, W.; Chen, Z. Modeling and Parameter Identification of Seated Human Body with the Reference Vector Guided Evolutionary Algorithm. Adv. Mech. Eng. 2021, 13, 16878140211062679. [Google Scholar] [CrossRef]
  18. 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]
  19. Yazdi, S.J.M.; Cho, K.S.; Kang, N. Characterization of the Viscoelastic Model of In Vivo Human Posterior Thigh Skin using Ramp-Relaxation Indentation Test. Korea-Aust. Rheol. J. 2018, 30, 293–307. [Google Scholar] [CrossRef]
  20. Zhang, X.; Yu, P.; Li, Y.; Qiu, Y.; Sun, C.; Wang, Z.; Liu, C. Dynamic Interaction Between the Human Body and the Seat During Vertical Vibration: Effect of Inclination of the Seat Pan and the Backrest on Seat Transmissibilities. Ergonomics 2022, 65, 691–703. [Google Scholar] [CrossRef]
  21. Xiu, W.; Zhang, L.; Ma, O. Experimental study of a momentum-based method for identifying the inertia barycentric parameters of a human body. Multibody Syst. Dyn. 2016, 36, 237–255. [Google Scholar] [CrossRef]
  22. Size Korea. The 8th Korean Human Dimensions Survey Project. Available online: https://sizekorea.kr/human-info/meas-report?measDegree=8 (accessed on 15 March 2021).
  23. Ma, Y.; Kwon, J.; Mao, Z.; Lee, K.; Li, L.; Chung, H. Segment Inertial Parameters of Korean Adults Estimated from Three-Dimensional Body Laser Scan Data. Int. J. Ind. Ergon. 2011, 41, 19–29. [Google Scholar] [CrossRef]
  24. Kim, K.S.; Kim, J.; Kim, K.J. Dynamic Modeling of Seated Human Body Based on Measurements of Apparent Inertia Matrix for Fore-and-Aft/Vertical/Pitch Motion. J. Sound Vib. 2011, 330, 5716–5735. [Google Scholar] [CrossRef]
  25. Liang, C.C.; Chiang, C.F. Modeling of a Seated Human Body Exposed to Vertical Vibrations in Various Automotive Postures. Ind. Health 2008, 46, 125–137. [Google Scholar] [CrossRef] [PubMed]
  26. Qiu, Y.; Griffin, M.J. Modelling the fore-and-aft apparent mass of the human body and the transmissibility of seat backrests. Veh. Syst. Dyn. 2011, 49, 703–722. [Google Scholar] [CrossRef]
  27. Wu, J.; Qiu, Y.; Sun, C. Modelling and analysis of coupled vibration of human body in the sagittal and coronal planes exposed to vertical, lateral and roll vibrations and the comparison with modal test. Mech. Syst. Signal Proc. 2022, 166, 108439. [Google Scholar] [CrossRef]
  28. Math Works, Inc. Genetic Algorithm and Direct Search Toolbox for Use with MATLAB: User’s Guide. Available online: http://cda.psych.uiuc.edu/matlab_pdf/gads_tb.pdf (accessed on 2 May 2023).
  29. Fairley, T.E.; Griffin, M.J. The apparent mass of the seated human body: Vertical vibration. J. Biomech. 1989, 22, 81–94. [Google Scholar] [CrossRef]
  30. Nawayseh, N.; Griffin, M.J. Non-linear dual-axis biodynamic response to vertical whole-body vibration. J. Sound Vibr. 2003, 268, 503–523. [Google Scholar] [CrossRef]
  31. Wu, J.; Qiu, Y. Modelling of Seated Human Body Exposed to Combined Vertical, Lateral and Roll Vibrations. J. Sound Vibr. 2020, 485, 115509. [Google Scholar] [CrossRef]
  32. Wong, J.Y. Terramechanics and Off-Road Vehicle Engineering: Terrain Behaviour, Off-Road Vehicle Performance and Design; Butterworth-Heinemann: Oxford, UK, 2009. [Google Scholar]
  33. Mansfield, N.J.; Griffin, M.J. Non-linearities in apparent mass and transmissibility during exposure to whole-body vertical vibration. J. Biomech. 2000, 33, 933–941. [Google Scholar] [CrossRef]
Figure 1. The proposed 7-DOF seated human model.
Figure 1. The proposed 7-DOF seated human model.
Applsci 13 10065 g001
Figure 2. The deformed and undeformed configurations of the 7-DOF model during vibrations.
Figure 2. The deformed and undeformed configurations of the 7-DOF model during vibrations.
Applsci 13 10065 g002
Figure 3. Detailed schematic of deflection of the spring damper.
Figure 3. Detailed schematic of deflection of the spring damper.
Applsci 13 10065 g003
Figure 4. Equipment used for vibration excitation and the mounted vehicle seat.
Figure 4. Equipment used for vibration excitation and the mounted vehicle seat.
Applsci 13 10065 g004
Figure 5. The NexGen monitoring system and locations of IMU sensors on the subject.
Figure 5. The NexGen monitoring system and locations of IMU sensors on the subject.
Applsci 13 10065 g005
Figure 6. Hybrid optimization method (HOM) flowchart.
Figure 6. Hybrid optimization method (HOM) flowchart.
Applsci 13 10065 g006
Figure 7. Transmissibility obtained from experiments at 0.5 and 2.0 m/s2 rms forcing magnitudes.
Figure 7. Transmissibility obtained from experiments at 0.5 and 2.0 m/s2 rms forcing magnitudes.
Applsci 13 10065 g007
Figure 8. Model fitting results obtained by GBA, GA, and HOM at 0.5 m/s2 rms forcing amplitude.
Figure 8. Model fitting results obtained by GBA, GA, and HOM at 0.5 m/s2 rms forcing amplitude.
Applsci 13 10065 g008
Figure 9. Model fitting results obtained by GBA, GA, and HOM at 2.0 m/s2 rms forcing amplitude.
Figure 9. Model fitting results obtained by GBA, GA, and HOM at 2.0 m/s2 rms forcing amplitude.
Applsci 13 10065 g009
Figure 10. Comparison of the objective function values by optimization methods.
Figure 10. Comparison of the objective function values by optimization methods.
Applsci 13 10065 g010
Figure 11. Target and mean value of FVAL by mGA with respect to the generations of evolution in the HOM process.
Figure 11. Target and mean value of FVAL by mGA with respect to the generations of evolution in the HOM process.
Applsci 13 10065 g011
Figure 12. Translational parameters (top: stiffness coefficients, bottom: damping coefficients).
Figure 12. Translational parameters (top: stiffness coefficients, bottom: damping coefficients).
Applsci 13 10065 g012
Figure 13. Torsional parameters (top: stiffness coefficients, bottom: damping coefficients).
Figure 13. Torsional parameters (top: stiffness coefficients, bottom: damping coefficients).
Applsci 13 10065 g013
Table 1. Specifications of the uniaxial exciter used in the experiment.
Table 1. Specifications of the uniaxial exciter used in the experiment.
ModelMax.
Dynamic Force
Max.
Displacement
Operating
Frequency
Vibration
Table
MTS 248.05
Hydraulic
±50 kN75 mm0.1~100 Hz1.2 × 1.2 m
Table 2. Identified human and geometric parameters.
Table 2. Identified human and geometric parameters.
Mass
(kg)
Moment of
Inertia (kg·m2)
Length
(m)
Thickness (m)Initial Angle
(deg.)
Spring Ends Distance (m)
m 1 5.66 J 1 0.037 L 1 0.263 θ 1 100 d 1 0.100
m 2 37.96 J 2 2.050 L 2 0.599 T 2 0.224 θ 2 111 d 2 0.479
m 3 18.74 J 3 0.324 L 3 0.495 T 3 0.156 θ 3 12 d 3 0.088
m 4 6.08 J 4 0.061 L 4 0.368 θ 4 306 d 4 0.459
m 5 1.56 J 5 0.005 L 5 0.248 T 5 0.076 θ 5 10 d 5 0.124
Table 3. Goodness of fit at 0.5 m/s2 rms.
Table 3. Goodness of fit at 0.5 m/s2 rms.
MethodHead (%)Trunk (%)Thigh (%)Leg (%)Average (%)
GBA95.892.594.682.891.4
GA98.497.390.983.992.6
HOM99.098.997.991.496.8
Table 4. Goodness of fit at 2.0 m/s2 rms.
Table 4. Goodness of fit at 2.0 m/s2 rms.
MethodHead (%)Trunk (%)Thigh (%)Leg (%)Average (%)
GBA90.591.490.792.991.4
GA99.999.593.494.296.7
HOM99.598.998.795.998.2
Table 5. FVAL obtained when GBA and GA were used as standalone algorithms.
Table 5. FVAL obtained when GBA and GA were used as standalone algorithms.
MethodObjective Function Value, FVAL
HeadTrunkThighLeg
GBA alone0.04780.03930.07190.0448
GA alone0.06270.05130.05930.0571
Table 6. FVAL obtained for the combined HOM steps.
Table 6. FVAL obtained for the combined HOM steps.
HOM StepsObjective Function Value, FVAL
HeadTrunkThighLeg
GA step0.03280.03330.03940.0348
GBA step0.02270.02930.02630.0301
MGA step0.01450.01750.01450.0181
Table 7. Computation time for each body segment when using GBA and GA alone.
Table 7. Computation time for each body segment when using GBA and GA alone.
MethodComputation Time (S)
HeadTrunkThighLeg
GBA alone1510162016101760
GA alone1310129013101390
Table 8. Computation time for each body segment for combined HOM steps.
Table 8. Computation time for each body segment for combined HOM steps.
HOM StepsComputation Time (S)
HeadTrunkThighLeg
GA step1066109210661100
GBA step210200206260
MGA step250240210290
Total1526153214821650
Table 9. Summary of sensitivity of key parameters to the appearance and magnitude of the resonance frequency at ±30% perturbation.
Table 9. Summary of sensitivity of key parameters to the appearance and magnitude of the resonance frequency at ±30% perturbation.
Examined ResponsesDecreasing Order of Sensitivity
First resonance frequency k 31 ,       k 51 ,       c 21
Second resonance frequency k 31 ,       k 51
Peak magnitude of resonance frequency c 31 ,       c 21 ,       c 3 T ,       c 51
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

Alabi, A.O.; Song, B.-G.; Bae, J.-J.; Kang, N. Development of a 7-DOF Biodynamic Model for a Seated Human and a Hybrid Optimization Method for Estimating Human-Seat Interaction Parameters. Appl. Sci. 2023, 13, 10065. https://doi.org/10.3390/app131810065

AMA Style

Alabi AO, Song B-G, Bae J-J, Kang N. Development of a 7-DOF Biodynamic Model for a Seated Human and a Hybrid Optimization Method for Estimating Human-Seat Interaction Parameters. Applied Sciences. 2023; 13(18):10065. https://doi.org/10.3390/app131810065

Chicago/Turabian Style

Alabi, Abeeb Opeyemi, Byoung-Gyu Song, Jong-Jin Bae, and Namcheol Kang. 2023. "Development of a 7-DOF Biodynamic Model for a Seated Human and a Hybrid Optimization Method for Estimating Human-Seat Interaction Parameters" Applied Sciences 13, no. 18: 10065. https://doi.org/10.3390/app131810065

APA Style

Alabi, A. O., Song, B. -G., Bae, J. -J., & Kang, N. (2023). Development of a 7-DOF Biodynamic Model for a Seated Human and a Hybrid Optimization Method for Estimating Human-Seat Interaction Parameters. Applied Sciences, 13(18), 10065. https://doi.org/10.3390/app131810065

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