Morphing the feature-based multi-blocks of normative/healthy vertebral geometries to scoliosis vertebral geometries: development of personalized finite element models

Personalized Finite Element (FE) models and hexahedral elements are preferred for biomechanical investigations. Feature-based multi-block methods are used to develop anatomically accurate personalized FE models with hexahedral mesh. It is tedious to manually construct multi-blocks for large number of geometries on an individual basis to develop personalized FE models. Mesh-morphing method mitigates the aforementioned tediousness in meshing personalized geometries every time, but leads to element warping and loss of geometrical data. Such issues increase in magnitude when normative spine FE model is morphed to scoliosis-affected spinal geometry. The only way to bypass the issue of hex-mesh distortion or loss of geometry as a result of morphing is to rely on manually constructing the multi-blocks for scoliosis-affected spine geometry of each individual, which is time intensive. A method to semi-automate the construction of multi-blocks on the geometry of scoliosis vertebrae from the existing multi-blocks of normative vertebrae is demonstrated in this paper. High-quality hexahedral elements were generated on the scoliosis vertebrae from the morphed multi-blocks of normative vertebrae. Time taken was 3 months to construct the multi-blocks for normative spine and less than a day for scoliosis. Efforts taken to construct multi-blocks on personalized scoliosis spinal geometries are significantly reduced by morphing existing multi-blocks.     

Analysis of pelvic strain in different gait configurations in a validated cohort of computed tomography based finite element models

The pelvis functions to transmit upper body loads to the lower limbs and is critical in human locomotion. Semi-automated, landmark-based finite element (FE) morphing and mapping techniques eliminate the need for segmentation and have shown to accelerate the generation of multiple specimen-specific pelvic FE models to enable the study of pelvic mechanical behaviour. The purpose of this research was to produce an experimentally validated cohort of specimen-specific FE models of the human pelvis and to use this cohort to analyze pelvic strain patterns during gait. Using an initially segmented specimen-specific pelvic FE model as a source model, four more specimen-specific pelvic FE models were generated from target clinical CT scans using landmark-based morphing and mapping techniques. FE strains from the five models were compared to the experimental strains obtained from cadaveric testing via linear regression analysis, (R² values ranging from 0.70 to 0.93). Inter-specimen variability in FE strain distributions was seen among the five specimen-specific pelvic FE models. The validated cohort of specimen-specific pelvic FE models was utilized to examine pelvic strains at different phases of the gait cycle. Each validated specimen-specific FE model was reconfigured into gait cycle phases representing heel-strike/heel-off and midstance/midswing. No significant difference was found in the double-leg stance and heel-strike/heel-off models (p = 0.40). A trend was observed between double-leg stance and midstance/midswing models (p = 0.07), and a significant difference was found between heel-strike/heel-off models and midstance/midswing models (p = 0.02). Significant differences were also found in comparing right vs. left models (heel-strike/heel-off p = 0.14, midstance/midswing p = 0.04).     

Morphing methods to parameterize specimen-specific finite element model geometries

Shape plays an important role in determining the biomechanical response of a structure. Specimen-specific finite element (FE) models have been developed to capture the details of the shape of biological structures and predict their biomechanics. Shape, however, can vary considerably across individuals or change due to aging or disease, and analysis of the sensitivity of specimen-specific models to these variations has proven challenging. An alternative to specimen-specific representation has been to develop generic models with simplified geometries whose shape is relatively easy to parameterize, and can therefore be readily used in sensitivity studies. Despite many successful applications, generic models are limited in that they cannot make predictions for individual specimens.We propose that it is possible to harness the detail available in specimen-specific models while leveraging the power of the parameterization techniques common in generic models. In this work we show that this can be accomplished by using morphing techniques to parameterize the geometry of specimen-specific FE models such that the model shape can be varied in a controlled and systematic way suitable for sensitivity analysis. We demonstrate three morphing techniques by using them on a model of the load-bearing tissues of the posterior pole of the eye. We show that using relatively straightforward procedures these morphing techniques can be combined, which allows the study of factor interactions. Finally, we illustrate that the techniques can be used in other systems by applying them to morph a femur. Morphing techniques provide an exciting new possibility for the analysis of the biomechanical role of shape, independently or in interaction with loading and material properties.     

