Apparent Modulus Research Articles

Anomalously High Apparent Young's Modulus of Ultrathin Freestanding PZT Films Revealed by Machine Learning Empowered Nanoindentation.

The mechanical properties at small length scales are not only significant for understanding the intriguing size-dependent behaviors but also critical for device applications. Nanoindentation via atomic force microscopy is widely used for small-scale mechanical testing, yet determining the Young's modulus of quasi-2D films from freestanding force-displacement curve of nanoindentation remains challenging, complicated by both bending and stretching that are highly nonlinear. To overcome these difficulties, a machine learning model is developed based on the back propagation (BP) neural network and finite element training to accurately determine the Young's modulus, pretension, and thickness of freestanding films from nanoindentation force-displacement curves simultaneously, improving the computational efficiency by two orders of magnitude over conventional brute force curve fitting. Using this technique, anomalously high apparent Young's modulus is discovered of 2.8nm thick PbZr0.2Ti0.8O3 (PZT) film as large as 229.4±12.1GPa, much higher than the bulk value. The enhancement can be attributed to the strain gradient-induced flexoelectric effect, and the corresponding flexoelectric coefficient is estimated to be ≈200nCm-1. The method is developed to enable an artificial intelligence (AI) system to study the mechanical properties of a wide range of low-dimensional materials.

Poroelastic Response of a Fractured Rock to Hydrostatic Pressure Oscillations

AbstractPoroelastic coupling between fractures and the surrounding rock is important to numerous applications in geosciences. We measure the in‐situ fluid pressure and local strain response of a fractured carbonate sample to hydrostatic pressure oscillations. A linear poroelastic model that represents the rock sample is parameterized using X‐ray imaging and ultrasonic wave transmission measurements. The numerical solution, based on Biot's quasistatic equations, is consistent with the measured frequency dependent dispersion of the apparent bulk modulus of the background matrix and the in‐situ pore pressure response, which is caused by fluid pressure diffusion from the compliant fractures into the stiffer matrix. The observed fluid pressure diffusion is causally related to the numerically quantified intrinsic attenuation at seismic frequencies, which is a major contributor to the dissipation of seismic waves. Our analysis supports the use of a simple approximation of fractures as compliant and planar inclusions in numerical simulations based on linear poroelasticity.

Insight of soy protein isolate to decrease the gel properties corn starch based binary system: Rheological and structural investigation

Dysphagia has become a growing concern for elderly individuals, and starch-based gel foods is expected to become one of the foods easy to swallow. This study investigated the effects of soy protein isolate (SPI) on the sol-gel behavior, rheological, structural and gel properties of corn starch (CS). SPI raised the gelatinization temperature but reduced enthalpy values, inhibiting the gelatinization of CS. There was a reduction of rheological properties in CS-SPI system, involving apparent viscosity, modulus, and shear structural recovery properties, while SPI increased the Tan δ, creep-recovery, and thixotropic properties. Structural analysis demonstrated that SPI introduction altered the lamellar structure, increased the relative crystallinity, and decreased the hydrogen bond strength and order degree in CS-SPI gels, forming a weaker gel network. As a result, the texture profiles such as hardness, gumminess, and chewiness were decreased from 308.60 g to 64.98 g, 285.44 g–50.62 g, and 279.66 g–50.15 g, respectively, whereas the gel strength of CS-SPI gel decreased from 123.91 g to 45.78 g. This study provided foundations for the production of gel food suitable for swallowing.

Fatigue analysis on flax-epoxy composites: insights and implications

This study investigates the performance of flax fiber and epoxy resin composites, focusing on fatigue testing, which is scarcely reported in the literature. Fatigue behavior is evaluated under full reverse bending loading conditions: R = −1. Apart of the definition of the classic S-N curve, the analysis deeply investigated the evolution of the hysteresis leaf generated during each cycle over the entire fatigue life. This approach allows evaluating interesting parameters such as the evolution of the loss factor and the apparent modulus over life, being these informations fundamental for practical application of natural-fiber composites in engineering applications.

Fabrication and characterization of electrospun core-shell nanofibers of bilayer zein/pullulan emulsions crosslinked by genipin

