Finite Element Numerical Simulation Software Research Articles

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.

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.

Simulation of Temperature Field in Railway Embankment Filling in Seasonal Frozen Areas

In order to examine the temperature field changes of railway roadbed in seasonal frozen areas under different filling heights, a typical cross-sectional model of Chagan Lake Station railway was established using the finite element numerical simulation software ANSYS. The temperature parameters of the roadbed were calculated based on the atmospheric temperature of local meteorological data, and the temperature fitting curve was used as the temperature boundary condition. Under different filling quantities, the temperature field of the roadbed was simulated for a freeze-thaw cycle, and the temperature curve changes were analyzed. The results indicate that as the filling height increases, the freezing depth line also increases accordingly. When the filling height is 2.5m, the maximum freezing depth at the roadbed is 1.7m. In a freeze-thaw cycle, the freezing process of soil develops from top to bottom, and the melting process develops in two directions: from top to bottom and from a certain depth of soil upwards, with the characteristics of "unidirectional freezing and bidirectional melting". The freezing period starts from early December of each year to the end of March of the following year, and completely melts until the end of April. Melting develops faster than freezing.

Stability analysis and support scheme design of the landslide behind Yang Erping's house in Zuopo Village

The stability of the landslide behind Yang'erpingfang, Zuopo Village, Luluo Town, Xingtai County is calculated and analyzed using the strength reduction method with MIDAS GTS finite element numerical simulation software, based on the treatment of Zhuwenping landslide in Bai'an Township, Xingtai County.

Study on Failure Energy per Unit Area of Concrete Specimens Based on Minimum Energy Dissipation Theory.

In order to study the strength change of concrete specimens under different loading conditions, based on the principle of minimum energy dissipation, the damage energy per unit area of concrete was studied. By using finite element numerical simulation software for concrete specimens with different failure modes of tension, pressure, bending and torsion, a double-broken line damage constitutive model is adopted. The failure forms of concrete specimens under different loading conditions, as well as the failure area and failure energy of each specimen during loading, are simulated and analyzed. The failure energy per unit area under different failure modes was quantitively calculated, the relationship between the failure area and failure energy consumption under different failure modes was analyzed. The results show that, under different failure modes, the failure area of concrete specimens is different, the energy consumed during failure is different, and the strength is different. However, no matter how the failure mode changes during the failure process, the failure energy W per unit area remains constant and fluctuates in the range of 2.0~6.0 mJ/cm2, which is related to the physical properties of concrete itself.

Ultra-broadband high-efficiency cross-polarization conversion and amplitude-phase manipulable multi-functional wavefront manipulation based on terahertz metasurface

Drawing upon the physical phenomenon of polarization transformation, this paper proposes an ultra-broadband, high-efficiency linear polarization converter composed of a metallic grating, an L-shaped metallic patch, and a dielectric substrate. The polarization conversion properties have been scrutinized using the finite element numerical simulation software CST. The computational outcomes reveal that the polarization converter operates within the frequency range of 0.5 THz to 1.8 THz, exhibiting a relative bandwidth of 113%, a transmission coefficient exceeding 0.87, a polarization conversion efficiency approaching 100%, and a phase coverage spanning 360°. Furthermore, a Fabry–Perot interference model was established utilizing Matlab to corroborate the concurrence between the theoretical analysis and the numerical findings. The polarization converter metasurface amalgamates both phase and transmission amplitude variations to accomplish not only a two-dimensional focusing lens operating between 1.55 THz and 1.65 THz, but also a spatial imaging capability utilizing transmission amplitude variation within the 0.5 THz to 1.15 THz range. The outcomes demonstrate that the devised metasurface exhibits ultra-broadband and high transmission efficacy, thus providing novel insights for the versatility of terahertz wave polarization and phase manipulation.

Study on the forming characteristics of polytetrafluoroethylene/copper jet with different preparation processes

In this paper, PTFE/Cu composite material for liner is taken as the research object, and the preparation process and jet forming characteristics of PTFE/Cu composite liner are studied. The liners were prepared by extrusion molding, molded sintering and hot-pressing sintering. Due to different preparation processes, different microstructures of the liner can occur, including defects such as pores and microcracks, resulting in different strength and density of the liner, leading to differences in the forming characteristics of the jet. Therefore, the forming process of the jet was simulated by the finite element numerical simulation software. It was found that there was obvious radial expansion effect in the head of the jet, but with the increase of density, the radial expansion effect was weakened, and the jet velocity decreased gradually. The strength and densification of the shaped charge liner prepared by different processes were different. The densification of the molded sintering liner was generally better than that of the other two kinds of shaped charge liners. As a result, the velocity of the jet formed by the molded sintering liner is always the highest, with a numerical simulation velocity of 6642 m/s and an experimental velocity of 6534.7 m/s. The second is the jet of the hot-pressing sintering liner and the lowest velocity is the jet of the extrusion molding cover, with a numerical simulation velocity of 6482 m/s, while the experimental velocity is only 6397.9 m/s. The jet velocity measured by the pulse X-ray experiment was compared with the velocity of the numerical simulation, and the error was within 2.96%, which verifies the accuracy of the numerical simulation.

