Finite Element Numerical Simulation Research Articles

Research on Mechanism of Non-Uniform In-Situ Stress Induced Casing Damage Based on Finite Element Analysis

Casing damage is a common problem encountered during oil and gas field development due to the complex stress state of the casing. Despite the large number of studies focusing on this problem, the mechanism of non-uniform in-situ stress-induced casing damage in a low-permeability reservoir is still unclear. In this study, casing damage due to non-uniform in-situ stress variations was investigated, and then the tectonic stress coefficients in the study area were determined by an in-situ stress inversion technique, which led to the derivation of formulas for calculating the maximum and minimum horizontal in-situ stresses. Subsequently, finite element numerical simulations were performed to assess the stress distribution during the formation of the casing cement sheath in a G155 block, a typical low-permeability reservoir. The results indicate that casing damage is caused not only by non-uniform in-situ stresses but also by various additional creep-induced loads. Subsequent finite element investigations into casing behavior under mudstone creep conditions indicated that immersion of mudstone in water instigated further shearing and deformation of the casing, culminating in premature well failure prior to water inundation. Notably, Von Mises stress levels exhibited a positive correlation with injection production ratios, with values exceeding critical thresholds leading to distinct modes of mechanical failure including shear-induced deformations, longitudinal tensile stress, and localized yielding near water wells. Maintenance of an optimal injection production ratio is identified as a key strategy for prolonging casing longevity in the region. To this end, recommendations include augmenting the casing wall thickness or enhancing the steel pressure specifications to mitigate casing damage progression, thereby extending the operational lifespan. This research serves as a pivotal theoretical framework for informing future development strategies aimed at mitigating and preempting casing failures in a low-permeability reservoir.

Research on Mechanism of Non-Uniform In-Situ Stress Induced Casing Damage Based on Finite Element Analysis

Casing damage is a common problem encountered during oil and gas field development due to the complex stress state of the casing. Despite the large number of studies focusing on this problem, the mechanism of non-uniform in-situ stress-induced casing damage in a low-permeability reservoir is still unclear. In this study, casing damage due to non-uniform in-situ stress variations was investigated, and then the tectonic stress coefficients in the study area were determined by an in-situ stress inversion technique, which led to the derivation of formulas for calculating the maximum and minimum horizontal in-situ stresses. Subsequently, finite element numerical simulations were performed to assess the stress distribution during the formation of the casing cement sheath in a G155 block, a typical low-permeability reservoir. The results indicate that casing damage is caused not only by non-uniform in-situ stresses but also by various additional creep-induced loads. Subsequent finite element investigations into casing behavior under mudstone creep conditions indicated that immersion of mudstone in water instigated further shearing and deformation of the casing, culminating in premature well failure prior to water inundation. Notably, Von Mises stress levels exhibited a positive correlation with injection production ratios, with values exceeding critical thresholds leading to distinct modes of mechanical failure including shear-induced deformations, longitudinal tensile stress, and localized yielding near water wells. Maintenance of an optimal injection production ratio is identified as a key strategy for prolonging casing longevity in the region. To this end, recommendations include augmenting the casing wall thickness or enhancing the steel pressure specifications to mitigate casing damage progression, thereby extending the operational lifespan. This research serves as a pivotal theoretical framework for informing future development strategies aimed at mitigating and preempting casing failures in a low-permeability reservoir.

A quantitative assessment of the rehabilitative impact of Kevlar liners on offshore pipelines

Offshore pipelines transport oil and natural gas that contain different corrosive contaminants, leading to internal corrosion. The presence of internal metal loss problems can significantly degrade the load bearing capacity and reliability of pipelines over time. To address this issue, a Kevlar liner, offering a more accessible and cost-effective solution, is employed for repairing corroded pipelines. However, a quantitative assessment of the rehabilitating effect of Kevlar liner has been lacking. This study performed an investigation of the burst pressure of internal corrosion pipelines reinforced with the Kevlar liner through a combination of experimental testing and finite element numerical simulations. While previous research has primarily concentrated on single external defect, this study investigated the burst pressure of pipelines in the presence of both single and dual internal corrosion defects. The results indicated that the presence of two corrosion defects further decreases the burst pressure of the pipelines. This study examined the reinforcement mechanism of Kevlar liner, revealing that it largely postpones the stress development in corrosion zone, and thereby enhances the ultimate internal pressure load-bearing capacity of the pipelines. Moreover, the contact stress resulting from the friction between the Kevlar liner and the steel pipeline contributes to the overall improvement in the load-bearing capacity of the corroded offshore pipelines.

