Intravoxel bone micromechanics for microCT-based finite element simulations

While micro-FE simulations have become a standard tool in computational biomechanics, the choice of appropriate material properties is still a relevant topic, typically involving empirical grey value-to-elastic modulus relations. We here derive the voxel-specific volume fractions of mineral, collagen, and water, from tissue-independent bilinear relations between mineral and collagen content in extracellular bone tissue (J. Theor. Biol. 287: 115, 2011), and from the measured X-ray attenuation information quantified in terms of grey values. The aforementioned volume fractions enter a micromechanics representation of bone tissue, as to deliver voxel-specific stiffness tensors. In order to check the relevance of this strategy, we convert a micro Computer Tomograph of a mouse femur into a regular Finite Element mesh, apply forces related to the dead load of a standing mouse, and then compare simulation results based on voxel-specific heterogeneous elastic properties to results based on homogeneous elastic properties related to the spatial average over the solid bone matrix compartment, of the X-ray attenuation coefficients. The element-specific strain energy density illustrates that use of homogeneous elastic properties implies overestimation of the organ stiffness. Moreover, the simulation reveals large tensile normal stresses throughout the femur neck, which may explain the mouse femur neck's trabecular morphology being quite different from the human case, where the femur neck bears compressive forces and bending moments.     

A penetration-based finite element method for hyperelastic 3D biphasic tissues in contact. Part II: finite element simulations

The penetration method allows for the efficient finite element simulation of contact between soft hydrated biphasic tissues in diarthrodial joints. Efficiency of the method is achieved by separating the intrinsically nonlinear contact problem into a pair of linked biphasic finite element analyses, in which an approximate, spatially and temporally varying contact traction is applied to each of the contacting tissues. In Part I of this study, we extended the penetration method to contact involving nonlinear biphasic tissue layers, and demonstrated how to derive the approximate contact traction boundary conditions. The traction derivation involves time and space dependent natural boundary conditions, and requires special numerical treatment. This paper (Part II) describes how we obtain an efficient nonlinear finite element procedure to solve for the biphasic response of the individual contacting layers. In particular, alternate linearization of the nonlinear weak form, as well as both velocity-pressure, v-p, and displacement-pressure, u-p, mixed formulations are considered. We conclude that the u-p approach, with linearization of both the material law and the deformation gradients, performs best for the problem at hand. The nonlinear biphasic contact solution will be demonstrated for the motion of the glenohumeral joint of the human shoulder joint.     

Nanoindentation testing and finite element simulations of cortical bone allowing for anisotropic elastic and inelastic mechanical response

Anisotropy is one of the most peculiar aspects of cortical bone mechanical behaviour, and the numerical approach can be successfully used to investigate aspects of bone tissue mechanics that analytical methods solve in approximate way or do not cover. In this work, nanoindentation experimental tests and finite element simulations were employed to investigate the elastic-inelastic anisotropic mechanical properties of cortical bone. The model allows for anisotropic elastic and post-yield behaviour of the tissue. A tension-compression mismatch and direction-dependent yield stresses are allowed for. Indentation experiments along the axial and transverse directions were simulated with the purpose to predict the indentation moduli and hardnesses along multiple orientations. Results showed that the experimental transverse-to-axial ratio of indentation moduli, equal to 0.74, is predicted with a ～3% discrepancy regardless the post-yield material behaviour; whereas, the transverse-to-axial hardness ratio, equal to 0.86, can be correctly simulated (discrepancy ～6% w.r.t. the experimental results) only employing an anisotropic post-elastic constitutive model. Further, direct comparison between the experimental and simulated indentation tests evidenced a good agreement in the loading branch of the indentation curves and in the peak loads for a transverse-to-axial yield stress ratio comparable to the experimentally obtained transverse-to-axial hardness ratio. In perspective, the present work results strongly support the coupling between indentation experiments and FEM simulations to get a deeper knowledge of bone tissue mechanical behaviour at the microstructural level. The present model could be used to assess the effect of variations of constitutive parameters due to age, injury, and/or disease on bone mechanical performance in the context of indentation testing.     

