Next Article in Journal
Impact of Treadmill Interval Running on the Appearance of Zinc Finger Protein FHL2 in Bone Marrow Cells in a Rat Model: A Pilot Study
Previous Article in Journal
Long-Chain Polyunsaturated Fatty Acid Status at Birth and Development of Childhood Allergy: A Systematic Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of Optimization Strategies for Musculoskeletal Modeling of the Wrist for Therapy Planning in Case of Total Wrist Arthroplasty

1
Department of Orthopaedics, Trauma and Reconstructive Surgery, RWTH Aachen University Hospital, 52074 Aachen, Germany
2
Department of Orthopaedic Surgery, Klinikum Wels-Grieskirchen, 4600 Wels, Austria
3
Department of Orthopaedic and Trauma Surgery, Eifelklinik St. Brigida, 52152 Simmerath, Germany
*
Author to whom correspondence should be addressed.
Life 2022, 12(4), 527; https://doi.org/10.3390/life12040527
Submission received: 11 February 2022 / Revised: 29 March 2022 / Accepted: 30 March 2022 / Published: 2 April 2022
(This article belongs to the Section Medical Research)

Abstract

:
The human wrist joint is an elegant mechanism. The wrist allows the positioning and orienting of the hand to the forearm. The computational modeling of the human hand, especially of the wrist joint, can reveal important information about biomechanical mechanisms and provide the basis for its dysfunction and pathologies. For instance, this could be used for therapy planning in total wrist arthroplasty (TWA). In this study, different optimization methods and sensitivity analyses of anatomical parameters for musculoskeletal modeling were presented. Optimization includes finding the best available value of an objective function, including a variety of different types of objective functions. In the simplest case, optimization consists of maximizing or minimizing a function by systematically choosing input values from within an allowed set and computing the value of the function. Optimization techniques are used in many facets, such as the model building of joints or joint systems such as the wrist. The purpose of this study is to show the variability and influence of the included information for modeling, investigating the biomechanical function and load situation of the joint in representative scenarios. These possibilities to take them into account by an optimization and seem essential for the application of computational modeling to joint pathologies.

1. Introduction

The wrist is used in many ways and in many different situations in our daily lives; due to that, it is often involved in injuries, diseases, or deformations that may affect the quality of life. In diseases such as rheumatoid arthritis (RA), the wrist joint is frequently painful with impaired motion. Furthermore, it is involved in approximately 25% of all sports-related injuries [1]. It is necessary to objectively quantify parameters that affect wrist function to evaluate functional impairment and the efficacy of therapeutic interventions, such as total wrist arthroplasty (TWA), general wrist surgery, or rehabilitation approaches/programs following acute or chronic wrist ailments.
“Knowledge of the magnitudes and directions of the forces passing through the human joint is essential in the design of many surgical and therapeutic treatment programs” [2]. This missing information is still a current problem and has an interest in musculoskeletal/biomechanical modeling (MBM).
In the case of focusing on the characterization of the wrist joint, one important point is the investigation of the relationship between body motion and biomechanical loads. Here, the intention is to determine optimal posture and exercises during physiological activities (e.g., activities of daily living, ergonomics, training, and also rehabilitation). Additionally, this information could be used to design implants and to prevent, evaluate, and manage wrist joint disorders (e.g., joint deformities because of RA and surgical fixation strategies).
MBM can give insights into the physiological/mechanical behaviour of the joints. Subject-specific biomechanical models are essential for future clinical applications, such as predicting the effects of orthopedic surgery. MBM and simulations may be useful tools to solve complex biomechanical problems, revolutionizing medical decision-making and surgical management strategies. The simulation of the wrist joint is integral to understanding the influence of soft tissue. For example, how muscles behave to move the skeleton in normal and pathological states [3]. Delp et al. stated that a biomechanical description of the wrist’s surrounding muscles is needed to help design, e.g., tendon transfer surgeries [4]. The MBM of the wrist can reveal facts about its biomechanical mechanisms and load application, providing the basis for investigating dysfunction [5].
An open issue in MBM remains to solve the muscular redundancy problem [6]. If the diversity of musculoskeletal function is considered, several different criteria for muscle selection may likely be utilized for different possible activities [7]. All linear or nonlinear optimization procedures assume that the body selects muscles for a given activity according to some criteria (e.g., minimization of muscle force) [7]. The potential of multi-objective optimization, when both musculo-tendon forces and joint reaction forces are minimized, is still unclear. Considering more than one objective function in the optimization problem introduces additional degrees of freedom (DOF) that can stand for the compromise between musculo-tendon and joint reaction forces. The choice of the objective function used to solve the muscular redundancy problem is the bottleneck.
Besides this problem, biomechanical models of the wrist/forearm system may provide an alternative approach to estimating the relative changes in strength and the applied load to the joint [8]. Building a complete biomechanical model of the wrist is very complicated because of its anatomy, including muscles, and the high number of DOF [5]. Sometimes, this complex approach can be very time-consuming.
This study aims to introduce, based on a simple planar model of the wrist joint, the ensemble of optimal solutions and the solution of a global objective method for the following emblematic three-objective problem: optimizing (and minimizing) musculo-tendon forces in a static situation. In this study, we propose possibilities for performing the concurrent muscle problem. We provide different illustrative examples to demonstrate the implementation of the optimization approach to the wrist. Furthermore, different anatomical parameters and their sensitivity in the case of modeling were presented and their possible influence on the results was demonstrated. The aim of the sensitivity analysis was to investigate inaccuracies in muscle modeling parameters and the development of concepts to reduce related error sources.

2. Material and Methods

2.1. Sensitivity Analysis

2.1.1. General Background of a Sensitivity Analysis

The sensitivity analysis (SA) evaluates the uncertainty in the output of a mathematical model. It can be divided and allocated into different sources of uncertainty in its inputs. However, the changing, e.g., the quality of X-ray images or other (imaging/geometrical) information, and the subjectivity of the included data; furthermore, the subjectivity of the users mandatorily of medical experts lead to inaccuracies in e.g., landmark identification. The extent and effect of these inaccuracies on the MBM output are of interest concerning the robustness of the idealization approach itself and the reliability in the case of the performed therapy planning.
For this approach, MATLAB (R2020a, Mathworks) was used, a well-established numerical computing environment.

