Sharif Digital Repository

Accurate estimation of muscle forces in various occupational tasks is critical for a reliable evaluation of spinal loads and subsequent assessment of risk of injury and management of back disorders. The majority of biomechanical models of multi-segmental spine estimate muscle forces and spinal loads based on the balance of net moments at a single level with no consideration for the equilibrium at remaining levels. This work aimed to quantify the extent of equilibrium violation and alterations in estimations when such models are performed at different levels. Results are compared with those of kinematics-driven model that satisfies equilibrium at all levels and EMG data. Regardless of the... 

Spinal load models have become an increasingly valuable tool for the evaluation of the stress placed on the spine. In order to get an accurate representation of these spinal loads, they must be compared to known tolerance values. Bone mineral content and density of the lumbar spine of 23 males and 21 females was measured using a dual energy x-ray absorptiometry. Compression tolerance values were predicted by previously published studies based upon bone mineral levels. Anthropometric measurements were recorded and related to the compression tolerance values through the use of multivariate linear regression techniques. Compression tolerance values based on the bone mineral content or density... 

Based on the force-deflection equation for a beam subjected to lateral point loads, a C2 continuous piecewise bicubic mathematical representation was proposed to model complicated geometrical surfaces, e.g. the articular surfaces of human joints. The method was then extended so that it could be used for mathematical modelling of incomplete nets of data points, as well as smoothing of noisy and/or filtering of erroneous data points. Mathematical techniques were also developed to calculate the required unknown parameters explicitly, with no need to solve the system of equations simultaneously. The performance of the proposed method was evaluated on a number of surface modelling problems,... 

This work aimed to evaluate trunk muscle forces, internal loads and stability margin under some simulated standing postures, with and without external loads, using a nonlinear finite element model of the T1-S1 spine with realistic nonlinear load-displacement properties. A novel kinematics-based algorithm was applied that exploited a set of spinal sagittal rotations, initially calculated to minimize balancing moments, to solve the redundant active-passive system. The loads consisted of upper body gravity distributed along the spine with or without 200 N held in the hands, either in the front of the body or on the sides. Nonlinear and linear stability/perturbation analyses at deformed,... 

A computational framework was developed to simulate the bone remodelling process as a structural topology optimisation problem. The mathematical formulation of the Level Set technique was extended and then implemented into a coronal plane model of the proximal femur to simulate the remodelling of internal structure and external geometry of bone into the optimal state. Results indicated that the proposed approach could reasonably mimic the major geometrical and material features of the natural bone. Simulation of the internal bone remodelling on the typical gross shape of the proximal femur, resulted in a density distribution pattern with good consistency with that of the natural bone. When... 

Traditional load-control musculoskeletal and finite element (FE) models of the spine fail to accurately predict in vivo intervertebral joint loads due mainly to the simplifications and assumptions when estimating redundant trunk muscle forces. An alternative powerful protocol that bypasses the calculation of muscle forces is to drive the detailed FE models by image-based in vivo displacements. Development of subject-specific models, however, both involves the risk of extensive radiation exposures while imaging in supine and upright postures and is time consuming in terms of the reconstruction of the vertebrae, discs, ligaments, and facets geometries. This study therefore aimed to introduce a... 

Spine musculoskeletal (MS) models make simplifying assumptions on the intervertebral joint degrees-of-freedom (rotational and/or translational), representation (spherical or beam-like joints), and properties (linear or nonlinear). They also generally neglect the realistic structure of the joints with disc nuclei/annuli, facets, and ligaments. We aim to develop a novel MS model where trunk muscles are incorporated into a detailed finite element (FE) model of the ligamentous T12-S1 spine thus constructing a gold standard coupled MS-FE model. Model predictions are compared under some tasks with those of our earlier spherical joints, beam joints, and hybrid (uncoupled) MS-FE models. The coupled... 

Epidemiological studies are divided over the causative role of body weight (BW) in low back pain. Biomechanical modeling is a valuable approach to examine the effect of changes in BW on spinal loads and risk of back pain. Changes in BW have not been properly simulated by previous models as associated alterations in model inputs on the musculature and moment arm of gravity loads have been neglected. A detailed, multi-joint, scalable model of the thoracolumbar spine is used to study the effect of BW (varying at five levels, i.e., 51, 68, 85, 102, and 119kg) on the L5-S1 spinal loads during various static symmetric activities while scaling moment arms and physiological cross-sectional areas of... 

Recent advances in medical imaging techniques have allowed pure displacement-control trunk models to estimate spinal loads with no need to calculate muscle forces. Sensitivity of these models to the errors in post-imaging evaluation of displacements (reported to be ∼0.4–0.9° and 0.2–0.3 mm in vertebral displacements) has not yet been investigated. A Monte Carlo analysis was therefore used to assess the sensitivity of results in both musculoskeletal (MS) and passive finite element (FE) spine models to errors in measured displacements. Six static activities in upright standing, flexed, and extended postures were initially simulated using a force-control hybrid MS-FE model. Computed vertebral... 

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... 

A wide range of loading conditions involving external forces with varying magnitudes, orientations and locations are encountered in daily activities. Here we computed the effect on trunk biomechanics of changes in force location (two levels) and orientation (5 values) in 4 subjects in upright standing while maintaining identical external moment of 15. Nm, 30. N. m or 45. Nm at the L5-S1. Driven by measured kinematics and gravity/external loads, the finite element models yielded substantially different trunk neuromuscular response with moderate alterations (up to 24% under 45 Nm moment) in spinal loads as the load orientation varied. Under identical moments, compression and shear forces at... 

