Traction-separation law parameters for the description of age-related changes in the delamination strength of the human descending thoracic aorta.

Aortic dissection is a life-threatening disease that consists in the development of a tear in the wall of the aorta. The initial tear propagates as a discontinuity leading to separation within the aortic wall, which can result in the creation of a so-called false lumen. A fatal threat occurs if the rupture extends through the whole thickness of the aortic wall, as blood may then leak. It is generally accepted that the dissection, which can sometime extend along the entire length of the aorta, propagates via a delamination mechanism. The aim of the present paper is to provide experimentally validated parameters of a mathematical model for the description of the wall's cohesion. A model of the peeling experiment was built in Abaqus. The delamination interface was described by a piecewise linear traction-separation law. The bulk behavior of the aorta was assumed to be nonlinearly elastic, anisotropic, and incompressible. Our simulations resulted in estimates of the material parameters for the traction-separation law of the human descending thoracic aorta, which were obtained by minimizing the differences between the FEM predictions and the delamination force given by the regression of the peeling experiments. The results show that the stress at damage initiation, Tc, should be understood as an age-dependent quantity, and under the assumptions of our model this dependence can be expressed by linear regression as Tc = - 13.03·10-4·Age + 0.2485 if the crack front advances in the axial direction, and Tc = - 7.58·10-4·Age + 0.1897 if the crack front advances in the direction of the aortic circumference (Tc [MPa], Age [years]). Other model parameters were the stiffness K and the separation at failure, δf-δc (K = 0.5 MPa/mm, δf-δc = 0.1 mm). The material parameters provided by our study can be used in numerical simulations of the biomechanics of dissection propagation through the aorta especially when age-associated phenomena are studied.

Inverse finite element analysis for an axisymmetric model of vertical tooth extraction.

BACKGROUND AND OBJECTIVE: Tooth extraction is a common clinical procedure with biomechanical factors that can directly influence patient outcomes. Recent development in atraumatic extraction techniques have endeavoured to improve treatment outcomes, but the characterization of extraction biomechanics is sparse. An axisymmetric inverse finite element (FE) approach is presented to represent the biomechanics of vertical atraumatic tooth extraction in an ex-vivo swine model. METHODS: Geometry and boundary conditions from the model are determined to match the extraction of swine incisors in a self-aligning ex vivo extraction experiment. Material parameters for the periodontal ligament (PDL) model are determined by solving an inverse FE problem using clusters of data obtained from 10 highly-controlled mechanical experiments. A seven-parameter visco-hyperelastic damage model, based on an Arruda-Boyce framework, is used for curve fitting. Three loading schemes were fit to obtain a common set of material parameters. RESULTS: The inverse FE results demonstrate good predictions for overall force-time curve shape, peak force, and time to peak force. The fit model parameters are sufficiently consistent across all three cases that a coefficient-averaged model was taken that compares well to all three cases. Notably, the initial modulus ,u, converged across trials to an average value of 0.472 MPa with an average viscoelastic constant g of 0.561. CONCLUSIONS: The presented model is found to have consistent parameters across loading cases. The capability of this model to represent the fundamental mechanical characteristics of the dental complex during vertical extraction loading is a significant advancement in the modelling of extraction procedures. Future work will focus on verifying the model as a predictive design tool for assessing new loading schemes in addition to investigating its applications to subject-specific problems.

Trabecular Bone Component Assessment under Orthodontic Loads and Movements during Periodontal Breakdown-A Finite Elements Analysis.

This numerical analysis, by employing Tresca and Von Mises failure criteria, assessed the biomechanical behavior of a trabecular bone component subjected to 0.6, 1.2, and 2.4 N orthodontic forces under five movements (intrusion, extrusion, tipping, rotation, and translation) and during a gradual horizontal periodontal breakdown (0-8 mm). Additionally, they assessed the changes produced by bone loss, and the ischemic and resorptive risks. The analysis employed eighty-one models of nine patients in 405 simulations. Both failure criteria showed similar qualitative results, with Tresca being quantitatively higher by 1.09-1.21. No qualitative differences were seen between the three orthodontic loads. Quantitatively, a doubling (1.2 N) and quadrupling (2.4 N) were visible when compared to 0.6 N. Rotation and translation followed by tipping are the most stressful, especially for a reduced periodontium, prone to higher ischemic and resorptive risks. In an intact periodontium, 1.2 N can be safely applied but only in a reduced periodontium for extrusion and intrusion. More than 0.6 N is prone to increasing ischemic and resorptive risks for the other three movements. In an intact periodontium, stress spreads in the entire trabecular structure. In a reduced periodontium, stress concentrates (after a 4 mm loss-marker for the stress change distribution) and increases around the cervical third of the remaining alveolar socket.

