Establishing finite element model credibility of a pedicle screw system under compression-bending: An end-to-end example of the ASME V&V 40 standard.

Computational modeling and simulation (CM&S) is a key tool in medical device design, development, and regulatory approval. For example, finite element analysis (FEA) is widely used to understand the mechanical integrity and durability of orthopaedic implants. The ASME V&V 40 standard and supporting FDA guidance provide a framework for establishing model credibility, enabling deeper reliance on CM&S throughout the total product lifecycle. Examples of how to apply the principles outlined in the ASME V&V 40 standard are important to facilitating greater adoption by the medical device community, but few published examples are available that demonstrate best practices. Therefore, this paper outlines an end-to-end (E2E) example of the ASME V&V 40 standard applied to an orthopaedic implant. The objective of this study was to illustrate how to establish the credibility of a computational model intended for use as part of regulatory evaluation. In particular, this study focused on whether a design change to a spinal pedicle screw construct (specifically, the addition of a cannulation to an existing non-cannulated pedicle screw) would compromise the rod-screw construct mechanical performance. This question of interest (?OI) was addressed by establishing model credibility requirements according to the ASME V&V 40 standard. Experimental testing to support model validation was performed using spinal rods and non-cannulated pedicle screw constructs made with medical grade titanium (Ti-6Al-4V ELI). FEA replicating the experimental tests was performed by three independent modelers and validated through comparisons of common mechanical properties such as stiffness and yield force. The validated model was then used to simulate F1717 compression-bending testing on the new cannulated pedicle screw design to answer the ?OI, without performing any additional experimental testing. This E2E example provides a realistic scenario for the application of the ASME V&V 40 standard to orthopedic medical device applications.

Asymmetry between the dorsal and ventral digging valves of the female locust: function and mechanics.

BACKGROUND: The female locust is equipped with unique digging tools, namely two pairs of valves-a dorsal and a ventral-utilized for excavating an underground hole in which she lays her eggs. This apparatus ensures that the eggs are protected from potential predators and provides optimal conditions for successful hatching. The dorsal and the ventral valves are assigned distinct roles in the digging process. Specifically, the ventral valves primarily function as anchors during propagation, while the dorsal valves displace soil and shape the underground tunnel. RESULTS: In this study, we investigated the noticeable asymmetry and distinct shapes of the valves, using a geometrical model and a finite element method. Our analysis revealed that although the two pairs of valves share morphological similarities, they exhibit different 3D characteristics in terms of absolute size and structure. We introduced a structural characteristic, the skew of the valve cross-section, to quantify the differences between the two pairs of valves. Our findings indicate that these structural variations do not significantly contribute to the valves' load-bearing capabilities under external forces. CONCLUSIONS: The evolutionary development of the form of the female locust digging valves is more aligned with fitting their respective functions rather than solely responding to biomechanical support needs. By understanding the intricate features of these locust valves, and using our geometrical model, valuable insights can be obtained for creating more efficient and specialized tools for various digging applications.

Measuring Electrical Resistivity at the Nanoscale in Phase-Change Materials.

Electrical resistivity is the key parameter in the active regions of many current nanoscale devices, from memristors to resistive random-access memory and phase-change memories. The local resistivity of the materials is engineered on the nanoscale to fit the performance requirements. Phase-change memories, for example, rely on materials whose electrical resistance increases dramatically with a change from a crystalline to an amorphous phase. Electrical characterization methods have been developed to measure the response of individual devices, but they cannot map the local resistance across the active area. Here, we propose a method based on operando electron holography to determine the local resistance within working devices. Upon switching the device, we show that electrical resistance is inhomogeneous on the scale of only a few nanometers.

Pre-intervention myocardial stress is a good predictor of aortic valvoluplasty outcome for fetal critical aortic stenosis and evolving HLHS.