2.1.2. Determination of Anatomical and Modeling Information

MBM of musculoskeletal structures requires accurate data on anatomical parameters e.g., anatomical landmarks, muscle lengths (MLs), moment arms (MAs), and those parameters describing the segment position.
The main functions of the wrist are to provide mobility and sustain the load. It is reasonable to assume that mechanical factors play an important role in the function of the wrist joint. To determine these mechanical factors, especially muscle loads, motions, and external loads have to be quantitated in precise mechanical terms. Therefore, the knowledge of anatomical and relevant input parameters is important for our example, which focuses primarily on muscles.

2.1.3. Physiological Cross-Sectional Area

The physiological cross-sectional area (PCSA) of a muscle is a measure of the number of sarcomeres in parallel with the angle of pull of the muscles [9]. In pennate muscles, the fibers act at an angle from the long axis. The angle between the orientation of the muscle fiber and the angle of the muscle’s line of action is called the pennation angle [9]. For each muscle with a parallel-fibered structure, the PCSA can be calculated by the following equation:
PCSA   ( mm 2 ) = mass   ( g ) density   ( g mm 3 ) · optimum   fiber   length   ( mm )
In pennate muscles, the PCSA becomes:
PCSA   ( mm 2 ) = mass   ( g ) · cos θ ( degree ) density   ( g mm 3 ) · optimum   fiber   length   ( mm )
where θ is the pennation angle that increases as the muscle shortens. As in the previous work of Borst et al. [10], a muscle density of 0.0010567 g mm 3 was assumed [5]. The PCSA was taken from the literature.

2.1.4. Peak Force Calculation (PFC)

The peak force is calculated as the product of PCSA and the specific tension σ as follows:
PFC   ( N ) = PCSA ( mm 2 ) · σ ( N mm 2 )
A wide range of tension values for skeletal muscles has been reported, e.g., [9,11]. Most of these tension values were measured during isometric conditions and ranged from 0.2 to 1.0 N mm 2 . These higher values were recorded in pennate muscles, which are those whose fibers lie at an angle from the main axis of the muscle [9]. Such an orientation effectively increases the PCSA above that measured and used in the stress calculation [9].

2.2. General Approach of MBM

MBM is a noninvasive way to investigate the relationship between body motion and internal biomechanical parameters in a wide range of conditions. Validated models reduce the need for experimental studies, providing useful information about joint kinematics, joint contact pressures/forces, soft tissue tensions, and range of motion (ROM) [12]. Knowledge of muscle and joint contact forces would be highly beneficial from a clinical perspective, but their direct measurement is not feasible in a clinical setting [13]. MBM potentially overcomes the limitations of current practice, and the treatments are based on subjective, static, and mostly qualitative assessments [14].
In the context of the wrist joint, this technique is a valuable approach to examine the effect of changes in muscle force on loads and the risk of wrist pain. Changes in muscle forces have not been properly simulated by previous models. The estimation of general and especially individual muscle forces and their contribution to joint moments is essential for understanding how humans solve the dynamics of movement [15].
We performed a static study for the (mechanical) system composed of the forearm and wrist, assumed with negligible weight (Figure 1).
The wrist is approximated as a perfect hinge joint for flexion and extension and radial and ulnar deviation according to the state of the art [16]. Six muscles realize the following main motion: M. flexor carpi radialis (FCR), M. palmaris longus (PL), M. flexor carpi ulnaris (FCU), M. extensor carpi ulnaris (ECU), Mm. extensors carpi radialis longus et brevis (ECRL and ECRB) [16,17,18,19,20]. The flexion is mainly produced by the PL, FCU, and FCR. The motion is supported by the flexor digitorum superficialis muscle. The extension is mainly produced by the ECRL, ECRB, and ECU. They will be assisted from the extensor digitorum muscle. The abduction muscles of the wrist (radial deviation) are the FCR, the ECRL, and the ECRB, and the abductor pollicis longus supports the motion. Muscles involved in the adduction of the wrist (ulnar deviation) are the FCU and the ECU.
For an exemplar computation, a maximum load of 100 N was applied to the wrist with a lever arm of 10 cm (=0.1 m), which resulted in a moment of 10 N·m (=1000 N·cm) (Figure 2). To hold the joint in static equilibrium in a linear case, the forces for each muscle must be calculated e.g., for flexion. That means, to ensure an equilibrium the extensor muscles ECU, ECRB, and ECRL must hold the hand with 10 Nm.
What makes this system indeterminate is the fact that there exist more unknowns than solvable equations. This means the system is (statically) indeterminate. There are many possible combinations of values for the muscle forces that will satisfy the equilibrium equation. The study and treatment of these conditions could greatly benefit from combined software tools that offer better insights into neuromuscular biomechanics, and predictive capabilities for optimal surgical treatment planning [21].
An open issue in MBM is solving the muscular redundancy problem [6]. This problem states that there are multiple ways to perform a movement with different muscle activations. In other words, having more muscles than DOFs [22]. Under normal circumstances, no simple one-to-one correspondence exists between a muscular redundancy problem and a solution to the problem.
There are the following two widely used approaches to determining a unique solution: (1) reduction methods and (2) optimization methods. The reduction methods try to reduce the number of unknowns via the introduction of additional boundary constraints to obtain an equal number of equations. For example, if two muscles have the same attachment and perform a similar task, they could be seen as one muscle. This reduces the number of unknown muscle forces/parameters by one. The attractiveness of the reduction method is that the same technique used to solve the muscle force problem could be used.

2.2.1. Calculating Static Muscle Forces Using a General Optimization Approach Based on the Weighted Sum Method

