In vivo friction properties of human skin

Extended finite element analysis for the 2nd and 3rd metatarsals stress fracture during push-off

Metatarsal stress fractures (MSF), particularly the 2nd and 3rd MSF, are common injuries among athletes. Although there are several practices to reduce foot and ankle injuries, there is no injury prevention program specifically designed to minimize MSF. This is mainly due to the lack of information about the loadings/postures that cause MSF. Therefore, this study aimed to investigate dangerous loadings/postures potentially causing MSF during push-off (PO). The analysis was conducted with Finite Element Modelling (FEM), calibrated with the three-point bending test, and validated with peak plantar pressure (PPP) and fracture force measurement. Extended Finite Element Method was used for MSF simulation such that ten different foot and ankle configurations were designed, with five for each of the 2nd and 3rd MSF under pure vertical loadings. A more complex loading, ankle eversion/inversion during PO, was also examined for the MSF. The average error percentage for the calibration of the model with the three-point bending test was 3.05%. The average error percentages for the validation of the model with PPP and fracture force measurements were 18% and 30%, respectively. The outcomes of pure vertical loadings indicated the higher potential for the 2nd and 3rd MSF at 30% PO and 70% PO, respectively. The results of ankle eversion/inversion loadings represented that the most dangerous posture for MSF was 30° ankle eversion for the 3rd metatarsal at 70% PO. These results provide a guide, including what postures to avoid for the 2nd and 3rd MSF among people who are at high risk of MSF.

Custom orthotic design by integrating 3D scanning and subject-specific FE modelling workflow

The finite element (FE) foot model can help estimate pathomechanics and improve the customized foot orthoses design. However, the procedure of developing FE models can be time-consuming and costly. This study aimed to develop a subject-specific scaled foot modelling workflow for the foot orthoses design based on the scanned foot surface data. Six participants (twelve feet) were collected for the foot finite element modelling. The subject-specific surface-based finite element model (SFEM) was established by incorporating the scanned foot surface and scaled foot bone geometries. The geometric deviations between the scaled and the scanned foot surfaces were calculated. The SFEM model was adopted to predict barefoot and foot-orthosis interface pressures. The averaged distances between the scaled and scanned foot surfaces were 0.23 ± 0.09 mm. There was no significant difference for the hallux, medial forefoot, middle forefoot, midfoot, medial hindfoot, and lateral hindfoot, except for the lateral forefoot region (p = 0.045). The SFEM model evaluated slightly higher foot-orthoses interface pressure values than measured, with a maximum deviation of 7.1%. These results indicated that the SFEM technique could predict the barefoot and foot-orthoses interface pressure, which has the potential to expedite the process of orthotic design and optimization.

Determining the changes in morphology and loading status following medial displacement calcaneal osteotomy for flatfoot using patient-specific finite element models

Medial displacement calcaneal osteotomy (MDCO) is the standard procedure for flatfoot. We investigated the effect of MDCO on the foot using a finite element analysis. Foot models were created from computed tomography data of 8 patients with flat feet. MDCO was performed on each model with bone translation distance of 4, 8, and 12 mm. The morphological changes, plantar pressures, and stress percentage on the talocrural and subtalar joints were evaluated before and after surgery. Morphological evaluation showed improvement in the medial longitudinal arch. The stress percentage of plantar pressure in the medial area decreased, and the stress percentage of plantar pressure in the mid- and lateral forefoot area increased. At the talocrural joint, the medial and middle stress percentage increased, while the lateral and posterior stress percentage decreased. In the subtalar joint, the stress percentage in the middle subtalar joint increased and that in the posterior subtalar joint decreased. Within the posterior subtalar joint, the anterior and medial stress percentage increased, while the posterior and lateral stress percentage decreased. Preoperative simulation using the finite element analysis may be useful in understanding postoperative morphological changes and loading conditions to perform patient-specific surgery.

Method for evaluating the fit of custom-fit Continuous Positive Airway Pressure masks using finite element analysis

Continuous Positive Airway Pressure (CPAP) is a common therapy used to treat breathing disorders such as obstructive sleep apnea. In previous work, we designed a custom-fit CPAP mask to address comfort and leakage issues patients often experience. This paper presents a method to create a finite element (FE) model to evaluate the fit of the custom-fit mask before fabrication. The model includes details such as a skull to represent the variable soft tissue thicknesses on the face, and two strap configurations, original and X. The model was tested on four subjects and results show that the X strap configuration results in a more even stress distribution, measured by standard deviation, on the face compared to the original strap, indicating better fit. The simulations also show gaps in the stress distribution that seem to correspond to areas of leakage based on two initial in vivo tests on two subjects. This simulation method proves to be a valuable tool for custom-fit mask development by allowing us to evaluate designs before fabrication.

A pneumatic reconfigurable socket for transtibial amputees

Many transtibial amputees rate the fit between their residual limb and prosthetic socket as the most critical factor in satisfaction with using their prosthesis. This study aims to address the issue of prosthetic socket fit by reconfiguring the socket shape at the interface of the residual limb and socket. The proposed reconfigurable socket shifts pressure from sensitive areas and compensates for residual limb volume fluctuations, the most important factors in determining a good socket fit. Computed tomography scan images are employed to create the phantom limb of an amputee and to manufacture the reconfigurable socket. The performance of the reconfigurable socket was evaluated both experimentally and numerically using finite element modelling. The study showed that the reconfigurable socket can reduce interface pressure at targeted areas by up to 61%.

Friction Blisters of the Feet: A New Paradigm to Explain Causation

Friction blisters on the feet commonly occur when individuals engage in active pursuits such as running, hiking, and military training. The high prevalence of blisters in active individuals underscores the fact that the pathomechanics of this condition are not fully understood. The traditional blister causation paradigm revolves around heat, moisture, and friction. In reality, foot friction blisters are caused by repetitive shear deformation. The 3 fundamental elements of blister-inducing shear deformation are (1) motion of bone, (2) high friction force, and (3) repetition of the resulting shear events. Rubbing at the skin surface is not a mechanism for friction blister formation. To that end, prevention of the friction blister continues to be an elusive quest for both the patient and the treating clinician. In this article, we aimed to highlight the limitations of the long-held blister-causation paradigm and offer a new explanation.

Body mapping of skin friction coefficient and tactile perception during the dynamic skin-textile interaction

The clothing fabric and skin interact continuously across the many regions of users' bodies during wear, which can lead to both physical skin damage and discomfort. Therefore, this investigation aimed to explore the regional differences in skin friction, tactile perception, and sensitivity in both females and males during the skin-textile interaction. The static and dynamic friction coefficient and textile perception (texture, stickiness, pleasantness, and discomfort) were measured across the 36 selected testing body areas by using a friction measurement device. The results revealed there was a significant difference in skin friction, tactile perceptions, and sensitivity across the various body regions. The anterior neck had the highest skin friction in both females and males, and participants generally rated higher texture perception in their anterior aspects compared to posterior and lateral regions. There was no significant difference in skin friction, tactile perception ratings, and sensitivity between females and males. Practitioner summary: This study sought to examine regional variations in skin friction, tactile perception, and sensitivity during the skin-textile interaction. There was a significant difference in skin friction, tactile perceptions, and sensitivity across the various body regions and no significant sex effect on skin friction, tactile perception ratings, and sensitivity.

A flexible omnidirectional rotating magnetic array for MRI-safe transdermal wireless energy harvesting through flexible electronics