Fetal critical aortic stenosis with evolving hypoplastic left heart syndrome (CAS-eHLHS) causes biomechanical and functional aberrations, leading to a high risk of progression to hypoplastic left heart syndrome (HLHS) at birth. Fetal aortic valvuloplasty (FAV) can resolve outflow obstruction and may reduce progression risk. However, it is currently difficult to accurately predict which patients will respond to the intervention and become functionally biventricular (BV) at birth, as opposed to becoming functionally univentricular (UV). This prediction is important for patient selection, parental counselling, and surgical planning. Therefore, we investigated whether biomechanics parameters from pre-FAV image-based computations could robustly distinguish between CAS-eHLHS cases with BV or UV outcomes in a retrospective cohort. To do so we performed image-based finite element biomechanics modelling of nine CAS-eHLHS cases undergoing intervention and six healthy fetal control hearts, and found that a biomechanical parameter, peak systolic myofibre stress, showed a uniquely large difference between BV and UV cases, which had a larger magnitude effect than echocardiography parameters. A simplified equation was derived for quick and easy estimation of myofibre stress from echo measurements via principal component analysis. When tested on a retrospective cohort of 37 CAS-eHLHS cases, the parameter outperformed other parameters in predicting UV versus BV outcomes, and thus has a high potential of improving outcome predictions, if incorporated into patient selection procedures. Physiologically, high myocardial stresses likely indicate a healthier myocardium that can withstand high stresses and resist pathological remodelling, which can explain why it is a good predictor of BV outcomes. KEY POINTS: Predicting the morphological birth outcomes (univentricular versus biventricular) of fetal aortic valvuloplasty for fetal aortic stenosis with evolving HLHS is important for accurate patient selection, parental counselling and management decisions. Computational simulations show that a biomechanics parameter, pre-intervention peak systolic myofibre stress, is uniquely robust in distinguishing between such outcomes, outperforming all echo parameters. An empirical equation was developed to quickly compute peak systolic myofibre stress from routine echo measurements and was the best predictor of outcomes among a wide range of parameters tested.

CutFEM-based MEG forward modeling improves source separability and sensitivity to quasi-radial sources: A somatosensory group study.

Source analysis of magnetoencephalography (MEG) data requires the computation of the magnetic fields induced by current sources in the brain. This so-called MEG forward problem includes an accurate estimation of the volume conduction effects in the human head. Here, we introduce the Cut finite element method (CutFEM) for the MEG forward problem. CutFEM's meshing process imposes fewer restrictions on tissue anatomy than tetrahedral meshes while being able to mesh curved geometries contrary to hexahedral meshing. To evaluate the new approach, we compare CutFEM with a boundary element method (BEM) that distinguishes three tissue compartments and a 6-compartment hexahedral FEM in an n = 19 group study of somatosensory evoked fields (SEF). The neural generators of the 20 ms post-stimulus SEF components (M20) are reconstructed using both an unregularized and a regularized inversion approach. Changing the forward model resulted in reconstruction differences of about 1 centimeter in location and considerable differences in orientation. The tested 6-compartment FEM approaches significantly increase the goodness of fit to the measured data compared with the 3-compartment BEM. They also demonstrate higher quasi-radial contributions for sources below the gyral crowns. Furthermore, CutFEM improves source separability compared with both other approaches. We conclude that head models with 6 compartments rather than 3 and the new CutFEM approach are valuable additions to MEG source reconstruction, in particular for sources that are predominantly radial.

Fluid-Dynamic Modeling of Flow in Embryonic Tissue Indicates That Lymphatic Valve Location Is Not Consistently Determined by the Local Fluid Shear or Its Gradient.

OBJECTIVE: Intravascular lymphatic valves often occur in proximity to vessel junctions. It is commonly held that disturbed flow at junctions is responsible for accumulation of valve-forming cells (VFCs) at these locations as the initial step in valve creation, and the one which explains the association with these sites. However, evidence in favor is largely limited to cell culture experiments. METHODS: We acquired images of embryonic lymphatic vascular networks from day E16.5, when VFC accumulation has started but the developing valve has not yet altered the local vessel geometry, stained for Prox1, which co-localizes with Foxc2. Using finite-element computational fluid mechanics, we simulated the flow through the networks, under conditions appropriate to this early development stage. Then we correlated the Prox1 distributions with the distributions of simulated fluid shear and shear stress gradient. RESULTS: Across a total of 16 image sets, no consistent correlation was found between Prox1 distribution and the local magnitude of fluid shear, or its positive or negative gradient. CONCLUSIONS: This, the first direct semi-empirical test of the localization hypothesis to interrogate the tissue from in vivo at the critical moment of development, does not support the idea that a feature of the local flow determines valve localization.

