Large-sliding contact elements accurately predict levels of bone-implant micromotion relevant to osseointegration

Primary stability is recognised as an important determinant in the aseptic loosening failure process of cementless implants. An accurate evaluation of the bone–implant relative micromotion is becoming important both in pre-clinical and clinical studies. If the biological threshold for micro-movements is in the range 100–200 μm then, in order to be discriminative, any method used to evaluate the primary stability should have an accuracy of 10–20 μm or better. Additionally, such method should also be able to report the relative micromotion at each point of the interface. None of the available experimental methods satisfies both requirements. Aim of the present study is to verify if any of the current finite element modelling techniques is sufficiently accurate in predicting the primary stability of a cementless prosthesis to be used to decide whether the micromotion may or may not jeopardise the implant osseointegration. The primary stability of an anatomic cementless stem, as measured in vitro, was used as a benchmark problem to comparatively evaluate different contact modelling techniques. Frictionless contact, frictional contact and press-fitted frictional contact conditions were modelled using alternatively node-to-node, node-to-face and face-to-face contact elements. The model based on face-to-face contact elements accounting for frictional contact and initial press-fit was able to predict the micromotion measured experimentally with an average (RMS) error of 10 μm and a peak error of 14 μm. All the other models presented errors higher than 20 μm assumed in the present study as an accuracy threshold.     

A new material mapping procedure for quantitative computed tomography-based, continuum finite element analyses of the vertebra

Inaccuracies in the estimation of material properties and errors in the assignment of these properties into finite element models limit the reliability, accuracy, and precision of quantitative computed tomography (QCT)-based finite element analyses of the vertebra. In this work, a new mesh-independent, material mapping procedure was developed to improve the quality of predictions of vertebral mechanical behavior from QCT-based finite element models. In this procedure, an intermediate step, called the material block model, was introduced to determine the distribution of material properties based on bone mineral density, and these properties were then mapped onto the finite element mesh. A sensitivity study was first conducted on a calibration phantom to understand the influence of the size of the material blocks on the computed bone mineral density. It was observed that varying the material block size produced only marginal changes in the predictions of mineral density. Finite element (FE) analyses were then conducted on a square column-shaped region of the vertebra and also on the entire vertebra in order to study the effect of material block size on the FE-derived outcomes. The predicted values of stiffness for the column and the vertebra decreased with decreasing block size. When these results were compared to those of a mesh convergence analysis, it was found that the influence of element size on vertebral stiffness was less than that of the material block size. This mapping procedure allows the material properties in a finite element study to be determined based on the block size required for an accurate representation of the material field, while the size of the finite elements can be selected independently and based on the required numerical accuracy of the finite element solution. The mesh-independent, material mapping procedure developed in this study could be particularly helpful in improving the accuracy of finite element analyses of vertebroplasty and spine metastases, as these analyses typically require mesh refinement at the interfaces between distinct materials. Moreover, the mapping procedure is not specific to the vertebra and could thus be applied to many other anatomic sites.     

Contact analysis of a posterior cervical spine plate using a three-dimensional canine finite element model

The development of a three-dimensional finite element model of a posteriorly plated canine cervical spine (C3-C6) including contact nonlinearities is described. The model was created from axial CT scans and the material properties were derived from the literature. The model demonstrated sufficient accuracy from the results of a mesh convergence test. Significant steps were taken toward establishing model validation by comparison of plate surface strains with a posteriorly plated canine cervical spine under three-point bending. This model was developed to better characterize the contact pressures at the various interfaces under average physiologic canine loading. The analysis showed that the screw-plate interfaces had the highest values of all the mechanical parameters evaluated.     

Some effects of global joint morphology on local stress aberrations near imprecisely reduced intra-articular fractures

Cadaver models of contact pressure aberration near displaced intra-articular fractures complement clinical experience, but inter-specimen variability often complicates interpretation of in vitro data. A contact finite element formulation is here used to study juxta-articular stress distributions in a plane strain model of tibial plateau step-off incongruity. Attention is focused on the influence of global morphologic parameters: intact joint surface curvatures, cartilage thickness, and cartilage stiffness. The computed stress distributions agreed well with experimental recordings for a typical 3 mm incongruity in an otherwise normal joint. Both decreased cartilage thickness and increased cartilage modulus led to elevations in the peak local contact stress, and to concentration of contact stress near the edge of the step-off incongruity. Similar effects were seen when reduction of global joint congruency was modelled by decreasing the concavity of the tibial plateau. While the observed degree of coupling between global morphology and local stress aberration was by no means negligible, the sensitivity of the stresses to variations in individual parameters was relatively mild. This suggests that the finite element results will be useful for experimental data interpretation.     

