Effects of internal stress concentrations in plantar soft-tissue—A preliminary three-dimensional finite element analysis

An exploratory in-situ dynamic mechanical analysis on the shearing stress-strain mechanism of human plantar soft tissue

A DMA (dynamic mechanical analysis)-like device based on the principle of classical viscoelasticity testing is invented to investigate the in-situ/in-vivo shear-bearing mechanism of plantar soft tissue. Forty-three volunteers were recruited for the shear-strain test in the longitudinal and transverse directions at five anatomical spots on the plantar surface. Several encouraging observations indicated significant variances among different spots and individuals, implying that the outer forefoot surrounding the second, fifth metatarsal head is a more intensive shear-bearing region on the plantar surface compared to the inner forefoot under the first metatarsal head, and drawing the hypothesis of a significant effect of BMI on the shear-bearing property. The speculations agree with our expectations and other previous research. The feasibility and practical value of this novel approach are substantiated, and these intriguing discoveries provide foundational underpinnings for further in-depth investigations.

In-shoe plantar pressure measurement technologies for the diabetic foot: A systematic review

Introduction

Loss of cutaneous protective sensation and high plantar pressures increase the risk for diabetic foot patients. Trauma and ulceration are imminent threats, making assessment and monitoring essential. This systematic review aims to identify systems and technologies for measuring in-shoe plantar pressures, focusing on the at-risk diabetic foot population.

Methods

A systematic search was conducted across four electronic databases (Scopus, Web of Science, PubMed, Oxford Journals) using PRISMA methodology, covering articles published in English from 1979 to 2024. Only studies addressing systems or sensors exclusively measuring plantar pressures inside the shoe were included.

Results

A total of 87 studies using commercially available devices and 45 articles proposing new systems or sensors were reviewed. The prevailing market offerings consist mainly of instrumented insoles. Emerging technologies under development often feature configurations with four, six or eight resistive sensors strategically placed within removable insoles. Despite some variability due to the inherent heterogeneity of human gait, these devices assess plantar pressure, although they present significant differences between them in measurement results. Individuals with diabetic foot conditions appears exhibit elevated plantar pressures, with reported peak pressures reaching approximately 1000 kPa. The results also showed significant differences between the diabetic and non-diabetic groups.

Conclusion

Instrumented insoles, particularly those incorporating resistive sensor technology, dominate the field. Systems employing eight sensors at critical locations represent a pragmatic approach, although market options extend to systems with up to 960 sensors. Differences between devices can be a critical factor in measurement and highlights the importance of individualized patient assessment using consistent measurement devices.

A combined musculoskeletal and finite element model of a foot to predict plantar pressure distribution

In this study, a combined subject-specific numerical and experimental investigation was conducted to explore the plantar pressure of an individual. The research utilized finite element (FE) and musculoskeletal modelling based on computed tomography (CT) images of an ankle-foot complex and three-dimensional gait measurements. Muscle forces were estimated using an individualized multi-body musculoskeletal model in five gait phases. The results of the FE model and gait measurements for the same subject revealed the highest stress concentration of 0.48 MPa in the forefoot, which aligns with previously-reported clinical observations. Additionally, the study found that the encapsulated soft tissue FE model with hyper-elastic properties exhibited higher stresses compared to the model with linear-elastic properties, with maximum ratios of 1.16 and 1.88 MPa in the contact pressure and von-Mises stress, respectively. Furthermore, the numerical simulation demonstrated that the use of an individualized insole caused a reduction of 8.3% in the maximum contact plantar pressure and 14.7% in the maximum von-Mises stress in the encapsulated soft tissue. Overall, the developed model in this investigation holds potential for facilitating further studies on foot pathologies and the improvement of rehabilitation techniques in clinical settings.

Pressure change at ankle joint in supramalleolar osteotomy with or without fibular osteotomy according to different types of varus ankle

The need for additional fibular osteotomy (FO) when performing supramalleolar osteotomy (SMO) in the varus ankle arthritis (VAA) is controversial. Some cadaveric studies have been performed to prove this; however, it is difficult to implement deformities including talar tilting and translation in cadavers. In this study, we created a model of VAA with the tilting and translation using three-dimensional (3-D) finite element (FE) analysis and analyzed the results of SMO with or without FO depending on the types of VAA. The validated normal foot and ankle 3-D FE model was constructed including the ankle cartilages of the talar dome and tibia plafond. The VAA models were determined and reconstructed by following the classification of VAA, VAA with medial translation for stage 3a, VAA with varus tilting (7.5°) for stage 3b. The postoperative SMO models (SMO with and without FO) were reconstructed by corresponding to each VAA models. The FE analysis conditions were commonly applied. The boundary condition of ankle joint was defined as "sliding condition" and applied 0.002 friction coefficient to realize lubricative property. Loading condition was assumed as a two-leg standing position and half of the subject body weight (325 N) was loaded on center of ground to vertical direction. Contact pressure changes were predicted at the medial ankle cartilage. As a result, in VAA with medial translation, isolated SMO may provide sufficient pressure reduction at the medial ankle joint. However, in VAA with varus tilting, SMO combined with FO could appropriately relieve concentrated pressure at the medial ankle joint.

Plantar pressure thresholds as a strategy to prevent diabetic foot ulcers: A systematic review

Background

The development of ulcers in the plantar region of the diabetic foot originates mainly from sites subjected to high pressure. The monitoring of these events using maximum allowable pressure thresholds is a fundamental procedure in the prevention of ulceration and its recurrence.

Objective

The aim of this review was to identify data in the literature that reveal an objective threshold of plantar pressure in the diabetic foot, where pressure is classified as promoting ulceration. The aim is not to determine the best and only pressure threshold for ulceration, but rather to clarify the threshold values most used in clinical practice and research, also considering the devices used and possible applications for offloading plantar pressure.

Design

A systematic review.

Methods

The search was performed in three electronic databases, by the PRISMA methodology, for studies that used a pressure threshold to minimize the risk of ulceration in the diabetic foot. The selected studies were subjected to eligibility criteria.

Results

