Porohyperelastic finite element model for the kangaroo humeral head cartilage based on experimental study and the consolidation theory

Solid-extracellular fluid interaction is believed to play an important role in the strain-rate dependent mechanical behaviors of shoulder articular cartilages. It is believed that the kangaroo shoulder joint is anatomically and biomechanically similar to human shoulder joint and it is easy to get in Australia. Therefore, the kangaroo humeral head cartilage was used as the suitable tissue for the study in this paper. Indentation tests from quasi-static (10-4/sec) to moderately high strain-rate (10-2/sec) on kangaroo humeral head cartilage tissues were conduced to investigate the strain-rate dependent behaviors. A finite element (FE) model was then developed, in which cartilage was conceptualized as a porous solid matrix filled with incompressible fluids. In this model, the solid matrix was modeled as an isotropic hyperelastic material and the percolating fluid follows Darcy’s law. Using inverse FE procedure, the constitutive parameters related to stiffness, compressibility of the solid matrix and permeability were obtained from the experimental results. The effect of solid-extracellular fluid interaction and drag force (the resistance to fluid movement) on strain-rate dependent behavior was investigated by comparing the influence of constant, strain dependent and strain-rate dependent permeability on FE model prediction. The newly developed porohyperelastic cartilage model with the inclusion of strain-rate dependent permeability was found to be able to predict the strain-rate dependent behaviors of cartilages.

Can measures of limb loading and dynamic stability during the squat maneuver provide an index of early functional recovery following unilateral total hip arthroplasty?

Objective: To investigate limb loading and dynamic stability during squatting in the early functional recovery of total hip arthroplasty (THA) patients. Design: Cohort study Setting: Inpatient rehabilitation clinic. Participants: A random sample of 61 THA patients (34♂/27♀; 62±9 yrs, 77±14 kg, 174±9 cm) was assessed twice, 13.2±3.8 days (PRE) and 26.6±3.3 days post-surgery (POST), and compared with a healthy reference group (REF) (22♂/16♀; 47±12yrs; 78±20kg; 175±10cm). Interventions: THA patients received two weeks of standard in-patient rehabilitation. Main Outcome Measure(s): Inter-limb vertical force distribution and dynamic stability during the squat maneuver, as defined by the root mean square (RMS) of the center of pressure in antero-posterior and medio-lateral directions, of operated (OP) and non-operated (NON)limbs. Self-reported function was assessed via FFb-H-OA 2.0 questionnaire. Results: At PRE, unloading of the OP limb was 15.8% greater (P<.001, d=1.070) and antero-posterior and medio-lateral center of pressure RMS were 30-34% higher in THA than REF P<.05). Unloading was reduced by 12.8% towards a more equal distribution from PRE to POST (P<.001, d=0.874). Although medio-lateral stability improved between PRE and POST (OP: 14.8%, P=.024, d=0.397; NON: 13.1%, P=.015, d=0.321), antero-posterior stability was not significantly different. Self-reported physical function improved by 15.8% (P<.001, d=0.965). Conclusion(s): THA patients unload the OP limb and are dynamically more unstable during squatting in the early rehabilitation phase following total hip replacement than healthy adults. Although loading symmetry and medio-lateral stability improved to the level of healthy adults with rehabilitation, antero-posterior stability remained impaired. Measures of dynamic stability and load symmetry during squatting provide quantitative information that can be used to clinically monitor early functional recovery from THA.

CORR Insights® : loss of cement-bone interlock in retrieved tibial components from total knee arthroplasties

Most surgeons cement the tibial component in total knee replacement surgery. Mid-term registry data from a number of countries, including those from the United Kingdom and Australia, support the excellent survivorship of cemented tibial components. In spite of this success, results can always be improved, and cementing technique can play a role. Cementing technique on the tibia is not standardized, and surgeons still differ about the best ways to deliver cement into the cancellous bone of the upper tibia. Questions remain regarding whether to use a gun or a syringe to inject the cement into the cancellous bone of the tibial plateau . The ideal cement penetration into the tibial plateau is debated, though most reports suggest that 4 mm to 10 mm is ideal. Thicker mantles are thought to be dangerous due to the risk of bone necrosis, but there is little in the literature to support this contention...

The effect of knee flexion contracture following total knee arthroplasty on the energy cost of walking.

This study evaluated the energy cost of walking (Cw) with knee flexion contractures (FC) simulated with a knee brace, in total knee arthroplasty (TKA) recipients (n=16) and normal controls (n=15), and compared it to baseline (no brace). There was no significant difference in Cw between the groups at baseline but TKA recipients walked slower (P=0.048) and with greater knee flexion in this condition (P=0.003). Simulated FC significantly increased Cw in both groups (TKA P=0.020, control P=0.002) and this occurred when FC exceeded 20° in the TKA group and 15° in the controls. Reported perceived exertion was only significantly increased by FC in the control group (control P<0.001, TKA P=0.058). Simulated knee FCs less than 20° do not increase Cw or perceived exertion in TKA recipients.