Assessment of the effect of mesh density on the material property discretisation within QCT based FE models: a practical example using the implanted proximal tibia

A three-dimensional, quantitative computed tomography based finite element model of a proximal implanted tibia was analysed in order to assess the effect of mesh density on material property discretisation and the resulting influence on the predicted stress distribution. The mesh was refined on the contact surfaces (matched meshes) with element sizes of 3, 2, 1.4, 1 and 0.8 mm. The same loading conditions were used in all models (bi-condylar load: 60% medial, 40% lateral). Significant variations were observed in the modulus distributions between the coarsest and finest mesh densities. Poor discretisation of the material properties also resulted in poor correlations of the stresses and risk ratios between the coarsest and finest meshes. Little difference in Young's modulus, von Mises stress and risk ratio distributions were observed between the three finest models; hence, it was concluded that for this particular case an element size of 1.4 mm on the contact surfaces was enough to properly describe the stiffness, stress and risk ratio distributions within the bone. Poor convergence of the material property distribution occurred when the element size was significantly larger than the pixel size of the source CT data. It was concluded that unless there is convergence in the Young's modulus distribution, convergence of the stress field or of other parameters of interest will not occur either.     

Quantitative computed tomography-based finite element models of the human lumbar vertebral body: effect of element size on stiffness,damage, and fracture strength predictions

This study investigated the numerical convergence characteristics of specimen-specific "voxel-based" finite element models of 14 excised human cadaveric lumbar vertebral bodies (age: 37-87; M = 6, F = 8) that were generated automatically from clinical-type CT scans. With eventual clinical applications in mind, the ability of the model stiffness to predict the experimentally measured compressive fracture strength of the vertebral bodies was also assessed. The stiffness of "low"-resolution models (3 x 3 x 3 mm element size) was on average only 4% greater (p = 0.03) than for "high"-resolution models (1 x 1 x 1.5 mm) despite interspecimen variations that varied over four-fold. Damage predictions using low- vs high-resolution models were significantly different (p = 0.01) at loads corresponding to an overall strain of 0.5%. Both the high (r2 = 0.94) and low (r2 = 0.92) resolution model stiffness values were highly correlated with the experimentally measured ultimate strength values. Because vertebral stiffness variations in the population are much greater than those that arise from differences in voxel size, these results indicate that imaging resolution is not critical in cross-sectional studies of this parameter. However, longitudinal studies that seek to track more subtle changes in stiffness over time should account for the small but highly significant effects of voxel size. These results also demonstrate that an automated voxel-based finite element modeling technique may provide an excellent noninvasive assessment of vertebral strength.     

Discretization error when using finite element models: Analysis and evaluation of an underestimated problem

Mesh convergence tests are often insufficiently performed in finite element analyses. There are many parameters which may have an effect on the mesh convergence behavior. The aim of this study was to identify the influence of different parameters on the mesh convergence behavior.For this purpose we used a simplified axis-symmetrical model of a single pedicle screw flank with surrounding bone to simulate a pull-out test. In parameter studies, the flank radii and the contact conditions at the bone–screw interface were varied. These parameter studies were carried out using an implicit and explicit solver. Thereby, the convergence criteria and the number of the substeps for the implicit nonlinear iteration process as well as the velocity and the material density for the explicit approach were considered.The mesh convergence behavior was influenced by varying the flank radii and the contact conditions. The implicit calculations led to a reaction force, which converged rapidly to a certain value with increasing mesh density, whereas the maximum von-Mises stress showed substantial convergence problems. The number of substeps and the convergence criteria of the iteration process strongly influenced the implicit solutions. In contrast, the maximum von-Mises stresses resulting from explicit calculations converged to a certain value after only a few refinement steps. Different pull-out velocities substantially affected the mesh convergence behavior, while the material density showed only a negligible influence.The results indicated the need to perform an appropriate mesh convergence test when using finite element methods. We were able to show that different parameters strongly influence the mesh convergence behavior and we demonstrated that convergence tests do not always lead to a satisfactory or acceptable solution.     

