Assessment of mechanical conditions in sub-dermal tissues during sitting: a combined experimental-MRI and finite element approach

A common but potentially severe malady afflicting permanent wheelchair users is pressure sores caused by elevated soft tissue strains and stresses over a critical prolonged period of time. Presently, there is paucity of information regarding deep soft tissue strains and stresses in the buttocks of humans during sitting. Strain and stress distributions in deep muscle and fat tissues were therefore calculated in six healthy subjects during sitting, in a double-donut Open-MR system, using a "reverse engineering" approach. Specifically, finite element (FE) models of the undeformed buttock were built for each subject using MR images taken at the coronal plane in a non-weight-bearing sitting posture. Using a second MR image taken from each subject during weight-bearing sitting we characterized the ischial tuberosity sagging toward the sitting surface in weight-bearing, and used these data as displacement boundary conditions for the FE models. These subject-specific FE analyses showed that maximal tissue strains and stresses occur in the gluteal muscles, not in fat or at the skin near the body-seat interface. Peak principal compressive strain and stress in the gluteus muscle were 74+/-7% and 32+/-9 kPa (mean+/-standard deviation), respectively. Peak principal compressive strain and stress in enveloping fat tissue were 46+/-7% and 18+/-4 kPa, respectively. Models were validated by comparing measured peak interface pressures under the ischial tuberosities (17+/-4 kPa) with those calculated by means of FE (18+/-3 kPa), for each subject. This is the first study to quantify sub-dermal tissue strain and stress distributions in sitting humans, in vivo. These data are essential for understanding the aetiology of pressure sores, particularly those that were recently termed "deep tissue injury" at the US National Pressure Ulcer Advisory Panel (NPUAP) 2005 Consensus Conference.     

The variability and complexity of sitting postural control are associated with discomfort

The present investigation examined the variability of sitting postural movement in relation to the development of perceived discomfort by means of linear and nonlinear analysis. Nine male subjects participated in this study. Discomfort ratings, kinetic and kinematics data were recorded during prolonged sitting. Body part discomfort index, displacement of the center of pressure (COP) in anterior–posterior and medial–lateral directions as well as lumbar curvature were calculated. Mean, standard deviation and sample entropy values were extracted from COP and lumbar curvature signals. Standard deviation and sample entropy were used to assess the degree of variability and complexity of sitting. A correlation analysis was performed to determine the correlation of each parameter with discomfort. There were no correlations between discomfort and any of the mean values. On the contrary, the standard deviations of the COP displacement in both directions and lumbar curvature were positively correlated to discomfort, whereas sample entropies were negatively correlated. The present study suggests that the increase in degree of variability and the decrease in complexity of sitting postural control are interrelated with the increase in perceived discomfort. Finally, the present study underlined the importance of quantifying motor variability for understanding the biomechanics of seated posture.     

Dynamics of seated computer work before and after prolonged constrained sitting

This laboratory study investigated seated computer work before and after prolonged constrained sitting. Discomfort ratings and kinetic and kinematics data were recorded in nine healthy males performing computer work for 5 min before and after 96 min of sitting. The displacement of the center of pressure (CoP) in anterior-posterior and medial-lateral directions and lumbar curvature (LC) were calculated. The root mean square, standard deviation, and sample entropy values were computed from the CoPs and LC signals to assess the magnitude, amount of variability, and regularity of sitting dynamics, respectively. The discomfort increased for the buttocks (p = .02).The standard deviation and sample entropy values of the CoPs and LC signals, respectively, increased (p < .04) and decreased (p < .004) whereas the root mean square remained unchanged (p > .15) after prolonged constrained sitting compared with before. This present study showed that during seated computer work, prolonged constrained sitting affected the amount of variability and the regularity of sitting postural control, whereas the magnitude was not affected. The importance of the dynamics of sitting control may challenge the idea of a static and ideal seated posture at work.     

Comfort evaluation as the example of anthropotechnical furniture design

Human health is becoming an increasingly important issue in contemporary hectic lifestyle imposed at work and by struggle to save time and money. Sitting comfort and quality of chairs which we use for the most of our time have, thus, become essential for healthy lifestyle. Sitting discomforts arise from prolonged sitting on the inappropriate chairs, which failing to provide sufficient support to the body cause discomfort and tiring. The studies of the office chair constructions have identified differences in perception of comfort provided by different types of seats. Four seat constructions and the comfort they provide to the sitters were compared by means of subjective indicators. After a two-day sitting on each of the studied chairs the subjects scored their perception of comfort and discomfort, using the questionnaire with 17 statements. Constructional forms and materials which contributed more to the sense of comfort by minimizing fatigue and pains developed by sitting were determined.     

