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,045 records · Page 58Linked to original sources

Effects of risedronate on trabecular microstructure and biomechanical properties in ovariectomized rat tibia.

We determined the effect of risedronate on the trabecular microstructure of ovariectomized rat tibiae, using micro-computed tomography, in order to investigate how changes in microstructure contribute to biomechanical properties. Fifty 18-week-old rats underwent sham operation (n=10) or ovariectomy (OVX) (n=40). The OVX rats were further divided into four groups (n=10 for each group) and treated with risedronate at doses of 0, 0.1, 0.5 or 2.5 mg/kg for 9 months. OVX caused deterioration of three-dimensional trabecular microstructure, notably structure model index (SMI) and connectivity density, while treatment of OVX rats with risedronate at 0.5 and 2.5 mg/kg improved those deleterious microstructural changes. Biomechanical property, as assessed by finite element analysis (FEA), correlated significantly with trabecular bone volume fraction (BV/TV), and the correlation further increased substantially when microstructural parameters were added, especially SMI and connectivity density, with risedronate therapy. Thus, it is suggested that, in addition to increasing bone mass, risedronate improves biomechanical property by maintaining a plate-like structure as well as connectivity of trabeculae.

Animals↗

An advanced computer-aided geometric modeling and fabrication method for human middle ear.

This paper presents a practical and systematic method for reconstructing accurate computer and physical models of the entire human middle ear. The proposed method starts with the histological section preparation of human temporal bone. Through tracing outlines of the middle ear components on the sections, a set of discrete points is obtained and employed to construct B-spline curves that represent the exterior contours of the components using a curve-fitting technique. The surface-skinning technique is then employed to quilt the B-spline curves for smooth boundary surfaces of the middle ear components using B-spline surfaces. The solid models of the middle ear components are constructed using these surfaces and then assembled to create the entire middle ear in a computer-aided design environment. This method not only provides an effective way to visualize and measure the three-dimensional structure of the middle ear, but also provides a detailed knowledge of middle ear geometry that is required for finite element analysis or multibody dynamic analysis of the human middle ear. In addition, the geometric model constructed using the proposed method is smooth and can be fabricated in various scales using solid freeform fabrication technology. The physical model of the human middle ear is extremely effective in realizing the middle ear anatomy and enhancing discussion and collaboration among researchers and physicians.

Anatomy, Cross-Sectional↗

Effect of material property representation on stresses in endoprostheses.

The present study sought to examine if, and to what extent, the representation of the elastic properties of material(s) in a model of an endoprosthesis per se or a construct that includes an endoprosthesis affects the stresses in the model, when the finite element analysis method is used to obtain these stresses. For this purpose, two models were examined. One was of an elastomeric metacarpophalangeal joint implant, with the material being assumed to have isotropic elastic properties in one case while it was represented as hyperelastic in the other. The model used in the second study case was that of a construct comprising the tibial component of a total knee replacement, the associated tibial insert, the resected bones, and the acrylic bone cement layer. In one analysis of this model, all the materials were taken to have isotropic elastic properties, while in the other analysis, the cortical bone was represented as transversely isotropic, with all other materials being considered to have isotropic elastic properties. The result show that, in both study cases considered, there is a strong effect of the representation of the elastic properties of material(s) on the stresses in the models.

Elasticity↗

Morphological changes in pedal phalanges through ornithopod dinosaur evolution: a biomechanical approach.

The evolution of ornithopod dinosaurs provides a well-documented example of the transition from digitigrady to subunguligrady. During this transition, the ornithopod pes was drastically altered from the plesiomorphic dinosaurian morphology (four digits, claw-shaped unguals, strongly concavo-convex joints, phalanges longer than wide, excavated collateral ligament fossae, presence of sagittal ridge, and prominent processes for the attachment of tendons) to a more derived condition (tridactyly, modification of the unguals into hooves, phalanges wider and thinner than long, lack of collateral ligament fossae, loss of sagittal ridge and tendon attachment processes, relatively flattened articular surfaces). These changes are particularly noteworthy given the overall conservatism in pedal morphology seen across Dinosauria. But what are the functional consequences of these specific morphological transitions? To study them, we examine a wide range of pedal morphologies in four non-avian dinosaurs and two birds. Our analyses of the external morphology, two-dimensional models (using Finite Element Analysis), and internal bone structure demonstrate that this evolutionary shift was accompanied by a loss of digit mobility and flexibility. In addition, pedal posture was modified to better align the pes with the main direction of the ground reaction force, thus becoming well suited to support high loads. These conclusions can be applied to other, parallel evolutionary changes (in both dinosaurs and mammals) that involved similar transitions to a subunguligrade posture.

