Multi-body simulation of a canine hind limb: model development, experimental validation and calculation of ground reaction forces

Numerical evaluation of internal femur osteosynthesis based on a biomechanical model of the loading in the proximal equine hindlimb

Femoral fractures are often considered lethal for adult horses because femur osteosynthesis is still a surgical challenge. For equine femur osteosynthesis, primary stability is essential, but the detailed physiological forces occurring in the hindlimb are largely unknown. The objective of this study was to create a numerical testing environment to evaluate equine femur osteosynthesis based on physiological conditions. The study was designed as a finite element analysis (FEA) of the femur using a musculoskeletal model of the loading situation in stance. Relevant forces were determined in the musculoskeletal model via optimization. The treatment of four different fracture types with an intramedullary nail was investigated in FEA with loading conditions derived from the model. The analyzed diaphyseal fracture types were a transverse (TR) fracture, two oblique fractures in different orientations (OB-ML: medial-lateral and OB-AP: anterior-posterior) and a "gap" fracture (GAP) without contact between the fragments. For the native femur, the most relevant areas of increased stress were located distally to the femoral head and proximally to the caudal side of the condyles. For all fracture types, the highest stresses in the implant material were present in the fracture-adjacent screws. Maximum compressive (-348 MPa) and tensile stress (197 MPa) were found for the GAP fracture, but material strength was not exceeded. The mathematical model was able to predict a load distribution in the femur of the standing horse and was used to assess the performance of internal fixation devices via FEA. The analyzed intramedullary nail and screws showed sufficient stability for all fracture types.

A three-dimensional musculoskeletal model of the dog

The domestic dog is interesting to investigate because of the wide range of body size, body mass, and physique in the many breeds. In the last several years, the number of clinical and biomechanical studies on dog locomotion has increased. However, the relationship between body structure and joint load during locomotion, as well as between joint load and degenerative diseases of the locomotor system (e.g. dysplasia), are not sufficiently understood. Collecting this data through in vivo measurements/records of joint forces and loads on deep/small muscles is complex, invasive, and sometimes unethical. The use of detailed musculoskeletal models may help fill the knowledge gap. We describe here the methods we used to create a detailed musculoskeletal model with 84 degrees of freedom and 134 muscles. Our model has three key-features: three-dimensionality, scalability, and modularity. We tested the validity of the model by identifying forelimb muscle synergies of a walking Beagle. We used inverse dynamics and static optimization to estimate muscle activations based on experimental data. We identified three muscle synergy groups by using hierarchical clustering. The activation patterns predicted from the model exhibit good agreement with experimental data for most of the forelimb muscles. We expect that our model will speed up the analysis of how body size, physique, agility, and disease influence neuronal control and joint loading in dog locomotion.

Development of a three-dimensional computer model of the canine pelvic limb including cruciate ligaments to simulate movement

Gait analysis as subjective visual assessment forms the foundation of the veterinarian's lameness examination. Pelvic limb lameness is frequently seen in dogs and the stifle joint with its cruciate ligaments, is a main cause of lameness due to cruciate ligament deficiency. In this study, we developed an open-source three-dimensional musculoskeletal pelvic limb model of a 30 kg Labrador Retriever including cruciate ligaments, simulating the gait cycle of the walking movement with the open-source programs NMSBuilder (Institutio Ortopedico Rizzoli, Bologna, Italy) and OpenSim (National Center for Simulation in Rehabilitation Research (NCSRR), Stanford, CA, USA). The computer model generated muscle activations based on motion data. The computed activations were similar to experimental electromyogram data. Highest joint torque was in extension/flexion in the stifle joint at 54 Nm at 14% of the gait cycle with cruciate ligaments. Highest stifle joint reaction force was 408 N at 16% of the gait cycle and was reduced after adding cruciate ligaments. Especially the cranial cruciate ligament loads up to 102 N (34% body weight). Cranial cruciate ligament forces increase with stifle extension and decrease with stifle flexion. On the contrary, the caudal cruciate ligament loads up to 27 N (9% body weight) during the swing phase with a flexed stifle joint. The model was validated with electromyogram data. The model's predictions are plausible because joint torques and forces match the applied ground reaction forces in curve progression and in timing. This model forms a basis for further investigations into stifle surgery after cruciate ligament deficiency.

Impact of wearing a service vest on three-dimensional truncal motion in dogs

OBJECTIVE

To develop a 3-D kinematic model to measure truncal motion in dogs and assess changes in truncal motion in dogs when wearing each of 2 service vests.

ANIMALS

5 adult mixed-breed dogs.

PROCEDURES

27 reflective markers were placed on the pelvis, trunk, and scapula of each dog. Six infrared cameras were placed around a treadmill to track the location of the markers within a calibrated space. Dogs were recorded during walking and trotting on the treadmill. Local and global coordinate systems were established, and a segmental rigid-body model of the trunk was created. Dogs were then recorded while wearing a custom vest and an adjustable vest during walking and trotting on the treadmill. Range of motion of the trunk when dogs were and were not wearing vests was compared by repeated-measures ANOVA.

RESULTS

An anatomic coordinate system was established by use of markers located at T1, T13, and the xiphoid process. Range of motion of the trunk during a gait cycle did not differ significantly regardless of the day of the test for both walking and trotting gaits. Trunk motion of dogs when walking and trotting was significantly reduced when dogs were wearing a vest, compared with trunk motion when not wearing a vest.