Research on the Optimal Design of Retaining Piles of a Wide Metro Tunnel Foundation Pit Based on Deformation Control

Engineers pay more and more attention to the economic benefits of foundation pit engineering. At present, the optimal design of the foundation pit supporting structure mainly focuses on strength and functional design, and there is no mature theoretical design method for deformation control. In this paper, a method for calculating the overall deformation of a foundation pit supporting structure based on the principle of minimum potential energy is proposed. Based on this method, the optimal design of the foundation pit of Guangzhou Baiyun District Comprehensive Transportation Hub Metro Station is realized. The deformation calculation results and optimization design scheme are validated by finite element numerical simulation and field monitoring data. The results show that the proposed theoretical algorithm predicts the pile deformation curves better than the finite element method, suggesting the proposed theoretical method is reasonable and the optimization scheme of the retaining pile is feasible. In the optimized design, the deformation of the foundation pit retaining pile is controlled by its push-back effect. The proposed deformation calculation method can realize the overall deformation calculation of the foundation pit supporting structure.

Influence patterns of shear-thickening gel on overload transmission: a numerical simulation study

To address signal aliasing challenges encountered by smart munition fuzes during penetration and layer counting processes, this study employs shear-thickening gel (STG) as a cushioning protective material for accelerometers. Based on the SHPB test platform, the test is designed and the finite element numerical simulation is carried out. The thickness of the shear-thickening gel specimen is treated as a research variable, enabling a systematic analysis of its mechanical filtering effect on overload acceleration. Furthermore, a spring-damping physical model for the cushioning structure is established, facilitating a theoretical exploration of the filtering mechanism on the acceleration signal. Research findings indicate that within a certain frequency bandwidth, STG materials can effectively filter out high-frequency signals received by the accelerometer. Additionally, augmenting the thickness of STG materials has a positive impact on achieving superior filtering performance within the high-frequency range.

Investigating parameters and ultimate load-bearing performance of novel anchor-box type cable–girder anchoring structure for large-span cable-stayed bridges

This study combined experimental research and finite element numerical simulations to investigate the ultimate bearing performance of a novel anchor-box-type cable–girder anchorage structure for large-span cable-stayed bridges. The stress distribution and deformation in the structure was assessed under varying loads, and the effects of the construction parameters on peak stress in key load-bearing welds were analyzed. The proposed design exhibits a safety margin of 3.2. The simulated and experimental stress distributions agree well with experimental data, although deformation and load-bearing capacity results are non-conservative. The maximum stress in the main load-bearing welds exhibits an exponential relationship with construction parameters.

Impact of Interfacial Transition Zone on Concrete Mechanical Properties: A Comparative Analysis of Multiphase Inclusion Theory and Numerical Simulations

With the increasing implementation of sustainable development strategies, recycled concrete (RC) has garnered attention in research circles due to its substantial environmental and economic advantages. The presence and properties of various interface transition zones (ITZs) in RC play a vital role in its mechanical properties. This research uses a combination of multiphase inclusion theory and finite element numerical simulation to investigate and compare the impact of ITZs on concrete’s mechanical properties. The multiphase inclusion theory offers a theoretical framework for understanding ITZ behavior in concrete, categorizing it into new mortar, old mortar, new ITZ, old ITZ, and natural aggregate based on meso-structure. With simplified RC at the mesoscale, the study accurately predicts the mechanical properties of RC by adjusting the elastic modulus, Poisson’s ratio, and thickness of new and old ITZ models. Through finite element simulation and theoretical validation, the study achieves a minimal error of 6.24% in predicting the elastic modulus and 1.75% in predicting Poisson’s ratio. These results highlight the effectiveness of multiphase inclusion theory in capturing the meso-structure characteristics of RC and forecasting its macroscopic mechanical behavior while comprehensively considering the complexity of ITZs.