Assessment of the Effect of Mesh Density on the Material Property Discretisation Within QCT Based FE Models: A Practical Example Using the Implanted Proximal Tibia

A three-dimensional, quantitative computed tomography based finite element model of a proximal implanted tibia was analysed in order to assess the effect of mesh density on material property discretisation and the resulting influence on the predicted stress distribution. The mesh was refined on the contact surfaces (matched meshes) with element sizes of 3, 2, 1.4, 1 and 0.8 mm. The same loading conditions were used in all models (bi-condylar load: 60% medial, 40% lateral). Significant variations were observed in the modulus distributions between the coarsest and finest mesh densities. Poor discretisation of the material properties also resulted in poor correlations of the stresses and risk ratios between the coarsest and finest meshes. Little difference in Young's modulus, von Mises stress and risk ratio distributions were observed between the three finest models; hence, it was concluded that for this particular case an element size of 1.4 mm on the contact surfaces was enough to properly describe the stiffness, stress and risk ratio distributions within the bone. Poor convergence of the material property distribution occurred when the element size was significantly larger than the pixel size of the source CT data. It was concluded that unless there is convergence in the Young's modulus distribution, convergence of the stress field or of other parameters of interest will not occur either.     

Development and validation of a weight-bearing finite element model for total knee replacement

Total knee arthroplasty (TKA) is a successful procedure for osteoarthritis. However, some patients (19%) do have pain after surgery. A finite element model was developed based on boundary conditions of a knee rig. A 3D-model of an anatomical full leg was generated from magnetic resonance image data and a total knee prosthesis was implanted without patella resurfacing. In the finite element model, a restarting procedure was programmed in order to hold the ground reaction force constant with an adapted quadriceps muscle force during a squat from 20° to 105° of flexion. Knee rig experimental data were used to validate the numerical model in the patellofemoral and femorotibial joint. Furthermore, sensitivity analyses of Young’s modulus of the patella cartilage, posterior cruciate ligament (PCL) stiffness, and patella tendon origin were performed. Pearson’s correlations for retropatellar contact area, pressure, patella flexion, and femorotibial ap-movement were near to 1. Lowest root mean square error for retropatellar pressure, patella flexion, and femorotibial ap-movement were found for the baseline model setup with Young’s modulus of 5 MPa for patella cartilage, a downscaled PCL stiffness of 25% compared to the literature given value and an anatomical origin of the patella tendon. The results of the conducted finite element model are comparable with the experimental results. Therefore, the finite element model developed in this study can be used for further clinical investigations and will help to better understand the clinical aspects after TKA with an unresurfaced patella.     

Cartilage thickness distribution affects computational model predictions of cervical spine facet contact parameters

With motion-sparing disk replacement implants gaining popularity as an alternative to anterior cervical discectomy and fusion (ACDF) for the treatment of certain spinal degenerative disorders, recent laboratory investigations have studied the effects of disk replacement and implant design on spinal kinematics and kinetics. Particularly relevant to cervical disk replacement implant design are any postoperative changes in solid stresses or contact conditions in the articular cartilage of the posterior facets, which are hypothesized to lead to adjacent-level degeneration. Such changes are commonly investigated using finite element methods, but significant simplification of the articular geometry is generally employed. The impact of such geometric representations has not been thoroughly investigated. In order to assess the effects of different models of cartilage geometry on load transfer and contact pressures in the lower cervical spine, a finite element model was generated using cadaver-based computed tomography imagery. Mesh resolution was varied in order to establish model convergence, and cadaveric testing was undertaken to validate model predictions. The validated model was altered to include four different geometric representations of the articular cartilage. Model predictions indicate that the two most common representations of articular cartilage geometry result in significant reductions in the predictive accuracy of the models. The two anatomically based geometric models exhibited less computational artifact, and relatively minor differences between them indicate that contact condition predictions of spatially varying thickness models are robust to anatomic variations in cartilage thickness and articular curvature. The results of this work indicate that finite element modeling efforts in the lower cervical spine should include anatomically based and spatially varying articular cartilage thickness models. Failure to do so may result in loss of fidelity of model predictions relevant to investigations of physiological import.     

Micro-mechanical modeling of the cement-bone interface: the effect of friction, morphology and material properties on the micromechanical response

