Stress changes in the femoral head due to porous ingrowth surface replacement arthroplasty

Finite element stress analyses were conducted of the canine femoral head before and after implantation of various surface replacement-type components. The femoral head was replaced by four implant geometries; (a) shell, (b) shell with peg, (c) shell with rod, and (d) a new epiphyseal replacement design. All implants were modelled to simulate bony ingrowth along the underside of the shell and along the surfaces of the peg and rod. The results indicated that in the normal femur the forces are transferred from the articular surface through the femoral head cancellous bone to the inferior cortical shell of the femoral neck. After shell-type surface replacement, forces were transferred more distally at the rim of the shell and at the end of the peg or rod, thereby reducing the stresses in the proximal head cancellous bone. Computer simulation of bone remodelling due to proximal bone stress reduction was shown to accentuate the abnormality of the stress fields. Surface replacement with a lower modulus material created a less abnormal redistribution of bone stresses. The new epiphyseal replacement design resulted in stress distributions similar to those in the normal femoral head and minimal shear stresses at the implant/bone interface. These findings suggest that the epiphyseal replacement concept may provide better initial mechanical integrity and create a more benign milieu for adaptive bone remodelling than conventional, shell-type surface replacement components.     

Stress amplifications in dental non-carious cervical lesions

This study aims to investigate the influence of the presence, shape and depth of NCCLs on the mechanical response of a maxillary second premolar subjected to functional and non-functional occlusal loadings using 3-D finite element (FE) analysis. A three-dimensional model of a maxillary second premolar and its supporting bone was constructed based on the contours of their cross-sections. From the sound model, cervical defects having either V- or U-shapes, as found clinically, were subtracted in three different depths. The models were loaded with 105 N to simulate normal chewing forces according to a functional occlusal loading (F1) vertically applied and two non-functional loadings (F2 and F3) obliquely oriented. Two alveolar bone crest heights were tested. Ansys™ FE software was used to compute stress distributions and maximum principal stress for each of the models. The presence of a lesion had no effect on the overall stress distribution of the system, but affected local stress concentrations. Non-functional loadings exhibited tensile stresses concentrating at the cervical areas and root surfaces, while the functional loading resulted in homogeneous stress distributions within the tooth. V-shaped lesions showed higher stress levels concentrated at the zenith of the lesion, whereas in U-shaped defect stresses concentrated over a wider area. As the lesions advanced in depth, the stress was amplified at their deepest part. A trend of stress amplification was observed with decreasing bone height. These results suggest a non-linear lesion progression with time, with the progression rate increasing with patient's age (deeper lesions and lower bone support).     

Stress analysis of the proximo-medial femur after total hip replacement

Stress-induced bone loss in the proximo-medial femur has been identified as a factor leading to loosening in the artificial hip joint. In an effort to develop a quantitative understanding of the stress distribution that causes bone loss, axial and hoop stresses in the medial calcar of the femur have been determined after total hip replacement, using finite element stress analysis. Stress distributions for a high and a low Young's modulus prosthesis material are compared for both collared and uncollared prosthesis designs. The use of a low-modulus material, and of a collar, are predicted to be advantageous, giving rise to proximo-medial stress patterns similar to those of the normal, intact femur.     

Evaluation of bone plate with low-stiffness material in terms of stress distribution

The aim of this study is to evaluate a newly developed bone plate with low-stiffness material in terms of stress distribution. In this numerical study, 3D finite element models of the bone plate with low-stiffness material and traditional bone plates made of stainless steel and Ti alloy have been developed by using the ANSYS software. Stress analyses have been carried out for all three models under the same loading and boundary conditions. Compressive stresses occurring in the intact portion of the bone (tibia) and at the fractured interface at different stages of bone healing have been investigated for all three types of bone-plate systems. The results obtained have been compared and presented in graphs. It has been seen that the bone plate with low-stiffness material offers less stress-shielding to the bone, providing a higher compressive stress at the fractured interface to induce accelerated healing in comparison with Ti alloy and stainless-steel bone plate. In addition, the effects of low-stiffness materials with different Young's modulus on stress distribution at the fractured interface have been investigated in the newly developed bone-plate system. The results showed that when a certain value of Young's modulus of low-stiffness material is exceeded, increase in stiffness of the bone plate does not occur to a large extent and stress distributions and micro-motions at the fractured interface do not change considerably.     

