Verification, validation and sensitivity studies in computational biomechanics

Computational techniques and software for the analysis of problems in mechanics have naturally moved from their origins in the traditional engineering disciplines to the study of cell, tissue and organ biomechanics. Increasingly complex models have been developed to describe and predict the mechanical behavior of such biological systems. While the availability of advanced computational tools has led to exciting research advances in the field, the utility of these models is often the subject of criticism due to inadequate model verification and validation (V&V). The objective of this review is to present the concepts of verification, validation and sensitivity studies with regard to the construction, analysis and interpretation of models in computational biomechanics. Specific examples from the field are discussed. It is hoped that this review will serve as a guide to the use of V&V principles in the field of computational biomechanics, thereby improving the peer acceptance of studies that use computational modeling techniques.     

Spinal facet joint biomechanics and mechanotransduction in normal, injury and degenerative conditions

The facet joint is a crucial anatomic region of the spine owing to its biomechanical role in facilitating articulation of the vertebrae of the spinal column. It is a diarthrodial joint with opposing articular cartilage surfaces that provide a low friction environment and a ligamentous capsule that encloses the joint space. Together with the disc, the bilateral facet joints transfer loads and guide and constrain motions in the spine due to their geometry and mechanical function. Although a great deal of research has focused on defining the biomechanics of the spine and the form and function of the disc, the facet joint has only recently become the focus of experimental, computational and clinical studies. This mechanical behavior ensures the normal health and function of the spine during physiologic loading but can also lead to its dysfunction when the tissues of the facet joint are altered either by injury, degeneration or as a result of surgical modification of the spine. The anatomical, biomechanical and physiological characteristics of the facet joints in the cervical and lumbar spines have become the focus of increased attention recently with the advent of surgical procedures of the spine, such as disc repair and replacement, which may impact facet responses. Accordingly, this review summarizes the relevant anatomy and biomechanics of the facet joint and the individual tissues that comprise it. In order to better understand the physiological implications of tissue loading in all conditions, a review of mechanotransduction pathways in the cartilage, ligament and bone is also presented ranging from the tissue-level scale to cellular modifications. With this context, experimental studies are summarized as they relate to the most common modifications that alter the biomechanics and health of the spine-injury and degeneration. In addition, many computational and finite element models have been developed that enable more-detailed and specific investigations of the facet joint and its tissues than are provided by experimental approaches and also that expand their utility for the field of biomechanics. These are also reviewed to provide a more complete summary of the current knowledge of facet joint mechanics. Overall, the goal of this review is to present a comprehensive review of the breadth and depth of knowledge regarding the mechanical and adaptive responses of the facet joint and its tissues across a variety of relevant size scales.     

Discrete element analysis in musculoskeletal biomechanics

This paper is written to honor Professor Y. C. Fung, the applied mechanician who has made seminal contributions in biomechanics. His work has generated great spin-off utility in the field of musculoskeletal biomechanics. Following the concept of the Rigid Body-Spring Model theory by T. Kawai (1978) for non-linear analysis of beam, plate, and shell structures and the soil-gravel mixture foundation, we have derived a generalized Discrete Element Analysis (DEA) method to determine human articular joint contact pressure, constraining ligament tension and bone-implant interface stresses. The basic formulation of DEA to solve linear problems is reviewed. The derivation of non-linear springs for the cartilage in normal diarthrodial joint contact problem was briefly summarized. Numerical implementation of the DEA method for both linear and non-linear springs is presented. This method was able to generate comparable results to the classic contact stress problem (the Hertzian solution) and the use of Finite Element Modeling (FEM) technique on selected models. Selected applications in human knee and hip joints are demonstrated. In addition, the femoral joint prosthesis stem/bone interface stresses in a non-cemented fixation were analyzed using a 2D plane-strain approach. The DEA method has the advantages of ease in creating the model and reducing computational time for joints of irregular geometry. However, for the analysis of joint tissue stresses, the FEA technique remains the method of choice.     

Occupational spine biomechanics: A journey to the spinal frontier

This paper provides a brief introduction to the variety of research areas focusing on spine biomechanics as it pertains to understanding and preventing low back injuries in the workplace. While certainly not a comprehensive review of the literature, some of the earliest, pioneering studies are presented from the following areas: (1) spine tissue testing, (2) estimating spine tissue loading, (3) manual materials handling studies, (4) prolonged or repetitive spine loading, (5) ergonomic assessment tools, (6) sudden/unexpected loading and (7) spine stability. Where possible, some of our own research contributions are integrated into the relevant sections. This paper concludes with a suggestion of some future research directions to continue and enhance the important impact of occupational spine biomechanics.     

Mobility impairment in the elderly: challenges for biomechanics research.