Coupling finite and boundary element methods to solve the Poisson-Boltzmann equation for electrostatics in molecular solvation.

The Poisson-Boltzmann equation is widely used to model electrostatics in molecular systems. Available software packages solve it using finite difference, finite element, and boundary element methods, where the latter is attractive due to the accurate representation of the molecular surface and partial charges, and exact enforcement of the boundary conditions at infinity. However, the boundary element method is limited to linear equations and piecewise constant variations of the material properties. In this work, we present a scheme that couples finite and boundary elements for the linearised Poisson-Boltzmann equation, where the finite element method is applied in a confined solute region and the boundary element method in the external solvent region. As a proof-of-concept exercise, we use the simplest methods available: Johnson-Nédélec coupling with mass matrix and diagonal preconditioning, implemented using the Bempp-cl and FEniCSx libraries via their Python interfaces. We showcase our implementation by computing the polar component of the solvation free energy of a set of molecules using a constant and a Gaussian-varying permittivity. As validation, we compare against well-established finite difference solvers for an extensive binding energy data set, and with the finite difference code APBS (to 0.5%) for Gaussian permittivities. We also show scaling results from protein G B1 (955 atoms) up to immunoglobulin G (20,148 atoms). For small problems, the coupled method was efficient, outperforming a purely boundary integral approach. For Gaussian-varying permittivities, which are beyond the applicability of boundary elements alone, we were able to run medium to large-sized problems on a single workstation. The development of better preconditioning techniques and the use of distributed memory parallelism for larger systems remains an area for future work. We hope this work will serve as inspiration for future developments that consider space-varying field parameters, and mixed linear-nonlinear schemes for molecular electrostatics with implicit solvent models.

Ultrathin Ionic Diodes with Electrostatically Heterogeneous Hybrid Interfaces of Nanoporous SiO2 Nanofilms and Polymer Layer-by-Layer Multilayers.

Nanofluidic ionic diodes have attracted much attention due to their unique functions as unidirectional ion transportation ability and promising applications from molecular sensing, and energy harvesting to emerging neuromorphic devices. However, it remains a challenge to fabricate diode-like nanofluidic systems with ultrathin film thickness <100 nm. Herein the formation of ultrathin ionic diodes from hybrid nanoassemblies of nanoporous (NP) SiO2 nanofilms and polyelectrolyte layer-by-layer (LbL) multilayers is described. Ultrathin ionic diodes are prepared by integrating polyelectrolyte multilayers onto photo-oxidized NP SiO2 nanofilms obtained from silsesquioxane-containing block copolymer thin films as a template. The obtained ultrathin ionic diodes exhibit ion current rectification (ICR) properties with high ICR factor = ≈20 under low ionic strength and asymmetric pH conditions. It is concluded that this ICR behavior arises from effective ion accumulation and depletion at the interface of NP SiO2 nanofilms and LbL multilayers attributed to high ion selectivity by combining the experimental data and theoretical calculations using finite element methods. These results demonstrate that the hybrid nano assemblies of NP SiO2 nanofilms and polyelectrolyte LbL multilayers have potential applications for (bio)sensing materials and integrated ionic circuits for seamless connection of human-machine interfaces.

Space-Confined Growth for Thickness-Controlled Cs3Bi2I9 Perovskite Single Crystal Wafers for X-Ray Detectors.

The Cs3Bi2I9 single crystal, as an all-inorganic non-lead perovskite, offers advantages such as stability and environmental friendliness. Its superior photoelectric properties, attributed to the absence of grain boundary influence, make it an outstanding X-ray detection material compared to polycrystals. In addition to material properties, X-ray detector performance is affected by the thickness of the absorption layer. Addressing this, a space-confined method is proposed. The temperature field is determined through finite element simulation, effectively guiding the design of the space-confined method. Through this innovative method, a series of thickness-controlled perovskite single crystal wafers (PSCWs) are successfully prepared. Corresponding X-ray detectors are then prepared, and the impact of single crystal thickness on device performance is investigated. With an increase in single crystal thickness, a rise followed by a decline in device sensitivity is observed, reaching an optimal value at 0.7 mm thickness at 40V mm-1 with a device performance of 11313.6µC Gy-1 cm-2. This space-confined method enables the direct growth of high-quality perovskite single crystals with specified thickness, eliminating the need for slicing or etching.

