Non-linear finite element model to assess the effect of tendon forces on the foot-ankle complex

Biomechanical role of bone grafting for calcaneal fracture fixation in the presence of bone defect: A finite element analysis

BACKGROUND

The purpose of this study was to compare the biomechanical stress and stability of calcaneal fixations with and without bone defect, before and after bone grafting, through a computational approach.

METHODS

A finite element model of foot-ankle complex was reconstructed, impoverished with a Sanders III calcaneal fracture without bone defect and with moderate and severe bone defects. Plate fixations with and without bone grafting were introduced with walking stance simulated. The stress and fragment displacement of the calcaneus were evaluated.

FINDINGS

Moderate and severe defect increased the calcaneus stress by 16.11% and 32.51%, respectively and subsequently decreased by 10.76% and 20.78% after bone grafting. The total displacement was increased by 3.99% and 24.26%, respectively by moderate and severe defect, while that of posterior joint facet displacement was 86.66% and 104.44%. The former was decreased by 25.73% and 35.96% after grafting, while that of the latter was reduced by 88.09% and 84.78% for moderate and severe defect, respectively.

INTERPRETATION

Our finite element prediction supported that bone grafting for fixation could enhance the stability and reduce the risk of secondary stress fracture in cases of bone defect in calcaneal fracture.

A finite element model of human hindfoot and its application in supramalleolar osteotomy

BACKGROUND

The majority of the ankle osteoarthritis cases are posttraumatic and affect younger patients with a longer projected life span. Hence, joint-preserving surgery, such as supramalleolar osteotomy becomes popular among young patients, especially those with asymmetric arthritis due to alignment deformities. However, there is a lack of biomechanical studies on postoperative evaluation of stress at ankle joints. We aimed to construct a verifiable finite element model of the human hindfoot, and to explore the effect of different osteotomy parameters on the treatment of varus ankle arthritis.

METHODS

The bones of the hindfoot are reconstructed using normal CT tomography data from healthy volunteers, while the cartilages and ligaments are determined from the literature. The finite element calculation results are compared with the weight-bearing CT (WBCT) data to validate the model. By setting different model parameters, such as the osteotomy height (L) and the osteotomy distraction distance (h), the effects of different surgical parameters on the contact stress of the ankle joint surface are compared.

FINDINGS

The alignment and the deformation of hindfoot bones as determined by the finite element analysis aligns closely with the data obtained from WBCT. The maximum contact stress of the ankle joint surface calculated by this model increases with the increase of the varus angle. The maximum contact stresses as a function of the L and h of the ankle joint surface are determined.

INTERPRETATION

The relationship between surgical parameters and stress at the ankle joint in our study could further help guiding the planning of the supramalleolar osteotomy according to the varus/valgus alignment of the patients.

Comparison of different internal fixation models in ankle arthrodesis using 3D finite-element analysis

BACKGROUND

The purpose of this study is to use three-dimensional finite-element analysis to better understand the biomechanical features of various internal fixators for ankle arthrodesis.

METHODS

We used finite-element analysis to compare four different types of internal fixations in ankle arthrodesis: Group A had three crossed screws (Ø6.5 mm); Group B had two crossed screws (Ø6.5 mm) and an anterior plate (Ø2.7 mm); Group C only had an anterior anatomical plate (Ø3.5 mm); Group D had one anterior anatomical plate (Ø3.5 mm) and one posterior-lateral screw (Ø6.5 mm). We adopted Ansys 21.0 software to analyze and compare the four types in terms of the displacement of the arthrodesis surface and the stress peak and stress distribution of these models under intorsion, extorsion, dorsiflexion torque, and neutral vertical load.

RESULTS

① Displacement of the arthrodesis surface: In Group A, the maximum displacement was larger than Group D under neutral vertical load and dorsiflexion torque but less than it under intorsion and extorsion torque. In Group B, the maximum displacement against dorsiflexion, neutral vertical load, intorsion, and extorsion was less than that in the other three fixation models. In Group C, the maximum displacement against the above four loading patterns were significantly higher than that in another three fixation models. ② Stress peak and stress distribution: based on the stress distribution of the four models, the peak von Mises stress was concentrated in the central sections of the compression screws, plate joints, and bending parts of the plates.

