Petrel2ANSYS: Accessible software for simulation of crustal stress fields using constraints provided by multiple 3D models employing different types of grids

Crustal stresses play an important role in both exploration and development in the oil and gas industry.However,it is difficult to simulate crustal stress distributions accurately,because of the incompatibilities that exist among different software.Here,a series of algorithms is developed and integrated in the Petrel2ANSYS to carry out two-way conversions between the 3D attribute models that employ corner-point grids used in Petrel and the 3D finite-element grids used in ANSYS.Furthermore,a modified method of simulating stress characteristics and analyzing stress fields using the finite-element method and multiple finely resolved 3D models is proposed.Compared to the traditional finite-element simulation-based approach,which involves describing the heterogeneous within a rock body or sedimentary facies in detail and simulating the stress distribution,the single grid cell-based approach focuses on a greater degree on combining the rock mechanics described by 3D corner-point grid models with the finely resolved material characteristics of 3D finite-element models.Different models that use structured and unstructured grids are verified in Petrel2ANSYS to assess the feasibility.In addition,with minor modifications,platforms based on the present algorithms can be extended to other models to convert corner-point grids to the finite-element grids constructed by other software.

Modeling and simulation of 3D thermal stresses of large-sized castings in solidification processes

When heavy machines and large scaled receiver system of communication equipment are manufactured, it always needs to produce large-sized steel castings, aluminum castings and etc. Some defects of hot cracking by thermal stress often appear during solidification process as these castings are produced, which results in failure of castings. Therefore predicting the effects of technological parameters for production of castings on the thermal stress during solidification process becomes an important means. In this paper, the mathematical models have been established and numerical calculation of temperature fields by using finite difference method (FDM) and then thermal stress fields by using finite element method (FEM) during solidification process of castings have been carried out. The technological parameters of production have been optimized by the results of calculation and the defects of hot cracking have been eliminated. Modeling and simulation of 3D thermal stress during solidification processes of large-sized castings provided a scientific basis, which promoted further development of advanced manufacturing technique.

NUMERICAL SIMULATION OF STRESS-STRAIN DISTRIBUTIONS FOR WELD METAL SOLIDIFICATION CRACKING IN STAINLESS STEEL

This paper analyzed the characteristics of welding solidification crack of stainless steels,and clearly re- vealed the the of the deformation in the molten - the pool and the solidification shrinkage on the stress - strain fields in the trail of molten - weld pool.Moreover, rheologic properties of the alloys in solid - liquid zone were also obtained by measuring the hading and unloading deform curves of the steels.As a result, a numerical model for simulation of stress - strain distributions of welding solidifi- cation crack was developed. On the basis of the model,the thesis simulated the driving force of solidifi- cation crack of stainless steels, that is, stress - strain fields in the trail of molten-weld pool with fi- nite element method.

Numerical simulation of the tropical Pacific response to equatorial wind stress anomalies

The results of the tropical Pacific response to the sudden onset of the equatorial wind stress anomalies are discussed. The ocean model is a barotropic, non-linearized one that includes reduced-gravity and an equation for the temperature of the ocean mixed-layer. The experiments are based on a state of equilibrium reached through a long running under the action of annual mean wind stress. There are two kinds of westward wind intensity regions: the whole tropical Pacific and the western tropical Pacific, which are all between latitude 6. 8癗 and 6. 8癝.In these cases, the results show that the positive sea surface temperature (SST) anomalies in the Eastern Pacific and the negative SST anomalies in the Western Pacific are produced, and the positive SST anomalies propagate eastward, just as those observed during the actual El Nino phenomena. The propagations of the Kelvin waves and Rossby waves in the ocean are discussed.Another experiment is also carried out in simulating the process of the decay of El Nino event after the weakened Equatorial Pacific eastern winds returned to normal. The results are similar to the observations, too.

Numerical simulation of cutting process of the slice components cutting machine

By the use of ANSYS/LS-DYNA FEA software,numerical simulation on the cutting process of cutting plates with a reamer was carried out in the paper. The logical improvement was brought forward and the phenomenon of stress concentration was deceased by weighted analysis. The effects of different cut velocities and cutting thickness on life-spans of reamers were investigated, and the cutting parameters were optimized to satisfy the cutting precision and cutting efficiency. The study will provide a guide for the practical production.

