Modeling and control of anterior-posterior and medial-lateral sways in standing posture.

To study essential anterior-posterior and medial-lateral sways of the stance caused by rotational movements about the ankle and hip joints, a mathematical model is developed for the 3D postural kinematics and dynamics. The model is in the form of nonlinear differential-algebraic equations corresponding to a biomechanical system with holonomic constraints. A nonlinear feedback control law is further derived for stabilizing the upright stance, whilst eliminating internal torques induced by the constraints on postural movements. Numerical simulations of the model parametrized with experimental data of human body segments illustrate the performance of postural balancing with the proposed control.

New insights into the biomechanics of Legg-Calvé-Perthes' disease: The Role of Epiphyseal Skeletal Immaturity in Vascular Obstruction.

OBJECTIVES: Legg-Calvé-Perthes' disease (LCP) is an idiopathic osteonecrosis of the femoral head that is most common in children between four and eight years old. The factors that lead to the onset of LCP are still unclear; however, it is believed that interruption of the blood supply to the developing epiphysis is an important factor in the development of the condition. METHODS: Finite element analysis modelling of the blood supply to the juvenile epiphysis was investigated to understand under which circumstances the blood vessels supplying the femoral epiphysis could become obstructed. The identification of these conditions is likely to be important in understanding the biomechanics of LCP. RESULTS: The results support the hypothesis that vascular obstruction to the epiphysis may arise when there is delayed ossification and when articular cartilage has reduced stiffness under compression. CONCLUSION: The findings support the theory of vascular occlusion as being important in the pathophysiology of Perthes disease.Cite this article: M. Pinheiro, C. A. Dobson, D. Perry, M. J. Fagan. New insights into the biomechanics of Legg-Calvé-Perthes' disease: The Role of Epiphyseal Skeletal Immaturity in Vascular Obstruction. Bone Joint Res 2018;7:148-156. DOI: 10.1302/2046-3758.72.BJR-2017-0191.R1.

Validation of a morphometric reconstruction technique applied to a juvenile pelvis.

Three-dimensional reconstructions of bone geometry from microCT (computed tomography) data are frequently used in biomechanical and finite element analyses. Digitization of bone models is usually a simple process for specimens with a complete geometry, but in instances of damage or disarticulation it can be very challenging. Subsequent to digitization, further imaging techniques are often required to estimate the geometry of missing bone or connecting cartilage. This paper presents an innovative approach to the reconstruction of incomplete scan data, to reproduce proper anatomical arrangements of bones, including absent connecting cartilaginous elements. Utilizing geometric morphometric tools, the reconstruction technique is validated through comparison of a reconstructed 9 year old pelvis, to the original CT data. A principal component analysis and an overlay of the two pelves provide a measure of the accuracy of the reconstructed model. Future work aims to investigate the biomechanical effects of any minor positional error on the bone's predicted structural properties through the use of finite element analysis.

Masticatory loading and bone adaptation in the supraorbital torus of developing macaques.

Research on the evolution and adaptive significance of primate craniofacial morphologies has focused on adult, fully developed individuals. Here, we investigate the possible relationship between the local stress environment arising from masticatory loadings and the emergence of the supraorbital torus in the developing face of the crab-eating macaque Macaca fascicularis. By using finite element analysis (FEA), we are able to evaluate the hypothesis that strain energy density (SED) magnitudes are high in subadult individuals with resulting bone growth in the supraorbital torus. We developed three micro-CT-based FEA models of M. fascicularis skulls ranging in dental age from deciduous to permanent dentitions and validated them against published experimental data. Applied masticatory muscle forces were estimated from physiological cross-sectional areas of macaque cadaveric specimens. The models were sequentially constrained at each working side tooth to simulate the variation of the bite point applied during masticatory function. Custom FEA software was used to solve the voxel-based models and SED and principal strains were computed. A physiological superposition SED map throughout the face was created by allocating to each element the maximum SED value from each of the load cases. SED values were found to be low in the supraorbital torus region throughout ontogeny, while they were consistently high in the zygomatic arch and infraorbital region. Thus, if the supraorbital torus arises to resist masticatory loads, it is either already adapted in each of our subadult models so that we do not observe high SED or a lower site-specific bone deposition threshold must apply.

Assessing mechanical function of the zygomatic region in macaques: validation and sensitivity testing of finite element models.