Twenty-six studies were included in this review. Seven thresholds were identified, five of which are intended for the inside of the shoe: a threshold of average peak pressure of 200 kPa; 25 % and 40-80 % reduction from initial baseline pressure; 32-35 mm Hg for a capillary perfusion pressure; and a matrix of thresholds based on patient risk, shoe size and foot region. Two other thresholds are intended for the barefoot, 450 and 750 kPa. The threshold of 200 kPa of pressure inside the shoe is the most agreed upon among the studies. Regarding the prevention of ulceration and its recurrence, the efficacy of the proposed threshold matrix and the threshold of reducing baseline pressure by 40-80 % has not yet been evaluated, and the evidence for the remaining thresholds still needs further studies.

Conclusions

Some heterogeneity was found in the studies, especially regarding the measurement systems used, the number of regions of interest and the number of steps to be considered for the threshold. Even so, this review reveals the way forward to obtain a threshold indicative of an effective steppingstone in the prevention of diabetic foot ulcer.

Mechanical characteristics of diabetic and non-diabetic plantar skin

Diabetic foot ulceration is linked to high amputation and mortality rates, with the substantial associated annual spend on the at-risk diabetic foot reflecting the intensive time and labour involved in treatment. Assessing plantar interactions and developing improved understanding of the formation pathways of diabetic ulceration is important to orthotic interventions and patient outcomes. Plantar skin surrogates which emulate the mechanical and tribological characteristics can help improve physical models of ulceration, reduce reliance on cadaveric use and inform more complex computational modelling approaches. The information available from existing studies to characterise plantar skin is limited, typically featuring ex-vivo representations of skin and subcutaneous tissue combined and given focus to shear studies with time dependency. The aim of this study is to improve understanding of plantar tissue mechanics by assessing the mechanical characteristics of plantar skin in two groups; (1) non-diabetic and (2) diabetic donors without the subcutaneous tissue attachment of previous work in this field. Digital image correlation was used to assess inherent skin pre-tension of the plantar rearfoot prior to dissection. Young's modulus, storage and loss moduli were tested for using tensile stress-strain failure analysis and tensile and compressive dynamic mechanical analysis, which was conducted on excised plantar rearfoot donor specimens for both disease state cohorts at frequencies reflecting those achieved in activities of daily living. Plantar skin thickness for donor specimens were comparable to values obtained using ultrasound acquired in vivo values. Median tensile storage and loss moduli, along with Young's modulus, was higher in the diabetic cohort. With a mean Young's modulus of 0.83 ± 0.49 MPa and 1.33 ± 0.43 MPa for non-diabetic and diabetic specimens respectively. Compressive studies showed consistency between cohorts for median storage and loss moduli. The outcomes from this study show mechanical characteristics of plantar skin without the involvement of subcuteanous tissues under reflective daily achieved loading regimes, showing differences in the non-diabetic and diabetic specimens trialled to support improved understanding of plantar tissue response under tribological interactions.

A fluoroscopic imaging-guided computational analyses to inform internal tissue loads within fat pad of the diabetic foot during gait

To accurately predict internal tissue loads for early diagnostics of diabetic foot ulcerations, a novel data-driven computational analysis was conducted. A dedicated dual fluoroscopic system was combined with a pressure mat to simultaneously characterize foot motions and soft tissue's material properties during gait. Finite element (FE) models of the heel pad of a diabetic patient were constructed with 3D trajectories of the calcaneus applied as boundary conditions to simulate gait events. The tensile and compressive stresses occurring in the plantar tissue were computed. Predictions of the layered tissue FE model with anatomically-accurate heel pad structures (i.e., fat and skin) were compared with those of the traditional lumped tissue (i.e., homogeneous) models. The influence of different material properties (patient-specific versus generic) on internal tissue stresses was also investigated. The results showed the peak tensile stresses in the layered tissue model were predominantly found in the skin and distributed towards the circumferential regions of the heel, while peak compressive stresses in the fat tissue-bone interface were up to 51.4% lower than those seen in the lumped models. Performing FE analyses at four different phases of walking revealed that ignorance of layered tissue structures resulted in an unphysiological increase of peak-to-peak value of stress fluctuation in the fat and skin tissue components. Thus, to produce more clinical-relevant predictions, foot FE models are suggested to include layered tissue structures of the plantar tissue for an improved estimation of internal stresses in the diabetic foot in gait.

Disease-Specific Finite element Analysis of the Foot and Ankle

Finite-element analysis is a computational modeling technique that can be used to quantify parameters that are difficult or impossible to measure externally in a geometrically complex structure such as the foot and ankle. It has been used to improve our understanding of pathomechanics and to evaluate proposed treatments for several disorders, including progressive collapsing foot deformity, ankle arthritis, syndesmotic injury, ankle fracture, plantar fasciitis, diabetic foot ulceration, hallux valgus, and lesser toe deformities. Parameters calculated from finite-element models have been widely used to make predictions about their biomechanical correlates.

The pressure time integral: An underused, clinically significant parameter as a determinant of neuropathic ulceration in diabetes

BACKGROUND

The purpose of this study was to evaluate plantar pressure dynamics in the occurrence of active diabetic neuropathic ulceration (DNU) and any changes in loadings occurring in individuals with a history of diabetic neuropathic ulceration (DHNU). Since current gold standard offloading strategies are not producing desirable healing outcomes and optimum ulcer prevention, this study aimed to better understand the effect of diabetic peripheral neuropathy (DPN) and ulceration on mean peak plantar pressure (MPPP) and pressure-time integral (PTI) changes.

RESEARCH QUESTION

Is there a redistribution of plantar pressure during gait in the presence of active and history of diabetic neuropathic ulceration?

METHODS

A prospective, cross-sectional study was conducted, in which, eighty adult participants were divided into four groups, namely, the DM, DPN, DNU and DHNU groups. Participants were instructed to walk barefoot over a Tekscan HR Mat™ at self-selected speed. MPPP and PTI data were analysed under five forefoot anatomical landmarks and compared between individuals with and without active neuropathic ulceration.

RESULTS

Minimal MPPP significant changes were observed between ulcerated and non-ulcerated groups, however, PTI values were significantly increased in the ulcerated groups under all plantar ulceration regions. No significant plantar pressure differences were observed between the DNU and DHNU groups. Logistic regression tests demonstrated that as PTIs under the hallux increase, the likelihood of an individual living with DPN developing ulceration increases.

