Finite element analysis of cancellous bone failure in the vertebral body of healthy and osteoporotic subjects.

The aim of this work is to assess the fracture risk prediction of the cancellous bone in the body of a lumbar vertebra when the mechanical parameters of the bone, i.e. stiffness, porosity, and strength anisotropy, of elderly and osteoporotic subjects are considered. For this purpose, a non-linear three-dimensional continuum-based finite element model of the lumbar functional spinal unit L4-L5 was created and strength analyses of the spongy tissue of the vertebral body were carried out. A fabric-dependent strength criterion, which accounts for the micro-architecture of the cancellous bone, based on histomorphometric analyses was used. The strength analyses have shown that the cancellous bone of none of the subject types undergoes failure under loading applied during normal daily life like axial compression; however, bone failure occurs for the osteoporotic segment, subjected to a combination of the compression preloading and moments in the sagittal or in the frontal plane, which are conditions that may not be considered to occur 'daily'. In particular, critical stress conditions are met because of the high porosity values in the horizontal direction within the cancellous bone. The computational approach presented in the paper can potentially predict the material fracture risk of the cancellous bone in the vertebral body and it may be usefully employed to draw failure maps representing, for a given micro-architecture of the spongy tissue, the critical loading conditions (forces and moments) that may lead to such a risk. This approach could be further developed in order to assess the effectiveness of biomedical devices within an engineering approach to the clinical problem of the spinal diseases.

Computational modeling of combined cell population dynamics and oxygen transport in engineered tissue subject to interstitial perfusion.

This work presents a computational model of tissue growth under interstitial perfusion inside a tissue engineering bioreactor. The model accounts both for the cell population dynamics, using a model based on cellular automata, and for the hydrodynamic microenvironment imposed by the bioreactor, using a model based on the Lattice-Boltzmann equation and the convection-diffusion equation. The conditions of static culture versus perfused culture were compared, by including the population dynamics along with oxygen diffusion, convective transport and consumption. The model is able to deal with arbitrary complex geometries of the spatial domain; in the present work, the domain modeled was the void space of a porous scaffold for tissue-engineered cartilage. The cell population dynamics algorithm provided results which qualitatively resembled population dynamics patterns observed in experimental studies, and these results were in good quantitative agreement with previous computational studies. Simulation of oxygen transport and consumption showed the fundamental contribution of convective transport in maintaining a high level of oxygen concentration in the whole spatial domain of the scaffold. The model was designed with the aim to be computationally efficient and easily expandable, i.e. to allow straightforward implementability of further models of complex biological phenomena of increasing scientific interest in tissue engineering, such as chemotaxis, extracellular matrix deposition and effect of mechanical stimulation.

Stress distribution in a post-restored tooth using the three-dimensional finite element method.

Clinicians are opting ever more frequently for restorative materials which have an elastic modulus similar to that of dentin when reconstructing endodontically treated teeth. Metallic posts, which are capable of causing dangerous and non-homogenous stresses in root dentin, are slowly being abandoned. Ideal posts may be those made of various types of fibre (carbon, mineral and glass) and which are adhesively luted into the canal. Among the different methods for evaluating the mechanical behaviour of posts in root canals (progressive loads and photo-elastic technique) the finite element method (FEM) presents many advantages. The aim of this paper is to evaluate, utilizing three-dimensional analysis of the finite elements, what the effect of material rigidity, depth of insertion and post diameter could be on the stress distribution in the different components of the single tooth-post-core reconstruction unit. The results of the FEM analyses, expressed as the distribution of Von Mises stress values, has allowed us to conclude that (i) fibreglass-reinforced composite distributes stress better than titanium alloy or stainless steel; (ii) fibreglass-reinforced composite posts should be inserted as deeply as possible (but maintaining 5-6 mm of gutta-percha apical seal); (iii) fibreglass-reinforced composite post diameter does not affect stress distribution, therefore, as much radicular dentin as possible should be preserved.

A comparative evaluation of chondrocyte/scaffold constructs for cartilage tissue engineering.