In this study, the electrospun core-shell nanofibers of zein/pullulan stabilized bilayer emulsions before and after genipin crosslinking were fabricated. The experimental results indicated that the addition of pullulan increased the apparent viscosity and elastic modulus of the bilayer emulsions, which was further increased after the chemical crosslinking of genipin. The nanofiber diameter increased from 102.9 nm to 169.9 nm with the increasing ratio of pullulan, but increased significantly to a range of 405.6–708.0 nm after genipin crosslinking. The electrospun nanofiber films of crosslinked emulsions had higher thermal stability and stronger water stability. The FTIR result proved the existence of hydrogen bond interaction between the zein, pullulan, and genipin molecules. In addition, before and after crosslinking, the encapsulation efficiency of electrospun fiber films for camellia oil was >77.68 %, and the maximum encapsulation efficiency could reach 87.94 %, and there was no significant change during the 7-day storage period, showing good stability. These research results can provide a theoretical basis for the encapsulation of hydrophobic active substances in zein-based emulsion electrospun core-shell nanofibers.

Rheological modelling of apple puree based on machine learning combined Monte Carlo simulation: Insight into the fundamental light- particle structure interaction processes

In this work, apple purees from different particle concentration and verifying in size were reconstituted to investigate their impacts on rheological behaviors, optical properties and light interaction at 900–1650 nm. The optical scattering of different puree particle concentration varied more intensively than particle size. Based on Monte Carol simulation (MCX), the highest proportion of diffuse reflection energy and absorption were observed in 1050 nm and 1480 nm, respectively. The averaged light propagation of all the photons in purees varied from 259.29 mm to 199.17 mm with the particle concentration from 25 % to 45 %. Finally, support vector machine regression based on optical scattering properties fused with MCX parameters at 1050 nm effectively evaluated puree apparent viscosity parameters (RPD > 2.39), viscous and elastic modulus (RPD > 2.47). These extended our knowledge of the fundamental light - particle structure interaction processes and provided a new solution for rheological modelling in pureed food.

On the influence of structural and chemical properties on the elastic modulus of woven bone under healing.

Woven bone, a heterogeneous and temporary tissue in bone regeneration, is remodeled by osteoblastic and osteoclastic activity and shaped by mechanical stress to restore healthy tissue properties. Characterizing this tissue at different length scales is crucial for developing micromechanical models that optimize mechanical parameters, thereby controlling regeneration and preventing non-unions. This study examines the temporal evolution of the mechanical properties of bone distraction callus using nanoindentation, ash analysis, micro-CT for trabecular microarchitecture, and Raman spectroscopy for mineral quality. It also establishes single- and two-parameter power laws based on experimental data to predict tissue-level and bulk mechanical properties. At the macro-scale, the tissue exhibited a considerable increase in bone fraction, controlled by the widening of trabeculae. The Raman mineral-to-matrix ratios increased to cortical levels during regeneration, but the local elastic modulus remained lower. During healing, the tissue underwent changes in ash fraction and in the percentages of Calcium and Phosphorus. Six statistically significant power laws were identified based on the ash fraction, bone fraction, and chemical and Raman parameters. The microarchitecture of woven bone plays a more significant role than its chemical composition in determining the apparent elastic modulus of the tissue. Raman parameters were demonstrated to provide more significant power laws correlations with the micro-scale elastic modulus than mineral content from ash analysis.

Oligo-cyclic loading-induced evolution of stress distribution and apparent amorphous modulus in lamellar stacks of high-density polyethylene

Assessing the resistance of high-density polyethylene (HDPE) against earthquake-like loads involves understanding the changes in structure and properties induced by oligo-cyclic loading at various length scales. To study the evolution of stress distribution and intrinsic properties within lamellar stacks from pristine to oligo-cyclic loading pre-conditioned materials, simultaneous in-situ SAXS/WAXS measurements were performed. During the elastic deformation of each pristine and preconditioned sample, crystal strain was tracked using the in-situ WAXS technique. Based on the established elastic tensor of the crystalline structure in polyethylene, we calculated the microscopic stress values within crystalline lamellae. In the pristine sample, lamellar stacks exhibit closely series-like coupling in the equatorial region and parallel-like coupling in the polar region of the spherulite. In the pre-conditioned sample, stress is primarily concentrated in the intra-fibrillar region, where the crystalline and amorphous phases are series-coupled, and strong strain concentration occurs in the inter-fibrillar region. By combining the local strain in the amorphous layer within the lamellar stacks in the equatorial region of the spherulite and the intra-fibrillar region with series-coupled lamellar stacks, measured by in-situ SAXS tests, the apparent amorphous modulus at the lamellar stack scale can be determined. This modulus changes from 71 to 106 MPa in the equatorial region of pristine spherulites to a notable 2000–7000 MPa in the intra-fibrillar region under the influence of oligo-cyclic pre-loading. Importantly, this apparent modulus is affected by both crystallization conditions and molecular structure, with molecular parameters exerting the primary influence.