Functional adaptation of the infant craniofacial system to mechanical loadings arising from masticatory forces.

The morphology and biomechanics of infant crania undergo significant changes between the pre- and post-weaning phases due to increasing loading of the masticatory system. The aims of this study were to characterize the changes in muscle forces, bite forces and the pattern of mechanical strain and stress arising from the aforementioned forces across crania in the first 48 months of life using imaging and finite element methods. A total of 51 head computed tomography scans of normal individuals were collected and analysed from a larger database of 217 individuals. The estimated mean muscle forces of temporalis, masseter and medial pterygoid increase from 30.9 to 87.0 N, 25.6 to 69.6 N and 23.1 to 58.9 N, respectively (0-48 months). Maximum bite force increases from 90.5 to 184.2 N (3-48 months). There is a change in the pattern of strain and stress from the calvaria to the face during postnatal development. Overall, this study highlights the changes in the mechanics of the craniofacial system during normal development. It further raises questions as to how and what level of changes in the mechanical forces during the development can alter the morphology of the craniofacial system.

Hidden functional complexity in the flora of an early land ecosystem.

The first land ecosystems were composed of organisms considered simple in nature, yet the morphological diversity of their flora was extraordinary. The biological significance of this diversity remains a mystery largely due to the absence of feasible study approaches. To study the functional biology of Early Devonian flora, we have reconstructed extinct plants from fossilised remains in silico. We explored the morphological diversity of sporangia in relation to their mechanical properties using finite element method. Our approach highlights the impact of sporangia morphology on spore dispersal and adaptation. We discovered previously unidentified innovations among early land plants, discussing how different species might have opted for different spore dispersal strategies. We present examples of convergent evolution for turgor pressure resistance, achieved by homogenisation of stress in spherical sporangia and by torquing force in Tortilicaulis-like specimens. In addition, we show a potential mechanism for stress-assisted sporangium rupture. Our study reveals the deceptive complexity of this seemingly simple group of organisms. We leveraged the quantitative nature of our approach and constructed a fitness landscape to understand the different ecological niches present in the Early Devonian Welsh Borderland flora. By connecting morphology to functional biology, these findings facilitate a deeper understanding of the diversity of early land plants and their place within their ecosystem.

Association of dorsal malunion in distal radius fractures with wrist osteoarthritis: Alterations of bone density and stress-distribution patterns in relation to deformation angles.

OBJECTIVE: Distal radius fractures (DRFs) with dorsal malunion increase the risk of osteoarthritis (OA), although the cause of post-DRF OA is yet to be elucidated. To clarify the abnormal effects of a post-DRF dorsal radius deformity, we evaluated the bone density (BD) and stress-distribution patterns of the articular surface in dorsally malunited DRFs. DESIGN: In 36 cases of dorsally malunited DRFs following extra-articular fractures, we generated three-dimensional computerized models of the malunited distal radius from computed tomography data and extracted the subchondral bones of the radiocarpal joint (RCJ) and distal radioulnar joint (DRUJ). Both BD and stress distribution in the subchondral bones were quantitatively evaluated by comparing the affected and normal sides. Correlations of alterations in high-BD distribution and deformation angles were analyzed. RESULTS: The center of high-BD distribution from the center of the RCJ in the volar(-)-dorsal(+) direction was dorsal (0.56 ± 0.72 mm) on the affected side compared with the normal side (-0.15 ± 0.63 mm) [95% CI: 0.43, 1.00, P < 0.0001]. The maximum stress distribution was also dorsal on the affected side (2.34 ± 3.52 mm) compared with the normal side (-2.49 ± 1.62 mm) [95% CI: 0.89, 1.79, P < 0.0001]. The alterations in BD and stress distribution correlated with the dorsiflexion and radial deviation angles. In the DRUJ, there was no significant difference in BD between the affected and normal sides. CONCLUSIONS: In dorsally malunited DRFs, the alignment change of the RCJ resulted in high BD-concentration areas and stress distribution on the dorsal side of the radius, which may constitute a precursor for OA.

