Achilles tendon stress is more sensitive to subject-specific geometry than subject-specific material properties: A finite element analysis

This study used subject-specific measures of three-dimensional (3D) free Achilles tendon geometry in conjunction with a finite element method to investigate the effect of variation in subject-specific geometry and subject-specific material properties on tendon stress during submaximal isometric loading. Achilles tendons of eight participants (Aged 25–35 years) were scanned with freehand 3D ultrasound at rest and during a 70% maximum voluntary isometric contraction. Ultrasound images were segmented, volume rendered and transformed into subject-specific 3D finite element meshes. The mean (±SD) lengths, volumes and cross-sectional areas of the tendons at rest were 62 ± 13 mm, 3617 ± 984 mm³ and 58 ± 11 mm² respectively. The measured tendon strain at 70% MVIC was 5.9 ± 1.3%. Subject-specific material properties were obtained using an optimisation approach that minimised the difference between measured and modelled longitudinal free tendon strain. Generic geometry was represented by the average mesh and generic material properties were taken from the literature. Local stresses were subsequently computed for combinations of subject-specific and generic geometry and material properties. For a given geometry, changing from generic to subject-specific material properties had little effect on the stress distribution in the tendon. In contrast, changing from generic to subject-specific geometry had a 26-fold greater effect on tendon stress distribution. Overall, these findings indicate that the stress distribution experienced by the living free Achilles tendon of a young and healthy population during voluntary loading are more sensitive to variation in tendon geometry than variation in tendon material properties.     

Mechanics of the right whale mandible: Full scale testing and finite element analysis

In an effort to better understand the mechanics of ship-whale collision and to reduce the associated mortality of the critically endangered North Atlantic right whale, a comprehensive biomechanical study has been conducted by the Woods Hole Oceanographic Institution and the University of New Hampshire. The goal of the study is to develop a numerical modeling tool to predict the forces and stresses during impact and thereby the resulting mortality risk to whales from ship strikes.Based on post-mortem examinations, jaw fracture was chosen as a fatal endpoint for the whales hit by a vessel. In this paper we investigate the overall mechanical behavior of a right whale mandible under transverse loading and develop a finite element analysis model of the bone. The equivalent elastic modulus of the cortical component of right whale mandible is found by comparing full-scale bending tests with the results of numerical modeling. The finite element model of the mandible can be used in conjunction with a vessel-whale collision event model to predict bone fracture for various ship strike scenarios.     

Three-dimensional finite element analysis of the composite and compomer onlays in primary molars

Resin onlay restoration is an esthetic alternative technique used for restoring extensively damaged primary molars. Understanding the behavior of materials under repeated functional stress and how the stress is transmitted to the remaining tooth structure is important. The aim of this study was to compare stresses in primary molars restored with indirect composite and compomer onlay. 3D frame models of the right mandibular and maxillary primary molars and the alveolar bone were created using computerized tomography images of a six-year-old girl. The enamel and dentine layers above the cement layer were unified to generate onlay restoration, and composite and compomer were used as restorative materials. The vertical occlusal load (100?N) was applied to the teeth in the occlusal contact areas. The von Mises stress distributions and normal stress distributions of the y-axis (parallel to the long axis of tooth) were evaluated. The occlusal stress is transmitted to the cervical part of healthy teeth by spreading it through the enamel layer. The composite and compomer restorative materials exhibited similar stress distribution patterns. An indirect technique creates a structure similar to the original morphological form, and it allows restorations to distribute high occlusal stresses and to minimize possible breakages.     

Multi-scale analysis of optic chiasmal compression by finite element modelling

The precise mechanism of bitemporal hemianopia (a type of partial visual field defect) is still not clear. Previous work has investigated this problem by studying the biomechanics of chiasmal compression caused by a pituitary tumour growing up from below the optic chiasm. A multi-scale analysis was performed using finite element models to examine both the macro-scale behaviour of the chiasm and the micro-scale interactions of the nerve fibres within it using representative volume elements. Possible effects of large deflection and non-linear material properties were incorporated. Strain distributions in the optic chiasm and optic nerve fibres were obtained from these models. The results of the chiasmal model agreed well with the limited experimental results available, indicating that the finite element modelling can be a useful tool for analysing chiasmal compression. Simulation results showed that the strain distribution in nasal (crossed) nerve fibres was much more nonuniform and locally higher than in temporal (uncrossed) nerve fibres. This strain difference between nasal and temporal nerve fibres may account for the phenomenon of bitemporal hemianopia.     

