Spatial deduction of mining-induced stress redistribution using an optimized non-negative matrix factorization model

Investigation of mining-induced stress is essential for the safety of coal production.Although the field monitoring and numerical simulation play a significant role in obtaining the structural mechanical behaviors,the range of monitoring is not sufficient due to the limits of monitoring points and the associated numerical result is not accurate.In this study,we aim to present a spatial deduction model to characterize the mining-induced stress distribution using machine learning algorithm on limited monitoring data.First,the framework of the spatial deduction model is developed on the basis of non-negative matrix factorization(NMF)algorithm and optimized by mechanical mechanism.In this framework,the spatial correlation of stress response is captured from numerical results,and the learned correlation is employed in NMF as a mechanical constrain to augment the limited monitoring data and obtain the overall mechanical performances.Then,the developed model is applied to a coal mine in Shandong,China.Experimental results show the stress distribution in one plane is derived by several monitoring points,where mining induced stress release is observed in goaf and stress concentration in coal pillar,and the intersection point between goaf and coal seam is a sensitive area.The indicators used to evaluate the property of the presented model indicate that 83%mechanical performances have been captured and the deduction accuracy is about 92.9%.Therefore,it is likely that the presented deduction model is reliable.

Numerical simulation of spatial distributions of mining-induced stress and fracture fields for three coal mining layouts

In this study, the spatial distributions of stress and fracture fields for three typical underground coal mining layouts, Le, non-pillar mining （NM）, top-coal caving mining （TCM） and protective coal-seam mining （PCM）, are modeled using discrete element software UDEC, The numerical results show that different mining layouts can lead to different mining-induced stress fields, resulting in diverse fracture fields, For the PCM, the mining influenced area in front of the mining faces is the largest, and the stress concentration factor in front of the mining faces is the lowest, The spatial shapes of the mining-induced fracture fields under NM, TCM and PCM differ, and they are characterized by trapezoidal, triangular and tower shapes, respectively, The fractal dimensions of mining-induced fractures of the three mining layouts decrease in the order of PCM, TCM and NM, It is also shown that the PCM can result in a better gas control effect in coal mines with high outburst potential, The numerical results are expected to provide a basis for understanding of mining-induced gas seepage fields and provide a reference for high- efficiency coal mining,

Stress distribution in roadbeds of slab tracks with longitudinal discontinuities

Stress concentration occurs in the foundations of railway tracks where discontinuous components are located.The exacerbated stress under the expansion joints in slab tracks may trigger foundation failures such as mud pumping.Although the higher stress due to the discontinuities of track structures has been discussed in past studies,few focused on the stress response of roadbeds in slab tracks and quantitatively characterized the stress pattern.In this paper,we performed a dynamic finite element analysis of a track-formation system,incorporating expansion joints as primary longitudinal discontinuities.The configurations of CRTS Ⅲ slab tracks and the contact conditions between concrete layers were considered.Numerical results show that longitudinal influencing length of induced stress on roadbed under wheel load relates to the contact conditions between concrete layers,increasing nonlinearly at a larger coefficient of friction.Given a measured coefficient of friction of 0.7,the calculated longitudinal influencing length(9.0 m) matches with field data.The longitudinal influencing length is not affected with the increasing velocity.As stress concentration arises with expansion joints,the worstcase scenario emerges when double-axle loads are exerted immediately above the expansion joints between concrete bases.A stress concentration factor Cvon the roadbed is proposed;it increases with the increasing velocity,changing from 1.33 to 1.52 at velocities between 5 and 400 km/h.The stress distribution on roadbeds transforms from a trapezoid pattern at continuous sections to a triangle pattern at points with longitudinal discontinuities.An explicit expression is finally proposed for the stress pattern on roadbed under expansion joints.Although structural discontinuities induce stress raiser,the extent of concentration is mitigated with increasing depth at different velocity levels.

Micro sliding friction model considering periodic variation stress distribution of contact surface and experimental verification

Micro sliding phenomenon widely exists in the operation process of mechanical systems,and the micro sliding friction mechanism is always a research hotspot.In this work,based on the total reflection method,a measuring device for interface contact behavior under two-dimensional(2D)vibration is built.The stress distribution is characterized by the light intensity distribution of the contact image,and the interface contact behavior in the 2D vibration process is studied.It is found that the vibration angle of the normal direction of the contact surface and its fluctuation affect the interface friction coefficient,the tangential stiffness,and the fluctuation amplitude of the stress distribution.Then they will affect the change of friction state and energy dissipation in the process of micro sliding.Further,an improved micro sliding friction model is proposed based on the experimental analysis,with the nonlinear change of contact parameters caused by the normal contact stress distribution fluctuation taken into account.This model considers the interface tangential stiffness fluctuation,friction coefficient hysteresis,and stress distribution fluctuation,whose simulation results are consistent well with the experimental results.It is found that considering the nonlinear effect of a certain contact parameter alone may bring a greater error to the prediction of friction behavior.Only by integrating multiple contact parameters can the accuracy of friction prediction is improved.

