Naomi Tsafnat

An experimentally validated micromechanical model of a rat vertebra under compressive loading.

In recent years, finite element analysis (FEA) has been increasingly applied to examine and predict the mechanical behaviour of craniofacial and other bony structures. Traditional methods used to determine material properties and validate finite element models (FEMs) have met with variable success, and can be time‐consuming. An implicit assumption underlying many FE studies is that relatively high localized stress/strain magnitudes identified in FEMs are likely to predict material failure. Here we present a new approach that may offer some advantages over previous approaches. Recently developed technology now allows us to both image and conduct mechanical tests on samples in situ using a materials testing stage (MTS) fitted inside the microCT scanner. Thus, micro‐finite element models can be created and validated using both quantitative and qualitative means. In this study, a rat vertebra was tested under compressive loading until failure using an MTS. MicroCT imaging of the vertebra before mechanical testing was used to create a high resolution finite element model of the vertebra. Load‐displacement data recorded during the test were used to calculate the effective Young’s modulus of the bone (found to be 128 MPa). The microCT image of the compressed vertebra was used to assess the predictive qualities of the FE model. The model showed the highest stress concentrations in the areas that failed during the test. Clearly, our analyses do not directly address biomechanics of the craniofacial region; however, the methodology adopted here could easily be applied to examine the properties and behaviour of specific craniofacial structures, or whole craniofacial regions of small vertebrates. Experimentally validated micro‐FE analyses are a powerful method in the study of materials with complex microstructures such as bone.

Modelling heating of liver tumours with heterogeneous magnetic microsphere deposition

Ferromagnetic embolization hyperthermia (FEH) is a novel treatment for liver cancer. Magnetic microspheres are injected into the hepatic artery and cluster in the periphery of tumours and are heated with externally applied magnetic fields. In order to more accurately simulate FEH, we modelled a three-dimensional heterogeneous distribution of heat sources. We constructed a fractal model of the vasculature in the periphery of a tumour. We used this model to compute the spatial distribution of the microspheres that lodge in capillaries. We used the distribution model as input to a finite-element heat transfer model of the FEH treatment. The overall appearance of the vascular tree is subjectively similar to that of the disorganized vascular network which encapsulates tumours. The microspheres are distributed in the tumour periphery in similar patterns to experimental observations. We expect the vasculature and microsphere deposition models to also be of interest to researchers of any targeted cancer therapies such as localized intra-arterial chemotherapy and selective internal radiotherapy. Our results show that heterogeneous microsphere distributions give significantly different results to those for a homogeneous model and thus are preferable when accurate results are required.

Micromechanics of Sea Urchin Spines

The endoskeletal structure of the Sea Urchin, Centrostephanus rodgersii, has numerous long spines whose known functions include locomotion, sensing, and protection against predators. These spines have a remarkable internal microstructure and are made of single-crystal calcite. A finite-element model of the spine’s unique porous structure, based on micro-computed tomography (microCT) and incorporating anisotropic material properties, was developed to study its response to mechanical loading. Simulations show that high stress concentrations occur at certain points in the spine’s architecture; brittle cracking would likely initiate in these regions. These analyses demonstrate that the organization of single-crystal calcite in the unique, intricate morphology of the sea urchin spine results in a strong, stiff and lightweight structure that enhances its strength despite the brittleness of its constituent material.

A three-dimensional fractal model of tumour vasculature

We constructed a three-dimensional fractal model of the vascular network in a tumour periphery. We model the highly disorganised structure of the neoplastic vasculature by using a high degree of variation in segment properties such as length, diameter and branching angle. The overall appearance of the vascular tree is subjectively similar to that of the disorganised vascular network which encapsulates tumours. The fractal dimension of the model is within the range of clinically measured values.

Pedicle screw-based posterior dynamic stabilizers for degenerative spine:In vitrobiomechanical testing and clinical outcomes: Pedicle Screw-based Posterior Dynamic Stabilizers for Degenerative Spine

Dynamic stabilization in a degenerate symptomatic spine may be advantageous compared with conventional fusion procedures, as it helps preserve motion and minimizes redistribution of loads at instrumented and adjacent segments. This article presents a systematic review of biomechanical and clinical evidence available on some of the pedicle screw based posterior dynamic stabilization (PDS) devices. Using Medline, Embase, and Scopus online databases, we identified four pedicle-screw-PDS devices for which both, biomechanical testing and clinical follow-up data are available: Graf artificial ligaments, Isobar TTL, Polyetheretherketone rods, and Dynesys. The current state-of-the-art of pedicle-screw-PDS devices is far from achieving its desired biomechanical efficacy, which has resulted in a weak support for the posited clinical benefits. Although pedicle-screw-PDS devices are useful in salvaging a moderately degenerate functionally suboptimal disc, for severe disc degeneration cases fusion is still the preferred choice. We conclude that a pedicle-screw-PDS device should aim at restoring load sharing amongst spinal elements while preserving the qualitative and quantitative nature of spinal motion, especially minimize posterior shift of the helical axis of motion. More precise and objective assessment techniques need to be standardized for in vivo evaluation of intervertebral motion and load sharing amongst spinal elements across different pedicle-screw-PDS devices.

Computational aerodynamic analysis of a Micro-CT based bio-realistic fruit fly wing

