A Patient-Specific Computer Tomography-Based Finite Element Methodology to Calculate the Six Dimensional Stiffness Matrix of Human Vertebral Bodies

2012 ◽  
Vol 134 (5) ◽  
Author(s):  
Yan Chevalier ◽  
Philippe K. Zysset
Keyword(s):  
Finite Element ◽  
Structural Properties ◽  
Vertebral Body ◽  
Computer Tomography ◽  
Stiffness Matrix ◽  
Anisotropic Elasticity ◽  
Patient Specific ◽  
Axial Stiffness ◽  
Element Analysis ◽  
Vertebral Bodies

In most finite element (FE) studies of vertebral bodies, axial compression is the loading mode of choice to investigate structural properties, but this might not adequately reflect the various loads to which the spine is subjected during daily activities or the increased fracture risk associated with shearing or bending loads. This work aims at proposing a patient-specific computer tomography (CT)-based methodology, using the currently most advanced, clinically applicable finite element approach to perform a structural investigation of the vertebral body by calculation of its full six dimensional (6D) stiffness matrix. FE models were created from voxel images after smoothing of the peripheral voxels and extrusion of a cortical shell, with material laws describing heterogeneous, anisotropic elasticity for trabecular bone, isotropic elasticity for the cortex based on experimental data. Validated against experimental axial stiffness, these models were loaded in the six canonical modes and their 6D stiffness matrix calculated. Results show that, on average, the major vertebral rigidities correlated well or excellently with the axial rigidity but that weaker correlations were observed for the minor coupling rigidities and for the image-based density measurements. This suggests that axial rigidity is representative of the overall stiffness of the vertebral body and that finite element analysis brings more insight in vertebral fragility than densitometric approaches. Finally, this extended patient-specific FE methodology provides a more complete quantification of structural properties for clinical studies at the spine.

Experimental Validation of Finite Element Analysis of Human Vertebral Collapse Under Large Compressive Strains

2014 ◽  
Vol 136 (4) ◽  
Author(s):  
Hadi S. Hosseini ◽  
Allison L. Clouthier ◽  
Philippe K. Zysset
Keyword(s):  
Finite Element Analysis ◽  
Finite Element ◽  
Trabecular Bone ◽  
Axial Compression ◽  
Vertebral Body ◽  
Computer Tomography ◽  
Strain Localization ◽  
Vertebral Collapse ◽  
Element Analysis ◽  
Custom Made

Osteoporosis-related vertebral fractures represent a major health problem in elderly populations. Such fractures can often only be diagnosed after a substantial deformation history of the vertebral body. Therefore, it remains a challenge for clinicians to distinguish between stable and progressive potentially harmful fractures. Accordingly, novel criteria for selection of the appropriate conservative or surgical treatment are urgently needed. Computer tomography-based finite element analysis is an increasingly accepted method to predict the quasi-static vertebral strength and to follow up this small strain property longitudinally in time. A recent development in constitutive modeling allows us to simulate strain localization and densification in trabecular bone under large compressive strains without mesh dependence. The aim of this work was to validate this recently developed constitutive model of trabecular bone for the prediction of strain localization and densification in the human vertebral body subjected to large compressive deformation. A custom-made stepwise loading device mounted in a high resolution peripheral computer tomography system was used to describe the progressive collapse of 13 human vertebrae under axial compression. Continuum finite element analyses of the 13 compression tests were realized and the zones of high volumetric strain were compared with the experiments. A fair qualitative correspondence of the strain localization zone between the experiment and finite element analysis was achieved in 9 out of 13 tests and significant correlations of the volumetric strains were obtained throughout the range of applied axial compression. Interestingly, the stepwise propagating localization zones in trabecular bone converged to the buckling locations in the cortical shell. While the adopted continuum finite element approach still suffers from several limitations, these encouraging preliminary results towardsthe prediction of extended vertebral collapse may help in assessing fracture stability in future work.

Patient-Specific 3D Finite-Element Analysis of Miniscrew Implants During Orthodontic Treatment

2009 ◽  
Author(s):  
Hussein H. Ammar ◽  
Victor H. Mucino ◽  
Peter Ngan ◽  
Richard J. Crout ◽  
Osama M. Mukdadi
Keyword(s):  
Finite Element ◽  
Orthodontic Treatment ◽  
3D Model ◽  
Patient Specific ◽  
Stress Concentrations ◽  
Element Analysis ◽  
Miniscrew Implants ◽  
Operative Planning ◽  
Miniscrew Implant ◽  
Maximum Stresses