Magnetic resonance imaging (MRI)-safe implantable wireless energy harvester offers substantial benefits to patients suffering from brain disorders, hearing impairment, and arrhythmias. However, rigid magnets in cutting-edge systems with limited numbers of rotation axis impose high risk of device dislodgement and magnet failure. Here, a flexible omnidirectional rotating magnetic array (FORMA) and a flexible MRI-safe implantable wireless energy-harvesting system have been developed. Miniaturized flexible magnetic balls 1 millimeter in diameter achieved by molding three-dimensional printed templates can rotate freely in elastomer cavities and supply a magnetic force of 2.14 Newtons at a distance of 1 millimeter between an implantable receiver and a wearable transceiver. The system can work stably under an acceleration of 9g and obtain a power output of 15.62 decibel milliwatts at a transmission frequency of 8 megahertz. The development of the FORMA may lead to life-long flexible and batteryless implantable systems and offers the potential to promote techniques for monitoring and treating acute and chronic diseases.

Developing failure criteria for laceration injury of dermal tissue

Despite its importance, there is a poor understanding of human injury tolerance to trauma generally, and more specifically understanding of the mechanics of skin penetration or laceration. The objective of this analysis is to determine the failure criteria that will allow the evaluation of the laceration risk of blunt-tipped edges within a computational modeling environment. An axisymmetric tissue finite element model was set up in Abaqus 2021 to match the experimental set-up from a previous study. The model simulated the pressing of penetrometer geometries into dermal tissue, and stress and strain outputs were evaluated at the experimental failure force. Two separate nonlinear hyperelastic material models were calibrated for the dermis to data from the literature (high and low stiffness models). For both the high-stiffness and low-stiffness skin models, the failure force appears to occur near a local maximum in the principal strain. All failures occurred after the maximum strain near or at the top surface is or above 59%, with mid-thickness strain at a similar level. The strain energy density is concentrated near the edge tip for each configuration, indicating highly localized material damage at the point of loading, and increases rapidly prior to the approximate failure force. As the edge is further compressed into the tissue, the stress triaxiality near the edge contacting point decreases towards zero. This study has identified general failure criteria for skin laceration which can be implemented in a computational model. A higher risk for laceration would be indicated with strain energy density larger than 60 mJ/mm3, dermal strain larger than 55%, and stress triaxiality below 0.1. These findings were largely insensitive to the dermal stiffness and broadly applicable across different indenter geometries. It is expected that this framework may be implemented to evaluate hazardous forces for product edges, interactions with robots, and interfaces with medical and drug delivery devices.

Three dimensional finite element analysis of biomechanics of osteotomy ends with three different fixation methods after hallux valgus minimally invasive osteotomy

PURPOSE

Biomechanical study of fixation methods post hallux valgus minimally invasive osteotomy using finite element technology hasn't been reported. This study aimed to compare maximum displacement and stress distribution of osteotomy ends after minimally invasive osteotomy fixed by bandage, Kirschner wire, Herbert screw.

METHODS

Foot CT images of a patient with mild-moderate hallux valgus were collected. Three-dimensional finite element model of hallux valgus was established through CT image. This study simulated bandage, Kirschner wire and Herbert screw fixation, and analyzed maximum displacement and stress distribution of osteotomy ends in plantar flexion position of foot after fixation.

RESULTS

Maximum equivalent stress of osteotomy end fixed with bandage, Kirschner wire, Herbert screw was 7.8615, 14.253, 8.3156 MPa, respectively. Total displacement of osteotomy end fixed by bandage, Kirschner wire, Herbert screw was 0.26,896, 0.022,779, 0.029,195 mm, respectively. Maximum stress of Kirschner wire and Herbert screw near osteotomy end was 154.7 and 46.404 MPa, respectively. Fixation strength and stability of Kirschner wire and Herber screw were better than bandage. Kirschner wire had stress concentration phenomenon, with potential fracture risk. Stress of Herbert screw was evenly distributed around osteotomy end, and there was a certain stress concentration, playing an important role in maintaining fracture end stability.

CONCLUSIONS

Herbert screw showed good fixation strength and stability, and stress distribution was uniform, which can well maintain stability of minimally invasive osteotomy ends. Findings of this study would provide a theoretical basis for selection of fixation methods after clinical minimally invasive osteotomy for hallux valgus.

Coupled Diffusion-Binding-Deformation Modelling for Phase-Transition Microneedles-Based Drug Delivery

Phase-transition microneedles (PTMNs)-based transdermal drug delivery (TDD) is gaining popularity due to its non-invasiveness and ability to deliver a wide range of drugs. PTMNs absorb interstitial skin fluid (ISF) and transport drugs from microneedle (MNs) domain to the skin without polymer dissolution. To establish PTMNs for practical use, one needs to understand and optimise the key parameters governing drug transport mechanisms to achieve controlled drug delivery. In addressing this point, we have developed a coupled diffusion-binding-deformation model to understand the effect of physicochemical parameters (e.g., swelling capacity, drug binding) of MN and skin mechanical properties on overall drug transport behaviour. The contact mechanics at the MN and skin interface is introduced to account for the resistive force exerted by the deformed skin to MN swelling. The model is validated with the reported data of in vitro insulin delivery using polyvinyl alcohol (PVA) MN. The drug binding parameters are estimated from the fitting of the cumulative release of insulin within 6 hours of MN insertion. To predict the in vivo data of insulin delivery using the PVA MN, one-compartment model of drug pharmacokinetics is incorporated. It is shown in the paper that the model is able to predict the final insulin concentration in blood and in good agreement with the reported experimental data. The proposed model is concluded to be a tool for the predictive design and development of PTMNs-based TDD systems.

The Impact of Different Overlay Materials on the Tactile Detection of Virtual Straight Lines

To improve the perception of haptic feedback, materials and sense-modifier effects should be examined. Teflon, Nylon mesh, and Silicone overlays were tested in combination with lateral vibrations to study their impact on the tactile sense. A feelable point moving along a line was implemented through the use of a dynamically moving interference maximum generated via the offset actuation of four haptic exciters affixed to corners of a Gorilla Glass surface. This feedback was presented to eight participants in a series of randomized experiments. Both the Nylon mesh and Teflon covering revealed a statistically significant (p < 0.05) impact of improvement to the user performance in the task of dynamic haptic virtual straight lines localization. While Silicone covering, having three times greater friction than Gorilla Glass, has less or no impact on both decision time, the number of task repetitions, and error rate (p > 0.05). The lateral vibration modifier (60 Hz) can also successfully be used with an increase in performance by about twofold, at least that was demonstrated for both the Nylon mesh and Teflon covering.

Numerical simulation of mechanical tests on a living skin using anisotropic hyperelastic law

The skin is a living tissue that behaves in a hyperelastic and anisotropic way. A constitutive law called HGO-Yeoh is proposed to model the skin by improving the classical HGO constitutive law. This model is implemented in a finite element code FER "Finite Element Research" to benefit from its tools, including the bipotential contact method, a very efficient function coupling contact and friction. Identifying the skin-related material parameters is done through an optimisation procedure using analytic and experimental data. A tensile test is simulated using the codes FER and ANSYS. Then, the results are compared with the experimental data. Finally, a simulation of an indentation test using a bipotential contact law is done.

A review of the characterizations of soft tissues used in human body modeling: Scope, limitations, and the path forward

