Construction and Simulation of Biomechanical Model of Human Hip Joint Muscle-Tendon Assisted by Elastic External Tendon by Hill Muscle Model

Luo, Xi; Cai, Guofeng; Ma, Kun; Cai, Aiqi

doi:https://doi.org/10.1155/2022/1987345

Computational Intelligence and Neuroscience

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Special Issue

Exploration of Human Cognition using Artificial Intelligence in Healthcare

View this Special Issue

Research Article | Open Access

Volume 2022 | Article ID 1987345 | https://doi.org/10.1155/2022/1987345

Construction and Simulation of Biomechanical Model of Human Hip Joint Muscle-Tendon Assisted by Elastic External Tendon by Hill Muscle Model

Xi Luo,¹Guofeng Cai,²Kun Ma,¹and Aiqi Cai³

Academic Editor: Arpit Bhardwaj

Received18 May 2022

Revised08 Jun 2022

Accepted14 Jun 2022

Published02 Aug 2022

Abstract

Based on the Hill muscle model (HMM), a biomechanical model of human hip muscle tendon assisted by elastic external tendon (EET) was preliminarily established to investigate and analyze the biomechanical transition between the hip joint (HJ) and related muscle tendons. Using the HMM, the optimal muscle fiber length and muscle force scaling variables were introduced by means of constrained optimization problems and were optimized. The optimized HMM was constructed with human parameters of 170 cm and 70 kg. The biomechanical model simulation test of the hip muscle tendon was performed in the automatic dynamic analysis of mechanical systems (ADAMS) software to analyze and optimize the changes in the root mean square error (RMSE), biological moment, muscle moment distribution coefficient (MDC), muscle moment, muscle force, muscle power, and mechanical work of the activation curves of the hip major muscle, iliopsoas muscle, rectus femoris muscle, and hamstring muscle under analyzing the optimized HMM and under different EET auxiliary stiffnesses from the joint moment level, joint level, and muscle level, respectively. It was found that the trends of the output joint moment of the optimized HMM and the biological moment of the human HJ were basically the same, r² = 0.883 and RMSE = 0.18 Nm/kg, and the average metabolizable energy consumption of the HJ was (243.77 ± 1.59) J. In the range of 35%∼65% of gait cycle (GC), the auxiliary moment showed a significant downward trend with the increase of EET stiffness, when the EET stiffness of the human body was less than 200 Nm/rad, the biological moment of the human HJ gradually decreased with the increase of EET stiffness, and the MDC of the iliopsoas and hamstring muscles gradually decreased; when the EET stiffness was greater than 200 Nm/rad, the increase of the total moment of the extensor muscles significantly increased, the MDC of the gluteus maximus and rectus muscles gradually increased, and the gluteus maximus and hamstring muscle moments and muscle forces gradually increased; the results show that the optimized muscle model based on Hill can reflect the law of human movement and complete the simulation test of HJ movements, which provides a new idea for the analysis of energy migration in the musculoskeletal system of the lower limb.

1. Introduction

All joints and muscles of the human lower limb are involved in the process of natural walking, and the dynamics between each joint and muscle are in the process of change during a complete change cycle, and the structural characteristics of the body’s joints and their corresponding movement characteristics are the basis of exoskeleton activity [1]. Elastic elements play an important role in the activity of the musculoskeletal system of the lower extremities, while they can have a significant effect on limb metabolic energy consumption (EC) [2]. The hip joint (HJ) is the largest positive work joint in the human lower limb, with stability and a large range of motion [3], and plays an important role in human walking and maintaining body balance [4]. The HJ is mainly composed of two parts: pelvis and femur, which are connected through femoral head and acetabulum. Approximately 21 muscles participated in the movement during the movement of the HJ. The extensor and flexor muscle groups of the HJ cooperate with each other to jointly control the activities of various joints of the human lower limb, and elastic tissues such as ligaments in the human. HJ can assist the muscle tissue to support the movement of the human body [5]. Previous findings have pointed out that during walking, when the body’s auxiliary joints show positive work, they can significantly reduce body EC [6]. In addition, it has been pointed out that even if the body’s auxiliary joints show negative work, the body EC can be similarly reduced by providing damping assistance [7]. Elastic external tendon (EET) tissue has significant characteristics of storing energy and releasing energy, and body energy can be effectively reduced by sharing part of the muscle force during walking or exercise in the lower limbs [8].