Nested structure role in the mechanical response of spicule inspired fibers.

Euplectella aspergillummarine sponge spicules are renowned for their remarkable strength and toughness. These spicules exhibit a unique concentric layering structure, which contributes to their exceptional mechanical resistance. In this study, finite element method simulations were used to comprehensively investigate the effect of nested cylindrical structures on the mechanical properties of spicules. This investigation leveraged scanning electron microscopy images to guide the computational modeling of the microstructure and the results were validated by three-point bending tests of 3D-printed spicule-inspired structures. The numerical analyses showed that the nested structure of spicules induces stress and strain jumps on the layer interfaces, reducing the load on critical zones of the fiber and increasing its toughness. It was found that this effect shows a tapering enhancement as the number of layers increases, which combines with a threshold related to the 3D-printing manufacturability to suggest a compromise for optimal performance. A comprehensive evaluation of the mechanical properties of these fibers can assist in developing a new generation of bioinspired structures with practical real-world applications.

Auxetics and FEA: Modern Materials Driven by Modern Simulation Methods.

Auxetics are materials, metamaterials or structures which expand laterally in at least one cross-sectional plane when uniaxially stretched, that is, have a negative Poisson's ratio. Over these last decades, these systems have been studied through various methods, including simulations through finite elements analysis (FEA). This simulation tool is playing an increasingly significant role in the study of materials and structures as a result of the availability of more advanced and user-friendly commercially available software and higher computational power at more reachable costs. This review shows how, in the last three decades, FEA proved to be an essential key tool for studying auxetics, their properties, potential uses and applications. It focuses on the use of FEA in recent years for the design and optimisation of auxetic systems, for the simulation of how they behave when subjected to uniaxial stretching or compression, typically with a focus on identifying the deformation mechanism which leads to auxetic behaviour, and/or, for the simulation of their characteristics and behaviour under different circumstances such as impacts.

Modelling Inductive Sensors for Arc Fault Detection in Aviation.

Modern aircraft are being equipped with high-voltage and direct current (HVDC) architectures to address the increase in electrical power. Unfortunately, the rise of voltage in low pressure environments brings about a problem with unexpected ionisation phenomena such as arcing. Series arcs in HVDC cannot be detected with conventional means, and finding methods to avoid the potentially catastrophic hazards of these events becomes critical to assure further development of more electric and all electric aviation. Inductive sensors are one of the most promising detectors in terms of sensitivity, cost, weight and adaptability to the circuit wiring in aircraft electric systems. In particular, the solutions based on the detection of the high-frequency (HF) pulses created by the arc have been found to be good candidates in practical applications. This paper proposes a method for designing series arc fault inductive sensors able to capture the aforementioned HF pulses. The methodology relies on modelling the parameters of the sensor based on the physics that intervenes in the HF pulses interaction with the sensor itself. To this end, a comparative analysis with different topologies is carried out. For every approach, the key parameters influencing the HF pulses detection are studied theoretically, modelled with a finite elements method and tested in the laboratory in terms of frequency response. The final validation tests were conducted using the prototypes in real cases of detection of DC series arcs.

Poroelastic Characterization and Modeling of Subcutaneous Tissue Under Confined Compression.

