Hole Matrix Mapping Model for Partitioned Sitting Surface Based on Human Body Pressure Distribution Matrix

(1) Objective: The objective of this study was to experimentally obtain the ideal pressure distribution model of buttock and thigh support for office workers in forward-leaning and upright sitting postures, reproduce the support provided by mesh materials with elastic materials, and propose an effective seat design scheme to improve the comfort of office workers. (2) Method: Based on the seven most popular mesh chairs on the market, pressure distribution experiments, and the fuzzy clustering algorithm, the relatively ideal body pressure distribution matrices were generated for office workers under two common sitting postures, and the corresponding partitioned sitting support surfaces were obtained. A prototype chair was created and validated by combining the ergonomics node coordinates and the physical properties of the materials. (3) Result: An ideal support model of four zones was constructed, and prototype pads were designed and produced according to this model. Subjects were recruited to test the ability of the prototypes to reproduce the ideal pressure distribution maps. (4) Conclusion: The four-zone ideal support model is capable of effectively representing the buttock and thigh support requirements in forward-leaning and upright sitting postures, and it is useful for the development of related products. Studying sitting postures and pressure values generated by different activities of office workers will help to refine the needs of office personnel and provide new ideas for the design of office chairs.


Introduction
In urban environments, heavy workloads force office workers to work at their desks continuously for long periods of time [1]. Over time, this will lead to muscle soreness in the waist, buttocks, and back, and even the deformation of the entire spine's physiological curvature [2][3][4][5]. Sitting for a long time has become increasingly recognized as the culprit behind chronic diseases for many office workers [6]. The physical damage caused by sitting for a long time has become a major public health risk. A good seat surface can provide reasonable support for the buttocks and thighs, helping ensure good spinal curvature in office sitting postures and reducing the harm to the human body caused by sitting for a long time. Therefore, research on the comfort of the seating surface is of great importance.
There are more related studies for seat comfort [1,[7][8][9][10], which provide significant theoretical bases for seat designs. The literature suggests that the support conditions for various contact areas between the seating surface and a human body should be varied for optimal comfort [11][12][13]. Vink and Lips found that the sensitivities of the buttocks and thighs are not the same, as the buttocks are less sensitive than the front of the thigh and the middle of the thigh, and the front of the thigh is more sensitive than the buttocks and the middle of the thigh [13]. Therefore, the support conditions of various parts and regions where the seat cushion and the human body come into contact should be different to optimize comfort.
In addition, the study of different sitting type states and their body pressure distribution is important for the study of seat comfort. For the definition of sitting posture, (memory foam) was mapped out, and the prototype was made. The specific process is shown in Figure 1. (3). Finally, a back-testing experiment is conducted on the prototype model to observe the sitting comfort of different groups so as to verify the relevant knowledge extracted from the model.

Experimental Apparatus and Parameter
In this experiment, the pressure distribution system made by Tactilus, Hangzhou City, Zhejiang Province, was used to collect the pressure distribution information of the subject's sitting surface. The pressure sensor pad consisted of 32 rows and 32 columns, for a total of 1024 sensors. The pressure sensor software could calculate and return objective indicators such as maximum pressure, average pressure, and contact area. The pressure distribution parameters in the experimental analysis were the maximum pressure, the average pressure, the maximum pressure gradient, and the average pressure gradient.
The experimental apparatus used mesh chairs. After comparing 40 mesh chairs available on the market, we select 7 mesh chairs with the highest user ratings as the samples. All 7 mesh chairs could adjust sitting height, sitting depth, sitting surface inclination, lumbar height, lumbar angle, backrest depth, backrest height, backrest angle, and pivot point height and depth. They also allowed a transition between an upright sitting posture and a forward-leaning sitting posture. Their shapes and basic parameters are shown in Table  1.
To eliminate the interference of muscle activation, it was necessary to ensure that the subjects' muscles were completely relaxed, so surface electromyography (sEMG) was used as an observable.  In this experiment, the pressure distribution system made by Tactilus, Hangzhou City, Zhejiang Province, was used to collect the pressure distribution information of the subject's sitting surface. The pressure sensor pad consisted of 32 rows and 32 columns, for a total of 1024 sensors. The pressure sensor software could calculate and return objective indicators such as maximum pressure, average pressure, and contact area. The pressure distribution parameters in the experimental analysis were the maximum pressure, the average pressure, the maximum pressure gradient, and the average pressure gradient.

Experiment on Human Body Pressure Distribution of the Buttocks and Thighs
The experimental apparatus used mesh chairs. After comparing 40 mesh chairs available on the market, we select 7 mesh chairs with the highest user ratings as the samples. All 7 mesh chairs could adjust sitting height, sitting depth, sitting surface inclination, lumbar height, lumbar angle, backrest depth, backrest height, backrest angle, and pivot point height and depth. They also allowed a transition between an upright sitting posture and a forward-leaning sitting posture. Their shapes and basic parameters are shown in Table 1.  According to previous research [42], the main group of people using office chairs in China has the following characteristics: the overall level of education is relatively high, those involved in design are the class of people using office chairs more often, the main age group is 18-25 years old, and a higher percentage of women. This cohort made the ideal design target group.
A total of 7 healthy design graduate students, 4 males and 3 females, were recruited for the experiment. To reduce data variation caused by different body shapes, the subjects' basic morphological parameters were as consistent as possible. The basic subject information is as follows: the male subjects were 24 ± 2 years old, 173.1 ± 3.3 cm tall, 70 ± 4.2 kg in weight, and the female subjects were 24 ± 2 years old, 162.1 ± 2.2. 5 cm tall, and 52.4 ±  According to previous research [42], the main group of people using office chairs in China has the following characteristics: the overall level of education is relatively high, those involved in design are the class of people using office chairs more often, the main age group is 18-25 years old, and a higher percentage of women. This cohort made the ideal design target group.
A total of 7 healthy design graduate students, 4 males and 3 females, were recruited for the experiment. To reduce data variation caused by different body shapes, the subjects' basic morphological parameters were as consistent as possible. The basic subject information is as follows: the male subjects were 24 ± 2 years old, 173.1 ± 3.3 cm tall, 70 ± 4.2 kg in weight, and the female subjects were 24 ± 2 years old, 162.1 ± 2.2. 5 cm tall, and 52.4 ±  According to previous research [42], the main group of people using office chairs in China has the following characteristics: the overall level of education is relatively high, those involved in design are the class of people using office chairs more often, the main age group is 18-25 years old, and a higher percentage of women. This cohort made the ideal design target group.
A total of 7 healthy design graduate students, 4 males and 3 females, were recruited for the experiment. To reduce data variation caused by different body shapes, the subjects' basic morphological parameters were as consistent as possible. The basic subject information is as follows: the male subjects were 24 ± 2 years old, 173.1 ± 3.3 cm tall, 70 ± 4.2 kg in weight, and the female subjects were 24 ± 2 years old, 162.1 ± 2.2. 5 cm tall, and 52.4 ±  According to previous research [42], the main group of people using office chairs in China has the following characteristics: the overall level of education is relatively high, those involved in design are the class of people using office chairs more often, the main age group is 18-25 years old, and a higher percentage of women. This cohort made the ideal design target group.
A total of 7 healthy design graduate students, 4 males and 3 females, were recruited for the experiment. To reduce data variation caused by different body shapes, the subjects' basic morphological parameters were as consistent as possible. The basic subject information is as follows: the male subjects were 24 ± 2 years old, 173.1 ± 3.3 cm tall, 70 ± 4.2 kg in weight, and the female subjects were 24 ± 2 years old, 162.1 ± 2.2. 5 cm tall, and 52.4 ±  According to previous research [42], the main group of people using office chairs in China has the following characteristics: the overall level of education is relatively high, those involved in design are the class of people using office chairs more often, the main age group is 18-25 years old, and a higher percentage of women. This cohort made the ideal design target group.
A total of 7 healthy design graduate students, 4 males and 3 females, were recruited for the experiment. To reduce data variation caused by different body shapes, the subjects' basic morphological parameters were as consistent as possible. The basic subject information is as follows: the male subjects were 24 ± 2 years old, 173.1 ± 3.3 cm tall, 70 ± 4.2 kg in weight, and the female subjects were 24 ± 2 years old, 162.1 ± 2.2. 5 cm tall, and 52.4 ±  According to previous research [42], the main group of people using office chairs in China has the following characteristics: the overall level of education is relatively high, those involved in design are the class of people using office chairs more often, the main age group is 18-25 years old, and a higher percentage of women. This cohort made the ideal design target group.
A total of 7 healthy design graduate students, 4 males and 3 females, were recruited for the experiment. To reduce data variation caused by different body shapes, the subjects' basic morphological parameters were as consistent as possible. The basic subject information is as follows: the male subjects were 24 ± 2 years old, 173.1 ± 3.3 cm tall, 70 ± 4.2 kg in weight, and the female subjects were 24 ± 2 years old, 162.1 ± 2.2. 5 cm tall, and 52.4 ±  According to previous research [42], the main group of people using office chairs in China has the following characteristics: the overall level of education is relatively high, those involved in design are the class of people using office chairs more often, the main age group is 18-25 years old, and a higher percentage of women. This cohort made the ideal design target group.
A total of 7 healthy design graduate students, 4 males and 3 females, were recruited for the experiment. To reduce data variation caused by different body shapes, the subjects' basic morphological parameters were as consistent as possible. The basic subject information is as follows: the male subjects were 24 ± 2 years old, 173.1 ± 3.3 cm tall, 70 ± 4.2 kg in weight, and the female subjects were 24 ± 2 years old, 162. To eliminate the interference of muscle activation, it was necessary to ensure that the subjects' muscles were completely relaxed, so surface electromyography (sEMG) was used as an observable.

Experimental Subjects
According to previous research [41], the main group of people using office chairs in China has the following characteristics: the overall level of education is relatively high, those involved in design are the class of people using office chairs more often, the main age group is 18-25 years old, and a higher percentage of women. This cohort made the ideal design target group.
A total of 7 healthy design graduate students, 4 males and 3 females, were recruited for the experiment. To reduce data variation caused by different body shapes, the subjects' basic morphological parameters were as consistent as possible. The basic subject information is as follows: the male subjects were 24 ± 2 years old, 173.1 ± 3.3 cm tall, 70 ± 4.2 kg in weight, and the female subjects were 24 ± 2 years old, 162.1 ± 2.2. 5 cm tall, and 52.4 ± 3.1 kg in weight.

Experimental Procedure
The test included two parts: a subjective evaluation test and an objective measurement. As described in Section 1, the experiment took two sitting postures, forward-leaning and upright, and the sitting posture holding time was specified with reference to the findings of Sun Xinxin et al. [21]. The average time to maintain a forward-leaning sitting posture was 11.6 min, after which the torso would alternately lean forward and backward to relieve body fatigue, and the average holding time for upright sitting posture was only 2.8 min. Considering the limited holding time of the two sitting postures, in the experiment, subjects were required to maintain two sitting postures for 2 min for upright sitting and 10 min for forward-leaning sitting, respectively, while the experimental sitting posture criteria were determined as follows: upright sitting posture was determined as the user's subjective perceived ideal sitting posture. The forward-leaning sitting posture was a forward tilt of the trunk centerline of less than 20 • [42].
Before starting the test, subjects were allowed to familiarize themselves with the seat for at least 2 min. When collecting body pressure distribution data, subjects were allowed to take a 2-min break in the process of position switching to reduce fatigue. Between tests with different seat types, subjects were allowed to stand up and move slightly for 5 min.
During the test, to ensure that the subjects' muscles were completely relaxed, EMG signals from the transverse erector spinae (located on both sides of the spine in the waist and in the pelvis) and the gluteus maximus (located on the posterolateral side of the pelvis) that are most closely related to the posture of the buttocks and thighs were monitored as observation data [43],. As shown in Table 2 and Figure 2. When the EMG showed electrical silence, the muscles were considered to be fully relaxed, so the pressure distribution data was valid [44]. Table 2. Instructions for the Location of the Muscle to be Tested and the Electrode Sheet to be Pasted.

Electrode location
Healthcare 2023, 11, x 5 of 20 At the end of the test, the subjects were asked to fill out an experience assessment form. The assessment was done with reference to six indicators mentioned in the paper: shoulder, back, lower back, hip, thigh and overall [11]. Each item was scored on a 5-point scale, with 1 being extremely uncomfortable and 5 being extremely comfortable, up to a maximum score of 30.

Name of the Target Muscle Psoas Gluteus Maximus
Electrode location

Standardization
The results of the study by kazushige et al. [46] showed that the cushion area could be divided as shown in Figure 3, with zone A being the 4 cm × 4 cm area under the sciatic bone, zone B being the 10 cm × 10 cm area around zone A, zone C being the entire buttock area, and zone D being the area outside zone C and including the leg area, and based on the results of this study, the size of the sciatic tuberosity area was obtained. This was to At the end of the test, the subjects were asked to fill out an experience assessment form. The assessment was done with reference to six indicators mentioned in the paper: shoulder, back, lower back, hip, thigh and overall [11]. Each item was scored on a 5-point scale, with 1 being extremely uncomfortable and 5 being extremely comfortable, up to a maximum score of 30.

Name of the Target Muscle Psoas Gluteus Maximus
Electrode location

Standardization
The results of the study by kazushige et al. [46] showed that the cushion area could be divided as shown in Figure 3, with zone A being the 4 cm × 4 cm area under the sciatic bone, zone B being the 10 cm × 10 cm area around zone A, zone C being the entire buttock area, and zone D being the area outside zone C and including the leg area, and based on the results of this study, the size of the sciatic tuberosity area was obtained. This was to At the end of the test, the subjects were asked to fill out an experience assessment form. The assessment was done with reference to six indicators mentioned in the paper: shoulder, back, lower back, hip, thigh and overall [11]. Each item was scored on a 5-point scale, with 1 being extremely uncomfortable and 5 being extremely comfortable, up to a maximum score of 30.

Standardization
The results of the study by kazushige et al. [46] showed that the cushion area could be divided as shown in Figure 3, with zone A being the 4 cm × 4 cm area under the sciatic bone, zone B being the 10 cm × 10 cm area around zone A, zone C being the entire buttock area, and zone D being the area outside zone C and including the leg area, and based on the results of this study, the size of the sciatic tuberosity area was obtained. This was to facilitate the acquisition and processing of body pressure distribution data, to predict the position of the sciatic tuberosity, to ensure that the sciatic tuberosity position was basically the same when the subjects were in the experiment, and to ensure that the data was not shifted when the body pressure distribution matrix was subsequently averaged to ensure accuracy.

Standardization
The results of the study by kazushige et al. [45] showed that the cushion area could be divided as shown in Figure 3, with zone A being the 4 cm × 4 cm area under the sciatic bone, zone B being the 10 cm × 10 cm area around zone A, zone C being the entire buttock area, and zone D being the area outside zone C and including the leg area, and based on the results of this study, the size of the sciatic tuberosity area was obtained. This was to facilitate the acquisition and processing of body pressure distribution data, to predict the position of the sciatic tuberosity, to ensure that the sciatic tuberosity position was basically the same when the subjects were in the experiment, and to ensure that the data was not shifted when the body pressure distribution matrix was subsequently averaged to ensure accuracy.

Standardization
The results of the study by kazushige et al. [46] showed that the cushion area could be divided as shown in Figure 3, with zone A being the 4 cm × 4 cm area under the sciatic bone, zone B being the 10 cm × 10 cm area around zone A, zone C being the entire buttock area, and zone D being the area outside zone C and including the leg area, and based on the results of this study, the size of the sciatic tuberosity area was obtained. This was to facilitate the acquisition and processing of body pressure distribution data, to predict the position of the sciatic tuberosity, to ensure that the sciatic tuberosity position was basically the same when the subjects were in the experiment, and to ensure that the data was not shifted when the body pressure distribution matrix was subsequently averaged to ensure accuracy.

Statistical Analysis
The acquisition of the ideal pressure distribution matrix needed to be determined by calculating the body pressure distribution parameters. The meaning of each parameter is described as follows.

Measurement Data Analysis
The contact area, maximum pressure, and average pressure were directly read out by BPMS. The ASCII matrix pressure distribution information saved during the experiment was imported into Excel. The average pressure gradient and maximum pressure gradient were calculated from the pressure distribution matrix.

Analysis of Subjective Evaluation Results
The subjective evaluation method of the experiment adopted the semantic differentiation method [46,47]. The five-degree adverbs from small to large were numbered 1-5, of which 1 was the least comfortable, 5 was the most comfortable, and the highest score was 20 points. The subjects were asked to subjectively evaluate the corresponding experimental items with numbers during each experiment.

Calculation of the Similarity of Body Pressure Distribution Maps
For matrix similarity, the cosine similarity equation below is used to calculate the similarity of the pressure distribution matrices.
Here Ai and Bi represent the values of pressure points of the ideal pressure distribution matrix and the measured prototype pressure distribution matrix, respectively. The closer the calculated similarity is to one, the more similar the two matrices are. The similarity of the body pressure distribution maps between males and between males and females was calculated separately, and the similarity values were averaged. The obtained results are shown in Table 3. The results show that the similarity of pressure distribution maps between the same sexes is greater than that between opposite sexes. In order to select the best among the best, the pressure distribution maps obtained from the experimental results were ranked according to the subjective evaluation of the subjects, and the top 10% of body pressure distribution maps were selected for comfort.
According to the results of 2.3.1, the pressure distribution maps were classified according to gender and divided into four categories: male upright and forward sitting, female upright and forward sitting. Figure 4 shows the human body pressure distribution of the top 10% of the samples in the comparative test.

Human Body Pressure Distribution Indicators
The key body pressure distribution indicators are shown in Table 4. The seats and subjects with the highest subjective evaluation of comfort in the forward-leaning and upright sitting postures are E/S7, A/S5, D/S7, and F/S6 for males and C/S1, D/S3, E/S3, and F/S2 for females as.

Human Body Pressure Distribution Indicators
The key body pressure distribution indicators are shown in Table 4. The seats and subjects with the highest subjective evaluation of comfort in the forward-leaning and upright sitting postures are E/S7, A/S5, D/S7, and F/S6 for males and C/S1, D/S3, E/S3, and F/S2 for females as.

Ideal Pressure Distribution of Sitting Posture
The top 10% body pressure distribution map of comfort in Section 2.3.2 was differentiated according to gender and averaged to obtain an approximate ideal pressure distribution map, as shown in Figure 5. In Darwin's theory of evolution by natural selection, the averageness hypothesis states [48] that the average value of a feature is better than the extreme values. Then, to a certain extent, the pressure distribution matrix with a higher comfort level can represent the approximate ideal pressure distribution matrix.

Ideal Pressure Distribution of Sitting Posture
The top 10% body pressure distribution map of comfort in Section 2.3.2 was differentiated according to gender and averaged to obtain an approximate ideal pressure distribution map, as shown in Figure 5. In Darwin's theory of evolution by natural selection, the averageness hypothesis states [49] that the average value of a feature is better than the extreme values. Then, to a certain extent, the pressure distribution matrix with a higher comfort level can represent the approximate ideal pressure distribution matrix.

Ideal Pressure Distribution Matrix and Its Zones
The above four approximate ideal pressure distribution matrices are averaged according to gender again to obtain an approximate ideal pressure distribution matrix and then carry out Fuzzy C-means Clustering (FCM) on this matrix, and the mathematical model is as follows, based on the obtained results, the ideal pressure distribution map can be partitioned intuitively and effectively.
where n represents the number of pressure points, c represents the number of zones, vk represents the cluster center of the k-th class, µik represents the degree of membership of the i-th sample belonging to the k-th class, ‖xi − vk‖2 represents the squared Euclidean distance from the sample xi to the cluster center vk, and m represents the fuzzy index, which is generally equal to two. The results of FCM clustering on the matrices are shown in Figure 6. Combined with the actual situation of sitting ergonomics, the ideal pressure distribution was mapped to the size of each area on the seat, as shown in Figure 7. The body pressure distribution indexes of each sub-area of the above pressure distribution were extracted and calculated as shown in Table 5, of which Area A was a non-pressure area and

Ideal Pressure Distribution Matrix and Its Zones
The above four approximate ideal pressure distribution matrices are averaged according to gender again to obtain an approximate ideal pressure distribution matrix and then carry out Fuzzy C-means Clustering (FCM) on this matrix, and the mathematical model is as follows, based on the obtained results, the ideal pressure distribution map can be partitioned intuitively and effectively.
where n represents the number of pressure points, c represents the number of zones, vk represents the cluster center of the k-th class, µik represents the degree of membership of the i-th sample belonging to the k-th class, xi − vk 2 represents the squared Euclidean distance from the sample xi to the cluster center vk, and m represents the fuzzy index, which is generally equal to two. The results of FCM clustering on the matrices are shown in Figure 6.

Determination of Key Ergonomic Node Coordinates
The anatomical software 3D body was used to display the subcutaneous tissue of the buttocks and thighs. The zone partition is made according to the ideal support zone sizes, Combined with the actual situation of sitting ergonomics, the ideal pressure distribution was mapped to the size of each area on the seat, as shown in Figure 7. The body pressure distribution indexes of each sub-area of the above pressure distribution were extracted and calculated as shown in Table 5, of which Area A was a non-pressure area and is not included in the statistics.

Determination of Key Ergonomic Node Coordinates
The anatomical software 3D body was used to display the subcutaneous tissue of the buttocks and thighs. The zone partition is made according to the ideal support zone sizes,

Determination of Key Ergonomic Node Coordinates
The anatomical software 3D body was used to display the subcutaneous tissue of the buttocks and thighs. The zone partition is made according to the ideal support zone sizes, as shown in Figure 6. In this study, 10 buttock and thigh histology experts were invited to use Figure 8 as an object to compare and evaluate the partition of the seven zones for the two groups of subjects. Analysis was conducted based on the evaluation and comparison results. The comparison between zones and corresponding key ergonomic nodes is shown in Table 6. as shown in Figure 6. In this study, 10 buttock and thigh histology experts were invited to use Figure 8 as an object to compare and evaluate the partition of the seven zones for the two groups of subjects. Analysis was conducted based on the evaluation and comparison results. The comparison between zones and corresponding key ergonomic nodes is shown in Table 6.  According to the human dimensions of Chinese adults provided by the China national standard GB10000-88, 573 ergonomic data samples of human body dimensions in sitting posture were collected using a shape ruler and a Martin measuring instrument. The samples were classified according to gender, and the contour curves of the buttocks and thighs of different genders were obtained. The ergonomic nodes corresponding to each zone in the curve are marked, and the median of the coordinates of each node in the group was calculated. The results are shown in Figure 9.   According to the human dimensions of Chinese adults provided by the China national standard GB10000-88, 573 ergonomic data samples of human body dimensions in sitting posture were collected using a shape ruler and a Martin measuring instrument. The samples were classified according to gender, and the contour curves of the buttocks and thighs of different genders were obtained. The ergonomic nodes corresponding to each zone in the curve are marked, and the median of the coordinates of each node in the group was calculated. The results are shown in Figure 9. as shown in Figure 6. In this study, 10 buttock and thigh histology experts were invited to use Figure 8 as an object to compare and evaluate the partition of the seven zones for the two groups of subjects. Analysis was conducted based on the evaluation and comparison results. The comparison between zones and corresponding key ergonomic nodes is shown in Table 6.  According to the human dimensions of Chinese adults provided by the China national standard GB10000-88, 573 ergonomic data samples of human body dimensions in sitting posture were collected using a shape ruler and a Martin measuring instrument. The samples were classified according to gender, and the contour curves of the buttocks and thighs of different genders were obtained. The ergonomic nodes corresponding to each zone in the curve are marked, and the median of the coordinates of each node in the group was calculated. The results are shown in Figure 9.

Protype of Ideal Cushion
By mapping the corresponding ergonomic node positions of each zone obtained above to the seat model, the shape and size of the seat after compression can be obtained, as shown in Figure 10.

Protype of Ideal Cushion
By mapping the corresponding ergonomic node positions of each zone ob above to the seat model, the shape and size of the seat after compression can be obt as shown in Figure 10. According to Equation (3), expressing the equivalent elastic coefficient in our ous study [50], the height and pressure distribution values of each area are used calculation to obtain the equivalent elastic coefficient required by each area.
where ki is the comprehensive stiffness coefficient of the material of the i-th zone, the total pressure of the i-th zone, and ΔHi is the height change.
According to the experimental method for studying the relationship between diameter and compression deformation [50], memory foam with an 8 cm thick cu pad, a density of 35D, and a staggered hole spacing of 20 mm was selected. The po ameter D is a variable. A rod with a 10 cm diameter disk end is used to simulate pressure. The results are obtained by linear fitting. The relationship between the po ameter and the equivalent elastic coefficient is shown in Table 7. The pore diameter is used as an independent variable to fit linearly to the equi elastic coefficient k value, as shown in Figure 11. The exponential function that has goodness of fit is shown below: According to Equation (3), expressing the equivalent elastic coefficient in our previous study [49], the height and pressure distribution values of each area are used in the calculation to obtain the equivalent elastic coefficient required by each area.
where k i is the comprehensive stiffness coefficient of the material of the i-th zone, ∆Fi is the total pressure of the i-th zone, and ∆Hi is the height change.
According to the experimental method for studying the relationship between pore diameter and compression deformation [49], memory foam with an 8 cm thick cushion pad, a density of 35D, and a staggered hole spacing of 20 mm was selected. The pore diameter D is a variable. A rod with a 10 cm diameter disk end is used to simulate head pressure. The results are obtained by linear fitting. The relationship between the pore diameter and the equivalent elastic coefficient is shown in Table 7. The pore diameter is used as an independent variable to fit linearly to the equivalent elastic coefficient k value, as shown in Figure 11. The exponential function that has good goodness of fit is shown below: Healthcare 2023, 11, x 12 of 19 Figure 11. Pore size-equivalent elastic coefficient regression curve.
The distribution of the mapped hole matrix is obtained using the RHINO parametric modeling tool, as shown in Figure 12a,b. According to these distributions, 35D polyester fiber is used as the material to make a standard ideal support pillow prototype, as shown in Figure 12c. The equivalent elastic coefficient values of the key ergonomic node mapping positions in each area are measured, as shown in Table 8. Figure 11. Pore size-equivalent elastic coefficient regression curve.
The distribution of the mapped hole matrix is obtained using the RHINO parametric modeling tool, as shown in Figure 12a,b. According to these distributions, 35D polyester fiber is used as the material to make a standard ideal support pillow prototype, as shown in Figure 12c. The equivalent elastic coefficient values of the key ergonomic node mapping positions in each area are measured, as shown in Table 8. The distribution of the mapped hole matrix is obtained using the RHINO parametric modeling tool, as shown in Figure 12a,b. According to these distributions, 35D polyester fiber is used as the material to make a standard ideal support pillow prototype, as shown in Figure 12c. The equivalent elastic coefficient values of the key ergonomic node mapping positions in each area are measured, as shown in Table 8.

Experiment on Pressure Distribution of the Model of Ideal Cushion Prototype
Five eligible subjects from each of the two groups participated in the verification experiment using the corresponding ideal cushion prototype models and following the procedure described in Section 2. The chair with the highest comprehensive score in the first experiment was selected for the control group. The subjects were asked to make subjective scores after the test, and the evaluation indicators included firmness, wrapping, support, and fit. Each item was rated on a five-point scale, with one being the least comfortable and five being the most comfortable, for a maximum possible score of 30. The resulting pressure profile is shown in Figure 13.

Experiment on Pressure Distribution of the Model of Ideal Cushion Prototype
Five eligible subjects from each of the two groups participated in the verification experiment using the corresponding ideal cushion prototype models and following the procedure described in Section 2. The chair with the highest comprehensive score in the first experiment was selected for the control group. The subjects were asked to make subjective scores after the test, and the evaluation indicators included firmness, wrapping, support, and fit. Each item was rated on a five-point scale, with one being the least comfortable and five being the most comfortable, for a maximum possible score of 30. The resulting pressure profile is shown in Figure 13.

Comparative Analysis of Data for Each Sex
The analysis method described in Section 2 is used for the comparative analysis of data on the two sitting postures. The average pressure, peak pressure, maximum pressure gradient, average pressure gradient, and ideal values of each zone are compared between the prototype and the control group ( Figure 14). The following information can be found, as shown by the mean pressure, the data range of each region of the prototype is significantly narrower and more concentrated in distribution, and the median value of the prototype is closer to the ideal value described in Table 4 when compared with the control group. The values of subdivisions B and D are closer to the ideal value, with only +0.09 and +0.11 for males and females in subdivision B, respectively, and −0.01 and-The peak pressure structure is similar to the mean pressure, with the median values in each region of the prototype being closer to the ideal, especially in subpart B, which shows a better approximation with a difference of +0.46 and −0.07. In terms of the maximum pressure

Comparative Analysis of Data for Each Sex
The analysis method described in Section 2 is used for the comparative analysis of data on the two sitting postures. The average pressure, peak pressure, maximum pressure gradient, average pressure gradient, and ideal values of each zone are compared between the prototype and the control group ( Figure 14). The following information can be found, as shown by the mean pressure, the data range of each region of the prototype is significantly narrower and more concentrated in distribution, and the median value of the prototype is closer to the ideal value described in Table 4 when compared with the control group. The values of subdivisions B and D are closer to the ideal value, with only +0.09 and +0.11 for males and females in subdivision B, respectively, and −0.01 and-The peak pressure structure is similar to the mean pressure, with the median values in each region of the prototype being closer to the ideal, especially in subpart B, which shows a better approximation with a difference of +0.46 and −0.07. In terms of the maximum pressure gradient, region B of the prototype has an advantage of +0.06 and +0.19 for males and females, respectively. On the contrary, in females, the values in regions C and D are significantly lower than the ideal values. In terms of mean pressure gradient, the median value distribution of each subdivision of the prototype was better than that of the control group, but the difference was not significant, with region C being the closest and almost equal to the ideal value. In summary, the data distribution of the prototype tends to be more concentrated. The data range of the prototype is significantly reduced, and the median values are closer to the ideal values. In particular, the median values of zones C and D are closer to the ideal value. The body pressure data in each zone of the prototype are superior to the data of the control group, among which the peak pressure and maximum pressure gradient are superior. The data groups of maximum pressure gradient and average pressure gradient have fewer data scattering and tighter distributions. However, for the maximum pressure values, the prototype has a certain deviation from the control group. This may be because the hole size distribution disperses the pressure, avoiding pressure concentration. Thus, it can be judged that each body shape of the office population can be included because the hole array memory foam pad can better restore the pressure distribution map of the mesh support. The ideal pressure distribution map is not limited to a single group of people because it can cover all body shapes within the office population.
In addition to restoring the ideal body pressure indicators in each partition, the similarity of the pressure distribution matrix is also an important basis for evaluating whether the prototype has restored the ideal support surface. The prototype, reference sample, ordinary sponge cushion and ideal pressure distribution matrix were calculated separately for similarity, and the results are shown in Figure 15 below. Among them, the ordinary sponge cushion is the average perforated sponge cushion (hole specification is the smallest hole diameter of 4.2), and the size specifications are all consistent with the sponge used in the prototype. The results show that the similarity between each sample and the ideal pressure distribution map is much higher than that of the ordinary sponge cushion, and the performance is also superior compared with the reference sample. The similarity of the prototypes to the ideal pressure distribution matrix was relatively consistent across all subjects in the male data, indicating that the ideal support for men was achieved. In The body pressure data in each zone of the prototype are superior to the data of the control group, among which the peak pressure and maximum pressure gradient are superior. The data groups of maximum pressure gradient and average pressure gradient have fewer data scattering and tighter distributions. However, for the maximum pressure values, the prototype has a certain deviation from the control group. This may be because the hole size distribution disperses the pressure, avoiding pressure concentration. Thus, it can be judged that each body shape of the office population can be included because the hole array memory foam pad can better restore the pressure distribution map of the mesh support. The ideal pressure distribution map is not limited to a single group of people because it can cover all body shapes within the office population.
In addition to restoring the ideal body pressure indicators in each partition, the similarity of the pressure distribution matrix is also an important basis for evaluating whether the prototype has restored the ideal support surface. The prototype, reference sample, ordinary sponge cushion and ideal pressure distribution matrix were calculated separately for similarity, and the results are shown in Figure 15 below. Among them, the ordinary sponge cushion is the average perforated sponge cushion (hole specification is the smallest hole diameter of 4.2), and the size specifications are all consistent with the sponge used in the prototype. The results show that the similarity between each sample and the ideal pressure distribution map is much higher than that of the ordinary sponge cushion, and the performance is also superior compared with the reference sample. The similarity of the prototypes to the ideal pressure distribution matrix was relatively consistent across all subjects in the male data, indicating that the ideal support for men was achieved. In contrast, the female data showed greater fluctuations, with no significant difference in the similarity of the prototypes compared to the reference sample. Only some (3/10) of the samples had better similarity than that of the reference sample, with the remaining others (7/20) having slightly lower similarity than the control pillows. The support surface requirements for women are more complex compared to men and require more experimental studies. Uniform pressure distribution can characterize better comfort, according to Ahmadian [51], where they characterized the ability of the cushion to produce uniform pressure distribution. The seat pressure distribution index (Seat Pressure Distribution SPD%) was proposed and calculated as follows: where Pi is the pressure on the i-th unit; Pm is the average pressure; n is the total number of pressure points with non-zero values; SPD% can be used for the calculation of pressure distribution in static and dynamic environments, SPD% characterizes the uniformity of the overall body pressure distribution on the cushion, the smaller the SPD%, the more uniform and the higher the overall comfort evaluation of the seat. The results are shown in Table 9 below. It can be seen that the SPD% value of the sample in upright sitting posture is lower than that of the reference sample, the values of the reference sample in forward sitting posture are higher than those of the reference sample, and the SPD% values of the male sample are similar to those of the reference sample, while the female data fluctuate relatively more. In the upright sitting posture, the sample shape surface and human buttocks and thighs fit better, SPD% is relatively small, the overall distribution of pressure is more uniform, to a certain extent, has been able to restore the overall comfort of the sample; forward sitting posture, male and female sample fit compared to the reference sample fit and a certain gap. Compared with the upright sitting posture, the forward sitting posture is more complex, and more factors need to be considered in the design.  Uniform pressure distribution can characterize better comfort, according to Ahmadian [50], where they characterized the ability of the cushion to produce uniform pressure distribution. The seat pressure distribution index (Seat Pressure Distribution SPD%) was proposed and calculated as follows: where Pi is the pressure on the i-th unit; Pm is the average pressure; n is the total number of pressure points with non-zero values; SPD% can be used for the calculation of pressure distribution in static and dynamic environments, SPD% characterizes the uniformity of the overall body pressure distribution on the cushion, the smaller the SPD%, the more uniform and the higher the overall comfort evaluation of the seat. The results are shown in Table 9 below. It can be seen that the SPD% value of the sample in upright sitting posture is lower than that of the reference sample, the values of the reference sample in forward sitting posture are higher than those of the reference sample, and the SPD% values of the male sample are similar to those of the reference sample, while the female data fluctuate relatively more. In the upright sitting posture, the sample shape surface and human buttocks and thighs fit better, SPD% is relatively small, the overall distribution of pressure is more uniform, to a certain extent, has been able to restore the overall comfort of the sample; forward sitting posture, male and female sample fit compared to the reference sample fit and a certain gap. Compared with the upright sitting posture, the forward sitting posture is more complex, and more factors need to be considered in the design.

Subjective Scoring Analysis
The overall subjective score shows a downward trend with the experiment time. The score of each body part is calculated by Equation (6): (6) where N is the mean of each subjective evaluation item, and Nmax is the maximum value of subjective evaluation. The larger the C value, the higher the subjects are evaluated. According to the equation, the total comfort value of the subjective evaluation obtained by the above questionnaire is calculated, as shown in Table 10. The group with the best subjective evaluation in the first human body pressure test experiment is selected as the reference group. Calculation shows that the subjective evaluation values of the prototype model are closer to the control group. The subjective score in the upright sitting position is higher than that in the control group. However, during working hours, the duration of the forward-leaning sitting posture is longer than that of the upright sitting posture, so the four major partitions of the seat surface for forward-leaning may not be enough; further partitions may be required.
The prototypes are sorted from low to high according to the similarity of the ideal pressure distribution matrices, and the subjective evaluation scores of the prototypes obtained are compared, as shown in Figure 16. The trend between subjective evaluation and similarity is somewhat consistent. The agreement between the subjective measurements and similarity was slightly greater for men than for women. It can be seen that the similarity with the ideal pressure distribution matrix characterizes the comfort evaluation to a large extent and also reflects that the ability to reproduce the ideal pressure distribution matrix is an important indicator for evaluating the comfort of the sitting surface.

Conclusions
The main conclusions of this study are as follows: (1) Through the analysis of the ideal pressure distribution map, the sitting contact surface can be divided into three zones: the ischial zone, the buttock zone, and the thigh zone. The ideal body pressure distribution index of each zone is shown in Table 4. (2) The similarity of the pressure distribution of a seat cushion to the ideal pressure distribution reflects comfort to some extent. The similarity of the pressure distribution matrices is compared with the subjective comfort scores. The two are highly consistent. Whether or not the ideal pressure distribution matrix can be reproduced more accurately is an important indicator for evaluating the comfort of the seat cushion. (3) The partitioned hole matrix mapping method can reproduce pressure distribution maps between different materials (e.g., mesh surface and memory foam). This method can be used to design memory foam chairs that perform competitively with existing mesh chairs and can be fabricated more easily for a lower cost. Therefore, this method is important for the design and development of ergonomic chairs.

Future Prospects
Through the experiment of the body pressure distribution of existing mesh chairs, the ideal pressure distribution map of office workers in the two most common sitting postures during working hours was studied. Based on this, the optimization scheme for the seat surface of elastic support materials is designed. However, we cannot rule out the possibility that we have missed a better support scheme, so the resultant ideal pressure distribution matrix can only be said to be relatively optimal to a certain extent. In addition, in the model prototype back-testing experiment, only gender is classified, and the differences in physique among populations have not been refined. Furthermore, in the course of the experiment, factors such as conversations with researchers could cause unnatural changes in posture, and the experimental posture is relatively simple and does not take into account the non-standard sitting postures of users. All these limitations should be addressed in future research.

Conclusions
The main conclusions of this study are as follows: (1) Through the analysis of the ideal pressure distribution map, the sitting contact surface can be divided into three zones: the ischial zone, the buttock zone, and the thigh zone. The ideal body pressure distribution index of each zone is shown in Table 4. (2) The similarity of the pressure distribution of a seat cushion to the ideal pressure distribution reflects comfort to some extent. The similarity of the pressure distribution matrices is compared with the subjective comfort scores. The two are highly consistent. Whether or not the ideal pressure distribution matrix can be reproduced more accurately is an important indicator for evaluating the comfort of the seat cushion. (3) The partitioned hole matrix mapping method can reproduce pressure distribution maps between different materials (e.g., mesh surface and memory foam). This method can be used to design memory foam chairs that perform competitively with existing mesh chairs and can be fabricated more easily for a lower cost. Therefore, this method is important for the design and development of ergonomic chairs.

Future Prospects
Through the experiment of the body pressure distribution of existing mesh chairs, the ideal pressure distribution map of office workers in the two most common sitting postures during working hours was studied. Based on this, the optimization scheme for the seat surface of elastic support materials is designed. However, we cannot rule out the possibility that we have missed a better support scheme, so the resultant ideal pressure distribution matrix can only be said to be relatively optimal to a certain extent. In addition, in the model prototype back-testing experiment, only gender is classified, and the differences in physique among populations have not been refined. Furthermore, in the course of the experiment, factors such as conversations with researchers could cause unnatural changes in posture, and the experimental posture is relatively simple and does not take into account the non-standard sitting postures of users. All these limitations should be addressed in future research.