In recent years, HJ has received much attention in orthopedic clinical and biomechanical research. The research on the transport biomechanical model of human HJ belongs to the research category of sports biomechanics, which has a wide range of research scope, mainly including biology and metrology, the establishment of biomechanical model, and the computer simulation of biological movement mechanism. Among them, the research of the human body is an important direction in sports biomechanics, which is mainly achieved by modeling. At present, biomechanical models of the HJ have been widely used in biology and biological motion machine simulation [9]. The human joint model is based on motion mechanics, and the research methods of motion biomechanics mainly include the application of multi-body theory to establish kinetic model and human simulation research. Muscle model is an important part of human biomechanical model to explore the musculoskeletal system, and the commonly used muscle motor system models are Hill theoretical, Huxley, and rheological models [10]. Based on the relationship between thigh muscle contraction and heat production in frogs, British physiologist Hill proposed the theory of Hill muscle contraction, which is the theoretical basis for the study of skeletal muscle dynamics and has been widely used in lower limb modeling [11]. The researchers applied Hill theory to skeletal muscle, tendon, and human movement control and established a Hill muscle model (HMM) [12]. The human muscle system is complex, and the biomechanics of the muscle system as well as the analysis of interactions between muscles are difficult. The traditional HMM lacks the measurement of physiological parameters such as the degree of muscle activation and is severely limited during use [13]. Although the optimization algorithm based on Hill theory can evaluate the muscle force in the movement of the musculoskeletal system.Its accuracy is limited, and forward dynamics calculation cannot be performed [14]. At present, scholars have established muscle mechanical models, but most of these models are test-based and mostly physiological models rather than biomechanical models [15]. Passive bone simulation based on energy conversion and migration in the human musculoskeletal system and assisted by EET is still a challenging problem in passive exoskeleton design. Therefore, it is of great significance to establish a biomechanical model of the HJ with relatively few parameters and clear mechanical meaning for orthopedic clinical and biomechanical studies.

Therefore, based on the HMM, the biomechanical model of human hip muscle tendon assisted by EET was preliminarily established by optimizing it, and the simulation experiment was performed based on the model to investigate the biomechanical transition among the HJ and related muscle tendons, providing a new idea for the analysis of hip biomechanics and energy migration of the lower limb musculoskeletal system.

2. Materials and Methods

2.1. HMM

The HMM was used to model a single muscle, which was often expressed as muscle-tendon unit (MTU) and was mainly composed of three parts: contractile element (CE), parallel element (PEE), and series element (SEE) [16]. CE refers to the active force generated by the action within the muscle; PEE refers to the passive tension generated within the muscle; SEE refers to the elastic tendon tissue that connects the muscle to the bone. The calculation method of the passive tendon force of SEE in the HMM can be expressed as follows: indicates the active force generated by CE when subjected to nerve stimulation; indicates the passive tension; is the pinna angle; is the muscle force of MTU; indicates the passive tendon force of SEE.

The MT length is calculated as follows: indicates the actual length of the muscle fiber; indicates the length of SEE; indicates the MT; is the pinna angle.

When , the tendon force and tendon length in the HMM can be expressed as follows:

Human movement is based on the mechanical properties of skeletal muscle [17]. The generation of follows the myofilament sliding theory. The generation of is related to the actual length , contraction speed , and activation degree of muscle fibers. is the maximum muscle force that can be generated during muscle isometric contraction, and is the normalized active muscle force after the maximum isometric normalization of . is the standardized muscle fiber length after normalized by the optimal muscle fiber length , and is the standardized muscle fiber contraction speed normalized by the maximum contraction speed of muscle fiber . is the normalized active muscle force-length relationship curve and is the normalized active muscle force-speed relationship curve.

The calculation method of and can be expressed as follows:

Passive muscle force is produced by passive stretching of non-contractile tissue within the muscle and is related to its own muscle fiber length , and when the muscle length exceeds the optimal muscle fiber length , can be expressed as follows: is the normalized passive muscle force after maximum isometric normalization of . is the normalized muscle fiber length after normalized by the optimal muscle fiber length , and is the normalized active muscle force-length relationship curve.

The generation of SEE passive tendon force is only related to the actual length of its own tendon at the current time. is the passive tendon force-length curve. is the normalized muscle fiber length after normalized by the optimal tendon relaxation length , and its normalized expression can be expressed as follows:

Muscle tissue produces muscle force by consuming metabolic energy. Tendon tissue has high energy conversion characteristics. Its energy conversion and metabolism are in balance, so it can be considered that metabolic energy is not consumed [18]. Based on the mathematical model of muscle EC of a single separated muscle, the metabolic EC of muscle fiber tissue in a gait cycle (GC) is estimated, and the specific metabolic energy power can be expressed as below. is the normalized contraction speed of muscle fibers.

The metabolic EC of muscle in a GC can be obtained by integrating the metabolic energy power .

2.2. Hill Muscle Optimization Model to Establish Data and Parameters