Maxillary expansion treatment using bone anchors: development and validation of a 3D finite element model

OBJECTIVE: Develop a finite element (FE) model of a skull to perform biomechanical studies of maxillary expansion using bone anchors (BA). MATERIALS AND METHODS: A skull model was developed and assigned material properties based on Hounsfield unit (HU) values of cone-beam computerized tomography (CBCT) images. A 3 mm diameter cylindrical BA was modelled and inserted in the palatal bone. A 4 mm transverse displacement was applied on the anchor. An evaluation on the effect on local stresses of BA implantation inclination angle was performed. RESULTS: Proper displacement results and strain-stress trends for the expansion process were present. Stress distribution patterns were similar as reported in the literature. No significant difference between BA inclination angles was found. CONCLUSION: This work leads to a better understanding and prediction of craniofacial and maxillary bone remodelling during ME with BA treatments and is a first step towards the development of patient specific treatments.     

Improved method for coliform verification

Modification of a method for coliform verification presented in Standard Methods for the Examination of Water and Wastewater is described. Modification of the method, which is based on beta-galactosidase production, involves incorporation of a lactose operon inducer in medium upon which presumptive coliform isolates are cultured prior to beta-galactosidase assay.     

Maxillary expansion treatment using bone anchors: development and validation of a 3D finite element model

Objective: Develop a finite element (FE) model of a skull to perform biomechanical studies of maxillary expansion using bone anchors (BA).

Materials and methods: A skull model was developed and assigned material properties based on Hounsfield unit (HU) values of cone-beam computerized tomography (CBCT) images. A 3 mm diameter cylindrical BA was modelled and inserted in the palatal bone. A 4 mm transverse displacement was applied on the anchor. An evaluation on the effect on local stresses of BA implantation inclination angle was performed.

Results: Proper displacement results and strain–stress trends for the expansion process were present. Stress distribution patterns were similar as reported in the literature. No significant difference between BA inclination angles was found.

Conclusion: This work leads to a better understanding and prediction of craniofacial and maxillary bone remodelling during ME with BA treatments and is a first step towards the development of patient specific treatments.     

Targeted quantitation of CVD-linked plasma proteins for biomarker verification and validation

Despite significant advances in treatment, cardiovascular disease (CVD) remains one of the leading causes of morbidity and mortality in developed and developing countries. Judicious monitoring of common risk factors has been unable to control this global epidemic, necessitating novel biomarkers for improved screening and earlier disease detection and management. Although numerous plasma proteins have been associated with CVD, only a few of these potential biomarkers have been validated for clinical use. Here we review the quantitative proteomic methods used to verify and validate new biomarker candidates in human plasma. These methods center on a bottom-up approach involving multiple or selected reaction monitoring, for targeted detection, with stable isotope-labeled standards, for peptide normalization. Also included are a discussion of future strategies for improved CVD protein biomarker verification and validation, recommendations for method translation to the clinic, and future projections for protein biomarker research.     

RETRACTED: Modelling subcortical bone in finite element analyses: A validation and sensitivity study in the macaque mandible

