- Open Access
A simulating analysis of the effects of increased joint stiffness on muscle loading in a thumb
© Wu et al; licensee BioMed Central Ltd. 2009
Received: 18 June 2009
Accepted: 16 December 2009
Published: 16 December 2009
The development of osteoarthritis (OA) in the hand results in increased joint stiffness, which in turn affects the grip strength. The goal of the present study is to theoretically analyze the muscle forces in a thumb in response to the increased joint stiffness.
The thumb was modeled as a linkage system consisting of a trapezium, a metacarpal bone, a proximal and a distal phalanx. Nine muscles were included in the model: flexor pollicis longus (FPL), extensor pollicis longus (EPL), extensor pollicis brevis (EPB), abductor pollicis longus (APL), flexor pollicis brevis (FPB), abductor pollicis brevis (APB), the transverse head of the adductor pollicis (ADPt), the oblique head of the adductor pollicis (ADPo), and opponens pollicis (OPP). Numerical tests were performed using an inverse dynamic approach. The joints were prescribed to an angular motion at one degree-of-freedom (DOF) each time with all other DOFs of the joints being mechanically constrained, while the muscle forces in response to the joint motions were predicted. The normal joint stiffness was assumed to be 0.05, 0.10, and 0.15 N m/rad for interphalangeal (IP), metacarpophalangeal (MCP), and carpometacarpal (CMC) joint, respectively. The joint stiffness was assumed to increase by 50% and 100%, simulating the biomechanical consequences of OA.
Our simulations indicated that the increase in joint stiffness induced substantial increases in muscle forces, especially in the EPL and FPL muscles in response to IP, MCP, or CMC extension/flexion motions.
Because the strength of the muscles in the fingers is limited, the muscles will not be able to overcome joint resistance if joint stiffness is increased to its limit due to OA. This may contribute to the reduced range of motion typically seen in OA.
The development of osteoarthritis (OA) in the hand is associated with difficulties in gripping activities . Previous studies indicated that the joint stiffness could increase by more than 100% due to OA in the hand . The increased joint stiffness in turn affects the grip strength [3–5]. The biomechanics underlaying the interactions between the muscular loading and joint stiffness variations due to OA has not been investigated. Since it is not convenient to experimentally measure the muscle forces in a finger under physiological conditions, biomechanical models of the hand and fingers are useful for studying such problems.
Multiple biomechanical models of the hands and fingers have been developed to simulate different problems; for example, the whole hand models by [6, 7] that simulated the muscle loading for static gripping and free movements, and the biomechanical model of the dynamics of the index finger by  that simulated the muscle forces in pinch grip and disc rotation. More practical biomechanical finger models were proposed by [9, 10] that included anatomically realistic tendon/muscle network connections in their models. Most of past simulation models and studies were developed for healthy normal hands. The effects of the altered joint stiffness due to pathological conditions, such as OA, on the musculoskeletal loading in a finger has not been analyzed to date.
Direct experimental determinations of the relationship between the muscle force and joint motions in the thumb have been performed by [11, 12]. In these experiments, the relationships between the joint motions and muscle forces were tested directly using cadaveric hands. The muscle in the thumb was pulled individually and all other muscles were either loaded by a small force  or left free , while the joint motions in response to the muscle forces were measured. Again, only healthy donors were considered in these experimental studies.
The goal of the present study is to theoretically analyze the muscle force in a thumb in response to increased joint stiffness. Specifically, we are going to examine: (1) the relationship between the joint motion and muscle force for a normal thumb, and (2) the effects of increased joint stiffness on the relationship between the joint motion and the muscle force. The analysis is to be conducted using an inverse dynamic method, i.e., the joint motion is prescribed while the muscle forces are predicted. We hypothesized that the relationship between the joint motion and muscle forces determined using inverse dynamics will be consistent with those observed in the experimental study , and that the elevated joint stiffness will increase the muscle force recruitment.
For the MCP and CMC joints, we assumed that the joint stiffness in extension/flexion was identical to that in abduction/adduction motion. The normal joint stiffness, which includes the joint resistance and the effects of the connective tissues, was assumed to be 0.05, 0.10, and 0.15 N m/rad for the IP, MCP, and CMC joint, respectively. The joint stiffness was assumed to increase by 50% and 100%, simulating the early stage of OA . The neutral position was considered to be 5 and 10 degrees of the IP and MCP flexion, respectively. All other joint angular components were zero at the neutral position.
Numerical tests were performed using an inverse dynamic approach. The joints were prescribed to an angular motion at one degree-of-freedom (DOF) each time with all other DOFs of the joints being mechanically constrained, while the muscles forces in response to the joint motions were predicted. The joints were moving from their prescribed extreme positions within a time period of 10 s and at constant speeds.
The recruitment of the muscle forces was calculated by using a min/max optimization procedure in AnyBody , in which the maximal normalized muscle force was minimized. The minimization of the cost function was subjected to the constraints: the muscle force was greater than or equal to zero and the maximal muscle force did not exceed its capacity, which is estimated by the physiological cross sectional area multiplied by a muscle strength factor of 35 N/cm2 . At any instance, the sums of the contributions of each individual muscle to joint moments were calculated and they were balanced with the external forces and the inertial forces of the segments.
Discussion and Conclusion
Our simulations indicated that the increase in the joint stiffness - a biomechanical consequence in the early stage of OA in the fingers - induces a substantial increase in muscle forces, especially in EPL and FPL muscles in response to the IP, MCP, or CMC extension/flexion motions. Because the strength of the muscles in the fingers is limited, the muscles will not be able to overcome the joint resistance and move the joints through the entire range of motion if the joint stiffness is increased significantly. This explains, in part, why the OA patients suffer a reduced range of motion [19, 4].
The predicted relationships between the joint motion and extrinsic muscle (FPL, EPL, APL, and EPB) activities are generally consistent with those observed in the experimental measurements . The relationship between the joint extension/flexion motions and muscle force responses predicted in our simulations agree well with those observed experimentally: the muscle forces in EPL and FPL are generated in response to the IP extension and flexion motion, respectively (Figure 2); in the MCP joint, the muscle forces in EPL and EPB are generated in response to the extension motion while those in the FPL are generated in response to flexion motion (Figure 3-left column); and in CMC joint, the muscle forces in EPL, APL, and EPB are generated in response to extension while that in FPL is generated in response to flexion motion (Figure 4). The relationship between the joint abduction/adduction motions and muscle force responses predicted in our simulations also agree in general with those observed experimentally: the abduction motions in the CMC and MCP joints induce the EPL muscle force (Figures 3-right column and 5); and adduction/abduction motions in CMC and MCP joints induce a small force response in the APL and EPB muscles (Figure 3-right column and Figure 5).
The only difference between the model predictions and the experimental observations  is in the FPL muscle in response to the adduction/abduction motions. Our predictions indicated that the MCP and CMC abduction/adduction motions generated force in the FPL (Figure 3-right column and Figure 5), which was, however, not observed experimentally . This does not mean that there is a conflict between the model predictions and the experimental observations. The current simulations indicated that, in order to generate the MCP and CMC abduction/adduction motions, the FPL muscle is required to maintain joint stability in the flexion/extension direction. It should be noted that the current study is different in nature from the experimental study : the simulations were performed using an inverse dynamic technique while the experiments  were virtually performed using a forward dynamic principle. In addition, only one muscle was pulled each time and all other muscles had no contributions in the experiments , whereas all nine muscles were considered to participate in the force sharing in the simulations.
Another factor of the current study, which could contribute to the difference between the theoretical predictions and experimental data, is the modeling of the MCP and CMC joint. These two joints, especially the CMC joint, are considered to consist of two "scoliotic saddle-shaped" joint surfaces. In the ideal case, there are only two DOFs (extension/flexion and abduction/adduction motions) and the motion along the axis, i.e., pronation/supination, is negligible in such a joint. In a real CMC joint, the pronation/supination motion is not completely constrained because of the flexibility of the cartilage layers and the space in the joint, as demonstrated in the experimental results .
In the current study, the joint stiffness is considered as linear and time-independent. Previous experimental data indicated, however, that the relationship of the joint-motion-moment in the fingers is nonlinear  and viscous . Because of the effects of the connective tissues, which are nonlinear and viscous, the stiffness property of the joint is also typically nonlinear and viscous. When the joint motion is very slow - the case simulated in the current study - the viscous effects of the joint moment response is negligible. Besides, joint viscosity is conventionally not evaluated in clinical studies. A relationship between the joint viscosity and joint OA conditions has not been established. Therefore, we did not consider the effects of the joint viscosity in the current parametric study.
The curves of the joint moment-angle of the fingers are typically characterized by a hysteresis loop: the joint moment-angle relationship becomes non-linear only towards the ends of the joint motion range, while it is nearly linear around the center of the neutral position. If we fit the joint-motion-moment using a linear model - considering only the joint motion around the neutral position - we need only one parameter (joint stiffness) to describe the characteristics of the joint. However, if the full feature of the joint-motion-moment is modeled using a nonlinear model, we will need more than three parameters. Therefore, the characteristics of the resistance in the joint are conventionally described by using "joint stiffness" in the clinical practice, implying that a linear model is applied to describe the joint-motion-moment around the neutral position. The linear joint stiffness applied in the current study reflects roughly the average joint stiffness as observed experimentally. Technically, it is possible to include the time-dependent and nonlinear characteristics of the joint stiffness into the finger model. However, in order to make it easier for the parametric studies, we felt it was more appropriate to assume a linear joint-motion-moment as in the current parametric study.
In summary, we theoretically analyzed the effects of the increased joint stiffness on muscle loading in a thumb in the current study. Our simulations indicated that the increase in joint stiffness induced substantial increases in muscle forces, especially in EPL and FPL muscles in response to the IP, MCP, or CMC extension/flexion motions. One of the potential applications of the proposed model is the estimation of the joint stiffness inversely using the test data of the gripping tests, which are conducted routinely in clinical diagnostics. The current simulation results might suggest that it is possible, theoretically, to improve the range of motion for OA patients in early stages by increasing muscle strength through exercise.
The findings and conclusions in this report are those of the authors and do not necessarily represent the views of the National Institute for Occupational Safety and Health.
- Chaisson CE, Zhang Y, Sharma L, Felson DT: Higher grip strength increases the risk of incident radiographic osteoarthritis in proximal hand joints. Osteoarthritis Cartilage 2000,8(Suppl A):S29–32.View ArticleGoogle Scholar
- Wright V, Johns RJ: Physical factors concerned with the stiffness of normal and diseased joints. Bull Johns Hopkins Hosp 1960, 106: 215–31.Google Scholar
- Altman R, Alarcon G, Appelrouth D, Bloch D, Borenstein D, Brandt K, Brown C, Cooke TD, Daniel W, Gray R, et al.: The American College of Rheumatology criteria for the classification and reporting of osteoarthritis of the hand. Arthritis Rheum 1990,33(11):1601–10. 10.1002/art.1780331101View ArticleGoogle Scholar
- Petersson IF, Jacobsson LT: Osteoarthritis of the peripheral joints. Best Pract Res Clin Rheumatol 2002,16(5):741–60. 10.1053/berh.2002.0266View ArticleGoogle Scholar
- Herzog W, Clark A, Longino D: Joint mechanics in osteoarthritis. Novartis Found Symp 2004, 260: 79–95. full_textView ArticleGoogle Scholar
- Sancho-Bru JL, Perez-Gonzalez A, Vergara-Monedero M, Giurintano D: A 3-D dynamic model of human finger for studying free movements. J Biomech 2001,34(11):1491–500. 10.1016/S0021-9290(01)00106-3View ArticleGoogle Scholar
- Sancho-Bru JL, Perez-Gonzalez A, Vergara M, Giurintano DJ: A 3D biomechanical model of the hand for power grip. J Biomech Eng 2003, 125: 78–83. 10.1115/1.1532791View ArticleGoogle Scholar
- Brook N, Mizrahi J, Shoham M, Dayan J: A biomechanical model of index finger dynamics. Med Eng Phys 1995, 17: 54–63. 10.1016/1350-4533(95)90378-OView ArticleGoogle Scholar
- Valero-Cuevas FJ, Johanson ME, Towles JD: Towards a realistic biomechanical model of the thumb: the choice of kinematic description may be more critical than the solution method or the variability/uncertainty of musculoskeletal parameters. J Biomech 2003,36(7):1019–30. 10.1016/S0021-9290(03)00061-7View ArticleGoogle Scholar
- Valero-Cuevas FJ: An integrative approach to the biomechanical function and neuromuscular control of the fingers. J Biomech 2005,38(4):673–84. 10.1016/j.jbiomech.2004.04.006View ArticleGoogle Scholar
- Imaeda T, Niebur G, Cooney WP, Linscheid RL, An KN: Kinematics of the normal trapeziometacarpal joint. J Orthop Res 1994,12(2):197–204. 10.1002/jor.1100120208View ArticleGoogle Scholar
- Li ZM, Tang J, Chakan M, Kaz R: Complex, multidimensional thumb movements generated by individual extrinsic muscles. J Orthop Res 2008,26(9):1289–95. 10.1002/jor.20641View ArticleGoogle Scholar
- An KN, Chao EY, Cooney WP, Linscheid RL: Normative model of human hand for biomechanical analysis. J Biomech 1979,12(10):775–88. 10.1016/0021-9290(79)90163-5View ArticleGoogle Scholar
- Wu JZ, An KN, Cutlip RG, Andrew ME, Dong RG: Modeling of the muscle/tendon excursions and moment arms in the thumb using the commercial software anybody. J Biomech 2009,42(3):383–8. 10.1016/j.jbiomech.2008.11.008View ArticleGoogle Scholar
- Smutz WP, Kongsayreepong A, Hughes RE, Niebur G, Cooney WP, An KN: Mechanical advantage of the thumb muscles. J Biomech 1998,31(6):565–70. 10.1016/S0021-9290(98)00043-8View ArticleGoogle Scholar
- Li ZM, Davis G, Gustafson NP, Goitz RJ: A robot-assisted study of intrinsic muscle regulation on proximal interphalangeal joint stiffness by varying metacarpophalangeal joint position. J Orthop Res 2006,24(3):407–15. 10.1002/jor.20046View ArticleGoogle Scholar
- Rasmussen J, Damsgaard M, Voigt M: Muscle recruitment by the min/max criterion - a comparative numerical study. J Biomech 2001,34(3):409–15. 10.1016/S0021-9290(00)00191-3View ArticleGoogle Scholar
- Epstein M, Herzog W: Theoretical models of skeletal muscle. 1st edition. Chichester-New York-Weinheim-Brisbance-Singapore-Toronto: John Wiley and Sons; 1998.Google Scholar
- Messier SP, Loeser RF, Hoover JL, Semble EL, Wise CM: Osteoarthritis of the knee: effects on gait, strength, and flexibility. Arch Phys Med Rehabil 1992, 73: 29–36.Google Scholar
- Esteki A, Mansour JM: An experimentally based nonlinear viscoelastic model of joint passive moment. J Biomech 1996,29(4):443–50. 10.1016/0021-9290(95)00081-XView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.