Adaptive finite element analysis of elliptic problems based on bubble-type local mesh generation

A new mesh adaptivity algorithm that combines a posteriori error estimation with bubble-type local mesh generation (BLMG) strategy for elliptic differential equations is proposed. The size function used in the BLMG is defined on each vertex during the adaptive process based on the obtained error estimator. In order to avoid the excessive coarsening and refining in each iterative step, two factor thresholds are introduced in the size function. The advantages of the BLMG-based adaptive finite element method, compared with other known methods, are given as follows: the refining and coarsening are obtained fluently in the same framework; the local a posteriori error estimation is easy to implement through the adjacency list of the BLMG method; at all levels of refinement, the updated triangles remain very well shaped, even if the mesh size at any particular refinement level varies by several orders of magnitude. Several numerical examples with singularities for the elliptic problems, where the explicit error estimators are used, verify the efficiency of the algorithm. The analysis for the parameters introduced in the size function shows that the algorithm has good flexibility.

Finite element analysis of CFRP strengthened steel hollow sections under tension

This paper presents a nonlinear finite element (FE) model for the analysis of very high strength (VHS) steel hollow sections wrapped by high modulus carbon fibre rein forced polymer (CFRP) sheets. The bond strength of CFRP wrapped VHS circular steel hollow section under tension is investigated using the FE model. The three dimensional FE model by Nonlinear static analysis has been carried out by Strand 7 finite element software. The model is validated by the experimental data obtained from Fawzia et al [1]. A detail parametric study has been performed to examine the effect of number of CFRP layers, different diameters of VHS steel tube and different bond lengths of CFRP sheet. The analytical model developed by Fawzia et al. [1] has been used to determine the load carrying capacity of different diameters of CFRP strengthened VHS steel tube by using the capacity from each layer of CFRP sheet. The results from FE model have found in reasonable agreement with the analytical model developed by Fawzia et al [1]. This validation was necessary because the analytical model by Fawzia et al [1] was developed by using only one diameter of VHS steel tube and fixed (five) number of CFRP layers. It can be concluded that the developed analytical model is valid for CFRP strengthened VHS steel tubes with diameter range of 38mm to 100mm and CFRP layer range of 3 to 5 layers. Based on the results it can also be concluded that the effective bond length is consistent for different diameters of steel tubes and different layers of CFRP. Three layers of CFRP is considered most effective wrapping scheme due to the cost effectiveness. Finally the distribution of longitudinal and hoop stress has been determined by the finite element model for different diameters of CFRP strengthened VHS steel tube.

Explicit finite element modelling of bridge girder bearing pedestals

Bridge girder bearings rest on pedestals to transfer the loading safely to the pier headstock. In spite of the existence of industry guidelines, due to construction complexities, such guidelines are often overlooked. Further, there is paucity of research on the performance of pedestals, although their failure could cause exorbitant maintenance costs. Although reinforced concrete pedestals are recommended in the industry design guidelines, unreinforced concrete and/ or epoxy glue pedestals are provided due to construction issues; such pedestals fail within a very short period of service. With a view to understanding the response of pedestals subject to monotonic loading, a three-dimensional nonlinear explicit finite element micro-model of unreinforced and reinforced concrete pedestals has been developed. Contact and material nonlinearity have been accounted for in the model. It is shown that the unreinforced concrete pedestals suffer from localised edge stress singularities, the failure of which was comparable to those in the field. The reinforced concrete pedestals, on the other hand, distribute the loading without edge stress singularity, again conforming to the field experience.

Comparison of eight published lumbar spine finite element models