Animals↗

Mechanical evaluation of a bio-active bone cement for total hip arthroplasty.

A new bio-active bone cement, known as CAP, has been developed as an alternative to acrylic bone cement. CAP has improved mechanical properties, with a high modulus that is over five times that of PMMA. The effects of this high modulus are examined by finite element analysis, when the CAP is used in place of PMMA to fix the femoral component in total hip prostheses. The results show a higher tensile stress of 8.76 MPa in the CAP cement, compared with 1.99 MPa in the PMMA cement. However, it is also shown that CAP has a superior fatigue strength of approximately 40 MPa, obtained from a cyclic loading test.

Arthroplasty, Replacement, Hip↗

Prediction of the strength and fracture location of the femoral neck by CT-based finite-element method: a preliminary study on patients with hip fracture.

The purpose of this study was to predict fracture load and fracture location of the femora by means of the originally developed CT-based finite-element method (FEM). The femora of ten patients with contralateral hip fracture were analyzed to estimate fracture strength and to investigate whether the predicted fracture locations were similar to those of contralateral hip fractures. FEM has been utilized to determine the stress or strain distribution in bones under a certain load. FEM analyses of the strength of the femora in cadavers and patients have been reported, but those of hip fracture patients have not been analyzed. The femora of ten patients with contralateral hip fracture and those of three volunteers were analyzed based on the axial CT images of the whole femora. Prediction of hip fracture load and failure locations was made using CT-based finite-element analysis software. The predicted strength of the patients was less than half that of volunteers, and the predicted fracture lines existed at the proximal femur in all patients. It can be concluded that the FEM analyses adopted in this study are able to predict the fracture locations and load of the femora in patients with hip fracture.

Aged↗

Quasi-automatic 3D finite element model generation for individual single-rooted teeth and periodontal ligament.

The paper demonstrates how to generate an individual 3D volume model of a human single-rooted tooth using an automatic workflow. It can be implemented into finite element simulation. In several computational steps, computed tomography data of patients are used to obtain the global coordinates of the tooth's surface. First, the large number of geometric data is processed with several self-developed algorithms for a significant reduction. The most important task is to keep geometrical information of the real tooth. The second main part includes the creation of the volume model for tooth and periodontal ligament (PDL). This is realized with a continuous free form surface of the tooth based on the remaining points. Generating such irregular objects for numerical use in biomechanical research normally requires enormous manual effort and time. The finite element mesh of the tooth, consisting of hexahedral elements, is composed of different materials: dentin, PDL and surrounding alveolar bone. It is capable of simulating tooth movement in a finite element analysis and may give valuable information for a clinical approach without the restrictions of tetrahedral elements. The mesh generator of FE software ANSYS executed the mesh process for hexahedral elements successfully.

Algorithms↗

Ultrasound velocity and attenuation in relation to maximum trabecula stress in the patella.

The ultrasound velocity and attenuation were examined in 16 sets of human patellae. The average ultrasound velocity of patella was shown to be greater in the superior/inferior direction than in the anterior/posterior and medial/lateral directions. The distribution of bone mineral density (BMD) was also examined. The BMD of the patella varied with location. BMD values were largest at the superior and lateral regions and decreased inferiorly and medially. A two-dimensional finite element analysis was performed on each patella. The maximum von Mises stress occurred along the cortical shell on the non-articular surface. The trabecular von Mises stress existed in the posterior region of the patella. Correlation study showed that patellar BMD was significantly associated with each of three directional ultrasound velocities. The relationship between BMD and ultrasound attenuation was not significant. It was also found that the ultrasound velocity and attenuation were not significantly correlated with the maximum von Mises stress.

Aged↗

An investigation of leakage tracts along stressed suture lines in phantom tissue.

Suturing is still one of the most important procedures used in trauma and surgery environments. It does however rely on piercing other wise healthy tissue. In doing so channels are opened which may act as conduits for bacteria to enter the body or to allow material such as blood or intestinal bacteria to leak from vessels. This work shows by means of finite element analysis and physical modelling through the use of a phantom that the tension in the suture has a direct bearing on the opening of these channels. It was found that for the phantom the channel cross-sectional area was approximately proportional to the applied tension and suture thread deflection. Leakage channels of up to 0.1mm(2) cross-section were achieved with suture tensions of 120 grammes force. Channels of this size are significantly larger than the bacteria suggesting transmission is theoretically possible.