The problems of mobility impairment in the elderly constitute new and major challenges for biomechanics research. This paper outlines what some of the important problems are, discusses the relevance of biomechanics research to these problems, and reviews some of the current state of knowledge about factors related to the biomechanics of mobility impairments in the elderly. The population of old adults is growing rapidly and the incidence of mobility impairments in old adults is high. Mobility impairment biomechanics research is needed to make the assessments of impairments more precise, to design therapeutic programs that are more effective and to learn more about how mobility impairments can be prevented.     

Interpreting principal components in biomechanics: Representative extremes and single component reconstruction

Principal component analysis is a powerful tool in biomechanics for reducing complex multivariate datasets to a subset of important parameters. However, interpreting the biomechanical meaning of these parameters can be a subjective process. Biomechanical interpretations that are based on visual inspection of extreme 5th and 95th percentile waveforms may be confounded when extreme waveforms express more than one biomechanical feature. This study compares interpretation of principal components using representative extremes with a recently developed method, called single component reconstruction, which provides an uncontaminated visualization of each individual biomechanical feature. Example datasets from knee joint moments, lateral gastrocnemius EMG, and lumbar spine kinematics are used to demonstrate that the representative extremes method and single component reconstruction can yield equivalent interpretations of principal components. However, single component reconstruction interpretation cannot be contaminated by other components, which may enhance the use and understanding of principal component analysis within the biomechanics community.     

Use of quadruped models in thoraco-abdominal biomechanics research

Pigs and dogs have become common models of human thoraco-abdominal impact response. This paper summarizes a comparative analysis of the dog and pig to the live human accomplished through a series of necropsies performed on pigs and dogs. The results are summarized below. Emphasis is placed on specific aspects which are felt to be important for impact biomechanics. In particular, emphasis is placed upon the effect of tethering structures because of their potential in explaining mechanisms of injury for specific types of trauma such as aortic and certain liver injuries. Some aspects of tethering in the pig and dog are significantly different from that of the live human so care should be taken when using these animals in thoraco-abdominal biomechanics experiments.     

Spine biomechanics

Current trends in spine research are reviewed in order to suggest future opportunities for biomechanics. Recent studies show that psychosocial factors influence back pain behaviour but are not important causes of pain itself. Severe back pain most often arises from intervertebral discs, apophyseal joints and sacroiliac joints, and physical disruption of these structures is strongly but variably linked to pain. Typical forms of structural disruption can be reproduced by severe mechanical loading in-vitro, with genetic and age-related weakening sometimes leading to injury under moderate loading. Biomechanics can be used to quantify spinal loading and movements, to analyse load distributions and injury mechanisms, and to develop therapeutic interventions. The authors suggest that techniques for quantifying spinal loading should be capable of measurement "in the field" so that they can be used in epidemiological surveys and ergonomic interventions. Great accuracy is not required for this task, because injury risk depends on tissue weakness as much as peak loading. Biomechanical tissue testing and finite-element modelling should complement each other, with experiments establishing proof of concept, and models supplying detail and optimising designs. Suggested priority areas for future research include: understanding interactions between intervertebral discs and adjacent vertebrae; developing prosthetic and tissue-engineered discs; and quantifying spinal function during rehabilitation. "Mechanobiology" has perhaps the greatest future potential, because spinal degeneration and healing are both mediated by the activity of cells which are acutely sensitive to their local mechanical environment. Precise characterisation and manipulation of this environment will be a major challenge for spine biomechanics.     

A dynamic model of jaw and hyoid biomechanics during chewing

Our understanding of human jaw biomechanics has been enhanced by computational modelling, but comparatively few studies have addressed the dynamics of chewing. Consequently, ambiguities remain regarding predicted jaw-gapes and forces on the mandibular condyles. Here, we used a new platform to simulate unilateral chewing. The model, based on a previous study, included curvilinear articular guidance, a mobile hyoid apparatus, and a compressible food bolus. Muscles were represented by Hill-type actuators with drive profiles tuned to produce target jaw and hyoid movements. The cycle duration was 732 ms. At maximum gape, the lower incisor-point was 20.1mm down, 5.8mm posterior, and 2.3mm lateral to its initial, tooth-contact position. Its maximum laterodeviation to the working-side during closing was 6.1mm, at which time the bolus was struck. The hyoid's movement, completed by the end of jaw-opening, was 3.4mm upward and 1.6mm forward. The mandibular condyles moved asymmetrically. Their compressive loads were low during opening, slightly higher on the working-side at bolus-collapse, and highest bilaterally when the teeth contacted. The model's movements and the directions of its condylar forces were consistent with experimental observations, resolving seeming discordances in previous simulations. Its inclusion of hyoid dynamics is a step towards modelling mastication.     

