Computational design and evaluation of the mechanical and electrical behavior of a piezoelectric scaffold: a preclinical study.

Piezoelectric scaffolds have been recently developed to explore their potential to enhance the bone regeneration process using the concept of piezoelectricity, which also inherently occurs in bone. In addition to providing mechanical support during bone healing, with a suitable design, they are supposed to produce electrical signals that ought to favor the cell responses. In this study, using finite element analysis (FEA), a piezoelectric scaffold was designed with the aim of providing favorable ranges of mechanical and electrical signals when implanted in a large bone defect in a large animal model, so that it could inform future pre-clinical studies. A parametric analysis was then performed to evaluate the effect of the scaffold design parameters with regard to the piezoelectric behavior of the scaffold. The designed scaffold consisted of a porous strut-like structure with piezoelectric patches covering its free surfaces within the scaffold pores. The results showed that titanium or PCL for the scaffold and barium titanate (BT) for the piezoelectric patches are a promising material combination to generate favorable ranges of voltage, as reported in experimental studies. Furthermore, the analysis of variance showed the thickness of the piezoelectric patches to be the most influential geometrical parameter on the generation of electrical signals in the scaffold. This study shows the potential of computer tools for the optimization of scaffold designs and suggests that patches of piezoelectric material, attached to the scaffold surfaces, can deliver favorable ranges of electrical stimuli to the cells that might promote bone regeneration.

A review on computer modeling of bone piezoelectricity and its application to bone adaptation and regeneration.

Bone is a hierarchical, multiphasic and anisotropic structure which in addition possess piezoelectric properties. The generation of piezoelectricity in bone is a complex process which has been shown to play a key role both in bone adaptation and regeneration. In order to understand the complex biological, mechanical and electrical interactions that take place during these processes, several computer models have been developed and used to test hypothesis on potential mechanisms behind experimental observations. This paper aims to review the available literature on computer modeling of bone piezoelectricity and its application to bone adaptation and healing. We first provide a brief overview of the fundamentals of piezoelectricity and bone piezoelectric effects. We then review how these properties have been used in computational models of bone adaptation and electromechanical behaviour of bone. In addition, in the last section, we summarize current limitations and potential directions for future work.

In Silico Optimization of Femoral Fixator Position and Configuration by Parametric CAD Model.

Structural analysis, based on the finite element method, and structural optimization, can help surgery planning or decrease the probability of fixator failure during bone healing. Structural optimization implies the creation of many finite element model instances, usually built using a computer-aided design (CAD) model of the bone-fixator assembly. The three most important features of such CAD models are: parameterization, robustness and bidirectional associativity with finite elements (FE) models. Their significance increases with the increase in the complexity of the modeled fixator. The aim of this study was to define an automated procedure for the configuration and placement of fixators used in the treatment of long bone fractures. Automated and robust positioning of the selfdynamisable internal fixator on the femur was achieved and sensitivity analysis of fixator stress on the change of major design parameters was performed. The application of the proposed methodology is considered to be beneficial in the preparation of CAD models for automated structural optimization procedures used in long bone fixation.

Dealing with Nap-Core Sandwich Composites: How to Predict the Effect of Symmetry.

The behavior of nap-core sandwiches was investigated with a special focus on the effect of symmetry in nap cores. A nap-core is, in general terms, a 3D-formed hollow structure made of knitted textile impregnated by a thermosetting resin. The molding process determines if the nap-core is double-sided (symmetric) or single-sided. The sandwich with nap-core owns various remarkable properties of a novel lightweight material, but the nap-core's complex structure makes the prediction of these properties a difficult task. While the analysis of a single-sided nap-core sandwich has been presented by the authors before, this study is focused on the simulation of symmetric nap-core sandwich. Overall, performance of the structure is examined with respect to several loading conditions. The simulation approach invokes a typical homogenization scheme to find the engineering properties of the nap-core's fabric with least computational time and memory resources. Results from experiments and simulations exhibit a good compatibility, which prove the fitness of the modeling method.

Selecting boundary conditions in physiological strain analysis of the femur: Balanced loads, inertia relief method and follower load.

Selection of boundary constraints may influence amount and distribution of loads. The purpose of this study is to analyze the potential of inertia relief and follower load to maintain the effects of musculoskeletal loads even under large deflections in patient specific finite element models of intact or fractured bone compared to empiric boundary constraints which have been shown to lead to physiological displacements and surface strains. The goal is to elucidate the use of boundary conditions in strain analyses of bones. Finite element models of the intact femur and a model of clinically relevant fracture stabilization by locking plate fixation were analyzed with normal walking loading conditions for different boundary conditions, specifically re-balanced loading, inertia relief and follower load. Peak principal cortex surface strains for different boundary conditions are consistent (maximum deviation 13.7%) except for inertia relief without force balancing (maximum deviation 108.4%). Influence of follower load on displacements increases with higher deflection in fracture model (from 3% to 7% for force balanced model). For load balanced models, follower load had only minor influence, though the effect increases strongly with higher deflection. Conventional constraints of fixed nodes in space should be carefully reconsidered because their type and position are challenging to justify and for their potential to introduce relevant non-physiological reaction forces. Inertia relief provides an alternative method which yields physiological strain results.