Numerical simulation of residual stress field in green power metallurgy compacts by modified Drucker-Prager Cap model

Compaction process simulation and residual stress prediction of green PM compact were carried out with elasto-plastic 3D FEA based on the modified Drueker-Prager Cap model in Abaqus. The model parameters of the investigated powder Distaloy AE were determined as functions of relative density through typical mechanical property tests of powder. The model was implemented as a user subroutine USDFLD. Single sided compaction of a d20 ram^5 mm disk green compact of Distaloy AE was simulated, and the residual stress of the disk after ejection was predicted with FEA. The FEA results of the compaction process and the residual stress of the disk show good agreement with compaction experiments and X-ray diffraction measurements, which validates the model and its parameters. The results indicate that the compressive residual stresses exist mainly in a thin layer on the side surface, but the residual stresses are very small on the top and bottom surfaces.

基于ABAQUS的螺纹钢树脂锚杆在软岩大变形隧洞中的应用研究

某输水工程隧洞围岩完整性较差,采用锚杆可以有效提高输水隧洞稳定性。文章利用ABAQUS软件对螺纹钢树脂锚杆在输水隧洞中的支护开展研究。结果表明,隧洞竖向位移大于水平位移,该处应为工程监测的重点部位;相比于初期锚杆支护方案,螺纹钢树脂锚杆支护方案的拱顶竖直位移较小(13.28mm),具有良好的支护效果;隧洞竖向位移中,拱底位移达到18.79mm,应采用加筋的方式防止拱底起鼓。研究结果为隧洞加固提供指导。

Numerical Simulation of the Whole Three-Dimensional Flow in a Stirred Tank with Anisotropic Algebraic Stress Model

In accordance to the anisotropic feature of turbulent flow, ananisotropic algebraic stress model is adopted to predict theturbulent flow field and turbulent characteristics generated by aRushton disc turbine with the improved inner-outer iterativeprocedure. The predicted turbulent flow is compared with experimentaldata and the simulation by the standard k-ε turbulence model. Theanisotropic algebraic stress model is found to give better predictionthan the standard k-ε turbulence model. The predicted turbulent flowfield is in accordance to experimental data and the trend of theturbulence intensity can be effectively reflected in the simulation.

Evaluation of coal pillar loads during longwall extraction using the numerical method and its application

It is very difficult to reasonably evaluate the loads acting on coal pillars in longwall panels during the planning of a new pillar system. The application of empirical equations is a common practice in calculating coal pillar loads while designing a new pillar. This paper proposes numerical models for evaluating coal pillar loads. The key of building a successful numerical model for calculating coal pillar loads lies in the fact that the model should represent the redistribution of stress all over the longwall panels and the surrounding areas, and it is especially important to include the characteristics of the stress rebuilding process in the gob areas, which are crucial for the building process of coal pillar loads. Based on the geo-mechanical background of the Baoshan Coal Mine, this paper details the procedures of applying numerical models to the evaluation of coal pillar loads and their local practices. The study results show it is feasible and reasonable to use numerical models to evaluate coal pillar loads.

Prediction of plastic zone size around circular tunnels in non-hydrostatic stress field

This paper discusses the calculation of plastic zone properties around circular tunnels to rock-masses that satisfy the Hoek–Brown failure criterion in non-hydrostatic condition,and reviews the calculation of plastic zone and displacement,and the basis of the convergence–confinement method in hydrostatic condition.A two-dimensional numerical simulation model was developed to gain understanding of the plastic zone shape.Plastic zone radius in any angles around the tunnel is analyzed and measured,using different values of overburden(four states)and stress ratio(nine states).Plastic zone radius equations were obtained from fitting curve to data which are dependent on the values of stress ratio,angle and plastic zone radius in hydrostatic condition.Finally validation of this equation indicate that results predict the real plastic zone radius appropriately.

Effective grouting area of jointed slope and stress deformation responses by numerical analysis with FLAC^(3D)

To study the grouting reinforcement mechanism in jointed rock slope, first, the theoretical deduction was done to calculate the critical length of slipping if the slope angle is larger than that of joint inclination; Second, the numerical calculation model was founded by FLAG^3D, so as to find the stress and deformation responses of rock mass in the state before and after grouting, the analysis results show that the range between the boundary of critical slipping block and the joint plane that passes the slope toe is the effective grouting area （EGA）. After excavation, large deformation occurs along the joint plane. After grouting, the displacements of rock particles become uniform and continuous, and large deformations along the joint plane are controlled; the dynamic displacement can re- flect the deformation response of slope during excavation in the state before and after grouting, as well as the shear location of potential slip plane. After grouting, the dynamic displacement of each monitoring point reaches the peak value with very few time steps, which indicate that the parameters of the joint plane, such as strength and stiffness, are improved; the stress field becomes uniform. Tensile area reduces gradually; whole stability of the slope and its ability to resist tensile and shear stress are improved greatly.

Prevention of rockburst by guide holes based on numerical simulations