Experimental-Numerical Assessment of Mechanical Behavior of Laboratory-Made Steel and NiTi Shape Memory Alloy Wire Ropes

The mechanical behaviors of laboratory-fabricated steel and superelastic shape memory alloy (SMA) wire ropes are assessed in this study through a comprehensive approach encompassing both experimental investigations and finite element (FE) numerical simulations. The assessment of steel wire ropes involves experimental scrutiny under sinusoidal cyclic loading and natural earthquake loading conditions. In parallel, SMA wire ropes’ behaviors are analyzed utilizing FE simulations employing the widely acknowledged ABAQUS software version 2020. The validation of all numerical simulations is undertaken against the experimentally observed behaviors. Moreover, full-scale steel wire ropes are subjected to shaking table tests to validate the simulations, facilitating a comparative analysis between the mechanical responses of SMA and steel wire ropes. The findings demonstrate that SMA wire ropes exhibit superelastic behavior akin to SMA wires, with marginal variations in overall response observed across distinct configurations, akin to steel wire ropes. Furthermore, augmenting the helix angle of SMA wire ropes results in reduced stress and increased strain when exposed to the El Centro earthquake scenario. Nevertheless, the mechanical response of SMA wire ropes closely mirrors that of a single wire.

Geomechanical Response Characteristics of Different Sedimentary Hydrodynamic Cycles—Exampled by Xujiahe Formation of Upper Triassic, Western Sichuan Basin

This study delves into the geomechanical responses of different sedimentary hydrodynamic cycles in deep tight sandstone formations. Employing core observation and thin section analysis, we quantitatively identified and characterized bedding planes, sedimentary microfacies, and tectonic fractures. Then, the intricate relationships between various architectural interfaces and geomechanical parameters were elucidated. Subsequently, utilizing finite element numerical simulation software, in situ stress and fracture parameters were derived. By identifying a fracture facies zone correlated with the sedimentary hydrodynamic cycle and production data, our findings unveil several key insights: (1) Geomechanical parameters (Young’s modulus, Poisson’s ratio, brittleness index) exhibited noteworthy variations within the T3x2−5 sand group, indicative of weak elasticity and robust plasticity. (2) The effective distance, influenced by diverse reservoir architecture interfaces, displayed variability, with each transition between peak-valley-peak or valley-peak-valley pinpointed as a distinct sedimentary hydrodynamic cycle. (3) In environments characterized by strong sedimentary hydrodynamics (between two level 3 architecture interfaces), fractures with larger strike angles and lower dip angles were observed to be more prevalent. (4) Three significant fracture faces—level I, level II, and level III—were discerned within the study area. Notably, reservoirs associated with level III exhibited characteristics suggestive of medium porosity and permeability, indicative of a gas layer. By thoroughly understanding the geomechanical response characteristics of formations such as the Xujiahe Formation, it is possible to guide the exploration and development of energy resources such as oil and natural gas. This helps to improve the efficiency and safety of resource extraction, promoting the sustainable utilization of energy.

Wave‐to‐Wire Model and Output Performance Analysis of a Point Absorption Wave Energy Converter Based on Triboelectric Nanogenerator

Herein, a point absorption wave energy converter (WEC) model combined with a contact–separation mode triboelectric nanogenerator (PAW TENG) is proposed. The model considers the interaction between the hydrodynamic forces captured by a point absorption conversion device of wave energy and the electrostatic force of the TENG. According to the established harmonic oscillation model of the buoy, the buoy's hydrodynamic properties and power generation characteristics are studied using finite element simulation and numerical simulation calculations. The best matching resistor circuit for PAW TENG in regular waves is investigated and analyzed. The electrical characteristics of the PAW TENG in the irregular wave spectrum are also solved. This work will provide design ideas for the experimental design of a novel point absorption WEC.

Numerical simulation and experimental study on the compression performance of a double layer cross-lattice-core sandwich structure