This article has been retracted: please see Elsevier Policy on Article Withdrawal (http://www.elsevier.com/locate/withdrawalpolicy).This article has been retracted at the request of the authors who unanimously wish to retract this paper because of errors in the validation protocol and data, and the finite element analysis data and results.     

The steady-state finite cable: Numerical method for non-linear membrane

An iterative numerical method is described for finding steady-state solutions to the one-dimensional cable equations for finite cable lengths and open-circuit termination. The method is suitable for any non-linear membrane with a single positive-slope crossing of the zero current axis, including those with regions of negative slope conductance. The method generates the current necessary to cause any desired voltage displacement at the input end of the cable as well as solutions for the transmembrane potential and axial current along the cable.     

Improved method for coliform verification.

Modification of a method for coliform verification presented in Standard Methods for the Examination of Water and Wastewater is described. Modification of the method, which is based on beta-galactosidase production, involves incorporation of a lactose operon inducer in medium upon which presumptive coliform isolates are cultured prior to beta-galactosidase assay.     

Use of micro-CT-based finite element analysis to accurately quantify peri-implant bone strains: a validation in rat tibiae

Although research has been addressed at investigating the effect of specific loading regimes on bone response around the implant, a precise quantitative understanding of the local mechanical response close to the implant site is still lacking. This study was aimed at validating micro-CT-based finite element (μFE) models to assess tissue strains after implant placement in a rat tibia. Small implants were inserted at the medio-proximal site of 8 rat tibiae. The limbs were subjected to axial compression loading; strain close to the implant was measured by means of strain gauges. Specimen-specific μFE models were created and analyzed. For each specimen, 4 different models were created corresponding to different representations of the bone–implant interface: bone and implant were assumed fully osseointegrated (A); a low stiffness interface zone was assumed with thickness of 40 μm (B), 80 μm (C), and 160 μm (D). In all cases, measured and computational strains correlated highly (R ² = 0.95, 0.92, 0.93, and 0.95 in A, B, C, and D, respectively). The averaged calculated strains were 1.69, 1.34, and 1.15 times higher than the measured strains for A, B, and C, respectively, and lower than the experimental strains for D (factor = 0.91). In conclusion, we demonstrated that specimen-specific FE analyses provide accurate estimates of peri-implant bone strains in the rat tibia loading model. Further investigations of the bone-implant interface are needed to quantify implant osseointegration.     

A finite element method for mechanical response of hair cell ciliary bundles

This paper describes the development of a methodology for performing a mechanical analysis of hair cell ciliary bundles. The cilia were modeled as shear deformable beams, and interconnections were modeled as two-force members. These models were incorporated into software, which performs a finite element analysis of a user-defined bundle. The algorithm incorporates aspects of the bundle such as geometric realignment and buckling of compressed side links. A sample bundle is introduced and results of modeling it are presented.     

Electrophoretic affinity chromatography: method validation

A new method for preparative-scale separation of biomolecules, electrophoretic affinity chromatography (EAC), is proposed in this paper. Separation by EAC is carried out in a long and ribbon-like multicompartment electrolyser separated by membranes, in which the two central compartments are used for packing the gel matrix and for sample loading respectively. Next to the central compartments are the elution compartments and electrode compartments. The electric field is applied perpendicular to the fluid flow in the compartments. Adsorption and desorption steps may both be carried out in the presence of an electric field, which transports the target components into the gel compartment for adsorption and the impurities into the elution compartments for washing. After the adsorption step an elution solution is introduced and the product is released from the gel matrix and washed out. Separation of human serum albumin (HSA) from human serum gives HSA product of high purity, as demonstrated by isoelectric focusing analysis. The characteristics of electrophoretic binding of HSA on Blue Sepharose Fast Flow are examined. The preliminary results show that this new method has advantages in terms of high rate of mass transfer and ease of scaling up, which are of particular interest when large-scale separation of biomolecules is considered.     

Homogenized stiffness matrices for mineralized collagen fibrils and lamellar bone using unit cell finite element models

Mineralized collagen fibrils have been usually analyzed like a two-phase composite material where crystals are considered as platelets that constitute the reinforcement phase. Different models have been used to describe the elastic behavior of the material. In this work, it is shown that when Halpin–Tsai equations are applied to estimate elastic constants from typical constituent properties, not all crystal dimensions yield a model that satisfy thermodynamic restrictions. We provide the ranges of platelet dimensions that lead to positive definite stiffness matrices. On the other hand, a finite element model of a mineralized collagen fibril unit cell under periodic boundary conditions is analyzed. By applying six canonical load cases, homogenized stiffness matrices are numerically calculated. Results show a monoclinic behavior of the mineralized collagen fibril. In addition, a 5-layer lamellar structure is also considered where crystals rotate in adjacent layers of a lamella. The stiffness matrix of each layer is calculated applying Lekhnitskii transformations, and a new finite element model under periodic boundary conditions is analyzed to calculate the homogenized 3D anisotropic stiffness matrix of a unit cell of lamellar bone. Results are compared with the rule-of-mixtures showing in general good agreement.     

Development and validation of subject-specific finite element models for blunt trauma study

This study developed and validated finite element (FE) models of swine and human thoraxes and abdomens that had subject-specific anatomies and could accurately and efficiently predict body responses to blunt impacts. Anatomies of the rib cage, torso walls, thoracic, and abdominal organs were reconstructed from X-ray computed tomography (CT) images and extracted into geometries to build FE meshes. The rib cage was modeled as an inhomogeneous beam structure with geometry and bone material parameters determined directly from CT images. Meshes of soft components were generated by mapping structured mesh templates representative of organ topologies onto the geometries. The swine models were developed from and validated by 30 animal tests in which blunt insults were applied to swine subjects and CT images, chest wall motions, lung pressures, and pathological data were acquired. A comparison of the FE calculations of animal responses and experimental measurements showed a good agreement. The errors in calculated response time traces were within 10% for most tests. Calculated peak responses showed strong correlations with the experimental values. The stress concentration inside the ribs, lungs, and livers produced by FE simulations also compared favorably to the injury locations. A human FE model was developed from CT images from the Visible Human project and was scaled to simulate historical frontal and side post mortem human subject (PMHS) impact tests. The calculated chest deformation also showed a good agreement with the measurements. The models developed in this study can be of great value for studying blunt thoracic and abdominal trauma and for designing injury prevention techniques, equipments, and devices.     

Singularity solution of Lanir's osmoelasticity: verification of discontinuity simulations in soft tissues

The literature characterizes cartilaginous tissues as osmoviscoelastic. Understanding the damage and failure of these tissues is essential for designing treatments. To determine tissue strength and local stresses, experimental studies—both clinical and animal—are generally supported by computational studies. Verification methods for computational studies of ionized porous media including cracks are hardly available. This study provides a method for verification and shows its performance. For this purpose, shear loading of a finite crack is addressed analytically and through a commercial finite element code. Impulsive shear loading by two-edge dislocation of a crack was considered in a 2D plane strain model for an ionized porous medium. To derive the analytical solution, the system of equation is decoupled by stress functions. The shear stress distribution at the plane of the crack is derived using Fourier and Laplace transformations. The analytical solution for the shear stress distribution is compared with computer simulations in ABAQUS version 6.4-5. Decoupling of the equations makes it possible to solve some boundary value problems in porous media taking chemical effects into account. The numerical calculations underestimate the shear stress at the crack-tips. Mesh refinement increases accuracy, but is still low in the neighborhood of the crack-tips.     

Microscale poroelastic metamodel for efficient mesoscale bone remodelling simulations

Bone functional tissue adaptation is a multiaspect physiological process driven by interrelated mechanical and biological stimuli which requires the combined activity of osteoclasts and osteoblasts. In previous work, the authors developed a phenomenological mesoscale structural modelling approach capable of predicting internal structure of the femur based on daily activity loading, which relied on the iterative update of the cross-sectional areas of truss and shell elements representative of trabecular and cortical bones, respectively. The objective of this study was to introduce trabecular reorientation in the phenomenological model at limited computational cost. To this aim, a metamodel derived from poroelastic microscale continuum simulations was used to predict the functional adaptation of a simplified proximal structural femur model. Clear smooth trabecular tracts are predicted to form in the regions corresponding to the main trabecular groups identified in literature, at minimal computational cost.     

The mechanics of heterotypic cell aggregates: insights from computer simulations

Finite element-based computer simulations are used to investigate a number of phenomena, including tissue engulfment, cell sorting, and checkerboard-pattern formation, exhibited by heterotypic cell aggregates. The simulations show that these phenomena can be driven by a single equivalent force, namely a surface (or interfacial) tension, that results from cytoskeletal components and cell-cell adhesions. They also reveal that tissue engulfment, cell sorting, and checkerboard-pattern formation involve several discernible mechanical features or stages. With the aid of analytical arguments, we identify the conditions necessary for each of these phenomena. These findings are consistent with previous experimental investigations and computer simulations, but pose significant challenges to current theories of cell sorting and tissue engulfment.