This study aimed to evaluate three biodegradable scaffolds as cell carriers for in vitro cartilage regeneration using mature human chondrocyte cells. We compared cell distribution, viability and morphology and we evaluated the mechanical properties of the constructs after 2 weeks of in vitro culture. The materials used as scaffolds were fibrin glue, a collagen sponge and a polyurethane foam (DegraPol(R)). Fibrin glue was found unsuitable as a chondrocyte carrier vehicle after culture times longer than a few days, probably due to significant barriers to nutrients and oxygen diffusion, and the material weakened rapidly. The collagen-based sponge was found to be unsuitable to support chondrocyte survival in vitro, although the presence of newly synthesized collagen was observed in these constructs. The synthetic biodegradable scaffold was more adequate in supporting cell survival and mechanical properties. After 2 weeks of static culture, the storage modulus obtained by dynamic shear testing was in the order of 0.7 kPa in fibrin constructs, 3.7 kPa in collagen constructs and 105 kPa in DegraPol(R) constructs. The better mechanical stability of the synthetic foam supports further investigation in the possible use of synthetic biomaterials as biodegradable scaffolds for in vitro cartilage regeneration. (Journal of Applied Biomaterials & Biomechanics 2004; 2: 55-64).

Elastic stabilization alone or combined with rigid fusion in spinal surgery: a biomechanical study and clinical experience based on 82 cases.

The authors report their experience with the treatment of lumbar instability by a kind of spine stabilization. The elastic stabilization, which follows a new philosophy, is obtained by an interspinous device, and should be used alone in degenerative disc disease, recurrent disc herniation and in very low grade instability, or in association with rigid fusion for the prevention of pathology of the border area. In collaboration with bioengineers, we carried out an experimental study on a lumbar spine model in order to calculate stresses and deformations of lumbar disc during simulation of motion, in physiological conditions and when elastic stabilization is combined with rigid fusion. Results suggest that elastic stabilization reduces stresses on the adjacent disc up to 28 degrees of flexion. Based on this preliminary result, we began to use elastic stabilization alone or combined with fusion in 1994. To date, we have performed 82 surgical procedures, 57 using stabilization alone and 25 combined with fusion, in patients affected by degenerative disc disease, disc herniation, recurrence of disc herniation or other pathologies. Clinical results are satisfactory, especially in the group of patients affected by recurrent disc herniation, in whom the elastic device was used alone.

Multiscale modelling as a tool to prescribe realistic boundary conditions for the study of surgical procedures.

This work was motivated by the problems of analysing detailed 3D models of vascular districts with complex anatomy. It suggests an approach to prescribing realistic boundary conditions to use in order to obtain information on local as well as global haemodynamics. A method was developed which simultaneously solves Navier-Stokes equations for local information and a non-linear system of ordinary differential equations for global information. This is based on the principle that an anatomically detailed 3D model of a cardiovascular district can be achieved by using the finite element method. In turn the finite element method requires a specific boundary condition set. The approach outlined in this work is to include the system of ordinary differential equations in the boundary condition set. Such a multiscale approach was first applied to two controls: (i) a 3D model of a straight tube in a simple hydraulic network and (ii) a 3D model of a straight coronary vessel in a lumped-parameter model of the cardiovascular system. The results obtained are very close to the solutions available for the pipe geometry. This paper also presents preliminary results from the application of the methodology to a particular haemodynamic problem: namely the fluid dynamics of a systemic-to-pulmonary shunt in paediatric cardiac surgery.

Mechanobiology of engineered cartilage cultured under a quantified fluid-dynamic environment.

Natural cartilage remodels both in vivo and in vitro in response to mechanical forces and hence mechanical stimulation is believed to have a potential as a tool to modulate extra-cellular matrix synthesis in tissue-engineered cartilage. Fluid-induced shear is known to enhance chondrogenesis on animal cells. A well-defined hydrodynamic environment is required to study the biochemical response to shear of three-dimensional engineered cell systems. We have developed a perfused-column bioreactor in which the culture medium flows through chondrocyte-seeded porous scaffolds, together with a computational fluid-dynamic model of the flow through the constructs' microstructure. A preliminary experiment of human chondrocyte growth under static versus dynamic conditions is described. The median shear stress imposed on the cells in the bioreactor culture, as predicted by the CFD model, is 3 x 10(-3) Pa (0.03 dyn/cm(2)) at a flow rate of 0.5 ml/min corresponding to an inlet fluid velocity of 44.2 mum/s. Providing a fluid-dynamic environment to the cells yielded significant differences in cell morphology and in construct structure.

