Marwan El-Rich

Finite element investigation of the loading rate effect on the spinal load-sharing changes under impact conditions

Sudden deceleration and frontal/rear impact configurations involve rapid movements that can cause spinal injuries. This study aimed to investigate the rotation rate effect on the L2-L3 motion segment load-sharing and to identify which spinal structure is at risk of failure and at what rotation velocity the failure may initiate? Five degrees of sagittal rotations at different rates were applied in a detailed finite-element model to analyze the responses of the soft tissues and the bony structures until possible fractures. The structural response was markedly different under the highest velocity that caused high peaks of stresses in the segment compared to the intermediate and low velocities. Under flexion, the stress was concentrated at the upper pedicle region of L2 and fractures were firstly initiated in this region and then in the lower endplate of L2. Under extension, maximum stress was located in the lower pedicle region of L2 and fractures started in the left facet joint, then they expanded in the lower endplate and in the pedicle region of L2. No rupture has resulted at the lower or intermediate velocities. The intradiscal pressure was higher under flexion and decreased when the endplate was fractured, while the contact forces were greater under extension and decreased when the facet surface was cracked. The highest ligaments stresses were obtained under flexion and did not reach the rupture values. The endplate, pedicle and facet surface represented the potential sites of bone fracture. Results showed that spinal injuries can result at sagittal rotation velocity exceeding 0.5 degrees /ms.

Effect of load position on muscle forces, internal loads and stability of the human spine in upright postures

A novel kinematics-based approach coupled with a non-linear finite element model was used to investigate the effect of changes in the load position and posture on muscle activity, internal loads and stability margin of the human spine in upright standing postures. In addition to 397 N gravity, external loads of 195 and 380 N were considered at different lever arms and heights. Muscle forces, internal loads and stability margin substantially increased as loads displaced anteriorly away from the body. Under same load magnitude and location, adopting a kyphotic posture as compared with a lordotic one increased muscle forces, internal loads and stability margin. An increase in the height of a load held at a fixed lever arm substantially diminished system stability thus requiring additional muscle activations to maintain the same margin of stability. Results suggest the importance of the load position and lumbar posture in spinal biomechanics during various manual material handling operations.

Calibration of Hyperelastic Material Properties of the Human Lumbar Intervertebral Disc under Fast Dynamic Compressive Loads

Under fast dynamic loading conditions (e.g. high-energy impact), the load rate dependency of the intervertebral disc (IVD) material properties may play a crucial role in the biomechanics of spinal trauma. However, most finite element models (FEM) of dynamic spinal trauma uses material properties derived from quasi-static experiments, thus neglecting this load rate dependency. The aim of this study was to identify hyperelastic material properties that ensure a more biofidelic simulation of the IVD under a fast dynamic compressive load. A hyperelastic material law based on a first-order Mooney-Rivlin formulation was implemented in a detailed FEM of a L2-L3 functional spinal unit (FSU) to represent the mechanical behavior of the IVD. Bony structures were modeled using an elasto-plastic Johnson-Cook material law that simulates bone fracture while ligaments were governed by a viscoelastic material law. To mimic experimental studies performed in fast dynamic compression, a compressive loading velocity of 1 m/s was applied to the superior half of L2, while the inferior half of L3 was fixed. An exploratory technique was used to simulate dynamic compression of the FSU using 34 sets of hyperelastic material constants randomly selected using an optimal Latin hypercube algorithm and a set of material constants derived from quasi-static experiments. Selection or rejection of the sets of material constants was based on compressive stiffness and failure parameters criteria measured experimentally. The two simulations performed with calibrated hyperelastic constants resulted in nonlinear load-displacement curves with compressive stiffness (7335 and 7079 N/mm), load (12,488 and 12,473 N), displacement (1.95 and 2.09 mm) and energy at failure (13.5 and 14.7 J) in agreement with experimental results (6551 ± 2017 N/mm, 12,411 ± 829 N, 2.1 ± 0.2 mm and 13.0 ± 1.5 J respectively). The fracture pattern and location also agreed with experimental results. The simulation performed with constants derived from quasi-static experiments showed a failure energy (13.2 J) and a fracture pattern and location in agreement with experimental results, but a compressive stiffness (1580 N/mm), a failure load (5976 N) and a displacement to failure (4.8 mm) outside the experimental corridors. The proposed method offers an innovative way to calibrate the hyperelastic material properties of the IVD and to offer a more realistic simulation of the FSU in fast dynamic compression.

Investigation of impact loading rate effects on the ligamentous cervical spinal load-partitioning using finite element model of functional spinal unit C2-C3.

