A new algorithm for volume mesh refinement on merging geometries: Application to liver and vascularisation

This paper presents an approach to calculate temperature in the electrical machines critical parts when operating in intermittent operations. Focus in this paper is particularly placed on the traction motors that should be run in duty cycles with thousands of operating points during several hours, e.g. in e-mobility and railway applications. Normally, application of numerical methods in duty-cycle calculations is limited due to the long computation time which, in turn, results in a lower accuracy in temperature estimation. Here, numerical models are simplified and used efficiently together with analytical models so that calculation time is reduced significantly. The numerical methods include computational fluid dynamics and finite element models. The reduction in calculation time is from hours to minutes. The developed numerical and analytical methods are verified by measurements on an open self-ventilated traction motor and a good agreement is achieved both in transient and at steady-state temperatures.

Adding magnetic resonance (MR)-derived information to standard transrectal ultrasound (TRUS) images for guiding prostate biopsy is of substantial clinical interest. A tumor visible on MR images can be projected on ultrasound (US) by using MR-US registration. A common approach is to use surface-based registration. The authors hypothesize that biomechanical modeling will better control deformation inside the prostate than a regular nonrigid surface-based registration method. The authors developed a novel method by extending a nonrigid surface-based registration algorithm with biomechanical finite element (FE) modeling to better predict internal deformations of the prostate. Data were collected from ten patients and the MR and TRUS images were rigidly registered to anatomically align prostate orientations. The prostate was manually segmented in both images and corresponding surface meshes were generated. Next, a tetrahedral volume mesh was generated from the MR image. Prostate deformations due to the TRUS probe were simulated using the surface displacements as the boundary condition. A three-dimensional thin-plate spline deformation field was calculated by registering the mesh vertices. The target registration errors (TREs) of 35 reference landmarks determined by surface and volume mesh registrations were compared. The median TRE of a surface-based registration with biomechanical regularization was 2.76 (0.81-7.96) mm. This was significantly different than the median TRE of 3.47 (1.05-7.80) mm for regular surface-based registration without biomechanical regularization. Biomechanical FE modeling has the potential to improve the accuracy of multimodal prostate registration when comparing it to a regular nonrigid surface-based registration algorithm and can help to improve the effectiveness of MR guided TRUS biopsy procedures.

Mesh refinement algorithms in an unstructured solver for multiphase flow simulation using discrete particles

Skull stiffness of children is of great significance to protect brain tissues during traffic accidents. In this paper, the effects of skull stiffness on traumatic brain injury (TBI) are investigated by using a validated 6-year-old child head finite element (FE) model according to Chinese standard GB/T 24550-2009 Results showed that the HIC value, intracranial pressure, Von Mises stress and shear strain of brain tissue increased with the decrease of skull stiffness during the impact of FE model with engine hood. Therefore, the injury risk of the child head decreases with the increase of skull stiffness. And the skull with higher stiffness can provide a better protection for brain tissues during traffic accidents.

In the search for possible causes of unfavorable results after stapes surgery, the study reported here focused on the anterior mallear ligament, since it has been previously reported that partial mallear fixation (PMF) leads to functional failure in 38% of cases of stapes revision surgery. The aims of the study were to identify effective methods for the diagnosis of PMF and experimentally assess the conductive hearing loss that results from PMF. The study included vibration amplitude measurements of the ossicles by laser Doppler interferometry (LDI) in 19 patients and 5 fresh human temporal bone (TB) specimens. Analysis of their dynamic behavior was performed by finite element modeling (FEM). Similar, significant changes of manubrium vibration patterns for PMF were found by FEM calculations, in TB experiments, and in patients. We could identify PMF either before operation, using LDI, or during operation, by manual palpation. In the TB experiments and FEM calculations, the attenuation of the stapes displacement due to an isolated PMF was approximately 10 dB and frequency-dependent. Untreated anterior mallear ligament fixation produced a persistent air-bone gap of approximately 10 dB after stapedioplasty.

