A nonlinear homogenized finite element analysis of the primary stability of the bone–implant interface

Stability of an implant is defined by its ability to undergo physiological loading–unloading cycles without showing excessive tissue damage and micromotions at the interface. Distinction is usually made between the immediate primary stability and the long-term, secondary stability resulting from the biological healing process. The aim of this research is to numerically investigate the effect of initial implantation press-fit, bone yielding, densification and friction at the interface on the primary stability of a simple bone–implant system subjected to loading–unloading cycles. In order to achieve this goal, human trabecular bone was modeled as a continuous, elasto-plastic tissue with damage and densification, which material constants depend on bone volume fraction and fabric. Implantation press-fit related damage in the bone was simulated by expanding the drilled hole to the outer contour of the implant. The bone–implant interface was then modeled with unilateral contact with friction. The implant was modeled as a rigid body and was subjected to increasing off-axis loading cycles. This modeling approach is able to capture the experimentally observed primary stability in terms of initial stiffness, ultimate force and progression of damage. In addition, it is able to quantify the micromotions around the implant relevant for bone healing and osseointegration. In conclusion, the computationally efficient modeling approach used in this study provides a realistic structural response of the bone–implant interface and represents a powerful tool to explore implant design, implantation press-fit and the resulting risk of implant failure under physiological loading.     

A bone preservation protocol that enables evaluation of osseointegration of implants with micro- and nanotextured surfaces

Development of surface treatments has enabled secure attachment of dental implants in less than 1 month. Consequently, it is necessary to characterize accurately the osseointegration of the implant surface in the region of the bone-implant contact (BIC). We developed a method for sample preparation that preserves both bone and BIC to permit analysis of the contact interface. We prepared eight nanotextured implants and implanted them in rabbit tibias. After healing for 30 days, outcomes were analyzed using both our bone preservation protocol and routine decalcification followed by preparation of histological sections stained by hematoxylin and eosin (H & E). Pull-out tests for implant osseointegration were performed after healing. Non-implanted samples of rabbit mandible were used as a control for assessing organic and mineralized bone characteristics and bone structure. Our bone preservation protocol enabled evaluation of many of the same bone characteristics as histological sections stained with H & E. Our protocol enables analysis of implant samples, implant surfaces and osseointegration without risk of BIC damage.     

氧化钛纳米管对即刻种植骨结合影响的实验研究

目的:评价氧化钛纳米管对犬即刻种植骨结合效果的影响。方法:犬拔牙后即刻将光滑表面(对照组)和氧化钛纳米管表面(实验组)种植体植入拔牙窝内,于12周后处死取材,进行显微CT扫描、组织学染色分析以及生物力学检测。结果:扫描电镜显示经过阳极氧化后,钛表面形成了直径为30-80纳米的纳米管状结构;12周后,显微CT扫描结果提示实验组骨体积分数(BV/TV)、骨小梁数目(Tb.N)、骨小梁厚度(Tb.Th)均显著高于对照组,骨小梁间隙(Tb.Sp)显著低于对照组,差异具有统计学意义(P0.05)。术后12周,实验组与对照组骨结合率分别为49.35±11.76%、31.79±13.07%,最大拔出力分别为105.28±27.87N、79.23±20.46N,实验组均显著高于对照组,差异均具有统计学意义(P0.05)。结论:氧化钛纳米管表面有利于促进即刻种植后骨结合的效果。     

The effect of surface roughness on the stress adaptation of trabecular architecture around a cylindrical implant

The effect of implant-bone bonding and the effect of implant surface roughness on bone remodeling near the bone-implant interface were studied by using a surface remodeling theory and the boundary element method. The study has shown that implant attachment plays an important role in bone remodeling near the implant. It has been observed in animal experiments and in clinical situations that the remodeled trabecular bone architecture around a cylindrical implant could vary, on one hand, from a hub surrounding the implant with a set of external spokes to, on the other hand, a hubless situation in which a set of spokes attach directly to the implant. It is shown here that the difference in these structures may be attributed to differences in implant attachment. The results show that the bone with perfect bonding or roller boundary condition without a gap remodeled to a hubless spoke trabecular bone architecture. On the other hand, the roller boundary condition with a specified gap yielded a spoke trabecular architecture with a hub or ring surrounding the implant. These quantitative results mirror the experimental and clinical observations. It is concluded that the hub is a consequence of the gap and not a consequence of the lack of friction between the implant and the bone.     