Free body analysis,beam mechanics,and finite element modeling of the mandible of Alligator mississippiensis

The mechanical behavior of mammalian mandibles is well‐studied, but a comprehensive biomechanical analysis (incorporating detailed muscle architecture, accurate material properties, and three‐dimensional mechanical behavior) of an extant archosaur mandible has never been carried out. This makes it unclear how closely models of extant and extinct archosaur mandibles reflect reality and prevents comparisons of structure–function relationships in mammalian and archosaur mandibles. We tested hypotheses regarding the mechanical behavior of the mandible of Alligator mississippiensis by analyzing reaction forces and bending, shear, and torsional stress regimes in six models of varying complexity. Models included free body analysis using basic lever arm mechanics, 2D and 3D beam models, and three high‐resolution finite element models of the Alligator mandible, incorporating, respectively, isotropic bone without sutures, anisotropic bone with sutures, and anisotropic bone with sutures and contact between the mandible and the pterygoid flange. Compared with the beam models, the Alligator finite element models exhibited less spatial variability in dorsoventral bending and sagittal shear stress, as well as lower peak values for these stresses, suggesting that Alligator mandibular morphology is in part designed to reduce these stresses during biting. However, the Alligator models exhibited greater variability in the distribution of mediolateral and torsional stresses than the beam models. Incorporating anisotropic bone material properties and sutures into the model reduced dorsoventral and torsional stresses within the mandible, but led to elevated mediolateral stresses. These mediolateral stresses were mitigated by the addition of a pterygoid‐mandibular contact, suggesting important contributions from, and trade‐offs between, material properties and external constraints in Alligator mandible design. Our results suggest that beam modeling does not accurately represent the mechanical behavior of the Alligator mandible, including important performance metrics such as magnitude and orientation of reaction forces, and mediolateral bending and torsional stress distributions. J.Morphol. 2011. © 2011 Wiley‐Liss, Inc.     

Collagen fiber orientation at the tendon to bone insertion and its influence on stress concentrations

The tendon to bone insertion serves the mechanical role of transferring loads from a relatively compliant tendon to a relatively rigid bone. The details of the mechanism of load transfer are of great importance, since current surgical procedures for tendon reattachment have high failure rates. We hypothesized that the microscopic structure of the insertion is optimized to minimize stress concentrations associated with this load transfer. To explore this, collagen fiber orientation distributions were measured in the supraspinatus tendons of rats. The angular deviation of fibers was fairly uniform across the insertion, and the mean angles of the local distributions deviated mildly from the tendon axis. To explore how these observed property distributions could influence load transfer, these distributions were used to derive material properties for an idealized two-dimensional mechanical model of an insertion. Comparison between stress concentrations in this idealized model and those in three comparison models suggests that the microstructure serves to (1) simultaneously reduce stress concentrations and material mass, and (2) shield the insertion's outward splay from the highest stresses.     

The influence of mechanical stimulus on the pattern of tissue differentiation in a long bone fracture--an FEM study