Research on key influencing factors of scrap tire-soil retaining wall

The emergence of STSW can reduce the amount of materials such as sand and gravel and improve its performance compared to conventional retaining walls. In order to investigate the key influencing factors of STSW, the compression modulus of Scrap tire-soil module (STSM) was obtained by designing a full-scale compression strength test. The shear strength was obtained through a scale-test of shear strength. On this basis, the gravity STSW was established by the finite element numerical simulation software, and the bearing capacity of STSW was analyzed. Further, the key influencing factors of the bearing capacity of STSW were analyzed by designing orthogonal tests, and suggestions for parameter selection for designing STSW were given. The results show that: the friction angle of STSW increases and becomes larger with the increase of area ratio, and the cohesion decreases and decreases with the increase of area ratio; with the increase of compression modulus and internal friction angle, the maximum shear stress inside the retaining wall increases, while the cohesive force has almost no effect on the maximum value of shear stress inside the retaining wall; when constructing retaining walls, it is recommended to use Internal filling material with large compression modulus, small area ratio and its in the range of [0.3–0.5] of STSM.

Case Study on the Effect of Delay-Time Differences between Columns during Blasting Demolition of RC Structures with a Small Height-to-Width Ratio

Blasting demolition plays an important role in building demolition, and numerical-simulation methods can facilitate the optimization of blasting-demolition schemes. To study the influence of the delay-time differences between columns on the blasting and demolition effect of small-aspect-ratio reinforced concrete structures, a separated common node model of small-aspect-ratio structures was established through the finite element numerical simulation software ANSYS/LS-DYNA. The results showed that the delay-time differences between columns influenced the position of the plastic hinge of the structure, and the plastic hinge was located at the end of the back span beam when the delay-time difference was 0.1~0.3 s. The plastic hinge moves to the top of the back row of the columns of the first floor at 0.4~0.8 s and to the end of each span of the beams at 0.9 s. When the plastic hinge is at the end of the beam, the orientation axis is located at the bottom of the back row of the columns of the first floor. Moreover, when the plastic hinge is at the top of the back row of the columns of the first floor, the orientation axis is at the same position as the plastic hinge. The delay-time difference also affects the recoil distance of the structure as well as the height of the burst pile. To ensure that the structure can be successfully collapsed under the preconditions of disintegration, for a small height-to-width ratio structure, the appropriate intercolumn delay-time difference is 0.4~0.5 s.

Dynamic Response Law and Failure Mechanism of Slope with Weak Interlayer under Combined Action of Reservoir Water and Seismic Force

The southwestern region of China is close to the Eurasian earthquake zone. Many engineering areas in southwestern China are affected by earthquakes and are close to the epicenter of earthquakes that occur in this region. During earthquakes, slopes with weak interlayers are more likely to cause large-scale landslides. In response to the low stability of slopes with weak interlayers in reservoir dam areas, the dynamic response law and failure mechanism of weak interlayered slopes under the combined action of reservoir water and seismic forces were studied through shaking table model tests and finite element numerical simulation software. The height of the water level and the size of the seismic waves were changed during these tests. The research results indicate that seismic waves are influenced by weak interlayers and are repeatedly superimposed between the weak interlayers and the slope surface, resulting in an acceleration amplification effect that increases by approximately 1.8 times compared to homogeneous slopes. Vertical earthquakes have a significant impact on the dynamic response of slopes, and their peak acceleration amplification coefficient can reach 0.83 times the horizontal peak acceleration. The stability of weak interlayers during earthquakes is the worst within the range of the direct action of reservoir water. The failure mode of a slope is as follows: earthquake action causes cracking in the upper part of the slope, and as the earthquake increases in intensity, and the infiltration of reservoir water intensifies, the cracks expand. The soft and muddy interlayer in the front section of the slope forms a sliding surface, and ultimately, the sliding failure forms an accumulation body at the foot of the slope. In reservoir dam areas, the stability of a slope is closely related to the engineering safety of the reservoir dam. Therefore, when a strong earthquake and the water level in a reservoir jointly affect a weak-interlayer slope, the slope is in the stage of plastic deformation and instability. The stability of the slope may be overestimated, and the slope is likely vulnerable to sliding instability, which needs to be monitored and treated.