FEBio finite element models of the human lumbar spine

Finite element analysis has proven to be a viable method for assessing many structure-function relationships in the human lumbar spine. Several validated models of the spine have been published, but they typically rely on commercial packages and are difficult to share between labs. The goal of this study is to present the development of the first open-access models of the human lumbar spine in FEBio. This modeling framework currently targets three deficient areas in the field of lumbar spine modeling: 1) open-access models, 2) accessibility for multiple meshing schemes, and 3) options to include advanced hyperelastic and biphasic constitutive models.     

A finite element analysis of the retinal hemorrhages accompanied by shaken baby syndrome/abusive head trauma

We aimed to elucidate the mechanism of the retinal hemorrhage (RH) accompanied by shaken baby syndrome or abusive head trauma (SBS/AHT) by analyses using a computational model. We focused on a hypothesis that the vitreoretinal traction due to acceleration and deceleration caused by abusive shaking leads to retinal hemorrhage. A finite element (FE) mechanical model with simple spherical geometry was constructed. When the FE mechanical model was virtually shaken, the intensity of the stress applied to the retinal plane agreed well with the results from an analysis using a physical model made of agar gel. Impacts due to falling events induced more intensive tensile stresses, but with shorter duration, than the shake did. By applying a mathematical theory on tackiness, we propose a hypothesis that the time integration of the stress, in the unit of Pa·s, would be a good predictor of the RH accompanied by SBS/AHT. A single cycle of abusive shake amounted to 101 Pa·s of time integration of inflicted stress, while a single impact event amounted to 36 Pa·s. This would explain the paradoxical observation that shaking induces RH while RH due to impact events is only seen in a major event such as a fatal motor vehicle accident.     

Stress analysis of irradiated human tooth enamel using finite element methods

The objectives of this project were to use finite element methods to determine how changes in the elastic modulus due to oral cancer therapeutic radiation alter the distribution of mechanical stresses in teeth and to determine if observed failures in irradiated teeth correlate with changes in mechanical stresses. A thin slice section finite element (FE) model was constructed from micro CT sections of a molar tooth using MIMICS and 3-Matic software. This model divides the tooth into three enamel regions, the dentin-enamel junction (DEJ) and dentin. The enamel elastic modulus was determined in each region using nano indentation for three experimental groups namely – control (non-radiated), in vitro irradiated (simulated radiotherapy following tooth extraction) and in vivo irradiated (extracted subsequent to oral cancer patient radiotherapy) teeth. Physiological loads were applied to the tooth models at the buccal and lingual cusp regions for all three groups (control, in vitro and in vivo). The principal tensile stress and the maximum shear stress were used to compare the results from different groups since it has been observed in previous studies that delamination of enamel from the underlying dentin was one of the major reasons for the failure of teeth following therapeutic radiation. From the FE data, we observed an increase in the principal tensile stress within the inner enamel region of in vivo irradiated teeth (9.97 ± 1.32 MPa) as compared to control/non-irradiated teeth (8.44 ± 1.57 MPa). Our model predicts that failure occurs at the inner enamel/DEJ interface due to extremely high tensile and maximum shear stresses in in vivo irradiated teeth which could be a cause of enamel delamination due to radiotherapy.     

Comparison of stress distribution in dental crown with different cusp angles: 3D finite element analysis

The purpose of this study is to measure the failure risk of a crown depending on the cusp angle. Three all-ceramic crown models consisting of C_H (high incline), C_M (middle incline), and C_L (low incline) are designed. Stress is applied to the crown with Loading case-1 (top of cusp tip) and Loading case-2 (middle of cusp ridge) with the use of FEA software. In Loading case-1 and case-2, the C_H showed the highest Maximum Principal Stress (MPS) while the CL showed the lowest MPS. The cusp angle is an influential factor affecting stress distribution in dental crowns.     

Three-dimensional finite element simulations of the mechanical response of the fingertip to static and dynamic compressions

