Experimental evidence of impingement induced strains at the interface and the periphery of an embedded acetabular cup implant

After total hip arthroplasty, impingement of implant components may occur during every-day patient activities causing increased shear stresses at the acetabular implant-bone interface. In the literature, impingement related lever-out moments were noted for a number of acetabular components. But there is little information about pelvic load transfer. The aim of the current study was to measure the three-dimensional strain distribution at the macrostructured hemispherical interface and in the periphery of a standard acetabular press-fit cup in an experimental implant-bone substitute model. An experimental setup was developed to simulate impingement loading via a lever arm representing the femoral component and the lower limb. In one experimental setup 12 strain gauges were embedded at predefined positions in the periphery of the acetabular cup implant inside a tray, using polyurethane composite resin as a bone substitute material. By incremental rotation of the implant tray in steps of 10 and 30 deg, respectively, the strains were measured at evenly distributed positions. With the described method 288 genuine strain values were measured in the periphery of an embedded acetabular cup implant in one experimental setup. In two additional setups the strains were evaluated at different distances from the implant interface. Both in radial and meridional interface directions strain magnitudes reach their peak near the rim of the cup below the impingement site. Values of equatorial strains vary near zero and reach their peaks near the rim of the cup on either side and in some distance from the impingement site. Interestingly, the maximum of averaged radial strains does not occur, as expected, close to the interface but at an interface offset of 5.6 mm. With the described experimental setup it is now possible to measure and display the three-dimensional strain distribution in the interface and the periphery of an embedded acetabular cup implant. The current study provides the first experimental proof of the high local stresses gradients in the direct vicinity of the impingement site. The results of the current study help for a better understanding of the impingement mechanism and its impact on acetabular cup stability.     

Finite element analysis of anterior tooth root stresses developed during endodontic treatment

Vertical tooth root fractures are diagnostically challenging, frustrating, and an increasingly common cause of failure of tooth restoration. These vertical root fractures have been associated with many causes, including the endodontic process itself. To investigate these endodontic causes, various phases of crown replacement for an anterior tooth were modeled using nonlinear, plane strain finite element (FE) analysis. Stresses developed during obturation, post positioning, crown placement, and masticatory and occlusal loading of the restored tooth were estimated using this analysis method. The minimum (compressive) principal stress was greatest during obturation of cones 1 and 2, and during mastication of the restored tooth. Although these stresses were significant (-150 to -280 MPa), they did not exceed the compressive strength of dentin. The maximum (tensile) principal stresses, 160 to 300 MPa, were also observed during obturation of cones 1 and 2. These peak stresses exceed the dentin tensile strength.     

Finite element analysis of a hemi-pelvis: the effect of inclusion of cartilage layer on acetabular stresses and strain

An appropriate method of application of the hip-joint force and stress analysis of the pelvic bone, in particular the acetabulum, is necessary to investigate the changes in load transfer due to implantation and to calculate the reference stimulus for bone remodelling simulations. The purpose of the study is to develop a realistic 3D finite element (FE) model of the hemi-pelvis and to assess stress and strain distribution during a gait cycle. The FE modelling approach of the pelvic bone was based on CT scan data and image segmentation of cortical and cancellous bone boundaries. Application of hip-joint force through an anatomical femoral head having a cartilage layer was found to be more appropriate than a perfectly spherical head, thereby leading to more accurate stress–strain distribution in the acetabulum. Within the acetabulum, equivalent strains varied between 0.1% and 0.7% strain in the cancellous bone. High compressive (15–30 MPa) and low tensile (0–5 MPa) stresses were generated within the acetabulum. The hip-joint force is predominantly transferred from the acetabulum through the lateral cortex to the sacroiliac joint and the pubic symphysis. The study is useful to understand the load transfer within the acetabulum and for further investigations on acetabular prosthesis.     

A parametric study of acetabular cup design variables using finite element analysis and statistical design of experiments.