Responses of the Anthocyanin and Osmolyte Contents of the Capsicum annuum Cultivars Planted in Wenshan Prefecture of Yunnan Province to the Drought Stress Simulated by PEG-6000

[Objective] The aim of this study was to compare the drought resistance difference of the main Capsicum annuum cultivars planted in Wenshan Prefecture of Yunnan Province. [Method] The total anthocyanin, soluble sugar, soluble protein and free proline contents of the leaves of the five main C. annuum cultivars planted in Wenshan Prefecture, i.e., ZSZ75-1, ZSZ49-1-1, 12WS-18-1, 12ZH01 and ZS130, under the drought stress simulated by PEG-6000 were studied by using spectrophotometry, and the total osmoregulation abilities of the leaves were evaluated by using subordinate function. [Result] Under the drought stress simulated by PEG-6000, the total anthocyanin contents of the veins, mesophylls and whole leaves of the five cultivars were all as： ZSZ75-1〉ZSZ49-1-1〉12WS-18-1〉12ZH01〉ZS130, the soluble sugar contents as： 12ZH01 〉12WS-18-1 〉ZS130〉ZSZ49-1-1 〉ZSZ75-1, the soluble protein contents as： ZSZ75-1〉12ZH01〉ZSZ49-1-1〉12WS-18-1〉ZS130, and the free proline contents as： ZS130〉ZSZ75-1 〉12WS-18-1 〉ZSZ49-1-1 〉12ZH01. Furthermore, the differences among the total anthocyanin contents of the veins, mesophylls and whole leaves of the five cultivars all reached the extremely significant levels, whereas the differences among the contents of the soluble sugar, soluble protein and free proline did not reach the significant levels, the correlation degrees among the contents of the total anthocyanin, soluble sugar, soluble protein and free proline of different cultivars were also inconsistent, and the total osmoregulation abilities of the five cultivar leaves were as： ZSZ75-1 〉7SZ49-1-1 〉12WS-18-1 〉ZS130〉12ZH01. [Conclusion] The anthocyanin accumulation in the specific organs of the five C. annuum cultivars in Wenshan Prefecture is beneficial to the strengthening of the leaf osmoregulation abilities under drought stress, favoring the drought tolerance of the cultivars.

Simulation of casting deformation based on mold surface element method

Deformation of casting during the solidification process has puzzled many engineers and scientists for years. In order to attain the goal of near-net forming by casting, numerical simulation is a powerful tool. Traditional methods compute the thermal stress of both the casting and the mold. This method suffers the problem of massive calculation and failure of convergence. This paper proposes an improved Mold Surface Element Method, the main idea of which is to use the surface elements instead of body elements to express the interactions between the casting and the mold. The proposed method shows a high computation efficiency and provides satisfactory precision for engineering. Two practical casting products were used to verify the proposed method. The simulated results agree well with those observed in practical products. The proposed method is believed to benefit production practice and to provide theoretical guidance.

Numerical analysis of the effects of rock bolts on stress redistribution around a roadway

Besides opening geometry, in situ stress and material properties, opening support also has significant effects on stress redistribution around a roadway. To investigate these effects of rock bolts on the stress redistribution around a roadway, a series of numerical studies were carried out using the finite difference method. Since the stress changes around a roadway caused by rock bolting is small relative to the in situ stress, they cannot obviously be observed in stress contour plots. To overcome this difficulty, a new result processing methodology was developed using the contouring program Surfer. With this methodology, the effects of rock bolts on stress redistribution can obviously be analyzed. Numerical results show that in the three patterns of rock bolts installed in the roof, in the roof and the two lateral sides, and in all the four sides of the rectangular roadway, the maximum stress magnitude of the increase is 0.931 MPa, 2.46 MPa,and 6.5 MPa, respectively; the bolt number of 5 can form an integrated ground arch; the appropriate length and pre-tensioned force of the rock bolt is 2.0 m and 60 k N, respectively. What is more, the ground arch action under the function of rock bolting is able to be effectively examined. The rock bolts dramatically increase the minor principal stress around a roadway which results in significant increase in material strength. Consequently, the major principal stress that the material can carry will greatly increase.With adequate supports, an integrated ground arch which is critical for the stability of roadway will be formed around the roadway.

MD simulation of a copper rod under thermal shock