We studied variations in the stress field around guide holes drilled during tunnel excavation to understand the mechanical mechanism by which these holes help prevent rockburst.The study used elasto-plastic analysis of a circular chamber under non-axisymmetrical loads.The results showed that the unloading of in-situ stresses, and the forming of a secondary stress field, leads to a severe change in the stress field around the guide holes.This causes the formation of an X-shaped area of plastic deformation, which prevents the rockburst.Adopting a sub-model finite element technique, we analyzed the factors that influence the distribution of the plastic area, such as the guide hole distribution and the in-situ stress state.The calculations showed that higher initial stresses result in greater adjustment to the stress field.When the stress concentration is greater the size of the plastic area surrounding the guide hole is larger.A multi-row distribution of the guide holes shaped like a quincunx can increase the interconnectivity of the plastic areas and allow the plastic area to extend from the tunnel wall deep into the surrounding rock.An optimized design was put forward based on the distribution of the plastic area around guide holes and the factors that influence it.

Application of preblasting to high-section top coal caving for steepthick coal sea

For mining extra-steep-thick coal seam, the sublevel top coal caving is a high efficient method in practical engineering. However, major challenges associated with mining high-section top-coal-caving （HSTCC） are related to the resulting high ground stresses. Inevitability, using the high-section sublevel top coal caving for extra-steep-thick coal seam, the large scale of mined-out area appears. If the prefracture blasting and hydraulic fracture techniques are utilized, the top coal damage and cracks will develop, and the mining complexity will increase, such as seam inclination, continuity, mechanical characteristics of roof and susceptibility of top coal, etc. First, the field conditions of B1＋2 seam were investigated at the ＋588 level of the Weihuliang Underground Mine of China. Subsequently, according to caving mechanism of strata response obtained from several special models including physical simulation tests and numerical simulation models, the prefracture process including blasting and injecting water were analyzed. Then, the prefracture blasting technique was successfully applied to the caving of 52 m-sublevel seam. Finally, the effects were verified by advanced detecting instruments, and the results show these methods and measurements are feasible and valid.

A case study of multi-seam coal mine entry stability analysis with strength reduction method

In this paper,the advantage of using numerical models with the strength reduction method(SRM) to evaluate entry stability in complex multiple-seam conditions is demonstrated.A coal mine under variable topography from the Central Appalachian region is used as a case study.At this mine,unexpected roof conditions were encountered during development below previously mined panels.Stress mapping and observation of ground conditions were used to quantify the success of entry support systems in three room-and-pillar panels.Numerical model analyses were initially conducted to estimate the stresses induced by the multiple-seam mining at the locations of the affected entries.The SRM was used to quantify the stability factor of the supported roof of the entries at selected locations.The SRM-calculated stability factors were compared with observations made during the site visits,and the results demonstrate that the SRM adequately identifies the unexpected roof conditions in this complex case.It is concluded that the SRM can be used to effectively evaluate the likely success of roof supports and the stability condition of entries in coal mines.

湖南印支期过铝质花岗岩的形成:岩浆底侵与地壳加厚热效应的数值模拟

在合理构建华南印支期地质-物理模型的基础上,利用FLAC软件,模拟了该区印支期过铝质富钾花岗岩形成与基性岩浆底侵,及陆壳变形叠置加厚两种动力学背景的可能联系.模拟结果表明,220 Ma±的基性岩浆底侵能导致地壳含水矿物相岩石的深熔,但除非印支期存在大规模基性岩浆的底侵作用,否则其热效应持续时间和热效应波及范围难以形成具大岩基规模的湖南印支期花岗岩.在陆壳叠置加厚模型中,地壳的叠置加厚可导致中下地壳界面温度升高到700℃以上,引起片麻质岩石熔融,当加厚因子达1.3,白云母矿物脱水熔融产生的熔体达到熔体流动临界比例（≥20％）,从而形成花岗岩基;结合印支期挤压逆冲推覆构造和同期基性火山岩极少出露的地质事实,认为陆壳变形加厚可能是湖南印支期构造岩浆作用形成的主导因素.

Investigation of Barree-Conway non-Darcy flow effects on coalbed methane production

Coalbed gas non-Darcy flow has been observed in high permeable fracture systems,and some mathematical and numerical models have been proposed to study the effects of non-Darcy flow using Forchheimer non-Darcy model.However,experimental results show that the assumption of a constant Forchheimer factor may cause some limitations in using Forchheimer model to describe non-Darcy flow in porous media.In order to investigate the effects of non-Darcy flow on coalbed methane production,this work presents a more general coalbed gas non-Darcy flow model according to Barree-Conway equation,which could describe the entire range of relationships between flow velocity and pressure gradient from low to high flow velocity.An expanded mixed finite element method is introduced to solve the coalbed gas non-Darcy flow model,in which the gas pressure and velocity can be approximated simultaneously.Error estimate results indicate that pressure and velocity could achieve first-order convergence rate.Non-Darcy simulation results indicate that the non-Darcy effect is significant in the zone near the wellbore,and with the distance from the wellbore increasing,the non-Darcy effect becomes weak gradually.From simulation results,we have also found that the non-Darcy effect is more significant at a lower bottom-hole pressure,and the gas production from non-Darcy flow is lower than the production from Darcy flow under the same permeable condition.