Crucial to the interpretation of the results of any finite element analysis of a skeletal system is a test of the validity of the results and an assessment of the sensitivity of the model parameters. We have therefore developed finite element models of two crania of Macaca fascicularis and investigated their sensitivity to variations in bone material properties, the zygomatico-temporal suture and the loading regimen applied to the zygomatic arch. Maximum principal strains were validated against data derived from ex vivo strain gauge experiments using non-physiological loads applied to the macaque zygomatic arch. Elastic properties of the zygomatic arch bone and the zygomatico-temporal suture obtained by nanoindentation resulted in a high degree of congruence between experimental and simulated strains. The findings also indicated that the presence of a zygomatico-temporal suture in the model produced strains more similar to experimental values than a completely separated or fused arch. Strains were distinctly higher when the load was applied through the modelled superficial masseter compared with loading an array of nodes on the arch. This study demonstrates the importance of the accurate selection of the material properties involved in predicting strains in a finite element model. Furthermore, our findings strongly highlight the influence of the presence of craniofacial sutures on strains experienced in the face. This has implications when investigating craniofacial growth and masticatory function but should generally be taken into account in functional analyses of the craniofacial system of both extant and extinct species.

Three dimensional stereolithography models of cancellous bone structures from muCT data: testing and validation of finite element results.

Stereolithography (STL) models of complex cancellous bone structures have been produced from three-dimensional micro-computed tomography data sets of human cancellous bone histological samples from four skeletal sites. The STL models have been mechanically tested and the derived stiffness compared with that predicted by finite element analysis. The results show a strong correlation (R2 = 0.941) between the predicted and calculated stiffnesses of the structures and show promise for the use of STL as an additional technique to complement the use of finite element models, for the assessment of the mechanical properties of complex cancellous bone structures.

Stereo visualization of 3D trabecular bone structures produced by bone remodelling simulation.

Adult human bone is constantly being renewed by a process known as remodelling. For cancellous bone this renewal process occurs at the interface between bone and marrow where bone is depleted by osteoclasts and rebuilt by osteoblasts. This remodelling process allows bone to repair itself. Software simulators for bone remodelling provide insight into the bone remodelling process; they allow investigation into bone form and structural properties, and they also allow the emulation of bone diseases and possible treatments for these diseases over long periods of time. BONESIM is a software that simulates bone remodelling in terms of Basic Multi-cellular Units (BMUs). 3D visualization of trabecular bone and its attributes is an essential tool in understanding this remodelling process for cancellous bone. It enables the bone researcher to quickly understand the dynamic behaviour of remodelling, the resulting geometry of the bone structure and it allows alternative remodelling scenarios to be compared. This paper presents the volume visualization technique that has been developed to provide this visualization tool.

Smoothing of pixelated finite element models of cancellous bone structures and the effect on the predicted structural properties of the bone.

Many areas of biomedical engineering involve the modelling of biological systems, often using data from medical scanning techniques such as computed microtomography (microCT), and the prediction of the mechanical properties of these systems via finite element models. These models, and also those produced from remodelling simulations on idealized bone structures, are inherently highly pixelated and therefore have a high degree of surface roughness. The purpose of this paper is to demonstrate that this surface roughness need not necessarily have an influence on the predicted properties of the object under examination. To demonstrate this, two-dimensional idealized models of cancellous bone structures were used that were initially depleted and then rebuilt stochastically. A hysteresis effect was observed such that a significant amount of rebuilding beyond the original density was required to regain the initial intact stiffness. To ensure that this effect was not an artefact of the high degree of surface roughness of the rebuilt structures, a two-stage smoothing procedure was applied to assess if this had any effect on the stiffness of the structures. The superpixelation of the structures appeared to have a more profound effect than the smoothing procedures, although the smoothed structures still had stiffness and density values similar to those of the original structures, with a hysteresis effect still evident. This proves that the pixelization of the structures does not have a significant effect on the predicted mechanical properties of the structures. This work has important implications for other models that exhibit a high degree of surface roughness.

Stochastically simulated assessment of anabolic treatment following varying degrees of cancellous bone resorption.