CONCLUSIONS AND CLINICAL RELEVANCE

A 3-D kinematic model for measuring truncal rotation was developed. Results indicated measurable differences in the gait of dogs when wearing each of the 2 service vests, compared with the gait when not wearing a vest.

Biomechanical comparison of a new expandable intramedullary nail and conventional intramedullary nails for femoral osteosynthesis in dogs

Intramedullary nailing of diaphyseal femoral fractures is a commonly used treatment method in dogs because of its biological and biomechanical advantages compared to bone plating. To achieve adequate resistance of the intramedullary nail against torsional and axial compressive forces, additional application of transcortical screws is needed. As these interlocking screws represent a frequent cause of post-operative complications, a new expandable intramedullary nail (EXPN) was developed, which was designed to provide adequate fracture stabilisation without the need for transcortical fixation. The evaluation of the biomechanical properties of the new EXPN with regard to torsional, compressive and bending stability as well as direct comparison to the biomechanical properties of conventional Steinmann (STMN)- and interlocking (ILN) nails was carried out with different biomechanical test arrangements. No significant statistical differences regarding the torsional and bending resistance between the EXPN and ILN group were seen, which indicates that rotatory as well as bending stability of the innovative EXPN is similar to the conventional ILN. Nevertheless, the percentage deviation between the attempted and successfully reached physiological compressive forces was significantly higher (p = 0.045) in the EXPN group compared to the ILN group, which indicates that the compressive stability of the innovative EXPN might be weaker compared to the ILN. In summary, the new EXPN represents an interesting alternative to conventional intramedullary nails. However, in direct comparison to conventional interlocking nails, the EXPN has shown weaknesses in the neutralization of axial compressive forces, which indicates that at least biomechanically the interlocking nail seems advantageous. Further in-vitro and in-vivo investigations are required before clinical use can be recommended.

Development of a Canine Rigid Body Musculoskeletal Computer Model to Evaluate Gait

Background

Kinematic and kinetic analysis have been used to gain an understanding of canine movement and joint loading during gait. By non-invasively predicting muscle activation patterns and forces during gait, musculoskeletal models can further our understanding of normal variability and muscle activation patterns and force profiles characteristic of gait.

Methods

Pelvic limb kinematics and kinetics were measured for a 2 year old healthy female Dachshund (5.4 kg) during gait using 3-D motion capture and force platforms. A computed tomography scan was conducted to acquire pelvis and pelvic limb morphology. Using the OpenSim modeling platform, a bilateral pelvic limb subject-specific rigid body musculoskeletal computer model was developed. This model predicted muscle activation patterns, muscle forces, and angular kinematics and joint moments during walking.

Results

Gait kinematics determined from motion capture matched those predicted by the model, verifying model accuracy. Primary muscles involved in generating joint moments during stance and swing were predicted by the model: at mid-stance the adductor magnus et brevis (peak activation 53.2%, peak force 64.7 N) extended the hip, and stifle flexor muscles (biceps femoris tibial and calcaneal portions) flexed the stifle. Countering vertical ground reaction forces, the iliopsoas (peak activation 37.9%, peak force 68.7 N) stabilized the hip in mid-stance, while the biceps femoris patellar portion stabilized the stifle in mid-stance and the plantar flexors (gastrocnemius and flexor digitorum muscles) stabilized the tarsal joint during early stance. Transitioning to swing, the iliopsoas, rectus femoris and tensor fascia lata flexed the hip, while in late swing the adductor magnus et brevis impeded further flexion as biceps femoris tibial and calcaneal portions stabilized the stifle for ground contact.

Conclusion

The musculoskeletal computer model accurately replicated experimental canine angular kinematics associated with gait and was used to predict muscle activation patterns and forces. Thus, musculoskeletal modeling allows for quantification of measures such as muscle forces that are difficult or impossible to measure in vivo.

Limb Kinematics, Kinetics and Muscle Dynamics During the Sit-to-Stand Transition in Greyhounds

Standing up from a prone position is a critical daily activity for animals: failing to do so effectively may cause an injurious fall or increase predation susceptibility. This sit-to-stand behaviour (StS) is biomechanically interesting because it necessitates transitioning through near-maximal joint motion ranges from a crouched (i.e., poor mechanical advantage) to a more upright posture. Such large joint excursions should require large length changes of muscle-tendon units. Here we integrate experimental and musculoskeletal simulation methods to quantify the joint motions, limb forces, and muscle fibre forces, activations and length changes during StS in an extreme athlete-the greyhound-which has large hindlimb muscles bearing short-fibred distal muscles and long tendons. Study results indicate that hindlimb anti-gravity muscle fibres operate near their ~50% limits of length change during StS; mostly by starting at highly lengthened positions. StS also requires high muscle activations (>50%), in part due to non-sagittal motions. Finally, StS movements require passive non-muscular support in the distal hindlimb where short-fibred muscles are incapable of sustaining StS themselves. Non-locomotor behaviours like StS likely impose important trade-offs between muscle fibre force capacity and length changes, as well as active and passive mechanisms of support, that have been neglected in locomotor biomechanics studies.