Soft tissue material properties are vital to human body models that evaluate interactions between the human body and its environment. Such models evaluate internal stress/strain responses in soft tissues to investigate issues like pressure injuries. Numerous constitutive models and parameters have been used to represent mechanical behavior of soft tissues in biomechanical models under quasi-static loading. However, researchers reported that generic material properties cannot accurately represent specific target populations due to large inter-individual variability. Two challenges that exist are experimental mechanical characterization and constitutive modeling of biological soft tissues and personalization of constitutive parameters using non-invasive, non-destructive bedside testing methods. It is imperative to understand the scope and appropriate applications for reported material properties. Thus, the goal of this paper was to compile studies from which soft tissue material properties were obtained and categorize them by source of tissue samples, methods used to quantify deformation, and material models used to describe tissues. The collected studies displayed wide ranges of material properties, and factors that affected the properties included whether tissue samples were in vivo or ex vivo, from humans or animals, the body region tested, body position during in vivo studies, deformation measurements, and material models used to describe tissues. Because of the factors that affected reported material properties, it is clear that much progress has been made in understanding soft tissue responses to loading, yet there is a need to broaden the scope of reported soft tissue material properties and better match reported properties to appropriate human body models.

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.

Determining frictional properties of pants and cushion cover materials using human soft tissue and a rigid sled and how they affect seated shear forces

Shear forces on the buttocks while seated are directly linked to friction, yet the frictional properties at the seat interface are unknown. Shear forces are one of the factors related to increase risk of pressure injury formation. The goals of this study included determining coefficients of friction between three cushion covers and two clothing fabrics using a mechanical system as well as human participants and to evaluate the impact of the cushion covers on shear loading on the buttocks while seated. A chair with separate seat pan tilt and back recline movements was built and instrumented with reflective markers and a load cell. A motion capture system and load cell were used to determine the angles of seat pan tilt at which the sled and participants started sliding, as well as shear forces at three recline angles for three cushion covers (vinyl, one-layer nylon, and two-layer nylon). Results showed the vinyl and two-layer nylon cushion covers respectively had the largest and smallest coefficients of friction for both pants materials. The coefficients of friction calculated with the human participants and rigid sled were within 10% of each other, demonstrating similar results. Further, increasing back recline increased shear load on the buttocks, while the two-layer nylon cover reduced shear forces seen on the buttocks. This work furthers the understanding of shear loading on the buttocks, will aid in the protocols for reducing pressure injuries, and suggests that coefficients of friction found using rigid bodies may be applied to deformable bodies.

An investigation into the hammer toe effects on the lower extremity mechanics and plantar fascia tension: A case for a vicious cycle and progressive damage

Hammer toes are one of the common deformities of the forefoot that can lead to compensatory changes during walking in individuals with this condition. Predicting the adverse effects of tissue damage on the performance of other limbs is very important in the prevention of progressive damage. Finite element (FE) and musculoskeletal modeling can be helpful by allowing such effects to be studied in a way where the internal stresses in the tissue could be investigated. Hence, this study aims to investigate the effects of the hammer toe deformity on the lower extremity, especially on the plantar fascia functions. To compare the joint reactions of the hammer toe foot (HTF) and healthy foot (HF), two musculoskeletal models (MSM) of the feet of a healthy individual and that of a participant with hammer toe foot were developed based on gait analysis. A previously validated 3D finite element model which was constructed using Magnetic Resonance Imaging (MRI) of the diabetic participant with the hammer toe deformity was processed at five different events during the stance phase of gait. It was found that the hammer toe deformity makes dorsiflexion of the toes and the windlass mechanism less effective during walking. Specifically, the FE analysis results showed that plantar fascia (PF) in HTF compared to HF played a less dominant role in load bearing with both medial and lateral parts of PF loaded. Also, the results indicated that the stored elastic energy in PF was less in HTF than the HF, which can indicate a higher metabolic cost during walking. Internal stress distribution shows that the majority of ground reaction forces are transmitted through the lateral metatarsals in hammer toe foot, and the probability of fifth metatarsal fracture and also progressive deformity was subsequently increased. The MSM results showed that the joint reaction forces and moments in the hammer toe foot have deviated from normal, where the metatarsophalangeal joint reactions in the hammer toe were less than the values in the healthy foot. This can indicate a vicious cycle of foot deformity, leading to changes in body weight force transmission line, and deviation of joint reactions and plantar fascia function from normal. These in turn lead to increased internal stress concentration, which in turn lead to further foot deformities. This vicious cycle cause progressive damage and can lead to an increase in the risk of ulceration in the diabetic foot.

Dynamic finite element analyses to compare the influences of customised total talar replacement and total ankle arthroplasty on foot biomechanics during gait

Objective, Total talar replacement (TTR) using a customised talus prosthesis is an emerging surgical alternative to conventional total ankle arthroplasty (TAA) for treating ankle problems. Upon satisfying clinical reports in the literature, this study explored the advantages of TTR in restoring foot biomechanics during walking compared with TAA through computational simulations.Methods, A dynamic finite element foot model was built from the MRIs of a healthy participant and modified into two implanted counterparts (TTR and TAA) by incorporating the corresponding prosthetic components into the ankle joint. Twenty bony parts, thirty-nine ligament/tendon units, nine muscle contractors, and bulk soft tissue were included in the intact foot model. The TTR prosthesis was reconstructed from the mirror image data of the participant's contralateral talus and the TAA prosthesis was modelled by reproducing the Scandinavian ankle replacement procedure in the model assembly. The model was meshed with explicit deformable elements and validated against existing experimental studies that have assessed specific walking scenarios. Simulations were performed using the boundary conditions (time-variant matrix of muscle forces, segment orientation, and ground reaction forces) derived from motion capture analyses and musculoskeletal modelling of the participant's walking gait. Outcome variables, including foot kinematics, joint loading, and plantar pressure were reported and compared among the three model conditions.

RESULTS

Linear regression indicated a better agreement between the TTR model and intact foot model in plots of joint motions and foot segment movements during walking (R² ​= ​0.721-0.993) than between the TAA and intact foot (R² ​= ​0.623-0.990). TAA reduced talocrural excursion by 21.36%-31.92% and increased (MTP) dorsiflexion by 3.03%. Compared with the intact foot, TTR and TAA increased the midtarsal joint contact force by 17.92% and 10.73% respectively. The proximal-to-distal force transmission within the midfoot was shifted to the lateral column in TTR (94.52% or 210.54 ​N higher) while concentrated on the medial column in TAA (41.58% or 27.55 ​N higher). The TTR produced a plantar pressure map similar to that of the intact foot. TAA caused the plantar pressure centre to drift medially and increased the peak forefoot pressure by 7.36% in the late stance.

CONCLUSION

The TTR better reproduced the foot joint motions, segment movements, and plantar pressure map of an intact foot during walking. TAA reduced ankle mobility while increasing movement of the adjacent joints and forefoot plantar pressure. Both implant methods changed force transmission within the midfoot during gait progression.The translational potential of this article Our work is one of the few to report foot segment movements and the internal loading status of implanted ankles during a dynamic locomotion task. These outcomes partially support the conjecture that TTR is a prospective surgical alternative for pathological ankles from a biomechanical perspective. This study paves the way for further clinical investigations and systematic statistics to confirm the effects of TTR on functional joint recovery.

Oblique impact responses of Hybrid III and a new headform with more biofidelic coefficient of friction and moments of inertia