Lattice sandwich structures are widely used in many engineering fields for their excellent mechanical properties, such as aerospace, transportation, construction and so on. In this paper, a new type of double layer cross-lattice-core sandwich structure is proposed. Firstly, based on the nonlinear finite element software ABAQUS, a finite element numerical simulation analysis is processed to the sandwich structure unite cell under flat compression. The deflection pattern, load-displacement curve and internal energy-displacement curve are obtained to discuss the compression performance of the unite cell with some specific geometric parameters, such as arm length ratio, clamped angle and thickness of the cross-lattice-core. It is found that the compression performance of the lattice sandwich structure with an arm length ratio of 1 and an included angle of 90° is optimal. The energy absorption effect decreases slightly with the increases of the arm length ratio increases, which will be greatly reduced when the core angle is lower than 90°. Then the unite cell model of the structure is prepared using 3D printing technology, and its compression performance is investigated by flat compression tests and verified in comparison with the numerical analysis results. Finally, the load-bearing property of the structure is further verified by compression analysis of a plate structure consisting of several unite cells. The results provide a reference for the design and engineering application of the structure.

Progressive Failure Mechanism of Shield Tunnel Face in Complex Urban Geological Environment

The construction of multiple tunnels across inland rivers has had a significant influence on the improvement of the transportation infrastructure. The technology for constructing tunnels is progressing towards the development of larger cross-sections, longer distances, and the ability to withstand high hydraulic pressure in complex hydrogeological conditions, including high-permeability strata. In order to ensure the face stability of shield tunnels under high hydraulic pressure that crosses a fault fracture zone, it is necessary to study the progressive failure mechanism of shield tunnel faces induced by high hydraulic pressure seepage. This paper employs finite element numerical simulation software to methodically examine the variation in the characteristics of the water seepage field, limiting support force, and face stability failure mode of shield tunnels passing through fault fracture zones with high hydraulic pressure under varying fault fracture width zones. The results show that the formation hydraulic gradient will progressively widen when the tunnel face is located within the undisturbed rock mass and is advanced towards the area of fault fracture. This will raise the likelihood of instability in the shield tunnel and progressively raise the limiting support force on the tunnel face. Moreover, as the tunnel face nears the region of fault fracture within the undisturbed rock mass, the damage range increases gradually. In addition, due to the increase in seepage force, the angle between the failure area and the horizontal plane becomes more and more gentle. On the contrary, as the tunnel’s face moves closer to the undisturbed rock mass from the region of the fault fracture, the damage range gradually decreases, and the dip angle between the damage area and the horizontal plane becomes steeper and steeper due to the decreasing seepage force in the process. The study findings presented in this work are highly significant, both theoretically and practically, for the design and management of safety.

Effects of jet grouting on adjacent ground through numerical modelling

When using jet grouting for soil stabilization, effects of the installation process, such as erosion or remolding of the surrounding soil, need to be evaluated. This could be particularly important for applications in challenging ground conditions and low initial ground stability. This study investigates such effects in soft marine clays (undrained shear strength of 20-30 kPa) using numerical finite element simulations for both flat and sloping terrains. Additionally, an analytical solution using the modified Hagen-Poiseuille equation was developed to i) assess the pressure induced by the grout at any installed depth, and ii) illustrate the importance of maintaining a proper amount of return flow to avoid build-up of excessive pressures in the ground during the jet grouting process. Combining the analytical and the finite element (FE) simulations, jet grouting induced soil displacements were determined. For flat terrain, the area affected by jet grouting was approximately three to four times the column radius. By contrast, for the studied slope (inclination of 1:3), the influenced area increased to about 10 times the column radius. The influenced area could increase further for steeper slopes. The obtained displacement magnitude was also considerably larger for sloping terrain than for flat terrain. The reported findings relate to jet-grouting in marine clays (which is not the normal application of this technology) and cannot be extrapolated to other soils.

Analysis of the stability of a novel assembled working shaft support structure during pipe jacking construction: Experiments and numerical simulations

