Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Finite Element Analysis”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,081 records · Page 60Linked to original sources

Estimation of bolus dispersion effects in perfusion MRI using image-based computational fluid dynamics.

Bolus tracking magnetic resonance imaging (MRI) is a powerful technique for measuring perfusion, and is playing an increasing role in the investigation of acute stroke. However, limitations have been reported when assessing patients with steno-occlusive disease. The presence of a steno-occlusive disease in the artery may cause bolus dispersion, which has been shown to introduce significant errors in cerebral blood flow (CBF) quantification. Bolus dispersion is commonly described by a vascular transport function, but the function that properly characterizes the dispersion is unknown. A novel method to quantify bolus dispersion errors on perfusion measurements is presented. A realistic patient-specific model is constructed from anatomical and physiologic MR data, and the arterial blood flow pattern and the transport of the bolus of contrast agent are computed using finite element analysis. The methodology presented was used also to evaluate the accuracy of three simple vascular models. The methodology was tested on MR data from two normal subjects and two subjects with mild carotid artery stenosis. The estimated CBF errors were of the order of 15% to 20%. However, the presence of stenosis did not necessarily introduce larger dispersion (not only the geometrical model but also the particular physiologic conditions influence the degree of bolus dispersion). The method described will contribute to a better understanding of errors introduced by dispersion effects, to the assessment and validation of vascular models, and to the development of new methods for the correction of dispersion errors in CBF quantification.

Blood Flow Velocity↗

Mechanical influences on tibial fracture healing.

Selected studies are summarized that measure interfragmentary fracture displacements in 6 degrees of freedom at intervals throughout healing in groups of patients with tibial diaphyseal fractures treated by external skeletal fixation. The results are compared with those obtained from experimental studies in which the ideal mechanical conditions for fracture healing were predicted. A finite element analysis model of the healing tibial fracture also was developed. Measured data were used for the analysis, and stress and strain patterns were defined for different stages of healing. Interfragmentary movement measured in the first 6 weeks after injury usually is a magnitude smaller in patients treated by external fixation than in patients treated with cast immobilization. This movement can be much smaller than that predicted to be optimal by experimental studies. A greater amplitude can be achieved, even in stable fractures, by ensuring patient activity. The interfragmentary movement is elastic during loading activity and is generally sinusoidal during steady walking. At the time of dynamization (the unlocking of the frame), a permanent set occurs at the fracture site in all planes. The cyclical movement range in each plane often decreases immediately after unlocking. The model analysis study of fracture healing predicts that tissue damage may occur in the later (hard callus) phase of healing, even while the fixation device is in place, because of abnormally high stresses and strains. This study indicates that fracture mechanics should be controlled more rigorously to provide amplitudes of movement in the first 4 to 6 weeks after fracture. The rigidity of fixation should be increased in the subsequent weeks until the fracture has healed and the frame is removed.

Bony Callus↗

Stress distributions in maxillary bone surrounding overdenture implants with different overdenture attachments.

In this study, effects of different overdenture attachments on the stress distributions in the maxillary bone surrounding the overdenture implants are studied. Four different types of attachment are considered. They are rigid Dalbo Stud, movable Dalbo Stress Broken, movable Dalro, and movable O-ring attachments. Three-dimensional finite element analysis was conducted with commercial package to obtain the stress distributions in the maxillary bone. Varying the attachment types and angle of inclination of load, the stress distributions in the portions of compact bone and trabecular bone were monitored separately. The analysis was conducted by assuming two different boundary conditions at the interface between cap and overdenture abutment in order to evaluate influence of interface boundary condition on stress distribution in the maxillary bone. They were perfect bonding condition and contact with friction at the interfaces. However, it is preferable to assume perfect bond condition at the interface for rigid type attachment systems and contact with friction at the interface for movable type attachment systems. From the numerical results, it was found that the load transfer mechanism of the implant system is altered significantly by the types of the overdenture attachment and also special care must be taken to assign proper boundary conditions at the interface for the analysis. The movable type Dalro attachment generated the highest maximum effective stress in the maxillary bone among the models under the same inclined loading condition for contact with friction. The rigid type Dalbo Stud attachment generated the smallest maximum effective stress in the maxillary bone among the models under the same inclined loading condition for perfect bonding condition.

Dental Cements↗

Altered left ventricular geometry changes the border zone temporal distribution of stress in an experimental model of left ventricular aneurysm: a finite element model study.