Artificial composite bone as a model of human trabecular bone: the implant-bone interface

The use of artificial bones in implant testing has become popular due to their low variability and ready availability. However, friction coefficients, which are critical to load transfer in uncemented implants, have rarely been compared between human and artificial bone, particularly for wet and dry conditions. In this study, the static and dynamic friction coefficients for four commercially used titanium surfaces (polished, Al(2)O(3) blasted, plasma sprayed, beaded) acting on the trabecular component of artificial bones (Sawbones) were compared to those for human trabecular bone. Artificial bones were tested in dry and wet conditions and normal interface stress was varied (0.25, 0.5, 1.0MPa). Friction coefficients were mostly lower for artificial bones than real bone. In particular, static friction coefficients for the dry polished surface were 20% of those for real bone and 42-61% for the dry beaded surface, with statistical significance (alpha<0.05). Less marked differences were observed for dynamic friction coefficients. Significant but non-systematic effects of normal stress or wet/dry condition on friction coefficients were observed within each surface type. These results indicate that the use of artificial bone models for pre-clinical implant testing that rely on interface load transfer with trabecular bone for mechanical integrity can be particularly sensitive to surface finish and lubrication conditions.     

Shear stability of an elastomeric disk spacer within an intervertebral joint: a parametric study.

A parametric study was conducted to compare the resistance to shear force provided by three surgical implantation techniques and three endplate designs for use with a polymeric lumbar intervertebral disk spacer. While under an axial load, the implant was pushed out laterally through the surgical window created for implantation. Force at 0.3 mm displacement, slope of the initial portion of the force-displacement graph, and maximum force were measured. The results indicate that implants with pegs, inclined planes, or domes on the surfaces of the disk will add significantly to the shear stability of the implant while maintaining the simplicity of a single part device.     

Primary stability of uncemented femoral resurfacing implants for varying interface parameters and material formulations during walking and stair climbing

Primary stability of uncemented resurfacing prosthesis is provided by an interference fit between the undersized implant and the reamed bone. Dependent on the magnitude of interference, the implantation process causes high shear forces and large strains which can exceed the elastic limit of cancellous bone. Plastification of the bone causes reduced stiffness and could lead to bone damage and implant loosening. The purpose in this study was to determine press-fit conditions which allow implantation without excessive plastic bone deformation and sufficient primary stability to achieve bone ingrowth. In particular, the influence of interference, bone quality and friction on the micromotion during walking and stair-climbing was investigated. Therefore elastic and plastic finite element (FE) models of the proximal femur were developed. Implantation was realized by displacing the prosthesis onto the femur while monitoring the contact pressure, plastic bone deformation as well as implantation forces. Subsequently a physiologic gait and stair-climbing cycle was simulated calculating the micromotion at the bone-implant interface. Results indicate that plastic deformation starts at an interference of 30 μm and the amount of plastified bone at the interface increases up to 90% at 150 μm interference. This effect did not reduce the contact pressure if interference was below 80 μm. The micromotion during walking was similar for the elastic and plastic FE models. A stable situation allowing bony ingrowth was achieved for both constitutive laws (elastic, plastic) for walking and stair climbing with at least 60 μm press-fit, which is feasible with clinically used implantation forces of 4 kN.     

An Improved Mechanical Testing Method to Assess Bone-implant Anchorage