To isolate the primary variables influencing acetabular cup and interface stresses, we performed an evaluation of cup loading and cup support variables, using a Statistical Design of Experiments (SDOE) approach. We developed three-dimensional finite element (FEM) models of the pelvis and adjacent bone. Cup support variables included fixation mechanism (cemented or noncemented), amount of bone support, and presence of metal backing. Cup loading variables included head size and cup thickness, cup/head friction, and conformity between the cup and head. Interaction between and among variables was determined using SDOE techniques. Of the variables tested, conformity, head size, and backing emerged as significant influences on stresses. Since initially nonconforming surfaces would be expected to wear into conforming surfaces, conformity is not expected to be a clinically significant variable. This indicates that head size should be tightly toleranced during manufacturing, and that small changes in head size can have a disproportionate influence on the stress environment. In addition, attention should be paid to the use of nonmetal backed cups, in limiting cup/bone interface stresses. No combination of secondary variables could compensate for, or override the effect of, the primary variables. Based on the results using the SDOE approach, adaptive FEM models simulating the wear process may be able to limit their parameters to head size and cup backing.     

Physical stresses at the air-wall interface of the human nasal cavity during breathing.

The nose is the front line defender of the respiratory system and is rich with mechanoreceptors, thermoreceptors, and nerve endings. A time-dependent computational model of transport through nasal models of a healthy human has been used to analyze the fields of physical stresses that may develop at the air-wall interface of the nasal mucosa. Simulations during quiet breathing revealed wall shear stresses as high as 0.3 Pa in the noselike model and 1.5 Pa in the anatomical model. These values are of the same order of those known to exist in uniform large arteries. The distribution of temperature near the nasal wall at peak inspiration is similar to that of wall shear stresses. The lowest temperatures occur in the vicinity of high stresses due to the narrow passageway in these locations. Time and spatial gradients of these stresses may have functional effects on nasal sensation of airflow and may play a role in the well-being of nasal breathing.     

The influence of friction and interference on the seating of a hemispherical press-fit cup: a finite element investigation.

The formation of gaps in the polar region of acetabular cups is seen as a drawback of press-fit fixation of non-cemented acetabular cups. Recent findings indicate a link between long-term polar gaps and the gaps present directly after implantation. In this study the process of press-fitting is simulated with a linear-elastic two-dimensional axisymmetric finite-element model. The aim of this paper is to investigate the possible importance of friction and interference on the formation of these gaps. A range of cup-bone friction coefficients (mu = 0.1-0.5) is assigned to the cup-bone interface in order to represent the unknown amount of friction occurring during press-fitting. The cup is modeled with a radius of 27 mm, whereas the radius of the cavity is varied between 26.50 and 26.75 mm, thus, creating 0.50 and 0.25 mm radial interference fits. The difference in cavity radius represents the discrepancy between the radius of the last-reamer-used and radius of the cavity it creates. The subchondral plate is considered as being completely removed during reaming. The effects of impact blows via the surgeon's mallet during surgery are modeled as a series of four load pulses, in which peak force is gradually increased from 0.5 to 4.0 kN. The effects of load removal as well as those of load application are investigated. On load application, the cup penetrates into the cavity, and on load removal, the cup rebounds. Depending on the friction, interference and load applied, the position of the cup after the load pulse is somewhere between its position at peak force and its position at the beginning of the pulse. Although the simplifications and conditions involved in the creation of the model necessitate caution when interpreting the results for all clinical cases, it is found that the seating of hemispherical cups in trabecular bone could be more satisfactory for intermediate values of friction (mu = 0.2-0.3) and smaller interference fits (0.25 mm).     

Residual stresses at the stem-cement interface of an idealized cemented hip stem

During the operation of total hip arthroplasty, when the cement polymerizes between the stem implant and the bone, residual stresses are generated in the cement. The purpose of this study was to determine whether including residual stresses at the stem-cement interface of cemented hip implants affected the cement stress distributions due to externally applied loads. An idealized cemented hip implant subjected to bending was numerically investigated for an early post-operative situation. The finite element analysis was three-dimensional and used non-linear contact elements to represent the debonded stem-cement interface. The results showed that the inclusion of the residual stresses at the interface had up to a 4-fold increase in the von Mises cement stresses compared to the case without residual stresses.     