Improved mathematical model of the wear of the cup articular surface in hip joint prostheses and comparison with retrieved components.

This paper presents an analytical model of the cobalt-based alloy-ultra-high molecular weight polyethylene (UHMWPE) wear coupling. Based on a previous model in which the cup wear volume over a gait cycle (WG) was calculated under the simplifying assumption of an ideal rigid coupling, the current version proposes a more realistic wear simulation. All three components of the hip loading force were considered for the contact pressure calculation and all three components of the hip motion were taken into account for the sliding distance calculation. The contact pressure distribution was calculated on the basis of the Hertzian theory for the elastic contact of two bodies with non-conforming geometrical shapes. The wear factor was taken from hip simulator wear tests. The calculated WG is 67 x 10(-6) mm3 for a standard reference patient. The parametric model simulations show that WG increases linearly with the patient weight, femoral head diameter and surface roughness. It increases non-linearly to a maximum and decreases to an asymptotic value with increasing cup/head clearance and with cup isotropic elastic modulus. The cup orientation in the pelvis affects only slightly the total amount of WG whereas it is the dominant factor affecting the shape of the wear distribution. The iso-wear maps show paracentral patterns at low cup inclination angles and marginal patterns at higher inclination angles. The maximum wear depth is supero-posterior when the cup is in neutral alignment and supero-anterior at increasing anteversion angles. Complex patterns with a combination of paracentral and marginal wear were obtained at specific clearance values and cup orientations. The results of the simulations are discussed in relation to the wear distribution measured on the articular surface of 12 UHMWPE components retrieved from failed hip joint prostheses, after a period of in situ functioning.

Quantitative evaluation of the prosthetic head damage induced by microscopic third-body particles in total hip replacement.

The increase of the femoral head roughness in artificial hip joints is strongly influenced by the presence of abrasive particulate entrapped between the articulating surfaces. The aim of the present study is to evaluate the dependence of such damage on the geometry of the particles entrapped in the joint, with reference to the UHMWPE/chrome-cobalt coupling. Five chrome-cobalt femoral heads and their coupled UHMWPE acetabular cups, retrieved at revision surgery after a short period of in situ functioning, have been investigated for the occurrence of third-body damage. This was found on all the prosthetic heads, where the peak-to-valley height of the scratches, as derived from profilometry evaluations, ranged from 0.3-1.3 microm. The observed damage has been divided into four classes, related to the particle motion while being embedded into the polymer. Two kinds of particle morphology have been studied, spherical and prismatic, with size ranging from 5-50 microm. In order to provide an estimation of the damage induced by such particles, a finite element model of the third-body interaction was set up. The peak-to-valley height of the impression due to the particle indentation on the chrome-cobalt surface is assumed as an index of the induced damage. The calculated values range from 0.1-0.5 microm for spherical particles of size ranging from 10-40 microm. In the case of prismatic particles, the peak-to-valley height can reach 1.3 microm and depends both on the size and width of the particle's free corner, indenting the chrome-cobalt. As an example, a sharp-edged particle of size 30 microm can induce on the chrome-cobalt an impression with peak-to-valley height of 0.75 microm, when embedded into the polyethylene with a free edge of 5 microm facing the metallic surface. Negligible damage is induced, if a free edge of 7.5 microm is indenting the counterface. Our findings offer new support to the hypothesis that microscopic third-body particles are capable of causing increased roughening of the femoral head and provide a quantitative evaluation of the phenomenon.

Modeling of the Norwood circulation: effects of shunt size, vascular resistances, and heart rate.

Hypoplastic left heart syndrome is the most common lethal cardiac malformation of the newborn. Its treatment, apart from heart transplantation, is the Norwood operation. The initial procedure for this staged repair consists of reconstructing a circulation where a single outlet from the heart provides systemic perfusion and an interpositioning shunt contributes blood flow to the lungs. To better understand this unique physiology, a computational model of the Norwood circulation was constructed on the basis of compartmental analysis. Influences of shunt diameter, systemic and pulmonary vascular resistance, and heart rate on the cardiovascular dynamics and oxygenation were studied. Simulations showed that 1) larger shunts diverted an increased proportion of cardiac output to the lungs, away from systemic perfusion, resulting in poorer O2 delivery, 2) systemic vascular resistance exerted more effect on hemodynamics than pulmonary vascular resistance, 3) systemic arterial oxygenation was minimally influenced by heart rate changes, 4) there was a better correlation between venous O2 saturation and O2 delivery than between arterial O2 saturation and O2 delivery, and 5) a pulmonary-to-systemic blood flow ratio of 1 resulted in optimal O2 delivery in all physiological states and shunt sizes.

