Numerical Finite Element Method Research Articles

Influence of Subsoil and Building Material Properties on Mine-Induced Soil–Structure Interaction Effect

Soil–structure interaction (SSI) refers to the dynamic interaction between a structure and the surrounding soil on which it rests. The behavior of the soil can significantly affect the response of the building structure. In the context of civil engineering and structural analysis, SSI becomes particularly important when considering the response of structures to dynamic loads such as earthquakes or so-called paraseismic loads, e.g., mining tremors. Several factors contribute to SSI. Soil and building structure material properties, foundation type, and loading conditions are the most important parameters. The article concerns SSI in the case of mining rock bursts in Poland. The influence of changes in site material conditions and building material properties on the SSI phenomenon was investigated. A few variants of different properties of typical construction materials (brick, reinforced concrete, and cellular concrete) in the case of selected representative building structure were considered. The subsoil material properties from the wide range were also taken into account. Numerical three-dimensional finite element method (FEM) analysis was applied. The adopted models of the soil-structure system were verified by data from in situ experimental vibration measurements. A significant influence of the subgrade material and the building structure material on the SSI was demonstrated.

Complex variable solution for ground deformation responses considering the grouting pressure of shallow shield

During shield tunneling in shallow strata, ground deformation is highly sensitive to grouting pressure response. Excessive grouting pressure can cause ground uplift, and an excessively low grouting pressure can cause excessive ground settlement, affecting the surrounding environment. Current studies seldom analyze the influence of shield excavation on ground deformation when the grouting pressure is considered. To ensure micro-disturbance of the surrounding environment during shield tunneling, this study presents a complex variable analysis method for analyzing the influence of grouting pressure on surface settlement. This solution maps the semi-infinite strata outside the segment into circular domains by introducing a complex function conformal mapping method. The nonuniform convergence deformation pattern BC-4 was introduced as the boundary condition for the displacement around the tunnel. The Laurent series expansion of the analytical function considering the grouting pressure in the circular domain was substituted into each boundary condition to solve the analytical function. Subsequently, the accuracy of the analytical method was verified by combining it with a finite element numerical method. Finally, the influence of grouting pressure on ground deformation under different working conditions was analyzed through parameterization. This analytical method can provide a theoretical basis for controlling grouting pressure parameters in shield tunnels.

6σ Robust Lightweight Design of the GIL/GIS Bus Based on Co‐Simulation of Insulation, Current Flow, and Mechanical Stress

GIL/GIS bus is widely used in the field of power transmission and distribution, and the demand for lightweight is becoming more and more prominent. To reduce the weight of the bus and ensure the safety and reliability of the structure, a new 6σ robust design method of the GIL/GIS bus based on electric‐thermal‐mechanical co‐simulation is proposed. Taking a 550 kV/5000A bus as the initial structure, considering the uncertainty of product size and performance, a 6σ robust lightweight model of bus insulation electric field strength, current temperature rise, and mechanical stress multi‐performance co‐simulation is established. The finite element numerical calculation method is used to calculate and analyze the performance parameters of the bus. The Kriging surrogate model is used to form the response surface, and the optimal bus structure is obtained by genetic algorithm. The results show that after the robustness optimization, the performance of the bus meets the design requirements, the weight per meter decreased from 65.11 to 39.12 kg, and the lightweight effect is remarkable. Compared with the traditional deterministic optimization design, the reliability of the structure is improved. The new method of the GIL/GIS bus lightweight design proposed in this study has high feasibility and application prospects, which can provide new ideas for the development of gas‐insulated electrical equipment. © 2024 Institute of Electrical Engineer of Japan and Wiley Periodicals LLC.

Thermal analysis with high accuracy of multi-beam mask fabrication

For a 7 nm technology node and beyond, multi-beam mask fabrication based on charged particles has attracted attention widely and shows great advantages in terms of throughput. However, the heating effect during mask writing is a serious problem and makes deformation error. To address this issue, an accurate analysis of heating with multi-beam writing is necessary. In this study, the thermal effects of electron beams on a mask during writing time (exposure time and nonexposure time) were simulated with a finite element numerical method. The variation in the temperature field with two writing paths (S-shaped and E-shaped) was analyzed. A comparative analysis of the mask’s deformation under different writing paths was conducted. Numerical research shows that the thermal analysis method in this study provides a guide for optimizing the process parameters of mask fabrication.

Numerical analysis of a pressure distribution field and fluid flow velocity vectors near cumulative perforation holes