Pipe jacking construction is a commonly employed trenchless technique for laying underground pipes in urban areas. However, the conventional support structure for pipe jacking shafts poses challenges, including difficulty in reserving pipe jacking holes, susceptibility to tilting, low bearing capacity, and the potential for failure during construction. Taking the 5# pipe jacking shaft of the pipe jacking construction for diverting the Yellow River through the river at Niukouyu in Zhengzhou City as the research background, investigates the soil settle deformation around the working shaft and the mechanical response law of the supporting structure during the short-distance, long-distance and second long-distance pipe jacking construction by combining the practical engineering field test and finite element numerical simulation, and carries out a sensitivity analysis on the main design parameters affecting the stability of the supporting structure through orthogonal test. The findings reveal that during the second pipe jacking construction, stress and deformation of the supporting structure are higher than those observed in the first pipe jacking. Notably, support piles and waist beams at the entrance of the pipe jacking experience greater force, and the back and side walls undergo increased force and deformation in the later stages of pipe jacking, and support pile spacing is the main control factor affecting the mechanical performance of the novel support structure. The study concludes that monitoring and protection measures should be reinforced, particularly in areas prone to failure and damage during construction. The insights gained from this research can serve as a reference for designing, optimizing, and safely monitoring novel assembled pipe jacking shaft support structures.

Rate-dependent energy dissipation of graded viscoelastic structures fabricated by grayscale vat photopolymerization

A major benefit of additive manufacturing technologies is precise control over structural topologies and material properties, which allows to tailor, for instance, energy absorption and dissipation. While vat photopolymerization is generally restricted to a single material, grayscale masked stereolithography (gMSLA) allows to customize material behavior by grading the light intensity within a structure. This study investigates the impact and opportunities of grayscale grading strategies on the rate-dependent mechanical behavior of structures fabricated by gMSLA. Considering the viscoelastic nature of polymers, rate-dependent energy dissipation is explored, introducing a parametric linear viscoelastic constitutive model for varying grayscales. The investigation includes the comprehensive characterization of mechanical properties, numerical finite element simulation, validation through experimental procedures, and exploration of dissipation energy under different strain rates. In this way, a rational function successfully determines the critical strain rate at which the maximum dissipation occurs. Overall, the research offers a comprehensive investigation of the mechanical dissipation behavior of graded 3D printed structures, laying the foundation for further studies and advancements aimed at optimizing these structures for enhanced energy absorption capabilities.

Multi-parameter terahertz metamaterial sensors based on single-layer quarter ring patterns

A terahertz metamaterial sensor based on a single-layer quarter ring pattern is proposed to realize the simultaneous measurement of temperature, refractive index, and other parameters through a metasurface consisting of a patterned quarter ring gold layer and a silicon substrate. Field finite element numerical simulations were performed to analyse the transmission and reflection spectra. A two-parameter sensor combined with electric field distribution was realized using multiple peaks, thus revealing the physical mechanisms of the interaction in metamaterial structures. In addition, the effect of each geometrical parameter on the sensor performance is also simultaneously analysed to verify the optimal selection of structural parameters. On this basis, the sensors are theoretically analysed at different incidence and polarization angles, and the results show that the transmittance is not sensitive to the polarization angle. The simulation results show that the transmittance at the resonant frequencies of 0.48 THz, 0.64 THz, 0.79 THz and 1.04 THz are 93.7 %, 79.1 %, 38.4 % and 50.8 %, respectively, and the refractive index sensitivities at the resonant frequencies are 6.84 GHz/RIU, 7.42 GHz/RIU, 8.22 GHz/RIU and 32.11 GHz/RIU, and quality factor Q of 6.76, 9.63, 15.51 and 26.28, respectively. The devices proposed in this study are simple in structure, easy to fabricate, easy to miniaturize and can be integrated for sensitivity, multi-parameter sensors. They can be functionally designed for chemical and biological sensing and other important applications.

Numerical Simulation and Deformation Prediction of Deep Pit Based on PSO-BP Neural Network Inversion of Soil Parameters.