SIGNIFICANCE

A significant increase in PTI values in the presence of ulceration highlights the importance of evaluating the duration of loads under forefoot regions during gait rather than just focusing on the magnitude of pressures during ulcer management and prevention. Moreover, results show that PTI values remain high in the presence of a history of neuropathic ulceration, possibly demonstrating the value of PTI as a clinical tool over MPPP in the assessment of the high-risk diabetic foot.

Footwear and insole design parameters to prevent occurrence and recurrence of neuropathic plantar forefoot ulcers in patients with diabetes: a series of N-of-1 trial study protocol

BACKGROUND

Foot complications occur in conjunction with poorly controlled diabetes. Plantar forefoot ulceration contributes to partial amputation in unstable diabetics, and the risk increases with concomitant neuropathy. Reducing peak plantar forefoot pressure reduces ulcer occurrence and recurrence. Footwear and insoles are used to offload the neuropathic foot, but the success of offloading is dependent on patient adherence. This study aims to determine which design and modification features of footwear and insoles improve forefoot plantar pressure offloading and adherence in people with diabetes and neuropathy.

METHODS

This study, involving a series of N-of-1 trials, included 21 participants who had a history of neuropathic plantar forefoot ulcers. Participants were recruited from two public hospitals and one private podiatry clinic in Sydney, New South Wales, Australia. This trial is non-randomised and unblinded. Participants will be recruited from three sites, including two high-risk foot services and a private podiatry clinic in Sydney, Australia. Mobilemat™ and F-Scan® plantar pressure mapping systems by TekScan® (Boston, USA) will be used to measure barefoot and in-shoe plantar pressures. Participants' self-reports will be used to quantify the wearing period over a certain period of between 2 and 4 weeks during the trial. Participant preference toward footwear, insole design and quality-of-life-related information will be collected and analysed. The descriptive and inferential statistical analyses will be performed using IBM SPSS Statistics (version 27). And the software NVivo (version 12) will be utilised for the qualitative data analysis.

DISCUSSION

This is the first trial assessing footwear and insole interventions in people with diabetes by using a series of N-of-1 trials. Reporting self-declared wearing periods and participants' preferences on footwear style and aesthetics are the important approaches for this trial. Patient-centric device designs are the key to therapeutic outcomes, and this study is designed with that strategy in mind.

TRIAL REGISTRATION

Australian New Zealand Clinical Trials Registry (ANZCTR) ACTRN12620000699965p. Registered on June 23, 2020.

Quantification of the in vivo stiffness and natural length of the human plantar aponeurosis during quiet standing using ultrasound elastography

This study aimed to identify the stiffness and natural length of the human plantar aponeurosis (PA) during quiet standing using ultrasound shear wave elastography. The shear wave velocity (SWV) of the PA in young healthy males and females (10 participants each) was measured by placing a probe in a hole in the floor plate. The change in the SWV with the passive dorsiflexion of the metatarsophalangeal (MP) joint was measured. The Young's modulus of the PA was estimated to be 64.7 ± 9.4 kPa, which exponentially increased with MP joint dorsiflexion. The PA was estimated to have the natural length when the MP joint was plantarflexed by 13.8°, indicating that the PA is stretched by arch compression during standing. However, the present study demonstrated that the estimated stiffness for the natural length in quiet standing was significantly larger than that in the unloaded condition, revealing that the PA during standing is stiffened by elongation and through the possible activation of intrinsic muscles. Such quantitative information possibly contributes to the detailed biomechanical modeling of the human foot, facilitating an improved understanding of the mechanical functions and pathogenetic mechanisms of the PA during movements.

A novel graded-stiffness footwear device for heel ulcer prevention and treatment: a finite element-based study

Diabetic heel ulceration is a serious, destructive, and costly complication of diabetes. In this study, a novel "graded-stiffness" offloading method was proposed. This method consists of heel support with multi-increasing levels of stiffness materials, to prevent and treat heel ulcers. A three-dimensional finite element model of the heel was used to evaluate the novel "graded-stiffness" orthotic device compared to two existing solutions: (1) an insole with a hole under the active ulcer and (2) an insole with a hole filled with a soft material (elastic modulus of 15 kPa). Volumetric exposure evaluation of internal tissues to stress was performed at two volume-of-interests: (1) the area of the heel soft tissues typically at high risk for ulceration, and (2) the soft tissues surrounding the high-risk area. The models predict that the "graded-stiffness" offloading solution is more effective than existing solutions in distributing and reducing heel internal loads, considering both volume-of-interests. Comparing different material gradient combinations for the offloading support reveals considerable variation of the heel stress distribution. In clinical practice, the "graded-stiffness" technological solution enables to form an adaptable and flexible system that can be customized to a specific patient, through adequate selection of the offloading materials, to fit the shape and size of the ulcer. This solution can be made as an off-the-shelf product or alternatively, be manufactured by-demand using 3D printing tools. The proposed novel practical offloading solution has the potential for streamlining and optimizing the prevention and treatment of diabetic heel ulcers.

Development and Analysis of Key Components of a Multi Motion Mode Soft-Bodied Pipe Robot

In order to enhance the environmental adaptability of peristaltic soft-bodied pipe robots, based on the nonlinear and hyperelastic characteristics of silicone rubber combined with the biological structure and motion characteristics of worms, a hexagonal prism soft-bodied bionic actuator is proposed. The actuator adopts different inflation patterns to produce different deformations, so that the soft-bodied robot can realize different motion modes in the pipeline. Based on the Yeoh binomial parameter silicone rubber constitutive model, the deformation analysis model of the hexagonal prism soft-bodied bionic actuator is established, and the numerical simulation algorithm is used to ensure both that the drive structure and deformation mode are reasonable, and that the deformation analysis theoretical model is accurate. The motion and dynamic characteristics of the prepared hexagonal prism soft-bodied bionic actuator are tested and analyzed, the motion and dynamic characteristic curves of the actuator are obtained, and the empirical deformation formula of the actuator is fitted. The experimental results are consistent with the deformation analysis model and numerical simulation result, which shows that the deformation analysis model and numerical simulation method are accurate and can provide design methods and reference basis for the development of a pneumatic soft-bodied body bionic actuator. The above research results can also prove that the hexagonal prism soft-bodied bionic actuator is reasonable and feasible.