Due to its ability to represent intricate systems with material nonlinearities as well as irregular loading, boundary, geometrical and material domains, the finite element (FE) method has been recognized as an important computational tool in spinal biomechanics. Current FE models generally account for a single distinct spinal geometry with one set of material properties despite inherently large inter-subject variability. The uncertainty and high variability in tissue material properties, geometry, loading and boundary conditions has cast doubt on the reliability of their predictions and comparability with reported in vitro and in vivo values. A multicenter study was undertaken to compare the results of eight well-established models of the lumbar spine that have been developed, validated and applied for many years. Models were subjected to pure and combined loading modes and their predictions were compared to in vitro and in vivo measurements for intervertebral rotations, disc pressures and facet joint forces. Under pure moment loading, the predicted L1-5 rotations of almost all models fell within the reported in vitro ranges; their median values differed on average by only 2° for flexion-extension, 1° for lateral bending and 5° for axial rotation. Predicted median facet joint forces and disc pressures were also in good agreement with previously published median in vitro values. However, the ranges of predictions were larger and exceeded the in vitro ranges, especially for facet joint forces. For all combined loading modes, except for flexion, predicted median segmental intervertebral rotations and disc pressures were in good agreement with in vivo values. The simulations yielded median facet joint forces of 0 N in flexion, 38 N in extension, 14 N in lateral bending and 60 N in axial rotation that could not be validated due to the paucity of in vivo facet joint forces. In light of high inter-subject variability, one must be cautious when generalizing predictions obtained from one deterministic model. This study demonstrates however that the predictive power increases when FE models are combined together. The median of individual numerical results can hence be used as an improved tool in order to estimate the response of the lumbar spine.

Deformation of a three-dimensional red blood cell in a stenosed micro-capillary

Red blood cells (RBCs) exhibit different types of motions and deformations when the blood flows through capillaries. Interestingly, due to the complex three-dimensional structure of the RBC membrane, RBCs show three-dimensional motions and deformations in the blood flow. These motions and deformations of the RBCs highly depend on the stiffness of the RBC membrane and on the geometrical parameters of the capillary through which blood flows. However, capillaries always do not have uniform cross sections and some capillaries have stenosed segments, where cross sectional area suddenly reduces. Further, some diseases can alter the stiffness of the RBC membrane drastically. In this study, the deformation behaviour of a single three-dimensional RBC is examined, when it moves through a stenosed capillary. A three-dimensional spring network is used to model the RBC membrane. The RBC’s inside and outside fluids are discretized into a finite number of mass points and treated by smoothed particle hydrodynamics (SPH) method. The capillary is considered as a rigid tube with a stenosed section. The deformation index, mean velocity and total energy of the RBC are analysed when it flows through the stenosed capillary. Further, motion and deformation of the RBCs with different membrane stiffness (KB) are compared when they flow through the stenosed segment of the capillary. The simulation results demonstrate the RBCs are subjected to a larger deformation when they move through the stenosed part of the capillary and the RBCs with lower KBvalues easily pass through the stenosed segment of the capillary. Further, RBCs having higher KBvalues have a lower mean velocity and it leads to slow down the overall blood flow rate

Finite element modeling of carbon fibre reinforced polymer (CFRP) strengthened steel tubes under axial impact

Steel hollow sections used in structures such as bridges, buildings and space structures involve different strengthening techniques according to their structural purpose and shape of the structural member. One such technique is external bonding of CFRP sheets to steel tubes. The performance of CFRP strengthening for steel structures has been proven under static loading while limited studies have been conducted on their behaviour under impact loading. In this study, a comprehensive numerical investigation is carried out to evaluate the response of CFRP strengthened steel tubes under dynamic axial impact loading. Impact force, axial deformation impact velocities are studied. The results of the numerical investigations are validated by experimental results. Based on the developed finite element (FE) model several output parameters are discussed. The results show that CFRP wrapping is an effective strengthening technique to increase the axial dynamic load bearing capacity by increasing the stiffness of the steel tube.

Soft tissue human thigh and buttock finite element model to simulate vehicle seat cushion indentation

Accurate modelling of automotive occupant posture is strongly related to the mechanical interaction between human body soft tissue and flexible seat components. This paper presents a finite-element study simulating the deflection of seat cushion foam and supportive seat structures, as well as human buttock and thigh soft tissue when seated. The thigh-buttock surface shell model was based on 95th percentile male subject scan data and made of two layers, covering thin to moderate thigh and buttock proportions. To replicate the effects of skin and fat, the neoprene rubber layer was modelled as a hyperelastic material with viscoelastic behaviour. The analytical seat model is based on a Ford production seat. The result of the finite-element indentation simulation is compared to a previous simulation of an indentation with a hard shell human model of equal geometry, and to the physical indentation result. We conclude that SAE composite buttock form and human-seat indentation of a suspended seat cushion can be validly simulated.

Functional adaptation and rehabilitation following total knee replacement

This research project evaluated the biomechanical and functional outcomes of patients following total knee replacement measured at 6 and 12 months following surgery. Using more objective measures, patients were examined to determine changes in biomechanical and neuromuscular function during performance of activities of daily living such as walking, stair climbing and turning. Adaptations in joint positioning and performance were identified and progressive improvements were made in some areas of locomotor function. The findings of the study provided important objective information to contribute to the design and evaluation of prostheses, new surgical and rehabilitation procedures and improved recovery of patients.