Rheological Characterization of Inorganic Curing Foam for Preventing Spontaneous Combustion of Coal

ABSTRACT This study examines the rheological characteristics of inorganic cured foam (ICF). The apparent viscosity, yield stress, and modulus of elasticity of fresh cement-fly ash composite slurries with varying fly ash blending ratios at water/ash ratios of 0.4, 0.5, 0.6, and 0.7 were determined. The findings indicate that, under identical testing conditions, the apparent viscosity, yield stress, and modulus of elasticity of fresh cement-fly ash composite slurries (FCFS) gradually decrease with increasing water/ash ratio. The apparent viscosity and modulus of elasticity of inorganic curing foam (ICF) decrease as the volume expansion rate increases, while the linear viscoelastic region (LVER) broadens with higher volume expansion rates. Furthermore, the fitted apparent viscosity-shear stress curve for FCFS aligns with the power law model, whereas the apparent viscosity-shear stress curve for ICF slurries is better represented by the Carreau-Yasuda function model. The results of the study will contribute to the industrial application of inorganic curing foams for the prevention of spontaneous coal combustion.

The impacts of different modification techniques on the gel properties and structural characteristics of fish gelatin

Fish gelatin hydrogels are typically characterized by weak texture and poor thermal stability. To address this issue, fish gelatin (FG) was modified using κ-carrageenan (κC) and transglutaminase (TG). FG-κC, FG-TG, and FG-κC-TG gels were prepared. The results demonstrated that the particle size and textural properties of FG gels were significantly enhanced by κC and κC-TG composite modifications, and their zeta potential was reduced (P < 0.05). As opposed to that, the improvement effects of TG modification were weaker. The highest particle size and hardness, along with the lowest zeta potential, were exhibited by the FG-κC-TG0.06% gel. Additionally, the FG-κC-TG gel system possessed higher apparent viscosity and elastic storage modulus. Fluorescence spectroscopy and UV absorption analysis indicated that fluorescence intensity was reduced by all three modification methods, and the conformational structure of FG was altered to varying degrees. Microstructure, FTIR, and XRD results revealed that the number of hydrogen bonds and isopeptide bonds in the composite gel systems was increased by κC and κC-TG modifications, enhancing the stability and orderliness of the gel's helical structure, and forming larger aggregates and denser network structures. Furthermore, due to the presence of covalent cross-linking, improved thermal stability was exhibited by the modified gelatin hydrogels. This study not only provides theoretical support for enhancing the performance of fish gelatin but also facilitates the application of modified fish gelatin in improving the sensory quality and flavor characteristics of food, as well as in serving as a matrix material in the biomedical field.

The role of collagen and crystallinity in the physicochemical properties of naturally derived bone grafts

Abstract Xenografts are commonly used for bone regeneration in dental and orthopaedic domains to repair bone voids and other defects. The first-generation xenografts were made through sintering, which deproteinises them and alters their crystallinity, while later xenografts are produced using cold-temperature chemical treatments to maintain the structural collagen phase. However, the impact of collagen and the crystalline phase on physicochemical properties have not been elucidated. We hypothesized that understanding these factors could explain why the latter provides improved bone regeneration clinically. In this study, we compared two types of xenografts, one prepared through a low-temperature chemical process (Treated) and another subsequently sintered at 1100 °C (Sintered) using advanced microscopy, spectroscopy, x-ray analysis, and compressive testing. Our investigation showed that the Treated bone graft was free of residual blood, lipids, or cell debris, mitigating the risk of pathogen transmission. Meanwhile, the sintering process removed collagen and the carbonate phase of the Sintered graft, leaving only calcium phosphate and increased mineral crystallinity. Micro computed tomography revealed that the Treated graft exhibited an increased high porosity (81%) and pore size compared to untreated bone, whereas the Sintered graft exhibited shrinkage, which reduced the porosity (72%), pore size, and strut size. Additionally, scanning electron microscopy (SEM) displayed crack formation around the pores of the Sintered graft. The Treated graft displayed median mechanical properties comparable to native cancellous bone and clinically available solutions, with an apparent modulus of 166 MPa, yield stress of 5.5 MPa, and yield strain of 4.9%. In contrast, the Sintered graft exhibited a lower median apparent modulus of 57 MPa. It failed in a brittle manner at a median stress of 1.7 MPa and strain level of 2.9%, demonstrating the structural importance of the collagen phase. This indicates why bone grafts prepared through cold-temperature processes are clinically favourable.