New oblique impact methods for evaluating head injury mitigation effects of helmets are emerging, which mandate measuring both translational and rotational kinematics of the headform. These methods need headforms with biofidelic mass, moments of inertia (MoIs), and coefficient of friction (CoF). To fulfill this need, working group 11 of the European standardization head protection committee (CEN/TC158) has been working on the development of a new headform with realistic MoIs and CoF, based on recent biomechanics research on the human head. In this study, we used a version of this headform (Cellbond) to test a motorcycle helmet under the oblique impact at 8 m/s at five different locations. We also used the Hybrid III headform, which is commonly used in the helmet oblique impact. We tested whether there is a difference between the predictions of the headforms in terms of injury metrics based on head kinematics, including peak translational and rotational acceleration, peak rotational velocity, and BrIC (brain injury criterion). We also used the Imperial College finite element model of the human head to predict the strain and strain rate across the brain and tested whether there is a difference between the headforms in terms of the predicted strain and strain rate. We found that the Cellbond headform produced similar or higher peak translational accelerations depending on the impact location (−3.2% in the front-side impact to 24.3% in the rear impact). The Cellbond headform, however, produced significantly lower peak rotational acceleration (−41.8% in a rear impact to −62.7% in a side impact), peak rotational velocity (−29.5% in a side impact to −47.6% in a rear impact), and BrIC (−29% in a rear-side impact to −45.3% in a rear impact). The 90th percentile values of the maximum brain strain and strain rate were also significantly lower using this headform. Our results suggest that MoIs and CoF have significant effects on headform rotational kinematics, and consequently brain deformation, during the helmeted oblique impact. Future helmet standards and rating methods should use headforms with realistic MoIs and CoF (e.g., the Cellbond headform) to ensure more accurate representation of the head in laboratory impact tests.

The contribution of the ligaments in progressive collapsing foot deformity: A comprehensive computational study

The contribution of each of the ligaments in preventing the arch loss, hindfoot valgus, and forefoot abduction seen in progressive collapsing foot deformity (PCFD) has not been well characterized. An improved understanding of the individual ligament contributions to the deformity would aid in selecting among available treatments, optimizing current surgical techniques, and developing new ones. In this study, we evaluated the contribution of each ligament to the maintenance of foot alignment using a finite element model of the foot reconstructed from computed tomography scan images. The collapsed foot was modeled by simulating the failure of all the ligaments involved in PCFD. The ligaments were removed one at a time to determine the impact of each ligament on foot alignment, and then restored one at a time to simulate isolated reconstruction. Our findings show that the failure of any one ligament did not immediately lead to deformity, but that combined failure of only a few (the plantar fascia, long plantar, short plantar, deltoid, and spring ligaments) could lead to significant deformity. The plantar fascia, deltoid, and spring ligaments were primarily responsible for the prevention of arch collapse, hindfoot valgus, and forefoot abduction, respectively. Moreover, to produce deformity, a considerable amount of attenuation in the spring, tibiocalcaneal, interosseous talocalcaneal, plantar naviculocuneiform, and first plantar tarsometatarsal ligaments, but only a small amount in the plantar fascia, long plantar, and short plantar ligaments was needed. The results of this study suggest that the ability of a ligament to prevent deformity may not correlate with its attenuation in a collapsed foot.

[Three-dimensional finite element analysis of Swanson prosthesis-arthroplasty of the first metatarsophalangeal joint combined with osteotomy and bone grafting of the first metatarsal bone for hallux valgus]

OBJECTIVE

To analyze the biomechanical changes of hallux valus after Swanson prosthesis-arthroplasty of the 1st metatarsophalangeal joint combined with osteotomy and bone grafting of the 1st metatarsal bone by three-dimensional finite element analysis, so as to provide data basis for studying the changes of foot morphology and physiological function after hallux valus correction surgery.

METHODS

A 65-year-old female patient with severe hallux valus admitted in January 2013 was selected as the research object. The CT data of the right foot was obtained, and the three-dimensional finite element models before and after Swanson prosthesis-arthroplasty of the 1st metatarsophalangeal joint combined with osteotomy and bone grafting of the 1st metatarsal bone were established by Mimics10.01, Geomagic Studio, and ANSYS12.0 software. ANSYS 12.0 software was used for nonlinear static stress analysis, and the hallux valgus angle (HVA), the intermetatarsal angle (IMA), and the von Mises stress distributions of the forefoot plantar surface and the 1st to 5th metatarsal bones were observed before and after operation.

RESULTS

The HVA and IMA were 56.3° and 16.3° before operation and 9.2° and 9.8° after operation, respectively. Before operation, the stress on the forefoot was the largest in the 4th metatarsal head zone and the smallest in the 1st metatarsal head zone; the stress on the medial side of the forefoot was significantly smaller than that on the lateral side, and the center of forefoot pressure was located on the lateral side. After operation, the stress on the forefoot was the largest in the 1st metatarsal head zone and the smallest in the 5th metatarsal head zone; the stress on the lateral side of the forefoot was significantly smaller than that on the medial side, and the center of forefoot pressure was located on the medial side. Before operation, the stress of the 5th metatarsal bone was the largest, and the 1st metatarsal bone was the smallest. After operation, the stress of the 1st metatarsal bone was the largest, and the 4th metatarsal bone was the smallest.

CONCLUSION

Swanson prosthesis-arthroplasty of the 1st metatarsophalangeal joint combined with osteotomy and bone grafting of the 1st metatarsal bone can effectively correct hallux valgus and make HVA, IMA, and plantar pressure distribution close to normal. However, postoperative stresses of the 1st to 5th metatarsal bones elevate, which may lead to associated complications.

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.

An investigation of the ankle contact forces in a foot with hammer toe deformity. A comparison of patient-specific approaches using finite element modeling and musculoskeletal simulation

The internal forces and stresses in the tissue are important as they are linked to the risk of mechanical trauma and injuries. Despite their value, the internal stresses and forces cannot be directly measured in-vivo. A previously validated 3D finite element model (FEM) was constructed using Magnetic Resonance Imaging (MRI) of a person with diabetes and hammer toe deformity. The foot model simulated at five different instances during the stance phase of gait. The internal stress distribution on the talus that was obtained using the FEM simulation, was used to calculate the joint reaction force at the ankle joint. In addition, the musculoskeletal model (MSM) of the participant with hammer toe foot was developed based on the gait analysis and was used to determine the muscle forces and joint reactions. The result showed that the vertical reaction forces obtained from the FEM and MSM follow a similar trend through the stance phase of gait cycle and are significantly correlated ( R=0.99). The joint reaction forces obtained through the two methods do not differ for the first 25% of the gait cycle, while the maximum difference was ∼0.7 Body weight that was observed at 50% of the stance phase. Clinical Relevance: Finite element modeling and musculoskeletal simulation can shed light on the internal forces at the ankle in pathological conditions such as hammer toe. The similarities and differences observed in the joint reaction forces calculated from the two methods can have implications in assessing the effect of clinical interventions.

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 Handheld Gyroscopic Device for Haptics and Hand Rehabilitation