The analysis of the mechanics of the contact interactions of fingers/handle and the stress/strain distributions in the soft tissues in the fingertip is essential to optimize design of tools to reduce many occupation-related hand disorders. In the present study, a three-dimensional (3D) finite element (FE) model for the fingertip is proposed to simulate the nonlinear and time-dependent responses of a fingertip to static and dynamic loadings. The proposed FE model incorporates the essential anatomical structures of a finger: skin layers (outer and inner skins), subcutaneous tissue, bone and nail. The soft tissues (inner skin and subcutaneous tissue) are considered to be nonlinearly viscoelastic, while the hard tissues (outer skin, bone and nail) are considered to be linearly elastic. The proposed model has been used to simulate two loading scenarios: (a) the contact interactions between the fingertip and a flat surface and (b) the indentation of the fingerpad via a sharp wedge. For case (a), the predicted force/displacement relationships and time-dependent force responses are compared with the published experimental data; for case (b), the skin surface deflection profiles were predicted and compared with the published experimental observations. Furthermore, for both cases, the time-dependent stress/strain distributions within the tissues of the fingertip were calculated. The good agreement between the model predictions and the experimental observations indicates that the present model is capable of predicting realistic time-dependent force/displacement responses and stress/strain distributions in the soft tissues for dynamic loading conditions.     

On the importance of 3D,geometrically accurate,and subject-specific finite element analysis for evaluation of in-vivo soft tissue loads

Pressure ulcers are a type of local soft tissue injury due to sustained mechanical loading and remain a common issue in patient care. People with spinal cord injury (SCI) are especially at risk of pressure ulcers due to impaired mobility and sensory perception. The development of load improving support structures relies on realistic tissue load evaluation e.g. using finite element analysis (FEA). FEA requires realistic subject-specific mechanical properties and geometries. This study focuses on the effect of geometry. MRI is used for the creation of geometrically accurate models of the human buttock for three able-bodied volunteers and three volunteers with SCI. The effect of geometry on observed internal tissue deformations for each subject is studied by comparing FEA findings for equivalent loading conditions. The large variations found between subjects confirms the importance of subject-specific FEA.     

Reliable simulations of the human proximal femur by high-order finite element analysis validated by experimental observations

Background: The mechanical response of patient-specific bone to various load conditions is of major clinical importance in orthopedics. Herein we enhance the methods presented in Yosibash et al. [2007. A CT-based high-order finite element analysis of the human proximal femur compared to in-vitro experiments. ASME Journal of Biomechanical Engineering 129(3), 297–309.] for the reliable simulations of the human proximal femur by high-order finite elements (FEs) and validate the simulations by experimental observations.

Method of approach: A fresh-frozen human femur was scanned by quantitative computed tomography (QCT) and thereafter loaded (in vitro experiments) by a quasi-static force of up to 1250 N. QCT scans were manipulated to generate a high-order FE bone model with distinct cortical and trabecular regions having inhomogeneous isotropic elastic properties with Young's modulus represented by continuous spatial functions. Sensitivity analyses were performed to quantify parameters that mostly influence the mechanical response. FE results were compared to displacements and strains measured in the experiments.

Results: Young moduli correlated to QCT Hounsfield Units by relations in Keyak and Falkinstein [2003. Comparison of in situ and in vitro CT scan-based finite element model predictions of proximal femoral fracture load. Medical Engineering and Physics 25, 781–787.] were found to provide predictions that match the experimental results closely. Excellent agreement was found for both the displacements and strains. The presented study demonstrates that reliable and validated high-order patient-specific FE simulations of human femurs based on QCT data are achievable for clinical computer-aided decision making.     

Structural analysis of two different stent configurations

Two different stent configurations (i.e. the well known Palmaz–Schatz (PS) and a new stent configuration) are mechanically investigated. A finite element model was used to study the two geometries under combining loads and a computational fluid dynamic model based on fluid structure interaction was developed investigating the plaque and the artery wall reactions in a stented arterial segment. These models determine the stress and displacement fields of the two stents under internal pressure conditions. Results suggested that stent designs cause alterations in vascular anatomy that adversely affect arterial stress distributions within the wall, which have impact in the vessel responses such as the restenosis. The hemodynamic analysis shows the use of new stent geometry suggests better biofluid mechanical response such as the deformation and the progressive amount of plaque growth.     

Implication of obesity on motion,posture and internal stress of the foot: an experimental and finite element analysis