CONCLUSIONS

The fixation model consisting of two crossed screws and an anterior outperformed the other three fixation models in terms of biomedical advantages; thus, this model can be deemed a safe and reliable internal fixation approach for ankle arthrodesis.

Single-cell RNA sequencing reveals differences between force application and bearing in ankle cartilage

Chondrocytes are the major functional elements of articular cartilage. Force has been demonstrated to influence the structure and function of articular cartilage and chondrocytes. Therefore, it is necessary to evaluate chondrocytes under different force conditions to gain deep insight into chondrocyte function. Six cartilage tissues from the distal tibia (referred to as the AT group) and five cartilage tissues from the trochlear surface of the talus (referred to as the ATa group) were obtained from 6 donors who had experienced fatal accidents. Single-cell RNA sequencing was used on these samples. A total of 149,816 cells were analyzed. Nine chondrocyte subsets were ultimately identified. Pseudotime analyses, enrichment analyses, cell-cell interaction studies, and single-cell regulatory network inference and clustering were performed for each cell type, and the differences between the AT and ATa groups were analyzed. Immunohistochemical staining was used to verify the existence of each chondrocyte subset and its distribution. The results suggested that reactive oxygen species related processes were active in the force-applied region, while tissue repair processes were common in the force-bearing region. Although the number of prehypertrophic chondrocytes was small, these chondrocytes seemed to play an important role in the ankle.

Clinically useful finite element models of the natural ankle - A review

BACKGROUND

Biomechanical simulation of the foot and ankle complex is a growing research area but compared to simulation of joints such as hip and knee, it has been under investigated and lacks consistency in research methodology. The methodology is variable, data is heterogenous and there are no clear output criteria. Therefore, it is very difficult to correlate clinically and draw meaningful inferences.

METHODS

The focus of this review is finite element simulation of the native ankle joint and we will explore: the different research questions asked, the model designs used, ways the model rigour has been ensured, the different output parameters of interest and the clinical impact and relevance of these studies.

FINDINGS

The 72 published studies explored in this review demonstrate wide variability in approach. Many studies demonstrated a preference for simplicity when representing different tissues, with the majority using linear isotropic material properties to represent the bone, cartilage and ligaments; this allows the models to be complex in another way such as to include more bones or complex loading. Most studies were validated against experimental or in vivo data, but a large proportion (40%) of studies were not validated at all, which is an area of concern.

INTERPRETATION

Finite element simulation of the ankle shows promise as a clinical tool for improving outcomes. Standardisation of model creation and standardisation of reporting would increase trust, and enable independent validation, through which successful clinical application of the research could be realised.

Model of calf muscle tear during a simulated eccentric contraction, comparison between ex-vivo experiments and discrete element model

The tearing of the muscle-tendon complex (MTC) is one of the common sports-related injuries. A better understanding of the mechanisms of rupture and its location could help clinicians improve the way they manage the rehabilitation period of patients. A new numerical approach using the discrete element method (DEM) may be an appropriate approach, as it considers the architecture and the complex behavior of the MTC. The aims of this study were therefore: first, to model and investigate the mechanical elongation response of the MTC until rupture with muscular activation. Secondly, to compare results with experimental data, ex vivo tensile tests until rupture were done on human cadavers {triceps surae muscle + Achilles tendon}. Force/displacement curves and patterns of rupture were analyzed. A numerical model of the MTC was completed in DEM. In both numerical and experimental data, rupture appeared at the myotendinous junction (MTJ). Moreover, force/displacement curves and global rupture strain were in agreement between both studies. The order of magnitude of rupture force was close between numerical (858 N for passive rupture and 996 N-1032 N for rupture with muscular activation) and experimental tests (622 N ± 273 N) as for the displacement of the beginning of rupture (numerical: 28-29 mm, experimental: 31.9 mm ± 3.6 mm). These differences could be explained by choices of DEM model and mechanical properties of MTC's components or their rupture strain values. Here we show that he MTC was broken by fibers' delamination at the distal MTJ and by tendon disinsertion at the proximal MTJ in agreement with experimental data and literature.

A computational model of force within the ligaments and tendons in progressive collapsing foot deformity