We propose a novel, gyroscopic device for haptics and hand rehabilitation, named Gymball. It consists of a fully actuated rotor-gimbal assembly encased in an easy-to-grip appealing design. When held, the device generates a gyroscopic torque which causes the user's hand to move about the wrist. Interviews with occupational therapists, simulations, and proof-of-concept models helped determine the design specifications of Gymball. Compared to the existing gyroscopic devices, Gymball has the following advantages. (i) A smaller form-factor with better user appeal while achieving 0.5 Nm torque. (ii) A wire entanglement-free design allowing complete rotations of the rotor-gimbal assembly. (iii) Negligible rotary imbalances owing to a symmetrical design, resulting in haptic signals with minimal vibratory noise. In this paper, we detail the design and analysis of the device. A feasibility study was conducted to validate prospect of using the device for haptic feedback or therapy. Specifically, the study focused on (i) whether the gyroscopic torque generated by the device can passively move the user's hand about the wrist and (ii) whether the produced hand motion can be controlled. The results show that Gymball can successfully generate about 7^° of hand oscillations. The amplitude and frequency of the hand oscillations can be controlled using the speed of rotor and gimbal.

Biomechanical comparison among five mid/hindfoot arthrodeses procedures in treating flatfoot using a musculoskeletal multibody driven finite element model

BACKGROUND AND OBJECTIVE

Mid/hindfoot arthrodesis could modify the misalignment of adult-acquired flatfoot and attenuate pain. However, the long-term biomechanical effects of these surgical procedures remain unclear, and the quantitative evidence is scarce. Therefore, we aimed to investigate and quantify the influences of five mid/hindfoot arthrodeses on the internal foot biomechanics during walking stance.

METHODS

A young participant with flexible flatfoot was recruited for this study. We reconstructed a subject-specific musculoskeletal multibody driven-finite element (FE) foot model based on the foot magnetic resonance imaging. The severe flatfoot model was developed from the flexible flatfoot through the attenuation of ligaments and the unloading of the posterior tibial muscle. The five mid/hindfoot arthrodeses simulations (subtalar, talonavicular, calcaneocuboid, double, and triple arthrodeses) and a control condition (no arthrodesis) were performed simultaneously in the detailed foot multibody dynamics model and FE model. Muscle forces calculated by a detailed multi-segment foot model and ground reaction force were used to drive the foot FE model. The internal foot loadings were compared among control and these arthrodeses conditions at the first and second vertical ground reaction force (VGRF) peak and VGRF valley instants.

RESULTS

The results indicated that the navicular heights in double and triple arthrodeses were higher than other surgical procedures, while the subtalar arthrodesis had the smallest values. Five mid/hindfoot arthrodeses reduced the peak plantar fascia stress compared to control. However, double and triple arthrodeses increased the peak medial cuneo-navicular joint contact pressures and peak foot pressures as well as the metatarsal bones stresses.

CONCLUSION

Although mid/hindfoot arthrodesis generally reduced the collapse of medial longitudinal arch and plantar fascia loading during the stance phase, the increased loading in the adjacent unfused joint and metatarsal bones for double and triple arthrodeses should be noted. These findings could account for some symptoms experienced by flatfoot patients after surgery, which may facilitate the optimization of surgical protocols.

Evaluation of facial tissue stresses under medical devices post application of a cyanoacrylate liquid skin protectant: An integrated experimental-computational study

Medical device‐related pressure ulcers (PUs) (injuries) are a subclass of PUs, associated with pressure and/or shear applied by a medical device onto the skin. Clinical application of a cyanoacrylate liquid skin protectant (CLSP) under the contours of skin‐contacting medical devices to shield an intact skin from the sustained mechanical loads that are applied by medical devices is a preventative option, but no computer modelling work has been reported to assess the biomechanical efficacy of such interventions. Here, we investigated the biomechanical protective effect of a polymerised cyanoacrylate coating using three‐dimensional, anatomically realistic finite element models of the ear with oxygen cannula and the mouth with endotracheal attachment device, informed by experimental studies. We have compared tissue stress exposures under the devices at these facial sites between conditions where the cyanoacrylate skin protectant has been applied or where the device was contacting the skin directly, without the shielding of the cyanoacrylate coating. The CLSP considerably reduced the skin stress concentration levels and overall tissue stress exposures under the aforementioned medical devices. This demonstrates strong biomechanical effectiveness of the studied cyanoacrylate‐based skin protectant in prevention of facial medical device‐related injuries at small, curved and thereby difficult to protect facial sites.

A numerical analysis of skin-PPE interaction to prevent facial tissue injury

The use of close-fitting PPE is essential to prevent exposure to dispersed airborne matter, including the COVID-19 virus. The current pandemic has increased pressure on healthcare systems around the world, leading to medical professionals using high-grade PPE for prolonged durations, resulting in device-induced skin injuries. This study focuses on computationally improving the interaction between skin and PPE to reduce the likelihood of discomfort and tissue damage. A finite element model is developed to simulate the movement of PPE against the face during day-to-day tasks. Due to limited available data on skin characteristics and how these vary interpersonally between sexes, races and ages, the main objective of this study was to establish the effects and trends that mask modifications have on the resulting subsurface strain energy density distribution in the skin. These modifications include the material, geometric and interfacial properties. Overall, the results show that skin injury can be reduced by using softer mask materials, whilst friction against the skin should be minimised, e.g. through use of micro-textures, humidity control and topical creams. Furthermore, the contact area between the mask and skin should be maximised, whilst the use of soft materials with incompressible behaviour (e.g. many elastomers) should be avoided.

The effect of frictional coefficients and sock material on plantar surface shear stress measurement

At present, there are no viable systems that can acquire in-shoe measurement of distributed shear forces. Foot-shoe interactions are such that skin shear is a notoriously difficult quantity to measure under the best of conditions. This is further complicated by the presence of forces normal to the skin surface that are large compared to the shear forces, which often results in crosstalk between pressure and shear signals. The present study used multibody dynamic simulations to investigate the combined effects of (i) coefficient of friction (COF) at skin-sock and sock-sensor interfaces, as well as (ii) sock stiffness on the accuracy of measured shear against the skin. These factors were systematically altered within a wide range (COF: 0.04, 0.34, 0.54, and 0.9; sock stiffness: 100, 250, 500, 1000, 1500 and 2000 N/m) to simulate a total of 96 scenarios. The correlation between the shear at the skin and at the sensor was used to compare each set of conditions. The results indicated that a high COF at the sock-sensor interface and a low sock stiffness would individually result in a significantly higher accuracy of shear measurements (p < 0.001). A low COF at the skin-sock interface was observed to reduce the occurred shear against the skin up to a factor of five, with very minimal effect on the accuracy of shear measurements (p = 0.98). These findings allow researchers to understand the potential effects of (i) sock stiffness, and (ii) coefficients of friction, on skin shear, and potentially correct for the effects of interface materials when trying to determine shear at the skin-sock interface.

Changes in the three-dimensional microscale topography of human skin with aging impact its mechanical and tribological behavior

Human skin enables interaction with diverse materials every day and at all times. The ability to grasp objects, feel textures, and perceive the environment depends on the mechanical behavior, complex structure, and microscale topography of human skin. At the same time, abrasive interactions, such as sometimes occur with prostheses or textiles, can damage the skin and impair its function. Previous theoretical and computational efforts have shown that skin's surface topography or microrelief is crucial for its tribological behavior. However, current understanding is limited to adult surface profiles and simplified two-dimensional simulations. Yet, the skin has a rich set of features in three dimensions, and the geometry of skin is known to change with aging. Here we create a numerical model of a dynamic indentation test to elucidate the effect of changes in microscale topography with aging on the skin's response under indentation and sliding contact with a spherical indenter. We create three different microrelief geometries representative of different ages based on experimental reports from the literature. We perform the indentation and sliding steps, and calculate the normal and tangential forces on the indenter as it moves in three distinct directions based on the characteristic skin lines. The model also evaluates the effect of varying the material parameters. Our results show that the microscale topography of the skin in three dimensions, together with the mechanical behavior of the skin layers, lead to distinctive trends on the stress and strain distribution. The major finding is the increasing role of anisotropy which emerges from the geometric changes seen with aging.