Subcutaneous tissue mechanics are important for drug delivery. Yet, even though this material is poroelastic, its mechanical characterization has focused on its hyperelastic response. Moreover, advancement in subcutaneous drug delivery requires effective tissue mimics such as hydrogels for which similar gaps of poroelastic data exist. Porcine subcutaneous samples and gelatin hydrogels were tested under confined compression at physiological conditions and strain rates of 0.01%/s in 5% strain steps with 2600 s of stress relaxation between loading steps. Force-time data were used in an inverse finite element approach to obtain material parameters. Tissues and gels were modeled as porous neo-Hookean materials with properties specified via shear modulus, effective solid volume fraction, initial hydraulic permeability, permeability exponent, and normalized viscous relaxation moduli. The constitutive model was implemented into an isogeometric analysis (IGA) framework to study subcutaneous injection. Subcutaneous tissue exhibited an initial spike in stress due to compression of the solid and fluid pressure buildup, with rapid relaxation explained by fluid drainage, and longer time-scale relaxation explained by viscous dissipation. The inferred parameters aligned with the ranges reported in the literature. Hydraulic permeability, the most important parameter for drug delivery, was in the range k 0 ∈ [ 0.142 , 0.203 ] mm 4 /(N s). With these parameters, IGA simulations showed peak stresses due to a 1-mL injection to reach 48.8 kPa at the site of injection, decaying after drug volume disperses into the tissue. The poro-hyper-viscoelastic neo-Hookean model captures the confined compression response of subcutaneous tissue and gelatin hydrogels. IGA implementation enables predictive simulations of drug delivery.

Effects of different contact angles during forefoot running on the stresses of the foot bones: a finite element simulation study.

Introduction: The purpose of this study was to compare the changes in foot at different sole-ground contact angles during forefoot running. This study tried to help forefoot runners better control and improve their technical movements by comparing different sole-ground contact angles. Methods: A male participant of Chinese ethnicity was enlisted for the present study, with a recorded age of 25 years, a height of 183 cm, and a body weight of 80 kg. This study focused on forefoot strike patterns through FE analysis. Results: It can be seen that the peak von Mises stress of M1-5 (Metatarsal) of a (Contact angle: 9.54) is greater than that of b (Contact angle: 7.58) and c (Contact angle: 5.62) in the three cases. On the contrary, the peak von Mises stress of MC (Medial Cuneiform), IC (Intermediate Cuneiform), LC (Lateral Cuneiform), C (Cuboid), N (Navicular), T (Tarsal) in three different cases is opposite, and the peak von Mises stress of c is greater than that of a and b. The peak von Mises stress of b is between a and c. Conclusion: This study found that a reduced sole-ground contact angle may reduce metatarsal stress fractures. Further, a small sole-ground contact angle may not increase ankle joint injury risk during forefoot running. Hence, given the specialized nature of the running shoes designed for forefoot runners, it is plausible that this study may offer novel insights to guide their athletic pursuits.

Solar fuels design: Porous cathodes modeling for electrochemical carbon dioxide reduction in aqueous electrolytes.

The reduction of carbon dioxide emissions is crucial to reduce the atmospheric greenhouse effect, fighting climate change and global warming. Electrochemical CO2 reduction is one of the most promising carbon capture and utilization technologies, that can be powered by solar energy and used to make added-value chemicals and green fuels, providing grid-stability, energy security, and environmental benefits. A two-dimensional finite-elements model for porous electrodes was developed and validated against experimental data, allowing the design and performance improvement of a porous zinc cathode morphology and its operational conditions for an electrolyzer producing syngas via the co-electrolysis of CO2 and water. Porosity, pore length, fiber geometric shape, inlet pressure, system temperature, and catholyte flow rate were explored, and these parameters were thoroughly tuned by using the smart-search Nelder-Mead's multi-parameter optimization algorithm to achieve pronouncedly higher, industrial-relevant current density values than those previously reported, up to 263.6 mA/cm2 at an applied potential of -1.1 V vs. RHE.

Comparative evaluation of methods for obtaining the stress-free geometry of keratoconus corneas with different levels of severity.

A biomechanical model that simulates the physiological pressure load on the cornea without considering the stresses in vivo will result in an overstress or underestimation of the stress field and incorrect deformation of the corneal structure. Therefore, it is essential to propose numerical models that consider the stress-free geometry of the cornea. In this study, the Displacement and Pre-stress methods were compared to obtain the stress-free geometry (S-FG) and the physiological estimated geometry (P-EG), based on the patient-specific geometric behavior and the computational time required to reach each geometry. The same shape and contour conditions were considered in the models obtained from both methods for each of the pathological scenarios analyzed. Both methods behaved differently to obtain the free geometry, and this difference increased with the severity grade of the disease. However, they behaved in a similar way to reach the physiological estimated geometry. The Displacement method required a lower computational cost to reach the free geometry, with both methods presenting a similar computational cost to obtain the physiological geometry. The stress-free geometries obtained by both methods allowed to characterize the existing biomechanical decompensation during the progression of the diseases. In conclusion, the calculation of the stress-free corneal geometry associated to the clinically measured intraocular pressure with the Displacement and Prestress Methods in keratoconus eyes allows the development of accurate and useable models in clinical practice in real time. This displacement method shows some benefits in terms of computational cost.