BACKGROUND: Left ventricular aneurysm (LVA) is a significant complication of myocardial infarction that may lead to global left ventricular (LV) dysfunction. However, the exact mechanism underlying the abnormal function has not been elucidated. In this study we tested the hypothesis that changes in LV geometry cause both an increase in wall stress and a change in the temporal distribution of stress in the LVA border zone (BZ) during systole. METHODS: Transmural anteroapical infarcts were created in adult Dorsett sheep (n=8) and were allowed to mature into LVAs for 10 weeks. Animals were imaged subsequently using MRI with simultaneous recording of intraventricular pressures. Cardiac models were constructed from the MRI images at end-diastole, isovolumic systole, peak-systole and end-systole. Two short-axis slices, 1 basal and 1 apical were analyzed. The apical slice included the septal and anterior component of the aneurysm as well as the corresponding BZs and normal myocardium. Regional wall stresses were calculated using finite element analysis and compared with stresses in corresponding regions from normal control sheep (n=7). RESULTS: In the LVA group, stress was significantly increased in the BZ at the end-diastolic, isovolumic, peak-systolic, and end-systolic instants (P<0.001 for all). In addition the temporal distribution of stress was significantly altered with maximum stress occurring at peak instead of isovolumic systole. CONCLUSIONS: Geometric changes in the LVA hearts increased wall stress and altered its temporal distribution in the BZ region. Correlation of this finding with the corresponding regional blood flow, oxygen consumption, and mechanical systolic performance may help elucidate the mechanism underlying the observed global LV dysfunction.

Animals↗

[Bone adaptation changes mechanical stress in the femur--a prospective two years follow up after Hüft-TEP implantation].

AIM: Periprosthetic adaptive bone remodeling after total hip arthroplasty (THA) has been frequently simulated in computer models, combining bone remodeling theory with finite element analysis. Unfortunately, there is still a lack of clinical validation data. The collection of prospective volumetric bone density data with a clinical computerized tomography study after THA was previously conducted. The objective of the study presented here is an additional evaluation of femoral strain. METHOD: In a prospective clinical trial with 7 female patients after THA computerized tomography (CT) and finite element (FE) modeling was patient specifically combined immediately after surgery, as well as at three and 24 months postoperatively. Mechanical strain was expressed by strain energy density. RESULTS: Corresponding to a bone density decrease, a decrease of the femoral strain energy density was observed during the two year follow-up after surgery (approximately 20 %). The predominant change occurred within the first three months and was found to be distally stronger than proximally. CONCLUSION: As far as we are aware, this is the first combination of fully prospective 3D CT density data in vivo with patient-specific finite element modeling. The assessment of mechanical strain data during a follow-up trial should be a new approach for analyzing hip stems in clinical biomechanics. A future confirmative study based on a statistical case number calculation would increase the evidence of the results presented here.

Absorptiometry, Photon↗

Evaluation of a testing method for the fatigue performance of total knee tibial trays.

Catastrophic failure of tibial baseplates due to fatigue fracture have occurred in patients. The International Standards Organisation (ISO) have proposed an endurance test to ensure a safety level for trays with respect to this problem. At present the magnitude of the applied load has not been defined. The purpose of this study is to evaluate the ISO test by comparing the results of two tibial trays: one with a clinical history of fracture of the tibial baseplate in some cases (Kinematic) and one without such a record (Kinemax). A load of 2000 N was selected, representing the peak force during the normal walking cycle. Finite element analysis (FEA) predicted that both the Kinematic and Kinemax trays would rapidly fail for two defined contact positions. Mechanical testing performed on both these components confirmed the FEA predictions. However, under service conditions, the Kinematic tray has been reported to have a failure of only 1-2%, the failure occurring at a mean of 54 months for early fracture and 7 years for late fracture. There are no published reports of failed Kinemax trays. It is concluded that the ISO test is an exaggeration of reality when a physiological load is used. As a result of our study we propose that a load of 500 N would effectively screen components with respect to the ISO test conditions for the likelihood of fracture.

Bone Plates↗

On the sensitivity of wall stresses in diseased arteries to variable material properties.