The effect of macropore morphology of actual anecic earthworm burrows on water infiltration: A COMSOL simulation

Soil macropores impact water infiltration. Due to the complexity of macropore structure, most macropore soil water dynamic simulations use simplified assumptions, such as cylindrical shape macropores, which do not well represent fluid flow in actual macropores. In this study, 3D digitalized structures of actual anecic earthworm burrows are obtained by combining tin casting with a 3D scanner. The in-situ earthworm burrows are imported into finite element numerical simulation software (COMSOL). Based on simulations, the effects of burrow spatial characteristics, such as number (N), length (L), average diameter (Davg), and tortuosity (τ) of partially penetrating (non-through) and fully penetrating (through) burrows, on water infiltration are clarified. The N value of non-through burrows correlates significantly with preferential flow infiltration rates (r = 1). In addition to the N value, the L value of non-through burrows is an important factor affecting preferential flow infiltration (r = 0.99). For the through earthworm burrows, macroporosity (εp) is the best predictor of Ksat. Unlike Davg and L, the τ of burrows has a low correlation coefficient with Ksat. The increase of adjacent infiltration is N-dependent for non-through burrows and is N-independent for through burrows. These findings demonstrate that actual macropore morphology can be used for numerical simulations, which provides a new pathway forward for pore-scale soil water dynamic research.

Adaptability Evaluation of Rotary Jet Grouting Pile Composite Foundation for Shallow Buried Collapsible Loess Tunnel

The deformation control effect of loess tunnel composite foundation plays an important role in optimization design and reinforcement effect evaluation. Systematically evaluate the adaptability of the composite foundation of jet grouting pile in shallow collapsible loess tunnel. Taking the shallow buried section of Fujiyao Tunnel with a buried depth of 20 m as an example, using MIDAS finite element numerical simulation software, the foundation deformation control during construction and settlement control after construction are systematically studied, the differential deformation control is analyzed, and the reinforcement effect of the tunnel bottom is evaluated. The results show that the uplift displacement can be controlled by changing the pile length and increasing the replacement ratio. The combination of long and short piles can significantly reduce the uneven settlement and plastic zone of the foundation. The uneven settlement of 9 mm can be used as the evaluation index of the composite foundation reinforcement effect in a shallow buried section of the loess tunnel.

Seismic effects of loess slopes using physical modeling and numerical simulation

A large shaking-table test of a loess slope with a geometric similarity ratio of 1:25 was established as a prototype model of the loess landslide in Yongguang Village, which occurred during the 2013 Minxian–Zhangxian M6.6 earthquake. Combined with finite-element–numerical-simulation software, the dynamic response characteristics of the loess slope were compared and analyzed under the action of the same Minxian horizontal seismic wave characteristics. The results show that the seismic waves observed in the model test and numerical simulation clearly exhibit non-linear amplification effects along the slope surface and in the loess vertical direction, reaching a maximum value at the top of the slope. The ground-motion acceleration waveforms observed in the model test and numerical simulation at each observation point of the slope were similar, and the horizontal peak ground acceleration (PGA) obtained from the numerical simulation at these locations was larger than that of the slope top. The horizontal aSimulated-max/aModel-max values ranged from 1.5 to 2.2 and the vertical aSimulated-max/aModel-max values were more discrete in the range of 0.5–2.7, while the horizontal and vertical-PGA-amplification factors at the slope surface were similar. The Fourier amplitude and main horizontal frequencies were mainly distributed between 0.1 and 12 Hz. The Fourier spectrum observed in the model test was approximately 0.5 Hz, while the main frequency of the vertical Fourier spectrum was approximately 4 Hz. Additionally, the main frequency of the horizontal and vertical Fourier spectra observed in the numerical simulation was approximately 2 and 1.2 Hz, respectively.

Numerical simulation analysis for the effect of water content on the intelligent compaction quality of roadbed