Epidermal differentiation governs engineered skin biomechanics

Engineered skin must be mechanically strong to facilitate surgical application and prevent damage during the early stages of engraftment. However, the evolution of structural properties during culture, the relative contributions of the epidermis and dermis, and any correlation with tissue morphogenesis are not well known. These aspects are investigated by assessing the mechanical properties of engineered skin (ES) and engineered dermis (ED) during a 21-day culture period, including correlations with cellular metabolism, cellular organization and epidermal differentiation. During culture, the epidermis differentiates and begins to cornify, as evidenced by immunostaining and surface electrical capacitance. Tensile testing reveals that the ultimate tensile strength and linear stiffness increase linearly with time for ES, but are relatively unchanged for ED. ES strength correlates significantly with epidermal differentiation (p<0.001) and a composite strength model indicates that strength is largely determined by the epidermis. These data suggest that strategies to improve ES biomechanics should target the dermis. Additionally, time-dependant changes in average ES strength and percent elongation can be used to set upper bound limits on mechanical stimulation profiles to avoid tissue damage.     

An upper bound computational model for investigation of fusion effects on adjacent segment biomechanics of the lumbar spine

Abstract

Prediction of the biomechanical effects of fusion surgery on adjacent segments is a challenge in computational biomechanics of the spine. In this study, a two-segment L3-L4-L5 computational model was developed to simulate the effects of spinal fusion on adjacent segment biomechanical responses under a follower load condition. The interaction between the degenerative segment (L4-5) and the adjacent segment (L3-4) was simulated using an equivalent follower spring. The spring stiffness was calibrated using a rigid fusion of a completely degenerated disc model at the L4-5 level, resulting in an upper bound response at the adjacent (L3-4) segment. The obtained upper bound equivalent follower spring was used to simulate the upper bound biomechanical responses of fusion of the disc with different degeneration grades. It was predicted that as the disc degeneration grade at the degenerative segment decreased, the effect on the adjacent segment responses decreased accordingly after fusion. The data indicated that the upper bound computational model can be a useful computational tool for evaluation of the interaction between segments and for investigation of the biomechanical mechanisms of adjacent segment degeneration after fusion.     

A case for spiking neural network simulation based on configurable multiple-FPGA systems

Recent neuropsychological research has begun to reveal that neurons encode information in the timing of spikes. Spiking neural network simulations are a flexible and powerful method for investigating the behaviour of neuronal systems. Simulation of the spiking neural networks in software is unable to rapidly generate output spikes in large-scale of neural network. An alternative approach, hardware implementation of such system, provides the possibility to generate independent spikes precisely and simultaneously output spike waves in real time, under the premise that spiking neural network can take full advantage of hardware inherent parallelism. We introduce a configurable FPGA-oriented hardware platform for spiking neural network simulation in this work. We aim to use this platform to combine the speed of dedicated hardware with the programmability of software so that it might allow neuroscientists to put together sophisticated computation experiments of their own model. A feed-forward hierarchy network is developed as a case study to describe the operation of biological neural systems (such as orientation selectivity of visual cortex) and computational models of such systems. This model demonstrates how a feed-forward neural network constructs the circuitry required for orientation selectivity and provides platform for reaching a deeper understanding of the primate visual system. In the future, larger scale models based on this framework can be used to replicate the actual architecture in visual cortex, leading to more detailed predictions and insights into visual perception phenomenon.     

Effect of loss of balance on biomechanics platform measures of sway: influence of stance and a method for adjustment

This paper describes a method for adjusting biomechanics platform measures of sway for loss of balance. Area and velocity measures of sway were determined in forty-seven elderly women, in double and single leg stance, first with their eyes open, then closed. Subjects were rarely able to complete 10 s trials during single leg stances. Therefore, a method was developed for eliminating data associated with loss of balance. Monitoring changes in vertical force and velocity by computer, those points exceeding trial specific thresholds associated with loss of balance were truncated. In double leg stances, loss of balance increased area measures by 0.3%, but did not effect velocity measures. In contrast, the loss of balance increased area measures by 0-3%, but did not effect velocity measures. In contrast, the loss of balance experienced by most subjects in single leg stance exaggerated area measures by 16-38%, and velocity measures by up to 10%. In double leg stances the correlations between unadjusted area measures and area measures adjusted for loss of balance ranged from 0.98 to 1.00. In single leg stances, the correlations for the area measures ranged from 0.69 to 0.89. The correlations between adjusted and unadjusted velocity measures were 1.00 and 0.93 for the double and single leg stances respectively. Although the question of which sway measure is best remains unanswered, this study provides useful data for future research. First, it demonstrates a method for modifying area representations of the center of pressure excursions for the effects of loss of balance.(ABSTRACT TRUNCATED AT 250 WORDS)     

