Finite element simulation of interactions between pelvic organs: Predictive model of the prostate motion in the context of radiotherapy

The setting up of predictive models of the pelvic organ motion and deformation may prove an efficient tool in the framework of prostate cancer radiotherapy, in order to deliver doses more accurately and efficiently to the clinical target volume (CTV). A finite element (FE) model of the prostate, rectum and bladder motion has been developed, investigating more specifically the influence of the rectum and bladder repletions on the gland motion. The required organ geometries are obtained after processing the computed tomography (CT) images, using specific softwares. Due to their structural characteristics, a 3D shell discretization is adopted for the rectum and the bladder, whereas a volume discretization is adopted for the prostate. As for the mechanical behavior modelling, first order Ogden hyperelastic constitutive laws for both the rectum and bladder are identified. The prostate is comparatively considered as more rigid and is accordingly modelled as an elastic tissue undergoing small strains. A FE model is then created, accounting for boundary and contact conditions, internal and applied loadings being selected as close as possible to available anatomic data.The order of magnitude of the prostate motion predicted by the FE simulations is similar to the measurements done on a deceased person, accounting for the delineation errors, with a relative error around 8%. Differences are essentially due to uncertainties in the constitutive parameters, pointing towards the need for the setting up of direct measurement of the organs mechanical behavior.     

髋骨术后仿真生物力学分析与临床研究

目的:探讨坐耻骨肿瘤切除术后髋骨生物力学变化及其与恢复期并发症产生的关系,指导术后骨盆重建.方法:构建单侧完整髋骨及坐耻骨肿瘤切除术后髋骨有限元模型,在相同约束和负载条件下,计算髋骨相同部位位移、应力及应变值,比较与分析模拟手术前后力学变化,结合临床评价意义.结果:手术前后髋骨节点位移、应力及应变分布区域相似,手术前后骶髂关节节点位移及应变、髋臼顶节点应力及应变有显著性差异;手术前后骶髂关节节点应力、坐骨大切迹应力及应变无显著性差异.结论:坐耻骨肿瘤切除术后主要影响同侧骶髂关节、髋臼顶的生物力学状态,恢复期并发症发生的原因由骶髂关节、髋臼顶生物力学改变及耻骨联合与耻骨上支连接与稳定作用的消失所致,三者相互作用,相互影响.     

Effect of bone fragment impact velocity on biomechanical parameters related to spinal cord injury: A finite element study

Several experimental and computational studies have investigated the effect of bone fragment impact on the spinal cord during trauma. However, the effect of the impact velocity of a fragment generated by a burst fracture on the stress and strain inside the spinal cord has not been computationally investigated, even though spinal canal occlusion and peak pressure at various impact velocities were provided in experimental studies. These stresses and strains are known factors related to clinical symptoms or injuries. In this study, a fluid-structure interaction model of the spinal cord, dura mater, and cerebrospinal fluid was developed and validated. The von-Mises stress distribution in the cord, the longitudinal strain, the cord compression and cross-sectional area at the impact center, and the obliteration of the cerebrospinal fluid layer were analyzed for three pellet sizes at impact velocities ranging from 1.5 m/s to 7.5 m/s. The results indicate that stress in the cord was substantially elevated when the initial impact velocity of the pellet exceeded a threshold of 4.5 m/s. Cord compression, reduction in cross-sectional area, and obliteration of the cerebrospinal fluid increased gradually as the velocity of the pellet increased, regardless of the size of the pellet. The present study provides insight into the mechanisms underlying spinal cord injury.     

Elevation and orientation of external loads influence trunk neuromuscular response and spinal forces despite identical moments at the L5–S1 level

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 the L5–S1 as well as forces in extensor thoracic muscles progressively decreased as orientation of external forces varied from downward gravity (90°) all the way to upward (−25°) orientation. In contrast, forces in local lumbar muscles followed reverse trends. Under larger horizontal forces at a lower elevation, lumbar muscles were much more active whereas extensor thoracic muscle forces were greater under smaller forces at a higher elevation. Despite such differences in activity pattern, the spinal forces remained nearly identical (<6% under 45 Nm moment). The published recorded surface EMG data of extensor muscles trend-wise agreed with computed local muscle forces as horizontal load elevation varied but were overall different from results in both local and global muscles when load orientation altered. Predictions demonstrate the marked effect of external force orientation and elevation on the trunk neuromuscular response and spinal forces and questions attempts to estimate spinal loads based only on consideration of moments at a spinal level.     

Simple finite element models for use in the design of therapeutic footwear

