Computational hip joint simulator for wear and heat generation

This paper presents a computational simulator for the hip to compute the wear and heat generation on artificial joints. The friction produced on artificial hip joints originates wear rates that can lead to failure of the implant. Furthermore, the frictional heating can increase the wear. The developed computational model calculates the wear in the joint and the temperature in the surrounding zone, allowing the use of different combinations of joint materials, daily activities and different individuals. The pressure distribution on the joint bearing surfaces is obtained with the solution of a contact model. The heat generation by friction and the volumetric wear is computed from the pressure distribution and the sliding distance. The temperature is obtained from the solution of a transient heat conduction problem that includes the time-dependent heat generated by friction. The contact and heat conduction problems are solved numerically with the Finite Element Method. The developed computational model performs a full simulation of the acetabular bearing surface behaviour, which is useful for acetabular cup design and material selection. The results obtained by the present model agree with experimental and clinical data, as well as other numerical studies.     

Research on the Interstitial Fluid Load Support Characteristics and Start-Up Friction Mechanisms of PVA-HA-Silk Composite Hydrogel

Hydrogel has been extensively studied as an articular cartilage repair and replacement material. PVA-HA-Silk composite hydrogel was prepared by freezing-thawing method in this paper. Mechanical properties were determined by experiments and the friction coefficient of PVA-HA-Silk composite hydrogel against steel ball was verified using micro-tribometer. Finite Element Method （FEM） was used to study the lubrication mechanism of PVA-HA-Silk composite hydrogel and the relation between the interstitial fluid load support and the start-up friction resistance. The results show that the elastic modulus and the permeability are 2.07 MPa and 10^-15m^4N^-1s^-1, respectively, and the start-up friction coefficients of PVA-HA-Silk composite hydrogel are in the range of 0.154）.2 at different contact loads, contact time and sliding speeds. The start-up friction resistance of PVA-HA-Silk composite hydrogel increases with the contact load and contact time. With the increase in sliding speed, the start-up friction resistance of PVA-HA-Silk composite hydrogel decreases. There is an inverse relation between the start-up friction resistance and the interstitial fluid load support. The change of fluid flow with the increase in sliding displacement has an important effect on the interstitial fluid load support and friction resistance. The interstitial fluid load support decreases with the increase in contact load and contact time, while the interstitial fluid load support reinforces with the increase in sliding speed. Moreover, PVA-HA-Silk composite hydrogel has mechanical properties of recovery and self-lubricating.     

The Influence of Normal Load and Sliding Speed on Frictional Properties of Skin

The study of frictional properties of human skin is important for medical research, skin care products and textile exploi- tation. In order to investigate the influence of normal load and sliding speed on the frictional properties of skin and its possible mechanism, tests were carded out on a multi-specimen friction tester. When the normal load increases from 0.1 N to 0.9 N, normal displacement and the friction coefficient of skin increase. The friction coefficient is dependent on the load, indicating that both adhesion and deformation contribute to the friction behaviour. The deformation friction was interpreted using the plough model of friction. When sliding speed increases from 0.5 mm·s^-1 to 4 mm·s^-1, the friction coefficient increases and ＂stick-slip＂ phenomena increase, indicating that hysteretic friction contributes to the friction. The hysteretic friction was in- terpreted using schematic of energy translation during the rigid spherical probe sliding on the soft skin surface, which provides an explanation for the influence of the sliding speed on the frictional characteristics of the skin.     

Minimal adhesion surface area in tangentially loaded digital contacts

The stick-to-slip transition of a fingertip in contact with a planar surface does not occur instantaneously. As the tangential load increases, portions of the skin adhere while others slip, giving rise to an evolution of the contact state, termed partial slip. We develop a quasi-static model that predicts that if the coefficient of kinetic friction is larger than the coefficient of static friction, then the stuck surface area diminishes as the tangential load increases until reaching a 'minimal adhesion surface area' where it vanishes abruptly. This phenomenon was observed in recently measured finger-slip image data (André et al., 2011) that were processed by an optic flow detection algorithm. We examined the results of 10 trials. Four of them exhibited the minimal adhesion surface area phenomenon, four of them did not, and two were inconclusive.     

Physical Factors Influencing Pleasant Touch during Tactile Exploration

Background