Relevance. The need to determine a reliable pressure distribution field and fluid filtration velocity vectors inside the perforation hole and the reservoir rock surrounding it. Aim. On the basis of numerical finite element modeling of the fluid flow inside the perforation holes and its filtration in the surrounding reservoir rock, to reveal the patterns of pressure distribution and fluid filtration vectors in the perforation hole, on its walls and in the near-wellbore zone. Objects. Near-wellbore zone of a limestone reservoir of one of the oil fields in the south of the Perm Region, including perforations. Methods. Numerical finite element method for calculating the flow and filtration of liquid in the near-wellbore zone, taking into account the geometry of perforation holes. Results. The paper considers the main relationships used in numerical simulation of fluid flow and filtration in the ANSYS finite element modeling software package. The authors have developed the finite element scheme of the near-wellbore zone, including cumulative perforation holes and taking into account their geometric parameters, as well as the fact that inside the holes the fluid flow is modeled in open space using the Navier–Stokes equations, and in the surrounding reservoir rock based on the filtration equations and Darcy's law. Numerical calculations were carried out, on the basis of which the distribution of pressure, flow velocities and fluid filtration inside the holes and in the near-wellbore zone as a whole was obtained. Calculations were made with varying pressure in the well (or pressure drawdown), as well as for different values of reservoir permeability. The calculation results showed that for the actual drawdown values of 10 MPa and the reservoir permeability of 50 mD, the value of pressure change inside the perforation holes will not exceed 0.01 MPa, i.e. it can be assumed that it practically does not change inside the hole. It is noted that the maximum value of the filtration rate corresponds to the top of the perforation holes and then its value decreases as it approaches the borehole wall. It is concluded that in further modeling of the stress-strain state of the near-wellbore zone, taking into account the holes of cumulative perforation on the surface of the holes, it is permissible to set a constant pressure value equal to the pressure in the well, and not logarithmic or any other distribution of it.

Mathematical Modeling and Econometric Analysis of Composite Fuel Combustion Process

In this article, the process of combustion of composite fuel is studied. Particles of hydrocarbon fuel were taken as the object of research. The analysis of the combustion process is carried out on the basis of mathematical modeling. The created mathematical model was solved by the numerical method of finite elements using the Elcut program, and a computer model was developed. The results obtained on the basis of the numerical method were compared with the experimental results, as well as an econometric analysis and forecast based on the regression model. A comparison of the results obtained by the finite element method using the Elcut program with experimental data shows that the error is on average 5% and indicates the physical validity of the equations and boundary conditions (1)-(7) for the mathematical modeling of the considered tasks and use of the Elcut program for numerical calculations.

SR-M−GAN: A generative model for high-fidelity stress fields prediction of the composite bolted joints

Predicting the stress fields induced by composite bolted joints constitutes a valuable research endeavor aimed at enhancing the structural integrity of equipment in the aerospace and maritime fields. The stress field characteristics of composite bolted joints may undergo changes due to potential mechanical behaviors, such as damage, constitutive relationships and so on, which exhibit complexity and pose challenges for prediction. Nevertheless, modeling the intricate stress field distribution directly from experimental data poses a formidable challenge, necessitating a substantial investment in materials and resources. Traditional finite element numerical methods have the capability to model the distribution of the stress fields. However, to obtain a high-fidelity solution of this distribution, it is imperative to generate a fine-grid mesh and allocate the significant computational resources. The incorporation of computational techniques in coarse −grid would enhance the efficiency of computations, but these approaches compromise the accuracy of the calculations. In this paper, an innovative deep learning framework is introduced. The network named as the SR-M−GAN is drew inspiration from the conditional generative adversarial network (cGAN) used in the image super-resolution. The primary aim of SR-M−GAN is to mitigate the intrinsic trade-off between result fidelity and computational complexity. The model proposed in this paper effectively reconstructs high-fidelity stress field data from calculations on coarse grids and exhibits excellent reconstruction metrics. Simultaneously, the model only incurs additional computational costs during the deployment of the network model and training. When compared to experimental data, the model continues to demonstrate superior performance in both accuracy and computational efficiency. Therefore, this method offers a high-fidelity data augmentation approach for stress analysis in composite bolted joints, laying the groundwork for large-scale analysis of stress field characteristics in composite bolted joints.

Ramalan Kestabilan Cerun yang Diperkukuhkan Menggunakan Cerucuk Menggunakan Model Sistem Inferens Neuro-Fuzzy Adaptif (ANFIS)