In order to gain insight into the micro-mechanical behavior of the cement-bone interface, the effect of parametric variations of frictional, morphological and material properties on the mechanical response of the cement-bone interface were analyzed using a finite element approach. Finite element models of a cement-bone interface specimen were created from micro-computed tomography data of a physical specimen that was sectioned from an in vitro cemented total hip arthroplasty. In five models the friction coefficient was varied (mu=0.0; 0.3; 0.7; 1.0 and 3.0), while in one model an ideally bonded interface was assumed. In two models cement interface gaps and an optimal cement penetration were simulated. Finally, the effect of bone cement stiffness variations was simulated (2.0 and 2.5 GPa, relative to the default 3.0 GPa). All models were loaded for a cycle of fully reversible tension-compression. From the simulated stress-displacement curves the interface deformation, stiffness and hysteresis were calculated. The results indicate that in the current model the mechanical properties of the cement-bone interface were caused by frictional phenomena at the shape-closed interlock rather than by adhesive properties of the cement. Our findings furthermore show that in our model maximizing cement penetration improved the micromechanical response of the cement-bone interface stiffness, while interface gaps had a detrimental effect. Relative to the frictional and morphological variations, variations in the cement stiffness had only a modest effect on the micro-mechanical behavior of the cement-bone interface. The current study provides information that may help to better understand the load-transfer mechanisms taking place at the cement-bone interface.     

HR-pQCT-based homogenised finite element models provide quantitative predictions of experimental vertebral body stiffness and strength with the same accuracy as μFE models

This study validated two different high-resolution peripheral quantitative computer tomography (HR-pQCT)-based finite element (FE) approaches, enhanced homogenised continuum-level (hFE) and micro-finite element (μFE) models, by comparing them with compression test results of vertebral body sections. Thirty-five vertebral body sections were prepared by removing endplates and posterior elements, scanned with HR-pQCT and tested in compression up to failure. Linear hFE and μFE models were created from segmented and grey-level CT images, and apparent model stiffness values were compared with experimental stiffness as well as strength results. Experimental and numerical apparent elastic properties based on grey-level/segmented CT images (N=35) correlated well for μFE (r2=0.748/0.842) and hFE models (r2=0.741/0.864). Vertebral section stiffness values from the linear μFE/hFE models estimated experimental ultimate apparent strength very well (r2=0.920/0.927). Calibrated hFE models were able to predict quantitatively apparent stiffness with the same accuracy as μFE models. However, hFE models needed no back-calculation of a tissue modulus or any kind of fitting and were computationally much cheaper.     

Assessment of factors influencing finite element vertebral model predictions

This study aimed to establish model construction and configuration procedures for future vertebral finite element analysis by studying convergence, sensitivity, and accuracy behaviors of semiautomatically generated models and comparing the results with manually generated models. During a previous study, six porcine vertebral bodies were imaged using a microcomputed tomography scanner and tested in axial compression to establish their stiffness and failure strength. Finite element models were built using a manual meshing method. In this study, the experimental agreement of those models was compared with that of semiautomatically generated models of the same six vertebrae. Both manually and semiautomatically generated models were assigned gray-scale-based, element-specific material properties. The convergence of the semiautomatically generated models was analyzed for the complete models along with material property and architecture control cases. A sensitivity study was also undertaken to test the reaction of the models to changes in material property values, architecture, and boundary conditions. In control cases, the element-specific material properties reduce the convergence of the models in comparison to homogeneous models. However, the full vertebral models showed strong convergence characteristics. The sensitivity study revealed a significant reaction to changes in architecture, boundary conditions, and load position, while the sensitivity to changes in material property values was proportional. The semiautomatically generated models produced stiffness and strength predictions of similar accuracy to the manually generated models with much shorter image segmentation and meshing times. Semiautomatic methods can provide a more rapid alternative to manual mesh generation techniques and produce vertebral models of similar accuracy. The representation of the boundary conditions, load position, and surrounding environment is crucial to the accurate prediction of the vertebral response. At present, an element size of 2x2x2 mm(3) appears sufficient since the error at this size is dominated by factors, such as the load position, which will not be improved by increasing the mesh resolution. Higher resolution meshes may be appropriate in the future as models are made more sophisticated and computational processing time is reduced.     

Finite element implementation of a generalized Fung-elastic constitutive model for planar soft tissues