2D, coronal plane, finite elements models (FEMs) were developed from orthogonal radiographs of a diaphyseal tibial fracture and its reparative tissue at four different time points during healing. Each callus was separated into regions of common tissue histology by computerised radiographic analysis. Starting point values of tissue material properties from the literature were refined by the model to simulate exactly the mechanical behaviour of the subject's callus and bone during loading. This was achieved by matching measured inter-fragmentary displacements with calculated inter-fragmentary forces. Stress and strain distributions in the callus and bone were calculated from peak inter-fragmentary displacements measured during natural walking activity, and were correlated with the subsequently observed pattern of tissue differentiation and maturation of the callus. The growth and stiffening of the external callus progressively reduced the inter-fragmentary gap strain. Partial maturation of the gap tissue was apparent only one week before fixator removal. Principal stresses in the callus were compared with 'yield stresses' in corresponding tissue from the literature. This indicated the presence of stress concentrations medial and lateral to the fracture gap, which probably caused tissue damage during normal activity levels. Tissue damage may also have precipitated partial structural failure of the callus, both of which were believed to have delayed healing during the middle third of the fixation period. Had the fixation device provided greater inter-fragmentary support during early healing, this may have prevented callus failure and the consequent delay in healing. A further benefit of this would have been the reduction of the initially high intra-gap tissue strains to a magnitude more conducive to earlier maturation of the bridging tissue that united the bone.     

Method for the design and analysis of a non-linear anisotropic lower limb prosthetic socket

A simplified computer-based structural analysis procedure has been developed for evaluating the stresses in a lower limb prosthetic socket when subjected to external loads. This technique has been tested by studying a simple linear axisymmetric shell subjected to uniform loads. It was found that the stress and strain distributions obtained are close to those determined by a finite element technique. The method that has been proposed involves an incremental procedure, which can be used to analyse the behaviour of a prosthetic socket of irregular shape made from non-linear anisotropic materials. By employing the proposed procedure to study an example, in which a simplified lower limb prosthetic socket is subject to internal pressure and frictional forces, the effects of the material non-linearity have been found to be significant.     

Morphological and numerical studies of the stress compatibility of patella implants]

Failure of total knee arthroplasty is relatively often caused by problems of the patellofemoral replacement. The purpose of this study was to analyze the distribution of stresses within an anatomical patella and the changes in stress distribution after patellar resurfacing with a Miller-Galante I patellar implant using two- and three dimensional finite element models (FEM). To assess validity, FEM results were compared with morphological findings from contact radiographs and densitographs. Internal orientation of bone trabeculae is in good agreement with the orientation of theoretically calculated principal stresses. Almost unchanged principal tensile stresses after implantation, together with the lack of extreme stress peaks within the cancellous bone ensure stress compatibility of the implant. In the case of a firmly seated implant with good bone ingrowth, increased von Mises stresses are found near the fixation peg/plate junction. Their relevance for improved bone ingrowth near this part of the interface is emphasized. At the same time, material failure at the peg/plate junction can be better understood. An analysis of the early postoperative period assuming nonlinear interface conditions failed to demonstrate an uniform distribution of normal and tangential interface forces.     

Finite element analysis of microelectrotension of cell membranes

Electric fields can be focused by micropipette-based electrodes to induce stresses on cell membranes leading to tension and poration. To date, however, these membrane stress distributions have not been quantified. In this study, we determine membrane tension, stress, and strain distributions in the vicinity of a microelectrode using finite element analysis of a multiscale electro-mechanical model of pipette, media, membrane, actin cortex, and cytoplasm. Electric field forces are coupled to membranes using the Maxwell stress tensor and membrane electrocompression theory. Results suggest that micropipette electrodes provide a new non-contact method to deliver physiological stresses directly to membranes in a focused and controlled manner, thus providing the quantitative foundation for micreoelectrotension, a new technique for membrane mechanobiology.     

Virtualisation of stress distribution in heart valve tissue

This study presents an image-based finite element analysis incorporating histological photomicrographs of heart valve tissues. We report stress fields inside heart valve tissues, where heterogeneously distributed collagen fibres are responsible for transmitting forces into cells. Linear isotropic and anisotropic tissue material property models are incorporated to quantify the overall stress distributions in heart valve tissues. By establishing an effective predictive method with new computational tools and by performing virtual experiments on the heart valve tissue photomicrographs, we clarify how stresses are transferred from matrix to cell. The results clearly reveal the roles of heterogeneously distributed collagen fibres in mitigating stress developments inside heart valve tissues. Moreover, most local peak stresses occur around cell nuclei, suggesting that higher stress may be mediated by cells for biomechanical regulations.     