An automated method to morph finite element whole-body human models with a wide range of stature and body shape for both men and women

Field data analyses have shown that small female, obese, and/or older occupants are at increased risks of death and serious injury in motor-vehicle crashes compared with mid-size young men. The current adult finite element (FE) human models represent occupants in the same three body sizes (large male, mid-size male, and small female) as those for the contemporary adult crash dummies. Further, the time needed to develop an FE human model using the traditional method is measured in months or even years. In the current study, an improved regional mesh morphing method based on landmark-based radial basis function (RBF) interpolation was developed to rapidly morph a mid-size male FE human model into different geometry targets. A total of 100 human models with a wide range of human attributes were generated. A pendulum chest impact condition was applied to each model as an initial assessment of the resulting variability in response. The morphed models demonstrated mesh quality similar to the baseline model. The peak impact forces and chest deflections in the chest pendulum impacts varied substantially with different models, supportive of consideration of population variation in evaluating the occupant injury risks. The method developed in this study will enable future safety design optimizations targeting at various vulnerable populations that cannot be considered with the current models.     

Effects of obesity on occupant responses in frontal crashes: a simulation analysis using human body models

The objective of this study is to investigate the effects of obesity on occupant responses in frontal crashes using whole-body human finite element (FE) models representing occupants with different obesity levels. In this study, the geometry of THUMS 4 midsize male model was varied using mesh morphing techniques with target geometries defined by statistical models of external body contour and exterior ribcage geometry. Models with different body mass indices (BMIs) were calibrated against cadaver test data under high-speed abdomen loading and frontal crash conditions. A parametric analysis was performed to investigate the effects of BMI on occupant injuries in frontal crashes based on the Taguchi method while controlling for several vehicle design parameters. Simulations of obese occupants predicted significantly higher risks of injuries to the thorax and lower extremities in frontal crashes compared with non-obese occupants, which is consistent with previous field data analyses. These higher injury risks are mainly due to the increased body mass and relatively poor belt fit caused by soft tissues for obese occupants. This study demonstrated the feasibility of using a parametric human FE model to investigate the obesity effects on occupant responses in frontal crashes.     

Subject specific finite element mesh generation of the pelvis from biplanar x-ray images: application to 120 clinical cases

Several Finite Element (FE) models of the pelvis have been developed to comprehensively assess the onset of pathologies and for clinical and industrial applications. However, because of the difficulties associated with the creation of subject-specific FE mesh from CT scan and MR images, most of the existing models rely on the data of one given individual. Moreover, although several fast and robust methods have been developed for automatically generating tetrahedral meshes of arbitrary geometries, hexahedral meshes are still preferred today because of their distinct advantages but their generation remains an open challenge. Recently, approaches have been proposed for fast 3D reconstruction of bones based on X-ray imaging. In this study, we adapted such an approach for the fast and automatic generation of all-hexahedral subject-specific FE models of the pelvis based on the elastic registration of a generic mesh to the subject-specific target in conjunction with element regularity and quality correction. The technique was successfully tested on a database of 120 3D reconstructions of pelvises from biplanar X-ray images. For each patient, a full hexahedral subject-specific FE mesh was generated with an accurate surface representation.     

Growing multiblock structures: a semi-automated approach to block placement for multiblock hexahedral meshing