Miniscrew implants have seen increasing clinical use as orthodontic anchorage devices with demonstrated stability. The focus of this study is to develop and simulate operative factors, such as load magnitudes and anchor locations to achieve desired motions in a patient-specific 3D model undergoing orthodontic treatment with miniscrew implant anchorage. A CT scan of a patient skull was imported into Mimics software (Materialise, 12.1). Segmentation operations were performed on the images to isolate the mandible, filter out noise, then reconstruct a smooth 3D model. A model of the left canine was reconstructed with the PDL modeled as a thin solid layer. A miniscrew was modeled with dimensions based on a clinical implant (BMK OAS-T1207) then inserted into the posterior mandible. All components were volumetrically meshed and optimized in Mimics software. Elements comprising the mandible bone and teeth were assigned a material based on their gray value ranges in HU from the original scan, and meshes were exported into ANSYS software. All materials were defined as linear and isotropic. A nonlinear PDL was also defined for comparison. For transverse forces applied on the miniscrew, maximum stresses increased linearly with loading and appeared at the neck or first thread and in the cortical bone. A distal tipping force was applied on the canine, and maximum stresses appeared in the tooth at the crown and apex and in the bone at the compression surface. Under maximum loading, stresses in bone were sufficient for resorption. The nonlinear PDL exhibited lower stresses and deflections than the linear model due to increasing stiffness. Numerous stress concentrations were seen in all models. Results of this study demonstrate the potential of patient-specific 3D reconstruction from CT scans and finite-element simulation as a versatile and effective pre-operative planning tool for orthodontists.

scholarly journals Nonlinear Axial Stiffness Characteristics of Axisymmetric Bolted Joints

1990 ◽  
Vol 112 (3) ◽  
pp. 442-449 ◽  
Author(s):  
I. R. Grosse ◽  
L. D. Mitchell
Keyword(s):  
Finite Element ◽  
Linear Theory ◽  
Design Theory ◽  
Joint Stiffness ◽  
Bolted Joints ◽  
Axial Stiffness ◽  
Outer Diameter ◽  
Bolted Joint ◽  
Element Analysis ◽  
Stiffness Characteristics

A critical assessment of the current design theory for bolted joints which is based on a linear, one-dimensional stiffness analysis is presented. A detailed nonlinear finite element analysis of a bolted joint conforming to ANSI standards was performed. The finite element results revealed that the joint stiffness is highly dependent on the magnitude of the applied load. The joint stiffness changes continuously from extremely high for small applied loads to the bolt stiffness during large applied loads, contrary to the constant joint stiffness of the linear theory. The linear theory is shown to be inadequate in characterizing the joint stiffness. The significance of the results in terms of the failure of bolted joints is discussed. A number of sensitivity studies were carried out to assess the effect of various parameters on the axial joint stiffness. The results revealed that bending and rotation of the joint members, interfacial friction, and the bolt/nut threading significantly influence the axial stiffness characteristics of the bolted joint. The two-dimensional, axisymmetric finite element model includes bilinear gap elements to model the interfaces. Special orthotropic elements were used to model the bolt/nut thread interaction. A free-body-diagram approach was taken by applying loads to the outer diameter of the joint model which correspond to internal, uniformly distributed line-shear and line-moment loads in the joint. A number of convergence studies were performed to validate the solution.

scholarly journals A Methodology for the Derivation of Unloaded Abdominal Aortic Aneurysm Geometry With Experimental Validation

2016 ◽  
Vol 138 (10) ◽  
Author(s):  
Santanu Chandra ◽  
Vimalatharmaiyah Gnanaruban ◽  
Fabian Riveros ◽  
Jose F. Rodriguez ◽  
Ender A. Finol
Keyword(s):  
Finite Element ◽  
Abdominal Aortic Aneurysm ◽  
Aortic Aneurysm ◽  
Stress Component ◽  
Wall Stress ◽  
Patient Specific ◽  
Element Analysis ◽  
Nodal Surface ◽  
Abdominal Aortic ◽  
Peak Wall Stress