Prosthetic hip failure: retrospective radiograph image analysis of the acetabular cup

A method has been developed for quantifying movement and wear of the acetabular component (cup) of total hip replacements (THR) from routine postoperative and review radiographs. The method uses both interactive and automatic computer image analysis techniques. Dimensions of the prosthesis are used to scale the measurements and so overcom variation in radiographic alignment. The application of the method is illustrated by retrospective investigations of cup migration and wear using review radiographs taken over a follow-up of at least 12 years.     

Finite element analysis of stresses developed in blood sacs of a pusherplate blood pump

Mechanical circulatory support (MCS) devices are blood pumps that support or replace the function of the native heart. It is important to minimize the material stresses in the flexing blood sac or diaphragm in order to increase the duration of support these devices can provide. An axisymmetric finite element model of a pusherplate blood pump was constructed to evaluate the effect of various design parameters on the material stresses in a segmented poly(ether polyurethane urea) seamless blood sac. The design parameters of interest were the sac thickness, pump case wall taper, and radius of the sac between the pusherplate and pump case wall. The analysis involved a quasi-static analysis of the systolic ejection phase of the pump. The finite element solution suggested that the principal stresses and strains increased almost linearly with sac thickness. The pump case wall taper had the largest effect; decreasing the peak principal stresses by approximately 35% when the pump case was straight versus tapered. Lastly, the model demonstrated that the radius of the blood sac between the pusherplate and pump case wall had little or no effect on the magnitude of the blood sac stresses. Therefore, this study suggests that in order to minimize the stresses in a blood sac of a pusherplate blood pump, a straight pump case should be chosen with the thinnest sac.     

Finite element analysis on mechanical state on the osteoclasts under gradient fluid shear stress

Biomechanics and Modeling in Mechanobiology - Mechanical loading, such as fluid shear stress (FSS), is regarded as the main factor that regulates the biological responses of bone cells. Our...     

Finite element evaluation of the AIA shear specimen for bone

An elastic-plastic finite element analysis is performed on the AIA shear specimen to evaluate its effectiveness to yield ultimate shear strength values. The effect of geometry, material properties, and yield criteria are discussed in the light of applications to human femoral cortical bone. Specimen dimensions are noted as follows: W, width, D, hole diameter and H, distance between holes. As the H/D ratio increases the stress distribution tends more toward pure shear at the same time the overshoot in the shear distribution increases. An H/D ratio equal to 1.2-1.5 is optimal. The H/W parameter does not affect the overshoot noticeably but it does slightly affect the purity of shear. The material parameters do affect the performance of the shear specimen. However, the effect of the material parameters are far more pronounced in the anisotropic case than it is in the isotropic case. In the isotropic case, the Young modulus does not affect the overshoot. The increase in Poisson's ratio does slightly decrease the overshoot. For the anisotropic case, the increase in the ratio of shear modulus to Young modulus in the transverse direction (G/E2) results in an increase in the overshoot (in the shear distribution). The increase in the ratio of the Young modulus in the transverse direction to that of the axial direction (E2/E1) also results in an increase in the overshoot. Creating a notch at the top of the hole is shown to have the effect of decreasing the overshoot. Its effect on the purity of the shear is rather slight. It is found that plasticity is initiated at the sides of the two holes where the tensile normal stresses are maximum. The plastic region first expands around the perimeter of the hole then radially outward; and finally, it expands into the significant region. If the W/H parameter is less than 5, a sizable portion of the width of the specimen around the hole can go plastic with the significant region still being in the elastic state. Such a situation can cause tearing of the specimen across the width. A W/H ratio of 6 or more can prevent that danger. It is also found that the onset of plasticity brings about higher overshoot and higher purity of shear. The notched shear specimen performs better in actual tests and is more reliable in producing shear failures. The shear strength results obtained from AIA shear tests tend to confirm those shear strength results obtained from torsion tests.     

Finite element model of stresses in the anterior lens capsule of the eye