Recent advances in material science have led to a substantial increase in the topographical complexity of implant surfaces, both on a micro- and a nano-scale. As such, traditional methods of describing implant surfaces - namely numerical determinants of surface roughness - are inadequate for predicting in vivo performance. Biomechanical testing provides an accurate and comparative platform to analyze the performance of biomaterial surfaces. An improved mechanical testing method to test the anchorage of bone to candidate implant surfaces is presented. The method is applicable to both early and later stages of healing and can be employed for any range of chemically or mechanically modified surfaces - but not smooth surfaces. Custom rectangular implants are placed bilaterally in the distal femora of male Wistar rats and collected with the surrounding bone. Test specimens are prepared and potted using a novel breakaway mold and the disruption test is conducted using a mechanical testing machine. This method allows for alignment of the disruption force exactly perpendicular, or parallel, to the plane of the implant surface, and provides an accurate and reproducible means for isolating an exact peri-implant region for testing.     

Sexual dimorphism in the attachment ability of the Colorado potato beetle Leptinotarsa decemlineata (Coleoptera: Chrysomelidae) to rough substrates

Many representatives of the beetle family Chrysomelidae exhibit a distinctive sexual dimorphism in the structure of adhesive tarsal setae. The present study demonstrates the influence of surface roughness on the friction force of Leptinotarsa decemlineata males and females. The maximum friction force of individual beetles was measured on epoxy resin surfaces (smooth and with asperities ranging from 0.3 to 12.0 microm) using a centrifugal force tester. On the smooth surface, no considerable differences between males and females were found, whereas on rough surfaces, females attached significantly (up to two times) stronger than males. Clawless beetles generated lower forces than intact ones, but demonstrated similar differences between males and females. The results indicate that the female adhesive system has its main functional trait in a stronger specialisation to rough plant surfaces whereas the adhesive system of males possess a certain trade-off between attachment to rough plant surfaces during locomotion on vegetation and to the smooth surface of the female elytra, while mating.     

Simulation of bone strain by orthodontic implants using the finite element method]

Load direction of applied forces, implant geometry and other biomechanical parameters lead to varying reactions in the surrounding bone structure. Three types of endosseous implant measuring 9 mm in length and 3.3 mm in diameter with and without superperiosteal step, and a threaded surface were investigated with the aid of a finite element method using the COSMOS/M 2.5 program. The load on the implant was investigated under vertical, horizontal, and diagonal forces of between 0.01 N and 100 N. Vertical loading of simple implants caused bone deformation of more than 600 mu eps. The application of the superperiosteal step clearly reduced the deformation. The largest deformations under vertical loading were observed in the trabecular bone with all 3 implant geometries. On horizontal loading the deformation shifted from the trabecular to the cortical bone and was particularly marked at the transition between the two. The smallest deformations, less than 300 mu eps, were measured at implants with a superperiosteal step under diagonal loading. The thread did not improve loading capacity. Implants with a superperiosteal step are recommended since they contribute to more rapid healing and strengthening of the bone.     

Influence of Different Abutment Designs on the Biomechanical Behavior of One-Piece Zirconia Dental Implants and Their Surrounding Bone: A 3D-FEA

BackgroundIn a dental implant/bone system, the design factors affect the value and distributions of stress and deformations that plays a pivotal role on the stability, durability and lifespan of the implant/bone system.ObjectiveThe aim of this study was to compare the influence of different abutment designs on the biomechanical behavior of one-piece zirconia dental implants and their surrounding bone tissues using three-dimensional finite element analysis.MethodsA three-dimensional geometrical model of a zirconia dental implant and its surrounding bone tissue were created. The occlusal loading force applied to the prosthetic abutments was a combination of 114.6 N in the axial direction, 17.1 N in the lingual direction and 23.4 N toward the mesial direction where these components represent masticatory force of 118.2 N in the angle of approximately 75° to the occlusal plane.ResultsThe system included implant abutment Model 01 showed a decrease of 9.58%, 9.92% and 3.62% at least in the average value of maximum von Mises stress compared to Model 02, Model 03 and Model 04 respectively. The results also showed that the system included implant abutment Model 01 decreases the average value of maximum deformation of 16.96%, 7.17% and 9.47% at least compared to Model 02, Model 03 and Model 04 respectively.ConclusionThe one-piece zirconia dental implant abutment Model 01 presents a better biomechanical behavior in the peri-implant bone than others. It can efficiently distribute the applied load and present more homogeneous behavior of stress distribution and has less deformation than others, which will enhance the stability of implant/bone system and prolong its lifespan.     