In this work, we present a novel method for the derivation of the unloaded geometry of an abdominal aortic aneurysm (AAA) from a pressurized geometry in turn obtained by 3D reconstruction of computed tomography (CT) images. The approach was experimentally validated with an aneurysm phantom loaded with gauge pressures of 80, 120, and 140 mm Hg. The unloaded phantom geometries estimated from these pressurized states were compared to the actual unloaded phantom geometry, resulting in mean nodal surface distances of up to 3.9% of the maximum aneurysm diameter. An in-silico verification was also performed using a patient-specific AAA mesh, resulting in maximum nodal surface distances of 8 μm after running the algorithm for eight iterations. The methodology was then applied to 12 patient-specific AAA for which their corresponding unloaded geometries were generated in 5–8 iterations. The wall mechanics resulting from finite element analysis of the pressurized (CT image-based) and unloaded geometries were compared to quantify the relative importance of using an unloaded geometry for AAA biomechanics. The pressurized AAA models underestimate peak wall stress (quantified by the first principal stress component) on average by 15% compared to the unloaded AAA models. The validation and application of the method, readily compatible with any finite element solver, underscores the importance of generating the unloaded AAA volume mesh prior to using wall stress as a biomechanical marker for rupture risk assessment.

scholarly journals Effect of coronal plane acetabular correction on joint contact pressure in Periacetabular osteotomy: a finite-element analysis

2022 ◽  
Vol 23 (1) ◽  
Author(s):  
Kenji Kitamura ◽  
Masanori Fujii ◽  
Miho Iwamoto ◽  
Satoshi Ikemura ◽  
Satoshi Hamai ◽  
...  
Keyword(s):  
Finite Element Analysis ◽  
Finite Element ◽  
Contact Pressure ◽  
Periacetabular Osteotomy ◽  
Anterior Wall ◽  
Coronal Plane ◽  
Patient Specific ◽  
Element Analysis ◽  
Joint Contact ◽  
Joint Contact Pressure

Abstract Background The ideal acetabular position for optimizing hip joint biomechanics in periacetabular osteotomy (PAO) remains unclear. We aimed to determine the relationship between acetabular correction in the coronal plane and joint contact pressure (CP) and identify morphological factors associated with residual abnormal CP after correction. Methods Using CT images from 44 patients with hip dysplasia, we performed three patterns of virtual PAOs on patient-specific 3D hip models; the acetabulum was rotated laterally to the lateral center-edge angles (LCEA) of 30°, 35°, and 40°. Finite-element analysis was used to calculate the CP of the acetabular cartilage during a single-leg stance. Results Coronal correction to the LCEA of 30° decreased the median maximum CP 0.5-fold compared to preoperatively (p <  0.001). Additional correction to the LCEA of 40° further decreased CP in 15 hips (34%) but conversely increased CP in 29 hips (66%). The increase in CP was associated with greater preoperative extrusion index (p = 0.030) and roundness index (p = 0.038). Overall, virtual PAO failed to normalize CP in 11 hips (25%), and a small anterior wall index (p = 0.049) and a large roundness index (p = 0.003) were associated with residual abnormal CP. Conclusions The degree of acetabular correction in the coronal plane where CP is minimized varied among patients. Coronal plane correction alone failed to normalize CP in 25% of patients in this study. In patients with an anterior acetabular deficiency (anterior wall index < 0.21) and an aspherical femoral head (roundness index > 53.2%), coronal plane correction alone may not normalize CP. Further studies are needed to clarify the effectiveness of multiplanar correction, including in the sagittal and axial planes, in optimizing the hip joint’s contact mechanics.

Incorporation of Spinal Flexibility Measurements Into Finite Element Analysis

1990 ◽  
Vol 112 (4) ◽  
pp. 481-483 ◽  
Author(s):  
Mack G. Gardner-Morse ◽  
Jeffrey P. Laible ◽  
Ian A. F. Stokes
Keyword(s):  
Experimental Data ◽  
Finite Element Analysis ◽  
Finite Element ◽  
Stiffness Matrix ◽  
Beam Element ◽  
Technical Note ◽  
Element Analysis ◽  
Spinal Motion ◽  
Spinal Flexibility

This technical note demonstrates two methods of incorporating the experimental stiffness of spinal motion segments into a finite element analysis of the spine. The first method is to incorporate the experimental data directly as a stiffness matrix. The second method approximates the experimental data as a beam element.

Sign in / Sign up

Export Citation Format

Share Document