Bite hard: Linking cranial loading mechanics to ecological differences in gnawing behavior in caviomorph rodents.

The mammalian skull is very malleable and has notably radiated into highly diverse morphologies, fulfilling a broad range of functional needs. Although gnawing is relatively common in mammals, this behavior and its associated morphology are diagnostic features for rodents. These animals possess a very versatile and highly mechanically advantageous masticatory apparatus, which, for instance, allowed caviomorph rodents to colonize South America during the Mid-Eocene and successfully radiate in over 200 extant species throughout most continental niches. Previous work has shown that differences in bite force within caviomorphs could be better explained by changes in muscle development than in mechanical advantages (i.e., in cranial overall morphology). Considering the strong bites they apply, it is interesting to assess how the reaction forces upon the incisors (compression) and the powerful adductor musculature pulling (tension) mechanically affect the cranium, especially between species with different ecologies (e.g., chisel-tooth digging). Thus, we ran finite element analyses upon crania of the subterranean Talas' tuco-tuco Ctenomys talarum, the semi-fossorial common degu Octodon degus, and the saxicolous long-tailed chinchilla Chinchilla lanigera to simulate: (A) in vivo biting in all species, and (B) rescaled muscle forces in non-ctenomyid rodents to match those of the tuco-tuco. Results show that the stress patterns correlate with the mechanical demands of distinctive ecologies, on in vivo-based simulations, with the subterranean tuco-tuco being the most stressed species. In contrast, when standardizing all three species (rescaled models), non-ctenomyid models exhibited a several-fold increase in stress, in both magnitude and affected areas. Detailed observations evidenced that this increase in stress was higher in lateral sections of the snout and, mainly, the zygomatic arch; between approximately 2.5-3.5 times in the common degu and 4.0-5.0 times in the long-tailed chinchilla. Yet, neither species, module, nor simulation condition presented load factor levels that would imply structural failure by strong, incidental biting. Our results let us conclude that caviomorphs have a high baseline for mechanical strength of the cranium because of the inheritance of a very robust "rodent" model, while interspecific differences are associated with particular masticatory habits and the concomitant level of development of the adductor musculature. Especially, the masseteric and zygomaticomandibular muscles contribute to >80% of the bite force, and therefore, their contraction is responsible for the highest strains upon their origin sites, that is, the zygomatic arch and the snout. Thus, the robust crania of the subterranean and highly aggressive tuco-tucos allow them to withstand much stronger forces than degus or chinchillas, such as the ones produced by their hypertrophied jaw adductor muscles or imparted by the soil reaction.

The impact of variations in subject geometry, respiration and coil repositioning on the specific absorption rate in parallel transmit abdominal imaging at 7 T.

Parallel transmit MRI at 7 T has increasingly been adopted in research projects and provides increased signal-to-noise ratios and novel contrasts. However, the interactions of fields in the body need to be carefully considered to ensure safe scanning. Recent advances in physically flexible body coils have allowed for high-field abdominal imaging, but the effects of increased variability on energy deposition need further exploration. The aim of this study was to assess the impact of subject geometry, respiration phase and coil positioning on the specific absorption rate (SAR). Ten healthy subjects (body mass index [BMI] = 25 ± 5 kg m-2 ) were scanned (at 3 T) during exhale breath-hold and images used to generate body models. Seven of these subjects were also scanned during inhale. Simplifications of the coil and body models were first explored, and then finite-difference time-domain simulations were run with a typical eight-channel parallel transmit coil positioned over the abdomen. Simulations were used to generate 10 g averaged SAR (SAR10g ) maps across 100,000 phase settings, and the worst-case scenario 10 g averaged SAR (wocSAR10g ) was identified using trigonometric maximisation. The average maximum SAR10g across the 10 subjects with 1 W input power per channel was 1.77 W kg-1 . Hotspots were always close to the body surface near the muscle wall boundary. The wocSAR10g across the 10 subjects ranged from 2.3 to 3.2 W kg-1 and was inversely correlated to fat volume percentage (R = 8) and BMI (R = 0.6). The coefficient of variation values in SAR10g due to variations in subject geometry, respiration phase and realistic coil repositioning were 12%, 4% and 12%, respectively. This study found that the variability due to realistic coil repositioning was similar to the variability due to differing healthy subject geometries for abdominal imaging. This is important as it suggests that population-based modelling is likely to be more useful than individual modelling in setting safe thresholds for abdominal imaging.