Effect of boundary conditions on yield properties of human femoral trabecular bone

Trabecular bone plays an important mechanical role in bone fractures and implant stability. Homogenized nonlinear finite element (FE) analysis of whole bones can deliver improved fracture risk and implant loosening assessment. Such simulations require the knowledge of mechanical properties such as an appropriate yield behavior and criterion for trabecular bone. Identification of a complete yield surface is extremely difficult experimentally but can be achieved in silico by using micro-FE analysis on cubical trabecular volume elements. Nevertheless, the influence of the boundary conditions (BCs), which are applied to such volume elements, on the obtained yield properties remains unknown. Therefore, this study compared homogenized yield properties along 17 load cases of 126 human femoral trabecular cubic specimens computed with classical kinematic uniform BCs (KUBCs) and a new set of mixed uniform BCs, namely periodicity-compatible mixed uniform BCs (PMUBCs). In stress space, PMUBCs lead to 7–72 % lower yield stresses compared to KUBCs. The yield surfaces obtained with both KUBCs and PMUBCs demonstrate a pressure-sensitive ellipsoidal shape. A volume fraction and fabric-based quadric yield function successfully fitted the yield surfaces of both BCs with a correlation coefficient \(R^{2} \ge 0.93\). As expected, yield strains show only a weak dependency on bone volume fraction and fabric. The role of the two BCs in homogenized FE analysis of whole bones will need to be investigated and validated with experimental results at the whole bone level in future studies.     

Surface porous fibre-reinforced composite bulk bone substitute

The aim of this study was to describe and evaluate the significance of a porous surface with bioactive glass granules (S53P4) covering an artificial bulk material based on polymethylmetacrylate (PMMA) and fibre-reinforced composite (FRC) technology. Effort was focused particularly on characters of the porous surface and biomechanical properties of the material in vitro, and test in vivo the implant in reconstruction in an experimental long bone segment defect model. The defect, 10 mm in length, created in the shaft of rabbit tibia, was reconstructed by the implant and fixed by intramedullary K-wires. The implant was incorporated within 4 weeks by new bone growth from the host bone covering particularly its posterior surface and cortex/implant junctions with bridging trabecular bone. Later, at 8 weeks, new bone was found also at the cortex/implant interface and in the medullary canal of the implant. Histometric measurements revealed direct bone/implant surface contact in 34% at the interface. Bioactive glass granules in the porous surface evoked the most direct contact with bone. The implants manufactured from PMMA only served as a control group, and showed significantly lower osteoconductive properties. Biomechanical measurements in vitro of fibre-reinforced PMMA specimens revealed values for bending strength and the flexural modulus to match them to human bone. This artificial bulk bone material based on PMMA/FRC technology seems to have proposing properties to be used as a bone substitute on load-bearing conditions. This revised version was published online in July 2006 with corrections to the Cover Date.     

Frictional heating of total hip implants. Part 1: measurements in patients

Hip implants heat up due to friction during long lasting, high loading activities like walking. Thermal damage in the surrounding soft and hard tissues and deteriorated lubrication of synovial fluid could contribute to implant loosening. The goal of this study was to determine the implant temperatures in vivo under varying conditions. Temperatures and contact forces in the joints were measured in seven joints of five patients using instrumented prostheses with alumina ceramic heads and telemetry data transmission. The peak temperature in implants with polyethylene cups rose up to 43.1 degrees C after an hour of walking but varied considerably individually. Even higher temperatures at the joints are probable for patients with higher body weight or while jogging. The peak temperature was lower with a ceramic cup, showing the influence of friction in the joint. During cycling the peak temperatures were lower than during walking, proving the effect of force magnitudes on the produced heat. However, no positive correlation was found between force magnitude and maximum temperature during walking. Other individual parameters than just the joint force influence the implant temperatures. Based on the obtained data and the available literature about thermal damage of biological tissues a detrimental effect of friction induced heat on the stability of hip implants cannot be excluded. Because the potential risk for an individual patient cannot be foreseen, the use and improvement of low friction implant materials is important.     