A safe rehabilitation exercise for anterior cruciate ligament (ACL) injuries needs to be compatible with the normal knee arthrokinematics to avoid abnormal loading on the joint structures. The objective of this study was to measure the amount of the anterior tibial translation (ATT) of the ACL-deficient knees during selective open and closed kinetic chain exercises. The intact and injured knees of fourteen male subjects with unilateral ACL injury were imaged using uniplanar fluoroscopy, while the subjects performed forward lunge and unloaded/loaded open kinetic knee extension exercises. The ATTs were measured from fluoroscopic images, as the distance between the tibial and femoral reference... 

In this study, a nonlinear poroelastic model of intervertebral disc as an infrastructure was developed. Moreover, a new element was defined consisting a disc (Viscoelastic Euler Beam Element) and a vertebra (Rigid Link) as a unit element. Using the new element, three different viscoelastic finite element models were prepared for lumbar motion segment (L4/L5). Prolonged loading (short-term and long-term creep) and cyclic loading were applied to the models and the results were compared with results of in vivo tests. Simplification of the models by using the new element leads to reduction of the runtime of the models in dynamic analyses to few minutes without losing the accuracy in the results  

A new wearable assistive device (WAD) was developed to decrease required force on the lumbar spine in static holding tasks. In order to obtain moments on lumbar spine in two conditions, with and without WAD, a biomechanical static model was used for estimation of external moments on lumbar spine. The results of biomechanical models indicated that there was a reduction in the lumbar moment ranging from 20% to 43% using WAD depending on the load and flexion angle. A total of 15 male healthy subjects were tested to experimentally verify the predicted reduction of external moments on the spine by wearing WAD. Normalized electromyography (EMG) of the right and left lumbar and thoracic erector... 

This study was designed to investigate the effect of changing the location of the center of mass (COM) of transfemoral prostheses on the spatiotemporal and kinematical characteristics of the amputee gait, while maintaining the prosthetic mass fixed. Ten men with unilateral traumatic transfemoral amputation participated in gait analysis, 2-min walk and subject preference tests. Weights were added to the original prosthetic legs in three conditions: 600. g added to the ankle, 600. g added at 10. cm distal to the prosthetic knee, and 300. g added to the ankle, and 300. g at 10. cm distal to the knee. For each prosthetic mass condition, the stride and step lengths, stepping speed, stance, swing... 

Multiple joint interactions are critical to produce stable coordinated movements and can be influenced by low back pain and task conditions. Inter-segmental coordination pattern and variability were assessed in subjects with and without chronic nonspecific low back pain (CNSLBP). Kinematic data were collected from 22 CNSLBP and 22 healthy volunteers during repeated trunk flexion-extension in various conditions of symmetry, velocity, and loading; each at two levels. Sagittal plane angular data were time normalized and used to calculate continuous relative phase for each data point. Mean absolute relative phase (MARP) and deviation phase (DP) were derived to quantify lumbar-pelvis and... 

Atrogenic injuries to paraspinal muscles during the posterior lumbar surgery (PLS) cause a reduction in their cross-sectional areas (CSAs) and contractile densities over time post-surgery. This study aims to quantify such alterations. Method: Pre- and postoperative CSAs (~6 months interval) of all paraspinal muscles were measured in six patients undergoing PLS using a 3-T magnetic resonance (MR) scanner to quantify the alterations in geometrical and tissue effective contractile (non-fatty) CSAs of these muscles at all lumbar levels. To examine the presence of any confounding effects on recorded changes within ~7-month period, measurements were also carried out on ten healthy volunteers.... 

Underlying mechanisms of obesity-related back pain remain unexplored. Thus, we aim to determine the effect of obesity and its shapes on the spinal loads and the associated risks of injury. Obesity shapes were initially constructed by principal component analysis based on datasets on 5852 obese individuals. Spinal loads, cycles to vertebral failure and trunk stability margin were estimated in a subject-specific trunk model taking account of personalized musculature, passive ligamentous spine, obesity shapes, segmental weights, spine kinematics and bone mineral density. Three obesity shapes (mean and extreme abdominal circumferences) at three body weights (BWs) of 86, 98 and 109 kg were... 

During pregnancy, traumas can threaten maternal and fetal health. Various trauma effects on a pregnant uterus are little investigated. In the present study, a finite element model of a uterus along with a fetus, placenta, amniotic fluid, and two most effective ligament sets is developed. This model allows numerical evaluation of various loading on a pregnant uterus. The model geometry is developed based on CT-scan data and validated using anthropometric data. Applying Ogden hyper-elastic theory, material properties of uterine wall and placenta are developed. After simulating the “rigid-bar” abdominal loading, the impact force and abdominal penetration are investigated. Findings are compared... 

Finite element analysis is an effective tool to evaluate the material properties of living tissue. For an interactive optimization procedure, the finite element analysis usually needs many simulations to reach a reasonable solution. The metamodel analysis of finite element simulation can be used to reduce the computation of a structure with complex geometry or a material with composite constitutive equations. The intervertebral disc is a complex, heterogeneous, and hydrated porous structure. A poroelastic finite element model can be used to observe the fluid transferring, pressure deviation, and other properties within the disc. Defining reasonable poroelastic material properties of the...