Stiffness and Shear Stress Distribution of Glulam Beams in Elastic-Plastic Stage:Theory,Experiments and Numerical Modelling

Traditional methods focus on the ultimate bending moment of glulam beams and the fracture failure of materials with defects,which usually depends on empirical parameters.There is no systematic theoretical method to predict the stiffness and shear distribution of glulam beams in elastic-plastic stage,and consequently,the failure of such glulam beams cannot be predicted effectively.To address these issues,an analytical method considering material nonlinearity was proposed for glulam beams,and the calculating equations of deflection and shear stress distribution for different failure modes were established.The proposed method was verified by experiments and numerical models under the corresponding conditions.Results showed that the theoretical calculations were in good agreement with experimental and numerical results,indicating that the equations proposed in this paper were reliable and accurate for such glulam beams with wood material in the elastic-plastic stage ignoring the influence of mechanic properties in radial and tangential directions of wood.Furthermore,the experimental results reported by the previous studies indicated that the method was applicable and could be used as a theoretical reference for predicting the failure of glulam beams.

Effects of Material Parameters on Stress Distribution in Casing-cement-formation(CCF)Multilayer Composite System

This work focus on the stress distribution of the casing-cement-formation(CCF)multilayer composite system,which is a borehole system with multiple casings and cement sheathes.Mostof the previous relevant studies are based on the traditional CCF system with the single casing and cement sheath,but these results are not adaptive to the CCF system multiple composite system.In this paper,the FEM numerical model of CCF multilayer composite system was constructed.Numerical simulations were calculated and compared with the system which consists of the single casing and cement sheath.Results show that the multilayer composite system possesses better performance.On this basis,the sensitivity analysis of main influence mechanical parameters such as in-situ stress,the elastic of cement sheathes and the elastic of formation are conducted.The cement sheath on the inside,namely cement sheath-1,is sensitive to its elastic modulus;meanwhile,the cement sheath on the outside,namely cement sheath-2,is not so sensitive to the elastic modulus of cement sheath-1.Cement sheath-1 and cement sheath-2 are all sensitive to the elastic modulus of cement sheath-2,and the mises stress of them has opposite trend to the elastic modulus of cement sheath-2.The proper values of elastic modulus of cement sheath-1 and cement sheath-2 are 5GPa and 5GPa to 30GPa,respectively.Under the in-situ stress ratio σh/σH=0.7,the maximum mises stress of cementsheath-1 and cement sheath-2 increase as the increase of σh,and they are nearly equal when σh=15GPa.This research can be helpful for the design and analysis of CCF multilayer composite system.

Finite element analysis of stress distribution and burst failure of SiC_f/Ti-6Al-4V composite ring

A three-dimensional cyclic symmetry finite element model of titanium-matrix composites（TMCs） ring was developed to investigate the stress distribution and burst failure. The effects of fiber volume fractions, reinforced areas, thermal residual stresses and two different temperatures on stress distribution were studied. The burst speed was obtained through analyzing the hoop tensile stresses under a series of rotating speeds. The results indicate that at the two different temperatures, the influences of fiber volume fractions and reinforced areas on stress level and distribution are different. Some proposals are provided for the structure design of the TMCs ring. With regard to thermal residual stresses, a larger reinforced area is an advisable choice for design of the ring at higher temperature.

Deformation characteristics and reinforcement technology for entry subjected to mining-induced stresses