Different Design Feature Combinations of Flatfoot Orthosis on Plantar Fascia Strain and Plantar Pressure: A Muscle-Driven Finite Element Analysis With Taguchi Method

Customized foot orthosis is commonly used to modify foot posture and relieve foot pain for adult acquired flexible flatfoot. However, systematic investigation of the influence of foot orthotic design parameter combination on the internal foot mechanics remains scarce. This study aimed to investigate the biomechanical effects of different combinations of foot orthoses design features through a muscle-driven flatfoot finite element model. A flatfoot-orthosis finite element model was constructed by considering the three-dimensional geometry of plantar fascia. The plantar fascia model accounted for the interaction with the bulk soft tissue. The Taguchi approach was adopted to analyze the significance of four design factors combination (arch support height, medial posting inclination, heel cup height, and material stiffness). Predicted plantar pressure and plantar fascia strains in different design combinations at the midstance instant were reported. The results indicated that the foot orthosis with higher arch support (45.7%) and medial inclination angle (25.5%) effectively reduced peak plantar pressure. For the proximal plantar fascia strain, arch support (41.8%) and material stiffness (37%) were strong influencing factors. Specifically, higher arch support and softer material decreased the peak plantar fascia strain. The plantar pressure and plantar fascia loading were sensitive to the arch support feature. The proposed statistics-based finite element flatfoot model could assist the insole optimization and evaluation for individuals with flatfoot.

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.

Comparative Functional Morphology of Human and Chimpanzee Feet Based on Three-Dimensional Finite Element Analysis

To comparatively investigate the morphological adaptation of the human foot for achieving robust and efficient bipedal locomotion, we develop three-dimensional finite element models of the human and chimpanzee feet. Foot bones and the outer surface of the foot are extracted from computer tomography images and meshed with tetrahedral elements. The ligaments and plantar fascia are represented by tension-only spring elements. The contacts between the bones and between the foot and ground are solved using frictionless and Coulomb friction contact algorithms, respectively. Physiologically realistic loading conditions of the feet during quiet bipedal standing are simulated. Our results indicate that the center of pressure (COP) is located more anteriorly in the human foot than in the chimpanzee foot, indicating a larger stability margin in bipedal posture in humans. Furthermore, the vertical free moment generated by the coupling motion of the calcaneus and tibia during axial loading is larger in the human foot, which can facilitate the compensation of the net yaw moment of the body around the COP during bipedal locomotion. Furthermore, the human foot can store elastic energy more effectively during axial loading for the effective generation of propulsive force in the late stance phase. This computational framework for a comparative investigation of the causal relationship among the morphology, kinematics, and kinetics of the foot may provide a better understanding regarding the functional significance of the morphological features of the human foot.

The compressive, shear, biochemical, and histological characteristics of diabetic and non-diabetic plantar skin are minimally different

Diabetes is associated with lower limb co-morbidities, including ulceration and subsequent amputation. As a systemic disease, diabetes affects the microstructure of soft tissues, and material microstructural changes are known to affect the macroscale mechanics. However, the associations between diabetes-related disruptions to essential microstructural components and mechanical changes in plantar skin with diabetes has not been thoroughly characterized. Plantar skin specimens were collected from four diabetic and eight non-diabetic donors at six plantar locations (hallux; first, third, and fifth metatarsals; lateral midfoot; calcaneus) from matched pairs. Mechanical testing was performed on fresh frozen specimens from one foot, and histomorphological measurement and biochemical quantification were performed on specimens from the other foot. Mechanical (compressive and shear moduli and viscoelastic slopes) and biochemical/histological (total quantity of collagen and elastin; dermal and epidermal thickness) parameters were correlated using linear mixed effects regression. There were no significant differences by disease state. Skin thicknesses were positively correlated with initial compression modulus and all three shear moduli. The final compressive modulus was significantly lower at the third metatarsal than the fifth metatarsal, lateral midfoot, and calcaneus, while the final shear modulus was significantly higher at the calcaneus than at the hallux, first, and third metatarsals. Epidermal thickness was significantly higher at the calcaneus compared to all other locations. While differences were not significant by disease state, the strong differences by locations and significant but weak correlations between skin thickness and mechanics can inform future research to understand the mechanism of ulcer formation in the diabetic foot.

Prediction of diabetic foot ulcer progression: a computational study

The development of foot ulcers is a common consequence of severe diabetes. Due to vascular disorders and impeded healing caused by the disease, most foot ulcers have been reported to be affected by body weight and progress with time. Also, abnormal distribution of plantar pressures has been observed to cause the formation of additional ulcers, which may collectively lead to traumatic amputations. While a study of such pathophysiology is not possible through experiments, a few computational modelling works have investigated diabetic foot ulcers. To date, ulcers with a few sizes and locations have been studied, and their effect on the plantar stresses has been quantified. In this work, we have attempted to study the effect of all possible ulcer locations on the generated plantar peak stresses and peak stress locations where additional ulcers may form. Also, the effect of ulcer location on the possible ulcer growth was investigated. A full-scale foot model was developed and a total of 52 ulcer locations were simulated separately, with standing and walking loads. The generated stresses were normalised with the foot size and statistically analysed to develop novel formulations for predicting peak plantar stresses and their locations for any known ulcer location. The results from this study are anticipated to provide important guidelines to doctors and medical practitioners for predicting foot ulcer progression in diabetic patients with existing ulcers and allow the administration of timely preventive interventions.

The effect of induced joint restriction on plantar pressure during gait - a pilot study

BACKGROUND

The aim of this study was to analyse the effect of induced lower limb joint restriction on plantar pressures during gait. Focusing on restricting a single joint, without the effect of other co-morbidities, would provide better understanding as to the resultant plantar loadings during gait, which would be especially beneficial in patients requiring offloading procedures.

RESEARCH QUESTION

Does induced lower limb joint restriction affect plantar pressure distribution during gait?

METHODS