A finite element model has been developed in the Deform-3D software package for calculating the thermal field of low-alloy steel coil rolled products during rolling in the finishing group of stands and cooling on the discharge roller table of the 2000 wide-strip hot rolling mill of PJSC Magnitogorsk Iron and Steel Works (PJSC MMK). The proposed model takes into account heat release during plastic deformation, cooling of the strip due to metal contact with air, heat transfer to the work rolls, cooling of the strip by water supplied in the process of interstand cooling and in the accelerated cooling unit on the discharge roller table. At the same time, the features of cooling by both upper and lower collectors are taken into account. To achieve the required coiling temperature, the possibility of differentiated activation of collectors of the upper and lower sections of the accelerated cooling unit was implemented, which was introduced into the finite element model. Cooling by water jets was updated by dividing the contact of the cooled strip with water into two zones - the area of water impact with the sheet surface and the area of film boiling. As a result of the analysis of the thermal state of the sheet, an assessment was made of the nature of its change depending on the characteristics of each stage of thermomechanical treatment (TMT). The use of the presented model in the development of TMT technological modes makes it possible to evaluate their correctness, as well as to provide the possibility of predicting the structure and properties formed over the thickness of rolled steel in the process of controlled rolling with accelerated cooling.P. P. Poletskov, S. V. Denisov, P. A. Stekanov, E. V. Braichev, and P. G. Adishchev, post-graduate student of the MTa-21-2 group, took part in the work.The work was carried out at the Nosov Magnitogorsk State Technical University with the financial support of the Ministry of Education and Science of Russia as part of the implementation of a comprehensive project to create high-tech production (Agreement with the Ministry of Education and Science of Russia No. 075-11-2021-063 dated 06/25/2021).

Endovascular coil embolization is now widely used to treat cerebral aneurysms (CA) as an alternative to surgical clipping. It involves filling the aneurysmal sac with metallic coils to reduce flow, induce clotting, and promote the formation of a coil/thrombus mass which protects the aneurysm wall from hemodynamic forces and prevents rupture. However, a significant number of aneurysms are incompletely coiled leading to aneurysm regrowth and/or recanalization. Computational models of aneurysm coiling may provide important new insights into the effects of intrasaccular coil and thrombus on aneurysm wall stresses. Porcine blood and platinum coils were used to construct an in vitro coil thrombus mass (CTM) for mechanical testing. A uniaxial compression test was performed with whole blood clots and CTM, with coil packing densities (CPDs) of 10%, 20%, and 30% to obtain compressive stress/strain responses. A fourth-order polynomial mechanical response function was fit to the experimentally obtained stress/strain responses for each CPD in order to represent their mechanical properties for computational simulations. Patient-specific three-dimensional (3D) geometries of three aneurysms with simple geometry and four with complex geometry were reconstructed from digital subtraction angiography (DSA) images. The CPDs were digitally inserted in the aneurysm geometries and finite element modeling was used to determine transmural peak/mean wall stress (MWS) with and without coil packing. Reproducible stress/strain curves were obtained from compression testing of CTM and the polynomial mechanical response function was found to approximate the experimental stress/strain relationship obtained from mechanical testing to a high degree. An exponential increase in the CTM stiffness was observed with increasing CPD. Elevated wall stresses were found throughout the aneurysm dome, neck, and parent artery in simulations of the CAs with no filling. Complete, 100% filling of the aneurysms with whole blood clot and CPDs of 10%, 20%, and 30% significantly reduced MWS in simple and complex geometry aneurysms. Sequential increases in CPD resulted in significantly greater increases in MWS in simple but not complex geometry aneurysms. This study utilizes finite element analysis to demonstrate the reduction of transmural wall stress following coil embolization in patient-specific computational models of CAs. Our results provide a quantitative measure of the degree to which CPD impacts wall stress and suggest that complex aneurysmal geometries may be more resistant to coil embolization treatment. The computational modeling employed in this study serves as a first step in developing a tool to evaluate the patient-specific efficacy of coil embolization in treating CAs.

A detailed understanding of the changes in load transfer due to implantation is necessary to identify potential failure mechanisms of orthopedic implants. Computational finite element (FE) models provide full field data on intact and implanted bone structures, but their validity must be assessed for clinical relevance. The aim of this study was to test the validity of FE predicted strain distributions for the intact and implanted pelvis using the digital image correlation (DIC) strain measurement technique. FE models of an in vitro hemipelvis test setup were produced, both intact and implanted with an acetabular cup. Strain predictions were compared to DIC and strain rosette measurements. Regression analysis indicated a strong linear relationship between the measured and predicted strains, with a high correlation coefficient (R = 0.956 intact, 0.938 implanted) and a low standard error of the estimate (SE = 69.53 με, 75.09 με). Moreover, close agreement between the strain rosette and DIC measurements improved confidence in the validity of the DIC technique. The FE model therefore was supported as a valid predictor of the measured strain distribution in the intact and implanted composite pelvis models, confirming its suitability for further computational investigations.