Finite element (FE) analysis is a cornerstone of orthopaedic biomechanics research. Three-dimensional medical imaging provides sufficient resolution for the subject-specific FE models to be generated from these data-sets. FE model development requires discretisation of a three-dimensional domain, which can be the most time-consuming component of a FE study. Hexahedral meshing tools based on the multiblock method currently rely on the manual placement of building blocks for mesh generation. We hypothesise that angular analysis of the geometric centreline for a three-dimensional surface could be used to automatically generate building block structures for the multiblock hexahedral mesh generation. Our algorithm uses a set of user-defined points and parameters to automatically generate a multiblock structure based on a surface's geometric centreline. This significantly reduces the time required for model development. We have applied this algorithm to 47 bones of varying geometries and successfully generated a FE mesh in all cases. This work represents significant advancement in automatically generating multiblock structures for a wide range of geometries.     

A critical assessment of standard processing methods for the preparation of palynological samples

Standard processing techniques for the isolation of organic walled dinoflagellate cysts from geological samples are examined, with particular attention to the size and type of sieve mesh used. Variations within the ‘standard’ processing techniques used by different laboratories are identified, and an assessment of the retention capacities of meshes of different sizes and different materials is carried out. Some dinoflagellate cysts and large numbers of Lycopodium spores, used for the calculations of absolute abundance data, were found to pass through 20 μm meshes. This is due to a combination of factors including: the diagonal aperture diameter of a 20 μm mesh measuring over 28 μm; the three-dimensional properties of different mesh weaves (nylon and polyester); and the non-spherical shape of the particles. Experiments demonstrate that the maximum mesh size that should be used in palynological processing is 15 μm. Nylon mesh is more practical to use than polyester as processing time is reduced, but nylon is degraded by contact with acid solutions. Meshes with apertures < 15 μm may be used, though this may be impractical for large samples containing significant quantities of fine siliciclastic or organic material.     

中国近海甲藻环沟藻属2个新记录种:纺锤环沟藻和莫氏环沟藻

甲藻环沟藻属于一类无色素体、表面有脊的裸甲藻, 因可捕食一些重要的赤潮生物而在海洋生态系统中扮演着重要的角色。有关中国近海环沟藻属的物种多样性信息非常有限。本文报道了2个新记录种——纺锤环沟藻(Gyrodinium fusiforme)和莫氏环沟藻(G. moestrupii)。纺锤环沟藻细胞呈纺锤形, 长48.0-58.0 μm, 宽18.0-23.0 μm, 长宽比为2.4-3.0, 和模式种相比体型和长宽比都较小。莫氏环沟藻细胞也呈纺锤形, 长约30 μm, 宽约15 μm。我们测定了纺锤环沟藻和莫氏环沟藻大亚基的部分序列, 并根据大亚基序列利用最大似然法和贝叶斯法建立了系统发育树。结果显示环沟藻属是单源的, 纺锤环沟藻和裂缝环沟藻(G. fissum)聚合在一起, 但是与螺旋环沟藻(G. spirale)分离。纺锤环沟藻和莫氏环沟藻分别可以摄食米氏凯伦藻(Karenia mikimotoi)和具齿原甲藻(Prorocentrum dentatum), 前者在米氏凯伦藻赤潮中的大量出现显示它可以促进赤潮的消退。     

Growing multiblock structures: a semi-automated approach to block placement for multiblock hexahedral meshing

Finite element (FE) analysis is a cornerstone of orthopaedic biomechanics research. Three-dimensional medical imaging provides sufficient resolution for the subject-specific FE models to be generated from these data-sets. FE model development requires discretisation of a three-dimensional domain, which can be the most time-consuming component of a FE study. Hexahedral meshing tools based on the multiblock method currently rely on the manual placement of building blocks for mesh generation. We hypothesise that angular analysis of the geometric centreline for a three-dimensional surface could be used to automatically generate building block structures for the multiblock hexahedral mesh generation. Our algorithm uses a set of user-defined points and parameters to automatically generate a multiblock structure based on a surface's geometric centreline. This significantly reduces the time required for model development. We have applied this algorithm to 47 bones of varying geometries and successfully generated a FE mesh in all cases. This work represents significant advancement in automatically generating multiblock structures for a wide range of geometries.     