When scanning surfaces, humans perceive some of their physical attributes. These percepts are frequently accompanied by a sensation of (un)pleasantness. We therefore hypothesized that aspects of the mechanical activity induced by scanning surfaces with fingertips could be objectively associated with a pleasantness sensation. Previously, we developed a unidimensional measure of pleasantness, the Pleasant Touch Scale, quantifying the pleasantness level of 37 different materials. Findings of this study suggested that the sensation of pleasantness was influenced by the average magnitude of the frictional forces brought about by sliding the finger on the surface, and by the surface topography. In the present study, we correlated (i) characteristics of the fluctuations of frictional forces resulting from the interaction between the finger and the surface asperities as well as (ii) the average friction with the sensation of pleasantness.

Results

Eight blindfolded participants tactually explored twelve materials of the Pleasant Touch Scale through lateral sliding movements of their index fingertip. During exploration, the normal and tangential interaction force components, f_N and f_T, as well as the fingertip trajectory were measured. The effect of the frictional force on pleasantness sensation was investigated through the analysis of the ratio f_T to f_N, i.e. the net coefficient of kinetic friction, μ. The influence of the surface topographies was investigated through analysis of rapid f_T fluctuations in the spatial frequency domain. Results showed that high values of μ were anticorrelated with pleasantness. Furthermore, surfaces associated with fluctuations of f_T having higher amplitudes in the low frequency range than in the high one were judged to be less pleasant than the surfaces yielding evenly distributed amplitudes throughout the whole spatial frequency domain.

Conclusion

Characteristics of the frictional force fluctuations and of the net friction taking place during scanning can reliably be correlated with the pleasantness sensation of surfaces.     

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.     

A contact model to simulate human–artifact interaction based on force optimization: implementation and application to the analysis of a training machine

Musculoskeletal multibody models are increasingly used to analyze and optimize physical interactions between humans and technical artifacts. Since interaction is conveyed by contact between the human body and the artifact, a computationally robust modeling approach for frictional contact forces is a crucial aspect. In this contribution, we propose a parametric contact model and formulate an associated force optimization problem to simultaneously estimate unknown muscle and contact forces in an inverse dynamic manner from a prescribed motion trajectory. Unlike existing work, we consider both the static and the kinetic regime of Coulomb’s friction law. The approach is applied to the analysis of a leg extension training machine with the objective to reduce the stress on the tibiofemoral joint. The uncertainty of the simulation results due to a tunable parameter of the contact model is of particular interest.     

Determinants of friction in soft elastohydrodynamic lubrication

Elastohydrodynamic lubrication (EHL) protects soft tissues from damage and wear in many biological systems (e.g. synovial joints, cornea of the eye, and pleural surfaces of the lung and chest wall). Among studies of lubrication of deformable solids, few have examined the effects of external loads, geometry, and material properties on EHL of soft tissues. To examine these effects, we studied the tribology of soft tissues in a two-dimensional finite element simulation of a thin layer of fluid separating a sliding rigid surface from a soft asperity or bump with an initial sinusoidal shape. We computed the frictional force, deformation of the solid, and change in fluid thickness as functions of independent variables: sliding velocity, normal load, material properties, and bump amplitude and length. Double-logarithmic regression was used to determine the exponents of the scaling relationships of friction coefficient and minimum fluid thickness to the independent variables. The analysis showed that frictional shear force is strongly dependent on velocity, viscosity, and load, moderately dependent on bump length and elasticity, and only weakly dependent on the bump amplitude. The minimum fluid thickness is strongly dependent on velocity and viscosity, and changes moderately with load, elasticity, amplitude, and length. The shape of the bump has little effect. The results confirm that the shear-induced deformation of an initially symmetrical shape, including generalizations to other symmetrical geometries such as quadratic or piecewise linear bumps, leads to load-supporting behavior.     

Method and apparatus for soft tissue material parameter estimation using tissue tagged Magnetic Resonance Imaging

We describe an experimental method and apparatus for the estimation of constitutive parameters of soft tissue using Magnetic Resonance Imaging (MRI), in particular for the estimation of passive myocardial material properties. MRI tissue tagged images were acquired with simultaneous pressure recordings, while the tissue was cyclically deformed using a custom built reciprocating pump actuator A continuous three-dimensional (3D) displacement field was reconstructed from the imaged tag motion. Cavity volume changes and local tissue microstructure were determined from phase contrast velocity and diffusion tensor MR images, respectively. The Finite Element Method (FEM) was used to solve the finite elasticity problem and obtain the displacement field that satisfied the applied boundary conditions and a given set of material parameters. The material parameters which best fit the FEM predicted displacements to the displacements reconstructed from the tagged images were found by nonlinear optimization. The equipment and method were validated using inflation of a deformable silicon gel phantom in the shape of a cylindrical annulus. The silicon gel was well described by a neo-Hookian material law with a single material parameter C1=8.71+/-0.06kPa, estimated independently using a rotational shear apparatus. The MRI derived parameter was allowed to vary regionally and was estimated as C1 =8.80+/-0.86kPa across the model. Preliminary results from the passive inflation of an isolated arrested pig heart are also presented, demonstrating the feasibility of the apparatus and method for isolated heart preparations. FEM based models can therefore estimate constitutive parameters accurately and reliably from MRI tagging data.     

