Orthotropic active strain models for the numerical simulation of cardiac biomechanics

An active strain formulation for orthotropic constitutive laws arising in cardiac mechanics modeling is introduced and studied. The passive mechanical properties of the tissue are described by the Holzapfel-Ogden relation. In the active strain formulation, the Euler-Lagrange equations for minimizing the total energy are written in terms of active and passive deformation factors, where the active part is assumed to depend, at the cell level, on the electrodynamics and on the specific orientation of the cardiac cells. The well-posedness of the linear system derived from a generic Newton iteration of the original problem is analyzed and different mechanical activation functions are considered. In addition, the active strain formulation is compared with the classical active stress formulation from both numerical and modeling perspectives. Taylor-Hood and MINI finite elements are employed to discretize the mechanical problem. The results of several numerical experiments show that the proposed formulation is mathematically consistent and is able to represent the main key features of the phenomenon, while allowing savings in computational costs.

High-precision satellite positioning system as a new tool to study the biomechanics of human locomotion.

New Global Positioning System (GPS) receivers allow now to measure a location on earth at high frequency (5Hz) with a centimetric precision using phase differential positioning method. We studied whether such technique was accurate enough to retrieve basic parameters of human locomotion. Eight subjects walked on an athletics track at four different imposed step frequencies (70-130steps/min) plus a run at free pace. Differential carrier phase localization between a fixed base station and the mobile antenna mounted on the walking person was calculated. In parallel, a triaxial accelerometer, attached to the low back, recorded body accelerations. The different parameters were averaged for 150 consecutive steps of each run for each subject (total of 6000 steps analyzed). We observed a perfect correlation between average step duration measured by accelerometer and by GPS (r=0.9998, N=40). Two important parameters for the calculation of the external work of walking were also analyzed, namely the vertical lift of the trunk and the velocity variation per step. For an average walking speed of 4.0km/h, average vertical lift and velocity variation were, respectively, 4.8cm and 0.60km/h. The average intra-individual step-to-step variability at a constant speed, which includes GPS errors and the biological gait style variation, were found to be 24. 5% (coefficient of variation) for vertical lift and 44.5% for velocity variation. It is concluded that GPS technique can provide useful biomechanical parameters for the analysis of an unlimited number of strides in an unconstrained free-living environment.

Wireless GPS-based phase-locked synchronization system for outdoor environment.

Synchronization of data coming from different sources is of high importance in biomechanics to ensure reliable analyses. This synchronization can either be performed through hardware to obtain perfect matching of data, or post-processed digitally. Hardware synchronization can be achieved using trigger cables connecting different devices in many situations; however, this is often impractical, and sometimes impossible in outdoors situations. The aim of this paper is to describe a wireless system for outdoor use, allowing synchronization of different types of - potentially embedded and moving - devices. In this system, each synchronization device is composed of: (i) a GPS receiver (used as time reference), (ii) a radio transmitter, and (iii) a microcontroller. These components are used to provide synchronized trigger signals at the desired frequency to the measurement device connected. The synchronization devices communicate wirelessly, are very lightweight, battery-operated and thus very easy to set up. They are adaptable to every measurement device equipped with either trigger input or recording channel. The accuracy of the system was validated using an oscilloscope. The mean synchronization error was found to be 0.39 μs and pulses are generated with an accuracy of <2 μs. The system provides synchronization accuracy about two orders of magnitude better than commonly used post-processing methods, and does not suffer from any drift in trigger generation.

Predicting trabecular bone microdamage initiation and accumulation using a non-linear perfect damage model

Studies evaluating the mechanical behavior of the trabecular microstructure play an important role in our understanding of pathologies such as osteoporosis, and in increasing our understanding of bone fracture and bone adaptation. Understanding of such behavior in bone is important for predicting and providing early treatment of fractures. The objective of this study is to present a numerical model for studying the initiation and accumulation of trabecular bone microdamage in both the pre- and post-yield regions. A sub-region of human vertebral trabecular bone was analyzed using a uniformly loaded anatomically accurate microstructural three-dimensional finite element model. The evolution of trabecular bone microdamage was governed using a non-linear, modulus reduction, perfect damage approach derived from a generalized plasticity stress-strain law. The model introduced in this paper establishes a history of microdamage evolution in both the pre- and post-yield regions

Superficial and deep changes of cellular mechanical properties following cytoskeleton disassembly.

The cytoskeleton, composed of actin filaments, intermediate filaments, and microtubules, is a highly dynamic supramolecular network actively involved in many essential biological mechanisms such as cellular structure, transport, movements, differentiation, and signaling. As a first step to characterize the biophysical changes associated with cytoskeleton functions, we have developed finite elements models of the organization of the cell that has allowed us to interpret atomic force microscopy (AFM) data at a higher resolution than that in previous work. Thus, by assuming that living cells behave mechanically as multilayered structures, we have been able to identify superficial and deep effects that could be related to actin and microtubule disassembly, respectively. In Cos-7 cells, actin destabilization with Cytochalasin D induced a decrease of the visco-elasticity close to the membrane surface, while destabilizing microtubules with Nocodazole produced a stiffness decrease only in deeper parts of the cell. In both cases, these effects were reversible. Cell softening was measurable with AFM at concentrations of the destabilizing agents that did not induce detectable effects on the cytoskeleton network when viewing the cells with fluorescent confocal microscopy. All experimental results could be simulated by our models. This technology opens the door to the study of the biophysical properties of signaling domains extending from the cell surface to deeper parts of the cell.