Simple model of arch support: Relevance to Charcot Neuroarthropathy

BACKGROUND

Charcot neuropathy is a common complication resulting from poorly controlled diabetes and peripheral neuropathy leading to the collapse, and ultimately the breakdown, of the midfoot. Mechanically, it is likely that a compromised arch support in this, or any other patient group that experiences foot flattening, would be associated with slippage at the distal and proximal interface regions of the plantar surface of the foot and the adjacent support surface. This slippage, although difficult to quantify with standard motion capture systems used in a gait laboratory, could potentially be assessed with systems for monitoring interface shear stresses. However, before investing in such systems, a correlation between arch flattening and interface shear stresses needs to be verified.

METHODS

For this purpose, a sagittal plane model of a foot was developed using a multi-body dynamics package (MSC Adams). This model mimicked a subject swaying back and forth, and was constructed to show the dependence of interface stresses on altered arch support.

FINDINGS

The model's predictions matched typical FootSTEPS data: lengthening of the arch of 1-2 mm, sway oscillations of 0.22-0.33 s and frictional force differences (calcaneus relative to forefoot) of 60 N. Of clinical relevance, when the stiffness of the plantar spring (representing aponeurosis and intrinsic muscles) was reduced by 10%, the frictional force difference increased by about 6.5%.

INTERPRETATION

The clinical implications of this study are that, while arch lengthening of less than 2 mm might be difficult to measure reliably in a gait lab, using shear sensors under the forefoot and hindfoot should allow arch support to be assessed in a repeatable manner.

Associations between changes in loading pattern, deformity, and internal stresses at the foot with hammer toe during walking; a finite element approach

Over the past decade, Finite Element (FE) modelling has been used as a method to understand the internal stresses within the diabetic foot. Foot deformities such as hammer toe have been associated with increased risk of foot ulcers in diabetic patients. Hence the aim of this study is to investigate the influence of hammer toe deformity on internal stresses during walking. A 3D finite element model of the human foot was constructed based on capturing Magnetic Resonance Imaging (MRI) of a diabetic neuropathic volunteer exhibiting hammer toe. 3D gait measurements and a multi-body musculoskeletal model for the same participant were used to define muscle forces. FE simulations were run at five different instances during the stance phase of gait. Peak plantar pressure and pressure distribution results calculated from the model showed a good agreement with the experimental measurement having less than 11% errors. Maximum von Mises internal stresses in the forefoot hard tissue were observed at the 3^rd and 5^th metatarsals and 4^th proximal phalanx. Moreover, presence of hammer toe deformity was found to shift the location of maximum internal stresses on the soft tissue to the forefoot by changing the location of centre of pressure with internal stress 1.64 times greater than plantar pressure. Hammer toe deformity also showed to reduce the involvement of the first phalanx in internal/external load-bearing during walking. The findings of this study support the association between changes in loading pattern, deformity, and internal stresses in the soft tissue that lead to foot ulceration.

A review of foot finite element modelling for pressure ulcer prevention in bedrest: Current perspectives and future recommendations

Pressure ulcers (PUs) are a major public health challenge, having a significant impact on healthcare service and patient quality of life. Computational biomechanical modelling has enhanced PU research by facilitating the investigation of pressure responses in subcutaneous tissue and skeletal muscle. Extensive work has been undertaken on PUs on patients in the seated posture, but research into heel ulcers has been relatively neglected. The aim of this review was to address the key challenges that exist in developing an effective FE foot model for PU prevention and the confusion surrounding the wide range of outputs reported. Nine FE foot studies investigating heel ulcers in bedrest were identified and reviewed. Six studies modelled the posterior part of the heel, two included the calf and foot, and one modelled the whole body. Due to the complexity of the foot anatomy, all studies involved simplification or assumptions regarding parts of the foot structure, boundary conditions and material parameters. Simulations aimed to understand better the stresses and strains exhibited in the heel soft tissues of the healthy foot. The biomechanical properties of soft tissue derived from experimental measurements are critical for developing a realistic model and consequently guiding clinical decisions. Yet, little to no validation was reported in each of the studies. If FE models are to address future research questions and clinical applications, then sound verification and validation of these models is required to ensure accurate conclusions and prediction of patient outcomes. Recommendations and considerations for future FE studies are therefore proposed.

Effects of Custom-Made Insole Materials on Frictional Stress and Contact Pressure in Diabetic Foot with Neuropathy: Results from a Finite Element Analysis

Offloading plantar pressure in a diabetic foot with neuropathy is challenging in conventional clinical practice. Custom-made insole (CMI) materials play an important role in plantar pressure reduction, but the assessment is costly and time-consuming. Finite element analysis (FEA) can provide an efficient evaluation of different insoles on the plantar pressure distribution. This study investigated the effect of CMI materials and their combinations on plantar pressure reduction for the diabetic foot with neuropathy using FEA. The study was conducted by constructing a three-dimensional foot model along with CMI to study the peak contact pressure between the foot and CMI. The softer material (E = 5 MPa) resulted in a better reduction of peak contact pressure compared with the stiffer material (E = 11 MPa). The plantar pressure was well redistributed with softer material compared with the stiffer material and its combination. In addition, the single softer material resulted in reduced frictional stress under the first metatarsal head compared with the stiffer material and the combination of materials. The softer material and its combination have a beneficial effect on plantar pressure reduction and redistribution for a diabetic foot with neuropathy. This study provided an effective approach for CMI material selection using FEA.

Real-Time Analysis of the Dynamic Foot Function: A Machine Learning and Finite Element Approach

Finite element analysis (FEA) has been widely used to study foot biomechanics and pathological functions or effects of therapeutic solutions. However, development and analysis of such foot modeling is complex and time-consuming. The purpose of this study was therefore to propose a method coupling a FE foot model with a model order reduction (MOR) technique to provide real-time analysis of the dynamic foot function. A generic and parametric FE foot model was developed and dynamically validated during stance phase of gait. Based on a design of experiment of 30 FE simulations including four parameters related to foot function, the MOR method was employed to create a prediction model of the center of pressure (COP) path that was validated with four more random simulations. The four predicted COP paths were obtained with a 3% root-mean-square-error (RMSE) in less than 1 s. The time-dependent analysis demonstrated that the subtalar joint position and the midtarsal joint laxity are the most influential factors on the foot functions. These results provide additionally insight into the use of MOR technique to significantly improve speed and power of the FE analysis of the foot function and may support the development of real-time decision support tools based on this method.

Wearable Artificial Fingers With Skin Vibration and Multi-Axis Force Sensors

Tactile sensations are based on stimulation elicited on the skin through mechanical interactions between the skin and an object. Hence, it is important to consider skin properties as well as objects. In this article, we aim to develop wearable artificial fingers for quantitative evaluations reflecting individual differences in human fingers. In a previous study, a wearable skin vibration sensor was attached to artificial fingers and it was demonstrated that the skin vibrations differed according to the dimension of surface ridge and the artificial finger is useful for roughness evaluation. This article improved the artificial finger to measure the contact force and friction in addition to the skin vibration. A small three-axis force sensor was embedded into the base of the finger, and normal and friction forces were estimated via a multi-regression method. Furthermore, artificial fingers with different hardness were prepared and six different textures were used to investigate tactile evaluation. Experimental results showed that the artificial fingers could measure normal and friction forces along with the skin vibration and were useful to evaluate textures. Resulting distributions of the vibration intensity and friction coefficient were different for the soft and hard artificial fingers, indicating the complex influence of skin properties on tactile sensations.