Factors affecting subject-specific finite element models of implant-fitted rat bone specimens: critical analysis of a technical protocol

The authors propose a protocol to derive finite element (FE) models from micro computer tomography scans of implanted rat bone. A semi-automatic procedure allows segmenting the images using specimen-specific bone mineral density (BMD) thresholds. An open-source FE model generator processes the segmented images to a quality tetrahedral mesh. The material properties assigned to each element are integrated from the BMD field. Piecewise, threshold-dependent density–elasticity relationships are implemented to limit the effects of metal artefacts. A detailed sensitivity study highlights the coherence of the generated models and quantifies the influence of the modelling parameters on the results. Two applications of the protocol are proposed. The stiffness of bare and implanted rat tibiae specimens is predicted by simulating three-point bending and inter-implant displacement, respectively. Results are compared with experimental tests. The mean value and the variability between the specimens are well captured in both tests.     

The influence of ligament modelling strategies on the predictive capability of finite element models of the human knee joint

In finite element (FE) models knee ligaments can represented either by a group of one-dimensional springs, or by three-dimensional continuum elements based on segmentations. Continuum models closer approximate the anatomy, and facilitate ligament wrapping, while spring models are computationally less expensive. The mechanical properties of ligaments can be based on literature, or adjusted specifically for the subject. In the current study we investigated the effect of ligament modelling strategy on the predictive capability of FE models of the human knee joint. The effect of literature-based versus specimen-specific optimized material parameters was evaluated. Experiments were performed on three human cadaver knees, which were modelled in FE models with ligaments represented either using springs, or using continuum representations. In spring representation collateral ligaments were each modelled with three and cruciate ligaments with two single-element bundles. Stiffness parameters and pre-strains were optimized based on laxity tests for both approaches. Validation experiments were conducted to evaluate the outcomes of the FE models.Models (both spring and continuum) with subject-specific properties improved the predicted kinematics and contact outcome parameters. Models incorporating literature-based parameters, and particularly the spring models (with the representations implemented in this study), led to relatively high errors in kinematics and contact pressures. Using a continuum modelling approach resulted in more accurate contact outcome variables than the spring representation with two (cruciate ligaments) and three (collateral ligaments) single-element-bundle representations. However, when the prediction of joint kinematics is of main interest, spring ligament models provide a faster option with acceptable outcome.     

Finite element modeling of the human kidney for probabilistic occupant models: Statistical shape analysis and mesh morphing

Statistical shape analysis was conducted on 15 pairs (left and right) of human kidneys. It was shown that the left and right kidney were significantly different in size and shape. In addition, several common modes of kidney variation were identified using statistical shape analysis. Semi-automatic mesh morphing techniques have been developed to efficiently create subject specific meshes from a template mesh with a similar geometry. Subject specific meshes as well as probabilistic kidney meshes were created from a template mesh. Mesh quality remained about the same as the template mesh while only taking a fraction of the time to create the mesh from scratch or morph with manually identified landmarks. This technique can help enhance the quality of information gathered from experimental testing with subject specific meshes as well as help to more efficiently predict injury by creating models with the mean shape as well as models at the extremes for each principal component.     

The influence of material property and morphological parameters on specimen-specific finite element models of porcine vertebral bodies

The use of finite element (FE) methods in spinal research is increasing, but there is only limited information available on the influence of different input parameters on the model predictions. The aim of this study was to investigate the role of these parameters in FE models of the vertebral body. Experimental tests were undertaken on porcine lumbar vertebral bodies and scans of the specimens were used to create specimen-specific FE models. Three models were created for each specimen with combinations of generic and specimen-specific parameters. Stiffness and strength predictions were also made directly from the specimen trabecular bone volume fraction (BVF) and cross-sectional area (CSA). The agreement between the experimental results and the FE models with generic morphology was poorer (concordance coefficients = 0.058, 0.125 for stiffness, strength) than those made from the BVF and CSA (concordance coefficients = 0.638, 0.609). The greatest levels of agreement were found with the morphologically specific models including element-specific material properties (concordance coefficients = 0.881, 0.752). This indicates that highly specific models, both in terms of morphology and bone quality, are necessary if the FE tool is to be used effectively for spinal research and clinical practice.     