Determination of displacement, stress- and strain-distribution in the human heart: a FE-model on the basis of MR imaging.

The determination of characteristic cardiac parameters, such as displacement, stress and strain distribution are essential for an understanding of the mechanics of the heart. The calculation of these parameters has been limited until recently by the use of idealised mathematical representations of biventricular geometries and by applying simple material laws. On the basis of 20 short axis heart slices and in consideration of linear and nonlinear material behaviour we have developed a FE model with about 100,000 degrees of freedom. Marching Cubes and Phong's incremental shading technique were used to visualise the three dimensional geometry. In a quasistatic FE analysis continuous distribution of regional stress and strain corresponding to the endsystolic state were calculated. Substantial regional variation of the Von Mises stress and the total strain energy were observed at all levels of the heart model. The results of both the linear elastic model and the model with a nonlinear material description (Mooney-Rivlin) were compared. While the stress distribution and peak stress values were found to be comparable, the displacement vectors obtained with the nonlinear model were generally higher in comparison with the linear elastic case indicating the need to include nonlinear effects.

Translation of biomechanical concepts in bone tissue engineering: from animal study to revision knee arthroplasty.

Bone defects in revision knee arthroplasty are often located in load-bearing regions. The goal of this study was to determine whether a physiologic load could be used as an in situ osteogenic signal to the scaffolds filling the bone defects. In order to answer this question, we proposed a novel translation procedure having four steps: (1) determining the mechanical stimulus using finite element method, (2) designing an animal study to measure bone formation spatially and temporally using micro-CT imaging in the scaffold subjected to the estimated mechanical stimulus, (3) identifying bone formation parameters for the loaded and non-loaded cases appearing in a recently developed mathematical model for bone formation in the scaffold and (4) estimating the stiffness and the bone formation in the bone-scaffold construct. With this procedure, we estimated that after 3 years mechanical stimulation increases the bone volume fraction and the stiffness of scaffold by 1.5- and 2.7-fold, respectively, compared to a non-loaded situation.

A musculoskeletal shoulder model based on pseudo-inverse and null-space optimization.

The goal of the present work was assess the feasibility of using a pseudo-inverse and null-space optimization approach in the modeling of the shoulder biomechanics. The method was applied to a simplified musculoskeletal shoulder model. The mechanical system consisted in the arm, and the external forces were the arm weight, 6 scapulo-humeral muscles and the reaction at the glenohumeral joint, which was considered as a spherical joint. The muscle wrapping was considered around the humeral head assumed spherical. The dynamical equations were solved in a Lagrangian approach. The mathematical redundancy of the mechanical system was solved in two steps: a pseudo-inverse optimization to minimize the square of the muscle stress and a null-space optimization to restrict the muscle force to physiological limits. Several movements were simulated. The mathematical and numerical aspects of the constrained redundancy problem were efficiently solved by the proposed method. The prediction of muscle moment arms was consistent with cadaveric measurements and the joint reaction force was consistent with in vivo measurements. This preliminary work demonstrated that the developed algorithm has a great potential for more complex musculoskeletal modeling of the shoulder joint. In particular it could be further applied to a non-spherical joint model, allowing for the natural translation of the humeral head in the glenoid fossa.

Diagnostic et traitement de l'instabilité rotulienne [Diagnosis and treatment of patellar instability].

Patellar instability concerns mostly young and active patients. Because of the condition's multifactorial origin, knowledge about the risk factors of patella instability and the biomechanics of the femoropatellar joint are necessary to understand this disorder. Most of the patellar dislocations are spontaneously reduced. Only the patient's history, clinical examination and appropriate imaging lead to the correct diagnosis. If haemarthrosis is present, MRI is recommended to rule out concomitant injuries. A traumatic patellar dislocation without history of prior dislocation is treated conservatively. Surgical treatment is indicated in presence of intra-articular injury or for recurrent dislocations. The surgical procedure should correct the different factors of patella instability identified by preoperative analysis.

Mechanical external work and recovery at preferred walking speed in obese subjects

PURPOSE: The aim of this study was to compare the mechanical external work (per kg) and pendular energy transduction at preferred walking speed (PWS) in obese versus normal body mass subjects to investigate whether obese adults adopt energy conserving gait mechanics. METHODS: The mechanical external work (Wext) and the fraction of mechanical energy recovered by the pendular mechanism (Rstep) were computed using kinematic data acquired by an optoelectronic system and were compared in 30 obese (OG; body mass index [BMI] = 39.6 +/- 0.6 kg m(-2); 29.5 +/- 1.3 yr) and 19 normal body mass adults (NG; BMI = 21.4 +/- 0.5 kg m(-2); 31.2 +/- 1.2 yr) walking at PWS. RESULTS: PWS was significantly lower in OG (1.18 +/- 0.02 m s(-1)) than in NG (1.33 +/- 0.02 m s(-1); P <or= 0.001). There was no significant difference in Wext per unit mass between groups (OG: 0.36 +/- 0.03 J kg(-1) m(-1); NG: 0.31 +/- 0.02 J kg(-1) m(-1); P = 0.12). Rstep was significantly lower in OG (68.4% +/- 2.0%) compared with NG (74.4% +/- 1.0%; P = 0.01). In OG only, Wext per unit mass was positively correlated with PWS (r = 0.57; P < 0.001). CONCLUSION: Obese adults do not appear to alter their gait to improve pendular energy transduction and may select slower PWS to reduce mechanical and metabolic work.