The cervical spine functions as a complex mechanism that responds to sudden loading in a unique manner, due to intricate structural features and kinematics. The spinal load-sharing under pure compression and sagittal flexion/extension at two different impact rates were compared using a bio-fidelic finite element (FE) model of the ligamentous cervical functional spinal unit (FSU) C2-C3. This model was developed using a comprehensive and realistic geometry of spinal components and material laws that include strain rate dependency, bone fracture, and ligament failure. The range of motion, contact pressure in facet joints, failure forces in ligaments were compared to experimental findings. The model demonstrated that resistance of spinal components to impact load is dependent on loading rate and direction. For the loads applied, stress increased with loading rate in all spinal components, and was concentrated in the outer intervertebral disc (IVD), regions of ligaments to bone attachment, and in the cancellous bone of the facet joints. The highest stress in ligaments was found in capsular ligament (CL) in all cases. Intradiscal pressure (IDP) in the nucleus was affected by loading rate change. It increased under compression/flexion but decreased under extension. Contact pressure in the facet joints showed less variation under compression, but increased significantly under flexion/extension particularly under extension. Cancellous bone of the facet joints region was the only component fractured and fracture occurred under extension at both rates. The cervical ligaments were the primary load-bearing component followed by the IVD, endplates and cancellous bone; however, the latter was the most vulnerable to extension as it fractured at low energy impact.

On the load-sharing along the ligamentous lumbosacral spine in flexed and extended postures: Finite element study.

A harmonic synergy between the load-bearing and stabilizing components of the spine is necessary to maintain its normal function. This study aimed to investigate the load-sharing along the ligamentous lumbosacral spine under sagittal loading. A 3D nonlinear detailed Finite Element (FE) model of lumbosacral spine with realistic geometry was developed and validated using wide range of numerical and experimental (in-vivo and in-vitro) data. The model was subjected to 500 N compressive Follower Load (FL) combined with 7.5 Nm flexion (FLX) or extension (EXT) moments. Load-sharing was expressed as percentage of total internal force/moment developed along the spine that each spinal component carried. These internal forces and moments were determined at the discs centres and included the applied load and the resisting forces in the ligaments and facet joints. The contribution of the facet joints and ligaments in supporting bending moments produced additional forces and moments in the discs. The intervertebral discs carried up to 81% and 68% of the total internal force in case of FL combined with FLX and EXT, respectively. The ligaments withstood up to 67% and 81% of the total internal moment in cases of FL combined with EXT and FLX, respectively. Contribution of the facet joints in resisting internal force and moment was noticeable at levels L4-S1 only particularly in case of FL combined with EXT and reached up 29% and 52% of the internal moment and force, respectively. This study demonstrated that spinal load-sharing depended on applied load and varied along the spine.

Ergonomic Analysis and the Need for Its Integration for Planning and Assessing Construction Tasks

AbstractThe execution of daily construction tasks exposes workers to one or multiple ergonomic risk factors (awkward postures, static force, vibration, repetition, environmental risk, contact stress) and thus varying risks of developing musculoskeletal disorders. As a result, musculoskeletal disorders are common issues in construction and result in costly delays and disability claims. Though there is recent research investigating the epidemiology and causal factors for musculoskeletal injury, the construction industry has not fully embraced this as part of its safety practices. This study presents state-of-the-art ergonomic techniques, Canadian ergonomic legislation, and work-related musculoskeletal disorder (WRMSD) lost-time claims (LTC) statistics to show the resultant economic (cash and productivity) losses and adverse social (occupational health and safety) impact of WRMSDs resulting from current practice and legislation. The potential short- and long-term productivity and cost merits of incorporating...

Symmetry analysis of talus bone: A Geometric morphometric approach

Objective The main object of this study was to use a geometric morphometric approach to quantify the left-right symmetry of talus bones. Methods Analysis was carried out using CT scan images of 11 pairs of intact tali. Two important geometric parameters, volume and surface area, were quantified for left and right talus bones. The geometric shape variations between the right and left talus bones were also measured using deviation analysis. Furthermore, location of asymmetry in the geometric shapes were identified. Results Numerical results showed that talus bones are bilaterally symmetrical in nature, and the difference between the surface area of the left and right talus bones was less than 7.5%. Similarly, the difference in the volume of both bones was less than 7.5%. Results of the three-dimensional (3D) deviation analyses demonstrated the mean deviation between left and right talus bones were in the range of -0.74 mm to 0.62 mm. It was observed that in eight of 11 subjects, the deviation in symmetry occurred in regions that are clinically less important during talus surgery. Conclusions We conclude that left and right talus bones of intact human ankle joints show a strong degree of symmetry. The results of this study may have significance with respect to talus surgery, and in investigating traumatic talus injury where the geometric shape of the contralateral talus can be used as control. Cite this article: Bone Joint Res 2014;3:139–45.

Effects of inter-individual lumbar spine geometry variation on load-sharing: Geometrically personalized Finite Element study