An eFTD-VP framework for efficiently generating patient-specific anatomically detailed facial soft tissue FE mesh for craniomaxillofacial surgery simulation

Accurate surgical planning and prediction of craniomaxillofacial surgery outcome requires simulation of soft tissue changes following osteotomy. This can only be achieved by using an anatomically detailed facial soft tissue model. The current state-of-the-art of model generation is not appropriate to clinical applications due to the time-intensive nature of manual segmentation and volumetric mesh generation. The conventional patient-specific finite element (FE) mesh generation methods are to deform a template FE mesh to match the shape of a patient based on registration. However, these methods commonly produce element distortion. Additionally, the mesh density for patients depends on that of the template model. It could not be adjusted to conduct mesh density sensitivity analysis. In this study, we propose a new framework of patient-specific facial soft tissue FE mesh generation. The goal of the developed method is to efficiently generate a high-quality patient-specific hexahedral FE mesh with adjustable mesh density while preserving the accuracy in anatomical structure correspondence. Our FE mesh is generated by eFace template deformation followed by volumetric parametrization. First, the patient-specific anatomically detailed facial soft tissue model (including skin, mucosa, and muscles) is generated by deforming an eFace template model. The adaptation of the eFace template model is achieved by using a hybrid landmark-based morphing and dense surface fitting approach followed by a thin-plate spline interpolation. Then, high-quality hexahedral mesh is constructed by using volumetric parameterization. The user can control the resolution of hexahedron mesh to best reflect clinicians’ need. Our approach was validated using 30 patient models and 4 visible human datasets. The generated patient-specific FE mesh showed high surface matching accuracy, element quality, and internal structure matching accuracy. They can be directly and effectively used for clinical simulation of facial soft tissue change.     

Automated muscle wrapping using finite element contact detection

Realistic muscle path representation is essential to musculoskeletal modeling of joint function. Algorithms predicting these muscle paths typically rely on a labor intensive predefinition of via points or underlying geometries to guide wrapping for given joint positions. While muscle wrapping using anatomically precise three-dimensional (3D) finite element (FE) models of bone and muscle has been achieved, computational expense and pre-processing associated with this approach exclude its use in applications such as subject-specific modeling. With the intention of combining advantageous features of both approaches, an intermediate technique relying on contact detection capabilities of commercial FE packages is presented. We applied the approach to the glenohumeral joint, and validated the method by comparison against existing experimental data. Individual muscles were modeled as a straight series of deformable beam elements and bones as anatomically precise 3D rigid bodies. Only the attachment locations and a default orientation of the undeformed muscle segment were pre-defined. The joint was then oriented in a static position of interest. The muscle segment free end was then moved along the shortest Euclidean path to its origin on the scapula, wrapping the muscle along bone surfaces by relying on software contact detection. After wrapping for a given position, the resulting moment arm was computed as the perpendicular distance from the line of action vector to the humeral head center of rotation.This approach reasonably predicted muscle length and moment arm for 27 muscle segments when compared to experimental measurements over a wide range of shoulder motion. Artificial via points or underlying contact geometries were avoided, contact detection and multiobject wrapping on the bone surfaces were automatic, and low computational cost permitted wrapping of individual muscles within seconds on a standard desktop PC. These advantages may be valuable for both general and subject-specific musculoskeletal modeling.     

Finite element analysis of trabecular bone structure: a comparison of image-based meshing techniques