Previous attempts by researchers to predict the fatigue behavior of bone cement have been capable of predicting the location of final failure in complex geometries but incapable of predicting cement fatigue life to the right order of magnitude of loading cycles. This has been attributed to a failure to model the internal defects present in bone cement and their associated stress singularities. In this study, dog-bone-shaped specimens of bone cement were micro-computed-tomography (microCT) scanned to generate computational finite element (FE) models before uniaxial tensile fatigue testing. Acoustic emission (AE) monitoring was used to locate damage events in real time during tensile fatigue tests and to facilitate a comparison with the damage predicted in FE simulations of the same tests. By tracking both acoustic emissions and predicted damage back to microCT scans, barium sulfate (BaSO(4)) agglomerates were found not to be significant in determining fatigue life (p=0.0604) of specimens. Both the experimental and numerical studies showed that diffuse damage occurred throughout the gauge length. A good linear correlation (R(2)=0.70, p=0.0252) was found between the experimental and the predicted tensile fatigue life. Although the FE models were not always able to predict the correct failure location, damage was predicted in simulations at areas identified as experiencing damage using AE monitoring.

ObjectivesThe objective of this study is to demonstrate that computational finite element models can be used to reliably simulate dynamic interaction between a pacifier, the palate, and the tongue during nonnutritive sucking (NNS). The interactions can be quantified by the results of finite element analyses which include deformation, strain, stress, contact force, and contact area.Materials and MethodsA finite element model was created based upon CAD solid models of an infant pacifier and palate. The silicone pacifier bulb is represented by a hyperelastic constitutive law. Contact surfaces are defined between the pacifier and palate. A time and spatially varying pressure load is applied to the bulb representing peristaltic interaction with the tongue. A second time‐varying, periodic pressure representing NNS is applied to the model simultaneously. A large displacement, nonlinear transient dynamic analysis is run over two NNS cycles.ResultsResults from the finite element analysis show the deformed shape of the bulb with maximum principal elastic strain of 0.23 and a range of maximum principal stress on the palate from 0.60 MPa (tensile) to −0.27 MPa (compressive) over the NNS cycles. The areas of contact between the pacifier and the palate are shown in surface contour plots.ConclusionsA nonlinear transient dynamic finite element model can simulate the mechanical behavior of a pacifier and its interaction with the tongue and contact with the palate subject to NNS. Quantitative results predicting deformation, strain, stress, contact force, and contact area can be used in comparative studies to provide insight on how pacifiers cause changes in dental, orthognathic, and facial development.

In this paper, energy efficient power control for the uplink two-tier networks where a macrocell tier with a massive multiple-input multiple-output (MIMO) base station is overlaid with a small cell tier is investigated. We propose a distributed energy efficient power control algorithm which allows each user in the two-tier network taking individual decisions to optimize its own energy efficiency (EE) for the multi-user and multi-cell scenario. The distributed power control algorithm is implemented by decoupling the EE optimization problem into two steps. In the first step, we propose to assign the users on the same resource into the same group and each group can optimize its own EE, respectively. In the second step, multiple power control games based on evolutionary game theory (EGT) are formulated for each group, which allows each user optimizing its own EE. In the EGT-based power control games, each player selects a strategy giving a higher payoff than the average payoff, which can improve the fairness among the users. The proposed algorithm has a linear complexity with respect to the number of subcarriers and the number of cells in comparison with the brute force approach which has an exponential complexity. Simulation results show the remarkable improvements in terms of fairness by using the proposed algorithm.

Comparison of rubbing induced vibration responses using varying-thickness-twisted shell and solid-element blade models

Development and evaluation of a finite element model of the THOR for occupant protection of spaceflight crewmembers

An approach for patient-specific multi-domain vascular mesh generation featuring spatially varying wall thickness modeling

More validity of finite element analysis in implant biomechanics requires element downsizing. However, excess downsizing needs computer memory and calculation time. To investigate the effectiveness of element downsizing on the construction of a three-dimensional finite element bone trabeculae model, with different element sizes (600, 300, 150 and 75 microm) models were constructed and stress induced by vertical 10 N loading was analysed. The difference in von Mises stress values between the models with 600 and 300 microm element sizes was larger than that between 300 and 150 microm. On the other hand, no clear difference of stress values was detected among the models with 300, 150 and 75 microm element sizes. Downsizing of elements from 600 to 300 microm is suggested to be effective in the construction of a three-dimensional finite element bone trabeculae model for possible saving of computer memory and calculation time in the laboratory.