The aim of this study was to investigate the recovery in cancellous bone stiffness resulting from anabolic treatment following varying degrees of resorption, using a stochastic simulation applied to a simplistic structure consisting of five vertical and five horizontal trabeculae. The structure was initially resorbed, and "bone" elements were stochastically removed until nominal resorptions of 10%, 15%, 20%, 25%, and 30% were achieved. A stochastic simulation of anabolic treatment was then applied where bone elements were added, continuing until the original stiffness had been regained, for example, simulating treatment of a patient with an anabolic agent after a period of postmenopausal resorption. The resorption and anabolic simulations were repeated three times for each of the nominal resorptions. The stiffness of the bone structure decreased linearly with resorption, with a slope of approximately -2 and an R(2) of 97.0%; hence, the stiffness fell at approximately twice the rate of the reduction in density. When the various structures regained their original density, the resultant stiffness also had a linear relationship with the original resorption, with a slope of -1 and a lower R(2) of 86.1%. This implies that the reduction in stiffness, when original density was regained, fell proportionately with the degree of initial resorption and, therefore, after a resorption of 30%, when original density was regained, the stiffness of the resultant structure was approximately 30% less than that of the original structure. The density required for the original stiffness to be regained increased linearly with the degree of initial resorption, with a slope of approximately 0.5 and an R(2) of 65.2%, lower than that observed for the previous relationships. This indicates a greater spread of data and suggests greater variability in the formation phase beyond the point of regained original density. Because irreversible connectivity reduction is widely considered to be one of the earliest manifestations of estrogen loss, these findings, although obtained on a simulation of a simplistic cancellous bone structure, support the concept of early intervention to prevent potentially irreversible deterioration of trabecular architecture after menopause.

Finite Element Analysis of Simulations of Cancellous Bone Resorption.

A stochastic simulation of the resorption of cancellous bone has been developed and integrated with a finite element model to predict the resultant change in structural properties of bone as bone density decreases. The resorption represents the net imbalance of osteoclast and osteoblast activity that occurs in osteoporosis. A simple lattice structure of trabecular bone is considered, with an examination of the lattice geometry and discretization indicating that just five trabeculae need to be modelled. The results from the analysis show how the mechanical properties of the cancellous bone degrade with osteoporosis and demonstrate how the method can be used to predict the relationships between stiffness and density or porosity.

Dynamic stochastic simulation of cancellous bone resorption.

A stochastic simulation of cancellous bone resorption was developed and applied to a simple two-dimensional lattice structure representing the vertebral body. The simulation is based upon the concept of a basic multicellular unit (BMU) where net resorption (-deltaB.BMU) is considered at bone/marrow surfaces. The cancellous bone structure is defined as a binary matrix with the size of the pixels corresponding to a square element of approximately 20 microm dimension. The simulation considers both the probability that any surface pixel will be activated into a BMU and, if activated, the length of the resorption cavity. The relationship between relative stiffness and density for the simulation was predicted by finite element analysis. The stochastic simulation was iterated eight times with the mechanical properties assessed after each stage. Perforation of a single trabeculae was first observed at step 2, the structure completely lacking connectivity and mechanical integrity by step 8. The slope of the stiffness-porosity graph was greater than unity for the first five steps, but thereafter approached zero because the structure had lost connectivity and effectively collapsed. The eight-step simulation was repeated five times and demonstrated that, although the stiffness/density relationships were similar at the extremes of density, the dependence of stiffness upon density varied. This clearly demonstrates the stochastic nature of the simulation upon cancellous bone structure, and is probably indicative of a significant dependence of mechanical integrity upon perforation effects.

Telescience at the University of California, Berkeley.

The University of California at Berkeley (UCB) is a member of a university consortium involved in telescience testbed activities under the sponsorship of NASA. Our Telescience Testbed Project consists of three experiments using flight hardware being developed for the Extreme Ultraviolet Explorer project at UCB's Space Sciences Laboratory. The first one is a teleoperation experiment investigating remote instrument control using a computer network such as the Internet. The second experiment is an effort to develop a system for operation of a network of remote workstations allowing coordinated software development, evaluation, and use by widely dispersed groups. The final experiment concerns simulation as a method to facilitate the concurrent development of instrument hardware and support software. We describe our progress in these areas.

A 10-microm infrared camera.

An IR camera has been built at the University of California at Berkeley for astronomical observations. The camera has been used primarily for high angular resolution imaging at mid-IR wavelengths. It has been tested at the University of Arizona 61- and 90-in. telescopes near Tucson and the NASA Infrared Telescope Facility on Mauna Kea, HI. In the observations the system has been used as an imager with interference coated and Fabry-Perot filters. These measurements have demonstrated a sensitivity consistent with photon shot noise, showing that the system is limited by the radiation from the telescope and atmosphere. Measurements of read noise, crosstalk, and hysteresis have been made in our laboratory.