Hip joint degeneration due to cam impingement: a finite element analysis

The goal of this study was to investigate the impact of cam impingement, a biomechanical risk factor, on hip joint degeneration and ultimately coxarthrosis. 3D finite element solid models of a healthy and a pathologic hip were developed based on clinical reports. The biphasic characteristics of cartilaginous tissues were considered to identify localised solid matrix overloading during normal walking and sitting down (SD). Localised femoral intrusion at the anterior-superior pelvic horn was revealed in the pathologic hip during SD, where the radial and meridional solid stresses in the acetabular cartilage and circumferential solid stresses within the acetabular labrum increased by 3.7, 1.5 and 2.7 times, respectively. The increased solid-on-solid stresses, reduction in fluid-load support and associated higher friction during articulation may result in joint wear and other degenerative changes in the hip.     

Heterogeneous linear elastic trabecular bone modelling using micro-CT attenuation data and experimentally measured heterogeneous tissue properties

High-resolution voxel-based finite element software, such as FEEBE developed at the NCBES, is widely used for studying trabecular bone at the micro-scale. A new approach to determine heterogeneous bone tissue material properties for computational models was proposed in this study. The specimen-specific range of tissue moduli across strut width was determined from nanoindentation testing. This range was mapped directly using linear interpolation to that specimen's micro-computed tomography (microCT) grey value range as input material properties for finite element analysis. The method was applied to cuboid trabecular bone samples taken from eight, 4-year-old (skeletally mature) ovine L5 vertebrae. Before undergoing experimental uniaxial compression tests, the samples were microCT scanned and 30 microm resolution finite element models were generated. The linear elastic finite element models were compressed to 1% strain. This material property assignment method for computational models accurately reproduced the experimentally determined apparent modulus and concentrations of stress at locations of failure.     

Mechanical compression-induced pressure sores in rat hindlimb: muscle stiffness, histology, and computational models.

Pressure sores affecting muscles are severe injuries associated with ischemia, impaired metabolic activity, excessive tissue deformation, and insufficient lymph drainage caused by prolonged and intensive mechanical loads. We hypothesize that mechanical properties of muscle tissue change as a result of exposure to prolonged and intensive loads. Such changes may affect the distribution of stresses in soft tissues under bony prominences and potentially expose additional uninjured regions of muscle tissue to intensified stresses. In this study, we characterized changes in tangent elastic moduli and strain energy densities of rat gracilis muscles exposed to pressure in vivo (11.5, 35, or 70 kPa for 2, 4, or 6 h) and incorporated the abnormal properties that were measured in finite element models of the head, shoulders, pelvis, and heels of a recumbent patient. Using in vitro uniaxial tension testing, we found that tangent elastic moduli of muscles exposed to 35 and 70 kPa were 1.6-fold those of controls (P < 0.05, for strains /=5%). Histological (phosphotungstic acid hematoxylin) evaluation showed that this stiffening accompanied extensive necrotic damage. Incorporating these effects into the finite element models, we were able to show that the increased muscle stiffness in widening regions results in elevated tissue stresses that exacerbate the potential for tissue necrosis. Interfacial pressures could not predict deep muscle (e.g., longissimus or gluteus) stresses and injuring conditions. We conclude that information on internal muscle stresses is required to establish new criteria for pressure sore prevention.     

Effects of a posture-sensing air seat device (PSASD) on kinematics and trunk muscle activity during continuous computer work

The prevalence and incidence of musculoskeletal disorders is high with computer workers, and poor sitting posture can be considered a factor contributing to low back discomfort. In the clinical literature, maintaining a neutral spinal curvature has been considered an optimal sitting posture. This study investigated the flexion and lateral flexion of trunk movements and trunk muscle activity during computer work with and without a posture-sensing air seat device (PSASD). By sensing a certain amount of increased pressure over the baseline, posture-related visual feedback was given to participants through the PSASD. Eleven regular computer workers participated in this study. PSASD had the function of alerting the subject to their poor posture by using visual feedback. Subjects performed 20 min of computer work with and without a PSASD. Surface electromyography was used to measure the activity of the erector spine and internal abdominal oblique. Kinematic data were obtained using an electrogoniometer. The results showed that the mean of trunk flexion and lateral flexion was significantly reduced with PSASD. The activity of the erector spine and internal oblique was significantly higher with the PSASD than without. Our findings indicated that the PSASD helps to prevent habitual poor posture by maintaining an erect sitting posture during prolonged computer work.     