In this paper, thermoelastic problem of onedimensional copper rod under thermal shock is simulated using molecular dynamics method by adopting embedded atom method potential. The rod is on axis x, the left outermost surface of which is traction free and the right outermost surface is fixed. Free boundary condition is imposed on the outermost surfaces in direction y and z. The left and right ends of the rod are subjected to hot and cold baths, respectively. Temperature, displacement and stress distributions are obtained along the rod at different moments, which are shown to be limited in the mobile region, indicating that the heat propagation speed is limited rather than infinite. This is consistent with the prediction given by generalized thermoelastic theory. From simulation results we find that the speed of heat conduction is the same as the speed of thermal stress wave. In the present paper, the simulations are conducted using the large-scale atomic/molecular massively parallel simulator and completed visualization software.

Simulation of hydrogen diffusion in welded joint of X80 pipeline steel

Hydrogen diffusion coefficients of different regions in the welded joint of X80 pipeline steel were measured using the electro-chemical permeation technique. Using ABAQUS software, hydrogen diffusion in X80 pipeline steel welded joint was studied in consideration of the inhomogeneity of the welding zone, and temperature-dependent thermo-physical and mechanical properties of the metals. A three dimensional finite element model was developed and a coupled thermo-mechanical-diffusion analysis was performed. Hydrogen concentration distribution across the welded joint was obtained. It is found that the postweld residual hydrogen exhibits a non-uniform distribution across the welded joint. A maximum equivalent stress occurs in the immediate vicinity of the weld metal. The heat affected zone has the highest hydrogen concentration level, followed by the weld zone and the base metal.Simulation results are well consistent with theoretical analysis.

DEVELOPMENT AND APPLICATION OF ADAPTIVE MESH TECHNIQUE IN THREE DIMENSIONAL NUMERICAL SIMULATION OF WELDING PROCESS

The adaptive mesh mesh technique is developed and applied in three dimensional numerical simulation of welding process on the base of the commercial software. Special user subroutine is worked out to accom- plish this function.This technique can make the dense mesh moving simultaneously with the heat source while the other area of the structure with much coarser mesh, greatly reducing the number of nodes and elements in the analysis.Temperature field,displacement and stress distributions during welding pro- cess me analyzed by FEM method with adaptive mesh and the analysis is also conducted with normal FEM method. The temperature field,displacement and stress distributions obtained with both methods are shown in contrast. The results show that the temperature fields and the displacement distributions of simulation on adaptive mesh correspond wery well with that of without adaptive mesh. Though the stress distributions have some difference,but the trends of the stress distribution are corresponding.The com- parison of the computing time of the two meshes indicates that the adaptive that the adaptive mesh can greatly reduce the calculation time when used for welding process.

Molecular dynamics study of thermal stress and heat propagation in tungsten under thermal shock

Using molecular dynamics （MD） simulation, we study the thermal shock behavior of tungsten （W）, which has been used for the plasma facing material （PFM） of tokamaks. The thermo-elastic stress wave, corresponding to the collective displacement of atoms, is analyzed with the Lagrangian atomic stress method, of which the reliability is also analyzed. The stress wave velocity corresponds to the speed of sound in the material, which is not dependent on the thermal shock energy. The peak pressure of a normal stress wave increases with the increase of thermal shock energy. We analyze the temperature evolution of the thermal shock region according to the Fourier transformation. It can be seen that the “obvious” velocity of heat propagation is less than the velocity of the stress wave; further, that the thermo-elastic stress wave may contribute little to the transport of kinetic energy. The heat propagation can be described properly by the heat conduction equation. These results may be useful for understanding the process of the thermal shock of tungsten.

Analysis of the stress intensity factor of welded joints under prior corrosion and complex stress fields

Taking cruciform welded joints under different corrosion degrees, together with 45° inclined angle and full penetration load-carrying fillet, as the research object, the stress intensity factor of a quarter-circular comer crack of welded joints is calculated based on FRANC3D with ABAQUS. Effects of different corrosion pit sizes, crack aspect ratios, and crack depths on stress intensity factor are analyzed. The results show that pit depth plays a major role in stress intensity factor, while the effect of pit radius is relatively small. The cracking modes of the surface and the deepest point are mode I, and mixed modes I and II, respectively. Effects of pit depths, crack aspect ratios, and crack depths on the stress intensity factor at the surface point are greater than at the deepest point.

