Extramuscular myofascial force transmission: experiments and finite element modeling.

The specific purpose of the present study was to show that extramuscular myofascial force transmission exclusively has substantial effects on muscular mechanics. Muscle forces exerted at proximal and distal tendons of the rat extensor digitorium longus (EDL) were measured simultaneously, in two conditions (1) with intact extramuscular connections (2) after dissecting the muscles' extramuscular connections to a maximum extent without endangering circulation and innervation (as in most in situ muscle experiments). A finite element model of EDL including the muscles' extramuscular connections was used to assess the effects of extramuscular myofascial force transmission on muscular mechanics, primarily to test if such effects lead to distribution of length of sarcomeres within muscle fibers. In condition (1), EDL isometric forces measured at the distal and proximal tendons were significantly different (F(dist) > F(prox), DeltaF approximates maximally 40% of the proximal force). The model results show that extramuscular myofascial force transmission causes distributions of strain in the fiber direction (shortening in the proximal, lengthening in the distal ends of fibers) at higher lengths. This indicates significant length distributions of sarcomeres arranged in series within muscle fibers. Stress distributions found are in agreement with the higher distal force measured, meaning that the muscle fiber is no longer the unit exerting equal forces at both ends. Experimental results obtained in condition (2) showed no significant changes in the length-force characteristics (i.e., proximo-distal force differences were maintained). This shows that a muscle in situ has to be distinguished from a muscle that is truly isolated in which case the force difference has to be zero. We conclude that extramuscular myofascial force transmission has major effects on muscle functioning.

Design and analysis of coupling methods for modular endoprosthetic systems as an alternative to the conical coupling.

Modular endoprostheses are often used in bone tumor management. However, the conical coupling that connects the various modules has several shortcomings. As an alternative, four new couplings have been developed. To find out if they have sufficient strength and show no movement during loading, each coupling was analysed using the finite element method. Bolt force and friction coefficient was varied to examine their influence. From the analysis it was concluded that coupling B, a dovetail coupling, meets all requirements and is the best alternative to the conical coupling. Sensitivity to bolt force and friction coefficient is very limited.

A Hill type model of rat medial gastrocnemius muscle that accounts for shortening history effects.

The aim of the present study was to develop a Hill type muscle model that accounts for the effects of shortening history. For this purpose, a function was derived that relates force depression to starting length, shortening amplitude and contraction velocity. History parameters were determined from short-range isokinetic experiments on rat medial gastrocnemius muscle (GM). Simulations of isokinetic as well as isotonic experiments were performed with the new model and a standard Hill type model. The simulation results were compared with experimental results of rat GM to evaluate if incorporation of history effects leads to improvements in model predictions. In agreement with the experimental results, the new model qualitatively described force reduction during and after isokinetic shortening as well as the experimental observation that isometric endpoints of isotonic contractions are attained at higher muscle lengths than is expected from the fully isometric length-force curve. Consequently, the new model gave a better quantitative prediction of the experimental results compared to the standard model. It was concluded that incorporation of history effects can improve the predictive power of a Hill type model considerably. The applicability of the model to conditions other than those described in the present paper is discussed.

Modelling functional effects of muscle geometry.

Muscle architecture is an important aspect of muscle functioning. Hence, geometry and material properties of muscle have great influence on the force-length characteristics of muscle. We compared experimental results for the gastrocnemius medialis muscle (GM) of the rat to model results of simple geometric models such as a planimetric model and three-dimensional versions of this model. The capabilities of such models to adequately calculate muscle geometry and force-length characteristics were investigated. The planimetric model with elastic aponeurosis predicted GM muscle geometry well: maximal differences are 6, 1, 4 and 6% for fiber length, aponeurosis length, fiber angle and aponeurosis angle respectively. A slanted cylinder model with circular fiber cross-section did not predict muscle geometry as well as the planimetric model, whereas the geometry results of a second slanted cylinder model were identical to the planimetric model. It is concluded that the planimetric model is capable of adequately calculating the muscle geometry over the muscle length range studied. However, for modelling of force-length characteristics more complex models are needed, as none of the models yielded results sufficiently close to experimental data. Modelled force-length characteristics showed an overestimation of muscle optimum length by 2 mm with respect to experimental data, and the force at the ascending limb of the length force curve was underestimated. The models presented neglect important aspects such as non-linear geometry of muscle, certain passive material properties and mechanical interactions of fibers. These aspects may be responsible for short-comings in the modelling. It is argued that, considering the inability to adequately model muscle length-force characteristics for an isolated maximally activated (in situ) muscle, it is to be expected that prediction will fail for muscle properties in conditions of complex movement with many interacting factors. Therefore, modelling goals should be limited to the heuristic domain rather than expect to be able to predict or even approach medical or biological reality. However, the increased understanding about muscular mechanisms obtained from heuristic use of such simple models may very well be used in creating progress in, for example, clinical applications.

