Finite Element Calculations Research Articles

Influence of natural hyoid bone position and surgical repositioning on upper airway patency: A computational finite element modeling study

The hyoid bone's inferior baseline position in obstructive sleep apnea (OSA) has led to surgical hyoid repositioning (SHR) treatment, yet outcomes vary widely. The influence of baseline hyoid position (BHP; phenotype) and SHR on upper airway (UA) function remains unclear. We aimed to investigate their impact on the UA using computational modeling. Methods: A validated finite element model of the rabbit UA was advanced and used to simulate changes in BHP and SHR, alone and in combination. The hyoid was displaced in cranial, caudal, anterior, anterior-cranial and anterior-caudal directions from 1-4mm. Model outcomes included UA collapsibility, measured using closing pressure (Pclose), cross-sectional area (CSA) and soft tissue mechanics (displacement, stress and strain). Results: Graded BHP increments increased Pclose for all directions, and up to 29-43% at 4mm (relative to the original BHP). Anterior-based SHR decreased Pclose (~-115% at 4mm) and increased ΔCSA (~+35% at 4mm). Cranial SHR decreased ΔPclose (-29%), minimally affecting CSA. Caudal SHR increased ΔPclose (+27%) and decreased ΔCSA (-7%). Anterior-cranial and anterior-caudal SHR produced the highest stresses and strains. SHR effects on UA outcomes were dependent on BHP, with more caudal BHPs leading to less effective surgeries. Conclusion: BHP (phenotype) and SHR both alter UA outcomes, with effects dependent on hyoid displacement direction and magnitude. BHP influences the effectiveness of SHR in reducing UA collapsibility. These findings provide further insights into the hyoid's role in UA patency and suggest that considering the hyoid's baseline position and surgical repositioning direction/increment may help improve hyoid surgeries for OSA treatment.

Integral equation based wellbore preconditioner for 3D electromagnetic response modeling

Modeling subsurface controlled-source electromagnetic (CSEM) responses using the finite-element (FE) method is challenging in the presence of highly conductive wellbore casing. The very large conductivity contrast between the casing and the host formation leads to increased computation time and potentially unstable solutions. We address this difficulty by preconditioning an FE solver with an integral equation (IE) primary solution that captures the CSEM response of a realistic-sized steel wellbore casing. Our hybrid IE-FE approach determines the primary field using 2D integral-equation forward modeling and then interpolates the IE-computed solution onto the nodes. Then using an existing FE simulator, we solve for the secondary electric and magnetic fields. This approach removes the need for an ultra-fine FE mesh around the wellbore, thereby improving FE solution stability while greatly reducing FE computation time. Our method is illustrated by modeling the CSEM responses of idealized fluid-bearing zones.

Optimization of dissipative mufflers

The method of selecting the configuration of dissipative mufflers with the required acoustic efficiency is considered. The peculiarity of the considered approach is the use of an integral indicator of acoustic efficiency and dimensionless geometric parameters. The studies were carried out using finite element calculations. In the finite element model of a dissipative mufflers, acoustic characteristics of a fibrous sound-absorbing material obtained from experimental studies were used.

Determination of the moment–rotation relation of continuous timber–concrete composite floors based on the component method

Timber–concrete composite (TCC) combines the advantages of timber and reinforced concrete construction. In multi-storey buildings, the design of TCC floors is often governed by the deformation serviceability requirements. One way to improve deformation serviceability is through the design of continuous floor systems. However, the current practice of TCC slab systems is often limited to load-carrying as single-span beams. This paper presents an analytical model based on the component method for determining the moment–rotation relation of continuous TCC slabs in the range with negative bending moments. The investigations involve the development of load–displacement curves for the individual force-transmitting components, including timber and reinforced concrete, followed by their assembly into a component model. The model has been verified using finite element calculations and the results demonstrate that the moment–rotation relation of continuous TCC systems can be accurately captured. As the variability of the material properties strongly influences the stiffness of the joint and therefore the deformation and stresses in the TCC slab, a probabilistic analysis of the material parameters was performed using the Monte Carlo method. The probabilistic investigations show that the variability of the input parameters has a significant impact on the joint stiffness, with notable scattering observed in the moment–rotation relation, especially in the range of the cracking phase of the concrete. The results from the probabilistic investigations enable initial statements about the joint stiffness of continuous timber–concrete composite slabs in ranges with negative bending moments.

Hierarchical meta-porous materials as sound absorbers