Progressive collapsing foot deformity results from degeneration of the ligaments and posterior tibial tendon (PTT). Our understanding of the relationship between their failures remains incomplete. We sought to improve this understanding through computational modeling of the forces in these soft tissues. The impact of PTT and ligament failures on force changes in the remaining ligaments was investigated by quantifying ligament force changes during simulated ligament and tendon cutting in a validated finite element model of the foot. The ability of the PTT to restore foot alignment was also evaluated by increasing the PTT force in a foot with attenuated ligaments and comparing the alignment angles to the intact foot. We found that failure of any one of the ligaments led to overloading the remaining ligaments, except for the plantar naviculocuneiform, first plantar tarsometatarsal, and spring ligaments, where removing one led to unloading the other two. The combined attenuation of the plantar fascia, long plantar, short plantar, and spring ligaments significantly overloaded the deltoid and talocalcaneal ligaments. Isolated PTT rupture had no effect on foot alignment but did increase the force in the deltoid and spring ligaments. Moreover, increasing the force within the PTT to 30% of body weight was effective at restoring foot alignment during quiet stance, primarily through reducing hindfoot valgus and forefoot abduction as opposed to improving arch collapse. Our findings suggest that early intervention might be used to prevent the progression of deformity. Moreover, strengthening the PTT through therapeutic exercise might improve its ability to restore foot alignment.

Biomechanical Analysis of Foot–Ankle Complex during Jogging with Rearfoot Strike versus Forefoot Strike

Background and Aim. In order to reduce foot and ankle injuries induced by jogging, two-foot strike patterns, rearfoot strike (RFS), and forefoot strike (FFS), were adopted and compared. First, RFS jogging and FFS jogging were experimentally studied, so as to acquire kinematic and kinetic data, including foot strike angle, knee flexion angle, and ground reaction force (GRF). Then, a 3D finite element model of foot–ankle complex was reconstructed from the scanned 2D-stacked images. Biomechanical characteristics, including plantar pressure, stress of metatarsals, midfoot bone, calcaneus and cartilage, and tensile force of plantar fascia and ligaments, were obtained. The results showed that RFS jogging and FFS jogging had a similar change trend and a close peak value of GRF. Since possessing more momentum in the push stage and less momentum in the brake stage, FFS jogging could be in favor of a higher jogging speed. However, FFS jogging produced larger metatarsal stress in the 5th metatarsal and much larger tensile force of plantar fascia, which might cause metatarsal fracture and heel pain. While RFS jogging produced larger plantar pressure in the hindfoot area, larger calcaneus stress, and much larger tarsal navicular stress, which might cause heel tissue injury, calcaneus damage, and stress fracture of naviculocuneiform joint. In addition, talocrural and talocalcaneal joint cartilage could bear jogging loads, as the peak contact pressure were both small in RFS jogging and FFS jogging. Therefore, jogging with rearfoot or FFS pattern should be chosen according to the health condition of foot–ankle parts.

Modeling of the Achilles Subtendons and Their Interactions in a Framework of the Absolute Nodal Coordinate Formulation

Experimental results have revealed the sophisticated Achilles tendon (AT) structure, including its material properties and complex geometry. The latter incorporates a twisted design and composite construction consisting of three subtendons. Each of them has a nonstandard cross-section. All these factors make the AT deformation analysis computationally demanding. Generally, 3D finite solid elements are used to develop models for AT because they can discretize almost any shape, providing reliable results. However, they also require dense discretization in all three dimensions, leading to a high computational cost. One way to reduce degrees of freedom is the utilization of finite beam elements, requiring only line discretization over the length of subtendons. However, using the material models known from continuum mechanics is challenging because these elements do not usually have 3D elasticity in their descriptions. Furthermore, the contact is defined at the beam axis instead of using a more general surface-to-surface formulation. This work studies the continuum beam elements based on the absolute nodal coordinate formulation (ANCF) for AT modeling. ANCF beam elements require discretization only in one direction, making the model less computationally expensive. Recent work demonstrates that these elements can describe various cross-sections and materials models, thus allowing the approximation of AT complexity. In this study, the tendon model is reproduced by the ANCF continuum beam elements using the isotropic incompressible model to present material features.

Evaluation of assumptions in foot and ankle biomechanical models

BACKGROUND