Predictive analysis using artificial intelligence (AI) has transformed the landscape of forecasting analysis in various research fields. The advancements in AI modelling algorithms have enhanced decision-making, trend identification, and process optimization. In geotechnical engineering, AI assists in predicting soil behaviour, structural stability, and slope stability. The AI model discussed in this paper is the Adaptive Neuro-Fuzzy Inference System (ANFIS). In this study, the ANFIS model predicts slope stability by examining the Factor of Safety (FOS) value. Slope stability analyses reinforced with continuous bored pile walls generated by the numerical computation of the finite element method (FEM) in two dimensions (2D) and three dimensions (3D) are compared with the predictions of the ANFIS model. The numerical FEM computations employ PLAXIS 2D and PLAXIS 3D software. Meanwhile, the ANFIS model is designed within the MATLAB software platform involving 112 data samples. With six input pile parameters and one output, the finding shows that the ANFIS model can learn complex non-linear data and accurately predict the output. This is supported by the R² values of 0.9771 and 0.9965 from comparing the forecasting output with the 2D and 3D FEM outputs, respectively. Meanwhile, the low RMSE values of 0.0187 and 0.0180 each confirm this.

FIRE DESIGN OF STEEL MEMBER BY COMPONENT‐BASED FINITE ELEMENT METHOD

AbstractThe Component‐Based Finite Element Method (CBFEM) is the quite commonly used method to analyse and design connections of steel structures. It is the combination of the analytical component method and the numerical Finite Element Method (FEM). In this paper, the CBFEM is applied to analyse steel beams and columns at elevated temperature. The main objective of this study is to verify the CBFEM for predicting the resistance of steel members at elevated temperatures.

Seismic Response Analysis of Multi Span Hyperbolic Brick Arch Thin Shell Structure

The research objective is to investigate the multi-span hyperbolic brick arch thin-shell structure located in Changleyuan, Baoji City, Shaanxi Province. The study utilizes on-site dynamic testing and finite element numerical simulation methods to analyze the dynamic characteristics of the thin-shell structure and its seismic response under various seismic waves and directions. The simulation results indicate the following findings: Vertical earthquakes cause more severe damage to structures compared to horizontal earthquakes. Bidirectional seismic excitation has a greater impact on the structure than unidirectional and three-dimensional seismic excitation. The obtained results provide a scientific basis for the seismic protection, reinforcement, and repair of this multi-span double-curved brick arch thin-shell structure.

Reshaping force for deformed casing repairing with hydraulic rolling reshaper and its influencing factors

Hydraulic rolling reshaper is an advanced reshaping tool to solve the problem of casing deformation, which has been widely used in recent years. When it is used for well repair operation, the reshaping force provided by ground devices is generally determined by experience. However, too large reshaping force may destroy the deformed casing, and too small reshaping force may also prolong the construction period and affect the repairing effect. In this paper, based on Hertz contact theory and elastic-plastic theory, combined with the process parameters of shaping, and considering the structural characteristics of the deformed casing and reshaper, we propose a mathematical model for calculating the reshaping force required for repairing deformed casing by hydraulic rolling reshaper. Meanwhile, the finite element model and numerical method of hydraulic rolling reshaper repairing deformed casing are established by using the finite element method, and the reliability of the mathematical model is verified by several examples. On this basis, the control variable method is used to investigate the influence of each parameter on the reshaping force, and the influence degree of each parameter is explored by orthogonal simulation test and Pearson correlation analysis. The research results not only provide an important theoretical basis for the prediction of reshaping force in on-site construction, but also provide a reference for the subsequent improvement of the shaping process. ©2024 The Authors. Publishing services by Elsevier B.V. on behalf of KeAi Communications Co. Ltd. This is an openaccess article under the CCBY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

Preliminary design and evaluation method for the seismic isolation layer of a shield tunnel

Theoretical analyses and numerical simulations, proposed to increase the efficiency of tunnel seismic isolation measures, are too complex for solving practical problems. Herein, a simplified design method for a tunnel isolation layer is developed for preliminary calculations of the required physical parameters based on the surrounding rock conditions and seismic isolation performance. The internal force response of a shield tunnel in homogeneous ground before and after the installation of a seismic isolation layer is determined using the proposed method, considering the coupled effect of the isolation layer and surrounding rock by calculating the equivalent elastic modulus. The effectiveness of the derived method is verified using wave function expansion and finite element numerical simulation methods. The results obtained using the simplified design formulae are consistent with those obtained using the wave function expansion method. Furthermore, the effects of the seismic wave frequency, surrounding rock type, elastic modulus, and thickness of the seismic isolation layer on the seismic isolation effect are examined via parametric analysis. The proposed method could be applied to preliminary design and evaluation of seismic isolation layers for shield tunnels with different surrounding rock types.

High forward thrust metasurface beam-riding sail.