Molecular friction in an actomyosin molecular machine

In muscle contraction, it has been widely recognized that a binding state exists between myosin and actin in the presence of Mg-ATP. To estimate the magnitude of binding strength, I introduce a concept of frictional phenomena which occurs between two sliding bodies in contact each other. In such cases, the sliding speed can be formulated as a function of the actin-myosin bond strength. In order to validate this, the present theory is applied for the two movement assay systems with no external load; one movement assay of Phalloidin Rhodamine bound F-actin on a myosin coated hydrophobic cover glass and another assay of myosin coated beads along actin cables of Nitella. If a coefficient of 0.005 is applied to the kinetic friction, 1pN for the sliding force per cross-bridge and 10 microns sec-1 for the sliding speed, it is found that the bond strength between actin and one myosin head is about 200 pN in the contracting state.     

The Frictional Coefficient of Bovine Knee Articular Cartilage

1 Introduction Cartilage has excellent biomechanical and tri- bological properties with low friction and minimum wear in diarthrodial joints throughout the lifetime of most people, and the lifetime of articular cartilage can be 40 years or longer. This has inspired material and bionic scientists to study the mechanism of such excellent tri- bological characteristics in order to develop artificial joints. Various mechanisms have been proposed to ex- plain the remarkable low friction behavior of…     

Three-dimensional finite element simulations of the mechanical response of the fingertip to static and dynamic compressions

The analysis of the mechanics of the contact interactions of fingers/handle and the stress/strain distributions in the soft tissues in the fingertip is essential to optimize design of tools to reduce many occupation-related hand disorders. In the present study, a three-dimensional (3D) finite element (FE) model for the fingertip is proposed to simulate the nonlinear and time-dependent responses of a fingertip to static and dynamic loadings. The proposed FE model incorporates the essential anatomical structures of a finger: skin layers (outer and inner skins), subcutaneous tissue, bone and nail. The soft tissues (inner skin and subcutaneous tissue) are considered to be nonlinearly viscoelastic, while the hard tissues (outer skin, bone and nail) are considered to be linearly elastic. The proposed model has been used to simulate two loading scenarios: (a) the contact interactions between the fingertip and a flat surface and (b) the indentation of the fingerpad via a sharp wedge. For case (a), the predicted force/displacement relationships and time-dependent force responses are compared with the published experimental data; for case (b), the skin surface deflection profiles were predicted and compared with the published experimental observations. Furthermore, for both cases, the time-dependent stress/strain distributions within the tissues of the fingertip were calculated. The good agreement between the model predictions and the experimental observations indicates that the present model is capable of predicting realistic time-dependent force/displacement responses and stress/strain distributions in the soft tissues for dynamic loading conditions.     

Frictional properties of regenerated cartilage in vitro

Although tribological function is the most important mechanical property of articular cartilage, few studies have examined this function in tissue-engineered cartilage. We investigated changes in the frictional properties of cartilage regenerated from the inoculation of rabbit chondrocytes into fibroin sponge. A reciprocating friction-testing apparatus was used to measure the friction coefficient of the regenerated cartilage under a small load. The specimen was slid against a stainless steel plate in a water vessel filled with physiological saline. The applied load was 0.03 N, the stroke length was 20 mm, and the mean sliding velocity was 0.8 mm/s. The friction coefficient of the regenerated cartilage decreased with increasing cultivation time, because a hydrophilic layer of synthesized extracellular matrix was formed on the fibroin sponge surface. The friction coefficient of the regenerated cartilage was as low as that of natural cartilage in the early stages of the sliding tests, but it increased with increasing duration of sliding owing to exudation of interstitial water from the surface layer.     

Three-dimensional finite element simulations of the mechanical response of the fingertip to static and dynamic compressions