The entry at Zhangcun coal mine in Lu＇an coal mining area in Shanxi Province suffered from severe mining-induced stresses with the heading face driven oppositely to an adjacent working face. In this paper, the characteristics of deformation and failure of the entry were investigated in terms of the tempo-spatial relations between heading and working faces through field study and numerical modeling. The three-dimensional （3D） finite difference models were built to investigate stresses, displacements and damages in the surrounding rocks of the entry and the working face. The field study includes selection of reinforcing methods and materials, design parameters, and determination of cable prestress. The monitoring data of entry deformation and stress along the cables during every stage were presented. The state of the reinforced entry was evaluated based on the monitoring data. The results demonstrate that before the heading face of the entry crosses the adjacent working face, the influence of advanced abutment pressure caused by adjacent working face upon the entry is not significant. After they cross each other, however, the lateral abutment pressure will have an evident impact on the entry. The displacement rate of the entry will be greatly increased and reaches a certain value within a certain distance between the heading face and the working face. Then, it will increase again with the presence of secondary mining-induced pressure on the entry when the present working face advances. The fully-grouted cable with short length, high strength and high prestress is an effective way to reinforce the entry suffering from severe mining-induced stresses, which greatly reduces the displacement and failure possibility of the entry. Finally, the principles and recommendations for reinforcing design of entries suffering from severe mining-induced stresses were proposed according to field study, numerical modeling and experiences from other coal mines. Problems encountered in field study and suggestions for reinforcement were also discussed.

Nonlinear Analysis of Axial-load and Stress Distribution for Threaded Connection

Analytical method for the distributions of axial-load and stress is based on elastic assumption, but the threaded connections are often in plastic deformation stage in practice. Meanwhile the strain in the threaded connection is difficult to measure. So it is necessary to study the reliable numerical method. At present neither the convergence analysis of the computational results nor the elastic-plastic analysis in the loading-unloading process are studied. In this paper, von Mises plasticity and kinematic hardening model is used to describe the material response. A new convergence criterion for nonlinear finite element analysis of the loading-unloading process is proposed. An axisymmetric finite element model according to the proposed convergence criterion is developed and used to analyze the distributions of axial-load and stress. It can be conclude that the stress distribution analysis is more dependent on the mesh density than the axial-load distribution analysis. The stress distribution result indicates that with increasing of applied load, the engaged threads close to the nut-bearing surface become plastic firstly. The axial-load distribution result reveals that the load percentage carried by single thread depends on the position of thread and load intensity. When the load is relatively small, the applied load is mainly carried by the engaged threads near the nut-bearing surface, when the load is larger, the differences of percentages for all threads become small. The proposed convergence analyzing procedure is applicable for other nonlinear analyses. The obtained distributions of axial-load and stress can be a reference of engineering application.

Fault-Induced Coal Burst Mechanism under Mining-Induced Static and Dynamic Stresses

Fault is a common geological structure that has been revealed in the process of underground coal excavation and mining.The nature of its discontinuous structure controls the deformation,damage,and mechanics of the coal or rock mass.The interaction between this discontinuous structure and mining activities is a key factor that dominates fault reactivation and the coal burst it can induce.This paper first summarizes investigations into the relationships between coal mining layouts and fault occurrences,along with relevant conceptual models for fault reactivation.Subsequently,it proposes mechanisms of fault reactivation and its induced coal burst based on the superposition of static and dynamic stresses,which include two kinds of fault reactivations from:mining-induced quasi-static stress(FRMSS)-dominated and seismic-based dynamic stress(FRSDS)-dominated.These two kinds of fault reactivations are then validated by the results of experimental investigations,numerical modeling,and in situ microseismic monitoring.On this basis,monitoring methods and prevention strategies for fault-induced coal burst are discussed and recommended.The results show that fault-induced coal burst is triggered by the superposition of high static stress in the fault pillar and dynamic stress from fault reactivation.High static stress comes from the interaction of the fault and the roof structure,and dynamic stress can be ascribed to FRMSS and FRSDS.The results in this paper could be of great significance in guiding the monitoring and prevention of fault-induced coal bursts.

Stress distribution and surrounding rock control of mining near to the overlying coal pillar in the working face

The occurrence of overlying coal pillar(OCP)exerts a strong effect on the stress and strain distribution of the surrounding rock in the stope.In this paper,the stress distribution characteristics are analyzed via the numerical calculation with the account of OCP presence or absence.In addition,this study revealed the joint effect of side pressure relief area of the goaf and stress concentration in OCP on the final stress distribution.Furthermore,the rules of abutment stress distribution affected by three influencing factors,namely horizontal-vertical distances between OCP and working face and buried depth of OCP,are analyzed.The functional model linking the peak stress of surrounding rock with the above influencing factors is developed.The field application of the above results proved that the rib spalling and deformation of a 2.95 m-high and 5.66 m-wide roadway could be efficiently controlled by rationally adjusting working states of the support,and adopting the hydraulic prop coordinated with the p type metal beam and anchor cable to strengthen the surrounding rock of working face and roadway,respectively.The proposed measures are considered appropriate to satisfy the safe operation requirements.

Measurements of in situ stress and mining-induced stress in Beiminghe Iron Mine of China