The optimized model based on the HMM was used to establish human parameters [19, 20], which were modeled with human parameters of 170 cm and 70 kg. The length of the thigh connecting rod was 434 mm; the distance of the single-joint muscle insertion point was 210 mm; the relative position of the iliopsoas pelvic attachment point was (26 mm, 22 mm, and 5 mm); the length of the moment arm at the joint was 28 mm, the physiological cross-sectional area of the muscle was 28.9 cm²; the muscle pinna angle was 13.9°; the relative position of the rectus femoris pelvic attachment point was (41 mm, −4 mm, and −36 mm); the length of the moment arm at the joint was 50 mm; the physiological cross-sectional area was 34.8 cm²; the muscle pinna angle was 12.4°; the relative position of gluteus maximus pelvic attachment points was (−88 mm, 66 mm, and −39 mm); the moment arm length at its joint was 60 mm, the physiological cross-sectional area was 46.8 cm², and the muscle pinna angle was 21.0°; the moment arm length at the hamstring joint was 28 mm, the moment arm length at the lumbar joint was 28 mm, the physiological cross-sectional area was 73.0 cm², and the muscle pinna angle was 11.3°.

2.3. Establishment of Hill Muscle Optimization Model

The movement of the human lower limb mainly occurs in the sagittal plane. The HJ MT model was constructed based on the sagittal plane. The specific results are shown in Figure 1. In order to make the model closer to the human movement mode and reduce its complexity, only the main muscles related to HJ movement were considered, including gluteus maximus, iliopsoas, rectus femoris, hamstring muscle, and passive elastic elements, of which hamstring muscle included semitendinosus, semimembranosus, and long head of biceps femoris. The mechanical properties of each muscle were represented in an HMM. The α of the HJ muscle was small, and considering it as 0, then the HMM simplified to the fusiform Hill model. The angle of HJ and knee joint movement was expressed by δ and η, respectively, and the angle of iliopsoas muscle and joint moment arm was β1, the angle of rectus femoris and joint moment arm was β2 and β3, the angle of gluteus maximus and joint moment arm was β4, and the angle of hamstring muscle and joint moment arm was β5 and β6.

(a)

(b)

In order to make the HMM consistent with the biomechanical characteristics of the HJ when the human body is walking, the unknown parameter vector of the model is solved by the way of constraint optimization problem, is the angle variable of the muscle force arm at the joint and the muscle force scaling variable. The calculation method is , and is the variable of the optimal muscle fiber length of the muscle. The calculation method is . The muscle moment calculation method of each muscle unit acting on the HJ of the optimized HMM can be expressed as follows: is the vertical action arm length of each muscle force at a certain time in the GC; and represent the unknown scaling variable.

The calculation method of net joint moment at HJ of optimized HMM is as follows: is the passive flexion moment generated for elastic tissue at the joint, and denotes the flexion moment generated by gluteus maximus, iliopsoas, rectus femoris, and hamstring muscle, respectively.

The average difference between the model joint moment and the experimentally obtained human HJ biological moment over a GC is expressed as root mean square error (RMSE). is the number of sampling points. Combining the constraints of each muscle parameter, the constrained optimization problem is expressed as follows:

The standardized passive muscles were calculated according to the known standardization combined with PEE kinetic equation. The standardized active muscle force was obtained according to the CE kinetic equation, and then the actual muscle force at a certain time of each muscle GC was calculated according to the maximum equal moment of the standardized and unknown scaled variables, was calculated using the geometric relationship of the musculoskeletal system; the net joint moment of the musculoskeletal model was obtained to obtain the moment error value; the model parameters were optimized using the RMSE between the output moment of the model and the experimental moment as the solution problem, and the optimal solution of H was obtained using MATLAB genetic algorithm.

The subset I constrained optimization problem for the unknown vector Q can be expressed as follows: denotes the minimum value of muscle activation during one GC.

The optimal variables m and n are substituted into the dynamic equation of model muscle to solve , the normalized muscle fiber length , and actual contraction speed under the corresponding tendon force are calculated according to the force-length dynamic equation of the tendon, and at the same time is substituted into the PEE dynamic equation to solve and . The muscle activation is obtained by introducing them into the CE dynamic equation according to the obtained parameters, and finally the optimal solution of I is obtained according to the constrained optimization problem. The specific flow for solving the scaling variable H and the optimal fiber length variable I is given in Figure 2.

EET can provide auxiliary flexion moment for HJ during lower limb walking, and based on its effect, the external elastic tendon is assumed to be a passive linear torsional spring parallel to the HJ. For an EET with a stiffness of J, its resulting flexion moment during walking of the lower limb can be expressed as follows:

EET acts directly with human HJ flexor muscle group [21]. During the auxiliary period of EET, the calculation method of residual total flexion biological moment provided by rectus femoris and iliopsoas muscle based on Hill muscle optimization model can be expressed as follows: represents the moment of iliopsoas muscle; represents the moment of rectus femoris; represents less than the total flexion moment of flexor at HJ during natural walking.

The additional calculation method of resistance moment provided by gluteus maximus and hamstring muscle based on Hill muscle optimization model is as follows:

The biological moment solution and biomechanical simulation framework of each muscle of the model under the action of EET are shown in Figure 3.

2.4. Simulation Test Conditions of HJ MT Biomechanical Model