Revised planimetric model of unipennate skeletal muscle: a mechanical approach.

OBJECTIVE: Planimetric models which are simple, in the sense that small numerical effort is needed, are used to study functional consequences of skeletal muscle architecture. This paper argues with the approach to derive force of a unipennate muscle based on only equilibrium of the aponeurosis (tendon-sheet). In such an approach intramuscular pressure gradients are neglected and no suitable aponeurosis force can be determined. METHOD: The approach presented in this paper is based on mechanical equilibrium of whole muscle. A volume-related force is introduced to keep muscle volume constant. Mechanical equilibrium of whole muscle yields a different relation between fiber and muscle force as well as length changes as a consequence of pennation, compared with relations derived when only equilibrium of aponeurosis is considered. RESULTS: The newly derived relation improved prediction of the rat gastrocnemius medialis muscle force-length characteristics. CONCLUSION: The prediction of muscle geometry and the prediction of force-length characteristics are very good with a simple model such as a planimetric model. This conclusion suggests that the influence of properties neglected in such a simple model are either small or are internally compensated for in the net effects.

Numerical modelling of blood flow behaviour in the valved catheter of the PUCA-pump, a LVAD.

Mechanical heart assistance, performed by the PUlsatile CAtheter (PUCA) pump, chronologically takes place by sucking blood from the left ventricle and ejecting it into the ascending aorta. Within the pump activity the problem of hemolysis and clotting is encountered. In this study the influence of valve geometry on blood cell damage and stagnant zones has been investigated. A variable valve length coupled to a catheter ejection gap and a variable valve angle have been studied. In case of the studied valve, optimal parameter values have been determined. Compared to small catheter ejection gaps with a corresponding valve length, blood damage is found to be less for large ejection gaps with corresponding valve dimensions. In systole a valve positioned in a 0 degree angle proves to be best, whereas in diastole a +20 degree angle is preferable. Because the system is operating in both systole and diastole, a 0 degree angle valve is applied.

First clinical experience with a noninvasively extendable endoprosthesis: a limb-saving procedure in children suffering from a malignant bone tumor.

A modular endoprosthetic system that can be extended noninvasively has been applied for the first time in a growing child who underwent a tumor resection in his leg. The main goal of the study was to test the extendable mechanism that noninvasively corrects leg length differences caused by growth disturbances of the affected leg. The use of this endoprosthetic system resulted in good restoration of function. Six extensions were performed resulting in 19.5 mm of prosthetic growth. Unfortunately, an ingrown toenail caused infection of the endoprosthesis, and the infection necessitated extirpation of the prosthesis 15 months postoperatively. Two months later the patient died of acute leukemia. Analysis of the endoprosthesis revealed some manufacturing shortcomings, none of which impaired the function of the endoprosthesis.

An inverse dynamics model for the analysis, reconstruction and prediction of bipedal walking.

Walking is a constrained movement which may best be observed during the double stance phase when both feet contact the floor. When analyzing a measured movement with an inverse dynamics model, a violation of these constraints will always occur due to measuring errors and deviations of the segments model from reality, leading to inconsistent results. Consistency is obtained by implementing the constraints into the model. This makes it possible to combine the inverse dynamics model with optimization techniques in order to predict walking patterns or to reconstruct non-measured rotations when only a part of the three-dimensional joint rotations is measured. In this paper the outlines of the extended inverse dynamics method are presented, the constraints which define walking are defined and the optimization procedure is described. The model is applied to analyze a normal walking pattern of which only the hip, knee and ankle flexions/extensions are measured. This input movement is reconstructed to a kinematically and dynamically consistent three-dimensional movement, and the joint forces (including the ground reaction forces) and joint moments of force, needed to bring about this movement are estimated.