Square Wave Electrogravimetry Monitors Transient Submonolayer Adsorption of Redox-Active Molecules with a Precision of Less Than 1% of a Monolayer

Square wave voltammetry on electrolytes containing the reversible redox pairs flavin adenine dinucleotide (FAD) and methyl viologen (MV) was complemented by a new form of acoustic microgravimetry. The analytes FAD and MV have promising applications in bacterial fuel cells and redox flow batteries. The instrument operates similarly to a quartz crystal microbalance with dissipation monitoring (EQCM-D). It reports shifts in resonance frequency and half bandwidth on several overtones. Compared to the existing instruments, the time resolution was improved by a factor of 40 down to 2.5 ms by multifrequency lockin amplification, which is important for the analysis of the transients in electrochemistry. The frequency noise after accumulation was reduced below 10 mHz ≙ 0.12 ng/cm2 ≙ 1.2 pm (assuming ρ = 1 g/cm3) by employing modulation and averaging. Square wave modulation of the electrode potential is superior to linear ramps (cyclic voltammetry) because it reduces the contribution of non-Faraday currents and, also, improves sensitivity. The shifts of frequency,Δf, and half bandwidth, ΔΓ, are in line with the Sauerbrey prediction (Δf/n ≈const. and ΔΓ/n < -Δf/nwith n the overtone order). Small deviations (Figure B and D) presumably are caused by the layer’s softness. The apparent shear modulus, G app, of this layer is in the GPa range. The frequency difference, Δf/n|CT (Figure E), amounts to an apparent mass transfer rate. Plots of Δf/n|CT versus potential sometimes show curves similar to the electric current as observed for e. g. MV. For FAD, the results obtained depended on pH and the apparent mass transfer rates are much different from the current (Figure C and E). This difference is most likely caused by adsorbed molecules that are intermediates of the underlying two-electron process. The explanation is further corroborated by the response time of the QCM signals being much larger than the RC time of double layer recharging as determined with EIS. An interpretation in terms of adsorption/desorption is more plausible than an interpretation in terms of changed viscosity in the diffuse double layer.[1] Caption: Currents and frequency shifts obtained on a solution containing FAD. Df/n|av evidences the preferred adsorption of the intermediate state (FADH).

Formulations and evaluations of structural and physico-chemical properties of soy yogurts: Effect of incorporating soy protein isolate/chitosan complexed microgels

Novel complexed microgels with enhanced pH and ion stability were constructed using heat-denatured soy protein isolate (SPI) and Chitosan (CS). Their characteristics and their effects on the textural, rheological, and tribological properties of soy yogurt were investigated. The complexed microgels showed larger particle sizes, lower surface hydrophobicity, and higher turbidity than pure SPI microgels. Incorporating the SPI-CS complexed microgels into the yogurts increased their protein content from 3.94% to 5.07% and reduced their fat content from 2.17% to 1.99%. Additionally, incorporating the microgels enhanced the texture of the yogurts, including increased firmness, elasticity, and viscosity, alongside improved water holding capacity (WHC). The incorporation of 0.8% SPI-CS microgels led to the optimum physicochemical characteristics. Rheological analyses demonstrated that incorporating SPI-CS microgels led to an increase in apparent viscosity and shear modulus. Furthermore, the presence of the SPI-CS microgels significantly improved the lubrication properties. Microstructural analysis showed subtle differences in the gel structure of the yogurts after microgel addition. The three-dimensional network structure within the microgel-fortified yogurts was more compact than that in the control. Thus, SPI-CS microgels potentially improve the nutritional profile and quality attributes of plant-based yogurt products.

Scale-inspired programmable robotic structures with concurrent shape morphing and stiffness variation.