Continuum biomechanics continues to aid in the design of clinical procedures and biomedical devices. Such design requires detailed information on the biomechanical properties of the tissues of interest as well as appropriate methods of analysis. In this paper, we use a fully nonlinear, virtual work based finite element method to study the normal stress field in the anterior lens capsule of the porcine eye. The analysis shows that recently measured regional variations in material symmetry may combine with regional variations in membrane thickness to yield a nearly uniform and equibiaxial stress field in normalcy. These findings are discussed in terms of potential implications to the underlying mechanobiology and in particular, in terms of perturbations in stress that result from cataract surgery, the most commonly performed surgical procedure in North America.     

Finite element model of stresses in the anterior lens capsule of the eye

Continuum biomechanics continues to aid in the design of clinical procedures and biomedical devices. Such design requires detailed information on the biomechanical properties of the tissues of interest as well as appropriate methods of analysis. In this paper, we use a fully nonlinear, virtual work based finite element method to study the normal stress field in the anterior lens capsule of the porcine eye. The analysis shows that recently measured regional variations in material symmetry may combine with regional variations in membrane thickness to yield a nearly uniform and equibiaxial stress field in normalcy. These findings are discussed in terms of potential implications to the underlying mechanobiology and in particular, in terms of perturbations in stress that result from cataract surgery, the most commonly performed surgical procedure in North America.     

Finite element stress analysis of a push-out test. Part 1: Fixed interface using stress compatible elements.

In this first part of a two-part paper, interelement stress compatible finite elements are developed and used to perform the stress analysis of a push-out test with a fixed interface. In the formulation, the required continuity of some of the stresses along either a specific interface or all interelement interfaces is enforced by a penalty procedure. The model is axisymmetric and consists of two cylinders attached to each other through the interface. Various relative material properties and boundary conditions are simulated in order to examine their effects on the interface stresses. Both loadings of axial compression force and axial torque are considered. The predicted results exhibit identical interelement stresses and displacements even when highly dissimilar materials are used. They also exhibit a complex state of interface stresses which depend on the geometry, material arrangement, boundary conditions, and loading. The variation of the shear stress is often highly nonuniform and the radial normal stresses are likely to be large. The present results, therefore, disagree with the common assumptions made in the pull-out tests in the orthopaedic applications. Finally, stress analysis of a number of possible testing configurations could lead to the design of an optimal pull-out test which maximizes the usefulness of the measured results in terms of the interface bond strength and factors affecting it.     

Finite element stress analysis of a push-out test. Part II: Free interface with nonlinear friction properties.

In this second part of a two-part paper, nonlinear frictional properties measured at the bone/porous-surfaced metal interface are used to perform the stress analysis of a push-out test assuming free interface. In this case, the friction at the interface is the only mechanism to resist the externally applied load. Similar to the part I, the model is axisymmetric and consists of two cylinders in contact with each other through the interface. Various relative material properties and boundary conditions are simulated in order to examine their effects on the interface stresses and overall push-out resistance. The role of the force-fit and the load direction (push-out versus pull-out) on the results is also investigated. The computed radial and shear stresses are found to markedly vary both with location along the interface and with the testing configuration. The ultimate push-out resistance is also found to significantly alter as the material arrangement and boundary conditions change. The predicted push-out load augments with an increase in the force-fit and diminishes to nil in the absence of a press-fit. For the cases studied here, there is a relative difference of as large as 13 percent between the push-out response and the pull-out response so far as the interface stresses and the maximum resistance are concerned. Therefore, any comparison between the results of push-out (or pull-out) tests performed with different design configurations appears to be invalid.     

Joint contact stresses calculated for acetabular dysplasia patients using discrete element analysis are significantly influenced by the applied gait pattern