Effects of sitting postures on risks for deep tissue injury in the residuum of a transtibial prosthetic-user: a biomechanical case study

Transtibial amputation prosthetic-users are at risk of developing deep tissue injury (DTI) while donning their prosthesis for prolonged periods; however, no study addresses the mechanical loading of the residuum during sitting with a prosthesis. We combined MRI-based 3D finite element modelling of a residuum with an injury threshold and a muscle damage law to study risks for DTI in one sitting subject in two postures: 30°-knee-flexion vs. 90°-knee-flexion. We recorded skin-socket pressures, used as model boundary conditions. During the 90°-knee-flexion simulations, major internal muscle injuries were predicted (>1000 mm³). In contrast, the 30°-knee-flexion simulations only produced minor injury ( < 14 mm³). Predicted injury rates at 90°-knee-flexion were over one order of magnitude higher than those at 30°-knee-flexion. We concluded that in this particular subject, prolonged 90°-knee-flexion sitting theoretically endangers muscle viability in the residuum. By expanding the studies to large subject groups, this research approach can support development of guidelines for DTI prevention in prosthetic-users.     

A new method to investigate in vivo knee behavior using a finite element model of the lower limb

Several finite element models have been developed for estimating the mechanical response of joint internal structures, where direct or indirect in vivo measurement is difficult or impossible. The quality of the predictions made by those models is largely dependent on the quality of the experimental data (e.g. load/displacement) used to drive them. Also numerical problems have been described in the literature when using implicit finite element techniques to simulate problems that involve contacts and large displacements. In this study, a unique strategy was developed combining high accuracy in vivo three-dimensional kinematics and a lower limb finite element model based on explicit finite element techniques. The method presents an analytical technique applied to a dynamic loading condition (impact during hopping on one leg). The validation of the lower limb model focused on the response of the whole model and the knee joint in particular to the imposed 3D femoral in vivo kinematics and ground reaction forces. The approach outlined in this study introduces a generic tool for the study of in vivo knee joint behavior.     

Computational and experimental models of the human torso for non-penetrating ballistic impact

Both computational finite element and experimental models of the human torso have been developed for ballistic impact testing. The human torso finite element model (HTFEM), including the thoracic skeletal structure and organs, was created in the finite element code LS-DYNA. The skeletal structure was assumed to be linear-elastic while all internal organs were modeled as viscoelastic. A physical human surrogate torso model (HSTM) was developed using biosimulant materials and the same anthropometry as the HTFEM. The HSTM response to impact was recorded with piezoresistive pressure sensors molded into the heart, liver and stomach and an accelerometer attached to the sternum. For experimentation, the HSTM was outfitted with National Institute of Justice (NIJ) Level I, IIa, II and IIIa soft armor vests. Twenty-six ballistic tests targeting the HSTM heart and liver were conducted with 22 caliber ammunition at a velocity of 329 m/s and 9 mm ammunition at velocities of 332, 358 and 430 m/s. The HSTM pressure response repeatability was found to vary by less than 10% for similar impact conditions. A comparison of the HSTM and HTFEM response showed similar pressure profiles and less than 35% peak pressure difference for organs near the ballistic impact point. Furthermore, the peak sternum accelerations of the HSTM and HTFEM varied by less than 10% for impacts over the sternum. These models provide comparative tools for determining the thoracic response to ballistic impact and could be used to evaluate soft body armor design and efficacy, determine thoracic injury mechanisms and assist with injury prevention.     

Search for critical loading condition of the spine--a meta analysis of a nonlinear viscoelastic finite element model

The relative vulnerability of spinal motion segments to different loading combinations remains unknown. The meta-analysis described here using the results of a validated L2-L3 nonlinear viscoelastic finite element model was designed to investigate the critical loading and its effect on the internal mechanics of the human lumbar spine. A Box-Behnken experimental design was used to design the magnitude of seven independent variables associated with loads, rotations and velocity of motion. Subsequently, an optimization method was used to find the primary and secondary variables that influence spine mechanical output related to facet forces, disc pressure, ligament forces, annulus matrix compressive/shear stresses and anulus fibers strain. The mechanical responses with respect to the two most-relevant variables were then regressed linearly using the response surface quadratic model. Axial force and sagittal rotation were identified as the most-relevant variables for mechanical responses. The procedure developed can be used to find the critical loading for finite element models with multi input variables. The derived meta-models can be used to predict the risk associated with various loading parameters and in setting safer load limits.     