A variety of biomechanical models have been used in studies of foot and ankle disorders. Assumptions about the element types, material properties, and loading and boundary conditions are inherent in every model. It was hypothesized that the choice of these modeling assumptions could have a significant impact on the findings of the model.

METHODS

We investigated the assumptions made in a number of biomechanical models of the foot and ankle and evaluated their effects on the results of the studies. Specifically, we focused on: (1) element choice for simulation of ligaments and tendons, (2) material properties of ligaments, cortical and trabecular bones, and encapsulating soft tissue, (3) loading and boundary conditions of the tibia, fibula, tendons, and ground support.

FINDINGS

Our principal findings are: (1) the use of isotropic solid elements to model ligaments and tendons is not appropriate because it allows them to transmit unrealistic bending and twisting moments and compressive forces; (2) ignoring the difference in elastic modulus between cortical and trabecular bones creates non-physiological stress distribution in the bones; (3) over-constraining tibial motion prevents anticipated deformity within the foot when simulating foot deformities, such as progressive collapsing foot deformity; (4) neglecting the Achilles tendon force affects almost all kinetic and kinematic parameters through the foot; (5) the axial force applied to the tibia and fibula is not equal to the ground reaction force due to the presence of tendon forces.

INTERPRETATION

The predicted outcomes of a foot model are highly sensitive to the model assumptions.

New Distinct Component Patterns for Plantar Pressure Variables by Using Principal Component Analysis

BACKGROUND

It is important to determine the plantar pressure distribution of schoolchildren by applying static and dynamic foot analyses using a pedobarography device. However, it is difficult to obtain clear interpretations from results that can be explained by a large number of plantar pressure variables. The aim of this study was to use principal component analysis (PCA) to predict the main components for reducing the size of big data sets, provide a practical overview, and minimize information loss on the subject of plantar pressure assessment in youths.

METHODS

In total, 112 schoolchildren were included in the study (mean ± SD: age, 10.58 ± 1.27 years; body mass index, 18.86 ± 4.33). During the research, a pedobarography device was used to obtain plantar pressure data. Each foot was divided into six anatomical regions and evaluated. Global and regional plantar pressure distributions, load and surface areas, pressure-time integrals, weight ratios, and geometric foot properties were calculated.

RESULTS

The PCA yielded ten principal components that together account for 81.88% of the variation in the data set and represent new and distinct patterns. Thus, 137 variables affecting the subject were reduced to ten components.

CONCLUSIONS

The numerous variables that affect static and dynamic plantar pressure distributions can be reduced to ten components by PCA, making the research results more concise and understandable.

Effects of badminton insole design on stress distribution, displacement and bone rotation of ankle joint during single-leg landing: a finite element analysis

Previous research has reported that up to 92% of injuries amongst badminton players consist of lower limb, whereby 35% of foot fractures occurred at the metatarsal bone. In sports, insoles are widely used to increase athletes' performance and prevent many injuries. However, there is still a lack of badminton insole analysis and improvements. Therefore, this study aimed to biomechanically analyse three different insole designs. A validated and converged three-dimensional (3D) finite element model of ankle-foot complex was developed, which consisted of the skin, talus, calcaneus, navicular, three cuneiform, cuboid, five metatarsals and five phalanges. Three existing insoles from the market, (1) Yonex Active Pro Truactive, (2) Victor VT-XD 8 and (3) Li-Ning L6200LA, were scanned using a 3D scanner. For the analysis, single-leg landing was simulated. On the superior surface of the skin, 2.57 times of the bodyweight was axially applied, and the inferior surface of the outsole was fixed. The results showed that Insole 3 was the most optimum design to reduce peak stress on the metatarsals (3.807 MPa). In conclusion, the optimum design of Insole 3, based on the finite element analysis, could be a justification of athletes' choices to prevent injury and other complications.

Workflow assessing the effect of Achilles tendon rupture on gait function and metatarsal stress: Combined musculoskeletal modeling and finite element analysis