Electroporation-Based Biopsy Treatment Planning with Numerical Models and Tissue Phantoms.

Molecular sampling with vacuum-assisted tissue electroporation is a novel, minimally invasive method for molecular profiling of solid lesions. In this paper, we report on the design of the battery-powered pulsed electric field generator and electrode configuration for an electroporation-based molecular sampling device for skin cancer diagnostics. Using numerical models of skin electroporation corroborated by the potato tissue phantom model, we show that the electroporated tissue volume, which is the maximum volume for biomarker sampling, strongly depends on the electrode's geometry, needle electrode skin penetration depths, and the applied pulsed electric field protocol. In addition, using excised human basal cell carcinoma (BCC) tissues, we show that diffusion of proteins out of human BCC tissues into water strongly depends on the strength of the applied electric field and on the time after the field application. The developed numerical simulations, confirmed by experiments in potato tissue phantoms and excised human cancer lesions, provide essential tools for the development of electroporation-based molecular markers sampling devices for personalized skin cancer diagnostics.

Finite element analysis in the optimization of posterior-stabilized total knee arthroplasty.

Posterior-stabilized total knee arthroplasty (PS-TKA) is associated with high rates of satisfaction and functional recovery. This is notably attributed to implant optimization in terms of design, choice of materials, positioning and understanding of biomechanics. Finite elements analysis (FEA) is an assessment technique that contributed to this optimization by ensuring mechanical results based on numerical simulation. By close teamwork between surgeons, researchers and engineers, FEA enabled testing of certain clinical impressions. However, the methodological features of the technique led to wide variations in the presentation and interpretation of results, requiring a certain understanding of numerical and biomechanical fields by the orthopedic community. The present study provides an up-to-date review, aiming to address the following questions: what are the principles of FEA? What is the role of FEA in studying PS design in TKA? What are the key elements in the literature for understanding the role of FEA in PS-TKA? What is the contribution of FEA for understanding of tibiofemoral and patellofemoral biomechanical behavior? What are the limitations and perspectives of digital simulation and FEA in routine practice, with a particular emphasis on the "digital twin" concept? LEVEL OF EVIDENCE: V, expert opinion.

The inclusion of the epithelium in numerical models of the human cornea.

We present a patient-specific finite element model of the human cornea that accounts for the presence of the epithelium. The thin anterior layer that protects the cornea from the external actions has a scant relevance from the mechanical point of view, and it has been neglected in most numerical models of the cornea, which assign to the entire cornea the mechanical properties of the stroma. Yet, modern corneal topographers capture the geometry of the epithelium, which can be naturally included into a patient-specific solid model of the cornea, treated as a multi-layer solid. For numerical applications, the presence of a thin layer on the anterior cornea requires a finer discretization and the definition of two constitutive models (including the corresponding properties) for stroma and epithelium. In this study, we want to assess the relevance of the inclusion of the epithelium in the model of the cornea, by analyzing the effects in terms of uncertainties of the mechanical properties, stress distribution across the thickness, and numerical discretization. We conclude that if the epithelium is modeled as stroma, the material properties should be reduced by 10%. While this choice represents a sufficiently good approximation for the simulation of in vivo mechanical tests, it might result into an under-estimation of the postoperative stress in the simulation of refractive surgery.

Effects of Increasing the Orthodontic Forces over Cortical and Trabecular Bone during Periodontal Breakdown-A Finite Elements Analysis.