An alternative to the reduction methods are optimization methods. One common approach to the force distribution problem is to formulate a static optimization model.
First, a weighted sum method is proposed. In this approach the following percentual part of the muscles for the flexion motion of the wrist in proportion to the main muscle cross-sectional was considered for the calculation:
M = i = 1 m ( d i · A i A F )
F 0
A = i = 1 m A i = 1
where F denotes the unknown muscle forces of ith forearm muscle, Ai denotes the cross-sectional area of the ith muscle, M the wrist joint moment, di the moment arm of the ith muscle of the wrist, m is the number of muscles involved in the motion (m = 3), and A is the main/the summation cross-sectional area of all involved muscles (A = 1 = 100%).

2.2.2. Calculating Static Muscle Forces Using Optimization Based on an Approach by Crowninshield and Brand

The approach proposed by Crowninshield and Brand is surely the most used optimization model in the literature [7]. Here, the cost function takes the form of a weighted sum of muscle forces raised to some non-linear power. The weighting is based on constant values for the PCSA. The basis is that every physical system reaches its equilibrium state by the energy minimization principle. Transferred to the body this means that the muscles are activated in such a way that fatigue is minimized.
The principle is to minimize the so-called objective function U [7]:
Minimize   { U }
The objective function U proposed by Crowninshield and Brand is given as follows [7]:
U = i = 1 m ( F i PCSA i ) n
subject to:
M i = 1 m ( d i · F i ) = 0 M = i = 1 m ( d i · F i ) F i 0 i = 1   to   m
where Fi denotes the unknown muscle forces of ith forearm muscle, PCSAi denotes the physiological cross-sectional area of the ith muscle, M the wrist joint moment, di the moment arm vector of the ith muscle at the wrist, m is the number of muscles crossing the joint, and n is the exponent in the muscle endurance-force relationship. This norm is minimized when the radicand (summation of the nth power of muscle stresses) is minimized.
Depending on the value of n, the Crowninshield and Brand criteria may yield a linear programming problem (n = 1), a quadratic programming problem (n = 2), or a general nonlinear programming problem (n > 2).
Considering the optimization formulation applied to the indeterminate problem of an applied force of 100 N in case of flexion to the hand and the distance to the rotational center is D = 10 cm (Figure 2). For this approach, MATLAB (R2020a, Mathworks) was used. It includes hundreds of general-purpose and specialized functions for numerical computation, organized into libraries called toolboxes that cover such areas as simulation, optimization, systems, and control and signal processing. If Equation (9) is taken and insert the resulting moment, it follows:
1000 i = 1 m ( d i · F i ) = 0 1000 = i = 1 m ( d i · F i ) F i 0 i = 1   to   3
In the case of following Crowninshield and Brand and assuming n = 1 and m = 3, the following optimization problem will be originated [7]:
Min ( F ECU   F ECRB   F ECRL ) = i = 1 3 ( F i PCSA i ) 1 = ( F ECU PCSA ECU ) + ( F ECRB PCSA ECRB ) + ( F ECRL PCSA ECRL )
subject to:
M = i = 1 3 ( d i · F i ) = ( d ECU · F ECU ) + ( d ECRB · F ECRB ) + ( d ECRL · F ECRL )
In addition,
{ A eq = b eq A ineq 0
The problem can be classified as a linear problem. The objective function and the constraints are linear functions of the minimizing variables FECU, FECRB, and FECRL. To utilize MATLAB for solving the linear programming problem, the following “linprog” solver was used:
[ 1 PCSA ECU 1 PCSA ECRB 1 PCSA ECRL ] T [ F ECU F ECRB F ECRL ]
subject to:
[ d ECU d ECRB d ECRL ] [ F ECU F ECRB F ECRL ] = M = 1000   N · cm
where:
A eq = [ d ECU d ECRB d ECRL ] x = [ F ECU F ECRB F ECRL ] b eq = M
If we assume n = 2 and m = 3, the following optimization problem will be originated:
Min ( F ECU   F ECRB   F ECRL ) = i = 1 3 ( F i PCSA i ) 2 = ( F ECU PCSA ECU ) 2 + ( F ECRB PCSA ECRB ) 2 + ( F ECRL PCSA ECRL ) 2
subject to:
M = i = 1 3 ( d i · F i ) = ( d ECU · F ECU ) + ( d ECRB · F ECRB ) + ( d ECRL · F ECRL )
For the case of a quadratic approach it follows:
1 2 [ F ECU   F ECRB   F ECRL ] T [ 2 PCSA ECU 2 0 0 0 2 PCSA ECRB 2 0 0 0 2 PCSA ECRL 2 ] [ F ECU F ECRB F ECRL ]
where:
x T = 1 2 [ F ECU   F ECRB   F ECRL ] T H = [ 2 PCSA ECU 2 0 0 0 2 PCSA ECRB 2 0 0 0 2 PCSA ECRL 2 ] x = [ F ECU F ECRB F ECRL ]
subject to:
[ d ECU d ECRB d ECRL ] [ F ECU F ECRB F ECRL ] = 1000   N · cm
In the case where n > 2, e.g., if n = 3, the optimization problem may be stated as follows:
Min ( F ECU   F ECRB   F ECRL ) = i = 1 3 ( F i PCSA i ) 3 = ( F ECU PCSA ECU ) 3 + ( F ECRB PCSA ECRB ) 3 + ( F ECRL PCSA ECRL ) 3
subject to the following:
M = i = 1 3 ( d i · F i ) = ( d ECU ·   F ECU ) + ( d ECRB ·   F ECRB ) + ( d ECRL ·   F ECRL )
Now it is a general nonlinear programming problem. The reason for that is the nonlinearity of the objective function. To solve this problem, the “fmincon” solver in MATLAB was used.

2.2.3. Calculating Static Muscle Forces Using an Individual Optimization Approach

An alternative optimization approach was published by An et al. [2]. Here, a linear approach was presented with additional constraints that only a few muscles will be active during a special movement.
Min   σ
Subject to the following:
F i PCSA i σ
in which PCSA represents the physiological cross-sectional area of the ith muscle, and the left-hand side of the equation represented the “stress” (σi) of the ith muscle.
The resultant joint forces and moments were then obtained by the remaining force and moment equilibrium equations. The flexion-extension moment equations are as follows:
d i · F i = M
in which M represented the intersegmental moments due to externally applied load in flexion-extension about the joint center. The moment arms in flexion-extension of the ith muscle were di. Fi represents the magnitude of the ith unknown muscle force.
An et al. justified their method based on the fact that nonlinear programming methods can be very sensitive to ill-conditioned data [2]. In MATLAB, this type of problem will be solved using the “fminimax” solver.

3. Results

3.1. Determination of Anatomical Parameters for Analyzing the Wrist Joint

3.1.1. Anatomical Parameters

The different anatomical parameters are shown in Table 1, which are used for analyzing the muscles strength and loads of the wrist joint.

3.1.2. Moment Arms of the Muscles

Figure 3 shows the technical information for the different muscles during wrist motion.

3.2. Sensitivity Analysis

Table 2 shows the parameters for the sensitivity analysis.
Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8 show the results of the sensitivity analysis of the wrist joint.

3.3. Optimization of the Muscle Redundancy Problem

3.3.1. Results for Using a General Optimization Approach Based on the Weighted Sum Method

If we know to take the mean lever arms for each muscle, e.g., from [31], we will obtain for ECU a lever arm of 6 mm (0.006 m), for ECRB a lever arm of 12 mm (0.012 m), and for ECRL a lever arm of 7 mm (0.007 m). To calculate the part of the force for each muscle, we can choose the PCSA and calculate the percentual part of the summary of the complete cross-sectional area.
For example, if we take the PCSA from [23], we find an area for the ECU of 480 mm2 (=4.8 cm2), for an ECRB of 450 mm2 (=4.5 cm2), and for ECRL of 370 mm2 (=3.7 cm2). The sum of all is then 1300 mm2 (=13 cm2 = 100%). That means the percentual part for the ECU is 36.9%, for the ECRB it is 34.6%, and for the ECRL it is 28.5%.
We obtained the following:
1000   N · cm = 0.6   cm · ECU + 1.2   cm · ECRB + 0.7   cm · ECRL
With the substitution we obtain the following:
1000   N · cm = 0.6   cm · 0.37   F + 1.2   cm · 0.35   F + 0.7   cm · 0.28   F
That means the following:
1000   N · cm = 0.6   cm · 0.37   F + 1.2   cm · 0.35   F + 0.7   cm · 0.28   F 1000   N · cm = 0.22   cm · F + 0.42   cm · F + 0.20   cm · F 1000   N · cm = 0.84   cm · F 1190.48   N = F
That means we have to ensure a muscle force in the sum of 1190.48 N, and for each muscle and a force of the following:
ECU = 440.48 N (=0.37 F)
ECRB = 416.67 N (=0.35 F)
ECRL = 333.33 N (=0.28 F)

3.3.2. Results for Calculating Static Muscle Forces Using Optimization Based on an Approach by Crowninshield and Brand

As these solutions are optimal in the sense that they minimize the utility function U as follows:
For n = 1:
ECU = 0 N
ECRB = 833.34 N
ECRL = 0 N
For n = 2:
ECU = 313.03 N
ECRB = 550.24 N
ECRL = 216.99 N
For n = 3:
ECU = 377.05 N
ECRB = 484.03 N
ECRL = 275.62 N
This example illustrates a concept in biomechanics that muscles generally must generate large forces to balance external loads because of their often very short moment arms. Crowninshield and Brand discussed how, in the case the exponent in the objective function increases, the resulting force is more distributed between the movement-involved muscles to reduce the onset of muscle fatigue [7].

3.3.3. Calculating Static Muscle Forces Using Optimization Based on an Approach by An et al.

As these solutions are optimal in the sense that they minimize the utility function U:
ECU = 441.58 N
ECRB = 413.98 N
ECRL = 340.39 N

4. Discussion

The carpus is a challenging and complex mechanical system because of its organization of bones and soft tissue structures. The aim of this in silico investigation is to put forward the idea of multi-objective optimization in MBM and the understanding of the wrist by the recognition that there is no standard carpal mechanical system and no typical wrist. This study investigated the effect of variations in muscle forces acting on the wrist using a simple MBM approach.

4.1. Sensitivity Analysis

Indeed, the method proposed in this paper allows comparative statements and is based on information taken out of the literature (Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8). However, differences in the impact of inaccuracies in the modeling parameters on the simulation results could be detected and should be further investigated. They showed a huge variation and could influence in the case of choosing them in a model environment. The extent of scattering observed verified the necessity to perform a sensitivity analysis of the available anatomical information. The strong scattering originates from the individuals who were investigated in the different references, e.g., because of the stress values, which range from 0.2 to 1.0 N mm 2 .

4.2. Optimization

We tried to introduce the different optimization approaches and the idea of MBM via an easily simplified model of the wrist because it is easier to analyze and interpret the effects of different optimizations on a simple model before moving to complex ones [6]. In this case, we are interested in examining the muscle forces using optimization methods; we have to consider some assumptions. The main one is that the solution to the problem is unique to each individual task. The more realistic the approximation is based on physiological or functional parameters, the better the estimation of realistic results [32].
The outcomes indicate that muscle forces had a marked impact on computed wrist loads. The knowledge of the forces passing through the wrist joint is essential in the design of many surgical and therapeutic treatments because normal carpal mechanics relies on the interactions between the muscles, ligaments, and carpal bone morphology. The novel feature of this work is the introduction of different optimization approaches that efficiently emulate the behavior of static load analyses within the optimization’s constraints. The example applications yielded reasonable load predictions. To obtain these results, we had to adjust the optimization constraints. The optimization function is one of the most important points in the optimization process because its minimization or maximization determines the solution to the problem [32]. The objective function must be related to the physical problem, e.g., the physiological and pathological situations, respectively.
We also made sure that the different optimization approaches were able to match the (experimentally/literature-based) obtained loads.
For example, when it was planned to include such an approach in the daily clinical process, the fmincon algorithm from the MATLAB Optimization Toolbox was simple to implement and very fast to be of much practical value. If a computer cluster is available, e.g., fmincon could be run in parallel using the MATLAB Parallel Computing Toolbox [33]. That implies that such an approach produces additional information without elongating the daily clinical workflow.
Contact surfaces and ligaments are not explicitly considered in this investigation. The interactions between the muscular and joint structures, and the possibility to take them into account in MBM, seem essential for applications to joint pathologies [6].
The next step should be that such multi-objective optimization approaches should be evaluated on more complex musculoskeletal models and, in the best case, validated by the results of an instrumented prosthesis. No data for the wrist joint in the context of the instrumented prosthesis have been reported so far in the literature.

4.3. Limitation

The replication of wrist motions in vitro using only six tendons was a limitation of this study. These tendons were selected since they insert at the base of the metacarpals, thereby primarily affecting the wrist. All muscles contributing to wrist torque, including the extrinsic muscles of the fingers and thumb, should be considered for future MBM approaches. The inclusion of one type of objective function, one joint without an angle and in the neutral position, respectively, and one external moment does not fully represent reality.
Because of the consideration of different optimization criteria, a comparison of the criteria and the resulting outcome is difficult because of the following: (1) the possible different formulation of the problem; (2) the different model assumptions used; (3) the different accuracy and reliability of the input data (e.g., PCSA, error in measurement data); (4) the different optimization algorithms used; (5) the use of different degrees of freedom according to the complexity of the model [32].
Furthermore, the planar model of the wrist used in this study seemed appropriate for an exploratory study of multi-objective optimization but is not appropriate for validation purposes. Therefore, future work should preferentially consist of multi-objective optimizations applied to more complex models of the wrist.

5. Conclusions

A comprehensive, unifying model of the functional kinematics of the carpus remains elusive. The complexity involved both in determining the input parameters and the interaction between the muscles involved is non-trivial.
The MATLAB options provide a powerful and flexible approach for generating optimal simulations of musculoskeletal parameters using the implemented optimization approaches. This should facilitate the use of these methods in developing therapies and devices for clinical conditions. The approach presented here shows that numerical simulations play an important role in solving complex engineering problems. Moreover, they have the potential to revolutionize medical decision-making and treatment strategies.

Author Contributions

Conceptualization, J.E., J.L., V.Q., B.R., M.P., F.H. and F.M.; methodology, J.E., J.L., V.Q., B.R., M.P., F.H. and F.M.; resources, J.E., J.L., V.Q., B.R., M.P., F.H. and F.M.; data curation, J.E. and M.P.; writing—original draft preparation, J.E.; writing—review and editing, J.E., J.L., V.Q., B.R., M.P., F.H. and F.M.; visualization, J.E.; supervision, F.M.; project administration, J.E. and F.M.; funding acquisition, J.E. All authors have read and agreed to the published version of the manuscript.

Funding

The research was partly funded by the German Research Foundation (Deutsche Forschungsgemeinschaft (DFG)) ES 442/1-1 and RA 2187/4-1, and partly funded by “NRW-Patent-Validierung”-EFRE-0400355.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no potential conflict of interest concerning the research, authorship, and/or publication of this article.

References

  1. Avery, D.M.; Rodner, C.M.; Edgar, C.M. Sports-related wrist and hand injuries: A review. J. Orthop. Surg. Res. 2016, 11, 99. [Google Scholar] [CrossRef] [Green Version]
  2. An, K.N.; Kwak, B.M.; Chao, E.Y.; Morrey, B.F. Determination of muscle and joint forces: A new technique to solve the indeterminate problem. J. Biomech. Eng. 1984, 106, 364–367. [Google Scholar] [CrossRef]
  3. Li, Z.; Mogk, J.P.M.; Lee, D.; Bibliowicz, J.; Agur, A.M. Development of an architecturally comprehensive database of forearm flexors and extensors from a single cadaveric specimen. Comput. Methods Biomech. Biomed. Eng. Imaging Vis. 2015, 3, 3–12. [Google Scholar] [CrossRef]
  4. Delp, S.L.; Grierson, A.E.; Buchanan, T.S. Maximumisometric moments generated by the wrist muscles in flexion-extension and radial-ulnar deviation. J. Biomech. 1996, 29, 1371–1375. [Google Scholar] [CrossRef]
  5. Mirakhorlo, M.; Visser, J.M.A.; Goislard de Monsabert, B.; van der Helm, F.C.; Maas, H.; Veeger, H.E. Anatomical parameters for musculoskeletal modeling of the hand and wrist. Int. Biomech. 2016, 3, 40–49. [Google Scholar] [CrossRef]
  6. Dumas, R.; Moissenet, F.; Lafon, Y.; Cheze, L. Multi-objective optimisation for musculoskeletal modelling: Application to a planar elbow model. Proc. Inst. Mech. Eng. Part H J. Eng. Med. 2014, 228, 1108–1113. [Google Scholar] [CrossRef]
  7. Crowninshield, R.D.; Brand, R.A. A physiologically based criterion of muscle force prediction in locomotion. J. Biomech. 1981, 14, 793–801. [Google Scholar] [CrossRef]
  8. La Delfa, N.J.; Potvin, J.R. A musculoskeletal model to estimate the relative changes in wrist strength due to interacting wrist and forearm postures. Comput. Methods Biomech. Biomed. Eng. 2017, 20, 1403–1411. [Google Scholar] [CrossRef]
  9. Winter, D.A. Biomechanics and Motor Control of Human Movement, 4th ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2009; ISBN 978-0-470-39818-0. [Google Scholar]
  10. Borst, J.; Forbes, P.A.; Happee, R.; Veeger, D.H. Muscle parameters for musculoskeletal modelling of the human neck. Clin. Biomech. 2011, 26, 343–351. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Holzbaur, K.R.S.; Murray, W.M.; Delp, S.L. A model of the upper extremity for simulating musculoskeletal surgery and analyzing neuromuscular control. Ann. Biomed. Eng. 2005, 33, 829–840. [Google Scholar] [CrossRef] [PubMed]
  12. Majors, B.J.; Wayne, J.S. Development and validation of a computational model for investigation of wrist biomechanics. Ann. Biomed. Eng. 2011, 39, 2807. [Google Scholar] [CrossRef] [PubMed]
  13. Marra, M.A.; Vanheule, V.; Fluit, R.; Koopman, B.H.; Rasmussen, J.; Verdonschot, N.; Andersen, M.S. A subject-specific musculoskeletal modeling framework to predict in vivo mechanics of total knee arthroplasty. J. Biomech. Eng. 2015, 137, 020904. [Google Scholar] [CrossRef] [PubMed]
  14. Fregly, B.J.; Besier, T.F.; Lloyd, D.G.; Delp, S.L.; Banks, S.A.; Pandy, M.G.; D’lima, D.D. Grand challenge competition to predict in vivo knee loads. J. Orthop. Res. 2012, 30, 503–513. [Google Scholar] [CrossRef] [Green Version]
  15. Pizzolato, C.; Lloyd, D.G.; Sartori, M.; Ceseracciu, E.; Besier, T.F.; Fregly, B.J.; Reggiani, M. CEINMS: A toolbox to investigate the influence of different neural control solutions on the prediction of muscle excitation and joint moments during dynamic motor tasks. J. Biomech. 2015, 48, 3929–3936. [Google Scholar] [CrossRef] [PubMed]
  16. Eschweiler, J.; Li, J.; Quack, V.; Rath, B.; Baroncini, A.; Hildebrand, F.; Migliorini, F. Anatomy, Biomechanics, and Loads of the Wrist Joint. Life 2022, 12, 188. [Google Scholar] [CrossRef] [PubMed]
  17. Berger, R.A. The Anatomy and Basic Biomechanics of the Wrist joint. J. Hand Ther. 1996, 9, 84–93. [Google Scholar] [CrossRef]
  18. Taleisnik, J. The Wrist; Churchill Livingstone: New York, NY, USA, 1985; ISBN 9780443081347. [Google Scholar]
  19. Schmitt, R. Funktionelle Anatomie und Biomechanik des Karpus. Radiologe 2006, 46, 638–648. [Google Scholar] [CrossRef]
  20. Kijima, Y.; Viegas, S.F. Wrist anatomy and biomechanics. J. Hand Surg. 2009, 34, 1555–1563. [Google Scholar] [CrossRef] [PubMed]
  21. Mansouri, M.; Reinbolt, J.A. A platform for dynamic simulation and control of movement based on OpenSim and MATLAB. J. Biomech. 2012, 45, 1517–1521. [Google Scholar] [CrossRef] [Green Version]
  22. Kutch, J.J.; Valero-Cuevas, F.J. Muscle redundancy does not imply robustness to muscle dysfunction. J. Biomech. 2011, 44, 1264–1270. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. An, K.-N.; Hui, F.C.; Morrey, B.F.; Linscheid, R.L.; Chao, E.Y. Muscles across the elbow joint: A biomechanical analysis. J. Biomech. 1981, 14, 659–669. [Google Scholar] [CrossRef]
  24. Brand, P.W.; Beach, R.B.; Thompson, D.E. Relative tension and potential excursion of muscles in the forearm and hand. J. Hand Surg. 1981, 6, 209–219. [Google Scholar] [CrossRef]
  25. Cutts, A.; Alexander, R.M.; Ker, R.F. Ratios of cross-sectional areas of muscles and their tendons in a healthy human forearm. J. Anat. 1991, 176, 133. [Google Scholar] [PubMed]
  26. Lemay, M.A.; Crago, P.E. A dynamic model for simulating movements of the elbow, forearm, and wrist. J. Biomech. 1996, 29, 1319–1330. [Google Scholar] [CrossRef]
  27. Lieber, R.L.; Jacobson, M.D.; Fazeli, B.M.; Abrams, R.A.; Botte, M.J. Architecture of selected muscles of the arm and forearm: Anatomy and implications for tendon transfer. J. Hand Surg. 1992, 17, 787–798. [Google Scholar] [CrossRef]
  28. Loren, G.J.; Shoemaker, S.D.; Burkholder, T.J.; Jacobson, M.D.; Fridén, J.; Lieber, R.L. Human wrist motors: Biomechanical design and application to tendon transfers. J. Biomech. 1996, 29, 331–342. [Google Scholar] [CrossRef]
  29. Garland, A.K.; Shah, D.S.; Kedgley, A.E. Wrist tendon moment arms: Quantification by imaging and experimental techniques. J. Biomech. 2018, 68, 136–140. [Google Scholar] [CrossRef]
  30. Brand, P.W. Mechanics of individual muscles at individual joints. Clin. Mech. Hand 1985, 278–311. [Google Scholar]
  31. Horii, E.; An, K.N.; Linscheid, R.L. Excursion of prime wrist tendons. J. Hand Surg. 1993, 18, 83–90. [Google Scholar] [CrossRef]
  32. Tsirakos, D.; Baltzopoulos, V.; Bartlett, R. Inverse optimization’ Functional and physiological considerations related to the force-sharing problem. Crit. Rev. Biomed. Eng. 1997, 25, 371–407. [Google Scholar] [CrossRef]
  33. Lee, L.-F.; Umberger, B.R. Generating optimal control simulations of musculoskeletal movement using OpenSim and MATLAB. PeerJ 2016, 4, e1638. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Muscles direction and wrist joint moment arms: The force direction is parallel to the forearm so that the physiological cross-sectional area of the muscle is the leading parameter. (1): M. flexor carpi radialis, (2): M. palmaris longus, (3): M. flexor carpi ulnaris, (4): M. extensor carpi ulnaris, (5–6): Mm. extensors carpi radialis longus et brevis; U: Ulna; R: Radius; M: Metacarpal bone; C: Capitatum; L: Lunatum.
Figure 1. Muscles direction and wrist joint moment arms: The force direction is parallel to the forearm so that the physiological cross-sectional area of the muscle is the leading parameter. (1): M. flexor carpi radialis, (2): M. palmaris longus, (3): M. flexor carpi ulnaris, (4): M. extensor carpi ulnaris, (5–6): Mm. extensors carpi radialis longus et brevis; U: Ulna; R: Radius; M: Metacarpal bone; C: Capitatum; L: Lunatum.
Life 12 00527 g001
Figure 2. Free body diagram: (A): lever arm of ECRB; (B–C): lever arm of ECU and ECRL; (D): lever arm (10 cm = 0.1 m) of an applied load of 100 N.
Figure 2. Free body diagram: (A): lever arm of ECRB; (B–C): lever arm of ECU and ECRL; (D): lever arm (10 cm = 0.1 m) of an applied load of 100 N.
Life 12 00527 g002
Figure 3. Muscle moment arms: information taken from [11] are shown as black diamonds, the parameters shown as blue quadrat were taken from [28], the parameters shown as green quadrat were taken from [29], the parameters shown as red triangle were taken from [30], the parameters shown as red diamonds were taken from [31]. (M. flexor carpi radialis (FCR), M. palmaris longus (PL), M. flexor carpi ulnaris (FCU), M. extensor carpi ulnaris (ECU), Mm. extensors carpi radialis longus et brevis (ECRL and ECRB)).
Figure 3. Muscle moment arms: information taken from [11] are shown as black diamonds, the parameters shown as blue quadrat were taken from [28], the parameters shown as green quadrat were taken from [29], the parameters shown as red triangle were taken from [30], the parameters shown as red diamonds were taken from [31]. (M. flexor carpi radialis (FCR), M. palmaris longus (PL), M. flexor carpi ulnaris (FCU), M. extensor carpi ulnaris (ECU), Mm. extensors carpi radialis longus et brevis (ECRL and ECRB)).
Life 12 00527 g003
Figure 4. Resulting wrist muscle force in case of varied different modeling parameters. (M. flexor carpi radialis (FCR), M. palmaris longus (PL), M. flexor carpi ulnaris (FCU), M. extensor carpi ulnaris (ECU), Mm. extensors carpi radialis longus et brevis (ECRL and ECRB); PCSA: physiological cross-sectional area; PFC: Peak force calculation).
Figure 4. Resulting wrist muscle force in case of varied different modeling parameters. (M. flexor carpi radialis (FCR), M. palmaris longus (PL), M. flexor carpi ulnaris (FCU), M. extensor carpi ulnaris (ECU), Mm. extensors carpi radialis longus et brevis (ECRL and ECRB); PCSA: physiological cross-sectional area; PFC: Peak force calculation).
Life 12 00527 g004
Figure 5. The figure shows the resulting wrist moment during extension in case for different parameters (muscle forces and moment arm). (M. extensor carpi ulnaris (ECU), Mm. extensors carpi radialis longus et brevis (ECRL and ECRB)).
Figure 5. The figure shows the resulting wrist moment during extension in case for different parameters (muscle forces and moment arm). (M. extensor carpi ulnaris (ECU), Mm. extensors carpi radialis longus et brevis (ECRL and ECRB)).
Life 12 00527 g005
Figure 6. The figure shows the resulting wrist moment during flexion for varying different parameters (muscle forces and moment arms). (M. flexor carpi radialis (FCR), M. palmaris longus (PL), M. flexor carpi ulnaris (FCU)).
Figure 6. The figure shows the resulting wrist moment during flexion for varying different parameters (muscle forces and moment arms). (M. flexor carpi radialis (FCR), M. palmaris longus (PL), M. flexor carpi ulnaris (FCU)).
Life 12 00527 g006
Figure 7. The figure shows the resulting wrist moment during radial deviation for different parameters. The adduction is produced by the ECU and FCU. The abduction is produced by the FCR, ECRL, and ECRB. (M. flexor carpi radialis (FCR), M. palmaris longus (PL), Mm. extensors carpi radialis longus et brevis (ECRL and ECRB)).
Figure 7. The figure shows the resulting wrist moment during radial deviation for different parameters. The adduction is produced by the ECU and FCU. The abduction is produced by the FCR, ECRL, and ECRB. (M. flexor carpi radialis (FCR), M. palmaris longus (PL), Mm. extensors carpi radialis longus et brevis (ECRL and ECRB)).
Life 12 00527 g007
Figure 8. The figure shows the resulting wrist moment during ulnar deviation for different parameters. (M. flexor carpi ulnaris (FCU), M. extensor carpi ulnaris (ECU)).
Figure 8. The figure shows the resulting wrist moment during ulnar deviation for different parameters. (M. flexor carpi ulnaris (FCU), M. extensor carpi ulnaris (ECU)).
Life 12 00527 g008
Table 1. The presented information was taken out of the literature. There exists a huge variation of the detected data.
Table 1. The presented information was taken out of the literature. There exists a huge variation of the detected data.
MuscleVolumeSDFiber LengthSDTendon LengthSDPennation AngleSDMassSDPCSASDReferencen
[cm3] [mm] [mm] [deg] [g] [mm2]
FCR0.000.0076.400.00188.000.000.000.007.400.0091.600.00[5]1
12.404.9058.009.000.000.000.000.000.000.00200.0060.00[23]6
0.000.0052.000.000.000.000.000.000.000.000.000.00[24]15
0.000.0049.000.000.000.000.000.0035.000.00681.000.00[25]1
0.000.0063.000.00244.000.003.000.000.000.00160.000.00[11]*
0.000.00164.000.000.000.003.100.000.000.00199.000.00[26]*
0.000.000.000.000.000.000.000.000.000.000.000.00[3]1
0.000.000.000.000.000.000.000.000.000.00199.000.00[27]18
PL0.000.0089.200.00247.000.000.000.003.400.0036.100.00[5]1
5.103.9057.0010.000.000.000.000.000.000.0090.0060.00[23]6
0.000.0050.000.000.000.000.000.000.000.000.000.00[24]15
0.000.000.000.000.000.000.000.000.000.000.000.00[25]1
0.000.0064.000.00269.000.004.000.000.000.0060.000.00[11]*
0.000.00134.000.000.000.003.500.000.000.0069.000.00[26]*
0.000.000.000.000.000.000.000.000.000.000.000.00[3]1
0.000.000.000.000.000.000.000.000.000.0069.000.00[27]18
FCU0.000.0072.800.00200.000.0010.000.0013.300.00173.000.00[5]1
15.200.0048.000.000.000.000.000.000.000.0047.000.00[23]6
0.000.0042.000.000.000.000.000.000.000.000.000.00[24]15
0.000.0056.000.000.000.000.000.0049.700.00842.000.00[25]1
0.000.0051.000.00265.000.0012.000.000.000.00290.000.00[11]*
0.000.00228.000.000.000.0012.100.000.000.00342.000.00[26]*
26.800.0040.3014.500.000.0015.507.000.000.00631.900.00[3]1
0.000.000.000.000.000.000.000.000.000.000.000.00[27]18
ECU0.000.0045.600.00137.000.0011.000.008.900.00184.800.00[5]1
14.905.7045.005.000.000.000.000.000.000.00480.00140.00[23]6
0.000.0045.000.000.000.000.000.000.000.000.000.00[24]15
0.000.0056.000.000.000.000.000.0038.200.00644.000.00[25]1
0.000.0062.000.00228.000.004.000.000.000.00210.000.00[11]*
0.000.00182.000.000.000.003.500.000.000.00260.000.00[26]*
0.000.000.000.000.000.000.000.000.000.000.000.00[3]1
0.000.000.000.000.000.000.000.000.000.00260.000.00[27]18
ECRB0.000.0035.100.00185.000.0012.000.009.600.00258.800.00[5]1
15.808.1053.006.000.000.000.000.000.000.00450.00140.00[23]6
0.000.0061.000.000.000.000.000.000.000.000.000.00[24]15
0.000.0049.000.000.000.000.000.0040.500.00781.000.00[25]1
0.000.0059.000.00222.000.009.000.000.000.00220.000.00[11]*
0.000.00127.000.000.000.008.900.000.000.00273.000.00[26]*
0.000.000.000.000.000.000.000.000.000.000.000.00[3]1
0.000.000.000.000.000.000.000.000.000.00273.000.00[27]18
ECRL0.000.0072.700.00210.000.000.000.0013.900.00171.100.00[5]1
18.307.5078.005.000.000.000.000.000.000.00370.0090.00[23]6
0.000.0093.000.000.000.000.000.000.000.000.000.00[24]15
0.000.0090.000.000.000.000.000.0047.400.00500.000.00[25]1
0.000.0081.000.00224.000.000.000.000.000.00220.000.00[11]*
0.000.0093.700.000.000.002.500.000.000.00146.000.00[26]*
0.000.000.000.000.000.000.000.000.000.000.000.00[3]1
0.000.000.000.000.000.000.000.000.000.000.000.00[27]18
* Simulation study: M. flexor carpi radialis (FCR), M. palmaris longus (PL), M. flexor carpi ulnaris (FCU), M. extensor carpi ulnaris (ECU), Mm. extensors carpi radialis longus et brevis (ECRL and ECRB); PCSA: physiological cross-sectional area; SD: Standard deviation; n: numbers of the included specimen.
Table 2. Parameters for the SA. (M. flexor carpi radialis (FCR), M. palmaris longus (PL), M. flexor carpi ulnaris (FCU), M. extensor carpi ulnaris (ECU), Mm. extensors carpi radialis longus et brevis (ECRL and ECRB); PCSA: physiological cross-sectional area; PFC: Peak force calculation).
Table 2. Parameters for the SA. (M. flexor carpi radialis (FCR), M. palmaris longus (PL), M. flexor carpi ulnaris (FCU), M. extensor carpi ulnaris (ECU), Mm. extensors carpi radialis longus et brevis (ECRL and ECRB); PCSA: physiological cross-sectional area; PFC: Peak force calculation).
MusclePCSAMoment ArmPFC
[mm2][mm][N/mm2]
FlexExtRDUD
MinMaxMinMaxMinMaxMinMaxMinMaxMinMax
FCR906901020./../.515./../.0.21.0
PL301502030./../.010./../.
FCU408501020./../../../.1525
ECU180650./../.010./../.2030
ECRB220790./../.10201020./../.
ECRL140500./../.5152030./../.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Eschweiler, J.; Praster, M.; Quack, V.; Li, J.; Rath, B.; Hildebrand, F.; Migliorini, F. Comparison of Optimization Strategies for Musculoskeletal Modeling of the Wrist for Therapy Planning in Case of Total Wrist Arthroplasty. Life 2022, 12, 527. https://doi.org/10.3390/life12040527

AMA Style

Eschweiler J, Praster M, Quack V, Li J, Rath B, Hildebrand F, Migliorini F. Comparison of Optimization Strategies for Musculoskeletal Modeling of the Wrist for Therapy Planning in Case of Total Wrist Arthroplasty. Life. 2022; 12(4):527. https://doi.org/10.3390/life12040527

Chicago/Turabian Style

Eschweiler, Jörg, Maximilian Praster, Valentin Quack, Jianzhang Li, Björn Rath, Frank Hildebrand, and Filippo Migliorini. 2022. "Comparison of Optimization Strategies for Musculoskeletal Modeling of the Wrist for Therapy Planning in Case of Total Wrist Arthroplasty" Life 12, no. 4: 527. https://doi.org/10.3390/life12040527

APA Style

Eschweiler, J., Praster, M., Quack, V., Li, J., Rath, B., Hildebrand, F., & Migliorini, F. (2022). Comparison of Optimization Strategies for Musculoskeletal Modeling of the Wrist for Therapy Planning in Case of Total Wrist Arthroplasty. Life, 12(4), 527. https://doi.org/10.3390/life12040527

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

Article Metrics

Back to TopTop