Experimental and computational approach for the evaluation of the biomechanical effects of dental bridge misfit.

Dental bridges supported by osseointegrated implants are commonly used to treat the partially or completely edentulous jaw. The bridges are manufactured in metal alloy using a sequence of technological steps which well match the requirement to get custom overstructures but does not guarantee geometrical and dimensional tolerances. Dentists often experience that a perfect fit of the bridge with the abutments is almost impossible to achieve. When a misfitting bridge is forced on the abutments, deformations may occur inducing a permanent preload at the fixture-bone interface and the greater the misfit the greater is the preload and the risk of implant failure. This work gives an evaluation of the biomechanical effects induced by a misfitting bridge when forced on two supporting dental implants. The strains induced in the bridge have been measured using two purposely designed and fabricated experimental devices allowing different types of misfit. FEM 3D models of the bridge and of the bridge anchored to the bone by implants have been developed. The former has been validated by simulating the same loading conditions as in the experimental tests and comparing the bridge strains. Both models have been used for the evaluation of the stress induced in the bridge and at the fixture-bone interface by bridge length errors. The results show that the method may help to estimate the stress distribution in the bridge and bone as a consequence of different dental bridge misfits.

The in-vivo wear performance of prosthetic femoral heads with titanium nitride coating.

This paper reports the study performed on four titanium nitride (TiN) coated prosthetic femoral heads collected at revision surgery together with patient data. Surface topology has been examined using Scanning Electron Microscopy (SEM) and elemental analysis of both coating and substrate have been evaluated using energy-dispersive X-ray spectrometry. Quantitative assessment of the surface topography is achieved using contacting profilometry. The average Ra roughness value is calculated at five different locations for each femoral head. The UHMWPE counterface worn volume has been measured directly on the acetabular components. TiN fretting and coating breakthrough occurred in two of the four components examined. In the damaged coating areas the surface profile is macroscopically saw-toothed with average tooth height 1.5 microm. The average Ra value is 0.02 microm on the undamaged surfaces and 0.37 microm on the damaged ones. Failure of the coating adhesion resulted in the release of TiN fragments and of metallic particulate from the substrate fretting corrosion and in the increase of the head surface roughness affecting counterface debris production. Our results suggest that TiN-coated titanium alloy femoral heads are inadequate in the task of resisting third body wear mechanisms in vivo.

Computational model of the fluid dynamics in systemic-to-pulmonary shunts.

A systemic-to-pulmonary shunt is a connection created between the systemic and pulmonary arterial circulations in order to improve pulmonary perfusion in children with congenital heart diseases. Knowledge of the relationship between pressure and flow in this new, surgically created, cardiovascular district may be helpful in the clinical management of these patients, whose survival is critically dependent on the blood flow distribution between the pulmonary and systemic circulations. In this study a group of three-dimensional computational models of the shunt have been investigated under steady-state and pulsatile conditions by means of a finite element analysis. The model is used to quantify the effects of shunt diameter (D), curvature, angle, and pulsatility on the pressure-flow (DeltaP-Q) relationship of the shunt. Size of the shunt is the main regulator of pressure-flow relationship. Innominate arterial diameter and angles of insertion have less influence. Curvature of the shunt results in lower pressure drops. Inertial effects can be neglected. The following simplified formulae are derived: DeltaP=(0. 097Q+0.521Q(2))/D(4) and DeltaP=(0.096Q+0.393Q(2))/D(4) for the different shunt geometries investigated (straight and curved shunts, respectively).

Use of mathematical model to predict hemodynamics in cavopulmonary anastomosis with persistent forward flow.