Achilles tendon rupture (ATR) incidence has increased among badminton players in recent years. The foot internal stress was hard to obtain through experimental testing. The purpose of the current research is to develop a methodology that could improve the finite element model derived foot internal stress prediction for ATR clinical and rehabilitation applications. A subject-specific musculoskeletal model was combined with a 3D finite element model to predict the metatarsal stress. The 80% point during the push-off phase of walking was selected for the comparing between injured and uninjured sides. The surgical repaired Achilles tendon (AT) after 12 months was elongated by 5.5% than the uninjured tendon. At 80% point of stance phase, the ankle plantarflexion angle and AT force decreased by 39.6% and 21.9% on the injured side, respectively. The foot inversion degree increased by 22.9% and was accompanied by the redistribution of metatarsals von Mises stress. The stresses on the fourth and fifth metatarsals were increased by 59.5% and 85.9% on the injured side. The workflow is available to assess musculoskeletal disorders and obtain foot internal stress after ATR. The decreased ankle plantar flexor force may be affected by triceps surae muscle atrophy and weakened force transmission ability of elongated AT. The increased von Mises stress on fourth and fifth metatarsals accompanied by higher foot inversion may increase the ankle lateral sprain injury risk.

The impact of ligament tears on joint contact mechanics in progressive collapsing foot deformity: A finite element study

BACKGROUND

Patients with longstanding progressive collapsing foot deformity often develop osteoarthritis of the ankle, midfoot, or hindfoot joints, which can be symptomatic or lead to fixed deformities that complicate treatment. The development of deformity is likely caused by ligament degeneration and tears. However, the effect of individual ligament tears on changes in joint contact mechanics has not been investigated.

METHODS

A validated finite element model of the foot was used to compare joint contact areas, forces, and pressures between the intact and collapsed foot, and to evaluate the effect of individual ligament tears on joint contact mechanics.

FINDINGS

Collapsing the foot resulted in an increase in contact pressure in the subtalar, calcaneocuboid, tibiotalar, medial naviculocuneiform, and first tarsometatarsal joints but a decrease in contact pressure in the talonavicular joint. Rupture of the spring ligament was the main contributor to increased calcaneocuboid and subtalar joint contact pressures and decreased medial naviculocuneiform and first tarsometatarsal joint contact pressures, as well as talonavicular subluxation. Deltoid ligament rupture was the primary source of increased contact pressure in the medial naviculocuneiform, first tarsometatarsal, and tibiotalar joints.

INTERPRETATION

Degenerative tearing of the ligaments in flatfoot deformity leads to increased joint contact pressures, primarily in the calcaneocuboid, subtalar, and tibiotalar joints, which has been implicated in the development of osteoarthritis in these joints. An improved understanding of the relationship between ligament tears and joint contact pressures could provide support for the use of ligament reconstructions to prevent the development of arthrosis.

Finite Element Modelling of Planus and Rectus Foot Types for the Study of First Metatarsophalangeal and First Metatarsocuneiform Joint Contact Mechanics

Evaluating the contact mechanics of human joints is an important element in understanding the pathomechanics of orthopaedic diseases. Although physical testing is essential in the evaluation process, reliable computational models can augment these experiments by non-invasive predictions of biomechanical or surgical variables. The objective of this study was to perform verification of a framework for developing a medial forefoot finite element. Verification was conducted by comparing computational predictions to experimental measurements of first metatarsophalangeal and first metatarsocuneiform joint contact mechanics. A custom-built force-controlled cadaveric test-rig was used to derive measurements of contact pressure, force, and area. A quasi-static finite element was developed and driven under the same boundary and loading conditions. Calibration of cartilage moduli and mesh sensitivity analyses were performed. Mean errors in contact pressures, forces, and areas were 24%, 4%, and 40% at the first metatarsophalangeal joint and 23%, 12%, and 19% at the first Metatarsocuneiform joint, respectively. Verification of a medial forefoot finite element model development framework was presented and found to be within 30% for contact pressure and contact force of both joints. This study presents a method to verify and simulate realistic physiological loading to investigate orthopaedic diseases of the medial forefoot.

A Review of Finite Element Models of Ligaments in the Foot and Considerations for Practical Application

PURPOSE

Finite element (FE) modeling has been used as a research tool for investigating underlying ligaments biomechanics and orthopedic applications. However, FE models of the ligament in the foot have been developed with various configurations, mainly due to their complex 3D geometry, material properties, and boundary conditions. Therefore, the purpose of this review was to summarize the current state of finite element modeling approaches that have been used in the ?eld of ligament biomechanics, to discuss their applicability to foot ligament modeling in a practical setting, and also to acknowledge current limitations and challenges.