In order to obtain the distribution rules of in situ stress and mining-induced stress of Beiminghe Iron Mine, the stress relief method by overcoring was used to measure the in situ stress, and the MC type bore-hole stress gauge was adopted to measure the mining-induced stress. In the in situ stress measuring, the technique of improved hollow inclusion cells was adopted, which can realize complete temperature compensation. Based on the measuring results, the distribution model of in situ stress was established and analyzed. The in situ stress measuring result shows that the maximum horizontal stress is 1.75-2.45 times of vertical stress and almost 1.83 times of the minimum horizontal stress in this mineral field. And the mining-induced stress measuring result shows that, according to the magnitude of front abutment pressure the stress region can be separated into stress-relaxed area, stress- concentrated area and initial stress area. At the -50 m mining level of this mine, the range of stress-relaxed area is 0-3 m before mining face; the range of stress-concentrated area is 3-55 m before mining face, and the maximum mining-induced stress is 16.5-17.5 MPa, which is 15-20 m from the mining face. The coefficient of stress concentration is 1.85.

A Modification to Vertical Distribution of Tidal Flow Reynolds Stress in Shallow Sea

Tidal flow is a periodic movement of unsteady and non-uniform, which has acceleration and deceleration process obviously, especially in coastal shallow waters. Many researches show that vertical distribution of tidal flow Reynolds stress deviated from linear distribution. The parabolic distribution of the tidal flow Reynolds stress was proposed by Song et al. （2009）. Although the model fills better with field observations and indoor experimental data, it has the lower truncated series expansion of tidal flow Reynolds stress, and the description of the distribution is not very comprehensive By introducing the motion equation of tidal flow and improving the parabolic distribution established by Song et al. （2009）, the cubic distribution of the tidal flow Reynolds stress is proposed. The cubic distribution is verified well by field data （Bowden and Fairbairn, 1952; Bowden et al., 1959; Rippeth et al., 2002） and experimental data （Anwar and Atkins, 1980）, is consistent with the numerical model results of Kuo et al. （1996）, and is compared with the parabolic distribution of the tidal flow Reynolds stress. It is shown that this cubic distribution is not only better than the parabolic distribution, but also can better reflect the basic features of Reynolds stress deviating from linear distribution downward with the tidal flow acceleration and upward with the tidal flow deceleration, for the foundation of further study on the velocity profile of tidal flow.

Study of Stress and Strain Distributions of First Pass Conventional Spinning Under Different Roller-trace

Based on simplified axisymmetrical forming model, a elasto-plastic FEM simulation system of multi-pass conventional spinning is developed. Taking the typical draw-spinning as the study object, and establishing reasonable mechanics model, research on the first pass of spinning process is carried out with FEM system developed. The distributions of the stress and strain are obtained by three types of roller-trace curves: straight line, involute curves and quadratic curves. The results are as follows: (1) The values of equivalent stress and strain are the lowest under involute curve compared to other two curves, and they change relatively small and decrease with the increase of radius. The values of equivalent stress and strain is the highest under quadratic curves, and increase with the increase of radius. (2) The value of radial stress is smallest under involute curve, and is the largest under straight line. Value of radial stress is often used as the criterion of cracking limit, so its distribution laws can provide references for studying the condition of cracking in multi-pass conventional spinning under different roller-trace. (3) Tangential stress is compressive stress. Absolute value of tangential stress is the smallest under involute curve, and values of tangential stress are close between other two curves. The distribution laws of tangential stress can serve as a significant guide to research the critical condition of wrinkling in multi-pass conventional spinning under different roller-trace. (4) The reduction of thickness is the smallest under involute curve. The distribution of the thickness strain is very unequal under quadratic curves. The results obtained can provide references for selecting reasonable roller-trace in multi-pass conventional spinning.

Stress distribution rule of roadway affected by overhead mining in gently inclined coal seams group

In light of the severe deformation and destruction of the district raise tunnel in the mining area at the northern part of the Lubanshan colliery, by the theoretic analysis and numerical simulation, both the mining stress distribution in seams group and the deformation and destruction mechanism of floor district raise were investigated. The results show that, at the maximum vertical distance of 40 m, the abutment stress has an influence on the recovery of 2# and 3# coal seam and 8# coal seam at distance of 30 m. As a result, the recovery of 8# is rather than those of 2# or 3# coal seam, which contributes to the deformation and destruction of the district raise surrounding rock. The major factors affecting the abutment stress include the mining depth, mining height, residual gob space, adjacent working faces and short spacing coal seam recovery.