A new hip fracture risk index derived from FEA-computed proximal femur fracture loads and energies-to-failure.

Hip fracture risk assessment is an important but challenging task. Quantitative CT-based patient-specific finite element (FE) analysis (FEA) incorporates bone geometry and bone density in the proximal femur. We developed a global FEA-computed fracture risk index to increase the prediction accuracy of hip fracture incidence. PURPOSE: Quantitative CT-based patient-specific finite element (FE) analysis (FEA) incorporates bone geometry and bone density in the proximal femur to compute the force (fracture load) and energy necessary to break the proximal femur in a particular loading condition. The fracture loads and energies-to-failure are individually associated with incident hip fracture, and provide different structural information about the proximal femur. METHODS: We used principal component analysis (PCA) to develop a global FEA-computed fracture risk index that incorporates the FEA-computed yield and ultimate failure loads and energies-to-failure in four loading conditions of 110 hip fracture subjects and 235 age- and sex-matched control subjects from the AGES-Reykjavik study. Using a logistic regression model, we compared the prediction performance for hip fracture based on the stratified resampling. RESULTS: We referred the first principal component (PC1) of the FE parameters as the global FEA-computed fracture risk index, which was the significant predictor of hip fracture (p-value < 0.001). The area under the receiver operating characteristic curve (AUC) using PC1 (0.776) was higher than that using all FE parameters combined (0.737) in the males (p-value < 0.001). CONCLUSIONS: The global FEA-computed fracture risk index increased hip fracture risk prediction accuracy in males.

CT image-based biomarkers for opportunistic screening of osteoporotic fractures: a systematic review and meta-analysis.

The use of opportunistic computed tomography (CT) image-based biomarkers may be a low-cost strategy for screening older individuals at high risk for osteoporotic fractures and populations that are not sufficiently targeted. This review aimed to assess the discriminative ability of image-based biomarkers derived from existing clinical routine CT scans for hip, vertebral, and major osteoporotic fracture prediction. A systematic search in PubMed MEDLINE, Embase, Cochrane, and Web of Science was conducted from the earliest indexing date until July 2023. The evaluation of study quality was carried out using a modified Quality Assessment Tool for Diagnostic Accuracy Studies (QUADAS-2) checklist. The primary outcome of interest was the area under the curve (AUC) and its corresponding 95% confidence intervals (CIs) obtained for four main categories of biomarkers: areal bone mineral density (BMD), image attenuation, volumetric BMD, and finite element (FE)-derived biomarkers. The meta-analyses were performed using random effects models. Sixty-one studies were included in this review, among which 35 were synthesized in a meta-analysis and the remaining articles were qualitatively synthesized. In comparison to the pooled AUC of areal BMD (0.73 [95% CI 0.71-0.75]), the pooled AUC values for predicting osteoporotic fractures for FE-derived parameters (0.77 [95% CI 0.72-0.81]; p < 0.01) and volumetric BMD (0.76 [95% CI 0.71-0.81]; p < 0.01) were significantly higher, but there was no significant difference with the pooled AUC for image attenuation (0.73 [95% CI 0.66-0.79]; p = 0.93). Compared to areal BMD, volumetric BMD and FE-derived parameters may provide a significant improvement in the discrimination of osteoporotic fractures using opportunistic CT assessments.

A self-stiffening compliant intracortical microprobe.