<abstract> <p>In the process of intelligent compaction of roadbeds, the water content of the roadbed is one of the important influencing factors of compaction quality. In order to analyze the effect of water content on the compaction quality of roadbeds, this paper is developed by secondary development of Abaqus finite element numerical simulation software. At the same time, the artificial viscous boundary was set to eliminate the influence of boundary conditions on the results in the finite element modeling process, so that the numerical simulation can be refined to model. On this basis, the dynamic response analysis of the roadbed compaction process is performed on the finite element numerical simulation results. This paper established the correlation between compaction degree and intelligent compaction index CMV (Compaction Meter Value) and then analyzed the effect of water content on the compaction quality for the roadbed. The results of this paper show that the amplitude of the vertical acceleration is almost independent of the moisture content, and the vertical displacement mainly occurs in the static compaction stage. The vertical displacement changes sharply in the first 0.5 s when the vibrating wheel is in contact with the roadbed. The main stage of roadbed compaction quality increase is before the end of the first compaction. At the end of the first compaction, the roadbed compaction degree increased rapidly from 80% to 91.68%, 95.34% and 97.41%, respectively. With the increase in water content, the CMV gradually increased. At the end of the second compaction, CMV increased slightly compared with that at the end of the first compaction and stabilized at the end of the second compaction. The water content of the roadbed should be considered to be set slightly higher than the optimal water content of the roadbed by about 1% during the construction of the roadbed within the assumptions of this paper.</p> </abstract>

STUDY ON PREFABRICATED CONNECTOR OF DOUBLE-LAYER RECIPROCAL FRAME

A kind of lattice-type steel member is presented and a kind of prefabricated connector suitable for the connection between lattice-type steel members are proposed. The mechanical properties of the connectors are analyzed by using the finite element numerical simulation software ABAQUS. The connectors meet the design goals of the stiffness of connection stronger than members. Parameterized analysis is carried out on the prefabricated connector, and the flexural stiffness expression of the connector is obtained. The suggested values of each component of the prefabricated connectors are given based on the size of connected lattice-type steel members.

Dynamic analysis of shallow embedded jacket offshore structures under wave forces

In this paper, the stress of temporary structures during the construction of offshore structures is analysed through the study of the temporary buried drilling platform. Firstly, The fluid flow is controlled by the modified RANS equation (the Reynolds-average equation) and Forchheimer saturation resistance model, the free surface is tracked by volume of fluid (VOF) method, and the κ–ε closed equation is solved. The finite element numerical simulation software is used to establish a multi-factor coupling calculation model of the stress of temporary buried offshore structures. The maximum wave force of temporary structures under different buried depths, different structural sizes, different load periods and amplitudes is simulated, and the stability of Long Yuan wind turbine pile foundation drilling platform under the maximum wave force is analyzed. The research results of this paper can be used for dynamic design and analysis of wave load conditions, structural design dimensions and load periods in the engineering design of embedded offshore structures in shallow water under wave loads, so as to calculate the wave load values of structures reasonably and provide concrete theoretical basis for the establishment of offshore systems.

Roof Fractures of Near-Vertical and Extremely Thick Coal Seams in Horizontally Grouped Top-Coal Drawing Method Based on the Theory of a Thin Plate

During the mining process of the near-vertical seam, there will be movement and collapse of the “roof side” rock layer and the “overlying coal seam,” as well as the emergence of the “floor side” rock layer roof which is more complicated than the inclined and gently inclined coal seams, which causes problems with slippage or overturning damage. With the increase of the inclination of the coal seam, the impact of the destruction of the immediate roof on the stope and roadway gradually becomes prominent, while the impact of the destruction of the basic roof on the stope and the roadway gradually weakens. The destruction of the immediate roof of the near-vertical coal seam will cause a large area of coal and rock mass to suddenly rush to the working face and the two lanes, resulting in rapid deformation of the roadway, overturning of equipment, overturning of personnel, and even severe rock pressure disaster accidents, all of which pose a serious threat to coal mine safety and production. It is necessary to carry out research on the mechanical response mechanism of the immediate roof of near-upright coal seams, to analyse the weighting process of steeply inclined thick coal seam sub-level mining. A four fixed support plate model and top three clamped edges simply supported plate model for roof stress distribution are established before the first weighting of the roof during the upper and lower level mining process. The bottom three clamped edges simply supported plate model and two adjacent edges clamped on the edge of a simply supported plate model are established for roof stress distribution before periodic weighting of the roof during the upper and lower level mining process. The Galerkin method is used to make an approximate solution of deflection equation under the effect of sheet normal stress, and then roof failure criterion is established based on the maximum tensile stress strength criterion and generalized Hooke law. This paper utilizes FLAC3D finite element numerical simulation software, considering the characteristics of steeply inclined thick coal seam sub-level mining. It undertakes orthogonal numerical simulation experiment in three levels with different depths, coal seam angles, lateral pressure coefficient, and orientation of maximum horizontal principal stress, and translates roof stress of corresponding 9 simulation experiment into steeply inclined roof normal stress. We conclude that the distribution law of normal stress along dip and dip direction of a roof under the circumstance of different advancing distances and different sub-levels. The caving pace of first weight and periodical weight were counted under the effect of the roof uniform normal stress. It can better predict the weighting situation of the working face and ensure the safe, efficient, and sustainable mining of coal mines.