Paclitaxel-eluting biodegradable synthetic vascular prostheses: a step towards reduction of neointima formation?

BACKGROUND: Clinical small-caliber vascular prostheses are unsatisfactory. Reasons for failure are early thrombosis and late intimal hyperplasia. We thus prepared biodegradable small-caliber vascular prostheses using electrospun polycaprolactone (PCL) with slow-releasing paclitaxel (PTX), an antiproliferative drug. METHODS AND RESULTS: PCL solutions containing PTX were used to prepare nonwoven nanofibre-based 2-mm ID prostheses. Mechanical morphological properties and drug loading, distribution, and release were studied in vitro. Infrarenal abdominal aortic replacement was carried out with nondrug-loaded and drug-loaded prostheses in 18 rats and followed for 6 months. Patency, stenosis, tissue reaction, and drug effect on endothelialization, vascular remodeling, and neointima formation were studied in vivo. In vitro prostheses showed controlled morphology mimicking extracellular matrix with mechanical properties similar to those of native vessels. PTX-loaded grafts with suitable mechanical properties and controlled drug-release were obtained by factorial design. In vivo, both groups showed 100% patency, no stenosis, and no aneurysmal dilatation. Endothelial coverage and cell ingrowth were significantly reduced at 3 weeks and delayed at 12 and 24 weeks in PTX grafts, but as envisioned, neointima formation was significantly reduced in these grafts at 12 weeks and delayed at 6 months. CONCLUSIONS: Biodegradable, electrospun, nanofibre, polycaprolactone prostheses are promising because in vitro they maintain their mechanical properties (regardless of PTX loading), and in vivo show good patency, reendothelialize, and remodel with autologous cells. PTX loading delays endothelialization and cellular ingrowth. Conversely, it reduces neointima formation until the end point of our study and thus may be an interesting option for small caliber vascular grafts.

3D gait assessment in young and elderly subjects using foot-worn inertial sensors.

This study describes the validation of a new wearable system for assessment of 3D spatial parameters of gait. The new method is based on the detection of temporal parameters, coupled to optimized fusion and de-drifted integration of inertial signals. Composed of two wirelesses inertial modules attached on feet, the system provides stride length, stride velocity, foot clearance, and turning angle parameters at each gait cycle, based on the computation of 3D foot kinematics. Accuracy and precision of the proposed system were compared to an optical motion capture system as reference. Its repeatability across measurements (test-retest reliability) was also evaluated. Measurements were performed in 10 young (mean age 26.1±2.8 years) and 10 elderly volunteers (mean age 71.6±4.6 years) who were asked to perform U-shaped and 8-shaped walking trials, and then a 6-min walking test (6MWT). A total of 974 gait cycles were used to compare gait parameters with the reference system. Mean accuracy±precision was 1.5±6.8cm for stride length, 1.4±5.6cm/s for stride velocity, 1.9±2.0cm for foot clearance, and 1.6±6.1° for turning angle. Difference in gait performance was observed between young and elderly volunteers during the 6MWT particularly in foot clearance. The proposed method allows to analyze various aspects of gait, including turns, gait initiation and termination, or inter-cycle variability. The system is lightweight, easy to wear and use, and suitable for clinical application requiring objective evaluation of gait outside of the lab environment.

Plantar pressures in the tennis serve.

In-shoe loading patterns were examined in each foot (back and front) separately during two types of tennis serve [first (or flat) and second (or twist) serve] and two service stance styles [foot-up (back foot is moved forward next to front foot for push-off) and foot-back (feet remain at the same relative level)]. Ten competitive tennis players completed five trials for each type of serve and service stance style in random order. Plantar pressure distribution was recorded using Pedar insoles divided into nine areas for analysis. Mean and peak pressures (+15.2%, P < 0.01 and +12.8%, P < 0.05) as well as maximal forces (+20.2%, P < 0.01) were higher under the lateral forefoot of the front foot in first than in second serves, while mean forces were higher (+17.2%, P < 0.05) under the lesser toes. Relative load was higher on the lateral forefoot (+20.4%, P < 0.05) but lower (-32.5%, P < 0.05) on the medial heel of the front foot with foot-up compared with foot-back stance. Using a foot-up stance, loading of the back foot was higher (+31.8%, P < 0.01) under the lateral mid-foot but lower (-29.9%, P < 0.01) under the medial forefoot. The type of serve and the stance style adopted have a significant effect on foot loading. Such findings might help improve mechanical efficiency of the serve.