Search for critical loading condition of the spine–A meta analysis of a nonlinear viscoelastic finite element model

The relative vulnerability of spinal motion segments to different loading combinations remains unknown. The meta-analysis described here using the results of a validated L2–L3 nonlinear viscoelastic finite element model was designed to investigate the critical loading and its effect on the internal mechanics of the human lumbar spine. A Box-Behnken experimental design was used to design the magnitude of seven independent variables associated with loads, rotations and velocity of motion. Subsequently, an optimization method was used to find the primary and secondary variables that influence spine mechanical output related to facet forces, disc pressure, ligament forces, annulus matrix compressive/shear stresses and anulus fibers strain. The mechanical responses with respect to the two most-relevant variables were then regressed linearly using the response surface quadratic model. Axial force and sagittal rotation were identified as the most-relevant variables for mechanical responses. The procedure developed can be used to find the critical loading for finite element models with multi input variables. The derived meta-models can be used to predict the risk associated with various loading parameters and in setting safer load limits.     

Effects of thresholding techniques on microCT-based finite element models of trabecular bone

Microimaging based finite element analysis is widely used to predict the mechanical properties of trabecular bone. The choice of thresholding technique, a necessary step in converting grayscale images to finite element models, can significantly influence the predicted bone volume fraction and mechanical properties. Therefore, we investigated the effects of thresholding techniques on microcomputed tomography (micro-CT) based finite element models of trabecular bone. Three types of thresholding techniques were applied to 16-bit micro-CT images of trabecular bone to create three different models per specimen. Bone volume fractions and apparent moduli were predicted and compared to experimental results. In addition, trabecular tissue mechanical parameters and morphological parameters were compared among different models. Our findings suggest that predictions of apparent mechanical properties and structural properties agree well with experimental measurements regardless of the choice of thresholding methods or the format of micro-CT images.     

Numerical investigation of the intravascular coronary stent flexibility

Nowadays stent therapy is widely adopted to treat atherosclerotic vessel diseases. The high commercial value of these devices and the high prototypation costs require the use of finite element analyses, instead of classical trial and error technique, to design and verify new models. In this paper, we explore the advantages of the finite element method (FEM) in order to investigate new generation stent performance in terms of flexibility. Indeed, the ability of the stent to bend in order to accommodate curvatures and angles of vessels during delivery is one of the most significant prerequisites for optimal stent performance. Two different FEM models, resembling two new generation intravascular stents, were developed. The main model dimensions were obtained by means of a stereo microscope, analyzing one Cordis BX-Velocity and one Carbostent Sirius coronary stent. Bending tests under displacement control in the unexpanded and expanded configuration were carried out. The curvature index, defined as the ratio between the sum of rotation angles at the extremes and the length of the stent, yielded comparative information about the capability of the device to be delivered into tortuous vessels and to conform to their contours. Results, expressed in terms of the bending moment-curvature index, demonstrated a different response for the two models. In particular the Cordis model showed a higher flexibility. Lower flexibility in the expanded configurations for both models was detected. However this flexibility depends on how the contact takes place between the different parts of the struts.     

Simulation-based assessment for construction helmets

In recent years, there has been a concerted effort for greater job safety in all industries. Personnel protective equipment (PPE) has been developed to help mitigate the risk of injury to humans that might be exposed to hazardous situations. The human head is the most vulnerable to impact as a moderate magnitude can cause serious injury or death. That is why industries have required the use of an industrial hard hat or helmet. There have only been a few articles published to date that are focused on the risk of head injury when wearing an industrial helmet. A full understanding of the effectiveness of construction helmets on reducing injury is lacking. This paper presents a simulation-based method to determine the threshold at which a human will sustain injury when wearing a construction helmet and assesses the risk of injury for wearers of construction helmets or hard hats. Advanced finite element, or FE, models were developed to study the impact on construction helmets. The FE model consists of two parts: the helmet and the human models. The human model consists of a brain, enclosed by a skull and an outer layer of skin. The level and probability of injury to the head was determined using both the head injury criterion (HIC) and tolerance limits set by Deck and Willinger. The HIC has been widely used to assess the likelihood of head injury in vehicles. The tolerance levels proposed by Deck and Willinger are more suited for finite element models but lack wide-scale validation. Different cases of impact were studied using LSTC's LS-DYNA.     