BACKGROUND: The bidirectional cavopulmonary anastomosis with additional pulmonary blood flow is used as a staged procedure or a definitive palliation of univentricular hearts. In this paper the flow competition occurring between the caval and the pulmonary flows is investigated. The hemodynamics in the superior vena cava and the blood flow distribution into the lungs, as well as the systemic arterial oxygen availability, are correlated with the severity of the right ventricle outflow tract obstruction and the pulmonary arteriolar resistance. MATERIALS AND METHODS: Computer models of the pre- and postoperative hemodynamics of univentricular hearts were developed. The effects of increasing severity of the right ventricle outflow tract obstruction, with a pulmonary arteriolar resistance ranging from 0.8 to 7.9 nonindexed Woods units, were simulated. RESULTS: The study indicates that the presence of an additional pulmonary blood flow from the native pulmonary artery may be beneficial. Since an excessive additional blood flow may cause central venous hypertension, its optimal value should be chosen according to the value of pulmonary arteriolar resistance. The model was utilized to simulate four clinical cases. CONCLUSIONS: The simulations show that the model can predict the postoperative hemodynamics and could therefore be usefully applied to predict quantitatively the effect of the native pulmonary blood flow following bidirectional cavopulmonary anastomosis.

An in vitro study on compensation of mismatch of screw versus cement-retained implant supported fixed prostheses.

In common practice a perfect fit of the prosthetic framework with the implant abutments is almost impossible to achieve. The mismatch, which is principally induced by the technological process adopted to manufacture the fixed prostheses, strains the framework thus generating constraint reactions. These are static forces that load the implant components and the bone at the implant-bone interface and may cause the bone remodelling. Depending on the magnitude of such forces, i.e. depending on the magnitude of the mismatch, the bone remodelling may lead to the loosening of the screws and of the implant-bone interface and hence cause implant failure. The present study shows an in vitro comparison of 3 different connecting abutments (standard, EsthetiCone and CerAdapt, Nobel Biocare AB, Göteborg, Sweden) with relation to the mechanical stresses induced by geometrical mismatches (technology induced errors). Two experimental devices were purposely realized and used to assess the ability of the different abutments to compensate errors. One was designed for translation errors and the other for rotation errors. The experimental apparatus set-up includes 2 freestanding implants supporting a prosthetic structure and the connecting abutments. The implants and the abutments were used as delivered by the manufacturer, while the prostheses were purposely realized and instrumented with strain gauges. The data obtained with the error devices do not give quantitative information on what happens in clinical applications where the implants are connected to living bone, which is a tissue much more deformable than the steel used for the error devices. Results allow direct comparisons of the behaviour of the different investigated abutments with respect to position errors. The CerAdapt system (cement retained ceramic abutments) showed the least strain in presence of translation errors. The standard system (screw retained abutments) showed the least strain in presence of rotation errors.

A method for the evaluation of the change in volume of retrieved acetabular cups.

Studies on retrieved hip prostheses are currently performed in order to assess the wear mechanisms and the overall wear rate of such artificial joints. Several reported studies on the survival of artificial hips have been based on the measurement of the amount of worn material directly on retrieved acetabular cups. The estimation of the change in volume, V, of the cup cavity is particularly difficult in the case of slight wear due to several factors of which the most critical is the lack of information on the unworn geometry of the cup. This paper presents a new measuring technique, which is described in detail and has been applied to estimate the wear of 65 acetabular cups harvested from revised hip arthroplasty. The coordinate data of several points on the articular surface are sampled using a coordinate measuring machine (CMM). The value of V is calculated mathematically from the measurements by using the hypothesis that the actual shape of both slightly worn and highly worn surfaces has a small departure from a truly spherical shape. The wear volume is estimated with reasonable accuracy mainly on those cups showing penetration depths greater than 0.2 mm, corresponding to an amount of wear greater than 100 mm3 for a 32 mm cup. The uncertainty in the results is estimated for each cup. The repeatability of the technique is studied for a case showing very slight wear. The advantages, disadvantages and limitations of the method are presented and discussed.

Modelling evaluation of the testing condition influence on the maximum stress induced in a hip prosthesis during ISO 7206 fatigue testing.