METHODS

A comprehensive literature search was performed. Each article was analyzed in terms of the methods used for: (a) ligament geometry, (b) material property, (c) boundary and loading condition related to its application, and (d) model verification and validation.

RESULTS

Of the reviewed studies, 80% of the studies used simplified representations of ligament geometry, the non-linear mechanical behavior of ligaments was taken into account in only 19.2% of the studies, 33% of included studies did not include any kind of validation of the FE model.

CONCLUSION

Further refinement in the functional modeling of ligaments, the micro-structure level characteristics, nonlinearity, and time-dependent response, may be warranted to ensure the predictive ability of the models.

In Silico Finite Element Analysis of the Foot Ankle Complex Biomechanics: A Literature Review

Computational approaches, especially finite element analysis (FEA), have been rapidly growing in both academia and industry during the last few decades. FEA serves as a powerful and efficient approach for simulating real-life experiments, including industrial product development, machine design, and biomedical research, particularly in biomechanics and biomaterials. Accordingly, FEA has been a "go-to" high biofidelic software tool to simulate and quantify the biomechanics of the foot-ankle complex, as well as to predict the risk of foot and ankle injuries, which are one of the most common musculoskeletal injuries among physically active individuals. This paper provides a review of the in silico FEA of the foot-ankle complex. First, a brief history of computational modeling methods and finite element (FE) simulations for foot-ankle models is introduced. Second, a general approach to build an FE foot and ankle model is presented, including a detailed procedure to accurately construct, calibrate, verify, and validate an FE model in its appropriate simulation environment. Third, current applications, as well as future improvements of the foot and ankle FE models, especially in the biomedical field, are discussed. Finally, a conclusion is made on the efficiency and development of FEA as a computational approach in investigating the biomechanics of the foot-ankle complex. Overall, this review integrates insightful information for biomedical engineers, medical professionals, and researchers to conduct more accurate research on the foot-ankle FE models in the future.

Development and investigation of the applicability of microstructural models for polymeric low density foams directly obtained from computed tomography data

Nowadays, there are several methods to obtain simulation models for foams which consider important microstructural features. This research study presents a method to obtain low density foam models directly from computed tomography (CT) data. Finite element meshes are created from CT measurement results of a polypropylene extrusion foam with two different densities. Sensitivity studies with regard to the tension behaviour are performed with the microstructural models. The study shows that the tension behaviour highly depends on the examined area of the foam because the microstructure and density vary through the foam. Furthermore, the simulation results are validated with experimental results. The validation shows that the tension behaviour of the investigated foams characterised by the simulation approach is in good agreement with the experimentally observed behaviour and that specific microstructural characteristics (e.g. anisotropic cell shapes) are captured in the model.

Modulating nonlinear elastic behavior of biodegradable shape memory elastomer and small intestinal submucosa(SIS) composites for soft tissue repair

Structural repair of soft tissue for regenerative therapies can be advanced by developing biocompatible and bioresorbable materials with mechanical properties similar to the tissue targeted for therapy. Developing new materials modeling soft tissue mechanics can mitigate many limitations of material based therapies, specifically concerning the mechanical stress and deformation the material imposes on surrounding tissue structures. However, many elastomeric materials used in soft tissue repair lack the ability to be delivered through minimally invasive surgical (MIS) or transcatheter routes and require open surgical approaches for placement and application. We have developed a biocompatible and fully biodegradable shape memory elastomer, poly-(glycerol dodecanedioate) (PGD), which fulfills the requirements for hyperelasticity and exhibits shape memory behavior to serve as a novel substrate material for regenerative therapy in minimally invasive clinical procedures. Our previous work demonstrated control over the tangent modulus at 12.5% compressive strain between 1 and 3 MPa by increasing the crosslinking density in the polymer. In order to improve control over a broader range of mechanical properties, nonlinear behavior, and toughness, we 1) varied PGD physical crosslink density, 2) incorporated sheets of porcine small intestinal submucosa (SIS, Cook Biotech, Inc.) with varying thickness, and 3) mixed lyophilized SIS particulates into PGD at different weight percentages. Tensile testing (ASTM D412a) revealed PGD containing SIS sheets of were stiffer than controls (p < 0.01). Incorporating lyophilized SIS particulates into PGD increased the strain to failure (p < 0.001) compared to PGD controls. Test specimens with 1 ply sheets had greater tear strength (ASTM D624c) compared to PGD tear specimens prepared control specimens (p < 0.001). However, incorporating SIS particulates decreased tear strength of PGD-SIS 0.5 wt% particulate composites (p < 0.01) compared to PGD controls. Incorporating 2 ply and 4 ply sheets and 0.5 wt% particulates into PGD decreased the fixity and recovery of composite materials compared to controls (p < 0.01). Nonlinear modeling of stress strain curves under uniaxial tension demonstrated tunability of PGD-SIS composite materials to model various nonlinear soft tissues. These findings support the use of shape memory PGD-SIS composite materials towards the design of implantable devices for a variety of soft tissue regeneration applications by minimally invasive surgery.