In this study, we investigate if finite element (FE) analyses of human trabecular bone architecture based on 168 μm images can provide relevant information about the bone mechanical characteristics. Three human trabecular bone samples, one taken from the femoral head, one from the iliac crest, and one from the lumbar spine, were imaged with micro-computed tomography (micro-CT) using a 28 μm resolution. After reconstruction the resolution was coarsened to 168 μm. First, all reconstructions were thresholded and directly converted to FE-models built of hexahedral elements. For the coarser resolutions of two samples, this resulted in a loss of trabecular connections and a subsequent loss of stiffness. To reduce this effect, a tetrahedral element meshing based on the marching cubes algorithm, as well as a modified hexahedron meshing, which thresholds the image such that load carrying bone mass is preserved, were employed. For each sample elastic moduli and tissue Von Mises stresses of the three different 168 μm models were compared to those from the hexahedron 28 μm model. For one sample the hexahedron meshing at 168 μm produced excellent results. For the other two samples the results obtained from the hexahedral models at 168 μm resolution were poor. Considerably better results were attained for these samples when using the mass-compensated or tetrahedron meshing techniques. We conclude that the accuracy of the FE-models at 168 μm strongly depends on the bone morphology, in particular its trabecular thickness. A substantial loss of trabecular connections during the hexahedron meshing process indicates that poor FE results will be obtained. In this case the tetrahedron or mass-compensated hexahedron meshing techniques can reduce the loss of connections and produce better results than the plain hexahedron meshing techniques.     

APPLICATION OF THIN‐PLATE SPLINE TRANSFORMATIONS TO FINITE ELEMENT MODELS,OR, HOW TO TURN A BOG TURTLE INTO A SPOTTED TURTLE TO ANALYZE BOTH

Finite element (FE) models are popular tools that allow biologists to analyze the biomechanical behavior of complex anatomical structures. However, the expense and time required to create models from specimens has prevented comparative studies from involving large numbers of species. A new method is presented for transforming existing FE models using geometric morphometric methods. Homologous landmark coordinates are digitized on the FE model and on a target specimen into which the FE model is being transformed. These coordinates are used to create a thin‐plate spline function and coefficients, which are then applied to every node in the FE model. This function smoothly interpolates the location of points between landmarks, transforming the geometry of the original model to match the target. This new FE model is then used as input in FE analyses. This procedure is demonstrated with turtle shells: a Glyptemys muhlenbergii model is transformed into Clemmys guttata and Actinemys marmorata models. Models are loaded and the resulting stresses are compared. The validity of the models is tested by crushing actual turtle shells in a materials testing machine and comparing those results to predictions from FE models. General guidelines, cautions, and possibilities for this procedure are also presented.     

3D characterization of bone strains in the rat tibia loading model

Bone strain is considered one of the factors inducing bone tissue response to loading. Nevertheless, where animal studies can provide detailed data on bone response, they only offer limited information on experimental bone strains. Including micro-CT-based finite element (micro FE) models in the analysis represents a potent methodology for quantifying strains in bone. Therefore, the main objective of this study was to develop and validate specimen-specific micro FE models for the assessment of bone strains in the rat tibia compression model. Eight rat limbs were subjected to axial compression loading; strain at the medio-proximal site of the tibiae was measured by means of strain gauges. Specimen-specific micro FE models were created and analyzed. Repeated measurements on each limb indicated that the effect of limb positioning was small (COV?= 6.45 ± 2.27 %). Instead, the difference in the measured strains between the animals was high (54.2%). The computational strains calculated at the strain gauge site highly correlated to the measured strains (R ²?=?0.95). Maximum peak strains calculated at exactly 25% of the tibia length for all specimens were equal to 435.11 ± 77.88 microstrains (COV?=?17.19%). In conclusion, we showed that strain gauge measurements are very sensitive to the exact strain gauge location on the bone; hence, the use of strain gauge data only is not recommended for studies that address at identifying reliable relationships between tissue response and local strains. Instead, specimen-specific micro FE models of rat tibiae provide accurate estimates of tissue-level strains.