Simple finite element models for use in the design of therapeutic footwear

Therapeutic footwear is frequently prescribed in cases of rheumatoid arthritis and diabetes to relieve or redistribute high plantar pressures in the region of the metatarsal heads. Few guidelines exist as to how these interventions should be designed and what effect such interventions actually have on the plantar pressure distribution. Finite element analysis has the potential to assist in the design process by refining a given intervention or identifying an optimal intervention without having to actually build and test each condition. However, complete and detailed foot models based on medical image segmentation have proven time consuming to build and computationally expensive to solve, hindering their utility in practice. Therefore, the goal of the current work was to determine if a simplified patient-specific model could be used to assist in the design of foot orthoses to reduce the plantar pressure in the metatarsal head region. The approach is illustrated by a case study of a diabetic patient experiencing high pressures and pain over the fifth metatarsal head. The simple foot model was initially calibrated by adjusting the individual loads on the metatarsals to approximate measured peak plantar pressure distributions in the barefoot condition to within 3%. This loading was used in various shod conditions to identify an effective orthosis. Model results for metatarsal pads were considerably higher than measured values but predictions for uniform surfaces were generally within 16% of measured values. The approach enabled virtual prototyping of the orthoses, identifying the most favorable approach to redistribute the patient’s plantar pressures.     

Consequences of forward translation of the point of force application for the mechanics of running

In running humans, the point of force application between the foot and the ground moves forwards during the stance phase. Our aim was to determine the mechanical consequences of this 'point of force translation' (POFT). We modified the planar spring-mass model of locomotion to incorporate POFT, and then compared spring-mass simulations with and without POFT. We found that, if leg stiffness is adjusted appropriately, it is possible to maintain very similar values of peak vertical ground reaction force (GRF), stance time, contact length and vertical centre of mass displacement, whether or not POFT occurs. The leg stiffness required to achieve this increased as the distance of POFT increased. Peak horizontal GRF and mechanical work per step were lower when POFT occurred. The results indicate that the lack of POFT in the traditional spring-mass model should not prevent it from providing good predictions of peak vertical GRF, stance time, contact length and vertical centre of mass displacement in running humans, if an appropriate spring stiffness is used. However, the model can be expected to overestimate peak horizontal GRF and mechanical work per step. When POFT occurs, the spring stiffness in the traditional spring-mass model is not equivalent to leg stiffness. Therefore, caution should be exercised when using spring stiffness to understand how the musculoskeletal system adapts to different running conditions. This can explain the contradictory results in the literature regarding the effect of running speed on leg stiffness.     

Predicting timing of foot strike during running,independent of striking technique,using principal component analysis of joint angles

As 3-dimensional (3D) motion-capture for clinical gait analysis continues to evolve, new methods must be developed to improve the detection of gait cycle events based on kinematic data. Recently, the application of principal component analysis (PCA) to gait data has shown promise in detecting important biomechanical features. Therefore, the purpose of this study was to define a new foot strike detection method for a continuum of striking techniques, by applying PCA to joint angle waveforms. In accordance with Newtonian mechanics, it was hypothesized that transient features in the sagittal-plane accelerations of the lower extremity would be linked with the impulsive application of force to the foot at foot strike. Kinematic and kinetic data from treadmill running were selected for 154 subjects, from a database of gait biomechanics. Ankle, knee and hip sagittal plane angular acceleration kinematic curves were chained together to form a row input to a PCA matrix. A linear polynomial was calculated based on PCA scores, and a 10-fold cross-validation was performed to evaluate prediction accuracy against gold-standard foot strike as determined by a 10 N rise in the vertical ground reaction force. Results show 89–94% of all predicted foot strikes were within 4 frames (20 ms) of the gold standard with the largest error being 28 ms. It is concluded that this new foot strike detection is an improvement on existing methods and can be applied regardless of whether the runner exhibits a rearfoot, midfoot, or forefoot strike pattern.     

Straightening of curved pattern of collagen fibers under load controls aortic valve shape