In vivo fatigue failure of hip prosthesis stems has been extensively reported in literature. The ISO 7206 international standard has been developed to assess the fatigue reliability of hip prostheses. It describes the fatigue testing apparatus and procedure and it is currently adopted by several testing laboratories throughout the world. In this work we evaluate the maximum stress in a titanium alloy commercial stem in different testing conditions, ranging within the standard specification, using the finite element method applied to a 3D model of the stem. The calculated maximum von Mises stress ranges from +4.5 to -1.5% (for different cement constraint levels) and from +6.7 to -6.8% (for different stem angular orientations) with respect to that calculated at the nominal testing conditions. The results suggest that the ISO 7206 testing specification will give experimental data of reasonable accuracy, with probably no more scatter than that found in typical specimen test results. This is particularly important in the case of components manufactured from materials showing a fatigue resistance highly sensitive to stress variations, such as the Ti6A14V alloy, for which a small increase of the maximum applied stress corresponds to a significant decrease of the statistical fatigue life.

Computational fluid dynamic simulations of cavopulmonary connections with an extracardiac lateral conduit.

Complex congenital heart defects due to the absence of a ventricular chamber can often be treated by the Fontan surgical procedure. The objective of this work was to quantify the haemodynamics in the Fontan operation (cavopulmonary connection) with extracardiac lateral conduit. Four different models based on the finite element method were constructed with different lengths of inferior anastomosis (range 18-25 mm) and inclinations of the conduit (33 and 47.5 degrees). Mass conservation and Navier-Stokes equations were solved by means of the FIDAP code, based on the finite element method. The left-to-right pulmonary flow ratio and percentage inferior caval blood to the left lung were the highest with the smallest anastomosis and highest inclination: 1.35 and 83.26%, respectively. Dissipated power percentage was higher with the largest anastomosis than with the smallest (19.4 vs 15.8%). It was concluded that, when performing a total cavopulmonary connection, an extracardiac lateral conduit: (i) diverts more flow to the left lung, and (ii) shows higher energy losses when compared with a connection with intra-atrial tunnel. This study could be useful to evaluate the incidence of pulmonary arteriovenous malformations.

Lumbar dura mater biomechanics: experimental characterization and scanning electron microscopy observations.

UNLABELLED: There is no consensus about the anatomical structure of human dura mater. In particular, the orientation of collagen fibers, which are responsible for biomechanical behavior, is still controversial. The aim of this work was to evaluate the mechanical properties and the microstructure of the lumbar dura mater. We performed experimental mechanical characterization in longitudinal and circumferential directions and a scanning electron microscopy observation of the tissue. Specimens of human dura mater were removed from the dorsal-lumbar region (T12-L4/L5) of six subjects at autopsy; specimens of bovine dorsal-lumbar dura mater were obtained from two animals at slaughter. Human and bovine tissues both exhibited stronger tensile strength and stiffness in the longitudinal than in the circumferential direction. Scanning electron microscopy observations of dura mater showed that the collagen fibers are mainly oriented in a longitudinal direction, which accounts for its stronger tensile strength in this direction. We conclude that dura mater has a different mechanical response in the two directions investigated because the fiber orientation is predominantly longitudinal. IMPLICATIONS: In this experimental work, we studied the structural and functional relationship of human lumbar dura mater. We performed mechanical tests and microscopic observations on dura mater samples. The results show that the dura mater is mainly composed of longitudinally oriented collagen fibers, which account for higher tissue resistance in this direction.

Computational fluid dynamic and magnetic resonance analyses of flow distribution between the lungs after total cavopulmonary connection.

Total cavopulmonary connection is a surgical procedure adopted to treat complex congenital malformations of the right heart. It consists basically in a connection of both venae cavae directly to the right pulmonary artery. In this paper a three-dimensional model of this connection is presented, which is based on in vivo measurements performed by means of magnetic resonance. The model was developed by means of computational fluid dynamics techniques, namely the finite element method. The aim of this study was to verify the capability of such a model to predict the distribution of the blood flow into the pulmonary arteries, by comparison with in vivo velocity measurements. Different simulations were performed on a single clinical case to test the sensitivity of the model to different boundary conditions, in terms of inlet velocity profiles as well as outlet pressure levels. Results showed that the flow distribution between the lungs is slightly affected by the shape of inlet velocity profiles, whereas it is influenced by different pressure levels to a greater extent.