Efficacy of Computer-Aided Design and Manufacturing Versus Computer-Aided Design and Finite Element Modeling Technologies in Brace Management of Idiopathic Scoliosis: A Narrative Review

The efficiency and design quality of scoliosis braces produced by the conventional casting method depends highly on the orthotist's experience. Recently, advanced engineering techniques have been used with the aim of improving the quality of brace design and associated clinical outcomes. Numerically controlled machine tools have provided enormous opportunities for reducing the manufacturing time and saving material. However, the effectiveness of computer-aided brace manufacturing for scoliosis curve improvement is controversial. This narrative review is aimed at comparing the efficacy of braces made by the conventional method with those made by two computer-aided methods: computer-aided design and manufacturing (CAD-CAM), and computer-aided design and finite element modeling (CAD-FEM). The comparison was performed on scoliosis parameters in coronal, sagittal, and transverse planes. Scientific databases were searched, and 11 studies were selected for this review. Because of the diversity of study designs, it was not possible to decisively conclude which brace-manufacturing method is most effective. Similar effectiveness in curve correction was found in the coronal plane for braces made by using advanced manufacturing and conventional methods. In the sagittal plane, modern braces seem to be more effective than traditional braces, but there is an ongoing debate among clinicians, about which CAD-CAM and CAD-FEM brace provides a better treatment outcome. The relative effectiveness of modern and conventional methods in correcting deformities in the transverse plane is also a controversial subject. Overall, advanced engineering design and production methods can be proposed as time- and cost-efficient approaches for scoliosis management. However, there is insufficient evidence yet to conclude that CAD-CAM, and CAD-FEM methods provide significantly better clinical outcomes than those of conventional methods in the treatment of scoliosis curve. Moreover, for some factors, such as molding and the patient's posture during the data acquisition, in brace curve-correction plan, the orthotist's experience and scoliosis curve flexibility should be explored to confidently compare the outcomes of conventional, CAD-CAM, and CAD-FEM methods.

Finite element-based method for determining an optimal offloading design for treating and preventing heel ulcers

Diabetic heel ulceration, a serious, destructive, and costly complication of diabetes, is often treated by custom-made offloading footwear. One common offloading device is a custom-made insole designed with a hole under the damaged site that is intended to reduce local mechanical loads on the ulcer. However, current devices do not take into account the increasing loads at the wound peripheries, and quantitative assessments and scientific guidelines for the optimal design of the offloading hole are lacking. Here, we develop a novel method to determine the volumetric exposure to mechanical loading of a human heel, at two volume of interests (VOIs) during walking in 150 different finite-element footwear configurations. We defined the two VOIs as (1) the area of the heel soft tissues typically at high risk of ulceration, and (2) the soft tissues surrounding the high risk area. For all model variants, three hole-geometry parameters were defined: (1) radius, (2) radius of curvature (ROC) and (3) depth. We found two combinations of the offloading parameters which minimize heel loads in both VOIs. The first is with a large offloading radius, large ROC and large depth, whereas the second is with a large offloading radius, large depth but relatively small ROC. Our novel practical scientific analysis method, that takes into account the ulcer site as well as the peripheral area, has the potential to optimize development of offloading solutions by streamlining the examination of their biomechanical efficiency, and thus may revolutionize prevention and treatment of diabetic ulcers at any foot location.

Customized Three-Dimensional-Printed Orthopedic Close Contact Casts for the Treatment of Stable Ankle Fractures: Finite Element Analysis and a Pilot Study

Ankle fracture is one of the most common traumatic fractures among the elderly population. The majority of ankle fractures are stable types with the typically conservative strategy of close contact casting treatment. The continuous use of unventilated standard cast immobilization severely affects patient's satisfaction and compliance and markedly increases the rates of various complications. Three-dimensional (3D) printing for casts has advantages of lightweight, ventilated, proper-fit, and esthetic improvements. In this work, this novel 3D-printed cast has been applied to individuals with stable ankle fractures, and its effectiveness can be successfully validated with finite element analysis and a pilot study. A 30% reduction of the volume was chosen as the optimal result in topology optimization. Both 3D-printed casts and conventional casts showed significant ankle function improvement after immobilization for 6 weeks (p = 0.000). The 3D-printed casts were superior to the traditional casts in Olerud-Molander Ankle Scores (OMAS), with the mean difference of 8.3 ± 8.57 OMAS points (95% CI -10.8 to 27.5; p = 0.354) for 6 weeks, implying that the 3D-printed casts possibly maintain the equal clinical efficacy as the traditional casts. The statistically significant difference between groups from the 3D-printed cast and the traditional one observed in C-QUEST 2.0 was 11.3 ± 1.5 points (95% CI 8.0-14.6; p = 0.000), indicating that the 3D-printed cast possesses outperforming satisfaction and compliance and has great potential in practical applications. There were no severe complications in the 3D-printed casts, but more moderate complications were observed in the traditional casts.

Bio-Numerical Analysis of the Human Ankle-Foot Model Corresponding to Neutral Standing Condition

Background:

The foot is the most complex body’s structure; it is highly susceptible to disorders because of its loading pattern. The complexity of the foot structure geometry implies the use of reverse engineering tools to obtain a model that can accurately mimic the biomechanical behavior of the foot.

Objective:

The objective of this study is to establish a state-of-the-art ankle-foot finite element (FE) model with anatomically realistic geometry and structure in order to get the model that will suit all cases for future studies on stress injuries and foot insole designs under different loading conditions.

Material and Methods:

In this analytical study, tomography images were imported in DICOM format, after that, the object was exported in the form of three-dimensional structures in STL file format to define and assemble the structures. After that, the computer simulation on numerical model was done. One-way Analysis of variance (ANOVA) test was performed, and a threshold (p&lt;0.05) was used to indicate the significance of results.

Results:

The results showed no significant differences (P>0.05) between the values of the plantar pressure corresponding to neutral standing condition with other foot models in literature. The stresses transferred to the bone structure show that the relatively higher stress was located in the fifth, fourth and third tarsometatarsal, where the maximum von Mises stress in the bone structure was 2155.4 kPa.

Conclusion:

The state-of-the-art ankle-foot FE model with anatomically realistic geometry and structure will be very helpful for future studies on stress injuries and foot insole designs under different loading conditions.

Expandable and implantable bioelectronic complex for analyzing and regulating real-time activity of the urinary bladder

Underactive bladder or detrusor underactivity (DUA), that is, not being able to micturate, has received less attention with little research and remains unknown or limited on pathological causes and treatments as opposed to overactive bladder, although the syndrome may pose a risk of urinary infections or life-threatening kidney damage. Here, we present an integrated expandable electronic and optoelectronic complex that behaves as a single body with the elastic, time-dynamic urinary bladder with substantial volume changes up to ~300%. The system configuration of the electronics validated by the theoretical model allows conformal, seamless integration onto the urinary bladder without a glue or suture, enabling precise monitoring with various electrical components for real-time status and efficient optogenetic manipulation for urination at the desired time. In vivo experiments using diabetic DUA models demonstrate the possibility for practical uses of high-fidelity electronics in clinical trials associated with the bladder and other elastic organs.