The absorption of sound has great significance in many scientific and engineering applications, from room acoustics to noise mitigation. In this context, porous materials have emerged as a viable solution towards high absorption performance and lightweight designs. However, their performance is somewhat limited in the low-frequency regime. Inspired by the concept of recursive patterns over multiple length scales, typical of many natural materials, here, we propose a hierarchical organization of multilayered porous media and investigate their performance in terms of sound absorption. Two types of designs are considered: a hierarchical periodic (HP) and a hierarchical gradient (HG). In both cases, it is found that in some frequency ranges the introduction of multiple levels of hierarchy simultaneously allows for: (i) an increase in the level of absorption compared with the corresponding bulk block of porous material (or to the previous hierarchical levels (HLs)), and (ii) a reduction in the quantity of porous material required. Another advantage of using this approach is on the fabrication, since only one porous material is required. The performances are examined for both normal and oblique incidences, as well as for different values of the static air-flow resistivity. The methodological approach is based on the transfer-matrix method, optimization algorithms such as the metaheuristic greedy randomized adaptive search procedure (GRASP), finite-element calculations and measurements performed in an impedance tube.

Finite element analysis method for the sliding cable analysis: The sliding variable method

A sliding cable structure is typically characterized by uncertainty of sliding and complexity of friction at the cable-strut joints, which makes accurate and efficient analysis difficult. In this study, we have proposed the sliding variable method, to simulate the cable sliding process considering friction equivalently and analyze the internal force and deformation response at the equilibrium state quickly and accurately. The sliding variable method takes the sliding vector of the cable at cable-strut joints as the basic variable, extracts the slippage stiffness matrix and establishes the cable force balance equation to solve the sliding vector, and simulates the change of cable original length and cable force by applying equivalent temperature load. In this paper, first, the basic principles and calculation procedures of the sliding variable method without considering friction, are discussed in detail to simulate frictionless slip of continuous cable, using with roller-type cable-stayed joints. as an example. On this basis, a calculation model incorporating the effects of friction is analyzed, a method to judge the sliding situation and balance relation of cable-forces is revealed, and the calculation principle and steps of the method are further deduced. Meanwhile, choosing a multi-joint cable–pulleys model as an example, the accuracy and rationality of the sliding variable method with and without considering friction are verified respectively. This method requires a conventional finite element calculation and second calculation alternately. Finally, to improve the practicability and computational efficiency of the sliding variable method, a macro document was compiled using the secondary development function of ANSYS parametric programming language (APDL) according to the basic framework of the sliding variable method. The frictional sliding behaviors in two typical engineering cases, a cable girder and a prestressed fish-belly girder on a foundation pit, are analyzed by using the compiled macro document in finite element program. The results show that the compiled macro document can be both accurate and efficient in calculating and analyzing cable sliding problems and can be easily popularized and applied in engineering practices.

Investigation of higher-order springing of a ship in regular waves by experimental analysis and two-way CFD-FEA coupled method

In this paper, the higher-order springing phenomenon is addressed for a segmented barge ship model through experimental and numerical measures. An efficient in-house two-way coupled numerical solver between Computational Fluid Dynamics (CFD) and Finite Element Analysis (FEA) is developed and validated against experimental results. The coupling method is based on a domain-separated approach and necessitates the resolution of individual boundary value problems in distinct domains. To ensure convergence within these individual domains, an implicit numerical scheme is employed and further facilitated exchange of variables for coupling. The current approach emphasizes the overall convergence between two solvers, maintaining a strongly coupled setup to comprehensively address fluid-structure interaction phenomena, including added mass and damping effects. A series of tank tests were conducted first to measure the wave-induced sectional loads and motions, during which the springing responses to very high-order harmonics of wave load were observed. By comparing the numerical prediction with the tank test results for rigid body motion and flexible vertical bending moment (VBM), the proposed numerical method demonstrated agreement with experimental results, affirming its validity and robustness. Finally, the springing response up to 14th order harmonics is discussed and investigated.

Impact resistance study of lime slag–stabilised steel slag mixes

ABSTRACT Given the transient nature of impact loading, exploring the microscopic damage characteristics and damage mechanisms of materials through indoor tests is difficult. In this study, the damage morphology and energy dissipation characteristics of lime slag–stabilised steel slag (LSSS) mixes under vertical impact loading is explored using ABAQUS software, a finite element model of LSSS for falling hammer impact test is established, the validity of finite element calculation is verified by comparing the results of indoor tests, and the strain characteristics, deflection changes, and energy transformation process of LSSS under vertical impact loading are analysed. Results show that the LSSS tensile and compressive strains under impact loading are roughly divided into low strain rate, precracking, and cracking phases with a rapid increase in strain. Under the same impact energy conditions, the LSSS bottom tensile strain, span deflection, drop hammer kinetic energy conversion time, and LSSS plastic energy consumption curve rising slope are sensitive to drop hammer mass. Thus, improvement in drop hammer mass will prolong the impact loading of the drop hammer–LSSS system internal energy conversion time, and the LSSS strain peak range to the central concentration, enhancing span deflection and bottom deformation.