Accurate estimates of stress in an atherosclerotic lesion require knowledge of the material properties of its components (e.g., normal wall, fibrous plaque, calcified regions, lipid pools) that can only be approximated. This leads to considerable uncertainty in these computational predictions. A study was conducted to test the sensitivity of predicted levels of stress and strain to the parameter values of plaque used in finite element analysis. Results show that the stresses within the arterial wall, fibrous plaque, calcified plaque, and lipid pool have low sensitivities for variation in the elastic modulus. Even a +/- 50% variation in elastic modulus leads to less than a 10% change in stress at the site of rupture. Sensitivity to variations in elastic modulus is comparable between isotropic nonlinear, isotropic nonlinear with residual strains, and transversely isotropic linear models. Therefore, stress analysis may be used with confidence that uncertainty in the material properties generates relatively small errors in the prediction of wall stresses. Either isotropic nonlinear or anisotropic linear models provide useful estimates, however the predictions in regions of stress concentration (e.g., the site of rupture) are somewhat more sensitive to the specific model used, increasing by up to 30% from the isotropic nonlinear to orthotropic model in the present example. Changes resulting from the introduction of residual stresses are much smaller.

Anisotropy↗

A natural tooth's stress distribution in occlusion with a dental implant.

The loss of one or more teeth is normally treated with conventional fixed or removable partial dentures or with implant supported fixed or removable dentures. This study investigated stresses formed around the implant and the antagonist natural tooth under occlusal force in the substitution of a missing lower first molar with a rigid or resilient IMZ (Intra Mobil Zylinder) implant, using the finite element stress analysis method. The results indicate that a bite force of 143 N resulted in high compressive stresses around the roots of a natural tooth opposing a restoration supported by an IMZ implant with rigid type abutment. It is speculated that these high compressive stresses may contribute to intrusion of the tooth.

Bite Force↗

Influence of fibril taper on the function of collagen to reinforce extracellular matrix.

Collagen fibrils provide tensile reinforcement for extracellular matrix. In at least some tissues, the fibrils have a paraboloidal taper at their ends. The purpose of this paper is to determine the implications of this taper for the function of collagen fibrils. When a tissue is subjected to low mechanical forces, stress will be transferred to the fibrils elastically. This process was modelled using finite element analysis because there is no analytical theory for elastic stress transfer to a non-cylindrical fibril. When the tissue is subjected to higher mechanical forces, stress will be transferred plastically. This process was modelled analytically. For both elastic and plastic stress transfer, a paraboloidal taper leads to a more uniform distribution of axial tensile stress along the fibril than would be generated if it were cylindrical. The tapered fibril requires half the volume of collagen than a cylindrical fibril of the same length and the stress is shared more evenly along its length. It is also less likely to fracture than a cylindrical fibril of the same length in a tissue subjected to the same mechanical force.

Biomechanical Phenomena↗

A three-dimensional finite element model for arterial clamping.

Clamp induced injuries of the arterial wall may determine the outcome of surgical procedures. Thus, it is important to investigate the underlying mechanical effects. We present a three-dimensional finite element model, which allows the study of the mechanical response of an artery-treated as a two-layer tube-during arterial clamping. The important residual stresses, which are associated with the load-free configuration of the artery, are also considered. In particular, the finite element analysis of the deformation process of a clamped artery and the associated stress distribution is presented. Within the clamping area a zone of axial tensile peak-stresses was identified, which (may) cause intimal and medial injury. This is an additional injury mechanism, which clearly differs from the commonly assumed wall damage occurring due to compression between the jaws of the clamp. The proposed numerical model provides essential insights into the mechanics of the clamping procedure and the associated injury mechanisms. It allows detailed parameter studies on a virtual clamped artery, which can not be performed with other methodologies. This approach has the potential to identify the most appropriate clamps for certain types of arteries and to guide optimal clamp design.

Anisotropy↗

A transversely isotropic viscoelastic constitutive equation for brainstem undergoing finite deformation.

The objective of this study was to define the constitutive response of brainstem undergoing finite shear deformation. Brainstem was characterized as a transversely isotropic viscoelastic material and the material model was formulated for numerical implementation. Model parameters were fit to shear data obtained in porcine brainstem specimens undergoing finite shear deformation in three directions: parallel, perpendicular, and cross sectional to axonal fiber orientation and determined using a combined approach of finite element analysis (FEA) and a genetic algorithm (GA) optimizing method. The average initial shear modulus of brainstem matrix of 4-week old pigs was 12.7 Pa, and therefore the brainstem offers little resistance to large shear deformations in the parallel or perpendicular directions, due to the dominant contribution of the matrix in these directions. The fiber reinforcement stiffness was 121.2 Pa, indicating that brainstem is anisotropic and that axonal fibers have an important role in the cross-sectional direction. The first two leading relative shear relaxation moduli were 0.8973 and 0.0741, respectively, with corresponding characteristic times of 0.0047 and 1.4538 s, respectively, implying rapid relaxation of shear stresses. The developed material model and parameter estimation technique are likely to find broad applications in neural and orthopaedic tissues.