Measurement of transient and residual stresses during polymerization of bone cement for cemented hip implants

The initial fixation of a cemented hip implant relies on the strength of the interface between the stem, bone cement and adjacent bone. Bone cement is used as grouting material to fix the prosthesis to the bone. The curing process of bone cement is an exothermic reaction where bone cement undergoes volumetric changes that will generate transient stresses resulting in residual stresses once polymerization is completed. However, the precise magnitude of these stresses is still not well documented in the literature. The objective of this study is to develop an experiment for the direct measurement of the transient and residual radial stresses at the stem-cement interface generated during cement polymerization. The idealized femoral-cemented implant consists of a stem placed inside a hollow cylindrical bone filled with bone cement. A sub-miniature load cell is inserted inside the stem to make a direct measurement of the radial compressive forces at the stem-cement interface, which are then converted to radial stresses. A thermocouple measures the temperature evolution during the polymerization process. The results show the evolution of stress generation corresponding to volumetric changes in the cement. The effect of initial temperature of the stem and bone as well as the cement-bone interface condition (adhesion or no adhesion) on residual radial stresses is investigated. A maximum peak temperature of 70 degrees C corresponds to a peak in transient stress during cement curing. Maximum radial residual stresses of 0.6MPa in compression are measured for the preheated stem.     

An in vitro study of implant--tooth-supported connections using a robot test system.

Many unsolved problems in dental implant research concern the interfacial stress distributions between the implant components, as well as between the implant surface and contacting bone. To obtain a mechanical understanding of how vertical and horizontal occlusal forces are distributed in this context, it is crucial to develop in vitro testing systems to measure the force transmission between dental implants and attached prostheses. A new approach to such testing, involving a robotic system, is described in this investigation. The system has been designed to produce simulated mandibular movements and occlusal contact forces so that various implant designs and procedures can be thoroughly tested and evaluated before animal testing or human clinical trials. Two commonly used fixed prosthesis designs used to connect an implant and a tooth, a rigid connection and a nonrigid connection, were fabricated and used for experimental verification. The displacement and force distributions generated during simulated chewing activities were measured in vitro. Force levels, potentially harmful to human bone surrounding the connected dental implant and tooth, were analyzed. These results are useful in the design of prostheses and connecting components that will reduce failures and limit stress transfer to the implant/bone interface.     

2D-finite element analyses and histomorphology of lag screws with and without a biconcave washer.

For osteosynthesis and for bone transplant fixation in particular, a lag screw with a biconcave washer, the so called "Anchor Screw" (AS) has been introduced in maxillo-facial surgery. Using 2D-finite element analysis (FEA), the v.Mises and the circumferential stresses induced in underlying bone by this AS are analysed and compared to those under a conventional lag screw. The stress distributions below the biconcave washer of the AS were correlated with histomorphological bone reactions after AS osteosynthesis in two tumor patients, retrieved 12 weeks and 19 months after tumor surgery, respectively. Depending on the thickness of cortical bone, the v.Mises stress concentrations below the biconcave washer were lower than under the head of the conventional lag screw (CLS), but with a higher stress maximum concentrated around the rim of the washer. The circumferential stresses were only half as high around the AS, and thus the deformation of bone was reduced. As predicted by FEA, histology showed microcrack formation, but then after minimal resorption, remodelling of bone below the biconcave washer. Stable osteosynthesis could be demonstrated by bony union already after 12 weeks, and, while bone remodelling continued in the healed osteotomy, it had decreased around the screws after 19 months. It can be concluded from the biomechanical principles and the histomorphological findings that the AS appears superior to the CLS.     

An approach for the stress analysis of transversely isotropic biphasic cartilage under impact load.