A prospective, quasi-experimental study was conducted, recruiting ten healthy, adult participants who were instructed to walk barefoot over a Tekscan HR Mat™. This procedure was repeated after separately inducing restriction of the hip, knee and ankle joints. Mean peak plantar pressure (MPP) and pressure-time integral (PTI) data were analysed to compare between unrestricted and restricted data.

RESULTS

Significant plantar pressure changes were observed in the heel and first metatarsal regions. Rearfoot PTIs were increased with restriction of the contralateral hip (left p <0.001) (right p =0.02) and knee joints (left p =0.01) (right p =0.04). Both MPPs (left p =0.01; right p =0.01) and PTIs (left p =0.004; right p =0.03) were increased in the first metatarsal when restricting the hip joint of the same limb. MPP was decreased in the left first metatarsal with induced knee (left p =0.01; right p =0.04) and ankle (left and right p <0.001) joint restriction. Finally, MPP was decreased in the right first metatarsal with knee (left and right p =0.01) and ankle (left p =0.04; right p =0.01) joint restriction.

SIGNIFICANCE

Limited joint mobility may have a direct effect on plantar pressure, particularly with restriction in the hip and knee joints, hence careful attention should be given especially in patients with conditions involving plantar loadings. Results in this study also show that PTI changes during gait should be equally evaluated clinically along with peak plantar pressure analysis.

Preventive strategy of flatfoot deformity using fully automated procedure

A non-invasive, no radiation, out-of-hospital automated system is proposed to identify low arch integrated in the design and manufacturing of personalized orthoses using parametric modelling. The aim of the design process is to integrate assistive technology with assessment and prevent low arch progressing to a more serious case - flatfoot. In the automated procedure, we developed an assessment method including reliable thresholds of foot type classification and test protocol to reduce interferences due to preceding activities, an automation to translate scanned data into parametric design for orthotic customization, finite element model evaluating effectiveness of the personalized design, and a personalized comparative test to evaluate the long-term improvement of foot arch shape. Our low arch threshold established by subject-specific 3D models reduced the misclassification rate from 55%, as previously reported to 6.9%. Individuals who engaged in sedentary activity (i.e. sitting) had the greater change in arch height compared to active activity (i.e. standing and walking), which is more likely to affect the obtained measure. Therefore, a test protocol now states that participants are not allowed to sit over 100 min prior the measurement to reduce such interference. We have proposed and tested an automated algorithm to translate scanned data including seven foot's parameters into customised parametric design of the insert. The method decreases the required time of orthotic computer-aided design from over 3 h to less than 2 min. A finite element analysis procedure was additionally developed to assess the performance of geometries and material of designed orthotic based on the distribution of plantar pressure and internal stress. The personalized comparative assessment based on midfoot contact area was carried out periodically for follow-up and the orthotic could be re-designed if necessary. The proposed automated procedure develops a pre-screening system to distinguish low arch and provide preventatives before it becomes symptomatic. Furthermore, non-symptom flatfoot can be detected at early stages and referred to medics for further diagnosis or treatment.

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.

Three-Dimensional Finite Element Analysis and Biomechanical Analysis of Midfoot von Mises Stress Levels in Flatfoot, Clubfoot, and Lisfranc Joint Injury

Background

Midfoot deformity and injury can affect the internal pressure distribution of the foot. This study aimed to use 3D finite element and biomechanical analyses of midfoot von Mises stress levels in flatfoot, clubfoot, and Lisfranc joint injury.

Material/Methods

Normal feet, flatfeet, clubfeet (30 individuals each), and Lisfranc injuries (50 individuals) were reconstructed by CT, and 3D finite element models were established by ABAQUS. Spring element was used to simulate the plantar fascia and ligaments and set hyperelastic coefficients in encapsulated bone and ligaments. The stance phase was simulated by applying 350 N on the top of the talus. The von Mises stress of the feet and ankle was visualized and analyzed.

Results

The von Mises stress on healthy feet was higher in the lateral metatarsal and ankle bones than in the medial metatarsal bone. Among the flatfoot group, the stress on the metatarsals, talus, and navicular bones was significantly increased compared with that on healthy feet. Among patients with clubfeet, stress was mainly concentrated on the talus, and stress on the lateral metatarsal and navicular bones was significantly lower. The von Mises stress on the fractured bone was decreased, and the stress on the bone adjacent to the fractured bone was higher in Lisfranc injury. During bone dislocation alone or fracture accompanied by dislocation, the von Mises stress of the dislocated bone tended to be constant or increased.

Conclusions

Prediction of von Mises stress distribution may be used clinically to evaluate the effects of deformity and injury on changes in structure and internal pressure distribution on the midfoot.

OPTIMAL MESH CRITERIA IN FINITE ELEMENT MODELING OF HUMAN FOOT: THE DEPENDENCE FOR MULTIPLE MODEL OUTPUTS ON MESH DENSITY AND LOADING BOUNDARY CONDITIONS

The use of finite element models has gained popularity in the field of foot and footwear biomechanics to predict the stress–strain distribution and the treatment effectiveness of therapeutic insoles for pathological foot conditions. However, a comprehensive evaluation of mesh quality is often ignored, meanwhile no golden standard exists for the mesh density and selection of element size at an acceptable accuracy. Here, we make a convergence test and established anatomically-realistic foot models at different mesh densities. The study compared the discrepancy in output variables to the changes of element type and mesh density under barefoot and footwear conditions with compressive and shear loads, which are commonly encountered in foot and footwear biomechanics simulations. For a range of loading conditions simulated in 125 finite element models, the peak plantar pressure consistently converged with optimal mesh size determined at 2.5[Formula: see text]mm. The convergence variable of principal strains and stress tensors, however, varies significantly. The max von-Mises stress showed strong sensitive behavior to the changes of the mesh density. The pattern for contact pressure distribution became less accurate when the element sizes increase to 6.0[Formula: see text]mm; in particular, the locations of the pressure peak do not show remarkable changes, but the size of the area of contact still changes. The current study could offer a general guideline when generating a reasonable accurate finite element models for the analysis of plantar pressure distributions and stress/strain states employed for foot and footwear biomechanics evaluations.

Influence of arch support heights on the internal foot mechanics of flatfoot during walking: A muscle-driven finite element analysis