Obesity causes increased loading on the foot which can damage the soft tissue and bone ultimately leading to foot problems. Experimental and computational methods were used to analyse the chain of biomechanical changes in the lower limb due to obesity. The experimental study shows some changes in foot posture and gait where obese subjects were more likely to have pronated feet, smaller joint angles in the sagittal and frontal planes, smaller cadence, and smaller stride length. Anatomically correct finite element models generated on obese subjects showed increased and altered internal and plantar stress. Altered foot posture was identified as a key indicator of increased internal stress indicating the importance of foot posture correction.     

Muscle-driven finite element simulation of human foot movements

This paper describes a finite element scheme for realistic muscle-driven simulation of human foot movements. The scheme is used to simulate human ankle plantar flexion. A three-dimensional anatomically detailed finite element model of human foot and lower leg is developed and the idea of generating natural foot movement based entirely on the contraction of the plantar flexor muscles is used. The bones, ligaments, articular cartilage, muscles, tendons, as well as the rest soft tissues of human foot and lower leg are included in the model. A realistic three-dimensional continuum constitutive model that describes the biomechanical behaviour of muscles and tendons is used. Both the active and passive properties of muscle tissue are accounted for. The materials for bones and ligaments are considered as homogeneous, isotropic and linearly elastic, whereas the articular cartilage and the rest soft tissues (mainly fat) are defined as hyperelastic materials. The model is used to estimate muscle tissue deformations as well as stresses and strains that develop in the lower leg muscles during plantar flexion of the ankle. Stresses and strains that develop in Achilles tendon during such a movement are also investigated.     

Brief communication: comparing loading scenarios in lower first molar supporting bone structure using 3D finite element analysis

Finite element analysis (FEA) is a widespread technique to evaluate the stress/strain distributions in teeth or dental supporting tissues. However, in most studies occlusal forces are usually simplified using a single vector (i.e., point load) either parallel to the long tooth axis or oblique to this axis. In this pilot study we show how lower first molar occlusal information can be used to investigate the stress distribution with 3D FEA in the supporting bone structure. The LM₁ and the LP²‐LM¹ of a dried modern human skull were scanned by μCT in maximum intercuspation contact. A kinematic analysis of the surface contacts between LM₁ and LP²‐LM¹ during the power stroke was carried out in the occlusal fingerprint analyzer (OFA) software to visualize contact areas during maximum intercuspation contact. This information was used for setting the occlusal molar loading to evaluate the stress distribution in the supporting bone structure using FEA. The output was compared to that obtained when a point force parallel to the long axis of the tooth was loaded in the occlusal basin. For the point load case, our results indicate that the buccal and lingual cortical plates do not experience notable stresses. However, when the occlusal contact areas are considered, the disto‐lingual superior third of the mandible experiences high tensile stresses, while the medio‐lingual cortical bone is subjected to high compressive stresses. Developing a more realistic loading scenario leads to better models to understand the relationship between masticatory function and mandibular shape and structures. Am J Phys Anthropol, 2012. © 2011 Wiley Periodicals, Inc.     

Mapping anisotropy of the proximal femur for enhanced image based finite element analysis

Finite element (FE) models of bone derived from quantitative computed tomography (QCT) rely on realistic material properties to accurately predict bone strength. QCT cannot resolve bone microarchitecture, therefore QCT-based FE models lack the anisotropy apparent within the underlying bone tissue. This study proposes a method for mapping femoral anisotropy using high-resolution peripheral quantitative computed tomography (HR-pQCT) scans of human cadaver specimens. Femur HR-pQCT images were sub-divided into numerous overlapping cubic sub-volumes and the local anisotropy was quantified using a ‘direct-mechanics’ method. The resulting directionality reflected all the major stress lines visible within the trabecular lattice, and provided a realistic estimate of the alignment of Harvesian systems within the cortical compartment. QCT-based FE models of the proximal femur were constructed with isotropic and anisotropic material properties, with directionality interpolated from the map of anisotropy. Models were loaded in a sideways fall configuration and the resulting whole bone stiffness was compared to experimental stiffness and ultimate strength. Anisotropic models were consistently less stiff, but no statistically significant differences in correlation were observed between material models against experimental data. The mean difference in whole bone stiffness between model types was approximately 26%, suggesting that anisotropy can still effect considerable change in the mechanics of proximal femur models. The under prediction of whole bone stiffness in anisotropic models suggests that the orthotropic elastic constants require further investigation. The ability to map mechanical anisotropy from high-resolution images and interpolate information into clinical-resolution models will allow testing of new anisotropic material mapping strategies.