Biological organisms often have remarkable multifunctionality through intricate structures, such as concurrent shape morphing and stiffness variation in the octopus. Soft robots, which are inspired by natural creatures, usually require the integration of separate modules to achieve these various functions. As a result, the whole structure is cumbersome, and the control system is complex, often involving multiple control loops to finish a required task. Here, inspired by the scales that cover creatures like pangolins and fish, we developed a robotic structure that can vary its stiffness and change shape simultaneously in a highly integrated, compact body. The scale-inspired layered structure (SAILS) was enabled by the inversely designed programmable surface patterns of the scales. After fabrication, SAILS was inherently soft and flexible. When sealed in an elastic envelope and subjected to negative confining pressure, it transitioned to its designated shape and concurrently became stiff. SAILS could be actuated at frequencies as high as 5 hertz and achieved an apparent bending modulus change of up to 53 times between its soft and stiff states. We further demonstrated both the versatility of SAILS by developing a soft robot that is amphibious and adaptive and tunable landing systems for drones with the capacity to accommodate different loads.

Non-Destructive Hardness Indentation Measurement of Residual Stress on Large Aerospace Forged Components at the Engineering Site Based on Impact Hardness Tester.

Large forgings are crucial in aerospace applications; however, the residual stresses generated during their forming and heat treatment seriously affect their serviceability. Therefore, the non-destructive detection of residual stresses in large forgings is of far-reaching significance for ensuring the quality of forgings and realising precision machining. Although a variety of detection methods are available, there is still a lack of a programme that can comprehensively, accurately and non-destructively measure the residual stresses in large forgings. This study is dedicated to exploring the application of the bouncing impact indentation method in the non-destructive testing of residual stresses in large forgings. Through in-depth finite element simulations and orthogonal scheme analyses, we found that the elastic modulus, yield strength and work hardening indexes have significant effects on the impact indentation process. Further, we establish the dimensionless function of residual stress and indentation parameters, and successfully obtain the inversion algorithm of residual stress. The relative error of the calculated values of the indentation curves hm and hr in the simulation with reference values is not more than 3%, and the relative error of the corrected Pm inversion values for most virtual materials is not more than 5%. The folding elastic modulus and apparent elastic modulus obtained by inversion are controlled within 10%, which demonstrates a high value for engineering applications. In addition, we innovatively express the research results in the form of 3D stress diagrams, realising the digital expression of 3D residual stresses in large forgings based on feature point measurements and contour surface configurations, which provides intuitive and comprehensive data support for engineering practice.

Computational comparison study of virtual compression and shear test for estimation of apparent elastic moduli under various boundary conditions.

Virtual compression tests based on finite element analysis are representative noninvasive methods to evaluate bone strength. However, owing to the characteristic porous structure of bones, the material obtained from micro-computed tomography images in the finite-element model is not uniformly distributed. These characteristics cause differences in the apparent elastic moduli depending on the boundary conditions and affect the accuracy of bone-strength evaluation. Therefore, this study aimed to evaluate and compare the apparent elastic moduli under various, virtual-compression and shear-test boundary conditions. Four, nonuniform models were constructed with increasing model complexity. For representative boundary conditions, two, different, testing directions, and constrained surfaces were applied. As a result, the apparent elastic moduli of the nonuniform model varied up to 55.2% based on where the constrained surface was located in the single-end-cemented condition. Additionally, when connectivity in the test direction was lost, the accuracy of the apparent elastic moduli was low. A graphical comparison showed that the equivalent-stress distribution was more advantageous for analyzing load transferability and physical behavior than the strain-energy distribution. These results clearly show that the prediction accuracy of the apparent elastic moduli can be guaranteed if the boundary condition on the constraint and loading surfaces of the nonuniform model are applied symmetrically and the connectivity of the elements in the testing direction is well maintained. This study will aid in precision improvement of bone-strength-indicator determination for osteoporosis prevention.

Mechanical characterization of Xenopus laevis oocytes using atomic force microscopy

Mechanical properties are essential for the biological activities of cells, and they have been shown to be affected by diseases. Therefore, accurate mechanical characterization is important for studying the cell lifecycle, cell-cell interactions, and disease diagnosis. While the cytoskeleton and actin cortex are typically the primary structural stiffness contributors in most live cells, oocytes possess an additional extracellular layer known as the vitelline membrane (VM), or envelope, which can significantly impact their overall mechanical properties. In this study, we utilized nanoindentation via an atomic force microscope to measure the Young's modulus of Xenopus laevis oocytes at different force setpoints and explored the influence of the VM by conducting measurements on oocytes with the membrane removed. The findings revealed that the removal of VM led to a significant decrease in the apparent Young's modulus of the oocytes, highlighting the pivotal role of the VM as the main structural component responsible for the oocyte's shape and stiffness. Furthermore, the mechanical behavior of VM was investigated through finite element (FE) simulations of the nanoindentation process. FE simulations with the VM Young's modulus in the range 20–60 MPa resulted in force-displacement curves that closely resemble experimental in terms of shape and maximum force for a given indentation depth.