Enhancing solar energy systems: The role of computational and automation techniques in achieving net zero by 2050

This review aims to show how computational and automation can be applied to optimize the solar power system toward net-zero emissions in 2050. It emphasizes the power of data analytics, machine learning and automated systems for optimizing business models and measuring the performance of solar technology. The role of solar power in decarbonizing our energy needs and lowering carbon emissions is discussed along the global climate agenda. It also discusses novel uses such as self-teaching robots for solar panel inspection and automated manufacturing to increase the efficiency and sustainability of solar power generation. Additionally, it addresses the need for sophisticated research tools such as computational fluid dynamics (CFD) and finite element analysis (FEA) to create new technologies in the solar market and ultimately meet the ambitious renewable energy targets of 2050.

Tunable magnonic crystal in a hybrid superconductor–ferrimagnet nanostructure

One of the most intriguing properties of magnonic systems is their reconfigurability, where an external magnetic field alters the static magnetic configuration to influence magnetization dynamics. In this paper, we present an alternative approach to tunable magnonic systems. We studied theoretically and numerically a magnonic crystal induced within a uniform magnetic layer by a periodic magnetic field pattern created by the sequence of superconducting strips. We showed that the spin-wave spectrum can be tuned by the inhomogeneous stray field of the superconductor in response to a small uniform external magnetic field. Additionally, we demonstrated that modifying the width of superconducting strips and separation between them leads to the changes in the internal field which are unprecedented in conventional magnonic structures. The paper presents the results of semi-analytical calculations for realistic structures, which are verified by finite-element method computations.

Simulation study on parameter optimization of transcranial direct current stimulation based on rat brain slices

Transcranial direct current stimulation (tDCS) is an important method for treating mental illnesses and neurodegenerative diseases. This paper reconstructed two ex vivo brain slice models based on rat brain slice staining images and magnetic resonance imaging (MRI) data respectively, and the current densities of hippocampus after cortical tDCS were obtained through finite element calculation. Subsequently, a neuron model was used to calculate the response of rat hippocampal pyramidal neuron under these current densities, and the neuronal responses of the two models under different stimulation parameters were compared. The results show that a minimum stimulation voltage of 17 V can excite hippocampal pyramidal neuron in the model based on brain slice staining images, while 24 V is required in the MRI-based model. The results indicate that the model based on brain slice staining images has advantages in precision and electric field propagation simulation, and its results are closer to real measurements, which can provide guidance for the selection of tDCS parameters and scientific basis for precise stimulation.

Damage analysis and assessment of concrete T-girder bridge based on fire scene numerical reconstruction

Fire is a sporadic disaster of concrete bridges, with diverse fire environments and complex damage mechanisms. Accurately evaluating the damage situation of concrete bridges after a fire is exceedingly challenging. This study formulates a damage analysis and assessment method based on the step-by-step and progressively deepening working principle. The method relies on fire scene numerical reconstruction and encompasses key technical aspects, including bridge detection and analysis during the fire incident, fire scene numerical reconstruction, and subsequent bridge damage assessment. Building upon these principles, the study utilizes results from the detection and analysis of the concrete T-girder bridge during a fire incident to establish Computational Fluid Dynamics (CFD) and Finite Element Method (FEM) models for the numerical reconstruction of the fire scene. These models enable the calculation of varying temperature distributions and the evolution of the bridge under fire. Compared with the parameters obtained through the ISO834 method, the numerical reconstruction approaches not only enhances the accuracy of replicating the bridge combustion process but also enables the extraction of temperature field distribution patterns within the bridge fire space and its concrete components.

Simulation analysis of force and fatigue life of circular wheel of crawler vehicle

During loading and driving, the wheels bear the vertical load from the body mass and the excitation load generated by uneven road surface on the one hand, and bear the driving torque on the other hand. The load-bearing methods of wheels are generally divided into bottom load-bearing and top load-bearing. This paper describes the structural characteristics of the track system of the articulated track vehicle and the interaction relationship between the main components of the track system. Finite element calculation is carried out based on ANSYS software to obtain the stress distribution of each key component under various loading methods. It can be seen from the results that all key components of the track system can meet the strength and rigidity requirements; although there are also areas with large local stress, they are all within the safe range, which is mainly caused by stress concentration. Safe life is obtained through fatigue analysis.