The simulation experiment of HJ muscle tendon biomechanical model was mainly carried out in the software of automatic dynamic analysis of mechanical systems (ADAMS). It mainly provided constraint base and force base centered on the user through ADAMS/View module, and gathers the functions of interactive image, simulation calculation, and result analysis together to effectively carry out modeling and simulation. Moreover, MATLAB program was used for programming. The program was mainly composed of data input, data calculation and solution, data output, and graphical interface.

3. Results

3.1. Evaluation of Optimized HMM Performance

Firstly, the biomechanical properties of human lower limb walking under the optimized HMM were evaluated from the joint level. During a GC, the change trends of the output joint moment of the optimized HMM and the biological moment of the human HJ were basically the same, showing a trend of first increase and then decrease. The similarity of the output joint moment of the optimized HMM and the biological moment curve of the human HJ was compared using the coefficient r² and RMSE. The results showed that the HJ moment r² = 0.883, RMSE = 0.18 Nm/kg (Figure 4(a)). The total metabolic EC power during one GC of the HJ under the optimized HMM was statistically analyzed using the muscle energetics equation, and it was in a fluctuating state during one GC, with an average metabolic EC of (243.77 ± 1.59) J (Figure 4(b)).

(a)

(b)

The validity of the optimized HMM was further evaluated from the muscle level, and the results showed that the muscle activation curves of gluteus maximus, iliopsoas, rectus femoris, and hamstring muscle under the optimized HMM tended to be similar to the test activation curve during one GC; r² = 0.875, RMSE = 0.19 Nm/kg for iliopsoas; r² = 0.890, RMSE = 0.15 Nm/kg for rectus femoris; r² = 0.892, RMSE = 0.16 Nm/kg for gluteus maximus; and r² = 0.879, RMSE = 0.18 Nm/kg for hamstring muscle (Figure 5).

(a)

(b)

(c)

(d)

3.2. Analysis of Moment Change of Human HJ under Different EET Stiffness

The changes of human auxiliary moment and human HJ biomechanical moment under different EET stiffness were analyzed. With the increase of EET stiffness during GC, the corresponding auxiliary moment showed a significant trend of first reduction and then increase, and finally tended to be stable. In the range of 35%∼65% of GC, the auxiliary moment showed a significant downward trend with the increase of EET stiffness. When the EET stiffness was 0 Nm/rad, the auxiliary flexion moment showed a stable state (Figure 6(a)). With the increase of EET stiffness, in the range of 35%∼65% of GC, when the EET stiffness of human body was less than 200 Nm/rad, the biological moment of human HJ gradually decreased with the increase of EET stiffness, and when the EET stiffness of human body was greater than 200 Nm/rad, the biological moment of human HJ gradually increased with the increase of EET stiffness (Figure 6(b)).

(a)

(b)

By analyzing the changes of the total flexor moment of the HJ under different EET stiffness (Figure 7(a)), the total flexor moment of the HJ increased significantly with the increase of the EET stiffness in the range of 35% ∼ 65% of the GC. When the EET stiffness was less than 200 Nm/rad, the total flexor moment of the HJ decreased first and then increased with the decrease of the EET stiffness (Figure 7(a)). With the increase of EET stiffness, the total moment of human HJ extensor increased significantly in the range of 35%∼65% of GC. When the EET stiffness was greater than 200 Nm/rad, the increase range of total moment of extensor increased significantly (Figure 7(b)).

(a)

(b)

3.3. Analysis on the Change of Muscle Moment Distribution Coefficient (MDC) under Different EET Stiffness

The change of muscle MDC under different EET stiffness was analyzed. With the increase of EET stiffness, the change trend of MDC of different muscles was significantly different. When the stiffness of the external tendon was less than 200 Nm/rad, the MDC of iliopsoas muscle gradually decreased with the increase of the stiffness of the external tendon. When the stiffness of the external tendon was greater than 350 Nm/rad, the MDC of iliopsoas muscle gradually tended to a stable state with the increase of the stiffness of the external tendon (Figure 8(a)). When the stiffness of the external tendon was less than 250 Nm/rad, the rectus femoris MDC gradually increased with the increase of the stiffness of the external tendon. When the stiffness of the external tendon was greater than 250 Nm/rad, the rectus femoris MDC gradually decreased to zero with the increase of the stiffness of the external tendon (Figure 8(b)). With the increase of external tendon stiffness (ETS), gluteus maximus MDC gradually increased (Figure 8(c)). With the increase of ETS, hamstring muscle MDC showed a gradual downward trend (Figure 8(d)).

(a)

(b)

(c)

(d)

3.4. Analysis of Muscle Moment Changes under Different EET Stiffness

The moment changes of different muscles under different EET stiffness were analyzed. With the increase of EET stiffness, the moment changes of different muscles were obviously different. In the range of 35% ∼ 65% of GC, with the increase of ETS, the moment of iliopsoas and rectus femoris muscle gradually decreased. When the ETS was greater than 600 Nm/rad, the moment of iliopsoas and rectus femoris gradually approached zero (Figures 9(a) and 9(b)). In the range of 35% ∼ 65% of GC, the muscle moment of gluteus maximus and hamstring muscle gradually increased with the increase of ETS (Figures 9(c) and 9(d)).