Contact Pressures Between the Rearfoot and a Novel Offloading Insole: Results From a Finite Element Analysis Study

Pressure offloading is critical to diabetic foot ulcer healing and prevention. A novel product has been proposed to achieve this offloading with an insole that can be easily modified for each user. This insole consists of pressurized bubbles that can be selectively perforated and depressurized to redistribute weight to the nonulcer region of the foot. However, the effect of the insole design parameters, for example, bubble height and stiffness, on offloading effectiveness is unknown. To this end, a 3-dimensional finite element model was developed to simulate contact between the rearfoot and insole. The geometry of the calcaneus bone and soft tissue was based on the medical images of an average male patient, and material properties and loading conditions based on the values reported in the literature were used. The model predicts that increasing bubble height and stiffness leads to a more effectively offloaded region. However, the model also predicts that increasing stiffness leads to increasing contact pressures on the surrounding soft tissue. Thus, a combination of insole design parameters was determined, which completely offloads the desired region, while simultaneously reducing the contact pressure on the surrounding soft tissue. This design is expected to aid in diabetic foot ulcer healing and prevention.

Finite Element Analysis of Generalized Ligament Laxity on the Deterioration of Hallux Valgus Deformity (Bunion)

Hallux valgus is a common foot problem affecting nearly one in every four adults. Generalized ligament laxity was proposed as the intrinsic cause or risk factor toward the development of the deformity which was difficult to be investigated by cohort clinical trials. Herein, we aimed to evaluate the isolated influence of generalized ligament laxity on the deterioration using computer simulation (finite element analysis). We reconstructed a computational foot model from a mild hallux valgus participant and conducted a gait analysis to drive the simulation of walking. Through parametric analysis, the stiffness of the ligaments was impoverished at different degrees to resemble different levels of generalized ligament laxity. Our simulation study reported that generalized ligament laxity deteriorated hallux valgus by impairing the load-bearing capacity of the first metatarsal, inducing higher deforming force, moment and malalignment at the first metatarsophalangeal joint. Besides, the deforming moment formed a deteriorating vicious cycle between hallux valgus and forefoot abduction and may result in secondary foot problems, such as flatfoot. However, the metatarsocuneiform joint did not show a worsening trend possibly due to the overriding forefoot abduction. Controlling the deforming load shall be prioritized over the correction of angles to mitigate deterioration or recurrence after surgery.

Determination of Three-Dimensional Corrective Force in Adolescent Idiopathic Scoliosis and Biomechanical Finite Element Analysis

AIMS

In this study we have considered the three dimensional corrective forces for correction of scoliosis by using a patient specific finite element model.

MATERIALS AND METHODS

An objective function of corrective forces in three-dimensional space was defined. Computed tomography images were used to reconstruct three dimensional model of scoliotic trunk. Computer aided engineering software Abaqus was used to establish finite element model of deformed spine and its biomechanical characteristics were analyzed. By adjusting magnitude and position of corrective forces, objective function was minimized to achieve best orthopedic effect. The proposed corrective conditions were divided into three groups: (1) thoracic deformity; (2) lumbar deformity; (3) both thoracic and lumbar deformities were considered.

RESULTS

In all three cases, the objective function was reduced by 58, 52, and 63%, respectively. The best correction forces point was located on convex side of maximum displacement of vertebral body.

CONCLUSION

Using minimum objective function method, spinal deformity in three-dimensional space can be sufficiently reduced. This study provides scientific basis for design of a new corrective brace for treatment of scoliosis.

Monitoring Contact Behavior During Assisted Walking With a Lower Limb Exoskeleton

Interaction forces between a robotic cuff and the skin during assisted walking with a lower limb exoskeleton may cause skin injuries over time, such as blisters and ulcers. This study proposes a sensing cuff that can monitor the contact behavior between the exoskeleton and skin during assisted walking, and a functional test and assisted walking experiments were conducted to evaluate the performance of the proposed device. The functional test of the sensing cuff showed good performance of capturing the contact behavior for safety evaluation. The walking experiment involved subjects walking on a treadmill with a lower limb exoskeleton under different conditions (i.e., walking speed and clothing), and the sensing cuff attached to the exoskeleton measured the interaction forces and slip velocity. The magnitude of shear force in the movement direction peaked near the beginning and within 40 - 50% of the gait cycle. The contact safety of the lower limb exoskeleton during assisted walking was then evaluated based on the calculated shear stress. The designed sensing cuff could provide sufficient information regarding contact behavior and contact safety during assisted walking.

Finite element analysis of subtalar joint arthroereisis on adult-acquired flexible flatfoot deformity using customised sinus tarsi implant

Background

Subtalar arthroereisis may cause sinus tarsi pain complications. In this study, we aimed to introduce a customised implant that facilitated treatment effect and less impingement. The biomechanical outcome between the intact and implant conditions was compared using finite element analysis.

Methods

A female patient with flatfoot (age: 36 years, height: 156 ​cm, body mass: 51 ​kg) was recruited as the model patient. The customised implant was designed from the extracted geometry. Boundary and loading conditions were assumed from the data of a normal participant. Four gait instants, including the ground reaction force first peak (25% stance), valley (45%), initial push-off (60%) and second peak (75%) were analyzed.

Results

The navicular height was elevated by 4.2% at 25% stance, whereas the strain of the spring, plantar cuneonavicular and plantar cuboideonavicular ligaments were reduced. The talonavicular joint force decreased and the calcaneocuboid joint increased by half and 67%, respectively, representing a lateralised load pathway. There was a stress concentration at the sulcus tali reaching 15.29 ​MPa

Conclusion

Subtalar arthroereisis using a customised implant may produce some positive treatment effects in terms of navicular height elevation, ligament strain relief and lateralised joint loading pathway. Although the concentrated stress at the sulcus tali did not exceed the threshold of bone breakdown, we could not rule out the potential of vascular disturbance owing to the remarkable elevation of stress. Future study may enlarge the contact area of the bone–implant interface by considering customisation based on the dynamic change of the sinus tarsi during walking gait.

The translational potential of this article

Geometry mismatch of prefabricated implants could be the reason for complications. With the advancement of 3D printing, customising implant becomes possible and may improve treatment outcome. This study implemented a theoretical model approach to explore its potential under a simulation of walking.

Effect of Material Hardness on Friction Between a Bare Finger and Dry and Lubricated Artificial Skin

Understanding the tribological phenomena when fingers slide over soft surfaces such as skin is important for many practical applications. Therefore, this article analyzed the coefficients of friction for a bare finger sliding over artificial skin with different hardness under dry and lubricated surface conditions. This article contrasts with previous research that predominantly analyzed the contact between skin and hard surfaces or probes. Under dry conditions, the coefficient of friction was constant for artificial skins that were harder than the finger pad, irrespective of the normal force of the finger. However, the coefficient of friction decreased with increasing normal force for softer artificial skins. When the surface of the artificial skin model was lubricated with mica, the coefficient of friction exhibited normal-force dependence only for soft artificial skins, similar to the observations under dry conditions. This effect was due to the deformation friction; thus, the coefficient of friction increased as the normal force increased. Conversely, when the model was lubricated with TiO₂, the coefficient of friction depended on the normal force for all hardness levels. These findings provide insights into the friction experienced during skin-skin or skin-soft material contact under dry and lubricated conditions that can easily occur in daily life. Thus, the results of this study can be useful for the development of skin care products or assistive robots involving human-robot contact.