Animals↗

A 6-year prospective study of 620 stress-diversion surface (SDS) dental implants.

A 6-year prospective study was conducted to assess the clinical success rates and crestal bone response of a dental implant system with a stress-diversion design. Mathematical modeling, digital radiography with applied isodensity, and finite element analysis were used to highlight the effect of the stress distribution design. A total of 386 hydroxyapatite-coated prototypes and 234 commercial grit-blasted external hex implants were placed in virgin bone as well as various grafted maxillary regions, with 36% of the posterior implants being immediately loaded. Prototypes achieved 96.6% survival over a 3-year period. The grit-blasted implant, placed from 2000 to 2003, showed a 95% survival rate. There were no significant changes in crestal bone levels after the first 12 months of prosthetic loading. Engineering evaluations suggested that undesirable stresses were distributed from the crest of the ridge down through the center of the implant body.

Air Abrasion, Dental↗

A comparative study of most suitable miniplate fixation for mandibular symphysis fracture using a finite element model.

The purpose of this study is to determine the most stable fixation method for mandibular symphysis fractures by comparing the mechanical characteristics of models fixed at different positions with different numbers of plates. Fractures were generated in 3-dimensional finite element models, and were fixed with a single miniplate, parallel double miniplates, or perpendicular double miniplates. A 300 N perpendicular load was then applied on the left molar region, and a finite element analysis was performed. We compared vertical gaps between the fractured surfaces, maximum stress within the screw/plating system, and maximum stress around screw holes in the bone. Compared to the single miniplate, both the parallel and perpendicular double miniplates demonstrated significantly less stress in the screw/plating system and screw holes in the bone. In addition, the perpendicular double miniplateshad significantly smaller vertical gaps between fracture surfaces when compared to the single miniplate. Comparing parallel and perpendicular double miniplate fixations, less stress was found around the screw holes of the perpendicular miniplate models than those of the parallel miniplatemodels. There were no differences in vertical gaps or maximum stress within the screw/plating systems between the 2 double miniplate fixations. These results suggest that perpendicular double miniplate fixation is more suitable for fixing mandibular symphysis fractures.

Bone Plates↗

Validation data for periprosthetic bone remodelling theories.

Periprosthetic adaptive bone remodelling after total hip arthroplasty (THA) has been frequently simulated in computer models, combining bone remodelling theory with finite element analysis. Unfortunately, there still subsist a lack of clinical data, which are necessary for validation of these simulation results. Therefore, the objective of the current project is to collect prospective volumetric bone density data with a clinical computerized tomography study in seven patients after THA. A retrospective study 12 years after implantation in 11 patients was added. A data set of about 100000 bone voxels for each femur was collected. In all prospective cases, the predominant change is seen during the first year. The average density reduction in the horizontal slices was between 50 and 150 Hounsfield units (HU) (approx. 10%; p<0.001) after 2 years. Loss of density is particularly strong distal of the minor trochanter and decreases from proximal to distal. For the 12 years retrospective study, the contralateral femur provided the control. Similar trends comparable to the prospective 2-year follow-up CT density values were seen in most cases with density reductions of up to 400 HU (30%). However, in one of these cases there was no difference between the operated and the control density. As far as we are aware, this is the first collection of fully prospective 3D validation data in vivo for periprosthetic adaptive bone remodelling theories. The data are also unique as they are suitable for direct patient-specific 3D finite element meshing and individual weight-related loading.

Absorptiometry, Photon↗

Wall shear stress gradient topography in the normal left coronary arterial tree: possible implications for atherogenesis.

OBJECTIVE: Wall shear stress gradient (WSSG) in vitro has shown its importance in atherogenesis, probably as a local modulator of endothelial gene expression. The purpose of this study is to numerically analyse the WSSG distribution over the normal human left coronary artery (LCA) tree. RESEARCH DESIGN AND METHODS: A three-dimensional computer generated model of the LCA tree, based on an averaged human data set extracted from angiographies, was adopted for finite-element analysis. The LCA tree includes the left main coronary artery (LMCA), the left anterior descending (LAD), the left circumflex artery (LCxA) and their major branches. RESULTS: In proximal LCA tree regions where at bifurcations in regions opposite the flow atherosclerosis frequently occurs, low WSSG appears. At distal segments, the WSSG increases substantially due to increased velocity resulting from increased vessel tapering. Low WSSG occurs dividers, which are anatomic sites predisposed for atherosclerotic development. CONCLUSIONS: This computational work determines, probably for the first time, the topography of the WSSG in the normal human LCA tree. Spatial WSSG differentiation indicates that low values of this parameter probably correlate to atherosclerosis localization. However, further studies are needed to clarify the role of WSSG in atherogenesis.