Finite element analysis of dental implant loading on atrophic and non-atrophic cancellous and cortical mandibular bone – a feasibility study

The first aim of this study was to assess displacements and micro-strain induced on different grades of atrophic cortical and trabecular mandibular bone by axially loaded dental implants using finite element analysis (FEA). The second aim was to assess the micro-strain induced by different implant geometries and the levels of bone-to-implant contact (BIC) on the surrounding bone. Six mandibular bone segments demonstrating different grades of mandibular bone atrophy and various bone volume fractions (from 0.149 to 0.471) were imaged using a micro-CT device. The acquired bone STL models and implant (Brånemark, Straumann, Ankylos) were merged into a three-dimensional finite elements structure. The mean displacement value for all implants was 3.1±1.2 µm. Displacements were lower in the group with a strong BIC. The results indicated that the maximum strain values of cortical and cancellous bone increased with lower bone density. Strain distribution is the first and foremost dependent on the shape of bone and architecture of cancellous bone. The geometry of the implant, thread patterns, grade of bone atrophy and BIC all affect the displacement and micro-strain on the mandible bone. Preoperative finite element analysis could offer improved predictability in the long-term outlook of dental implant restorations.     

Experimental and Numerical Analyses of the Pull-out Response of a Steel Post/Bovine Bone Cementless Fixation

Effect of initial interference fit on pull-out strength in cementless fixation between bovine tibia and smooth stainless steel post was investigated in this study.Compressive behavior of bovine spongious bone was studied using mechanical testing in order to evaluate the elastic-plastic properties in different regions of the proximal tibia.Friction tests were carried out in the aim to evaluate the friction behavior of the contact between bovine spongious bone and stainless steel.A cylindrical stainless steel post inserted in a pre-drilled bovine tibia with an initial interference fit was taken as an in vitro model to assess the contribution of post fixation to the initial stability of the Total Knee Arthroplasty (TKA) tibial component.Pull-out experiments were carried out for different initial interference fits.Finite Element Models (FEM) using local elastic-plastic properties of the bovine bone were developed for the analysis of the experimental ultimate pull-out force results.At the post/bone interface,Coulomb friction was considered in the FEM calculations with pressure-dependent friction coefficient.It was found that the FEM results of the ultimate force are in good agreement with the experimental results.The analysis of the FEM interfacial stresses indicates that the micro-slip initiation depends on the local bone properties.     

Adaptive diffuse domain approach for calculating mechanically induced deformation of trabecular bone

Remodelling of trabecular bone is essentially affected by the mechanical load of the trabeculae. Mathematical modelling and simulation of the remodelling process have to include time-consuming calculations of the displacement field within the complex trabecular structure under loading. We present an adaptive diffuse domain approach for calculating the elastic bone deformation based on micro computer tomogram data of real trabecular bone structures and compared it with a conventional voxel-based finite element method. In addition to allowing for higher computational efficiency, the adaptive approach is characterised by a very smooth representation of the bone surface, which suggests that this approach would be suitable as a basis for future simulations of bone resorption and formation processes within the trabecular structure.     

Factors contributing to the integration of textural qualities: evidence from virtual surfaces