Use of an overhead goal alters vertical jump performance and biomechanics

This study examined whether an extrinsic motivator, such as an overhead goal, during a plyometric jump may alter movement biomechanics. Our purpose was to examine the effects of an overhead goal on vertical jump height and lower-extremity biomechanics during a drop vertical jump and to compare the effects on female (N = 18) versus male (N = 17) athletes. Drop vertical jump was performed both with and without the use of an overhead goal. Greater vertical jump height (p = 0.002) and maximum takeoff external knee flexion (quadriceps) moment (p = 0.04) were attained with the overhead goal condition versus no overhead goal. Men had significantly greater vertical jump height (p < 0.001), maximum takeoff vertical force (p = 0.009), and maximum takeoff hip extensor moment (p = 0.02) compared with women. A significant gender x overhead goal interaction was found for stance time (p = 0.02) and maximum ankle (p = 0.04) and knee flexion angles (p = 0.04), with shorter stance times and lower angles in men during overhead goal time. These results indicate that overhead goals may be incorporated during training and testing protocols to alter lower-extremity biomechanics and can increase performance.     

Modeling the biomechanics of fetal movements

Fetal movements in the uterus are a natural part of development and are known to play an important role in normal musculoskeletal development. However, very little is known about the biomechanical stimuli that arise during movements in utero, despite these stimuli being crucial to normal bone and joint formation. Therefore, the objective of this study was to create a series of computational steps by which the forces generated during a kick in utero could be predicted from clinically observed fetal movements using novel cine-MRI data of three fetuses, aged 20–22 weeks. A custom tracking software was designed to characterize the movements of joints in utero, and average uterus deflection of \(6.95 \pm 0.41\) mm due to kicking was calculated. These observed displacements provided boundary conditions for a finite element model of the uterine environment, predicting an average reaction force of \(0.52 \pm 0.15\) N generated by a kick against the uterine wall. Finally, these data were applied as inputs for a musculoskeletal model of a fetal kick, resulting in predicted maximum forces in the muscles surrounding the hip joint of approximately 8 N, while higher maximum forces of approximately 21 N were predicted for the muscles surrounding the knee joint. This study provides a novel insight into the closed mechanical environment of the uterus, with an innovative method allowing elucidation of the biomechanical interaction of the developing fetus with its surroundings.     

The art and science of multi-scale citizen science support

Citizen science and community-based monitoring programs are increasing in number and breadth, generating volumes of scientific data. Many programs are ill-equipped to effectively manage these data. We examined the art and science of multi-scale citizen science support, focusing on issues of integration and flexibility that arise for data management when programs span multiple spatial, temporal, and social scales across many domains. Our objectives were to: (1) briefly review existing citizen science approaches and data management needs; (2) propose a framework for multi-scale citizen science support; (3) develop a cyber-infrastructure to support citizen science program needs; and (4) describe lessons learned. We find that approaches differ in scope, scale, and activities and that the proposed framework situates programs while guiding cyber-infrastructure system development. We built a cyber-infrastructure support system for citizen science programs (www.citsci.org) and show that carefully designed systems can be adept enough to support programs at multiple spatial and temporal scales across many domains when built with a flexible architecture. The advantage of a flexible, yet controlled, cyber-infrastructure system lies in the ability of users with different levels of permission to easily customize the features themselves, while adhering to controlled vocabularies necessary for cross-discipline comparisons and meta-analyses. Program evaluation tied to this framework and integrated into cyber-infrastructure support systems will improve our ability to track effectiveness. We compare existing systems and discuss the importance of standards for interoperability and the challenges associated with system maintenance and long-term support. We conclude by offering a vision of the future of citizen science data management and cyber-infrastructure support.     

Biomechanics of Musculoskeletal System and Its Biomimetic Implications： A Review

Biological musculoskeletal system （MSK）, composed of numerous bones, cartilages, skeletal muscles, tendons, ligaments etc., provides form, support, movement and stability for human or animal body. As the result of million years of selection and evolution, the biological MSK evolves to be a nearly perfect mechanical mechanism to support and transport the human or animal body, and would provide enormously rich resources to inspire engineers to innovate new technology and methodology to develop robots and mechanisms as effective and economical as the biological systems. This paper provides a general review of the current status of musculoskeletal biomechanics studies using both experimental and computational methods. This includes the use of the latest three-dimensional motion analysis systems, various medical imaging modalities, and also the advanced rigid-body and continuum mechanics musculoskeletal modelling techniques. Afterwards, several representative biomimetic studies based on ideas and concepts inspired from the structures and biomechanical functions of the biological MSK are dis- cussed. Finally, the major challenges and also the future research directions in musculoskeletal biomechanics and its biomimetic studies are proposed.