Background and Objectives: Herein we used numerical analysis to study different biomechanical behaviors of mandibular bone subjected to 0.6 N, 1.2 N, and 2.4 N orthodontic loads during 0-8 mm periodontal breakdown using the Tresca failure criterion. Additionally, correlations with earlier FEA reports found potential ischemic and resorptive risks. Materials and Methods: Eighty-one models (nine patients) and 243 simulations (intrusion, extrusion, rotation, tipping, and translation) were analyzed. Results: Intrusion and extrusion displayed after 4 mm bone loss showed extended stress display in the apical and middle third alveolar sockets, showing higher ischemic and resorptive risks for 0.6 N. Rotation, translation, and tipping displayed the highest stress amounts, and cervical-third stress with higher ischemic and resorptive risks after 4 mm loss for 0.6 N. Conclusions: Quantitatively, rotation, translation, and tipping are the most stressful movements. All three applied forces produced similar stress-display areas for all movements and bone levels. The stress doubled for 1.2 N and quadrupled for 2.4 N when compared with 0.6 N. The differences between the three loads consisted of the stress amounts displayed in color-coded areas, while their location and extension remained constant. Since the MHP was exceeded, a reduction in the applied force to under 0.6 N (after 4 mm of bone loss) is recommended for reducing ischemic and resorptive risks. The stress-display pattern correlated with horizontal periodontal-breakdown simulations.

Orthodontic Internal Resorption Assessment in Periodontal Breakdown-A Finite Elements Analysis (Part II).

This finite elements analysis (FEA) assessed the accuracy of maximum shear stress criteria (Tresca) in the study of orthodontic internal surface resorption and the absorption-dissipation ability of dental tissues. The present study was conducted over eighty-one models totaling 324 simulations with various bone loss levels (0-8 mm), where 0.6 N and 1.2 N were applied in the intrusion, extrusion, rotation, tipping, and translation movements. Tresca criteria displayed localized high-stress areas prone to resorption for all situations, better visible in the dentine component. The internal resorptive risks are less than external ones, seeming to increase with the progression of the periodontal breakdown, especially after 4 mm. The internal and external surface high-stress areas are strictly correlated. The qualitative stress display for both forces was almost similar. The rotation and tipping displayed the highest resorptive risks for the pulp chamber, decreasing with bone loss. The resorptive risks seem to increase along with the progression of periodontal breakdown if the same applied force is kept. The dentine resemblance to ductile based on its high absorption-dissipation ability seems correct. Tresca seems to supply a better predictability of the prone-to-resorption areas than the other failure criteria.

A simulation-driven supervised learning framework to estimate brain microstructure using diffusion MRI.

We propose a framework to train supervised learning models on synthetic data to estimate brain microstructure parameters using diffusion magnetic resonance imaging (dMRI). Although further validation is necessary, the proposed framework aims to seamlessly incorporate realistic simulations into dMRI microstructure estimation. Synthetic data were generated from over 1,000 neuron meshes converted from digital neuronal reconstructions and linked to their neuroanatomical parameters (such as soma volume and neurite length) using an optimized diffusion MRI simulator that produces intracellular dMRI signals from the solution of the Bloch-Torrey partial differential equation. By combining random subsets of simulated neuron signals with a free diffusion compartment signal, we constructed a synthetic dataset containing dMRI signals and 40 tissue microstructure parameters of 1.45 million artificial brain voxels. To implement supervised learning models we chose multilayer perceptrons (MLPs) and trained them on a subset of the synthetic dataset to estimate some microstructure parameters, namely, the volume fractions of soma, neurites, and the free diffusion compartment, as well as the area fractions of soma and neurites. The trained MLPs perform satisfactorily on the synthetic test sets and give promising in-vivo parameter maps on the MGH Connectome Diffusion Microstructure Dataset (CDMD). Most importantly, the estimated volume fractions showed low dependence on the diffusion time, the diffusion time independence of the estimated parameters being a desired property of quantitative microstructure imaging. The synthetic dataset we generated will be valuable for the validation of models that map between the dMRI signals and microstructure parameters. The surface meshes and microstructures parameters of the aforementioned neurons have been made publicly available.

A Unified Determinant-Preserving Formulation for Compressible/Incompressible Finite Viscoelasticity.