The network of collagen fibers in the aortic valve leaflet is believed to play an important role in the strength and durability of the valve. However, in addition to its stress-bearing role, such a fiber network has the potential to produce functionally important shape changes in the closed valve under pressure load. We measured the average pattern of the collagen network in porcine aortic valve leaflets after staining for collagen. We then used finite element simulation to explore how this collagen pattern influences the shape of the closed valve. We observed a curved or bent pattern, with collagen fibers angled downward from the commissures toward the center of the leaflet to form a pattern that is concave toward the leaflet free edge. Simulations showed that these curved fiber trajectories straighten under pressure load, leading to functionally important changes in closed valve shape. Relative to a pattern of straight collagen fibers running parallel to the leaflet free edge, the concave pattern of curved fibers produces a closed valve with a 40% increase in central leaflet coaptation height and with decreased leaflet billow, resulting in a more physiological closed valve shape. Furthermore, simulations show that these changes in loaded leaflet shape reflect changes in leaflet curvature due to modulation of in-plane membrane stress resulting from straightening of the curved fibers. This effect appears to play an important role in normal valve function and may have important implications for the design of prosthetic and tissue engineered replacement valves.     

Three-dimensional Orthotropic Viscoelastic Finite Element Model of a Human Ligament

Ligaments undergo finite strain displaying hyperelastic behaviour as the initially tangled fibrils present straighten out, combined with viscoelastic behaviour (strain rate sensitivity). In the present study the anterior cruciate ligament of the human knee joint is modelled in three dimensions to gain an understanding of the stress distribution over the ligament due to motion imposed on the ends, determined from experimental studies. A three dimensional, finite strain material model of ligaments has recently been proposed by Pioletti in Ref. [2]. It is attractive as it separates out elastic stress from that due to the present strain rate and that due to the past history of deformation. However, it treats the ligament as isotropic and incompressible. While the second assumption is reasonable, the first is clearly untrue. In the present study an alternative model of the elastic behaviour due to Bonet and Burton (Ref. [4]) is generalized. Bonet and Burton consider finite strain with constant modulii for the fibres and for the matrix of a transversely isotropic composite. In the present work, the fibre modulus is first made to increase exponentially from zero with an invariant that provides a measure of the stretch in the fibre direction. At 12% strain in the fibre direction, a new reference state is then adopted, after which the material modulus is made constant, as in Bonet and Burton's model. The strain rate dependence can be added, either using Pioletti's isotropic approximation, or by making the effect depend on the strain rate in the fibre direction only.

A solid model of a ligament is constructed, based on experimentally measured sections, and the deformation predicted using explicit integration in time. This approach simplifies the coding of the material model, but has a limitation due to the detrimental effect on stability of integration of the substantial damping implied by the nonlinear dependence of stress on strain rate. At present, an artificially high density is being used to provide stability, while the dynamics are being removed from the solution using artificial viscosity. The result is a quasi-static solution incorporating the effect of strain rate. Alternate approaches to material modelling and integration are discussed, that may result in a better model.     

The effects of morphological irregularity on the mechanical behavior of interdigitated biological sutures under tension

Irregular interdigitated morphology is prevalent in biological sutures in nature. Suture complexity index has long been recognized as the most important morphological parameter to govern the mechanical properties of biological sutures. However, the suture complexity index alone does not reflect all aspects of suture morphology. The goal of this investigation was to determine that besides suture complexity index, whether the degree of morphological irregularity of biological sutures has influences on the mechanical properties, and if there is any, how to quantify these influences. To explore these issues, theoretical and finite element (FE) suture models with the same suture complexity index but different levels of morphological irregularity were developed. The quasi-static stiffness, strength for damage initiation and post-failure process of irregular sutures were studied. It was shown that for the same suture complexity index, when the level of morphological irregularity increases, the overall strain to failure will increase while tensile stiffness is retained; also, the total energy to fracture increases with a sacrifice in strength to damage initiation. These results reveal that morphological irregularity is another important independent parameter to govern and balance the mechanical properties of biological sutures. Therefore, from the mechanics point of view, the prevalence of irregular suture morphology in nature is a merit, not a defect.     

口腔颌面部爆炸伤的研究现状

口腔颌面部爆炸伤是指由致伤物爆炸所造成的口腔颌面部组织损伤,两种主要致伤因素是冲击波和高速破片,较一般火器伤而言,爆炸伤的致伤机制及致伤特点都有不同之处。对爆炸伤害的物理机制和病理生理反应方面研究有助于改进防护及改善治疗策略。本文在简述口腔颌面部爆炸伤的致伤因素和损伤特点的基础上,着重综述了动物模型和有限元模型的研究方法及结果,旨在为以后的模型研究提供思路和参考。