Blood Flow Velocity↗

[Effect of implant geometry on strain distribution in peri-implant bone].

Recent developments and results in the field of mechanical bone stimulation have made it possible to define limits between physiological and non-physiological transmission of strength to the bone. Now it is necessary to investigate the different geometries of dental implants regarding biomechanical reliability. The aim of this study was to examine the distribution of strain on the bone along different basic forms of implants during loading using finite element analysis (FEA). The following implant designs were included in the study: cylinder, threaded cylinder, cylinder with steps, threaded cylinder with steps, and double-disk implants. All implants had a length of 12 mm and a diameter of 4 mm, the axial loading was defined as 300 N. Threaded implants showed more homogeneous distribution of strain than stepped implants, however, the maximum values were below the physiological range. Shortening of the implants led to an increase of the values, however, there was a zone of hypophysiological strain apically. Disk implants showed extremely low values at the disk margins, possibly due to the lack of physiological bone stimulus by the disks. Overall, none of the implants showed optimal distribution of strain, even though homogeneous strain distribution is decisive for long-term implant stability.

Alveolar Process↗

The effect of striking angle on the buckling mechanism in blowout fracture.

BACKGROUND: The buckling mechanism is widely accepted as a mechanism of blowout fractures, along with the hydraulic mechanism. Although many studies have been performed related to the buckling mechanism, none of them have taken the direction of the striking force into consideration. As the orbital floor is not parallel to the horizontal plane, a difference in the striking force direction might affect resultant fracture patterns. The present study aims to investigate whether fracture patterns in the orbital floor were influenced by the striking force direction in terms of the buckling mechanism. METHODS: The authors produced three-dimensional models on a workstation simulating eight dry skulls and applied striking forces on the orbital rim of each model from three different angles (0, 15, and 30 degrees in the upward direction). Using finite element analysis, the authors calculated the width of the area where the resultant stresses exceed the bone's yielding criterion. The width was termed the "theoretical fracture width" because, theoretically, fracture was expected to occur in the area. Then, the authors compared the theoretical fracture width in groups with the three different striking force angles. Finally, the validity of the theoretical width was verified with an experiment on actual skull models. RESULTS: The theoretical fracture width was the greatest when the striking force was directed at 30 degrees in the upward direction. CONCLUSIONS: For the buckling mechanism, fracture would occur in a wider area of the orbital floor when striking force was directed upward than when the force was horizontally directed. This finding would be helpful in predicting fracture width in blowout fractures.

Biomechanical Phenomena↗

Ultrasonic backscattering from porcine whole blood of varying hematocrit and shear rate under pulsatile flow.

It was shown previously that ultrasonic scattering from whole blood varies during a flow cycle under pulsatile flow both in vitro and in vivo. It has been postulated that this cyclic variation may be associated with the dynamics of red cell aggregation because the shearing force acting on the red cell aggregates across the lumen is a function of time during a flow cycle. In all studies, the local shear rate variation as a function of time is unknown. The effect of shear rate on the red cell aggregation and, thus, on ultrasonic scattering from blood can only be merely speculated. One solution to this problem is to estimate the shear rate in a flow conduit by finite element analysis (FEA). An FEA computational fluid dynamics (CFD) tool was used to calculate local shear rate in a series of experiments in which ultrasonic backscattering from porcine whole blood under pulsatile flow was measured as a function of hematocrit and shear rate intravascularly with a 10-MHz catheter-mounted transducer in a mock flow loop. The results show that, at 20 beats per min (BPM), the magnitudes of the cyclic variation for hematocrits at 30, 40, and 50% were approximately 4 dB. However, at 60 BPM, the magnitude of cyclic variation was found to be minimal. The results also confirm previous findings that the amplitude and the timing of the peak of ultrasonic backscattering from porcine whole blood under pulsatile flow during a flow cycle are dependent upon the shear rate and hematocrit in a complicated way.

Animals↗