The finite element numerical simulation results of deep pit deformation are greatly influenced by soil layer parameters, which are crucial in determining the accuracy of deformation prediction results. This study employs the orthogonal experimental design to determine the combinations of various soil layer parameters in deep pits. Displacement values at specific measurement points were calculated using PLAXIS 3D under these varying parameter combinations to generate training samples. The nonlinear mapping ability of the Back Propagation (BP) neural network and Particle Swarm Optimization (PSO) were used for sample global optimization. Combining these with actual onsite measurements, we inversely calculate soil layer parameter values to update the input parameters for PLAXIS 3D. This allows us to conduct dynamic deformation prediction studies throughout the entire excavation process of deep pits. The results indicate that the use of the PSO-BP neural network for inverting soil layer parameters effectively enhances the convergence speed of the BP neural network model and avoids the issue of easily falling into local optimal solutions. The use of PLAXIS 3D to simulate the excavation process of the pit accurately reflects the dynamic changes in the displacement of the retaining structure, and the numerical simulation results show good agreement with the measured values. By updating the model parameters in real-time and calculating the pile displacement under different working conditions, the absolute errors between the measured and simulated values of pile top vertical displacement and pile body maximum horizontal displacement can be effectively reduced. This suggests that inverting soil layer parameters using measured values from working conditions is a feasible method for dynamically predicting the excavation process of the pit. The research results have some reference value for the selection of soil layer parameters in similar areas.

Model test and numerical analysis of shield tunnel face instability in sandy soil with strain-softening properties

During excavation of a shield tunnel in a sandy stratum, variations in the state of the sand affect the stratum stability of the shield tunnel face. When dense sand is loaded, its strength decreases gradually owing to dilatancy and it exhibits strain-softening behavior. Model tests and finite element numerical simulations were performed to clarify the failure characteristics of the tunnel face induced by strain-softening. In the model tests, a dense sand stratum was constructed; the movement of the supporting plate of the excavation face was controlled to simulate the progressive failure of the shield tunnel face. Based on a material state-dependent plasticity model, a numerical analysis model of the tunnel face was established to capture the strain-softening behavior of dense sand by including a nonlocal regularization technique. Comparing the model test and numerical simulation results, the evolution of the material state of the sand with progressive failure of the tunnel face was analyzed, and the failure characteristics caused by the strain-softening of sand were revealed.

Seepage process understanding at long-existing landslide dams through numerical analysis and hydrological measurements

Landslide dams consist of unconsolidated heterogeneous material and lack engineering measures to drain water and control pore water pressure. They may be porous and seepage through them could potentially lead to piping failure. In this research, the internal processes within a long-existing landslide dam are assessed under transient seepage force. The implemented approach includes a 3D finite element numerical simulation executing fully coupled flow-deformation and consolidation methods based on hydraulic data measurements and geotechnical laboratory tests. The nonlinear constitutive model ‘Hardening Soil’ is applied to accurately calculate the stress-induced pore water pressure, effective stress, deformation, and flow. Further, the possibility of slope failure due to seepage force is investigated through the strength reduction method. The results highlight the dependency of the seepage flow on the corresponding variation of the relative permeability and saturation in the soil mediums under different rates of seepage force. Small rates of seepage force, however, impose deformation at the dam's crown. High effective stress is obtained at negative small rates of seepage force where the long duration of fluctuation is modeled. In the drawdown simulation, there is a reverse relation between effective stress and the rate of the seepage force. Through the modeling process and based on the measured data, two seepage paths are detected within the landslide dam, while their activation depends on the lake level. The modeling approach and the required data analysis are suggested for utilization in further studies regarding the seepage process understanding at the long-existing landslide dams and their hazard assessments in addition to the common geomorphological approaches.

Study on the Ionic Transport Properties of 3D Printed Concrete

Three-dimensional printed concrete (3DPC) is an anisotropic heterogeneous material composed of a concrete matrix and the interfaces between layers and filaments that form during printing. The overall ion transport properties can be characterized by the equivalent diffusion coefficient. This paper first establishes a theoretical model to calculate the equivalent diffusion coefficient of 3DPC. Verification through numerical calculations shows that this theoretical model is highly precise. Based on this, the model was used to analyze the effects of dimensionless interface parameters on the equivalent diffusion coefficients in different directions of 3DPC. Finally, the dynamic ionic transport properties of 3DPC were investigated through finite element numerical simulation. The results of the dynamic study indicate that interfaces have a significant impact on the ion distribution and its evolution within 3DPC. The product of the interface diffusion coefficient and interface size can represent the ionic transport capacity of an interface. The stronger the ionic transport capacity of an interface, the higher the ion concentration at that interface. Due to the “drainage” effect of lateral interfaces, the ion concentration in the middle of 3DPC with a smaller equivalent diffusion coefficient is higher than that in 3DPC with a larger equivalent diffusion coefficient.