The analysis of the mechanics of the contact interactions of fingers/handle and the stress/strain distributions in the soft tissues in the fingertip is essential to optimize design of tools to reduce many occupation-related hand disorders. In the present study, a three-dimensional (3D) finite element (FE) model for the fingertip is proposed to simulate the nonlinear and time-dependent responses of a fingertip to static and dynamic loadings. The proposed FE model incorporates the essential anatomical structures of a finger: skin layers (outer and inner skins), subcutaneous tissue, bone and nail. The soft tissues (inner skin and subcutaneous tissue) are considered to be nonlinearly viscoelastic, while the hard tissues (outer skin, bone and nail) are considered to be linearly elastic. The proposed model has been used to simulate two loading scenarios: (a) the contact interactions between the fingertip and a flat surface and (b) the indentation of the fingerpad via a sharp wedge. For case (a), the predicted force/displacement relationships and time-dependent force responses are compared with the published experimental data; for case (b), the skin surface deflection profiles were predicted and compared with the published experimental observations. Furthermore, for both cases, the time-dependent stress/strain distributions within the tissues of the fingertip were calculated. The good agreement between the model predictions and the experimental observations indicates that the present model is capable of predicting realistic time-dependent force/displacement responses and stress/strain distributions in the soft tissues for dynamic loading conditions.     

Modeling of time-dependent force response of fingertip to dynamic loading

An extended exposure to repeated loading on fingertip has been associated to many vascular, sensorineural, and musculoskeletal disorders in the fingers, such as carpal tunnel syndrome, hand-arm vibration syndrome, and flexor tenosynovitis. A better understanding of the pathomechanics of these sensorineural and vascular diseases in fingers requires a formulation of a biomechanical model of the fingertips and analyses to predict the mechanical responses of the soft tissues to dynamic loading. In the present study, a model based on finite element techniques has been developed to simulate the mechanical responses of the fingertips to dynamic loading. The proposed model is two-dimensional and incorporates the essential anatomical structures of a finger: skin, subcutaneous tissue, bone, and nail. The skin tissue is assumed to be hyperelastic and viscoelastic. The subcutaneous tissue was considered to be a nonlinear, biphasic material composed of a hyperelastic solid and an invicid fluid, while its hydraulic permeability was considered to be deformation dependent. Two series of numerical tests were performed using the proposed finger tip model to: (a) simulate the responses of the fingertip to repeated loading, where the contact plate was assumed to be fixed, and the bone within the fingertip was subjected to a prescribed sinusoidal displacement in vertical direction; (b) simulate the force response of the fingertip in a single keystroke, where the keyboard was composed of a hard plastic keycap, a rigid support block, and a nonlinear spring. The time-dependent behavior of the fingertip under dynamic loading was derived. The model predictions of the time-histories of force response of the fingertip and the phenomenon of fingertip separation from the contacting plate during cyclic loading agree well with the reported experimental observations.     

Analysis of the Biotribological Behavior for the Stainless Steel/Polyethylene Contact Using a Knee Prosthesis Simulator

In this work,a friction and wear simulator was used to reproduce the Anterior-Posterior (AP) sliding and the Flexion-Extension (FE) rotation generated in the knee joint during human gait cycle.We chose to simplify the contact geometry between the Total Knee Arthroplasty (TKA) femoral component and tibial insert.A 304L stainless steel cylinder which replaces the femoral component was loaded onto a flat High Density Polyethylene (HDPE) block which replaces the tibial insert.The tribological behavior of the considered contact was analyzed by tracking the number of cycles,the friction coefficient,the roughness of the wear track on HDPE,the HDPE weight loss and the damage mechanisms.The friction coefficient shows a gradual increase with the number of cycles for both AP and FE kinematics.The evolution of friction coefficient with the number of cycles is not affected by the value of the imposed normal load in the case of AP sliding.For the FE rotation,decreased friction coefficient is obtained when the imposed normal load increases.For both considered AP and FE kinematics,the roughness of the wear track on the HDPE is not affected by the imposed normal load.It shows a progressive decrease when the number of cycles increases.The wear of HDPE obeys the Archard law and the wear coefficient increases with the normal force.For a given value of normal load,the obtained wear coefficient for the AP sliding is larger than that obtained for FE rotation.A predominant adhesive wear mechanism was identified for both AP and FE kinematics.Under the same normal load,damage development in terms of plastic deformation,micro-cracking and debonding is more pronounced for the AP sliding if compared with the FE rotation.For a given kinematics,the damage severity increases with the normal load.This finding is in good agreement with the predicted values of the wear coefficient according to the Archard law.     

Biomechanical analysis of the three-dimensional foot structure during gait: a basic tool for clinical applications

