Computational modeling to predict mechanical function of joints: application to the lower leg with simulation of two cadaver studies.

Computational models of musculoskeletal joints and limbs can provide useful information about joint mechanics. Validated models can be used as predictive devices for understanding joint function and serve as clinical tools for predicting the outcome of surgical procedures. A new computational modeling approach was developed for simulating joint kinematics that are dictated by bone/joint anatomy, ligamentous constraints, and applied loading. Three-dimensional computational models of the lower leg were created to illustrate the application of this new approach. Model development began with generating three-dimensional surfaces of each bone from CT images and then importing into the three-dimensional solid modeling software SOLIDWORKS and motion simulation package COSMOSMOTION. Through SOLIDWORKS and COSMOSMOTION, each bone surface file was filled to create a solid object and positioned necessary components added, and simulations executed. Three-dimensional contacts were added to inhibit intersection of the bones during motion. Ligaments were represented as linear springs. Model predictions were then validated by comparison to two different cadaver studies, syndesmotic injury and repair and ankle inversion following ligament transection. The syndesmotic injury model was able to predict tibial rotation, fibular rotation, and anterior/posterior displacement. In the inversion simulation, calcaneofibular ligament extension and angles of inversion compared well. Some experimental data proved harder to simulate accurately, due to certain software limitations and lack of complete experimental data. Other parameters that could not be easily obtained experimentally can be predicted and analyzed by the computational simulations. In the syndesmotic injury study, the force generated in the tibionavicular and calcaneofibular ligaments reduced with the insertion of the staple, indicating how this repair technique changes joint function. After transection of the calcaneofibular ligament in the inversion stability study, a major increase in force was seen in several of the ligaments on the lateral aspect of the foot and ankle, indicating the recruitment of other structures to permit function after injury. Overall, the computational models were able to predict joint kinematics of the lower leg with particular focus on the ankle complex. This same approach can be taken to create models of other limb segments such as the elbow and wrist. Additional parameters can be calculated in the models that are not easily obtained experimentally such as ligament forces, force transmission across joints, and three-dimensional movement of all bones. Muscle activation can be incorporated in the model through the action of applied forces within the software for future studies.

[1]  Stephen M Belkoff,et al.  A Biomechanical Evaluation of the Tibiofibular and Tibiotalar Ligaments of the Ankle , 2003, Foot & ankle international.

[2]  Kai-Nan An,et al.  Consequences of Partial and Total Plantar Fascia Release: A Finite Element Study , 2006, Foot & ankle international.

[3]  F.E. Zajac,et al.  An interactive graphics-based model of the lower extremity to study orthopaedic surgical procedures , 1990, IEEE Transactions on Biomedical Engineering.

[4]  E Y Chao,et al.  Simulation and animation of musculoskeletal joint system. , 1993, Journal of biomechanical engineering.

[5]  H. Grootenboer,et al.  Articular contact in a three-dimensional model of the knee. , 1991, Journal of Biomechanics.

[6]  S. Delp,et al.  Preserving plantar flexion strength after surgical treatment for contracture of the triceps surae: A computer simulation study , 1995, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[7]  Edmund Y S Chao Graphic-based musculoskeletal model for biomechanical analyses and animation. , 2003, Medical engineering & physics.

[8]  S. Siegler,et al.  The Mechanical Characteristics of the Collateral Ligaments of the Human Ankle Joint* , 1988, Foot & ankle.

[9]  J. Mcelhaney,et al.  Biomechanical Characteristics of Human Ankle Ligaments , 1985, Foot & ankle.

[10]  S. Nielsen,et al.  Lateral Talocalcaneal Instability Following Section of the Calcaneofibular Ligament: A Kinesiologic Study , 1987, Foot & ankle.

[11]  J Wismans,et al.  A three-dimensional mathematical model of the knee-joint. , 1980, Journal of biomechanics.

[12]  S J Piazza,et al.  Three-dimensional dynamic simulation of total knee replacement motion during a step-up task. , 2001, Journal of biomechanical engineering.

[13]  E. Chao,et al.  Three-dimensional dynamic hip contact area and pressure distribution during activities of daily living. , 2006, Journal of biomechanics.

[14]  J. Schuberth Orthopaedic surgery essentials: Foot and ankle , 2004 .

[15]  J. Owen,et al.  Comparison of the Syndesmotic Staple to the Transsyndesmotic Screw: A Biomechanical Study , 2005, Foot & ankle international.

[16]  Scott L. Delp,et al.  A Model of the Upper Extremity for Simulating Musculoskeletal Surgery and Analyzing Neuromuscular Control , 2005, Annals of Biomedical Engineering.

[17]  Ming Zhang,et al.  Three-dimensional finite element analysis of the foot during standing--a material sensitivity study. , 2005, Journal of biomechanics.

[18]  S. Delp,et al.  How muscle architecture and moment arms affect wrist flexion-extension moments. , 1997, Journal of biomechanics.

[19]  James H. Herndon,et al.  Tensile properties of the interosseous membrane of the human forearm , 1996 .

[20]  S L Woo,et al.  Measurement of changes in ligament tension with knee motion and skeletal maturation. , 1990, Journal of biomechanical engineering.

[21]  R. Adelaar,et al.  Elongation Behavior of Calcaneofibular and Cervical Ligaments during Inversion Loads Applied in an Open Kinetic Chain , 1998, Foot & ankle international.

[22]  S. Hirokawa,et al.  Three-dimensional mathematical model analysis of the patellofemoral joint. , 1991, Journal of biomechanics.

[23]  L Blankevoort,et al.  Ligament-bone interaction in a three-dimensional model of the knee. , 1991, Journal of biomechanical engineering.

[24]  Fong-Chin Su,et al.  Use of virtual, interactive, musculoskeletal system (VIMS) in modeling and analysis of shoulder throwing activity. , 2005, Journal of biomechanical engineering.

[25]  B Reggiani,et al.  Finite element analysis of a total ankle replacement during the stance phase of gait. , 2006, Journal of biomechanics.

[26]  S. D. Kwak,et al.  A Mathematical Formulation for 3D Quasi-Static Multibody Models of Diarthrodial Joints , 2000, Computer methods in biomechanics and biomedical engineering.

[27]  Loyd M. Myers,et al.  A Modeling and Simulation System for the Human Hand , 1990 .

[28]  F. Netter Atlas of Human Anatomy , 1967 .

[29]  S L Delp,et al.  A graphics-based software system to develop and analyze models of musculoskeletal structures. , 1995, Computers in biology and medicine.

[30]  Lawrence M. Elson,et al.  The Anatomy Coloring Book , 1977 .

[31]  M. S. Downey,et al.  McGlamry's Comprehensive Textbook of Foot and Ankle Surgery , 2001 .

[32]  S. Delp,et al.  Evaluation of a Deformable Musculoskeletal Model for Estimating Muscle–Tendon Lengths During Crouch Gait , 2001, Annals of Biomedical Engineering.