BACKGROUND

Different arch support heights of the customized foot orthosis could produce different effects on the internal biomechanics of the foot. However, quantitative evidence is scarce. Therefore, we aimed to investigate and quantify the influence of arch support heights on the internal foot biomechanics during walking stance.

METHODS

We reconstructed a foot finite element model from a volunteer with flexible flatfoot. The model enabled a three-dimensional representation of the plantar fascia and its interactions with surrounding osteotendinous structures. The volunteer walked in foot orthosis with different arch heights (low, neutral, and high). Muscle forces during gaits were calculated by a multibody model and used to drive a foot finite element model. The foot contact pressures and plantar fascia strains in different regions were compared among the insole conditions at the first and second vertical ground reaction force (VGRF) peak and VGRF valley instants.

RESULTS

The results indicated that peak foot pressures decreased in balanced standing and second VGRF as the arch support height increased. However, peak midfoot pressures increased during all simulated instants. Meanwhile, high arch support decreased the plantar fascia loading by 5%-15.4% in proximal regions but increased in the middle and distal regions.

CONCLUSION

Although arch support could generally decrease the plantar foot pressure and plantar fascia loading, the excessive arch height may induce high midfoot pressure and loadings at the central portion of the plantar fascia. The consideration of fascia-soft tissue interaction in modeling could improve the prediction of plantar fascia strains towards design optimization for orthoses.

Development and Performance Analysis of Pneumatic Soft-Bodied Bionic Actuator

The design of a pneumatic soft-bodied bionic actuator derives from the structural characteristics and motion mechanism of biological muscles, combined with the nonlinear hyperelasticity of silica gel, which can improve the mobility and environmental adaptability of soft-bodied bionic robots. Based on Yeoh's second-order constitutive model of silica gel, the deformation analysis model of the actuator is established, and the rationality of the structure design and motion forms of the actuator and the accuracy of the deformation analysis model are verified by using the numerical simulation algorithm. According to the physical model of the pneumatic soft-bodied bionic actuator, the motion and dynamic characteristics of the actuator are tested and analyzed, the curves of motion and dynamic characteristics of the actuator are obtained, and the empirical formula of the bending angle and driving torque of the actuator is fitted out. The results show that the deformation analysis model and numerical simulation method are accurate, and the pneumatic soft-bodied bionic actuator is feasible and effective, which can provide a design method and reference basis for the research and implementation of soft-bodied bionic robot actuator.

Development and Performance Analysis of Pneumatic Soft-Bodied Bionic Basic Execution Unit

This paper studies the design of pneumatic soft-bodied bionic basic execution unit with soft-rigid combination, which can be used as an actuator for pneumatic soft-bodied robots and soft-bodied manipulators. This study is inspired by structural characteristics and motion mechanism of biological muscles, combined with the nonlinear hyperelasticity of silica gel and the insertion of thin leaf spring structure in the nonretractable layer. Response surface analysis and numerical simulation algorithm are used to determine the optimal combination of structural dimension parameters by taking the maximum output bending angle of the basic executing unit as the optimization objective. Based on Odgen’s third-order constitutive model, the deformation analysis model of the basic execution unit is established. The physical model of pneumatic soft-bodied bionic basic execution unit is prepared through 3D printing, shape deposition, soft lithography, and other processing methods. Finally, the motion and dynamic characteristics of the physical model are tested through experiments and result analysis, thus obtaining curves and empirical formulas describing the motion and dynamic characteristics of the basic execution unit. The relevant errors are compared with the deformation analysis model of the execution unit to verify the feasibility and effectiveness of the design of the pneumatic soft-bodied bionic basic execution unit. The above research methods, research process, and results can provide a reference for the research and implementation of pneumatic and hydraulic driven soft-bodied robots and grasping actuators of soft-bodied manipulators.

Developing an Equivalent Solid Material Model for BCC Lattice Cell Structures Involving Vertical and Horizontal Struts

In this study, a body-centered cubic (BCC) lattice unit cell occupied inside a frame structure to create a so-called “InsideBCC” is considered. The equivalent quasi-isotropic properties required to describe the material behavior of the InsideBCC unit cell are equivalent Young’s modulus ( E e ) , equivalent shear modulus ( G e ) , and equivalent Poisson’s ratio ( ν e ) . The finite element analysis (FEA) based computational approach is used to simulate and calculate the mechanical responses of InsideBCC unit cell, which are the mechanical responses of the equivalent solid. Two separates finite element models are then developed for samples under compression: one with a 6 × 6 × 4 cell InsideBCC lattice cell structure (LCS) and one completely solid with equivalent solid properties obtained from a unit cell model. In addition, 6 × 6 × 4 cell specimens are fabricated on a fused deposition modeling (FDM) uPrint SEplus 3D printer using acrylonitrile butadiene styrene (ABS) material and tested experimentally under quasi-static compression load. Then, the results extracted from the finite element simulation of both the entire lattice and the equivalent solid models are compared with the experimental data. A good agreement between the experimental stress–strain behavior and that obtained from the FEA models is observed within the linear elastic limit.

Footwear and insole design features that reduce neuropathic plantar forefoot ulcer risk in people with diabetes: a systematic literature review

BACKGROUND

In people with diabetes, offloading high-risk foot regions by optimising footwear, or insoles, may prevent ulceration. This systematic review aimed to summarise and evaluate the evidence for footwear and insole features that reduce pathological plantar pressures and the occurrence of diabetic neuropathy ulceration at the plantar forefoot in people with diabetic neuropathy.

METHODS

Six electronic databases (Medline, Cinahl, Amed, Proquest, Scopus, Academic Search Premier) were searched in July 2019. The search period was from 1987 to July 2019. Articles, in English, using footwear or insoles as interventions in patients with diabetic neuropathy were reviewed. Any study design was eligible for inclusion except systematic literature reviews and case reports. Search terms were diabetic foot, physiopathology, foot deformities, neuropath*, footwear, orthoses, shoe, footwear prescription, insole, sock*, ulcer prevention, offloading, foot ulcer, plantar pressure.

RESULTS