The radiation pressure force and torque on a one-dimensional bi-grating composed of a Si-SiO2 high contrast binary metagrating is analyzed for the purpose of stable beam riding whereupon a high power laser having an expanding Gaussian irradiance distribution propels the grating in outer space, free from gravitational forces. The binary metagrating structure has been simultaneously optimized to afford high forward thrust, and corrective restoring forces and torques in the event of small linear and angular disturbances. We demonstrate that stability may be enhanced at the expense of forward thrust. The validity of our metamaterial findings is reinforced owing to good agreements between finite-difference time-domain and finite element numerical methods. To reduce mass and enhance forward acceleration this laser-driven sail was designed to be free of a stabilizing boom.

A standard calibration method based on a symmetric resistance network matrix for galvanic logging instruments.

An increasing number of measurement electrodes have been designed to satisfy the demand for high-resolution detection using galvanic logging technology in complex formations. The forward modeling response analysis of logging tools has important guiding significance in the design of galvanic logging tools. Based on a three-dimensional finite element numerical simulation method, we established a forward model of galvanic multi-electrodes in a complex formation. We also designed a symmetrical resistance network model of the formation with equivalent resistance between two electrodes. A symmetrical resistance network was derived using the balanced bridge method. The asymmetrical admittance matrix was extended to a symmetrical extended admittance matrix to realize a convenient calculation of the equivalent symmetrical resistance network in complex formations. Verification of the microcolumn-focused logging tool, with nine electrodes in a simulated standard well, and an evaluation of the degree of invasion in an actual oil well indicate that this calibration method can improve the measurement accuracy of galvanic logging instruments.

Investigation of the thermal regime of a linear DC motor

Linear DC motors (LDPT) with a limited armature stroke are used in mechanisms that reciprocate the working body or stabilize its position in tracking systems. The use of LDPT instead of rotary motion motors with a converter in the form of a crank mechanism, a screw pair or other similar mechanical systems leads to a significant simplification of the kinematic scheme, increased efficiency and, as a rule, increased accuracy of the complex and its reliability. There are a number of problems in the creation of LDPT, which are not fully reflected in the scientific and technical literature. In particular, with parametric optimization of the engine according to weight, size and energy criteria, one of the main limiters is overheating of the winding, and in magnetoelectric machines, also overheating of permanent magnets, which are very sensitive to the magnitude of the operating temperature and thermal shocks. The article presents an analytical method for calculating the thermal state of an engine in steady state, based on the method of thermal substitution schemes. The technique is designed for use in a parametric optimization program. To clarify the values of the a priori coefficients in the analytical expressions of the proposed algorithm, the authors modeled the LDPT using the numerical finite element method. Based on the simulation results, a picture of the temperature field of the machine was obtained and the critical heating areas were determined. Verification of the theoretical provisions was carried out during experimental studies of the prototype engine. The tests showed good convergence of the calculated and experimental parameters.

Paleotectonic Stress and Present Geostress Fields and Their Implications for Coalbed Methane Exploitation: A Case Study from Dahebian Block, Liupanshui Coalfield, Guizhou, China

The macroscopic structural fractures (joints) and geostress distribution characteristics of coal reservoirs are important factors affecting the exploitation of coalbed methane (CBM). In this study, the joints in the sedimentary strata of the Dahebian block in Liupanshui area, Guizhou Province were investigated. Directional coal samples were collected for observation and statistical analysis of coal microfractures, the paleotectonic stress fields of the study area were reconstructed, and the tectonic evolution was elucidated. The geostress distribution characteristics of the target coal seam (coal seam No. 11, P3l) in the study area were analyzed using the finite element numerical simulation method. The results indicate that the structural evolution of the Dahebian syncline in the study area can be divided into two stages. The Late Jurassic–Early Cretaceous stage (Early Yanshanian) is the first stage. Affected by the sinistral strike slip of the Weining–Ziyun–Luodian (WZL) fault zone, the derived stress field in the study area exhibits maximum principal stress (σ1) in the NEE–SWW direction. The Late Cretaceous stage (Late Yanshanian) is the second stage. Affected by the dextral strike slip of the WZL fault zone, the derived stress field exhibits σ1 in the NNW–SSE direction. The folds and faults formed in the first stage were modified by the structural deformation in the second stage. The dominant strikes of joints in the sedimentary strata are found to be in the NW–NNW (300°–360°) and NE (30°–60°) directions, with dip angles mostly ranging from 60° to 90°. The dominant strikes of coal microfractures are in the NW (285°–304°) and NE (43°–53°) directions. The distribution of geostress in the study area is characterized by high levels of geostress in the syncline center, decreasing towards the surrounding periphery. The overall trend of the geostress contour line is similar to the shape of the syncline and is influenced by folds and faults. The σ1 of coal seam No. 11 is vertical stress. The prediction results show that the joint density of coal seam No. 11 in the block is 36–50 joints/m, and the shape of the joint density contour line is also affected by the axial direction of the Dahebian syncline and the surrounding faults. The variation in coal seam joint density and the control effect of geostress on joints opening or closing affects the permeability of coal reservoirs. The study results provide significant guidance for the exploitation of CBM.