Effects of fit and bonding characteristics of femoral stems on adaptive bone remodeling.

Bone atrophy caused by stress-shielding may cause serious complications for the long-term fixation of hip stems. In particular, uncemented total hip arthroplasty is threatened by this problem, because the stems are usually larger and, as a consequence, stiffer than those of cemented implants. In the present study, the effects of fit and bonding characteristics of femoral hip stems were investigated, using the (nonlinear) finite element method in combination with adaptive bone remodeling theory to predict the bone density distribution in a bone or bone/implant configuration. Unknown parameters used in the theory, such as a reference equilibrium loading stimulus and a threshold (dead) zone of this stimulus, were established (triggered) by using the method to predict the density distributions in the natural femur and around fully coated uncemented implants. The computer simulation method can provide long term predictions of remodeling patterns around various implant configurations. Several cases were analyzed, whereby the coating conditions (fully, partly, or noncoated) and the fit characteristics (press fitted or overreamed) were varied. The computer predictions showed that partly coating can only significantly reduce bone atrophy relative to fully coated stems, when the coating is applied at a small region at the utmost proximal part of the stem. For smooth press-fit stems the predicted amount of bone loss (35 percent in the proximal medial region) was less than for a one-third proximally coated or a fully coated stem (50 to 54 percent predicted bone loss in the proximal medial region).(ABSTRACT TRUNCATED AT 250 WORDS)

Development and test of an extendable endoprosthesis for bone reconstruction in the leg.

A malignant bone tumour may develop in the femur of a child. In the majority of cases it will be necessary to resect the bone involved, growth plate and adjacent tissues. A modular endoprosthetic system has been developed which can be extended non-invasively to bridge the defect resulting from such a resection. Elongation is achieved by using an external magnetic field. In vitro tests with a prototype showed that the lengthening element met all requirements. Six animal experiments showed that the lengthening element also functioned in vivo.

Quantitative analysis of bone reactions to relative motions at implant-bone interfaces.

Connective soft tissues at the interface between implants and bone, such as in human joint replacements, can endanger the stability of the implant fixation. The potential of an implant to generate interface bone resorption and form soft tissue depends on many variables, including mechanical ones. These mechanical factors can be expressed in terms of relative motions between bone and implant at the interface or deformation of the interfacial material. The purpose of this investigation was to determine if interface debonding and subsequent relative interface motions can be responsible for interface degradation and soft tissue interposition as seen in experiments and clinical results. A finite element computer program was augmented with a mathematical description of interface debonding, dependent on interface stress criteria, and soft tissue interface interposition, dependent on relative interface motions. Three simplified models of orthopaedic implants were constructed: a cortical bone screw for fracture fixation plates, a femoral resurfacing prosthesis and a straight stem model, cemented in a bone. The predicted computer configurations were compared with clinical observations. The computer results showed how interface disruption and fibrous tissue interposition interrelate and possibly enhance each other, whereby a progressive development of the soft tissue layer can occur. Around the cortical bone screw, the predicted resorption patterns were relatively large directly under the screw head and showed a pivot point in the opposite cortex. The resurfacing cup model predicted some fibrous tissue formation under the medial and lateral cup rim, whereby the medial layer developed first because of higher initial interface stresses. The straight stem model predicted initial interface failure at the proximal parts. After proximal resorption and fibrous tissue interposition, the medial interface was completely disrupted and developed an interface layer. The distal and mid lateral side maintained within the strength criterion. Although the applied models were relatively simple, the results showed reasonable qualitative agreement with resorption patterns found in clinical studies concerning bone screws and the resurfacing cup. The hypothesis that interface debonding and subsequent relative (micro)motions could be responsible for bone resorption and fibrous tissue propagation is thereby sustained by the results.

The behavior of adaptive bone-remodeling simulation models.