Twenty-five studies were reviewed. The included articles used repeated measure (n = 12), case-control (n = 3), prospective cohort (n = 2), randomised crossover (n = 1), and randomised controlled trial (RCT) (n = 7) designs. This involved a total of 2063 participants. Eleven studies investigated footwear, and 14 studies investigated insoles as an intervention. Six studies investigated ulcer recurrence; no study investigated the first occurrence of ulceration. The most commonly examined outcome measures were peak plantar pressure, pressure-time integral and total contact area. Methodological quality varied. Strong evidence existed for rocker soles to reduce peak plantar pressure. Moderate evidence existed for custom insoles to offload forefoot plantar pressure. There was weak evidence that insole contact area influenced plantar pressure.

CONCLUSION

Rocker soles, custom-made insoles with metatarsal additions and a high degree of contact between the insole and foot reduce plantar pressures in a manner that may reduce ulcer occurrence. Most studies rely on reduction in plantar pressure measures as an outcome, rather than the occurrence of ulceration. There is limited evidence to inform footwear and insole interventions and prescription in this population. Further high-quality studies in this field are required.

Experimental and finite element investigation of total ankle replacement: A review of literature and recommendations

This paper briefly reviews the different methodology, technology, challenges, and outcomes of various studies related to TAR prosthesis based on numerical and experimental techniques. Very less in-vitro experimental studies on TAR have been found than finite element (FE) studies. Due to the invasive nature of the experimental approach, inadequacy and less clinical information, computational modelling has been widely used by the researchers. This paper critically examines the part related to FE modelling and experimental analysis. Some recommendation related to modelling of bones, cartilages, ligaments, muscles, and implant-bone interface condition were discussed for better understanding the results and better clinical significance.

Implementation of controlling strategy in a biomechanical lower limb model with active muscles for coupling multibody dynamics and finite element analysis

Computational biomechanics for human body modeling has generally been categorized into two separated domains: finite element analysis and multibody dynamics. Combining the advantages of both domains is necessary when tissue stress and physical body motion are both of interest. However, the method for this topic is still in exploration. The aim of this study is to implement unique controlling strategies in finite element model for simultaneously simulating musculoskeletal body dynamics and in vivo stress inside human tissues. A finite element lower limb model with 3D active muscles was selected for the implementation of controlling strategies, which was further validated against in-vivo human motion experiments. A unique feedback control strategy that couples together a basic Proportion-Integration-Differentiation (PID) controller and generic active signals from Computed Muscle Control (CMC) method of the musculoskeletal model or normalized EMG singles was proposed and applied in the present model. The results show that the new proposed controlling strategy show a good correlation with experimental test data of the normal gait considering joint kinematics, while stress distribution of local lower limb tissue can be also detected in real-time with lower limb motion. In summary, the present work is the first step for the application of active controlling strategy in the finite element model for concurrent simulation of both body dynamics and tissue stress. In the future, the present method can be further developed to apply it in various fields for human biomechanical analysis to monitor local stress and strain distribution by simultaneously simulating human locomotion.

Functional gradient structural design of customized diabetic insoles

Diabetic foot is a common and serious complication of diabetes, largely due to sensory neuropathy and excessive mechanical stresses. Studies have shown that reducing the contact pressure can effectively lower the incidence of diabetic foot. A new design method is proposed in this study for optimizing the stress distribution of the contact surfaces between the foot and the insole by applying functional gradient structural properties to the insole. Finite element analysis was employed for studying the contact mechanics, which laid the foundation for modulus readjustment during the optimization process. The moduli of the materials were correlated to the properties of the structural porous units. The customized insoles were manufactured using additive manufacturing technology and put into mechanical test. Results show that the designed insole helps in increasing the foot contact area by approximately 30% and reducing the peak contact pressure by 35%. Hence, the proposed method can be used to design customized insoles, particularly diabetic insoles, by offering better contact mechanics and good potential for reducing the severity of diabetic foot. The methodology is equally applicable to other designs involving optimization of material properties.

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.

A methodological framework for detecting ulcers' risk in diabetic foot subjects by combining gait analysis, a new musculoskeletal foot model and a foot finite element model

Diabetic foot is one of the most debilitating complications of diabetes and may lead to plantar ulcers. In the last decade, gait analysis, musculoskeletal modelling (MSM) and finite element modelling (FEM) have shown their ability to contribute to diabetic foot prevention and suggested that the origin of the plantar ulcers is in deeper tissue layers rather than on the plantar surface. Hence the aim of the current work is to develop a methodology that improves FEM-derived foot internal stresses prediction, for diabetic foot prevention applications. A 3D foot FEM was combined with MSM derived force to predict the sites of excessive internal stresses on the foot. In vivo gait analysis data, and an MRI scan of a foot from a healthy subject were acquired and used to develop a six degrees of freedom (6 DOF) foot MSM and a 3D subject-specific foot FEM. Ankle kinematics were applied as boundary conditions to the FEM together with: 1. only Ground Reaction Forces (GRFs); 2. OpenSim derived extrinsic muscles forces estimated with a standard OpenSim MSM; 3. extrinsic muscle forces derived through the (6 DOF) foot MSM; 4. intrinsic and extrinsic muscles forces derived through the 6 DOF foot MSM. For model validation purposes, simulated peak pressures were extracted and compared with those measured experimentally. The importance of foot muscles in controlling plantar pressure distribution and internal stresses is confirmed by the improved accuracy in the estimation of the peak pressures obtained with the inclusion of intrinsic and extrinsic muscle forces.

Bone cyst formation after ankle arthroplasty may be caused by stress shielding. A numerical simulation of the strain adaptive bone remodelling

BACKGROUND

The history of total ankle arthroplasty (TAA) has different evolution steps to improve the outcome. The third generation implants show an overall 8-year survival rate up to 93%. The main reported reason for early failure of TAA is aseptic loosening, cyst formation is also frequently reported. The aim of the present study is to use the finite element (FE) method to analyze the adaptive bone remodeling processes, including cyst formation after TAA.

METHODS

Bone characteristics applied to the model corresponded to information obtained from computed tomography. Finite element models for the tibia and the talus were developed and implant components were virtually implanted.

RESULTS

The calculated total bone loss is 2% in the tibia and 17% in the talus. Cysts and areas of increased bone density were detectable dependent on prosthesis design in the tibia and talus.