Numerical simulations of the anisotropic mechanical properties of soft tissues and tissue-derived biomaterials using accurate constitutive models remain an important and challenging research area in biomechanics. While most constitutive modeling efforts have focused on the characterization of experimental data, only limited studies are available on the feasibility of utilizing those models in complex computational applications. An example is the widely utilized exponential constitutive model proposed by Fung. Although present in the biomechanics literature for several decades, implementation of this model into finite element (FE) simulations has been limited. A major reason for limited numerical implementations are problems associated with inherent numerical instability and convergence. To address this issue, we developed and applied two restrictions for a generalized Fung-elastic constitutive model necessary to achieve numerical stability. These are (1) convexity of the strain energy function, and (2) the condition number of material stiffness matrix set lower than a prescribed value. These constraints were implemented in the nonlinear regression used for constitutive model parameter estimation to the experimental biaxial mechanical data. We then implemented the generalized Fung-elastic model into a commercial FE code (ABAQUS, Pawtucket, RI, USA). Single element and multi-element planar biaxial test simulations were conducted to verify the accuracy and robustness of the implementation. Results indicated that numerical convergence and accurate FE implementation were consistently obtained. The present study thus presents an integrated framework for accurate and robust implementation of pseudo-elastic constitutive models for planar soft tissues. Moreover, since our approach is formulated within a general FE code, it can be straightforwardly adopted across multiple software platforms.     

How the stiffness of meniscal attachments and meniscal material properties affect tibio-femoral contact pressure computed using a validated finite element model of the human knee joint

In an effort to prevent degeneration of articular cartilage associated with meniscectomies, both meniscal allografts and synthetic replacements are subjects of current interest and investigation. The objectives of the current study were to (1) determine whether a transversely isotropic, linearly elastic, homogeneous material model of the meniscal tissue is necessary to achieve a normal contact pressure distribution on the tibial plateau, (2) determine which material and boundary condition (attachments) parameters affect the contact pressure distribution most strongly, and (3) set tolerances on these parameters to restore the contact pressure distribution to within a specified error. To satisfy these objectives, a finite element model of the tibio-femoral joint of a human cadaveric knee (including both menisci) was used to study the contact pressure distribution on the tibial plateau. To validate the model, the contact pressure distribution on the tibial plateau was measured experimentally in the same knee used to create the model. Within physiologically reasonable bounds on five material parameters and four attachment parameters associated with a meniscal replacement, an optimization was performed under 1200 N of compressive load on the set of nine parameters to minimize the difference between the experimental and model results. The error between the experimental and model contact variables was minimized to 5.4%. The contact pressure distribution of the tibial plateau was sensitive to the circumferential modulus, axial/radial modulus, and horn stiffness, but relatively insensitive to the remaining six parameters. Consequently, both the circumferential and axial/radial moduli are important determinants of the contact pressure distribution, and hence should be matched in the design and/or selection of meniscal replacements. In addition, during surgical implantation of a meniscal replacement, the horns should be attached with high stiffness bone plugs, and the attachments of the transverse ligament and deep medial collateral ligament should be restored to minimize changes in the contact pressure distribution, and thereby possibly prevent the degradation of articular cartilage.     

The Penn State Safety Floor: Part II--Reduction of fall-related peak impact forces on the femur.

The goal of this study was to develop and validate a finite element model (FEM) for use in the design of a flooring system that would provide a stable walking surface during normal locomotion but would also deform elastically under higher loads, such as those resulting from falls. The new flooring system is designed to reduce the peak force on the femoral neck during a lateral fall onto the hip. The new flooring system is passive in nature and exhibits two distinct stiffnesses. During normal activities, the floor remains essentially rigid. Upon impact, the floor collapses and becomes significantly softer. The flooring system consists of a multitude of columns supporting a continuous walking surface. The columns were designed to remain stiff up to a specific load and, after exceeding this load, to deform elastically. The flooring returns to its original shape after impact. Part I of this study presented finite element and experimental results demonstrating that the floor deflection during normal walking remained less than 2 mm. To facilitate the floor's development further, a nonlinear finite element model simulating the transient-impact response of a human hip against various floor configurations was developed. Nonlinearities included in the finite element models were: changing topology of deformable-body-to-deformable-body contact, snap-through buckling, soft tissue stiffness and damping, and large deformations. Experimental models developed for validating the finite element model included an anthropomorphic hip, an impact delivery mechanism, a data collection system, and four hand-fabricated floor tiles. The finite element model discussed in this study is shown to capture experimentally observed trends in peak femoral neck force reduction as a function of flooring design parameters. This study also indicates that a floor can be designed that deflects minimally during walking and reduces the peak force on the femoral neck during a fall-related impact by 15.2 percent.     