Contraindications↗

A finite element model of apical force distribution from orthodontic tooth movement.

This study was undertaken to determine the types of orthodontic forces that cause high stress at the root apex. A 3-dimensional finite element model of a maxillary central incisor, its periodontal ligament (PDL), and alveolar bone was constructed on the basis of average anatomic morphology. The maxillary central incisor was chosen for study because it is one of the teeth at greatest risk for apical root resorption. The material properties of enamel, dentin, PDL, and bone and 5 different load systems (tipping, intrusion, extrusion, bodily movement, and rotational force) were tested. The finite element analysis showed that purely intrusive, extrusive, and rotational forces had stresses concentrated at the apex of the root. The principal stress from a tipping force was located at the alveolar crest. For bodily movement, stress was distributed throughout the PDL; however, it was concentrated more at the alveolar crest. We conclude that intrusive, extrusive, and rotational forces produce more stress at the apex. Bodily movement and tipping forces concentrate forces at the alveolar crest, not at the apex.

Alveolar Process↗

A finite-element approach for Young's modulus reconstruction.

Modulus imaging has great potential in soft-tissue characterization since it reveals intrinsic mechanical properties. A novel Young's modulus reconstruction algorithm that is based on finite-element analysis is reported here. This new method overcomes some limitations in other Young's modulus reconstruction methods. Specifically, it relaxes the force boundary condition requirements so that only the force distribution at the compression surface is necessary, thus making the new method more practical. The validity of the new method is demonstrated and the performance of the algorithm with noise in the input data is tested using numerical simulations. Details of how to apply this method under clinical conditions is also discussed.

Algorithms↗

Focal cryogen insults for inducing segmental osteonecrosis: computational and experimental assessments of thermal fields.

Orthopaedic management of femoral head osteonecrosis is a common clinical problem for which there presently is no good solution. Current animal models are inappropriate to study potential new solutions, since it has been difficult to replicate the natural history of structural collapse seen in the human disorder. Recently, progression to collapse was obtained for cryogenically induced osteonecrosis in emus, although the lesions involved were imprecisely controlled in terms of size or location. A new cryo-insult probe is here reported for the purpose of delivering well-prescribed local thermal insults in this new animal model, while minimizing damage to non-targeted regions. Finite element analysis was used to elucidate the influence of operator-controlled parameters upon the temporal/spatial variation of the thermal field. The numerical formulation includes convective heat transfer attributable to tissue bed perfusion. The computational results agreed closely with the results of thermocouple recordings in a companion bench-top experiment. The cryo-insult probe successfully produced segmental lesions in the emu model of sizes comparable to the computed freeze front diameters.

Animals↗

Micromechanically based poroelastic modeling of fluid flow in Haversian bone.

To explore the hypothesis that load-induced fluid flow in bone is a mechano-transduction mechanism in bone adaptation, unit cell micro-mechanical techniques are used to relate the microstructure of Haversian cortical bone to its effective poroelastic properties. Computational poroelastic models are then applied to compute in vitro Haversian fluid flows in a prismatic specimen of cortical bone during harmonic bending excitations over the frequency range of 10(0) to 10(6) Hz. At each frequency considered, the steady state harmonic response of the poroelastic bone specimen is computed using complex frequency-domain finite element analysis. At the higher frequencies considered, the breakdown of Poisueille flow in Haversian canals is modeled by introduction of a complex fluid viscosity. Peak bone fluid pressures are found to increase linearly with loading frequency in proportion to peak bone stress up to frequencies of approximately 10 kHz. Haversian fluid shear stresses are found to increase linearly with excitation frequency and loading magnitude up until the breakdown of Poisueille flow. Tan delta values associated with the energy dissipated by load-induced fluid flow are also compared with values measured experimentally in a concurrent broadband spectral analysis of bone. The computational models indicate that fluid shear stresses and fluid pressures in the Haversian system could, under physiologically realistic loading, easily reach the level of a few Pascals, which have been shown in other works to elicit cell responses in vitro.

Adaptation, Physiological↗

Predicting the risk of rupture of abdominal aortic aneurysms by utilizing various geometrical parameters: revisiting the diameter criterion.