The aerodynamic features of a bio-realistic 3D fruit fly wing in steady state (snapshot) flight conditions were analyzed numerically. The wing geometry was created from high resolution micro-computed tomography (micro-CT) of the fruit fly Drosophila virilis. Computational fluid dynamics (CFD) analyses of the wing were conducted at ultra-low Reynolds numbers ranging from 71 to 200, and at angles of attack ranging from -10° to +30°. It was found that in the 3D bio-realistc model, the corrugations of the wing created localized circulation regions in the flow field, most notably at higher angles of attack near the wing tip. Analyses of a simplified flat wing geometry showed higher lift to drag performance values for any given angle of attack at these Reynolds numbers, though very similar performance is noted at -10°. Results have indicated that the simplified flat wing can successfully be used to approximate high-level properties such as aerodynamic coefficients and overall performance trends as well as large flow-field structures. However, local pressure peaks and near-wing flow features induced by the corrugations are unable to be replicated by the simple wing. We therefore recommend that accurate 3D bio-realistic geometries be used when modelling insect wings where such information is useful.

The Role of Sacral Slope in the Progression of a Bilateral Spondylolytic Defect at L5 to Spondylolisthesis: A Biomechanical Investigation Using Finite Element Analysis

Study Design: A biomechanical study using finite element analysis. Objectives: The main objective of this study was to investigate the role of sacral slope in the progression of a L5 bilateral spondylolytic defect to spondylolisthesis. Methods: A 3-dimensional model of lumbosacral spine was built using computed tomography (CT) data procured from an anonymized healthy male subject. The segmented CT data was manipulated to generate 3 more models representing L5 bilateral spondylolytic defect with normal sacral slope (SS), sacral slope increased by 10° (SS+10), and sacral slope decreased by 10° (SS-10). The 3D models were imported into finite element modelling software Strand7 for preprocessing, running nonlinear static solves, and postprocessing of the results. Results: Directional biomechanical instabilities were induced in the lumbosacral spine as a result of changes in the L5-S1 disc shape secondary to the changes in sacral slope. Compared with the normal L5 lytic model, wedging of the L5-S1 disc (SS+10) resulted in a significantly greater range of motion in flexion (18% ↑) but extension motion characteristics were similar. Conversely, flattening of the L5-S1 disc (SS-10) resulted in a significantly greater range of motion in extension (16% ↑) but flexion motion characteristics were similar to that of the normal L5 lytic model. Conclusions: Variations in sacral slope while preserving the L5-S1 mid-disc height and orientation of the L5 vertebra resulted in variations in the L5-S1 disc shape. The results suggest that for such extremities in the L5-S1 disc shape different pathomechanisms exist for the progression of the L5 lytic defect to spondylolisthesis.

Global and segmental kinematic changes following sequential resection of posterior osteoligamentous structures in the lumbar spine: An in vitro biomechanical investigation using pure moment testing protocols

Lumbar spinal surgeries may compromise the integrity of posterior osteoligamentous structures implicating mechanical stability. Circumstances necessitating a concomitant surgery to achieve restabilisation are not well understood. The main objective of this in vitro study was to quantify global and segmental (index and adjacent levels) kinematic changes in the lumbar spine following sequential resection of the posterior osteoligamentous structures using pure moment testing protocols. Six fresh frozen cadaveric kangaroo lumbar spines (T12–S1) were tested under a bending moment in flexion–extension, bilateral bending, and axial torsion in a 6-degree-of-freedom Kinematic Spine Simulator. Specimens were tested in the following order: intact state (D0), after interspinous and supraspinous ligaments transection between L4 and L5 (D1), further after a total bilateral facetectomy between L4 and L5 (D2). Segmental motions at the cephalad, damaged, and caudal levels were recorded using an infrared-based motion tracking device. Following D1, no significant change in the global range of motion was observed in any of the bending planes. Following D2, a significant increase in the global range of motion from the baseline (D0) was observed in axial torsion (median normalised change +20%). At the damaged level, D2 resulted in a significant increase in the segmental range of motion in flexion–extension (+77%) and axial torsion (+492%). Additionally, a significant decrease in the segmental range of motion in axial torsion (−35%) was observed at the caudal level following D2. These results suggest that a multi-segment lumbar spine acts as a mechanism for transmitting motions, and that a compromised joint may significantly alter motion transfer to adjacent segments. We conclude that the interspinous and supraspinous ligaments play a modest role in restricting global spinal motions within physiologic limits. Following interspinous and supraspinous ligaments transection, a total bilateral facetectomy resulted in a significant increase in axial torsion motion, both at global and damaged levels, accompanied with a compensatory decrease in motion at the caudal level.

Modeling Dental Implants With Finite Element Analysis and Micro-Computed Tomography

X-ray micro-computed tomography (microCT) allows us to construct three-dimensional images of specimens at the micron scale in a non-destructive manner. The digital nature of the microCT images, which are in voxel form, make them ideal candidates for use in numerical modeling and simulation [1]. Finite element analysis (FEA) is a well-known technique for modeling the structural response of a system to mechanical loading, and is most useful in modeling complex systems which cannot be analyzed analytically. MicroCT datasets can be converted into finite element models, directly incorporating both the geometry of the specimen and information about the different materials in it. This method is known as micro-finite element analysis (microFEA). It is especially useful in the study of materials with complex microstructures.Copyright