Off-axis loads cause failure of the distal radius at lower magnitudes than axial loads: a finite element analysis

Distal radius (Colles') fractures are a common fall-related injury in older adults and frequently result in long-term pain and reduced ability to perform activities of daily living. Because the occurrence of a fracture during a fall depends on both the strength of the bone and upon the kinematics and kinetics of the impact itself, we sought to understand how changes in bone mineral density (BMD) and loading direction affect the fracture strength and fracture initiation location in the distal radius. A three-dimensional finite element model of the radius, scaphoid, and lunate was used to examine changes of +/-2% and +/-4% BMD, and both axial and physiologically relevant off-axis loads on the radius. Changes in BMD resulted in similar percent changes in fracture strength. However, modifying the applied load to include dorsal and lateral components (assuming a dorsal view of the wrist, rather than an anatomic view) resulted in a 47% decrease in fracture strength (axial failure load: 2752N, off-axis: 1448N). Loading direction also influenced the fracture initiation site. Axially loaded radii failed on the medial surface immediately proximal to the styloid process. In contrast, off-axis loads, containing dorsal and lateral components, caused failure on the dorsal-lateral surface. Because the radius appears to be very sensitive to loading direction, the results suggest that much of the variability in fracture strength seen in cadaver studies may be attributed to varying boundary conditions. The results further suggest that interventions focused on reducing the incidence of Colles' fractures when falls onto the upper extremities are unavoidable may benefit from increasing the extent to which the radius is loaded along its axis.     

Strain-rate sensitivity of the rabbit MCL diminishes at traumatic loading rates

This study was performed to determine whether the viscoelastic behavior of ligaments persists at high rates of loading, such as those associated with sports-related trauma or motor vehicle accidents. Medial collateral ligaments (MCLs) from 22 skeletally mature New Zealand White rabbits were tensile tested quasi-statically and via two impact conditions at displacement rates of 0.17 mm/s (n=22), 640+/-160 mm/s (n=10) and 2500+/-270 mm/s (n=12) (corresponding to strain rates of approximately 1.0%/s, 3660%/s and 14,000%/s, respectively). Despite dramatic increases in displacement rate, only a modest strain-rate effect was observed when the specimens tested quasi-statically were compared to those tested via impact (24% and 37% increases in stiffness and failure load, respectively). There were no differences in the structural (e.g. 145+/-30 and 136+/-29 N/mm stiffness values, respectively) or failure properties (e.g. 434+/-91 and 443+/-154 N failure load values, respectively) of the two impact-tested groups. Our findings suggest that the rabbit MCL is not viscoelastic at loading rates approximating those associated with high-energy trauma.     

Effect of boundary conditions,impact loading and hydraulic stiffening on femoral fracture strength

Patient specific quantitative CT (QCT) imaging data together with the finite element (FE) method may provide an accurate prediction of a patient's femoral strength and fracture risk. Although numerous FE models investigating femoral fracture strength have been published, there is little consent on the effect of boundary conditions, dynamic loading and hydraulic strengthening due to intra-medullary pressure on the predicted fracture strength. We developed a QCT-derived FE model of a proximal femur that included node-specific modulus assigned based on the local bone density. The effect of three commonly used boundary conditions published in literature were investigated by comparing the resulting strain field due to an applied fracture load. The models were also augmented with viscoelastic material properties and subject to a realistic impact load profile to determine the effect of dynamic loads on the strain field. Finally, the effect of hydraulic strengthening was investigated by including node specific permeability and performing a coupled pore diffusion and stress analysis of the FE model. Results showed that all boundary conditions yield the same strain field patterns, but peak strains were 22% lower and fracture load was 18% higher when loaded at the greater trochanter than when loaded at the femoral head. Comparison of the dynamic models showed that material viscoelasticity was important, but inertial effects (vibration and shock) were not. Finally, pore pressure changes did not cause significant hydraulic strengthening of bone under fall impact loading.     

Rate-independent characteristics of an arthroscopically implantable force probe in the human achilles tendon

The effect of loading rate on specimen calibration was investigated for an implantable force sensor of the two-point loading variety. This variety of sensor incorporates a strain gage to measure the compressive load applied to the sensor due to tensile loading in a soft tissue specimen. The Achilles tendon in each of four human cadaveric lower extremities was instrumented with a force sensor and then loaded in tension using a materials testing machine. Each specimen was tensile tested at three different displacement rates, 0.25, 2.5 and 12.7 cm s(-1), corresponding with mean loading rates of 33.8, 513.2, and 2838.6 N s(-1), respectively. A calibration curve relating the force sensor signal and applied tendon tension was generated for each specimen/ displacement rate combination. For each specimen, calibration curves were compared by calculating an RMS error for the entire data set (eRMS = 1.6% of the full load value) and a coefficient of determination, R2, of a curve fit through all of the data (R2 = 99.6%). Over the range of rates tested, no measurable change in sensor sensitivity due to loading rate was observed. Hysteresis for all displacement rates was on the order of 2.4%.     

Investigation of impact loading rate effects on the ligamentous cervical spinal load-partitioning using finite element model of functional spinal unit C2–C3

The cervical spine functions as a complex mechanism that responds to sudden loading in a unique manner, due to intricate structural features and kinematics. The spinal load-sharing under pure compression and sagittal flexion/extension at two different impact rates were compared using a bio-fidelic finite element (FE) model of the ligamentous cervical functional spinal unit (FSU) C2–C3. This model was developed using a comprehensive and realistic geometry of spinal components and material laws that include strain rate dependency, bone fracture, and ligament failure. The range of motion, contact pressure in facet joints, failure forces in ligaments were compared to experimental findings. The model demonstrated that resistance of spinal components to impact load is dependent on loading rate and direction. For the loads applied, stress increased with loading rate in all spinal components, and was concentrated in the outer intervertebral disc (IVD), regions of ligaments to bone attachment, and in the cancellous bone of the facet joints. The highest stress in ligaments was found in capsular ligament (CL) in all cases. Intradiscal pressure (IDP) in the nucleus was affected by loading rate change. It increased under compression/flexion but decreased under extension. Contact pressure in the facet joints showed less variation under compression, but increased significantly under flexion/extension particularly under extension. Cancellous bone of the facet joints region was the only component fractured and fracture occurred under extension at both rates. The cervical ligaments were the primary load-bearing component followed by the IVD, endplates and cancellous bone; however, the latter was the most vulnerable to extension as it fractured at low energy impact.     

Foot strike pattern,step rate,and trunk posture combined gait modifications to reduce impact loading during running

Elevated impact loading can be detrimental to runners as it has been linked to the increased risk of tibial stress fracture and plantar fasciitis. The objective of this study was to investigate the combined effects of foot strike pattern, step rate, and anterior trunk lean gait modifications on impact loading in runners. Nineteen healthy runners performed 12 separate gait modification trials involving: three foot strike patterns (rearfoot, midfoot, and forefoot strike), two step rates (natural and 10% increased), and two anterior trunk lean postures (natural and 10-degree increased flexion). Overall, forefoot strike combined with increased step rate led to the lowest impact loading rates, and rearfoot strike combined with anterior trunk lean led to the highest impact loading rates. In addition, there were interaction effects between foot strike pattern and step rate on awkwardness and effort, such that it was both more natural and easier to transition to a combined gait modification involving forefoot strike and increased step rate than to an isolated gait modification involving either forefoot strike or increased step rate. These findings could help to inform gait modifications for runners to reduce impact loading and associated injury risks.     

Perimenopausal risk of falling and incidence of distal forearm fracture.

A postal survey of 2000 women and 2000 men sampled from the electoral roll in Oxford was undertaken to ascertain whether changes with age in the risk of falling might explain the stepwise increases in age specific incidence rates of distal forearm fracture which occur in women at around the age of 50. Corrected response rates were 83% for women and 72% for men. In women, but not in men, there was a rise in the risk of falling from 45 years, peaking in the 55-59 year age group, and sinking to a nadir at ages 70-74. In both sexes rates rose in extreme old age. These variations were not attributable to preferential response from people who had suffered a fracture. It is concluded that changes in the risk of falling interact with osteoporosis to produce a perimenopausal rise in the incidence of forearm fractures and contribute to the fluctuations in incidence of these fractures in old age.     

Load transfer analysis of the distal radius from in-vivo high-resolution CT-imaging.

Prevention of osteoporotic bone fractures requires accurate diagnostic methods to detect the increase in bone fragility at an early stage of osteoporosis. However, today's bone fracture risk prediction, primarily based on bone density measurement, is not sufficiently precise. There is increasing evidence that, in addition to bone density, also the bone microarchitecture and its mechanical loading conditions are important factors determining the fracture risk. Recently, it has been shown that new high-resolution imaging techniques in combination with new computer modeling techniques based on the finite-element (FE) method can account for these additional factors. These techniques might provide information that is more relevant for the prediction of bone fracture risk. So far, however, these new imaged-based FE techniques have not been feasible in-vivo. The objectives of this study were to quantify the load transfer through the trabecular network in a distal radius using a computer model based on in-vivo high-resolution images and to determine if common regions of fractures can be explained as a result of high tissue loading in these regions. The left distal radius and the two adjacent carpal bones of a healthy volunteer were imaged using a high-resolution three-dimensional CT system providing an isotropic resolution of 165 microm. The bone representation was converted into a FE-model that was used to calculate stresses and strains in the trabecular network. The two carpal bones were loaded using different load ratios (for each load case 1000 N in total) representing impact forces on the hand either in near-neutral position or ulnar/radial deviation. The load transfer through the trabecular network of the radius was characterized by the tissue strain energy density (SED) distribution for all load cases. It was found that the distribution of the tissue loading depends on the ratio of the forces acting on the carpal bones. For all load cases the higher SED values (on average: 0.02 +/- 0.08 (S.D.) N mm(-2)) are found in a 10 mm region adjacent to the articular surface which corresponds well with the region where Colles- or Chauffeur-fractures occur. We expect that, eventually, this new approach can lead to a better prediction of the fracture risk than methods based on bone density alone since it accounts for the bone microstructure as well as its loading conditions.     

An ex vivo experiment to reproduce a forward fall leading to fractured and non-fractured radii

Forward falls represent a risk of injury for the elderly. The risk is increased in elderly persons with bone diseases, such as osteoporosis. However, half of the patients with fracture were not considered at risk based on bone density measurement (current clinical technique). We assume that loading conditions are of high importance and should be considered. Real loading conditions in a fall can reach a loading speed of 2 m/s on average. The current study aimed to apply more realistic loading conditions that simulate a forward fall on the radius ex vivo. Thirty radii from elderly donors (79 y.o. ± 12 y.o., 15 males, 15 females) were loaded at 2 m/s using a servo-hydraulic testing machine to mimic impact that corresponds to a fall. Among the 30 radii, 14 had a fracture after the impact, leading to two groups (fractured and non-fractured). Surfacic strain fields were measured using stereovision and allow for visualization of fracture patterns. The average maximum load was 2963 ± 1274 N. These experimental data will be useful for assessing the predictive capability of fracture risk prediction methods such as finite element models.     

Mechanical behaviour of in-situ chondrocytes subjected to different loading rates: a finite element study

Experimental findings indicate that in-situ chondrocytes die readily following impact loading, but remain essentially unaffected at low (non-impact) strain rates. This study was aimed at identifying possible causes for cell death in impact loading by quantifying chondrocyte mechanics when cartilage was subjected to a 5% nominal tissue strain at different strain rates. Multi-scale modelling techniques were used to simulate cartilage tissue and the corresponding chondrocytes residing in the tissue. Chondrocytes were modelled by accounting for the cell membrane, pericellular matrix and pericellular capsule. The results suggest that cell deformations, cell fluid pressures and fluid flow velocity through cells are highest at the highest (impact) strain rate, but they do not reach damaging levels. Tangential strain rates of the cell membrane were highest at the highest strain rate and were observed primarily in superficial tissue cells. Since cell death following impact loading occurs primarily in superficial zone cells, we speculate that cell death in impact loading is caused by the high tangential strain rates in the membrane of superficial zone cells causing membrane rupture and loss of cell content and integrity.     

Colles' fracture; a method of maintaining reduction

A major difficulty in the treatment of Colles' fracture is maintenance of reduction. Wedging of the cast, a procedure used in dealing with other orthopedic conditions, was adapted to the treatment of Colles' fracture and was employed in 23 cases. In most of them the method was effective in preventing displacement.     

Specimen-specific modeling of hip fracture pattern and repair

Hip fracture remains a major health problem for the elderly. Clinical studies have assessed fracture risk based on bone quality in the aging population and cadaveric testing has quantified bone strength and fracture loads. Prior modeling has primarily focused on quantifying the strain distribution in bone as an indicator of fracture risk. Recent advances in the extended finite element method (XFEM) enable prediction of the initiation and propagation of cracks without requiring a priori knowledge of the crack path. Accordingly, the objectives of this study were to predict femoral fracture in specimen-specific models using the XFEM approach, to perform one-to-one comparisons of predicted and in vitro fracture patterns, and to develop a framework to assess the mechanics and load transfer in the fractured femur when it is repaired with an osteosynthesis implant. Five specimen-specific femur models were developed from in vitro experiments under a simulated stance loading condition. Predicted fracture patterns closely matched the in vitro patterns; however, predictions of fracture load differed by approximately 50% due to sensitivity to local material properties. Specimen-specific intertrochanteric fractures were induced by subjecting the femur models to a sideways fall and repaired with a contemporary implant. Under a post-surgical stance loading, model-predicted load sharing between the implant and bone across the fracture surface varied from 59%:41% to 89%:11%, underscoring the importance of considering anatomic and fracture variability in the evaluation of implants. XFEM modeling shows potential as a macro-level analysis enabling fracture investigations of clinical cohorts, including at-risk groups, and the design of robust implants.     

Damage accumulation in vertebral trabecular bone depends on loading mode and direction

Osteoporotic vertebral fractures constitute a major clinical problem in ageing societies. A third of all vertebral fractures is caused by falls, 15% by lifting heavy loads or traffic accidents and over 50% are not relatable to a traumatic event. In the latter case vertebrae show sinter processes which indicate the accumulation of damage and permanent deformation. Accumulated damage may not be visible on radiographs but increases the risk of fracture and could lead to vertebral collapse. Clear understanding of the accumulation of damage and residual strains and their dependence on loading mode and direction is important for understanding vertebral fractures. Altogether, 251 cylindrical samples (8×18-25mm) were obtained from 50 male and 54 female fresh frozen human vertebrae (T1-L3) of 65 (21-94) years. Vertebrae were randomly assigned to three groups cranial-caudal, anterior-posterior and latero-lateral. Specimens were mechanically loaded in compression, tension or torsion in five load steps at a strain rate of 0.2%/s. Three conditioning cycles were driven per load step. Stress-strain curves were reconstructed from the force-displacement or from the moment-twist angle curves. Damage accumulated from 0 to 86% in compression, from 0 to 76% in tension and from 0 to 86% in torsion through the five load steps. Residual strains accumulated from 0 to -0.008mm/mm in compression, 0 to 0.006mm/mm in tension and 0 to 0.026rad/rad in torsion. Significantly less damage (p<0.05) but not residual strains accumulated in transverse directions. This study provides detailed experimental insights into the damage behaviour of vertebral trabecular bone under various loads occurring in vivo. Damage but not residual strain evolution seems to be anisotropic. Both seem to evolve differently under different loading modes. The results could be of importance in understanding vertebral fractures.     

Relationships between femoral fracture loads for two load configurations

Studies of proximal femoral strength usually involve one of two types of loading conditions, loading similar to joint loading during single-limb stance or loading simulating impact from a fall. When interpreting the results of studies involving only one of these load configurations, the question arises as to their applicability to the other configuration. In addition, it is desirable to know whether, for an individual bone, fracture load for one load configuration is indicative of fracture load for the other configuration. In this study, the relationship between proximal femoral fracture loads for single-limb stance loading and loading simulating impact from a type of fall was determined from mechanical testing of 17 matched pairs of human proximal femora. Fracture loads for these two configurations were found to be linearly related (r = 0.901, p < 0.001). However, the correlation between fracture loads is not notably stronger than correlations currently available between fracture load and measures of bone density and geometry. In addition, the regression results indicate that 81% of the variance in fracture load for one loading condition is accounted for by fracture load for the other loading condition. Thus, 19% of the variance remains unexplained, indicating that the results of studies involving only one load configuration are not necessarily indicative of those that would be found for another configuration.     

股骨近端抗旋髓内钉法治疗老年股骨转子下粉碎性骨折的效果分析

目的:探讨利用股骨近端抗旋髓内钉法治疗老年股骨转子下粉碎性骨折的临床效果,为临床提供参考。方法:对我院2009年6月-2013年1月收治的38例老年股骨转子下粉碎性骨折患者行股骨近端抗旋髓内钉法进行手术治疗,分析手术方法、效果及患者的预后效果。结果:28例行闭合复位,7例行骨折端切开复位,患者术后7-14d可下地活动,3 m内扶拐部分负重行走,31例(88.6%)获得随访,随访时间9-24个月,骨折愈合时间为3-6个月,平均3.7个月;Harris髋功能评分标准:优18例,良10例,中3例,优良率90.32%(28/31)。结论:PFNA是一种治疗老年股骨转子下粉碎性骨折的有效装置,能够减少骨折不愈合、髋内翻畸形愈合及内固定断裂、切割股骨头等并发症,在治疗老年股骨转子下粉碎性骨折时可达到较高的骨愈合率、较快的功能恢复。     

Proximal femur elastic behaviour is the same in impact and constant displacement rate fall simulation

Understanding proximal femur fracture may yield new targets for fracture prevention screening and treatment. The goal of this study was to characterize force–displacement and failure behaviours in the proximal femur between displacement control and impact loading fall simulations. Twenty-one human proximal femurs were tested in two ways, first to a sub-failure load at a constant displacement rate, then to fracture in an impact fall simulator. Comparisons of sub-failure energy and stiffness were made between the tests at the same compressive force. Additionally, the impact failure tests were compared with previous, constant displacement rate failure tests (at 2 and 100 mm/s) in terms of energy, yield force, and stiffness. Loading and displacement rates were characterized and related to specimen stiffness in the impact tests. No differences were observed between the sub-failure constant displacement and impact tests in the aforementioned metrics. Comparisons between failure tests showed that the impact group had the lowest absorbed energy, 24% lower maximum force and 160% higher stiffness than the 100 mm/s group (p<0.01$p<0.01$ for all), but suffered from low statistical power to differentiate the donor age and specimen BMD. Loading and displacement rates for the specimens tested using impact varied during each test and between specimens and did not show appreciable viscoelasticity. These results indicate that constant displacement rate testing may help understand sub-failure mechanical behaviour, but may not elucidate failure behaviours. The differences between the impact and constant displacement rate fall simulations have important ramifications for interpreting the results of previous experiments.     

Local and global measurements show that damage initiation in articular cartilage is inhibited by the surface layer and has significant rate dependence

Cracks in articular cartilage are a common sign of joint damage, but failure properties of cartilage are poorly understood, especially for damage initiation. Cartilage failure may be further complicated by rate-dependent and depth-dependent properties, including the compliant surface layer. Existing blunt impact methods do not resolve local cartilage inhomogeneities and traditional fracture mechanics tests induce crack blunting and may violate underlying assumptions of linear elasticity. To address this knowledge gap, we developed and applied a method to indent cartilage explants with a sharp blade and initiate damage across a range of loading rates (strain rates 0.5%/s–500%/s), while recording local sample deformation and strain energy fields using confocal elastography. To investigate the importance of cartilage’s compliant surface, we repeated the experiment for samples with the surface removed. Bulk data suggest a critical force at which the tissue cuts, but local strains reveals that the deformation the sample can sustain before reaching this force is significantly higher in the surface layer. Bulk and local results also showed significant rate dependence, such that samples were easier to cut at faster speeds. This result highlights the importance of rate for understanding cracks in cartilage and parallels recent studies of rate-dependent failure in hydrogels. Notably, local sample deformation fields were well fit by classical Hookean elasticity. Overall, this study illustrates how local and global measurements surrounding the initiation of damage in articular cartilage can be combined to reveal the importance of cartilage’s zonal structure in protecting against failure across physiologically relevant loading rates.     

Inelastic strain accumulation in cortical bone during rapid transient tensile loading

An experimental study examined the tensile stress-strain behavior of cortical bone during rapid load cycles to high strain amplitudes. Machined bovine and human cortical bone samples were subjected to loading cycles at a nominal load/unload rate of +/- 420 MPa/s. Loads were reversed at pre selected strain levels such that load cycles were typically completed in 0.5-0.7 seconds. Axial strain behavior demonstrated considerable nonlinearity in the first load cycle, while transverse strain behavior was essentially linear. For the human bone 29.1 percent (S.D. = 4.7 percent), and for the bovine bone 35.1 percent (S.D. = 10.8 percent) of the maximum nonlinear strain accumulated after load reversal, where nonlinear strain was defined as the difference between total strain and strain corresponding to linear elastic behavior. Average residual axial strain on unloading was 35.4 percent (S.D. = 1.2 percent) for human bone and 35.1 percent (S.D. = 2.9 percent) of maximum nonlinear strain. Corresponding significant volumetric strains and residual volumetric strains were found. The results support the conclusions that the nonlinear stress-strain behavior observed during creep loading also occurs during transient loading at physiological rates. The volume increases suggest that damage accumulation, i.e., new internal surfaces and voids, plays a major role in this behavior. The residual volume increases and associated disruptions in the internal structure of bone provide a potential stimulus for a biological repair response.     

Accuracy of finite element predictions in sideways load configurations for the proximal human femur

Subject-specific finite element models have been used to predict stress-state and fracture risk in individual patients. While many studies analysed quasi-axial loading configurations, only few works simulated sideways load configurations, such as those arising in a fall. The majority among these latter directly predicted bone strength, without assessing elastic strain prediction accuracy. The aim of the present work was to evaluate if a subject-specific finite element modelling technique from CT data that accurately predicted strains in quasi-axial loading configurations is suitable to accurately predict strains also when applying low magnitude loads in sideways configurations. To this aim, a combined numerical-experimental study was performed to compare finite element predicted strains with strain-gauge measurements from three cadaver proximal femurs instrumented with sixteen strain rosettes and tested non-destructively under twelve loading configurations, spanning a wide cone (0-30° for both adduction and internal rotation angles) of sideways fall scenarios. The results of the present study evidenced a satisfactory agreement between experimentally measured and predicted strains (R(2) greater than 0.9, RMSE% lower than 10%) and displacements. The achieved strain prediction accuracy is comparable to those obtained in state of the art studies in quasi-axial loading configurations. Still, the presence of the highest strain prediction errors (around 30%) in the lateral neck aspect would deserve attention in future studies targeting bone failure.     

The axial injury tolerance of the human foot/ankle complex and the effect of Achilles tension

Axial loading of the foot/ankle complex is an important injury mechanism in vehicular trauma that is responsible for severe injuries such as calcaneal and tibial pilon fractures. Axial loading may be applied to the leg externally, by the toepan and/or pedals, as well as internally, by active muscle tension applied through the Achilles tendon during pre-impact bracing. The objectives of this study were to investigate the effect of Achilles tension on fracture mode and to empirically model the axial loading tolerance of the foot/ankle complex. Blunt axial impact tests were performed on forty-three (43) isolated lower extremities with and without experimentally simulated Achilles tension. The primary fracture mode was calcaneal fracture in both groups. However, fracture initiated at the distal tibia more frequently with the addition of Achilles tension (p < 0.05). Acoustic sensors mounted to the bone demonstrated that fracture initiated at the time of peak local axial force. A survival analysis was performed on the injury data set using a Weibull regression model with specimen age, gender, body mass, and peak Achilles tension as predictor variables (R2 = 0.90). A closed-form survivor function was developed to predict the risk of fracture to the foot/ankle complex in terms of axial tibial force. The axial tibial force associated with a 50% risk of injury ranged from 3.7 kN for a 65 year-old 5th percentile female to 8.3 kN for a 45 year-old 50th percentile male, assuming no Achilles tension. The survivor function presented here may be used to estimate the risk of foot/ankle fracture that a blunt axial impact would pose to a human based on the peak tibial axial force measured by an anthropomorphic test device.