Improvement of apparent IFSS and specific modulus of CNT yarns

Due to the high strength of carbon nanotubes (CNTs) significant effort has been devoted to incorporating them into fibers and composites. To maximize their potential in these areas, there must be strong interactions between the CNT structures and the matrix of the composite. For CNT yarn, progress in this area has been limited. In this work, we present a method for modifying CNT yarn which improves the apparent interfacial shear strength (IFSS) by up to 45 %, as measured by single yarn pullout testing. An ammonium peroxide mixture (APM) treatment method is used to modify the yarn, and this is followed by treatment with commercial sizing agents. A tensioned heat treatment method is also employed to limit the effect of APM treatment on strength and modulus. This method increased the specific modulus of CNT yarns by as much as 16 %. The combined heat and chemical treatments yield increases as high as 43 % in IFSS and 6 % in specific tensile modulus while tempering the decreases in specific tensile strength. The yarn is characterized by SEM, EDS and Raman spectroscopy to show surface carbon being removed from the yarn, seemingly without damaging the CNT structures in the yarn. The main mechanism for IFSS increase appears to be removal of non-graphitic material from the yarn and slight oxidation of the yarn surface. The increase seen in IFSS has the potential to improve the performance of CNT yarn/polymer matrix composites for aerospace applications.

Nonlinear Poisson effect in affine semiflexible polymer networks.

Stretching an elastic material along one axis typically induces contraction along the transverse axes, a phenomenon known as the Poisson effect. From these strains, one can compute the specific volume, which generally either increases or, in the incompressible limit, remains constant as the material is stretched. However, in networks of semiflexible or stiff polymers, which are typically highly compressible yet stiffen significantly when stretched, one instead sees a significant reduction in specific volume under finite strains. This volume reduction is accompanied by increasing alignment of filaments along the strain axis and a nonlinear elastic response, with stiffening of the apparent Young's modulus. For semiflexible networks, in which entropic bending elasticity governs the linear elastic regime, the nonlinear Poisson effect is caused by the nonlinear force-extension relationship of the constituent filaments, which produces a highly asymmetric response of the constituent polymers to stretching and compression. The details of this relationship depend on the geometric and elastic properties of the underlying filaments, which can vary greatly in experimental systems. Here, we provide a comprehensive characterization of the nonlinear Poisson effect in an affine network model and explore the influence of filament properties on essential features of both microscopic and macroscopic response, including strain-driven alignment and volume reduction.

Elastic response of trabecular bone under compression calculated using the firm and floppy boundary lattice element method

Micro-Finite Element analysis (μFEA) has become widely used in biomechanical research as a reliable tool for the prediction of bone mechanical properties within its microstructure such as apparent elastic modulus and strength. However, this method requires substantial computational resources and processing time. Here, we propose a computationally efficient alternative to FEA that can provide an accurate estimation of bone trabecular mechanical properties in a fast and quantitative way. A lattice element method (LEM) framework based on the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) open-source software package is employed to calculate the elastic response of trabecular bone cores. A novel procedure to handle pore-material boundaries is presented, referred to as the Firm and Floppy Boundary LEM (FFB-LEM). Our FFB-LEM calculations are compared to voxel- and geometry-based FEA benchmarks incorporating bovine and human trabecular bone cores imaged by micro Computed Tomography (μCT). Using 14 computer cores, the apparent elastic modulus calculation of a trabecular bone core from a μCT-based input with FFB-LEM required about 15 min, including conversion of the μCT data into a LAMMPS input file. In contrast, the FEA calculations on the same system including the mesh generation, required approximately 30 and 50 min for voxel- and geometry-based FEA, respectively. There were no statistically significant differences between FFB-LEM and voxel- or geometry-based FEA apparent elastic moduli (+24.3% or +7.41%, and +0.630% or -5.29% differences for bovine and human samples, respectively).