Numerical Simulations of Upper Plenum Thermal-Hydraulics of Monju Reactor Vessel Using High Resolution Mesh Models

In order to evaluate the effects of the FHs （flow holes） on the inner barrel, which were installed in the upper plenum of the Monju reactor vessel, a high resolution meshes around the FHs was constructed. Using this model, it was mainly clear that in the 40% rated operational conditions, the shape of the FHs on the inner barrel did not change largely to the upper plenum thermal-hydraulics. The effect of the FHs on the honeycomb structure in the upper structure was also investigated in these calculations. The results indicated that the height of thermal stratification interface became lower than that evaluated from the test data.

Impacts of wind stress on saltwater intrusion in the Yangtze Estuary

The observation at the Chongxi gauging station indicated the salinity of saltwater spilling over from the North Branch to the South Branch increased abnormally from November 10 to 12 in 2009 (during neap tide) and from February 11 to 12 in 2010 (during moderate tide).We found for the first time that the strong northerly wind was responsible for the above abnormal salinity increase.Previous studies indicated that the saltwater intrusion in the Yangtze Estuary is influenced mainly by the river discharge,the tide,and the wind stress,but the impacts of variations of wind speed and direction on it have not been investigated.In this study the impacts of wind stress on the saltwater intrusion were numerically simulated and the associated mechanisms were analyzed.The model results were consistent with the observed data obtained at six gauging stations during February and March in 2007 and four gauging stations in March 2008,and the abnormal salinity risings were well captured.Meanwhile,if the wind speed is reduced by half,the salinity there will be significantly decreased.Driven by the monthly mean river discharge of 11000 m 3 /s and northerly wind of 5 m/s from January to February,the model simulated the temporal and spatial variation of saltwater intrusion.The wind-driven circulation,as well as the net water and salt fluxes from the North Branch into the South Branch,was calculated and analyzed in the cases of different wind speeds and directions.The results indicated that the intensity of the saltwater intrusion in the Yangtze Estuary is significantly influenced by the wind speeds and directions.

Investigation on Combustion Characteristics and NO Formation of Methane with Swirling and Non-Swirling High Temperature Air

Combustion characteristics of methane jet flames in an industrial burner working in high temperature combustion regime were investigated experimentally and numerically to clarify the effects of swirling high temperature air on combustion.Speziale-Sarkar-Gatski(SSG) Reynolds stress model,Eddy-Dissipation Model(EDM),Discrete Ordinates Method(DTM) combined with Weighted-Sum-of-Grey Gases Model(WSGG) were employed for the numerical simulation.Both Thermal-NO and Prompt-NO mechanism were considered to evaluate the NO formation.Temperature distribution,NO emissions by experiment and computation in swirling and non-swirling patterns show combustion characteristics of methane jet flames are totally different.Non-swirling high temperature air made high NO formation while significant NO prohibition were achieved by swirling high temperature air.Furthermore,velocity fields,dimensionless major species mole fraction distributions and Thermal-NO molar reaction rate profiles by computation interpret an inner exhaust gas recirculation formed in the combustion zone in swirling case.

Assessment of alternative scale-providing variables in a Reynolds-stress model using high-order methods

This paper is set in the high-order finite-difference discretization of the Reynolds-averaged Navier-Stokes(RANS)equations,which are coupled with the turbulence model equations.Three alternative scale-providing variables for the specific dissipation rate(o)are implemented in the framework of the Reynolds stress model(RSM)for improving its robustness.Specifically,g(=1/√ω)has natural boundary conditions and reduced spatial gradients,and a new numerical constraint is imposed on itω(=lnω)can preserve positivity and also has reduced spatial gradients;the eddy viscosity v,also has natural boundary conditions and its equation is improved in this work.The solution polynomials of the mean-flow and turbulence-model equations are both reconstructed by the weighted compact nonlinear scheme(WCNS).Moreover,several numerical techniques are introduced to improve the numerical stability of the equation system.A range of canonical as well as industrial turbulent flows are simulated to assess the accuracy and robustness of the scale-transformed models.Numerical results show that the scale-transformed models have significantly improved robustness compared to the w model and still keep the characteristics of RSM.Therefore,the high-order discretization of the RANS and RSM equations,which number 12 in total,can be successfully achieved.

Meandering Spiral Waves Induced by Time-Periodic Coupling Strength

Effects of time-periodic coupling strength （TPCS） on spiral waves dynamics are studied by numerical computations and mathematical analyses. We find that meandering or drifting spirals waves, which are not observed for the case of constant coupling strength, can be induced by TPCS. In particular, a transition between outward petal and inward petal meandering spirals is observed when the period of TPCS is varied. These two types of meandering spirals are separated by a drifting spiral, which can be induced by TPCS when the period of TPCS is very close to that of rigidly rotating spiral. Similar results can be obtained if the coupling strength is modulated by a rectangle wave. Furthermore, a kinetic model for spiral movement suggested by Diet al., [Phys. Rev. E 85 （2012） 046216] is applied for explaining the above findings. The theoretical results are in good qualitative agreement with numerical simulations.