Quantitative prediction of structural fractures in the Paleocene lower Wenchang formation reservoir of the Lufeng Depression

Currently, the lower Paleogene Wenchang formation in the Lufeng Depression is the primary focus of reservoir development. The structural fractures that have formed inside of it not only serve as the principal path for oil migration, but also as oil storage space. As a result, the distribution features of structural fractures are crucial for future petroleum exploration and development in the Lufeng Depression. At the same time, with the quantity of conventional reservoirs in the Lufeng Depression on the decline, it is critical to determine the fracture distribution criteria for deep unconventional reservoirs. In this work, the lower Paleogene Wenchang formation in the Lufeng Depression is used as the research stratum. Then, based on existing logging data for the research region, the distinct physical properties of different rock kinds are calculated. The simulation results of the paleotectonic stress field in the study area using the finite element numerical simulation software ANSYS show that the high-value areas of maximum principal stress are the high-value areas of the uplift belt and low uplift, and the areas with low maximum principal stress are the low-value areas of Lufeng 13 Sag and the gentle slope belt in the north of Lufeng middle-low uplift. The fracture density is quantitatively predicted after the stress field simulation, which shows good agreement between the anticipated and actual observed values, and an average error of 13.61%. The predicted findings may provide new ideas for future petroleum exploration. cited as: Li, H., Yu, F., Wang, M., Wang, Y., Liu, Y. Quantitative prediction of structural fractures in the Paleocene lower Wenchang formation reservoir of the Lufeng Depression. Advances in Geo-Energy Research, 2022, 6(5): 375-387. https://doi.org/10.46690/ager.2022.05.03

Research on Optimal Design of Pressure Relief Blasting Parameters of Middle Rock Pillar in Near Vertical Coal Seam

In order to prevent the occurrence of rockburst accident caused by the stress concentration of the middle rock pillar in the near-vertical coal seam and affect the mining of the lower coal seam in the known goaf, taking Wudong coal mine in Xinjiang as an example, based on the special geological and mining conditions of the mine, through the simulation experiment and numerical analysis of similar materials, the variation law of the stress distribution of the middle rock pillar with the increase in mining depth is analyzed. A presplitting blasting pressure relief technology of drilling a construction chamber in the middle of the rock column is proposed to form a buffer layer in the rock wall and effectively prevent the continuous downward movement of stress. In order to achieve the optimal blasting effect, the finite element numerical simulation software GTS NX is used to carry out the chamber rock mass hole blasting simulation experiment, the three influencing factors of blasting effect, such as charge quantity, blast hole diameter, and blast hole row spacing, are comprehensively considered, and the optimal pressure relief blasting parameters with the radius of borehole blasting crack circle and the width of plastic zone as the index are determined. The experimental results show that under the condition of determining the physical and mechanical parameters of coal and rock, the blast hole diameter has become the main control factor affecting the blasting pressure relief effect, followed by the charge amount. Through 9 groups of orthogonal tests, the optimal combination of the three indexes is 6 kg charge, 150 mm hole diameter and 10 m hole row spacing. The practical engineering application results of the combination are tested by downhole transient electromagnetic detection (TEM), and it is found that the pressure relief effect is obvious. The preliminary conclusions of the study have guiding significance for the pressure relief and erosion prevention of the middle rock pillar in Wudong mine and provide experimental methods and theoretical basis for the determination of blasting pressure relief parameters in other mines.

Analysis and Simulation of Natural Frequencies of Slightly Curved Pipes

For the transverse vibration of slightly curved pipe, a dynamic mechanical model based on Timoshenko beam theory is established for the first time. The natural vibration characteristic of slightly curved pipe under the influence of the fluid flow is analyzed. With generalized Hamiltonian principle, the governing equation of the transverse vibration of slightly curved pipe considering the effect of fluid-structure coupling is derived. Based on Galerkin truncation, the natural frequencies of slightly curved pipe are obtained by generalized eigenvalue method. Effects of the fluid velocity and the initial deflection on the natural vibration characteristic of the pipe are studied. The equivalent stiffness and equivalent damping method base on finite element simulation for the natural vibration of the slightly curved pipe are developed. Then through the finite element numerical simulation software, the results of Galerkin truncation method and the effectiveness of the Timoshenko model are verified. The results show that fluid velocity and the initial bending degree both have significant effects on the natural frequencies of slightly curved pipe.