There is a large, at times contradictory, body of research relating spinal curvature to Low Back Pain (LBP). Mechanical load is considered as important factor in LBP etiology. Geometry of the spinal structures and sagittal curvature of the lumbar spine govern its mechanical behavior. Thus, understanding how inter-individual geometry particularly sagittal curvature variation affects the spinal load-sharing becomes of high importance in LBP assessment. This study calculated and compared kinematics and load-sharing in three ligamentous lumbosacral spines: one hypo-lordotic (Hypo-L) with low lordosis, one normal-lordotic (Norm-L) with normal lordosis, and one hyper-lordotic (Hyper-L) with high lordosis in flexed and extended postures using 3D nonlinear Finite Element (FE) modeling. These postures were simulated by applying Follower Load (FL) combined with flexion or extension moment. The Hypo-L spine demonstrated stiffer behavior in flexion but more flexible response to extension compared to the Norm-L spine. The excessive lordosis stiffened response of the Hyper-L spine to extension but did not affect its resistance to flexion compared to the Norm-L spine. Despite the different resisting actions of the posterior ligaments to flexion moment, the increase of disc compression was similar in all the spines leading to similar load-sharing. However, resistance of the facet joints to extension was more important in the Norm- and Hyper-L spines which reduced the disc compression. The spinal curvature strongly influenced the magnitude and location of load on the spinal components and also altered the kinematics and load-sharing particularly in extension. Consideration of the subject-specific geometry and sagittal curvature should be an integral part of mechanical analysis of the lumbar spine.

The influence of cell micro-structure on the in-plane dynamic crushing of honeycombs with negative Poisson’s ratio

The in-plane dynamic crushing behaviors and energy-absorbed characteristics of honeycombs with negative Poisson’s ratio (NPR) have been studied by means of explicit dynamic finite element analysis (DFEA) using ANSYS/LS-DYNA. First, the honeycomb models filled with different reentrant cells by the variation of micro-cell configuration parameters (cell-wall angle and shape ratio) are established. The respective influences of micro-structure and impact velocities on the deformation behaviors, the dynamic plateau stresses and the absorbed energy of reentrant honeycombs are explored in detail. It is shown that owing to the variation of cell micro-structure, reentrant honeycombs display different macro-/micro- deformation properties during the crushing. For the given impact velocity, the dynamic plateau stresses are related to the shape ratio by a power law and to the cell-wall angle by least-square curves. And they are also proportional to the square of impact velocities for a high impact velocity. Based on the finite element simulated results and one-dimensional shock wave theory, an empirical formula for auxetic honeycomb to predict the dynamic plateau stress is derived in terms of relative density and impact velocity.

Monitoring for idiopathic scoliosis curve progression using surface topography asymmetry analysis of the torso in adolescents

BACKGROUND CONTEXT At first visit and each clinical follow-up session, patients with adolescent idiopathic scoliosis (AIS) undergo radiographic examination, from which the Cobb angle is measured. The cumulative exposure to X-ray radiation justifies efforts in developing noninvasive methods for scoliosis monitoring. PURPOSE To determine the capability of the three-dimensional markerless surface topography (ST) asymmetry analysis to detect ≥5° progression in the spinal curvature in patients with AIS over 1-year follow-up interval. STUDY DESIGN/SETTING Cross-sectional study in a specialized scoliosis clinic. PATIENT SAMPLE In this study, baseline and 1-year follow-up full torso ST scans of 100 patients with AIS were analyzed using three-dimensional markerless asymmetry analysis. OUTCOME MEASURES Patients with ΔCobb≥5° and ΔCobb<5° were categorized into progression and nonprogression groups, respectively. METHODS The ST scan of each full torso was analyzed to calculate the best plane of symmetry by minimizing the distances between the torso and its reflection about the plane of symmetry. Distance between the torso and its reflection was measured and displayed as deviation color maps. The difference of ST measurements between two successive acquisitions was used to determine if the scoliosis has progressed at least 5° or not. The classification tree technique was implemented using the local deformity of the torso in the thoracic-thoracolumbar (T-TL) and lumbar (L) regions to categorize curves into progression and nonprogression groups. The change in maximum deviation and root mean square of the deviations in the torso were the parameters effective in capturing the curve progression. Funding for this research is provided by the Scoliosis Research Society, and Women and Childrens Health Research Institute. RESULTS The classification model detected 85.7% of the progression and 71.6% of the nonprogression cases. The resulting false-negative rate of 4% for T-TL curves, representing the proportion of undetected progressions, confirmed that the technique shows promise to monitor the progression of T-TL scoliosis curves. Although 100% L curves with progression were detected using the deviation color maps of the torsos, because of the small number of analyzed L curves, further research is needed before the efficiency of the method in capturing the L curves with progression is confirmed. CONCLUSIONS Using the developed classification tree for the patients analyzed in this study, 43% of nonprogression cases between two visits would not have to undergo an X-ray examination.