(a)

(b)

(c)

(d)

3.5. Analysis of Muscle Force Change under Different EET Stiffness

The force changes of each muscle under different EET stiffness were analyzed, and the moment changes of different muscles were significantly different with the increase of EET stiffness. In the range of 35% ∼ 65% of GC, the iliopsoas muscle force gradually decreased with the increase of ETS (Figure 10(a)); in the range of 15% ∼ 65% of GC, the rectus femoris muscle force gradually decreased with the increase of ETS (Figure 10(b)); and in the range of 35% ∼ 60% of GC, the gluteus maximus and hamstring muscle forces gradually increased with the increase of ETS (Figures 10(c) and 10(d)).

(a)

(b)

(c)

(d)

3.6. Analysis of Muscle Power Changes under Different EET Stiffness

The change of muscle power under different EET stiffness was analyzed. With the increase of EET stiffness, the change trend of different muscle power was clearly different. In a GC, the muscle power of iliopsoas muscle decreased first and then increased. In the range of 55% ∼ 65% of the GC, with the increase of ETS, the muscle power of iliopsoas muscle gradually increased (Figure 11(a)); in one GC, the muscle power of rectus femoris decreased first and then increased. In the range of 60% ∼ 85% of GC, the muscle power of rectus femoris gradually decreased and tended to zero with the increase of ETS (Figure 11(b)). In the range of 35% ∼ 65% of GC, the muscle power of gluteus maximus gradually increased with the increase of ETS (Figure 11(c)). In the range of 35% ∼ 65% of GC, with the increase of ETS, hamstring muscle power gradually increased. When the ETS was 0, hamstring muscle power was zero (Figure 11(d)).

(a)

(b)

(c)

(d)

The changes of muscle mechanical work under different ETS were analyzed. With the increase of EET stiffness, the mechanical work of HJ decreased gradually. When the EET stiffness was less than 200 Nm/rad, the mechanical work of muscle unit showed a downward trend, and when the EET stiffness was greater than 200 Nm/rad, the mechanical work of muscle unit and muscle fiber showed an upward trend. With the increase of EET stiffness, the mechanical work of tendon did not change. With the increase of EET stiffness, the mechanical work of knee joint gradually decreased (Figure 12).

4. Discussion

The results showed that the changing trends of the output joint moment of the optimized HMM and the biological moment of the human HJ were basically the same, showing a trend of first increase and then decrease, with the HJ moment r² = 0.883 and RMSE = 0.18 Nm/kg. For the general musculoskeletal model, when the moment curve r² was greater than 0.883 and the RMSE value was less than 0.2 Nm/kg, it could be considered that the moment curve fitting had a higher goodness, and the established model was close to the mechanical characteristics of the human musculoskeletal [22]. Optimized HMM can effectively reproduce the biological characteristics of human HJ at joint moment level. The characteristics of the established model were further evaluated from the muscle level, and the results showed that the muscle activation curves of gluteus maximus, iliopsoas, rectus femoris, and hamstring muscles under the optimized HMM had similar trends to the test activation curves, and the shape of the activation curves of different muscle tests was similar to the results of Bogey and Barnes (2017) [23]. The total metabolic energy consumption power curve of the model muscle in a gait cycle was calculated by muscle energy equation. The result revealed the average metabolic EC of the HJ was (243.77 ± 1.59) J. It has been pointed out that the total metabolic EC rate produced during normal walking in the human body was (5.20 ± 0.91) W/kg [24], the net metabolic EC rate was about (3.39 ± 0.58) W/kg [25], and the net metabolic rate was (1.81 ± 0.33) W/kg [26]. During walking in the human body, the energy consumed by the HJ muscle accounted for 9/20 of the total metabolic energy. Therefore, the total metabolic EC rate of walking in the human body based on the optimized HMM was 6.22 W/kg, which was consistent with the metabolic EC level of the HJ muscle found in the current relevant studies [27, 28]. In conclusion, the optimized HMM can effectively reproduce the biological characteristics of human HJ at joint moment level, joint level, and muscle level.