This paper presents a formulation alongside a numerical solution algorithm to describe the mechanical response of bodies made of a large class of viscoelastic materials undergoing arbitrary quasistatic finite deformations. With the objective of having a unified formulation that applies to a wide range of highly compressible, nearly incompressible, and fully incompressible soft organic materials in a numerically tractable manner, the viscoelasticity is described within a Lagrangian setting by a two-potential mixed formulation. In this formulation, the deformation field, a pressure field that ensues from a Legendre transform, and an internal variable of state Fv that describes the viscous part of the deformation are the independent fields. Consistent with the experimental evidence that viscous deformation is a volume-preserving process, the internal variable Fv is required to satisfy the constraint det Fv=1. To solve the resulting initial-boundary-value problem, a numerical solution algorithm is proposed that is based on a finite-element (FE) discretization of space and a finite-difference discretization of time. Specifically, a Variational Multiscale FE method is employed that allows for an arbitrary combination of shape functions for the deformation and pressure fields. To deal with the challenging non-convex constraint det Fv=1, a new time integration scheme is introduced that allows to convert any explicit or implicit scheme of choice into a stable scheme that preserves the constraint det Fv=1 identically. A series of test cases is presented that showcase the capabilities of the proposed formulation.

Computational mechanobiology model evaluating healing of postoperative cavities following breast-conserving surgery.

Breast cancer is the most commonly diagnosed cancer type worldwide. Given high survivorship, increased focus has been placed on long-term treatment outcomes and patient quality of life. While breast-conserving surgery (BCS) is the preferred treatment strategy for early-stage breast cancer, anticipated healing and breast deformation (cosmetic) outcomes weigh heavily on surgeon and patient selection between BCS and more aggressive mastectomy procedures. Unfortunately, surgical outcomes following BCS are difficult to predict, owing to the complexity of the tissue repair process and significant patient-to-patient variability. To overcome this challenge, we developed a predictive computational mechanobiological model that simulates breast healing and deformation following BCS. The coupled biochemical-biomechanical model incorporates multi-scale cell and tissue mechanics, including collagen deposition and remodeling, collagen-dependent cell migration and contractility, and tissue plastic deformation. Available human clinical data evaluating cavity contraction and histopathological data from an experimental porcine lumpectomy study were used for model calibration. The computational model was successfully fit to data by optimizing biochemical and mechanobiological parameters through Gaussian process surrogates. The calibrated model was then applied to define key mechanobiological parameters and relationships influencing healing and breast deformation outcomes. Variability in patient characteristics including cavity-to-breast volume percentage and breast composition were further evaluated to determine effects on cavity contraction and breast cosmetic outcomes, with simulation outcomes aligning well with previously reported human studies. The proposed model has the potential to assist surgeons and their patients in developing and discussing individualized treatment plans that lead to more satisfying post-surgical outcomes and improved quality of life.

Evaluation of a multi-scale discrete fiber model for analyzing arterial failure.

So far, the prevalent rupture risk quantification of aortic aneurysms does not consider information of the underlying microscopic mechanisms. Uniaxial tension tests were performed on imaged aorta samples oriented in circumferential and longitudinal directions. To account for local heterogeneity in collagen fiber architecture, SHG imaging was performed on tissues at several locations prior to mechanical testing. This enabled the quantification of micro-scale information including organization of collagen fibers using relevant probability density functions. Two different modeling approaches are presented in this study for the sake of comparison. A multi-scale mechanical model was developed using this micro-structural information with collagen fibers as main components. accounting for non-affine fiber kinematics. Simultaneously, an embedded element model that accounts for affine fiber kinematics was developed in Abaqus using the same micro-structural information. Numerical simulations emulating uniaxial tension experiments were performed on the developed models. Global mechanical response of both models agreed well with the experimental data, although leading to mismatched material properties. The models present a rudimentary yet better than before representation of structure based description of aortic-tissue failure mechanics. reinforcing the importance of structural organization of micro-scale constituents and their kinematics in determining tissue failure.

FEM Simulation of AlSi10Mg Artifact for Additive Manufacturing Process Calibration with Industrial-Computed Tomography Validation.

Dimensional accuracy of selective laser melting (SLM) parts is one of manufacturers' major concerns. The additive manufacturing (AM) process is characterized by high-temperature gradients, consolidation, and thermal expansion, which induce residual stress on the part. These stresses are released by separating the part from the baseplate, leading to plastic deformation. Thermo-mechanical finite elements (FE) simulation can be adopted to determine the effect of process parameters on final geometrical accuracy and minimize non-compliant parts. In this research, a geometry for process parameter calibration is presented. The part has been manufactured and then analyzed with industrial computed tomography (iCT). An FE process simulation has been performed considering material removal during base plate separation, and the computed distortions have been compared with the results of the iCT, revealing good accordance between the final product and its digital twin.