Optimizing Transformation-Induced Plasticity to Resist Microvoid Softening

AbstractMany high-performance steels that are critical for energy-efficient, lightweight designs rely on transformation-induced plasticity (TRIP) to achieve superior combinations of strength and ductility/toughness. Further development of these alloys will require greater optimization of the metastable (retained) austenite phase responsible for TRIP. Considering the complex nature of TRIP and its effects on ductile fracture, an integrated computational materials engineering (ICME) approach to materials optimization is desired. In this work, we report the results of a large series of micromechanical finite element calculations that probe the interaction of TRIP and void-mediated ductile fracture mechanisms. The simulations identify the optimal austenite stability for maximizing the benefit of TRIP across a wide range of stress states. The applied stress triaxiality significantly influences the microvoid growth rate and the computationally determined optimal stability. The simulation results are compared with existing experimental data, demonstrating good agreement.

On the role of the retained porosity on the shock response of additively manufactured high-performance steel: Experiments, constitutive model and finite-element predictions

Experiments have shown that for quasi-static and moderate strain-rates (of the order of 102–103/s) the mechanical response of additively manufactured (AM) and traditionally processed high-strength steels is similar whereas the impact behavior is markedly different. In this paper, we reveal that the main reason for this difference is the retained porosity in the AM material. Fully-implicit finite element calculations are presented in which we simulate both the launching of the impact plate and the impact between the two plates. The constitutive model used is the elastic/plastic model for porous ductile materials with matrix displaying tension-compression asymmetry and Johnson-Cook hardening law that accounts for both strain-rate effects and plastic history. It is shown that even a very small initial porosity changes the wave front, decreases the Hugoniot while increasing the shock rise time, when compared to a void free material. Furthermore, quantitative comparisons between simulation results and plate impact data for both the AM and the wrought AF9628 steel are provided. The good agreement show that the model captures the impact response and illustrates the model capabilities to provide information on field variables that cannot be directly measured.Additive manufacturing (AM) of metals is rapidly advancing as a robust method for production of geometrically complex parts. To enhance understanding of material performance and open up additional application opportunities, dynamic characterization of newly printed alloys is required to validate their effectiveness. In this paper, we present results from plate impact testing of AF9628 steel, a newly developed high-strength low alloy martensitic steel for structural applications which require resistance to high-rate deformation. We put into evidence differences in the shock structure between the AM and the traditionally processed material. To gain understanding, we conduct fully-implicit finite element (FE) calculations in which we model both the launching of the impact plate and the impact between the two plates, respectively. An elastic/plastic damage model that accounts for the effects of the tension-compression asymmetry in plastic deformation and its influence on porosity evolution is used. The FE results reveal that even a very small amount of initial porosity leads to an increase in the shock rise time, explaining the observed trends. Furthermore, quantitative comparisons between simulation results and plate impact data for both the AM and the wrought AF9628 are provided. The good agreement show that the model captures the impact response and illustrates the model capabilities to provide information on field variables that cannot be directly measured.

A review of human cornea finite element modeling: geometry modeling, constitutive modeling, and outlooks.

The cornea is a vital tissue of the human body. The health status of the cornea has a great impact on the quality life of person. There has been a great deal of research on the human cornea biomechancis. However, the difficulty in obtaining the human cornea has greatly limited the research of cornea biomechancis. Using finite element modelling has become a very effective and economical means for studying mechanical properties of human cornea. In this review, the geometrical and constitutive models of the cornea are summarised and analysed, respectively. Some factors affecting of the finite element calculation are discussed. In addition, prospects and challenges for the finite element model of the human cornea are presented. This review will be helpful to researchers performing studies in the relevant fields of human cornea finite element analysis.

Strain-gradient and damage failure behavior in particle reinforced heterogeneous matrix composites

A strain gradient plasticity model with strengthening effect and damage effect (SGSED) is developed under the continuum medium framework for particle-reinforced heterogeneous matrix composites (PRHMCs). Boron carbide (B4C) particle-reinforced ultrafine-grained (UFG)/fine-grained (FG) heterogeneous matrix composites (B4Cp/(UFG/FG)) are fabricated, in which the UFG region consisted of carbon nanotubes (CNTs)/6061 aluminum (Al) flakes with grains in the ultrafine range and the FG region is processed via 6061Al. The SGSED model is written into the user subroutines using commercial finite element (FE) calculation software, and the three-dimensional (3D) FE representative volume element (RVE) for B4Cp/(UFG/FG) composites is established, from which the distribution of the interface-affected-zone (IAZ) formed of the strain gradient caused by the uncoordinated deformation of the UFG-FG heterogeneous matrix and reinforced phase-matrix is calculated. The evolution of the strain gradient in the deformation process of composites and the influence of the strain gradient on the progressive damage and crack evolution of composites are analyzed, and the strain gradient strengthening-toughening mechanism of composites is revealed. It is found that the IAZ has a considerable strengthening-toughening effect on the composites, which can reduce stress concentration at the interface between the reinforced phase and the matrix, and slow down the crack propagation of the matrix.