Stress analysis of contact models for isotropic articular cartilage under impacting loads shows high shear stresses at the interface with the subchondral bone and normal compressive stresses near the surface of the cartilage. These stress distributions are not consistent, with lesions observed on the cartilage surface of rabbit patellae from blunt impact, for example, to the patello-femoral joint. The purpose of the present study was to analyze, using the elastic capabilities of a finite element code, the stress distribution in more morphologically realistic transversely isotropic biphasic contact models of cartilage. The elastic properties of an incompressible material, equivalent to those of the transversely isotropic biphasic material at time zero, were derived algebraically using stress-strain relations. Results of the stress analysis showed the highest shear stresses on the surface of the solid skeleton of the cartilage and tensile stresses in the zone of contact. These results can help explain the mechanisms responsible for surface injuries observed during blunt insult experiments.     

Computational modeling of flow-induced shear stresses within 3D salt-leached porous scaffolds imaged via micro-CT

Flow-induced shear stresses have been found to be a stimulatory factor in pre-osteoblastic cells seeded in 3D porous scaffolds and cultured under continuous flow perfusion. However, due to the complex internal structure of porous scaffolds, analytical estimation of the local shear forces is impractical. The primary goal of this work is to investigate the shear stress distributions within Poly(l-lactic acid) scaffolds via computation. Scaffolds used in this study are prepared via salt leeching with various geometric characteristics (80–95% porosity and 215–402.5 μm average pore size). High resolution micro-computed tomography is used to obtain their 3D structure. Flow of osteogenic media through the scaffolds is modeled via lattice Boltzmann method. It is found that the surface stress distributions within the scaffolds are characterized by long tails to the right (a positive skewness). Their shape is not strongly dependent on the scaffold manufacturing parameters, but the magnitudes of the stresses are. Correlations are prepared for the estimation of the average surface shear stress experienced by the cells within the scaffolds and of the probability density function of the surface stresses. Though the manufacturing technique does not appear to affect the shape of the shear stress distributions, presence of manufacturing defects is found to be significant: defects create areas of high flow and high stress along their periphery. The results of this study are applicable to other polymer systems provided that they are manufactured by a similar salt leeching technique, while the imaging/modeling approach is applicable to all scaffolds relevant to tissue engineering.     

Epiphyseal-based designs for tibial plateau components--II. Stress analysis in the sagittal plane

A two-dimensional, finite element study was undertaken to establish the stresses in the proximal tibia before and after total knee arthroplasty. Equivalent-thickness models in a sagittal plane were created for the natural, proximal tibia and for the proximal tibia with two different types of tibial plateau components. All components simulated bony ingrowth fixation, i.e. no cement layer existed between component and bone. In addition, the interface between component and bone was assumed to be intimately connected, representing complete bony ingrowth and a rigid state of fixation. Two load cases were considered: a joint reaction force acting in conjunction with a patellar ligament force, simulating the knee at 40 degrees of flexion; and a joint reaction force directed along the long axis of the tibia. For the natural tibia model, the pattern of principal stresses for loadcase 1 more closely corresponds to the epiphyseal plate geometry and trabecular morphology than do the principal stress patterns for loadcase 2. Judging from the distribution of principal stresses, loadcase 1 represents a more severe test of implant design than does loadcase 2. The model of the component with a peg predicted that the trabecular bone near the tip of the peg will experience higher than normal stresses, while the bone stresses near the posterior aspect adjacent to the metal tray will be reduced. A component without pegs that incorporates a posterior chamfer and an anterior lip lead to stress distributions closer to those existing in the natural tibia. The interface geometry for this design is based upon the contour of the epiphyseal plate.     

Cortical bone thickness can adapt locally to muscular loading while changing with age