Gait modifications in acetabular dysplasia patients may influence cartilage contact stress patterns within the hip joint, with serious implications for clinical outcomes and the risk of developing osteoarthritis. The objective of this study was to understand how the gait pattern used to load computational models of dysplastic hips influences computed joint mechanics. Three-dimensional pre- and post-operative hip models of thirty patients previously treated for hip dysplasia with periacetabular osteotomy (PAO) were developed for performing discrete element analysis (DEA). Using DEA, contact stress patterns were calculated for each pre- and post-operative hip model when loaded with an instrumented total hip, a dysplastic, a matched control, and a normal gait pattern. DEA models loaded with the dysplastic and matched control gait patterns had significantly higher (p = 0.012 and p < 0.001) average pre-operative maximum contact stress than models loaded with the normal gait. Models loaded with the dysplastic and matched control gait patterns had nearly significantly higher (p = 0.051) and significantly higher (p = 0.008) average pre-operative contact stress, respectively, than models loaded with the instrumented hip gait. Following PAO, the average maximum contact stress for DEA models loaded with the dysplastic and matched control patterns decreased, which was significantly different (p < 0.001) from observed increases in maximum contact stress calculated when utilizing the instrumented hip and normal gait patterns. The correlation between change in DEA-computed maximum contact stress and the change in radiographic measurements of lateral center-edge angle were greatest (R² = 0.330) when utilizing the dysplastic gait pattern. These results indicate that utilizing a dysplastic gait pattern to load DEA models may be a crucial element to capturing contact stress patterns most representative of this patient population.     

Finite element analysis of the temporomandibular joint during lateral excursions of the mandible

One of the most significant characteristics of the temporomandibular joint (TMJ) is that it is in fact composed of two joints. Several finite element simulations of the TMJ have been developed but none of them analysed the different responses of its two sides during nonsymmetrical movement. In this paper, a lateral excursion of the mandible was introduced and the biomechanical behaviour of both sides was studied. A three-dimensional finite element model of the joint comprising the bone components, both articular discs, and the temporomandibular ligaments was used. A fibre-reinforced porohyperelastic model was introduced to simulate the behaviour of the articular discs, taking into account the orientation of the fibres in each zone of these cartilage components. The mandible movement during its lateral excursion was introduced as the loading condition in the analysis. As a consequence of the movement asymmetry, the discs were subjected to different load distributions. It was observed that the maximal shear stresses were located in the lateral zone of both discs and that the lateral attachment of the ipsilateral condyle-disc complex suffered a large distortion, due to the compression of this disc against the inferior surface of the temporal bone. These results may be related with possible consequences of a common disorder called bruxism. Although it would be necessary to perform an exhaustive analysis of this disorder, including the contact forces between the teeth during grinding, it could be suggested that a continuous lateral movement of the jaw may lead to perforations of both discs in their lateral part and may damage the lateral attachments of the disc to the condyle.     

Finite element analysis of active Eustachian tube function.

The inability to open the collapsible Eustachian tube (ET) has been related to the development of chronic otitis media. Although ET dysfunction may be due to anatomic and/or mechanical abnormalities, the precise mechanisms by which these structural properties alter ET opening phenomena have not been investigated. Previous investigations could only speculate on how these structural properties influence the tissue deformation processes responsible for ET opening. We have, therefore, developed a computational technique that can quantify these structure-function relationships. Cross-sectional histological images were obtained from eight normal adult human subjects, who had no history of middle ear disease. A midcartilaginous image from each subject was used to create two-dimensional finite element models of the soft tissue structures of the ET. ET opening phenomena were simulated by applying muscle forces on soft tissue surfaces in the appropriate direction and were quantified by calculating the resistance to flow (R(v)) in the opened lumen. A sensitivity analysis was conducted to determine the relative importance of muscle forces and soft-tissue elastic properties. Muscle contraction resulted in a medial-superior rotation of the medial lamina, stretching deformation in the Ostmann's fatty tissue, and lumen dilation. Variability in baseline R(v) values correlated with tissue size, whereas the functional relationship between R(v) and a given mechanical parameter was consistent in all subjects. ET opening was found to be highly sensitive to the applied muscle forces and relatively insensitive to cartilage elastic properties. These computational models have, therefore, identified how different tissue elements alter ET opening phenomena, which elements should be targeted for treatment, and the optimal mechanical properties of these tissue constructs.