A Porous Media Approach to Finite Deformation Behaviour in Soft Tissues.

The present work presents a porous medium formulation for the biomechanical analysis of soft tissues. An updated Lagrangian approach is developed to study the coupled effects of low speed flows of fluid phases, in partially or fully saturated conditions, and the finite deformation occurring in the solid matrix. The procedure developed allows both for the evaluation of coupled geometric and material non-linearities. The main theoretical and computational aspects of this multiphase formulation are discussed. The finite element method is used for the numerical solution of the resulting coupled system of equations. A reference case is reported with regard to healthy and degenerative phases of intervertebral segment. Results reported allow for a detailed interpretation of the formulation reliability, also by comparison with existing experimental data. In particular, the role played by the fluid on the load carrying mechanism is pointed out, thus stressing the importance of a multiphase approach to the overall behaviour of the spinal motion segment in time.

Investigation of the integration process of dental implants by means of a numerical analysis of dynamic response.

OBJECTIVES: The aim of this work was to present a preliminary numerical analysis of the integration process of dental implants using a finite element simulation of the dynamic response following impulse excitation. Assessment of the osseointegration process has been previously examined using a numerical approach by calculating the natural frequency of a cantilever attached to the implant. The methodology adopted in this work allows a direct measurement of the implant response following impulse loading and avoids the addition of a bulky cantilever set-up. METHODS: The geometric configuration was obtained by averaging the coordinate data from tomographic scans of 14 mandibles. The materials properties were approximated from experimental analysis performed on trabecular and cortical bone tissue. A load was applied to the top of the implant in one direction resulting in an initial displacement. The implant was then freed and allowed to vibrate over approximately 10 cycles. Three fixity conditions were assumed by changing the properties of the surrounding bone ranging from full integration to a poorly integrated implant typical of the situation during bone healing following surgery. The results of the three fixity conditions were compared by calculating the fundamental displacement amplitudes and frequencies of the vibrating impact. RESULTS: The calculated results indicated that the implant vibrated at a predominant frequency when partially integrated with a displacement principally in the direction of the applied impulse. However, when the implant was fully integrated a more complex vibration pattern ensued, suggesting the superposition of two or more fundamentals. SIGNIFICANCE: Attention has been paid to the formulation of the numerical model for validation purposes as well as a reliable reference for the optimum interpretation of the experimental data. In this way it was possible to establish a simulation procedure to investigate the response of the tissues surrounding the implant and their properties at different stages of healing. It should be pointed out that the numerical procedures represented a valid preliminary approach to the problem and were capable of indicating a guide to the optimum design of the experimental apparatus for measurement of displacement and frequency in vivo.

Adaptive finite-element approach for analysis of bone/prosthesis interaction.

The study uses the finite-element method to analyse the stress field in a perfectly bonded hip prosthesis arising from loading through body weight. Special attention is paid to the accuracy of the numerical analysis, and adaptive mesh refinement is introduced to reduce the discretisation error. The finite-element procedure developed is especially well suited to analyse the behaviour of a bonded interface as it is capable of calculating accurately the stress at the nodal positions while satisfying the natural discontinuity in the stress field at this location. An orthotropic material model is used for the representation of the behaviour of the bone, and an axisymmetric geometry with non-symmetrical loading is adopted for the analysis. The results demonstrate the usefulness of adaptive mesh refinement and the significance of adopting anisotropic material modelling in the context of tissue/prosthesis interaction.

The mechanical behaviour of bony endplate and annulus in prolapsed disc configuration.

The mechanical response of intervertebral joints is deeply influenced by disc degeneration. The phenomenon is expressed in terms of variations in the biomechanical properties of the material, whose compressibility characteristics change because of the liquid content loss in the tissue and, what is even more important, to prolapse. In this work, the problem is investigated by means of a computational mechanics approach; a coupled material and geometric non-linear model is developed, representing vertebra, annulus and nucleus submitted to an axial load. A transversely isotropic law is assumed for cortical bone in the vertebral body and an isotropic law for the cancellous portion; a hyperelastic formulation is assumed for the disc, allowing effective interpretation of the mechanical characteristics of degeneration. The results obtained are reported with regard to bony endplate and annulus behaviour; interaction phenomena between bony endplate and nucleus are emphasized.

A numerical approach to the biomechanical analysis of bone fracture healing.

Investigations are reported in the literature, by means of experimental, analytical and numerical methods, concerning the biomechanical properties of bone. However, the evolutionary phenomena of bone fracture healing does not have a large reference literature. This work investigates and describes the behaviour of inclined human femur fractures with external fixation up to complete healing. A numerical formulation based on the finite element method has been adopted. Geometric configuration is defined using data from a magnetic resonance process applied to a femur in vivo. A three dimensional model has been developed by adopting an orthotropic material law for cortical bone and an isotropic law for the fracture gap zone. Stress and strain responses of the bone and fixation device are investigated with reference to the evolutionary behaviour of the healing tissue.

A review of the biomechanical properties of bone as a material.

Experimental studies on bone all reveal important difficulties in data interpretation. This paper proposes an analysis of experimental studies performed so far, with particular attention to the anisotropic characteristics of bone, its behaviour in the post-elastic phase, and its dependence on viscoelastic phenomena. Mechanical properties are also correlated with variables such as moisture, deformation rate during testing, density, variations in different regions of the bone, and the most relevant strength criteria are recalled. The investigation performed is also intended to provide an evaluation of the degree of refinement of biomechanical experimental data for use in a numerical approach to bone mechanics.