With the increase of EET stiffness in GC, the corresponding auxiliary moment showed a significant trend of first reduction and then increased and finally tended to be stable, and in the range of 35% ∼ 65% of GC the auxiliary moment showed a significant downward trend with the increase of EET stiffness. When the EET stiffness of the human body was less than 200 Nm/rad, the biological moment of the human HJ gradually decreased with the increase of the EET stiffness, and when the EET stiffness of the human body was greater than 200 Nm/rad, the biological moment of the human HJ gradually increased with the increase of the EET stiffness, and the total moment of the extensor muscle increased significantly. There was a significant antagonistic effect among human flexor and extensor muscle groups [29], and the flexion moment generated by the HJ flexor muscles was obviously greater than the joint moment during human walking [30]. When the elastic external stiffness increased to 200 Nm/rad, the HJ moment value was significantly smaller than the auxiliary flexion moment value but the HJ extensor moment did not change during the reduction of the flexor moment. When the EET stiffness was greater than 200 Nm/rad, the human HJ flexor muscles can only provide energy through the EET auxiliary interval, which in turn leaded to a decrease in the flexion moment [31], and as the EET stiffness gradually became larger, the HJ moment gradually changed from the flexion moment to the extension moment [32]. In the range of 35% ∼ 65% of GC, the moment of iliopsoas and rectus femoris muscles gradually decreased with the increase of ETS, and the moment of iliopsoas and rectus femoris muscles gradually approached to zero when the ETS was greater than 600 Nm/rad. With the auxiliary of EET, the degree of auxiliary of the iliopsoas was more significant, and gluteus maximus appeared as a resistant muscle when the force generated by EET stiffness gradually increased. This was because monoarticular muscles such as iliopsoas and gluteus maximus were mainly driving joint movements [33], while biarticular muscles such as rectus femoris and hamstring muscle were mainly responsible for energy transfer between joints [34]. With the auxiliary of EET, the human control system tended to be more stable control than motor control [35]. In order to maintain the stability of the human motor system, single-joint muscles play an important role in it [36]. Therefore, the biomechanical parameters of iliopsoas and gluteus maximus changed most significantly with the auxiliary of EET. In the range of 35% ∼ 65% of GC, with the increase of ETS, the iliopsoas muscle force gradually increased, the rectus femoris muscle force gradually decreased, and the gluteus maximus and hamstring muscle forces gradually increased. The force of both iliopsoas and rectus femoris muscles decreased as the stiffness of EET auxiliary increased, and their corresponding muscle fiber lengths showed an increasing trend. The compliance of the tendon can have a significant effect on the length of the muscle fibers [37], and the tendon compliance of the iliopsoas muscle was small and short in length, so its muscle fiber length increased less; the tendon compliance of rectus femoris was long and flexible, but its muscle force decreased less, which in turn leaded to a smaller increase in its muscle fiber length.

5. Conclusion

A biomechanical model of human HJ MT assisted by EET was established based on the HMM, and it was found by simulation experiment that muscle activation had a significant effect on HJ muscle moment, muscle force as well as mechanical work changes with EET assistance, and metabolic EC during human walking could be reduced with EET assistance. However, there are still some shortcomings, only the biomechanical changes of human HJ MT assisted by EET are preliminarily explored and they are not compared with the ankle joint. In the future, the biomechanical properties of the two will be further compared. The optimized HMM can reflect the law of human motion, and complete the simulation test of typical movements of HJ, which provides a new idea for the analysis of energy migration of lower limb musculoskeletal system.

Data Availability

The datasets used during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

Xi Luo and Guofeng Cai contributed equally to this work and are co-1st author of this paper.

Acknowledgments

This work was supported by the Analysis of Testing Project of Kunming University of Science and Technology, Research on Biomechanical Analysis for Experimental of Hip Femoral Impingement Syndrome (FAI) (Project no. 2021P20191110002).

References