A novel three-dimensional numerical model of the foot, incorporating, for the first time in the literature, realistic geometric and material properties of both skeletal and soft tissue components of the foot, was developed for biomechanical analysis of its structural behavior during gait. A system of experimental methods, integrating the optical Contact Pressure Display (CPD) method for plantar pressure measurements and a Digital Radiographic Fluoroscopy (DRF) instrument for acquisition of skeletal motion during gait, was also developed in this study and subsequently used to build the foot model and validate its predictions. Using a Finite Element solver, the stress distribution within the foot structure was obtained and regions of elevated stresses for six subphases of the stance (initial-contact, heel-strike, midstance, forefoot-contact, push-off, and toe-off) were located. For each of these subphases, the model was adapted according to the corresponding fluoroscopic data, skeletal dynamics, and active muscle force loading. Validation of the stress state was achieved by comparing model predictions of contact stress distribution with respective CPD measurements. The presently developed measurement and numerical analysis tools open new approaches for clinical applications, from simulation of the development mechanisms of common foot disorders to pre- and post-interventional evaluation of their treatment.     

Influence of the stiffness of bone defect implants on the mechanical conditions at the interface--a finite element analysis with contact

The study focused on the influence of the implant material stiffness on stress distribution and micromotion at the interface of bone defect implants. We hypothesized that a low-stiffness implant with a modulus closer to that of the surrounding trabecular bone would yield a more homogeneous stress distribution and less micromotion at the interface with the bony bed. To prove this hypothesis we generated a three-dimensional, non-linear, anisotropic finite element (FE) model. The FE model corresponded to a previously developed animal model in sheep. A prismatic implant filled a standardized defect in the load-bearing area of the trabecular bone beneath the tibial plateau. The interface was described by face-to-face contact elements, which allow press fits, friction, sliding, and gapping. We assumed a physiological load condition and calculated contact pressures, shear stresses, and shear movements at the interface for two implants of different stiffness (titanium: E=110GPa; composite: E=2.2GPa). The FE model showed that the stress distribution was more homogeneous for the low-stiffness implant. The maximum pressure for the composite implant (2.1 MPa) was lower than for the titanium implant (5.6 MPa). Contrary to our hypothesis, we found more micromotion for the composite (up to 6 microm) than for the titanium implant (up to 4.5 microm). However, for both implants peak stresses and micromotion were in a range that predicts adequate conditions for the osseointegration. This was confirmed by the histological results from the animal studies.     

Design and numerical implementation of a 3-D non-linear viscoelastic constitutive model for brain tissue during impact

Finite Element (FE) head models are often used to understand mechanical response of the head and its contents during impact loading in the head. Current FE models do not account for non-linear viscoelastic material behavior of brain tissue. We developed a new non-linear viscoelastic material model for brain tissue and implemented it in an explicit FE code. To obtain sufficient numerical accuracy for modeling the nearly incompressible brain tissue, deviatoric and volumetric stress contributions are separated. Deviatoric stress is modeled in a non-linear viscoelastic differential form. Volumetric behavior is assumed linearly elastic. Linear viscoelastic material parameters were derived from published data on oscillatory experiments, and from ultrasonic experiments. Additionally, non-linear parameters were derived from stress relaxation (SR) experiments at shear strains up to 20%. The model was tested by simulating the transient phase in the SR experiments not used in parameter determination (strains up to 20%, strain rates up to 8s(-1)). Both time- and strain-dependent behavior were predicted accurately (R2>0.96) for strain and strain rates applied. However, the stress was overestimated systematically by approximately 31% independent of strain(rate) applied. This is probably caused by limitations of the experimental data at hand.     

Biomechanics of smooth adhesive pads in insects: influence of tarsal secretion on attachment performance

Many insects possess smooth adhesive pads on their legs, which adhere by thin films of a two-phasic secretion. To understand the function of such fluid-based adhesive systems, we simultaneously measured adhesion, friction and contact area in single pads of stick insects (Carausius morosus). Shear stress was largely independent of normal force and increased with velocity, seemingly consistent with the viscosity-effect of a continuous fluid film. However, measurements of the remaining force 2 min after a sliding movement show that adhesive pads can sustain considerable static friction. Repeated sliding movements and multiple consecutive pull-offs to deplete adhesive secretion showed that on a smooth surface, friction and adhesion strongly increased with decreasing amount of fluid. In contrast, pull-off forces significantly decreased on a rough substrate. Thus, the secretion does not generally increase attachment but does so only on rough substrates, where it helps to maximize contact area. When slides were repeated at one position so that secretion could accumulate, sliding shear stress decreased but static friction remained clearly present. This suggests that static friction which is biologically important to prevent sliding is based on non-Newtonian properties of the adhesive emulsion rather than on a direct contact between the cuticle and the substrate.