The process of adaptive bone remodeling can be described mathematically and simulated in a computer model, integrated with the finite element method. In the model discussed here, cortical and trabecular bone are described as continuous materials with variable density. The remodeling rule applied to simulate the remodeling process in each element individually is, in fact, an objective function for an optimization process, relative to the external load. Its purpose is to obtain a constant, preset value for the strain energy per unit bone mass, by adapting the density. If an element in the structure cannot achieve that, it either turns to its maximal density (cortical bone) or resorbs completely. It is found that the solution obtained in generally a discontinuous patchwork. For a two-dimensional proximal femur model this patchwork shows a good resemblance with the density distribution of a real proximal femur. It is shown that the discontinuous end configuration is dictated by the nature of the differential equations describing the remodeling process. This process can be considered as a nonlinear dynamical system with many degrees of freedom, which behaves divergent relative to the objective, leading to many possible solutions. The precise solution is dependent on the parameters in the remodeling rule, the load and the initial conditions. The feedback mechanism in the process is self-enhancing, denser bone attracts more strain energy, whereby the bone becomes even more dense. It is suggested that this positive feedback of the attractor state (the strain energy field) creates order in the end configuration. In addition, the process ensures that the discontinuous end configuration is a structure with a relatively low mass, perhaps a minimal-mass structure, although this is no explicit objective in the optimization process. It is hypothesized that trabecular bone is a chaotically ordered structure which can be considered as a fractal with characteristics of optimal mechanical resistance and minimal mass, of which the actual morphology depends on the local (internal) loading characteristics, the sensor-cell density and the degree of mineralization.

Effects of material properties of femoral hip components on bone remodeling.

Bone loss around femoral hip stems is one of the problems threatening the long-term fixation of uncemented stems. Many believe that this phenomenon is caused by reduced stresses in the bone (stress shielding). In the present study the mechanical consequences of different femoral stem materials were investigated using adaptive bone remodeling theory in combination with the finite element method. Bone-remodeling in the femur around the implant and interface stresses between bone and implant were investigated for fully bonded femoral stems. Cemented stems (cobalt-chrome or titanium alloy) caused less bone resorption and lower interface stresses than uncemented stems made from the same materials. The range of the bone resorption predicted in the simulation models was from 23% in the proximal medial cortex surrounding the cemented titanium alloy stem to 76% in the proximal medial cortex around the uncemented cobalt-chrome stem. Very little bone resorption was predicted around a flexible, uncemented "iso-elastic" stem, but the proximal interface stresses increased drastically relative to the stiffer uncemented stems composed of cobalt-chrome or titanium alloy. However, the proximal interface stress peak was reduced and shifted during the adaptive remodeling process. The latter was found particularly in the stiffer uncemented cobalt-chrome-molybdenum implant and less for the flexible iso-elastic implant.

Articular contact in a three-dimensional model of the knee.

This study is aimed at the analysis of articular contact in a three-dimensional mathematical model of the human knee-joint. In particular the effect of articular contact on the passive motion characteristics is assessed in relation to experimentally obtained joint kinematics. Two basically different mathematical contact descriptions were compared for this purpose. One description was for rigid contact and one for deformable contact. The description of deformable contact is based on a simplified theory for contact of a thin elastic layer on a rigid foundation. The articular cartilage was described either as a linear elastic material or as a non-linear elastic material. The contact descriptions were introduced in a mathematical model of the knee. The locations of the ligament insertions and the geometry of the articular surfaces were obtained from a joint specimen of which experimentally determined kinematic data were available, and were used as input for the model. The ligaments were described by non-linear elastic line elements. The mechanical properties of the ligaments and the articular cartilage were derived from literature data. Parametric model evaluations showed that, relative to rigid articular contact, the incorporation of deformable contact did not alter the motion characteristics in a qualitative sense, and that the quantitative changes were small. Variation of the elasticity of the elastic layer revealed that decreasing the surface stiffness caused the ligaments to relax and, as a consequence, increased the joint laxity, particularly for axial rotation. The difference between the linear and the non-linear deformable contact in the knee model was very small for moderate loading conditions. The motion characteristics simulated with the knee model compared very well with the experiments. It is concluded that for simulation of the passive motion characteristics of the knee, the simplified description for contact of a thin linear elastic layer on a rigid foundation is a valid approach when aiming at the study of the motion characteristics for moderate loading conditions. With deformable contact in the knee model, geometric conformity between the surfaces can be modelled as opposed to rigid contact which assumed only point contact.

The cause of axial rotation of the scoliotic spine.