Characteristics of evolution of mining-induced stress field in the longwall panel:insights from physical modeling

The evolution of mining-induced stress field in longwall panel is closely related to the fracture field and the breaking characteristics of strata.Few laboratory experiments have been conducted to investigate the stress field.This study investigated its evolution by constructing a large-scale physical model according to the in situ conditions of the longwall panel.Theoretical analysis was used to reveal the mechanism of stress distribution in the overburden.The modelling results showed that:(1)The major principal stress field is arch-shaped,and the strata overlying both the solid zones and gob constitute a series of coordinated load-bearing structures.The stress increasing zone is like a macro stress arch.High stress is especially concentrated on both shoulders of the arch-shaped structure.The stress concentration of the solid zone in front of the gob is higher than the rear solid zone.(2)The characteristics of the vertical stress field in different regions are significantly different.Stress decreases in the zone above the gob and increases in solid zones on both sides of it.The mechanical analysis show that for a given stratum,the trajectories of principal stress are arch-shaped or inverselyarched,referred to as the‘‘principal stress arch’’,irrespective of its initial breaking or periodic breaking,and determines the fracture morphology.That is,the trajectories of tensile principal stress are inversely arched before the first breaking of the strata,and cause the breaking lines to resemble an inverted funnel.In case of periodic breaking,the breaking line forms an obtuse angle with the advancing direction of the panel.Good agreement was obtained between the results of physical modeling and the theoretical analysis.

Thermal residual stresses and stress distributions under tensile and compressive loadings of short fiber reinforced metal matrix composites

The thermal residual stresses and the stress distributions of short fiber reinforced metal matrix composite under tensile and compressive loadings were studied using large strain axisymmetric elasto plastic finite element method. It is demonstrated that the thermal residual stresses can result in asymmetrical stress distributions and matrix plasticity. The thermal residual stresses decrease the stress transfer in tension and enhance the stress transfer in compression. The fiber volume fraction has more important effects on the thermal residual stresses and the stress distributions under tensile and compressive loadings than the fiber aspect ratio and the fiber end distance. [

Analytical model and application of stress distribution on mining coal floor

Given the analysis of underground pressure, a stress calculation model of coal floor stress has been established based on a theory of elasticity. The model presents the law of stress distribution on the relatively fixed position of the mining coal floor： the extent of stress variation in a fixed floor position decreases gradually along with depth, the decreasing rate of the vertical stress is clearly larger than that of the horizontal stress at a specific depth. The direction of the maximum principal stress changes gradually from a vertical direction to a horizontal direction with the advance of the working face. The deformation and permeability of the rock mass of the coal floor are obtained by contrasting the difference of the principal stress established from theoretical calculations with curves of stress-strain and permeability-strain from tests, which is an important mechanical basis for preventing water inrush from confined aauifers.

Analytical modeling and simulation of porous electrodes: Li-ion distribution and diffusion-induced stress

A new model of porous electrodes based on the Gibbs free energy is developed, in which lithium-ion(Liion) diffusion, diffusion-induced stress(DIS), Butler–Volmer(BV) reaction kinetics, and size polydispersity of electrode particles are considered. The influence of BV reaction kinetics and concentration-dependent exchange current density(ECD) on concentration profile and DIS evolution are numerically investigated. BV reaction kinetics leads to a decrease in Li-ion concentration and DIS. In addition, concentrationdependent ECD results in a decrease in Li-ion concentration and an increase in DIS. Size polydispersity of electrode particles significantly affects the concentration profile and DIS.Optimal macroscopic state of charge(SOC) should consider the influence of the microscopic SOC values and mass fractions of differently sized particles.

Simulation of Stress Distribution around Tunnels and Interaction between Tunnels Using an Elasto–plastic Model

This article presents a computer simulation of stress distribution around tunnels and interaction between tunnels using an elasto-plastic model. A finite element method using ANSYS software has been used for the analyses of one and two tunnels at different overburden depths with different separating distances between the tunnels. The results of numerical analyses indicate that stress distribution and stress concentration around the tunnels vary with the overburden depths. It is found that the coefficients of stress concentration for elasto-plastic medium are smaller than those for elastic one by 1.9%. Furthermore, the interaction between the two tunnels rapidly decreases with the increase of separation distance between them. In addition, for quantitatively describing the interaction between the two tunnels, a critical separation distance is introduced. The critical separation distances between the two tunnels at different overburden depths are 8 m, 12 m, and 14 m respectively. This fact is very important and essential for the design of mining tunnels and to ensure safety in tunnel engineering.