The authors estimated noninvasively the wall stress distribution for actual abdominal aortic aneurysms (AAAs) in vivo on a patient-to-patient basis and correlated the peak wall stress (PWS) with various geometrical parameters. They studied 39 patients (37 men, mean age 73.7 +/- 8.2 years) with an intact AAA (mean diameter 6.3 +/- 1.7 cm) undergoing preoperative evaluation with spiral computed tomography (CT). Real 3-dimensional AAA geometry was obtained from image processing. Wall stress was determined by using a finite-element analysis. The aorta was considered isotropic with linear material properties and was loaded with a static pressure of 120.0 mm Hg. Various geometrical parameters were used to characterize the AAAs. PWS and each of the geometrical characteristics were correlated by use of Pearson's rank correlation coefficients. PWS varied from 10.2 to 65.8 N/cm2 (mean value 37.1 +/- 9.9 N/cm2). Among the geometrical parameters, the PWS was well correlated with the mean centerline curvature, the maximum centerline curvature, and the maximum centerline torsion of the AAAs. The correlation of PWS with maximum diameter was nonsignificant. Multiple regression analysis revealed that the mean centerline curvature of the AAA was the only significant predictor of PWS and subsequent rupture risk. This noninvasive computational approach showed that geometrical parameters other than the maximum diameter are better indicators of AAA rupture.

Aged↗

A three-dimensional finite analysis of adaptive remodelling in the proximal femur.

A finite element analysis of adaptive bone remodelling in the proximal femur is presented. The use of a three-dimensional model permits a realistic representation of femur geometry, and also allows the possibility of examining the effects of fully three-dimensional loading situations. The long-term pattern of remodelling shows a realistic evolution of density distribution, with a tendency towards a steady state, though the simplified load cases used to model gait are not sufficient to predict the formation of the cortical shell.

Adaptation, Physiological↗

Semiautomated finite element mesh generation methods for a long bone.

The objective of this work was to develop and test a semi-automated finite element mesh generation method using computed tomography (CT) image data of a canine radius. The present study employs a direct conversion from CT Hounsfield units to elastic moduli. Our method attempts to minimize user interaction and eliminate the need for mesh smoothing to produce a model suitable for finite element analysis. Validation of the computational model was conducted by loading the CT-imaged canine radius in four-point bending and using strain gages to record resultant strains that were then compared to strains calculated with the computational model. Geometry-based and uniform modulus voxel-based models were also constructed from the same imaging data set and compared. The nonuniform voxel-based model most accurately predicted the axial strain response of the sample bone (R(2)=0.9764).

Animals↗

Mandibular stiffness in humans: numerical predictions.

The chin is a feature unique to humans. This study evaluates the effect of mandibular symphyseal design on biomechanical masticatory effectiveness as determined by structural stiffness and stress developed under flexural and torsional loading. A simple model of three symphyseal shapes (chin, flat symphysis and lingual buttress), was built to represent human, Neanderthal and higher primate symphyses and these were subjected to wishboning and torsional forces. Additionally, an anatomically detailed reconstruction was made of the CT scan of an actual human mandible, which was then also morphed into a chinless model. The results of a 3-D finite element analysis show firstly, that none of the three different symphyseal shapes is biomechanically more advantageous than the others for the given loading condition. Secondly, we show in a CT-derived model, that the presence of a chin does not confer significantly improved stiffness to torsional or flexural loading. These results indicate that the acquisition of a chin in modern humans is not related to the functional demands placed upon the mandible during mastication, but suggest that it may have developed in response to other biomechanical demands.

Animals↗

Fracturing rigid materials.

We propose a novel approach to fracturing (and denting) brittle materials. To avoid the computational burden imposed by the stringent time step restrictions of explicit methods or with solving nonlinear systems of equations for implicit methods, we treat the material as a fully rigid body in the limit of infinite stiffness. In addition to a triangulated surface mesh and level set volume for collisions, each rigid body is outfitted with a tetrahedral mesh upon which finite element analysis can be carried out to provide a stress map for fracture criteria. We demonstrate that the commonly used stress criteria can lead to arbitrary fracture (especially for stiff materials) and instead propose the notion of a time averaged stress directly into the FEM analysis. When objects fracture, the virtual node algorithm provides new triangle and tetrahedral meshes in a straightforward and robust fashion. Although each new rigid body can be rasterized to obtain a new level set, small shards can be difficult to accurately resolve. Therefore, we propose a novel collision handling technique for treating both rigid bodies and rigid body thin shells represented by only a triangle mesh.

Compressive Strength↗