To explain the cause of axial rotation in a scoliotic vertebral column, the influence of the gravitation force on a spine with a C-scoliosis has been investigated by means of a mechanical model. In this model the gravitation force takes hold of the three-dimensionally curved vertebral column eccentrically. From these reflections it appears that the axial rotation in the scoliotic spine can be explained by the moment distribution caused by this eccentrical gravitation force. The moment distribution, necessary for correction of the spine, is supposed to be opposite to the moments caused by the gravitation force. The moment distribution caused by the Harrington and the Luque spinal correction systems are compared to the calculated optimum correction moments. It appears that the moment distribution for the Harrington and Luque methods, necessary for the correction of the lateral deviation, are almost the same as the calculated correction moments. But the axial rotation appears to be increasing instead of decreasing in both correction systems.

Stiffness and hysteresis properties of some prosthetic feet.

A prosthetic foot is an important element of a prosthesis, although it is not always fully recognized that the properties of the foot, along with the prosthetic knee joint and the socket, are in part responsible for the stability and metabolic energy cost during walking. The stiffness and the hysteresis, which are the topics of this paper, are not properly prescribed, but could be adapted to improve the prosthetic walking performance. The shape is strongly related to the cosmetic appearance and so can not be altered to effect these improvements. Because detailed comparable data on foot stiffness and hysteresis, which are necessary to quantify the differences between different types of feet, are absent in literature, these properties were measured by the authors in a laboratory setup for nine different prosthetic feet, bare and with two different shoes. One test cycle consisted of measurements of load deformation curves in 66 positions, representing the range from heel strike to toe-off. The hysteresis is defined by the energy loss as a part of the total deformation energy. Without shoes significant differences in hysteresis between the feet exist, while with sport shoes the differences in hysteresis between the feet vanish for the most part. Applying a leather shoe leads to an increase of hysteresis loss for all tested feet. The stiffness turned out to be non-constant, so mean stiffness is used.(ABSTRACT TRUNCATED AT 250 WORDS)

Accelerations due to impact at heel strike using below-knee prosthesis.

The acceleration in the sagittal plane of the prosthetic tube at heel strike in normal walking was measured in five healthy amputees with their definitive below-knee prosthesis, every subject using six different prosthetic feet, wearing sport shoes as well as leather shoes. The experiments were carried out in the rehabilitation centre "Het Roessingh", Enschede, The Netherlands. Maximum accelerations were extracted from the acceleration-time-signal. Mean acceleration maxima of all subjects were calculated for each foot-shoe combination to eliminate the individual influence of the subjects. In the axial direction the maximal accelerations demonstrate a clear difference among the prosthetic feet and the shoes, while in dorsoventral (tangential) direction the inter-individual variation in the acceleration extremes dominates the difference between the types of footwear. In comparison with non-amputees the magnitude of the maximal axial acceleration at heel strike does not differ significantly.

An extendable modular endoprosthetic system for bone tumour management in the leg.

A modular endoprosthetic system has been developed at the Groningen University Hospital and the University of Twente. The system can bridge a defect resulting from the resection of a malignant bone tumour which has developed around the knee joint of a child. Since the other healthy leg continues to grow, the system includes an element whose length can be adjusted non-invasively by using an external magnetic field. In addition to this lengthening element, there are one hip and two knee components, connectors of various lengths, and fixation elements. The paper describes the elements of the modular endoprosthetic system. Tables are created by means of which the elemental composition of such an endoprosthesis can be determined for each individual patient.

Trends of mechanical consequences and modeling of a fibrous membrane around femoral hip prostheses.

In the present study, the effects of a fibrous membrane between cement and bone in a femoral total hip replacement were investigated. The study involved the problem of modeling this fibrous membrane in finite-element analyses, and its global consequences for the load-transfer mechanism and its resulting stress patterns. A finite-element model was developed, suitable to describe nonlinear contact conditions in combination with nonlinear material properties of the fibrous membrane. The fibrous tissue layer was described as a highly compliant material with little resistance against tension and shear. The analysis showed that the load transfer mechanism from stem to bone changes drastically when such a membrane is present. These effects are predominantly caused by tensile loosening and slip at the interface, and are enhanced by the nonlinear membrane characteristics. Using parametric analysis, it was shown that these effects on the load-transfer mechanism cannot be described satisfactorily with linear elastic models. Most importantly, the fibrous tissue interposition causes excessive stress concentrations in bone and cement, and relatively high relative displacements between these materials.