Mechanical loading of muscle action is concentrated at muscle attachment sites; thus there may be a potential for site-specific variation in cortical bone thickness. Humeri from an early 20th-century Finnish (Helsinki) and two medieval English (Newcastle, Blackgate and York, Barbican) populations were subjected to pQCT scanning to calculate site-specific cross-sectional cortical bone area (CA) for four locations and to measure cortical thickness at muscle attachment sites and non-attachment sites. We found that CA at 80% of humerus length was significantly reduced compared to more distal cross-sections, which can be due to reduced stresses at the proximal shaft. The principal direction of loading at 80% humerus length was towards mediolateral plane, likely due to fixing the humerus close to the torso. At 35% the main direction of loading was towards anteroposterior plane, reflecting elbow flexing forces. The principal direction of loading varied between populations, sides and sexes at 50% humerus length due to preference between elbow and shoulder joint; thus this location might be useful when trying to infer differences in activity. These changes are likely due to overall shaft adaptation to forces acting at the humerus. In addition, we found a potential for site-specific variation in cortical thickness; cortical bone at muscle attachment sites was significantly thicker compared to non-attachment sites. Lastly, CA at 35% of humerus length and cortical thickness at non-attachment sites decreased with age. These results underline the importance of muscle loading for bone mass preservation as well as indicate that a site-specific variation of bone mass is possible.     

Numerical study of tension and strain distribution around rat molars]

A knowledge of the mechanical processes triggered in the bone and periodontal ligament (PDL) by orthodontic forces applied to a tooth is of decisive importance for an understanding of the subsequent remodelling around the tooth. To investigate these mechanical relationships, three-dimensional finite element (FE) models of the first lower molar in the rat were established. On the basis of digitized serial histological sections, these FE models were generated semi-automatically. Using various simplified geometrical variations, an appropriate FE model for the analysis of the stress and strain distributions was established. The numerical analyses were carried out under a mesially directed force of 0.1 N. Stress distributions in the bone and PDL showed a similar pattern, while strains in the bone were lower than in the PDL by a factor of 10-5. The data confirm the assumption that strain patterns in the PDL may be the key stimulus of bone remodelling.     

Tissue stresses and strain in trabeculae of a canine proximal femur can be quantified from computer reconstructions

A quantitative assessment of bone tissue stresses and strains is essential for the understanding of failure mechanisms associated with osteoporosis, osteoarthritis, loosening of implants and cell- mediated adaptive bone-remodeling processes. According to Wolff's trajectorial hypothesis, the trabecular architecture is such that minimal tissue stresses are paired with minimal weight. This paradigm at least suggests that, normally, stresses and strains should be distributed rather evenly over the trabecular architecture. Although bone stresses at the apparent level were determined with finite element analysis (FEA), by assuming it to be continuous, there is no data available on trabecular tissue stresses or strains of bones in situ under physiological loading conditions. The objectives of this project were to supply reasonable estimates of these quantities for the canine femur, to compare trabecular-tissue to apparent stresses, and to test Wolff's hypothesis in a quantitative sense. For that purpose, the newly developed method of large-scale micro-FEA was applied in conjunction with micro-CT structural measurements. A three-dimensional high-resolution computer reconstruction of a proximal canine femur was made using a micro-CT scanner. This was converted to a large-scale FE-model with 7.6 million elements, adequately refined to represent individual trabeculae. Using a special-purpose FE-solver, analyses were conducted for three different orthogonal hip-joint loading cases, one of which represented the stance-phase of walking. By superimposing the results, the tissue stress and strain distributions could also be calculated for other force directions. Further analyses of results were concentrated on a trabecular volume of interest (VOI) located in the center of the head. For the stance phase of walking an average tissue principal strain in the VOI of 279 strain was found, with a standard deviation of 212 microstrain. The standard deviation depended not only on the hip-force magnitude, but also on its direction. In more than 95% of the tissue volume the principal stresses and strains were in a range from zero to three times the averages, for all hip-force directions. This indicates that no single load creates even stress or strain distributions in the trabecular architecture. Nevertheless, excessive values occurred at few locations only, and the maximum tissue stress was approximately half the value reported for the tissue fatigue strength. These results thus indicate that trabecular bone tissue has a safety factor of approximately two for hip-joint loads that occur during normal activities.