An anisotropic elastic-viscoplastic damage model for bone tissue

A new anisotropic elastic-viscoplastic damage constitutive model for bone is proposed using an eccentric elliptical yield criterion and nonlinear isotropic hardening. A micromechanics-based multiscale homogenization scheme proposed by Reisinger et al. is used to obtain the effective elastic properties of lamellar bone. The dissipative process in bone is modeled as viscoplastic deformation coupled to damage. The model is based on an orthotropic ecuntric elliptical criterion in stress space. In order to simplify material identification, an eccentric elliptical isotropic yield surface was defined in strain space, which is transformed to a stress-based criterion by means of the damaged compliance tensor. Viscoplasticity is implemented by means of the continuous Perzyna formulation. Damage is modeled by a scalar function of the accumulated plastic strain ${D(\kappa)}$ , reducing all element s of the stiffness matrix. A polynomial flow rule is proposed in order to capture the rate-dependent post-yield behavior of lamellar bone. A numerical algorithm to perform the back projection on the rate-dependent yield surface has been developed and implemented in the commercial finite element solver Abaqus/Standard as a user subroutine UMAT. A consistent tangent operator has been derived and implemented in order to ensure quadratic convergence. Correct implementation of the algorithm, convergence, and accuracy of the tangent operator was tested by means of strain- and stress-based single element tests. A finite element simulation of nano- indentation in lamellar bone was finally performed in order to show the abilities of the newly developed constitutive model.     

A visco-hyperelastic-damage constitutive model for the analysis of the biomechanical response of the periodontal ligament

The periodontal ligament (PDL), as other soft biological tissues, shows a strongly non-linear and time-dependent mechanical response and can undergo large strains under physiological loads. Therefore, the characterization of the mechanical behavior of soft tissues entails the definition of constitutive models capable of accounting for geometric and material non-linearity. The microstructural arrangement determines specific anisotropic properties. A hyperelastic anisotropic formulation is adopted as the basis for the development of constitutive models for the PDL and properly arranged for investigating the viscous and damage phenomena as well to interpret significant aspects pertaining to ordinary and degenerative conditions. Visco-hyperelastic models are used to analyze the time-dependent mechanical response, while elasto-damage models account for the stiffness and strength decrease that can develop under significant loading or degenerative conditions. Experimental testing points out that damage response is affected by the strain rate associated with loading, showing a decrease in the damage limits as the strain rate increases. These phenomena can be investigated by means of a model capable of accounting for damage phenomena in relation to viscous effects. The visco-hyperelastic-damage model developed is defined on the basis of a Helmholtz free energy function depending on the strain-damage history. In particular, a specific damage criterion is formulated in order to evaluate the influence of the strain rate on damage. The model can be implemented in a general purpose finite element code. The accuracy of the formulation is evaluated by using results of experimental tests performed on animal model, accounting for different strain rates and for strain states capable of inducing damage phenomena. The comparison shows a good agreement between numerical results and experimental data.     

Automatic generation of accurate subject-specific bone finite element models to be used in clinical studies

Most of the finite element models of bones used in orthopaedic biomechanics research are based on generic anatomies. However, in many cases it would be useful to generate from CT data a separate finite element model for each subject of a study group. In a recent study a hexahedral mesh generator based on a grid projection algorithm was found very effective in terms of accuracy and automation. However, so far the use of this method has been documented only on data collected in vitro and only for long bones. The present study was aimed at verifying if this method represents a procedure for the generation of finite element models of human bones from data collected in vivo, robust, accurate, automatic and general enough to be used in clinical studies. Robustness, automation and numerical accuracy of the proposed method were assessed on five femoral CT data sets of patients affected by various pathologies. The generality of the method was verified by processing a femur, an ileum, a phalanx, a proximal femur reconstruction, and the micro-CT of a small sample of spongy bone. The method was found robust enough to cope with the variability of the five femurs, producing meshes with a numerical accuracy and a computational weight comparable to those found in vitro. Even when the method was used to process the other bones the levels of mesh conditioning remained within acceptable limits. Thus, it may be concluded that the method presents a generality sufficient to cope with almost any orthopaedic application.