The Virtual Family computational whole-body anatomical human models were originally developed for electromagnetic (EM) exposure evaluations, in particular to study how absorption of radiofrequency radiation from external sources depends on anatomy. However, the models immediately garnered much broader interest and are now applied by over 300 research groups, many from medical applications research fields. In a first step, the Virtual Family was expanded to the Virtual Population to provide considerably broader population coverage with the inclusion of models of both sexes ranging in age from 5 to 84 years old. Although these models have proven to be invaluable for EM dosimetry, it became evident that significantly enhanced models are needed for reliable effectiveness and safety evaluations of diagnostic and therapeutic applications, including medical implants safety. This paper describes the research and development performed to obtain anatomical models that meet the requirements necessary for medical implant safety assessment applications. These include implementation of quality control procedures, re-segmentation at higher resolution, more-consistent tissue assignments, enhanced surface processing and numerous anatomical refinements. Several tools were developed to enhance the functionality of the models, including discretization tools, posing tools to expand the posture space covered, and multiple morphing tools, e.g., to develop pathological models or variations of existing ones. A comprehensive tissue properties database was compiled to complement the library of models. The results are a set of anatomically independent, accurate, and detailed models with smooth, yet feature-rich and topologically conforming surfaces. The models are therefore suited for the creation of unstructured meshes, and the possible applications of the models are extended to a wider range of solvers and physics. The impact of these improvements is shown for the MRI exposure of an adult woman with an orthopedic spinal implant. Future developments include the functionalization of the models for specific physical and physiological modeling tasks.

Epidural electrical stimulation (EES) targeting the dorsal roots of lumbosacral segments restored walking in people with spinal cord injury (SCI). However, EES was delivered with multielectrode paddle leads that were originally designed to target the dorsal column of the spinal cord. Here, we hypothesized that an arrangement of electrodes targeting the ensemble of dorsal roots involved in leg and trunk movements would result in superior efficacy, restoring more diverse motor activities after the most severe SCI. To test this hypothesis, we established a computational framework that informed the optimal arrangement of electrodes on a new paddle lead and guided its neurosurgical positioning. We also developed a software supporting the rapid configuration of activity-specific stimulation programs that reproduced the natural activation of motor neurons underlying each activity. We tested these neurotechnologies in three individuals with complete sensorimotor paralysis, as part of an ongoing clinical trial (clinicaltrials.gov, NCT02936453). Within a single day, activity-specific stimulation programs enabled the three individuals to stand, walk, cycle, swim, and control trunk movements. Neurorehabilitation mediated sufficient improvement to restore these activities in community settings, opening a realistic path to support everyday mobility with EES in people with SCI.

Computational modeling and simulations are increasingly being used to complement experimental testing for analysis of safety and efficacy of medical devices. Multiple voxel- and surface-based whole- and partial-body models have been proposed in the literature, typically with spatial resolution in the range of 1–2 mm and with 10–50 different tissue types resolved. We have developed a multimodal imaging-based detailed anatomical model of the human head and neck, named “MIDA”. The model was obtained by integrating three different magnetic resonance imaging (MRI) modalities, the parameters of which were tailored to enhance the signals of specific tissues: i) structural T1- and T2-weighted MRIs; a specific heavily T2-weighted MRI slab with high nerve contrast optimized to enhance the structures of the ear and eye; ii) magnetic resonance angiography (MRA) data to image the vasculature, and iii) diffusion tensor imaging (DTI) to obtain information on anisotropy and fiber orientation. The unique multimodal high-resolution approach allowed resolving 153 structures, including several distinct muscles, bones and skull layers, arteries and veins, nerves, as well as salivary glands. The model offers also a detailed characterization of eyes, ears, and deep brain structures. A special automatic atlas-based segmentation procedure was adopted to include a detailed map of the nuclei of the thalamus and midbrain into the head model. The suitability of the model to simulations involving different numerical methods, discretization approaches, as well as DTI-based tensorial electrical conductivity, was examined in a case-study, in which the electric field was generated by transcranial alternating current stimulation. The voxel- and the surface-based versions of the models are freely available to the scientific community.

Mass spring models are frequently used to simulate deformable objects because of their conceptual simplicity and computational speed. Unfortunately, the model parameters are not related to elastic material constitutive laws in an obvious way. Several methods to set optimal parameters have been proposed, but so far only with limited success. We analyze the parameter identification problem and show the difficulties, which have prevented previous work from reaching wide usage. Our main contribution is a new method to derive analytical expressions for the spring parameters from an isotropic linear elastic reference model. The method is described and expressions for several mesh topologies are derived. These include triangle, rectangle and tetrahedron meshes. The formulae are validated by comparing the static deformation of the MSM with reference deformations simulated with the finite element method.

Over the past decades, significant improvements have been made in the field of computational human phantoms (CHPs) and their applications in biomedical engineering. Their sophistication has dramatically increased. The very first CHPs were composed of simple geometric volumes, e.g., cylinders and spheres, while current CHPs have a high resolution, cover a substantial range of the patient population, have high anatomical accuracy, are poseable, morphable, and are augmented with various details to perform functionalized computations. Advances in imaging techniques and semi-automated segmentation tools allow fast and personalized development of CHPs. These advances open the door to quickly develop personalized CHPs, inherently including the disease of the patient. Because many of these CHPs are increasingly providing data for regulatory submissions of various medical devices, the validity, anatomical accuracy, and availability to cover the entire patient population is of utmost importance. The article is organized into two main sections: the first section reviews the different modeling techniques used to create CHPs, whereas the second section discusses various applications of CHPs in biomedical engineering. Each topic gives an overview, a brief history, recent developments, and an outlook into the future.