Two experiments involving indirect touch were carried out to explore the relationships among perceptual dimensions of haptically examined surfaces. Subjects in both experiments used a stylus to evaluate the properties of virtual surfaces created by a force-feedback device; four surface properties (resistance to normal force, coefficient of friction, texture scale, and vibration amplitude) were manipulated in various combinations. In Experiment 1, the extent to which there was a one-to-one relationship between specific stimulus properties and perceptual qualities ("perceptual separability") was evaluated. A substantial failure of separability was demonstrated, with friction tending to be more separable from the other properties than they were from one another. The pattern of results suggests that the amount of measured separability depends crucially on the way stimulus properties are defined (e.g., force versus displacement). In Experiment 2, surfaces with known perceptual properties were used to study the metric(s) of the relevant perceptual space. By specifying the perceptual, rather than the stimulus, coordinates of the surfaces, it was possible to bypass issues of perceptual separability. For surfaces of equal friction, a Euclidean metric captured the results (r(2) = 0.75) more effectively than a city-block metric did; neither metric did well when differences in friction were involved. The fact that-unlike stickiness-hardness, roughness, and perceived vibration intensity are all increasing functions of surface-normal forces may facilitate their integration into a Euclidean space, in both direct and indirect touch.     

Trabecular bone strains around a dental implant and associated micromotions—A micro-CT-based three-dimensional finite element study

The first objective of this computational study was to assess the strain magnitude and distribution within the three-dimensional (3D) trabecular bone structure around an osseointegrated dental implant loaded axially. The second objective was to investigate the relative micromotions between the implant and the surrounding bone. The work hypothesis adopted was that these virtual measurements would be a useful indicator of bone adaptation (resorption, homeostasis, formation).In order to reach these objectives, a μCT-based finite element model of an oral implant implanted into a Berkshire pig mandible was developed along with a robust software methodology. The finite element mesh of the 3D trabecular bone architecture was generated from the segmentation of μCT scans. The implant was meshed independently from its CAD file obtained from the manufacturer. The meshes of the implant and the bone sample were registered together in an integrated software environment. A series of non-linear contact finite element (FE) analyses considering an axial load applied to the top of the implant in combination with three sets of mechanical properties for the trabecular bone tissue was devised. Complex strain distribution patterns are reported and discussed. It was found that considering the Young’s modulus of the trabecular bone tissue to be 5, 10 and 15 GPa resulted in maximum peri-implant bone microstrains of about 3000, 2100 and 1400. These results indicate that, for the three sets of mechanical properties considered, the magnitude of maximum strain lies within an homeostatic range known to be sufficient to maintain/form bone. The corresponding micro-motions of the implant with respect to the bone microstructure were shown to be sufficiently low to prevent fibrous tissue formation and to favour long-term osseointegration.     

Rib fractures under anterior–posterior dynamic loads: Experimental and finite-element study

The purpose of this study was to investigate whether using a finite-element (FE) mesh composed entirely of hexahedral elements to model cortical and trabecular bone (all-hex model) would provide more accurate simulations than those with variable thickness shell elements for cortical bone and hexahedral elements for trabecular bone (hex–shell model) in the modeling human ribs. First, quasi-static non-injurious and dynamic injurious experiments were performed using the second, fourth, and tenth human thoracic ribs to record the structural behavior and fracture tolerance of individual ribs under anterior–posterior bending loads. Then, all-hex and hex–shell FE models for the three ribs were developed using an octree-based and multi-block hex meshing approach, respectively. Material properties of cortical bone were optimized using dynamic experimental data and the hex–shell model of the fourth rib and trabecular bone properties were taken from the literature. Overall, the reaction force–displacement relationship predicted by both all-hex and hex–shell models with nodes in the offset middle-cortical surfaces compared well with those measured experimentally for all the three ribs. With the exception of fracture locations, the predictions from all-hex and offset hex–shell models of the second and fourth ribs agreed better with experimental data than those from the tenth rib models in terms of reaction force at fracture (difference <15.4%), ultimate failure displacement and time (difference <7.3%), and cortical bone strains. The hex–shell models with shell nodes in outer cortical surfaces increased static reaction forces up to 16.6%, compared to offset hex–shell models. These results indicated that both all-hex and hex–shell modeling strategies were applicable for simulating rib responses and bone fractures for the loading conditions considered, but coarse hex–shell models with constant or variable shell thickness were more computationally efficient and therefore preferred.