Utilising a flexible intracortical microprobe to record/stimulate neurons minimises the incompatibility between the implanted microprobe and the brain, reducing tissue damage due to the brain micromotion. Applying bio-dissolvable coating materials temporarily makes a flexible microprobe stiff to tolerate the penetration force during insertion. However, the inability to adjust the dissolving time after the microprobe contact with the cerebrospinal fluid may lead to inaccuracy in the microprobe positioning. Furthermore, since the dissolving process is irreversible, any subsequent positioning error cannot be corrected by re-stiffening the microprobe. The purpose of this study is to propose an intracortical microprobe that incorporates two compressible structures to make the microprobe both adaptive to the brain during operation and stiff during insertion. Applying a compressive force by an inserter compresses the two compressible structures completely, resulting in increasing the equivalent elastic modulus. Thus, instant switching between stiff and soft modes can be accomplished as many times as necessary to ensure high-accuracy positioning while causing minimal tissue damage. The equivalent elastic modulus of the microprobe during operation is ≈ 23 kPa, which is ≈ 42% less than the existing counterpart, resulting in ≈ 46% less maximum strain generated on the surrounding tissue under brain longitudinal motion. The self-stiffening microprobe and surrounding neural tissue are simulated during insertion and operation to confirm the efficiency of the design. Two-photon polymerisation technology is utilised to 3D print the proposed microprobe, which is experimentally validated and inserted into a lamb's brain without buckling.

Strain in Hybrid Organic-Inorganic Metal Halide Perovskites Microstructures by Numerical Simulations.

Hybrid organic-inorganic metal halide perovskites (HOIPs) are promising materials for optoelectronics applications. Their optical and electrical properties can be controlled by strain engineering, that results from application of local elastic deformation or deposition on pre-patterned substrates acquiring a conformal 3D shape. Most interesting, their mechanical properties depend on their crystal structure, composition and dimensionality. We explore by numerical simulations the deformation of a selection of HOIPs comprising a broad range of elastic properties. We consider an axial symmetry with the formation of microdomes on flakes. Radial and vertical forces are considered, finding that the radial force is more effective to obtain large deformation. Large vertical displacement and strain is obtained for HOIPs with low stiffness. The layered nature of HOIPs, that are formed by inorganic layers of different thickness and organic spacers, is also investigated, revealing a non-monotonous trend with the proportion of inorganic to organic part.

A new biomechanical model of the mammal jaw based on load path analysis.

The primary function of the tetrapod jaw is to transmit jaw muscle forces to bite points. The routes of force transfer in the jaw have never been studied, but can be quantified using load paths--the shortest, stiffest routes from regions of force application to support constraints. Here we use load path analysis to map force transfer from muscle attachments to bite point and jaw joint, and to evaluate how different configurations of trabecular and cortical bone affect load paths. We created three models of the mandible of the Virginia opossum, Didelphis virginiana, each with a cortical bone shell, but with different material properties for the internal spaces: a cortical-trabecular model, in which the interior space is modeled with bulk properties of trabecular bone; a cortical-hollow model, in which trabeculae and mandibular canal are modeled as hollow; and a solid-cortical model, in which the interior is modeled as cortical bone. The models were compared with published in vivo bite force and bone strain data, and the load paths calculated for each model. The cortical-trabecular model, which most closely approximates the actual morphology, was best validated by in vivo data. In all three models the load path was confined to cortical bone, although its route within the cortex varied depending on the material properties of the inner model. Our analysis shows that most of the force is transferred through the cortical, rather than trabecular bone, and highlights the potential of load path analysis for understanding form-function relationships in the skeleton.

Tibial strains are sensitive to speed perturbations, but not grade perturbations, during running.

A fatigue-failure process is hypothesized to govern the development of tibial stress fractures, where bone damage is highly dependent on the peak strain magnitude. To date, much of the work examining tibial strain during running has ignored uphill and downhill running despite the prevalence of this terrain. This study examined the sensitivity of tibial strain to changes in running grade and speed using a combined musculoskeletal-finite element modelling routine. Seventeen participants ran on a treadmill at ±10, ±5 and 0 deg; at each grade, participants ran at 3.33 m s-1 and at a grade-adjusted speed of 2.50 and 4.17 m s-1 for uphill and downhill grades, respectively. Force and motion data were recorded in each grade and speed combination. Muscle and joint contact forces were estimated using inverse-dynamics-based static optimization. These forces were applied to a participant-adjusted finite element model of the tibia. None of the strain variables (50th and 95th percentile strain and strained volume ≥4000 µÎµ) differed as a function of running grade; however, all strain variables were sensitive to running speed (F1≥9.59, P≤0.03). In particular, a 1 m s-1 increase in speed resulted in a 9% (∼260 µÎµ) and 155% (∼600 mm3) increase in peak strain and strained volume, respectively. Overall, these findings suggest that faster running speeds, but not changes in running grade, may be more deleterious to the tibia.