Prediction of rock-burst-threatened areas in an island coal face and its prevention:A case study

The island coal face arises in coal mines with the purpose of preventing gas explosion or maintaining the balance between mining and tunneling. However, its particular stress conditions in the surrounding rock may increase the difficulty of stress control in the coal face and in its mining roadways, especially when the coal seam, the roof, and the floor have rock-burst propensities, The high energy accumulated in the island coal face and in its roof and floor will intensify rock-burst propensity or even induce rock burst, which further result in great casualties and financial losses. Taking island coal face 2321 in Jinqiao coal mine as a case, we propose a method for the prediction of rock-burst-threatened areas in an island coal face with weak rock-burst propensity. Based on the anaHysis of the movement of the overlying roof and characteristics of stress distribution, this method combined numerical simulation with drilling bits to ensure the prediction accuracy. The effects of coal pillars with different widths on the mitigation of stress concentration in the coal face and on the prevention of rock burst are analyzed together with the mech- anism behind. Finally, corresponding measures against the rock burst in the island coal face are proposed.

Technology of back stoping from level floors in gateway and pillar mining areas of extra-thick seams

According to the special requirements of secondary mining of resources in gateway-and-pillar goal in extra-thick seams of Shanxi, this paper presents a technical proposal of back stoping from level floors. Numerical simulation and theoretical analysis are ccsed to investigate the compaction characteristics of cavities under stress as well as an appropriate mining height of the primary-mining layer based on dif- ferent mining widths and pillar widths. For Yangjian coal mine, the mining thickness of the first seam during back stoping from level floor is determined to be 3 m, which meets the relevant requirements. Gateway-and-pillar goaf of a single layer has a range of influence of 9 m vertically. If gateway-and-pillar goaf occurs both in 9-1 and 9-5 layers, the range is extended to within 11.2 m. When the mining width of a gateway is less than 2 m or larger than 5 m, the gateway-and-pillar goal in the upper layer of the primary-mining seam can be filled in and compacted after stoping. When the working face is 2 m away from the gateway and pillar before entering into it and after passing through it, the coal body under the gateway and pillar is subjected to relatively high stress. During mining of the upper layer, moreover, the working face should interlock the goaf in primary-mining layer for 20 m.

Movement and deformation laws of the overlying strata in paste filling stope

We combined the similar simulation with numerical simulation to analyze the movement and deforma- tion features of overlying strata caused by paste backfill mining, study the movement and deformation laws of the overlying strata in paste backfill mining, structural movement of the stope strata as well as the stope stress distribution laws. Furthermore, authors also explored the key factors to the movement and deformation of the overlying strata in paste backfill mining. The results indicate that a caving zone existed in the bending zone only in the overlying strata of the paste backfill mining. Compared with the roof caving mining, the degree of stress concentration and area of influence in the paste filling stope were apparently smaller. And the degree of destruction and area of the overlying strata decreased prominently. Also, there was no apparent strata behavior in the working face. Lastly, the filling ratio was the key to control the movement and deformation of the overlying strata. Combined with a specific engineering example, the author proved the reliability of the simulation results and provided a theoretical basis for the further extension of the paste backfill mining.

The estimation of yearly probability gain for seismic statistical model

Based on the calculation method of information gain in the stochastic process presented by Vere-Jones, the rela tion between information gain and probability gain is studied, which is very common in earthquake prediction, and the yearly probability gain for seismic statistical model is proposed. The method is applied to the non stationary Poisson model with whole-process exponential increase and stress release model. In addition, the prediction method of stress release model is obtained based on the inverse function simulation method of stochastic variable.

Combination of direct-forcing fictitious domain method and sharp interface method for dielectrophoresis of particles

In this paper, we combine the direct-forcing fictitious domain （DF/FD） method and the sharp interface method to resolve the problem of particle dielectrophoresis in two dimensions. The flow field and the motion of particles are solved with the DF/FD method, the electric field is solved with the sharp inter- face method, and the electrostatic force on the particles is computed using the Maxwell stress tensor method. The proposed method is validated via three problems： effective conductivity of particle compos- ite between two planar plates, cell trapping in a channel, and motion of particles due to both conventional and traveling wave dielectrophoretic forces.