Numerical Simulation of Impact Effect on Stability of Transmission Tower Foundation

The impact effect can cause structural damage to a transmission tower’s foundation and affect its overall stability. In order to study the influence of the impact effect on the stability of transmission tower foundations, a three-dimensional finite element numerical simulation method was used to investigate the variations in the extent of damage, displacement, and inclination degree of a transmission tower foundation under different impact velocities, impact durations, impactor shapes, and impact locations. The results show that as the impact velocity increases, the damage value of the transmission tower foundation continuously increases, and the damaged area expands. The lateral displacement value increases continuously with the duration of the impact effect, and the variation in lateral displacement follows a linear function distribution. The inclination degree of the transmission tower foundation increases continuously with increased impact duration and can lead to overturning failure. A smaller impact contact area results in a larger compressive damage value for the transmission tower foundation, and different impact contact areas lead to different modes of failure for the transmission tower foundation. The damage value and damaged area of the transmission tower foundation vary with the location of the impact. By comparing the deformation of the transmission tower foundation before and after reinforcement, it is evident that the reinforcement design can significantly improve the deformation resistance and anti-overturning capacity of the transmission tower foundation.

Investigation of Ductile Point-Contact Yield Strength and Bilinear Reinforced Elastic–Plastic Mechanical Behavior Model

Ductile point-contact structures are widely used in engineering applications, with their yield strength and mechanical behavior directly related to the application safety of relevant engineering structures. The Hertz elastic contact theory can effectively express the mechanical behavior of a point-contact structure at the initial stage of load bearing. However, it cannot be applied when the structure has a characteristic yield. In this study, by analyzing the stress distribution characteristics of a point-contact structure prior to characteristic yield, the maximum contact stress on the contact surface was numerically calculated and summarized. Based on the assumption of simplified bilinear constitutive relation of materials, an elastic–plastic analytical model was derived to characterise the mechanical behavior of ductile point-contact structures before and after yield, with its validity verified using the finite element numerical calculations and experimental method.

The Stress Strain Analysis of Forward Extrusion Process of Hollow Elements

This paper presents an analysis focused on the stress-strain dynamics within the forward extrusion process for hollow elements. The forward extrusion of hollow elements is characterized by extensive plastic deformation and complex material flow. This process involves splitting the metal continuum, subjecting it to plastic deformation, and merging it into a seamless product with a predetermined cross-sectional shape as it exits the tool. Employing modern numerical finite element methods is essential for designing such processes, as it enables the simulation of metal behavior within the tool. These simulations help predict mechanical parameters within the deformable material continuum, such as stress, strain, temperatures, and strain rates. Additionally, they allow for the optimization of tool designs tailored to the specific product, improving process efficiency and quality. Understanding the stress-strain distribution at individual nodes enables proactive measures to prevent mechanical disruptions within the continuum, ensuring the integrity and quality of the final product—an essential criterion in the forward extrusion process for hollow elements.

An Efficient Parallel Iterative Solver for Controlled-source Electromagnetic 3-D Adaptive Forward Modeling in General Anisotropic Media

Controlled-source electromagnetic (CSEM) surveys are moving increasingly toward realistic and complicated scenarios where the survey region may contain undulated topography, complex geometries, and electrical anisotropy media. In this paper, the adaptive finite element numerical method is employed to discrete the total electric field equation, which can provide precise electromagnetic (EM) responses even with a coarse initial mesh. The unstructured tetrahedral grid is employed to effectively address arbitrary irregular geometry and mountain terrain. Then, the flexible generalized minimum residual solver (FGMRES), auxiliary-space Maxwell pre-conditioner, and grid division technique were used to solve the large-scale linear system of equations, which can stably solve ill-conditioned problems using fewer computing resources. Finally, we conducted a numerical experiment via our newly proposed forward modeling scheme on a synthetic model with multi-contrast electrical anisotropy. It validated that accurate EM fields could be obtained against the semi-analytic solutions, and this iterative solver has good robustness for various anisotropic media. As a result, we have developed a state-of-the-art 3D CSEM anisotropic forward modeling engine, which can quickly and accurately deal with large-scale and complex geo-models.