Finite element analysis of secondary effect of midfoot fusions on the spring ligament in the management of adult acquired flatfoot

BACKGROUND

Surgical treatment of adult acquired flatfoot deformity can involve arthrodesis of the midfoot to stabilize the medial column. Few experimental studies have assessed the biomechanical effects of these fusions, because of the difficulty of measuring these parameters in cadavers. Our objective was to quantify the biomechanical stress caused by various types of midfoot arthrodesis on the Spring ligament. To date this is not known.

METHODS

An innovative finite element model was used to evaluate flatfoot scenarios treated with various combinations of midfoot arthrodesis. All the bones, cartilages and tissues related to adult acquired flatfoot deformity were included, respecting their biomechanical characteristics. The stress changes on the Spring ligament were quantified. Both foot arch lengthening and falling were measured for each of the midfoot arthrodeses evaluated.

FINDINGS

Arthrodesis performed for stabilization of the talonavicular joint leads to a higher decrease in stress on the Spring ligament. Talonavicular fusion generated a Spring ligament stress decrease of about 61% with respect to the reference case (without any fusion). However, fusing the naviculocuneiform joints leads to an increase in the stress on the Spring ligament.

INTERPRETATION

This important finding has been unknown to date. We advocate caution regarding fusion of the naviculocuneiform joint as it leads to increased stresses across the Spring ligament and therefore accelerates the development of planovalgus.

Cartilage Stiffness Effect on Foot Biomechanics of Chinese Bound Foot: A Finite Element Analysis

The purpose of this study is to investigate the effect of cartilage stiffness on inner foot biomechanics of Chinese bound foot while balanced standing using finite element method. A three-dimensional FE model of bound foot involving 28 bones, 72 ligaments, 5 plantar fascia, cartilages, and encapsulated soft tissue was constructed and validated. To conduct the sensitivity analysis of cartilage stiffness, the incremental Young’s modulus of 1, 5, 10, and 15 MPa were assigned to the cartilage. 25% of the body weight was applied to the Achilles tendon to adjust the anterior- posterior displacement of center of pressure agreeable with the measured result. As the Young’s modulus of cartilage increased, the peak von Mises stress in the fifth metatarsal increased obviously, while that in the calcaneus remains unchanged. The plantar fascia experienced reduced total tension with stiffer cartilage. The cartilage stiffening also caused a general increase of contact pressure at mid- and forefoot joints. Cartilage stiffening due to foot binding gave rise to risks of foot pain and longitude arch damage. Knowledge of this study contributes to the understanding of bound foot biomechanical behavior and demonstrating the mechanism of long-term injury and function damage in terms of weight-bearing due to foot binding.

Foot internal stress distribution during impact in barefoot running as function of the strike pattern

The aim of the present study is to examine the impact absorption mechanism of the foot for different strike patterns (rearfoot, midfoot and forefoot) using a continuum mechanics approach. A three-dimensional finite element model of the foot was employed to estimate the stress distribution in the foot at the moment of impact during barefoot running. The effects of stress attenuating factors such as the landing angle and the surface stiffness were also analyzed. We characterized rear and forefoot plantar sole behavior in an experimental test, which allowed for refined modeling of plantar pressures for the different strike patterns. Modeling results on the internal stress distributions allow predictions of the susceptibility to injury for particular anatomical structures in the foot.