Finite Element Analysis and Computational Fluid Dynamics for the Flow Control of a Non-Return Multi-Door Reflux Valve

This paper presents a comprehensive analysis of a multi-door check valve using computational fluid dynamics (CFD) and finite element analysis (FEA) to evaluate flow performance under pressure test conditions, with an emphasis on its ability to prevent backflow. Check valves are essential components in various industries, ensuring fluid flow in one direction only while preventing reverse flow. The non-return multi-door reflux valve is increasingly preferred due to its superior backflow prevention, fluid control, and effective flow regulation. Rigorous testing under varying pressure conditions is essential to ensure that these valves perform optimally. This study uses CFD and FEA simulations to evaluate the structural integrity and flow characteristics of the valve, including pressure drop, flow velocity, backflow prevention effectiveness, and flow coefficient. A high-fidelity 3D model was created to simulate the valve’s behavior under various conditions, analyzing the effects of parameters such as the number of doors, their orientation, geometry, and operating conditions. The CFD results demonstrated a significant reduction in backflow and pressure drop across the valve. However, localized turbulence and flow separation near the valve doors, particularly under partially open conditions, have raised concerns about potential wear. The velocity profiles indicated a uniform distribution at full opening with laminar velocity profiles and minimal resistance to flow. The results of the FEA showed that the stresses induced by the fluid forces were below critical levels, with the highest stress concentrations observed around the hinge points of the valve doors. Although the valve structure remained intact under normal operating conditions, some areas may have required reinforcement to ensure long-term durability. Combined CFD and FEA analyses demonstrated that the valve effectively preserves system integrity, prevents backflow, and maintains consistent performance under various pressure and flow conditions. These findings provide valuable insights into design improvements, performance optimization, and enhancing the efficiency and reliability of reflux valve systems in industrial applications.

Dynamic reliability assessment for high CFRDs incorporating ground motion and material parameter variability via GPDEM

The randomness of earthquake and material parameters significantly impacts the seismic response and dynamic reliability of geotechnical engineering structures. In this study, a dynamic reliability analysis method for high concrete face rockfill dams (CFRDs) under the influence of random factors is proposed by combining random sample generation methods with the generalized probability density evolution method (GPDEM). Firstly, the spectral representation-random function method is employed to achieve dimensionality reduction expression for intensity-frequency non-stationary ground motion. Sensitivity analysis is then used to identify the random variables of the rockfill constitutive model, and the improved generalized F-discrepancy (GF-discrepancy) method is utilized to generate random material parameters. Then, the random samples of earthquake and material parameter are introduced into the finite element calculation of CFRD, conducting a series of dynamic finite element calculations. Finally, GPDEM is used to derive probability information and dynamic reliability from the deterministic time history responses. The results indicate that random factors contribute to the variability of seismic response and dynamic reliability. Stochastic ground motion significantly affects extreme values, while variability in material parameters has a more pronounced influence on the overall process of seismic response. These findings provide valuable references for the seismic design of CFRDs.

Development of a new soil–structure contact stress sensor for underground construction applications

This paper describes the design, development, calibration and validation of a novel soil–structure contact stress sensor. The new sensor design combines a novel operating principle, fibre Bragg grating (FBG) strain sensing and data-driven mapping techniques to create a multi-axis contact stress sensor that is both economical and suitably robust for deployment in underground construction applications. The instrumentation process is informed using a ‘virtual twin’ of the sensor in which synthetic data is generated by extracting and interpolating virtual FBG strains obtained from a large number of 3D finite-element calculations. A physical prototype is subsequently developed to demonstrate proof of concept. Results from laboratory validation tests give confidence in the sensor's ability to provide accurate contact stress measurements in typical soil–structure interface shear applications. In particular, the novel sensor structure and operating principle was shown to achieve excellent measurement of effective normal stress. The new sensor design harnesses many of the inherent benefits of FBGs including immunity to electromagnetic noise and water ingress, and the use a single lightweight cable and connector, which significantly simplifies installation on site compared to electrical multi-axis sensors.