Therapeutic footwear is frequently prescribed in cases of rheumatoid arthritis and diabetes to relieve or redistribute high plantar pressures in the region of the metatarsal heads. Few guidelines exist as to how these interventions should be designed and what effect such interventions actually have on the plantar pressure distribution. Finite element analysis has the potential to assist in the design process by refining a given intervention or identifying an optimal intervention without having to actually build and test each condition. However, complete and detailed foot models based on medical image segmentation have proven time consuming to build and computationally expensive to solve, hindering their utility in practice. Therefore, the goal of the current work was to determine if a simplified patient-specific model could be used to assist in the design of foot orthoses to reduce the plantar pressure in the metatarsal head region. The approach is illustrated by a case study of a diabetic patient experiencing high pressures and pain over the fifth metatarsal head. The simple foot model was initially calibrated by adjusting the individual loads on the metatarsals to approximate measured peak plantar pressure distributions in the barefoot condition to within 3%. This loading was used in various shod conditions to identify an effective orthosis. Model results for metatarsal pads were considerably higher than measured values but predictions for uniform surfaces were generally within 16% of measured values. The approach enabled virtual prototyping of the orthoses, identifying the most favorable approach to redistribute the patient’s plantar pressures.     

A Visco-hyperelastic Model With Damage for the Knee Ligaments Under Dynamic Constraints

The aim of this study was to identify the behaviour laws governing the knee ligaments, accounting for the damage incurred by the structure under dynamic constraints. The model is developed using a thermodynamic formulation based on the coupling between a viscoelastic model and a damage model. Identification is carried out using the results of dynamic traction tests performed on a bone ligament/bone complex to which traction velocities of around 1.98 m/s were applied. The results show the ability of the model to account for the brittle and ductile failure processes occurring in the cruciate and lateral ligaments, respectively.     

Impact of material and morphological parameters on the mechanical response of the lumbar spine – A finite element sensitivity study

Finite element models are frequently used to study lumbar spinal biomechanics. Deterministic models are used to reflect a certain configuration, including the means of geometrical and material properties, while probabilistic models account for the inherent variability in the population. Because model parameters are generally uncertain, their predictive power is frequently questioned. In the present study, we determined the sensitivities of spinal forces and motions to material parameters of intervertebral discs, vertebrae, and ligaments and to lumbar morphology. We performed 1200 model simulations using a generic model of the human lumbar spine loaded under pure moments. Coefficients of determination and of variation were determined for all parameter and response combinations. Material properties of the vertebrae displayed the least impact on results, whereas those of the discs and morphology impacted most. The most affected results were the axial compression forces in the vertebral body and in several ligaments during flexion and the facet-joint forces during extension. Intervertebral rotations were considerably affected only when several parameters were varied simultaneously. Results can be used to decide which model parameters require careful consideration in deterministic models and which parameters might be omitted in probabilistic studies. Findings allow quantitative estimation of a model׳s precision.     

The shape and mobility of the thoracic spine in asymptomatic adults – A systematic review of in vivo studies

A comprehensive knowledge of the thoracic shape and kinematics is essential for effective risk prevention, diagnose and proper management of thoracic disorders and assessment of treatment or rehabilitation strategies as well as for in silico and in vitro models for realistic applications of boundary conditions.After an extensive search of the existing literature, this study summarizes 45 studies on in vivo thoracic kyphosis and kinematics and creates a systematic and detailed database. The thoracic kyphosis over T1–12 determined using non-radiological devices (34°) was relatively less than measured using radiological devices (40°) during standing. The majority of kinematical measurements are based on non-radiological devices. The thoracic range of motion (RoM) was greatest during axial rotation (40°), followed by lateral bending (26°), and flexion (21°) when determined using non-radiological devices during standing. The smallest RoM was identified during extension (13°). The lower thoracic level (T8–12) contributed more to the RoM than the upper (T1–4) and middle (T4–8) levels during flexion and lateral bending. During axial rotation and extension, the middle level (T4–8) contributed the most. Coupled motion was evident, mostly during lateral bending and axial rotation. With aging, the thoracic kyphosis increased by about 3° per decade, whereas the RoM decreased by about 5° per decade for all load directions. These changes with aging mainly occurred in the lower region (T6–12). The influence of sex on thoracic kyphosis and the RoM has been described as partly contradictory. Obesity was found to decrease the thoracic RoM. Studies comparing standing, sitting and lying reported the effect of posture as significant.     

Comparison of eight published static finite element models of the intact lumbar spine: Predictive power of models improves when combined together

Finite element (FE) model studies have made important contributions to our understanding of functional biomechanics of the lumbar spine. However, if a model is used to answer clinical and biomechanical questions over a certain population, their inherently large inter-subject variability has to be considered. Current FE model studies, however, generally account only for a single distinct spinal geometry with one set of material properties. This raises questions concerning their predictive power, their range of results and on their agreement with in vitro and in vivo values.