Determination of stress levels and profiles in the periodontal ligament by means of an improved three-dimensional finite element model for various types of orthodontic and natural force systems

This study was undertaken to investigate the stress-strain levels and distribution within the periodontal ligament for various types of physiological and orthodontic force systems, assuming that the bone resorption process, leading to tooth movements, is partly controlled by those conditions. Two finite element models were developed, simulating a full and partial mandibular morphology, respectively. Both models were based on morphology and physical parameters of human autopsy material. The effect of changing material parameters and structure, type of boundary conditions, calculation method and fineness of the model on the stress levels and profiles in the periodontal ligament was evaluated by a series of tests. A structure optimization technique was used to investigate the load bearing characteristics of the mandible and the influence of the anisotropic material properties of both the mandible and the segment. A ‘multiple modelling’ technique based on both the mandible and the segment was developed to test various types of boundary conditions in the analysis of the segment. Results presented as ‘stress profiles’ showing the correlation between the applied force system and the stress distribution in the periodontal ligament, based on the improved finite element models, were established.     

Exposure to internal muscle tissue loads under the ischial tuberosities during sitting is elevated at abnormally high or low body mass indices

Deep tissue injury (DTI) is a severe pressure ulcer characteristic of chairfast or bedfast individuals, such as those with impaired mobility or neurological disorders. A DTI differs from superficial pressure ulcers in that the onset of DTI occurs under intact skin, in skeletal muscle tissue overlying bony prominences, and progression of the wound continues subcutaneously until skin breakdown. Due to the nature of this silently progressing wound, it is highly important to screen potentially susceptible individuals for their risk of developing a DTI. Abnormally low and high values of the body mass index (BMI) have been proposed to be associated with pressure ulcers, but a clear mechanism is lacking. We hypothesize that during sitting, exposure to internal muscle tissue loads under the ischial tuberosities (IT) is elevated at abnormally high or low body mass indices. Our aims in this study were: (a) to develop biomechanical models of the IT region in the buttocks that represent an individual who is gaining or losing weight drastically. (b) To determine changes in internal tissue load measures: principal compression strain, strain energy density (SED), principal compression stress and von Mises stress versus the BMI. (c) To determine percentage volumes of muscle tissue exposed to critical levels of the above load measures, which were defined based on our previous animal and tissue engineered model experiments: strain≥50%, stress≥2 kPa, SED≥0.5 kPa. A set of 21 finite element models, which represented the same individual, but with different BMI values within the normal range, above it and below it, was solved for the outcome measures listed above. The models had the same IT shape, size, distance between the IT, and (non-linear) mechanical properties for all soft tissues, but different thicknesses of gluteus muscles and fat tissue layers, corresponding to the BMI level. The resulted data indicated a trend of progressive increase in internal tissue loading, particularly in volumetric exposure to critical loading for BMI values outside the 17≤BMI≤22 kg/m² range, supporting our hypothesis for this study. We concluded that exposure to internal muscle tissue loads under the IT during sitting is optimally reduced at the low-normal BMI range, which is important not only in the context of DTI research, but also for understanding general sitting biomechanics.     

A computationally efficient strategy to estimate muscle forces in a finite element musculoskeletal model of the lower limb

Concurrent multiscale simulation strategies are required in computational biomechanics to study the interdependence between body scales. However, detailed finite element models rarely include muscle recruitment due to the computational burden of both the finite element method and the optimization strategies widely used to estimate muscle forces. The aim of this study was twofold: first, to develop a computationally efficient muscle force prediction strategy based on proportional-integral-derivative (PID) controllers to track gait and chair rise experimental joint motion with a finite element musculoskeletal model of the lower limb, including a deformable knee representation with 12 degrees of freedom; and, second, to demonstrate that the inclusion of joint-level deformability affects muscle force estimation by using two different knee models and comparing muscle forces between the two solutions. The PID control strategy tracked experimental hip, knee, and ankle flexion/extension with root mean square errors below 1°, and estimated muscle, contact and ligament forces in good agreement with previous results and electromyography signals. Differences up to 11% and 20% in the vasti and biceps femoris forces, respectively, were observed between the two knee models, which might be attributed to a combination of differing joint contact geometry, ligament behavior, joint kinematics, and muscle moment arms. The tracking strategy developed in this study addressed the inevitable tradeoff between computational cost and model detail in musculoskeletal simulations and can be used with finite element musculoskeletal models to efficiently estimate the interdependence between muscle forces and tissue deformation.