CONCLUSION

Our FE simulation provides a theoretical explanation for cyst formation and increasing bone density depending on implant design. However, cysts are not mono-causal, histo-chemical reactions should also be considered. Further clinical studies are necessary to evaluate the relevance of cyst formation and therapeutic strategies.

Subject-specific finite element modelling of the human foot complex during walking: sensitivity analysis of material properties, boundary and loading conditions

The objective of this study was to develop and validate a subject-specific framework for modelling the human foot. This was achieved by integrating medical image-based finite element modelling, individualised multi-body musculoskeletal modelling and 3D gait measurements. A 3D ankle-foot finite element model comprising all major foot structures was constructed based on MRI of one individual. A multi-body musculoskeletal model and 3D gait measurements for the same subject were used to define loading and boundary conditions. Sensitivity analyses were used to investigate the effects of key modelling parameters on model predictions. Prediction errors of average and peak plantar pressures were below 10% in all ten plantar regions at five key gait events with only one exception (lateral heel, in early stance, error of 14.44%). The sensitivity analyses results suggest that predictions of peak plantar pressures are moderately sensitive to material properties, ground reaction forces and muscle forces, and significantly sensitive to foot orientation. The maximum region-specific percentage change ratios (peak stress percentage change over parameter percentage change) were 1.935-2.258 for ground reaction forces, 1.528-2.727 for plantar flexor muscles and 4.84-11.37 for foot orientations. This strongly suggests that loading and boundary conditions need to be very carefully defined based on personalised measurement data.

Validation of plantar pressure simulations using finite and discrete element modelling in healthy and diabetic subjects

Plantar pressure simulation driven by integrated 3D motion capture data, using both a finite element and a discrete element model, is compared for ten healthy and ten diabetic neuropathic subjects. The simulated peak pressure deviated on average between 16.7 and 34.2% from the measured peak pressure. The error in the position of the peak pressure was on average smaller than 4.2 cm. No method was more accurate than the other although statistical differences were found between them. Both techniques are thus complementary and useful tools to better understand the alteration of diabetic foot biomechanics during gait.

Strategies towards rapid generation of forefoot model incorporating realistic geometry of metatarsals encapsulated into lumped soft tissues for personalized finite element analysis

Use of finite element (FE) foot model as a clinical diagnostics tool is likely to improve the specificity of foot injury predictions in the diabetic population. Here we proposed a novel workflow for rapid construction of foot FE model incorporating realistic geometry of metatarsals encapsulated into lumped forefoot's soft tissues. Custom algorithms were implemented to perform unsupervised segmentation and mesh generation to directly convert CT data into a usable FE model. The automatically generated model provided higher efficiency and comparable numerical accuracy when compared to the model constructed using a traditional solid-based mesh process. The entire procedure uses MATLAB as the main platform, and makes the present approach attractive for creating personalized foot models to be used in clinical studies.

Footwear and insole design features to prevent foot ulceration in people with diabetes: a systematic review protocol

REVIEW QUESTION/OBJECTIVE

The aim of this systematic review is to identify the key design features of footwear and insoles that are used to offload the plantar surface of the foot to prevent foot ulceration in adults with diabetes.More specifically, the objectives are to identify the key design features of footwear and insoles to offload the plantar surface of the foot with regard to.

Biomechanical study on surgical fixation methods for minimally invasive treatment of hallux valgus

Hallux valgus (HV) was one of the most frequent female foot deformities. The aim of this study was to evaluate mechanical responses and stabilities of the Kirschner, bandage and fiberglass fixations after the distal metatarsal osteotomy in HV treatment. Surface traction of different forefoot regions in bandage fixation and the biomechanical behavior of fiberglass bandage material were measured by a pressure sensor device and a mechanical testing, respectively. A three-dimensional foot finite element (FE) model was developed to simulate the three fixation methods (Kirschner, bandage and fiberglass fixations) in weight bearing. The model included 28 bones, sesamoids, ligaments, plantar fascia, cartilages and soft tissue. The peak Von Mises stress (MS) and compression stress (CS) of the distal fragment were predicted from the three fixation methods: Kirschner fixation (MS=6.71MPa, CS=1.232MPa); Bandage fixation (MS=14.90MPa, CS=9.642MPa); Fiberglass fixation (MS=15.83MPa, CS=19.70MPa). Compared with the Kirschner and bandage fixation, the fiberglass fixation reduced the relative movement of osteotomy fragments and obtained the maximum CS. We concluded that fiberglass fixation in HV treatment was helpful to the bone healing of distal fragment. The findings were expected to guide further therapeutic planning of HV patient.

Quantifying Dynamic Changes in Plantar Pressure Gradient in Diabetics with Peripheral Neuropathy

Diabetic foot ulcers remain one of the most serious complications of diabetes. Peak plantar pressure (PPP) and peak pressure gradient (PPG) during walking have been shown to be associated with the development of diabetic foot ulcers. To gain further insight into the mechanical etiology of diabetic foot ulcers, examination of the pressure gradient angle (PGA) has been recently proposed. The PGA quantifies directional variation or orientation of the pressure gradient during walking and provides a measure of whether pressure gradient patterns are concentrated or dispersed along the plantar surface. We hypothesized that diabetics at risk of foot ulceration would have smaller PGA in key plantar regions, suggesting less movement of the pressure gradient over time. A total of 27 participants were studied, including 19 diabetics with peripheral neuropathy and 8 non-diabetic control subjects. A foot pressure measurement system was used to measure plantar pressures during walking. PPP, PPG, and PGA were calculated for four foot regions – first toe (T1), first metatarsal head (M1), second metatarsal head (M2), and heel (HL). Consistent with prior studies, PPP and PPG were significantly larger in the diabetic group compared with non-diabetic controls in the T1 and M1 regions, but not M2 or HL. For example, PPP was 165% (P = 0.02) and PPG was 214% (P < 0.001) larger in T1. PGA was found to be significantly smaller in the diabetic group in T1 (46%, P = 0.04), suggesting a more concentrated pressure gradient pattern under the toe. The proposed PGA may improve our understanding of the role of pressure gradient on the risk of diabetic foot ulcers.