L. D. Mendonça, J. M. Ocarino, N. F. N. Bittencourt, L. G. Macedo, and S. T. Fonseca, “Association of hip and foot factors with patellar tendinopathy (Jumper’s Knee) in Athletes,” Journal of Orthopaedic & Sports Physical Therapy, vol. 48, no. 9, pp. 676–684, 2018.
View at: Publisher Site | Google Scholar
H. Kainz and M. H. Schwartz, “The importance of a consistent workflow to estimate muscle-tendon lengths based on joint angles from the conventional gait model,” Gait & Posture, vol. 88, pp. 1–9, 2021.
View at: Publisher Site | Google Scholar
E. A. Audenaert, V. Khanduja, C. Bauwens, T. Van Hoof, C. Pattyn, and G. Steenackers, “A discrete element model to predict anatomy of the psoas muscle and path of the tendon: Design implications for total hip arthroplasty,” Clinical Biomechanics, vol. 70, pp. 186–191, 2019.
View at: Publisher Site | Google Scholar
W. Chen, S. Wu, T. Zhou, and C. Xiong, “On the biological mechanics and energetics of the hip joint muscle-tendon system assisted by passive hip exoskeleton,” Bioinspiration & Biomimetics, vol. 14, no. 1, Article ID 016012, 2018.
View at: Publisher Site | Google Scholar
D. Chicoine, M. Bouchard, S. Laurendeau, G. Moisan, E. L. Belzile, and P. Corbeil, “Biomechanical effects of three types of foot orthoses in individuals with posterior tibial tendon dysfunction,” Gait & Posture, vol. 83, pp. 237–244, 2021.
View at: Publisher Site | Google Scholar
L. Wade, G. A. Lichtwark, and D. J. Farris, “Joint and muscle-tendon coordination strategies during submaximal jumping,” Journal of Applied Physiology, vol. 128, no. 3, pp. 596–603, 1985.
View at: Publisher Site | Google Scholar
E. Martín-Sosa, J. Martínez-Reina, J. Mayo, and J. Ojeda, “Influence of musculotendon geometry variability in muscle forces and hip bone-on-bone forces during walking,” PLoS One, vol. 14, no. 9, Article ID e0222491, 2019.
View at: Publisher Site | Google Scholar
S. Burzyński, A. Sabik, W. Witkowski, and P. Łuczkiewicz, “Influence of the femoral offset on the muscles passive resistance in total hip arthroplasty,” PLoS One, vol. 16, no. 5, Article ID e0250397, 2021.
View at: Publisher Site | Google Scholar
O. Hapa, N. D. Demirkıran, B. Hüsemoğlu, M. Edizer, and H. Havitçioğlu, “Anatomic implications of lesser trochanterplasty,” Acta Orthopaedica et Traumatologica Turcica, vol. 52, no. 1, pp. 54–57, 2018.
View at: Publisher Site | Google Scholar
A. Cerda-Lugo, A. González, A. Cardenas, and D. Piovesan, “Modeling the neuro-mechanics of human balance when recovering from a fall: a continuous-time approach,” BioMedical Engineering Online, vol. 19, no. 1, p. 67, 2020.
View at: Publisher Site | Google Scholar
Y. Tan, Z. Fu, L. Duan et al., “Hill-based musculoskeletal model for a fracture reduction robot,” International Journal of Medical Robotics and Computer Assisted Surgery, vol. 17, no. 3, 14 pages, 2021.
View at: Publisher Site | Google Scholar
T. Siebert, N. Stutzig, and C. Rode, “A hill-type muscle model expansion accounting for effects of varying transverse muscle load,” Journal of Biomechanics, vol. 66, pp. 57–62, 2018.
View at: Publisher Site | Google Scholar
R. Rockenfeller, M. Günther, N. Stutzig et al., “Exhaustion of Skeletal Muscle Fibers Within Seconds: Incorporating Phosphate Kinetics Into a Hill-Type Model,” Frontiers in Physiology, vol. 11, p. 306, 2020.
View at: Publisher Site | Google Scholar
B. Xiong, N. Zeng, Y. Li et al., “Determining the Online Measurable Input Variables in Human Joint Moment Intelligent Prediction Based on the Hill Muscle Model,” Sensors, vol. 20, no. 4, p. 1185, 2020.
View at: Publisher Site | Google Scholar
P. Bujalski, J. Martins, and L. Stirling, “A Monte Carlo analysis of muscle force estimation sensitivity to muscle-tendon properties using a Hill-based muscle model,” Journal of Biomechanics, vol. 79, pp. 67–77, 2018.
View at: Publisher Site | Google Scholar
F. Heinen, S. R. Sørensen, M. King et al., “Muscle-Tendon Unit Parameter Estimation of a Hill-Type Musculoskeletal Model Based on Experimentally Obtained Subject-Specific Torque Profiles,” Journal of Biomechanical Engineering, vol. 141, no. 6, 2019.
View at: Google Scholar
J. L. Palladino, “Functional Requirements of a Mathematical Model of Muscle Contraction,” in Proceedings of the 2019 41st Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), pp. 6956–6959, Berlin, Germany, July 2019.
View at: Publisher Site | Google Scholar
J. Guo, Y. Sun, Y. Hao, L. Cui, and G. Ren, “A mass-flowing muscle model with shape restrictive soft tissues: correlation with sonoelastography,” Biomechanics and Modeling in Mechanobiology, vol. 19, no. 3, pp. 911–926, 2020.
View at: Publisher Site | Google Scholar
C. Kleinbach, O. Martynenko, J. Promies, D. F. B. Haeufle, J. Fehr, and S. Schmitt, “Implementation and validation of the extended Hill-type muscle model with robust routing capabilities in LS-DYNA for active human body models,” BioMedical Engineering Online, vol. 16, no. 1, p. 109, 2017.
View at: Publisher Site | Google Scholar
Z. Shao, Q. Wu, B. Chen, and H. Wu, “Force and deformation transmission characteristics of a compliant tendon-sheath actuation system based on Hill-type muscle model,” Proceedings of the Institution of Mechanical Engineers - Part H: Journal of Engineering in Medicine, vol. 233, no. 7, pp. 695–705, 2019.
View at: Publisher Site | Google Scholar
R. Rockenfeller and M. Günther, “Hill equation and Hatze’s muscle activation dynamics complement each other: enhanced pharmacological and physiological interpretability of modelled activity-pCa curves,” Journal of Theoretical Biology, vol. 431, pp. 11–24, 2017.
View at: Publisher Site | Google Scholar
A. Mohammadi Nejad Rashty, M. Grimmer, and A. Seyfarth, “Hopping frequency influences elastic energy reuse with joint series elastic actuators,” Journal of Biomechanics, vol. 119, Article ID 110319, 2021.
View at: Publisher Site | Google Scholar
R. A. Bogey and L. A. Barnes, “Estimates of individual muscle power production in normal adult walking,” Journal of NeuroEngineering and Rehabilitation, vol. 14, no. 1, p. 92, 2017.
View at: Publisher Site | Google Scholar
J. C. Martin and J. A. Nichols, “Simulated work-loops predict maximal human cycling power,” Journal of Experimental Biology, vol. 221, no. Pt 13, 2018.
View at: Publisher Site | Google Scholar
Y. P. Lim, Y. C. Lin, and M. G. Pandy, “Effects of step length and step frequency on lower-limb muscle function in human gait,” Journal of Biomechanics, vol. 57, pp. 1–7, 2017.
View at: Publisher Site | Google Scholar
T. Zumbrunn, R. Patel, M. P. Duffy et al., “Cadaver-Specific Models for Finite-Element Analysis of Iliopsoas Impingement in Dual-Mobility Hip Implants,” The Journal of Arthroplasty, vol. 33, no. 11, pp. 3574–3580, 2018.
View at: Publisher Site | Google Scholar
V. R. Allen, R. E. Kambic, S. M. Gatesy, and J. R. Hutchinson, “Gearing effects of the patella (knee extensor muscle sesamoid) of the helmeted guineafowl during terrestrial locomotion,” Journal of Zoology, vol. 303, no. 3, pp. 178–187, 2017.
View at: Publisher Site | Google Scholar
M. Zhu, D. Musson, M. Oliver, E. Firth, J. Cornish, and J. Munro, “Modelling gluteus medius tendon degeneration and repair in a large animal model,” Archives of Orthopaedic and Trauma Surgery, vol. 142, no. 1, pp. 1–12, 2022.
View at: Publisher Site | Google Scholar
J. C. Sedlmayr, K. T. Bates, J. J. Wisco, and E. R. Schachner, “Revision of hip flexor anatomy and function in modern humans, and implications for the evolution of hominin bipedalism,” The Anatomical Record, vol. 305, no. 5, pp. 1147–1167, 2022.
View at: Publisher Site | Google Scholar
T. R. Souza, W. Schallig, K. Veerkamp et al., “Muscle actions on crossed and non-crossed joints during upright standing and gait: A comprehensive description based on induced acceleration analysis,” Journal of Biomechanics, vol. 130, Article ID 110874, 2022.
View at: Publisher Site | Google Scholar
L. C. Tsai, Y. A. Ko, K. E. Hammond, J. W. Xerogeanes, G. L. Warren, and C. M. Powers, “Increasing hip and knee flexion during a drop-jump task reduces tibiofemoral shear and compressive forces: implications for ACL injury prevention training,” Journal of Sports Sciences, vol. 35, no. 24, pp. 2405–2411, 2017.
View at: Publisher Site | Google Scholar
T. J. M. Dick, L. K. Punith, and G. S. Sawicki, “Humans falling in holes: adaptations in lower-limb joint mechanics in response to a rapid change in substrate height during human hopping,” Journal of The Royal Society Interface, vol. 17, no. 163, Article ID 20200047, 2019.
View at: Publisher Site | Google Scholar
T. J. M. Dick, C. J. Clemente, L. K. Punith, and G. S. Sawicki, “Series elasticity facilitates safe plantar flexor muscle-tendon shock absorption during perturbed human hopping,” Proceedings of the Royal Society B: Biological Sciences, vol. 288, no. 1947, Article ID 20210201, 2021.
View at: Publisher Site | Google Scholar
R. W. Nuckols, K. Z. Takahashi, D. J. Farris, S. Mizrachi, R. Riemer, and G. S. Sawicki, “Mechanics of walking and running up and downhill: A joint-level perspective to guide design of lower-limb exoskeletons,” PLoS One, vol. 15, no. 8, Article ID e0231996, 2020.
View at: Publisher Site | Google Scholar
D. Tokur, M. Grimmer, and A. Seyfarth, “Review of balance recovery in response to external perturbations during daily activities,” Human Movement Science, vol. 69, Article ID 102546, 2020.
View at: Publisher Site | Google Scholar
S. E. Miller, A. D. Segal, G. K. Klute, and R. R. Neptune, “Hip recovery strategy used by below-knee amputees following mediolateral foot perturbations,” Journal of Biomechanics, vol. 76, pp. 61–67, 2018.
View at: Publisher Site | Google Scholar
C. E. Shell, A. D. Segal, G. K. Klute, and R. R. Neptune, “The effects of prosthetic foot stiffness on transtibial amputee walking mechanics and balance control during turning,” Clinical Biomechanics, vol. 49, pp. 56–63, 2017.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Xi Luo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies