Based on the analysis method for tailings dam in upstream raising method presently used in metallurgy and nonferrous metals tailings depository in the world, an effective stress analysis method of seismic response for...Based on the analysis method for tailings dam in upstream raising method presently used in metallurgy and nonferrous metals tailings depository in the world, an effective stress analysis method of seismic response for high tailings dam was developed according to the results of engineering geological exploration, static and dynamic test and stability analysis on Baizhishan tailing dam 113.5 m high. The law of generation, diffusion and dissipation of seismic pore water pressure during and after earthquake was investigated, and the results of tailings dam’s acceleration, seismic dynamic stress and pore water pressure were obtained. The results show that the seismic stability and liquefaction resistance of high tailings dam are strengthened remarkably, and the scope and depth of liquefaction area at the top of dam are reduced greatly. The interior stress is compressive stress, the stress level of every element is less than 1.0 and the safety coefficient of every element is greater than 1.0. The safety coefficient against liquefaction of every element of tailing dam is greater than 1.5 according to the effective stress analysis of seismic response by finite element method. The calculated results prove that liquefaction is the main reason of seismic failure of high tailing dams, and the effect of seismic inertia forces on high tailing dams’ stability during earthquake is secondary reason.展开更多
The size-dependent effect on the biaxial and shear nonlinear buckling analysis of an isotropic and orthotropic micro-plate based on the surface stress, the modified couple stress theory （MCST）, and the nonlocal elas...The size-dependent effect on the biaxial and shear nonlinear buckling analysis of an isotropic and orthotropic micro-plate based on the surface stress, the modified couple stress theory （MCST）, and the nonlocal elasticity theories using the differential quadrature method （DQM） is presented. Main advantages of the MCST over the classical theory （CT） are the inclusion of the asymmetric couple stress tensor and the consideration of only one material length scale parameter. Based on the nonlinear von Karman assumption, the governing equations of equilibrium for the micro-classical plate consid- ering midplane displacements are derived based on the minimum principle of potential energy. Using the DQM, the biaxial and shear critical buckling loads of the micro-plate for various boundary conditions are obtained. Accuracy of the obtained results is validated by comparing the solutions with those reported in the literature. A parametric study is conducted to show the effects of the aspect ratio, the side-to-thickness ratio, Eringen＇s nonlocal parameter, the material length scale parameter, Young＇s modulus of the surface layer, the surface residual stress, the polymer matrix coefficients, and various boundary conditions on the dimensionless uniaxial, biaxial, and shear critical buckling loads. The results indicate that the critical buckling loads are strongly sensitive to Eringen＇s nonlocal parameter, the material length scale parameter, and the surface residual stress effects, while the effect of Young＇s modulus of the surface layer on the critical buckling load is negligible. Also, considering the size dependent effect causes the increase in the stiffness of the orthotropic micro-plate. The results show that the critical biaxial buckling load increases with an increase in G12/E2 and vice versa for E1/E2. It is shown that the nonlinear biaxial buckling ratio decreases as the aspect ratio increases and vice versa for the buckling amplitude. Because of the most lightweight micro-composite materials with high strength/weight and stiffness/weight ratios, it is anticipated that the results of the present work are useful in experimental characterization of the mechanical properties of micro-composite plates in the aircraft industry and other engineering applications.展开更多
Due to space constraints in urban areas,metro tunnels are typically constructed in pairs,with a small clearance.The interaction between twin tunnels leads to a significantly more complex ground deformation and stress ...Due to space constraints in urban areas,metro tunnels are typically constructed in pairs,with a small clearance.The interaction between twin tunnels leads to a significantly more complex ground deformation and stress distribution than that observed in a single tunnel scenario,particularly if the tunnels are excavated in sequence.A series of physical model tests were conducted to investigate soil deformation and stress disturbances caused by the excavation of twin tunnels.The test results indicate that the interaction between the twin tunnels was observed.Due to the soil arching effect,the excavation of Tunnel 2 increases the soil stress acting on Tunnel 1.An analytical method was proposed to determine soil stress considering the soil arching effect and the interaction between twin tunnels.The method categorized the relative locations between twin tunnels into non-influenced,partially influenced,and fully influenced scenarios.For non-influenced and fully influenced scenarios,the soil stresses above twin tunnels were calculated based on a symmetric major principal stress trace.For the partially influenced scenario,however,the soil arch above Tunnel 2 was asymmetric due to the interaction,and the stress distribution was obtained based on a new asymmetric major principal stress trace.The soil stress on Tunnel 1 was influenced by the load transferred from Tunnel 2 and calculated based on the force equilibrium.A comparison of the analytical and test results indicates that the proposed method effectively predicts the soil stress in the cover layer above twin tunnels excavated sequentially,considering the interaction and soil arching effects.展开更多
Pre-driven longwall retracement roadway(PLRR)is commonly used in large mine shaft.The support crushing disasters occur frequently during the retracement,and roof management is necessary.Taking the 31107 panel as resea...Pre-driven longwall retracement roadway(PLRR)is commonly used in large mine shaft.The support crushing disasters occur frequently during the retracement,and roof management is necessary.Taking the 31107 panel as research background,the roof breaking structure of PLRR is analyzed.It is concluded that the roof cutting with vertical hydraulic fracture(HF)at a specified position,that is,fixed-length roof cutting,can reduce support load and keep immediate roof intact.The extended finite element method(XFEM)is applied to simulate hydraulic fracturing.The results show that both the axial and transverse hydraulic fracturing cannot effectively create vertical HFs.Therefore,a novel construction method of vertical HF based on the stress shadow effect(SSE)is proposed.The stress reversal region and HF orientation caused by the prefabricated hydraulic fracture(PF)are verified in simulation.The sub-vertical HFs are obtained between two PFs,the vertical extension range of which is much larger than that of directional hydraulic fracturing.The new construction method was used to determine the field plan for fixed-length roof cutting.The roof formed a stable suspended structure and deformation of the main PLRR was improved after hydraulic fracturing.展开更多
The fatigue life evaluation of the girth butt weld within the welded cast steel joint was studied based on the extrapolation notch stress method.Firstly,the mesh sensitivity of the finite element model of the welded c...The fatigue life evaluation of the girth butt weld within the welded cast steel joint was studied based on the extrapolation notch stress method.Firstly,the mesh sensitivity of the finite element model of the welded cast steel joint was analyzed to determine the optimal mesh size.Based on the stress field analysis of the finite element model of the welded cast steel joint at the weld toe and weld root,the sharp model of the extrapolation notch stress method was applied to derive the effective notch stress of the rounded model belonging to the effective notch stress method,in which the key problem is to calculate the extrapolation point C,and the extrapolation point C has an exponential function relationship with the geometric parameters of the welded cast steel joint.By setting different values of geometric parameters,the corresponding value of parameter C is calculated,and then the functional relationship between the extrapolation point C and the geometric parameters can be obtained by the multiple linear regression analysis.Meanwhile,the fatigue life evaluation of the girth butt weld within welded cast steel joints based on the effective notch stress was performed according to the guideline recommended by the IIW(International Institute of Welding).The results indicate that the extrapolation notch stress method can effectively simplify the process of calculating the effective notch stress and accurately evaluate the fatigue life of the girth butt weld within welded cast steel joints.展开更多
Based on elastoplastic model, 2D and 3D finite element method （FEM） are used to calculate the stress and displacement distribution in the soft clay slope under gravity and uniform load at the slope top. Stability an...Based on elastoplastic model, 2D and 3D finite element method （FEM） are used to calculate the stress and displacement distribution in the soft clay slope under gravity and uniform load at the slope top. Stability analyses indicate that 3D boundary effect varies with the stress level of the slope. When the slope is stable, end effect of 3D space is not remarkable. When the stability decreases, end effect occurs; when the slope is at limit state, end effect reaches maximum. The energy causing slope failure spreads preferentially along y-z section, and when the failure resistance capability reaches the limit state, the energy can extend along x-axis direction. The 3D effect of the slope under uniform load on the top is related to the ratio of load influence width to slope height, and the effect is remarkable with the decrease of the ratio.展开更多
In wall-bounded turbulent flow calculations, the past focus has been directed to the modelling of the Reynolds-stress gradients. Not much attention has been paid to the effects of the numerical methods used to calcula...In wall-bounded turbulent flow calculations, the past focus has been directed to the modelling of the Reynolds-stress gradients. Not much attention has been paid to the effects of the numerical methods used to calculate these terms and the modelled equations. Discrepancies between model calculations and measurements are quite often attributed to incorrect modelling, while the suitability and accuracy of the numerical methods used are seldom scrutinized. Instead, alternate near-wall and Reynolds-stress models are proposed to remedy the incorrect turbulent flow calculations. On the other hand, if care is not taken in the numerical treatment of the Reynolds-stress gradient terms, physically unrealistic results and solution instability could occur. Previous studies by the author and his collaborators on the effects of numerical methods have shown that some of the more commonly used numerical methods could enhance numerical stability in the solution procedure but would introduce considerable inaccuracy to the results. The flow cases chosen to demonstrate these inaccuracies are a backstep flow and flow in a square duct, where flow complexities are present. The current investigation attempts to show that the above-mentioned effects of numerical methods could also occur in the calculation of a developing plane channel flow, where flow complexities are absent. In addition, this study shows that the results thus obtained lead to a predicted skin friction coefficient that is influenced more by the numerical method used than by the turbulence model invoked. Together, these results show that numerical treatment of the Reynolds-stress gradients in the equations play an important role, even for a developing plane channel flow.展开更多
The seepage property of low-permeability rock is of significant importance for the design and safety analysis of underground cavities. By using a self-developed test system, both permeability and porosity of granite f...The seepage property of low-permeability rock is of significant importance for the design and safety analysis of underground cavities. By using a self-developed test system, both permeability and porosity of granite from an underground oil storage depot were measured. In order to study the influence of rock types on permeability, a tight sandstone was selected as a contrast. The experimental results suggested that the porosity of this granite is less than 5% and permeability is low to 10–20 m^2 within the range of effective stress. During the loading process, both exponential relationship and power law can be utilized to describe the relationship between effective stress and permeability. However, power law matches the experimental data better during the unloading condition. The stress dependent porosity of granite during loading process can be described via an exponential relationship while the match between the model and experimental data can be improved by a power law in unloading paths. The correlation of permeability and porosity can be described in a power law form. Besides, granite shows great different evolution rules in permeability and porosity from sandstone. It is inferred that this difference can be attributed to the preparing of samples and different movements of microstructures subjected to effective stress.展开更多
Finite element method(FEM) simulations were employed to investigate the quenching residual stress distributions of 7085 aluminum alloy plates.The effect of dimensional variation on the quenching residual stress distri...Finite element method(FEM) simulations were employed to investigate the quenching residual stress distributions of 7085 aluminum alloy plates.The effect of dimensional variation on the quenching residual stress distributions was studied and discussed by using models with different dimensions(length,width,and thickness).The accuracy and efficiency of the models were verified by other numerical examples.The order of the dimension effects on the quenching residual stress distributions is:thickness> width=length.The maximum tensile stress and compressive stress increase from 33 to 190 and 39 to 270 MPa,respectively,as the thickness increases from 30 to 150 mm.The ultimate maximum tensile stress(about190 MPa) is equivalent to half of the quenching yield strength at 20℃,while the ultimate maximum compressive stress(about 300 MPa) is equivalent to 80 % of the quenching yield strength at 20℃.There are stress fluctuations at the edge of the large plate both in rolling and in transverse directions.The ratio of the fluctuation region along the rolling direction and transverse direction increases as the thickness increases,while it decreases as the length or width increases.The actual length of the fluctuation region is almost a constant value for the plates with a thickness of 115 mm(about 500 mm in length and 300 mm in width).展开更多
In this study, a size-dependent composite laminated skew Mindlin plate model is proposed based on a new modified couple stress theory. This plate model can be viewed as a simplified couple stress theory in engineering...In this study, a size-dependent composite laminated skew Mindlin plate model is proposed based on a new modified couple stress theory. This plate model can be viewed as a simplified couple stress theory in engineering mechanics. Governing equations and related boundary conditions are derived based on the principle of minimum potential energy. The Rayleigh–Ritz method is employed to obtain the numerical solutions of the center deflections of simply supported plates with different ply orientations. Numerical results show that the normalized center deflections obtained by the proposed model are always smaller than those obtained by the classical one, i.e. the present model can capture the scale effects of microstructures. Moreover, a phenomenon reveals that the ply orientation would make a significant influence on the magnitude of scale effects of composite laminated plates at micro scale. Additionally, the present model of thick skew plate can be degenerated to the model of Kirchhoff plate based on the modified couple stress theory by adopting the assumptions in Bernoulli–Euler beam and material isotropy.展开更多
In the simulation of rupture processes of seismic sources by using either numerical method or rock mechanics experiments, improper setting of the specimen size will influence the stress field near the faults. In this ...In the simulation of rupture processes of seismic sources by using either numerical method or rock mechanics experiments, improper setting of the specimen size will influence the stress field near the faults. In this study, 2D finite element method (FEM) was used to calculate the stress field of rock specimens in different sizes with fixed-size elliptic holes. The calculated stress field was compared with analytic solution for elliptic-hole problem in an infinite medium. Numerical results showed that boundary effect of a rock specimen with an elliptic hole on stress field under uniaxial compression cannot be neglected. Critical aspect ratio of the specimen is about 3:2, and critical ratio of distance between the tip of the hole and the border of specimen (d) to the major axis of the elliptic hole (l) is about 7.3. Numerical analysis on rock specimen size can provide theoretical reference for rock specimen experiments, and it is also helpful for setting of model sizes in numerical simulations of fault movement.展开更多
Time-dependent behaviors due to various mismatch strains are very important to the reliability of micro-/nano-devices.This paper aims at presenting an analytical model to study the viscoelastic stress relaxation of th...Time-dependent behaviors due to various mismatch strains are very important to the reliability of micro-/nano-devices.This paper aims at presenting an analytical model to study the viscoelastic stress relaxation of the laminated microbeam caused by mismatch strain.Firstly,Zhang’s two-variable method is used to establish a mechanical model for predicting the quasi-static stress relaxation of the laminated microbeam.Secondly,the related analytical solutions are obtained by combining the differential method and the eigenvalue method in the temporal domain.Finally,the influence of the substrateto-film thickness/modulus ratio on the relaxation responses of the laminated microbeam subject to a step load of the mismatch strain is studied.The results show that the present predictions are consistent with the previous theoretical studies.Furthermore,the thickness dependence of stress relaxation time of the laminated microbeam is jointly determined by the intrinsic structural evolution factors and tension-bending coupling state;the stress relaxation time can be controlled by adjusting the substrate-to-film thickness/modulus ratio.展开更多
An integrated approach is presented using field data input from measured geological information into numerical simulation for understanding effects of induced stress on geometries of multiple fracture propagation.We e...An integrated approach is presented using field data input from measured geological information into numerical simulation for understanding effects of induced stress on geometries of multiple fracture propagation.We establish a benchmark study based on comparison of field result with numerical computations.The comparison then acts as reference measures for studying effects of changing in-situ stress,fracturing fluid viscosity and fracture spacing on propagation and geometries of multiple fractures between neighbouring wellbores in an undisclosed gas field.This leverages understanding of more complexities associated with inter-well multiple fracture growth that are idealized as straight from certain perspectives.Although some studies focus on stress interference from pre-existing fractures,actual fracture propagation geometries may be far from the idealized scenarios.Therefore,the stress shadow effects between growing multiple hydraulic fractures,if not taken into account,can lead to unrealistic estimation of hydraulic fracture trajectories.Consequently,more attention should be paid to the actual propagation of hydraulic fractures.Actual field geologic information is provided through well-logging and field mapping data.Very short fractures were propagated in these wells before the operation was terminated due to technical problems.The reservoir depth in the area is about 2170 m.At such depth,quasi-brittleness of shale should be accounted for by using relevant methods that capture rock ductility such as traction separation law.Abaqus commercial software is used to conduct the numerical computation using extended finite element method based on cohesive zone modeling.Application of this technique is further validated using Kristianovich-Geertsma-de Klerk analytical solution.This study is important in field implementation of infill drilling with well-informed mechanism of fracture interference in the inter-well region.展开更多
The infiltration,evaporation and variation of the groundwater table have significant effects on the suction stress of the soils and the supporting earth pressures of the foundation excavation.The distribution of the s...The infiltration,evaporation and variation of the groundwater table have significant effects on the suction stress of the soils and the supporting earth pressures of the foundation excavation.The distribution of the suction stresses above the ground water table is derived under different fluxes at the ground surface,according to the soil-water characteristic parameters and the effective degree of saturation.In consideration of the cohesive stress formed from the soil suction stress and the relevant anti sliding effect,the calculation model of supporting earth pressures for foundation excavation is established by the variational limit equilibrium method under the steady flow condition.The evolution of the supporting earth pressures is studied in detail for foundation excavation under different fluxes at the ground surface.The effects of the soil-water characteristic parameters,the ground water table and the internal friction angle on the supporting earth pressures are discussed.The results show that the suction stress is reduced because of the infiltration,and thus the supporting earth pressure increases.The larger the air-entry pressures and the pore size are,the smaller the supporting earth pressures are.The higher the ground water table is,the larger the supporting earth pressures are.In order to reduce the construction risk,the effects of the suction stress and the evolution of the potential critical sliding surface should be considered during the calculation of the supporting earth pressures.展开更多
In micro-electro-mechanical systems,interface expansion issues are commonly encountered,and due to their small size,they often exist at the micro-or nano-scale.The influence of the micro-structural effect on interface...In micro-electro-mechanical systems,interface expansion issues are commonly encountered,and due to their small size,they often exist at the micro-or nano-scale.The influence of the micro-structural effect on interface mechanics cannot be ignored.This paper focuses on studying the impact of micro-structural effect on interface crack propagation.Modified couple stress theory(MCST)is used to study the buckling delamination of ultra-thin film-substrate systems.The equivalent elastic modulus(EEM)and equivalent flexural rigidity(EFR)are derived based on MCST.Substituting EEM and EFR into the classical Kirchhoff plate theory,the governing equations of ultra-thin film-substrate system with micro-structural effect can be obtained.The finite element method(FEM)was used to calculate the critical strain energy release rate for crack extension.Differences between the three theoretical approaches of MCST,classical theory(CT),and FEM were compared.The effects of stress ratioσ/σ_(c),initial crack length,film thickness,and micro-structural effect parameters on crack extension were analyzed.The results show that the FEM calculations coincide with the CT calculations.The stress ratioσ/σ_(c),initial crack length,film thickness,and micro-structural effect parameters have significantly influence crack extension.展开更多
The current national criteria cannot accurately reflect the multi-crack interference effect of defective gas pipelines,and thus result in conservative assessment results.In order to improve the safety assessment accur...The current national criteria cannot accurately reflect the multi-crack interference effect of defective gas pipelines,and thus result in conservative assessment results.In order to improve the safety assessment accuracy of defective gas pipelines,we compared the singular element method with the three-dimensional virtual crack closure technology(3D-VCCT),then,selected the 3D-VCCT to build a multi-crack interference model of gas pipelines,and finally,studied the interference effect of subsidiary cracks on the front points of the main crack by analyzing the variation of the interference factors between cracks.The following research results were obtained:first,as the subsidiary crack size changes,the strongest interference effect of a parallel collinear crack occurs near the surface of the crack,exhibiting an enhances effect;second,under certain conditions,the interference effect of a parallel coaxial crack is more significant at the deepest point and the surface point of the crack,and exhibits an enhanced effect at the deepest point and a weakened effect at the surface point;third,when the horizontal spacing between cracks is 6 times greater than the major semi-axis of the main crack,the interference effect of parallel collinear cracks and parallel bias cracks is very little,and multi-cracks can be simplified as a single crack for analysis.The research results not only can provide technical references for fracture analysis of multi-crack gas pipelines,but also can provide theoretical basis for the safety assessment of multi-cracks of the gas pipelines.展开更多
基金Projects(03JJY3078, 04JJ40032) supported by the Natural Science Foundation of Hunan Province, China project(03A006) supported by Scientific Research Fund of Hunan Provincial Education Department, China
文摘Based on the analysis method for tailings dam in upstream raising method presently used in metallurgy and nonferrous metals tailings depository in the world, an effective stress analysis method of seismic response for high tailings dam was developed according to the results of engineering geological exploration, static and dynamic test and stability analysis on Baizhishan tailing dam 113.5 m high. The law of generation, diffusion and dissipation of seismic pore water pressure during and after earthquake was investigated, and the results of tailings dam’s acceleration, seismic dynamic stress and pore water pressure were obtained. The results show that the seismic stability and liquefaction resistance of high tailings dam are strengthened remarkably, and the scope and depth of liquefaction area at the top of dam are reduced greatly. The interior stress is compressive stress, the stress level of every element is less than 1.0 and the safety coefficient of every element is greater than 1.0. The safety coefficient against liquefaction of every element of tailing dam is greater than 1.5 according to the effective stress analysis of seismic response by finite element method. The calculated results prove that liquefaction is the main reason of seismic failure of high tailing dams, and the effect of seismic inertia forces on high tailing dams’ stability during earthquake is secondary reason.
基金supported by the Iranian Nanotechnology Development Committee and the University of Kashan(No.363452/10)
文摘The size-dependent effect on the biaxial and shear nonlinear buckling analysis of an isotropic and orthotropic micro-plate based on the surface stress, the modified couple stress theory （MCST）, and the nonlocal elasticity theories using the differential quadrature method （DQM） is presented. Main advantages of the MCST over the classical theory （CT） are the inclusion of the asymmetric couple stress tensor and the consideration of only one material length scale parameter. Based on the nonlinear von Karman assumption, the governing equations of equilibrium for the micro-classical plate consid- ering midplane displacements are derived based on the minimum principle of potential energy. Using the DQM, the biaxial and shear critical buckling loads of the micro-plate for various boundary conditions are obtained. Accuracy of the obtained results is validated by comparing the solutions with those reported in the literature. A parametric study is conducted to show the effects of the aspect ratio, the side-to-thickness ratio, Eringen＇s nonlocal parameter, the material length scale parameter, Young＇s modulus of the surface layer, the surface residual stress, the polymer matrix coefficients, and various boundary conditions on the dimensionless uniaxial, biaxial, and shear critical buckling loads. The results indicate that the critical buckling loads are strongly sensitive to Eringen＇s nonlocal parameter, the material length scale parameter, and the surface residual stress effects, while the effect of Young＇s modulus of the surface layer on the critical buckling load is negligible. Also, considering the size dependent effect causes the increase in the stiffness of the orthotropic micro-plate. The results show that the critical biaxial buckling load increases with an increase in G12/E2 and vice versa for E1/E2. It is shown that the nonlinear biaxial buckling ratio decreases as the aspect ratio increases and vice versa for the buckling amplitude. Because of the most lightweight micro-composite materials with high strength/weight and stiffness/weight ratios, it is anticipated that the results of the present work are useful in experimental characterization of the mechanical properties of micro-composite plates in the aircraft industry and other engineering applications.
基金supported by the National Natural Science Foundation of China(Grant No.52308463)the Shanghai Rising-Star Program(Grant No.23YF1449100)the Fundamental Research Funds for the Central Universities(Grant No.2023-2-ZD08).
文摘Due to space constraints in urban areas,metro tunnels are typically constructed in pairs,with a small clearance.The interaction between twin tunnels leads to a significantly more complex ground deformation and stress distribution than that observed in a single tunnel scenario,particularly if the tunnels are excavated in sequence.A series of physical model tests were conducted to investigate soil deformation and stress disturbances caused by the excavation of twin tunnels.The test results indicate that the interaction between the twin tunnels was observed.Due to the soil arching effect,the excavation of Tunnel 2 increases the soil stress acting on Tunnel 1.An analytical method was proposed to determine soil stress considering the soil arching effect and the interaction between twin tunnels.The method categorized the relative locations between twin tunnels into non-influenced,partially influenced,and fully influenced scenarios.For non-influenced and fully influenced scenarios,the soil stresses above twin tunnels were calculated based on a symmetric major principal stress trace.For the partially influenced scenario,however,the soil arch above Tunnel 2 was asymmetric due to the interaction,and the stress distribution was obtained based on a new asymmetric major principal stress trace.The soil stress on Tunnel 1 was influenced by the load transferred from Tunnel 2 and calculated based on the force equilibrium.A comparison of the analytical and test results indicates that the proposed method effectively predicts the soil stress in the cover layer above twin tunnels excavated sequentially,considering the interaction and soil arching effects.
基金financially supported by the Postgraduate Research&Practice Innovation Program of Jiangsu Province(KYCX21_2358)the National Key Research and Development Program of China(2020YFB1314204)National Natural Science Foundation of China(No.52074239)。
文摘Pre-driven longwall retracement roadway(PLRR)is commonly used in large mine shaft.The support crushing disasters occur frequently during the retracement,and roof management is necessary.Taking the 31107 panel as research background,the roof breaking structure of PLRR is analyzed.It is concluded that the roof cutting with vertical hydraulic fracture(HF)at a specified position,that is,fixed-length roof cutting,can reduce support load and keep immediate roof intact.The extended finite element method(XFEM)is applied to simulate hydraulic fracturing.The results show that both the axial and transverse hydraulic fracturing cannot effectively create vertical HFs.Therefore,a novel construction method of vertical HF based on the stress shadow effect(SSE)is proposed.The stress reversal region and HF orientation caused by the prefabricated hydraulic fracture(PF)are verified in simulation.The sub-vertical HFs are obtained between two PFs,the vertical extension range of which is much larger than that of directional hydraulic fracturing.The new construction method was used to determine the field plan for fixed-length roof cutting.The roof formed a stable suspended structure and deformation of the main PLRR was improved after hydraulic fracturing.
基金The National Key Research and Development Program of China(No.2017YFC0805100),the National Natural Science Foundation of China(No.51578137)the Priority Academic Program Development of Jiangsu Higher Education Institutions,the Open Research Fund Program of Jiangsu Key Laboratory of Engineering Mechanics.
文摘The fatigue life evaluation of the girth butt weld within the welded cast steel joint was studied based on the extrapolation notch stress method.Firstly,the mesh sensitivity of the finite element model of the welded cast steel joint was analyzed to determine the optimal mesh size.Based on the stress field analysis of the finite element model of the welded cast steel joint at the weld toe and weld root,the sharp model of the extrapolation notch stress method was applied to derive the effective notch stress of the rounded model belonging to the effective notch stress method,in which the key problem is to calculate the extrapolation point C,and the extrapolation point C has an exponential function relationship with the geometric parameters of the welded cast steel joint.By setting different values of geometric parameters,the corresponding value of parameter C is calculated,and then the functional relationship between the extrapolation point C and the geometric parameters can be obtained by the multiple linear regression analysis.Meanwhile,the fatigue life evaluation of the girth butt weld within welded cast steel joints based on the effective notch stress was performed according to the guideline recommended by the IIW(International Institute of Welding).The results indicate that the extrapolation notch stress method can effectively simplify the process of calculating the effective notch stress and accurately evaluate the fatigue life of the girth butt weld within welded cast steel joints.
文摘Based on elastoplastic model, 2D and 3D finite element method （FEM） are used to calculate the stress and displacement distribution in the soft clay slope under gravity and uniform load at the slope top. Stability analyses indicate that 3D boundary effect varies with the stress level of the slope. When the slope is stable, end effect of 3D space is not remarkable. When the stability decreases, end effect occurs; when the slope is at limit state, end effect reaches maximum. The energy causing slope failure spreads preferentially along y-z section, and when the failure resistance capability reaches the limit state, the energy can extend along x-axis direction. The 3D effect of the slope under uniform load on the top is related to the ratio of load influence width to slope height, and the effect is remarkable with the decrease of the ratio.
文摘In wall-bounded turbulent flow calculations, the past focus has been directed to the modelling of the Reynolds-stress gradients. Not much attention has been paid to the effects of the numerical methods used to calculate these terms and the modelled equations. Discrepancies between model calculations and measurements are quite often attributed to incorrect modelling, while the suitability and accuracy of the numerical methods used are seldom scrutinized. Instead, alternate near-wall and Reynolds-stress models are proposed to remedy the incorrect turbulent flow calculations. On the other hand, if care is not taken in the numerical treatment of the Reynolds-stress gradient terms, physically unrealistic results and solution instability could occur. Previous studies by the author and his collaborators on the effects of numerical methods have shown that some of the more commonly used numerical methods could enhance numerical stability in the solution procedure but would introduce considerable inaccuracy to the results. The flow cases chosen to demonstrate these inaccuracies are a backstep flow and flow in a square duct, where flow complexities are present. The current investigation attempts to show that the above-mentioned effects of numerical methods could also occur in the calculation of a developing plane channel flow, where flow complexities are absent. In addition, this study shows that the results thus obtained lead to a predicted skin friction coefficient that is influenced more by the numerical method used than by the turbulence model invoked. Together, these results show that numerical treatment of the Reynolds-stress gradients in the equations play an important role, even for a developing plane channel flow.
基金Projects(11172090,51479049,11272113,11572110,51209075)supported by the National Natural Science Foundation of ChinaProject(BK2012809)supported by the Natural Science Foundation of Jiangsu Province,ChinaProject(201406710042)supported by China Scholarship Council
文摘The seepage property of low-permeability rock is of significant importance for the design and safety analysis of underground cavities. By using a self-developed test system, both permeability and porosity of granite from an underground oil storage depot were measured. In order to study the influence of rock types on permeability, a tight sandstone was selected as a contrast. The experimental results suggested that the porosity of this granite is less than 5% and permeability is low to 10–20 m^2 within the range of effective stress. During the loading process, both exponential relationship and power law can be utilized to describe the relationship between effective stress and permeability. However, power law matches the experimental data better during the unloading condition. The stress dependent porosity of granite during loading process can be described via an exponential relationship while the match between the model and experimental data can be improved by a power law in unloading paths. The correlation of permeability and porosity can be described in a power law form. Besides, granite shows great different evolution rules in permeability and porosity from sandstone. It is inferred that this difference can be attributed to the preparing of samples and different movements of microstructures subjected to effective stress.
基金financially supported by the National Basic Research Program of China(No.2012CB619504)the National Natural Science Foundation of China(No.51274046)
文摘Finite element method(FEM) simulations were employed to investigate the quenching residual stress distributions of 7085 aluminum alloy plates.The effect of dimensional variation on the quenching residual stress distributions was studied and discussed by using models with different dimensions(length,width,and thickness).The accuracy and efficiency of the models were verified by other numerical examples.The order of the dimension effects on the quenching residual stress distributions is:thickness> width=length.The maximum tensile stress and compressive stress increase from 33 to 190 and 39 to 270 MPa,respectively,as the thickness increases from 30 to 150 mm.The ultimate maximum tensile stress(about190 MPa) is equivalent to half of the quenching yield strength at 20℃,while the ultimate maximum compressive stress(about 300 MPa) is equivalent to 80 % of the quenching yield strength at 20℃.There are stress fluctuations at the edge of the large plate both in rolling and in transverse directions.The ratio of the fluctuation region along the rolling direction and transverse direction increases as the thickness increases,while it decreases as the length or width increases.The actual length of the fluctuation region is almost a constant value for the plates with a thickness of 115 mm(about 500 mm in length and 300 mm in width).
基金supported by the National Natural Sciences Foundation of China(No.11572204)
文摘In this study, a size-dependent composite laminated skew Mindlin plate model is proposed based on a new modified couple stress theory. This plate model can be viewed as a simplified couple stress theory in engineering mechanics. Governing equations and related boundary conditions are derived based on the principle of minimum potential energy. The Rayleigh–Ritz method is employed to obtain the numerical solutions of the center deflections of simply supported plates with different ply orientations. Numerical results show that the normalized center deflections obtained by the proposed model are always smaller than those obtained by the classical one, i.e. the present model can capture the scale effects of microstructures. Moreover, a phenomenon reveals that the ply orientation would make a significant influence on the magnitude of scale effects of composite laminated plates at micro scale. Additionally, the present model of thick skew plate can be degenerated to the model of Kirchhoff plate based on the modified couple stress theory by adopting the assumptions in Bernoulli–Euler beam and material isotropy.
基金National Natural Science Foundation of China (40234042).
文摘In the simulation of rupture processes of seismic sources by using either numerical method or rock mechanics experiments, improper setting of the specimen size will influence the stress field near the faults. In this study, 2D finite element method (FEM) was used to calculate the stress field of rock specimens in different sizes with fixed-size elliptic holes. The calculated stress field was compared with analytic solution for elliptic-hole problem in an infinite medium. Numerical results showed that boundary effect of a rock specimen with an elliptic hole on stress field under uniaxial compression cannot be neglected. Critical aspect ratio of the specimen is about 3:2, and critical ratio of distance between the tip of the hole and the border of specimen (d) to the major axis of the elliptic hole (l) is about 7.3. Numerical analysis on rock specimen size can provide theoretical reference for rock specimen experiments, and it is also helpful for setting of model sizes in numerical simulations of fault movement.
基金Project supported by the National Natural Science Foundation of China(Nos.12172204,11772182,11272193,and 10872121)the Program of Shanghai Municipal Education Commission(No.2019-01-07-00-09-E00018)。
文摘Time-dependent behaviors due to various mismatch strains are very important to the reliability of micro-/nano-devices.This paper aims at presenting an analytical model to study the viscoelastic stress relaxation of the laminated microbeam caused by mismatch strain.Firstly,Zhang’s two-variable method is used to establish a mechanical model for predicting the quasi-static stress relaxation of the laminated microbeam.Secondly,the related analytical solutions are obtained by combining the differential method and the eigenvalue method in the temporal domain.Finally,the influence of the substrateto-film thickness/modulus ratio on the relaxation responses of the laminated microbeam subject to a step load of the mismatch strain is studied.The results show that the present predictions are consistent with the previous theoretical studies.Furthermore,the thickness dependence of stress relaxation time of the laminated microbeam is jointly determined by the intrinsic structural evolution factors and tension-bending coupling state;the stress relaxation time can be controlled by adjusting the substrate-to-film thickness/modulus ratio.
基金the support of the project on“research on the basic theory of ultra-deep gas reservoir development,key Project of National Natural Science Foundation of China.”。
文摘An integrated approach is presented using field data input from measured geological information into numerical simulation for understanding effects of induced stress on geometries of multiple fracture propagation.We establish a benchmark study based on comparison of field result with numerical computations.The comparison then acts as reference measures for studying effects of changing in-situ stress,fracturing fluid viscosity and fracture spacing on propagation and geometries of multiple fractures between neighbouring wellbores in an undisclosed gas field.This leverages understanding of more complexities associated with inter-well multiple fracture growth that are idealized as straight from certain perspectives.Although some studies focus on stress interference from pre-existing fractures,actual fracture propagation geometries may be far from the idealized scenarios.Therefore,the stress shadow effects between growing multiple hydraulic fractures,if not taken into account,can lead to unrealistic estimation of hydraulic fracture trajectories.Consequently,more attention should be paid to the actual propagation of hydraulic fractures.Actual field geologic information is provided through well-logging and field mapping data.Very short fractures were propagated in these wells before the operation was terminated due to technical problems.The reservoir depth in the area is about 2170 m.At such depth,quasi-brittleness of shale should be accounted for by using relevant methods that capture rock ductility such as traction separation law.Abaqus commercial software is used to conduct the numerical computation using extended finite element method based on cohesive zone modeling.Application of this technique is further validated using Kristianovich-Geertsma-de Klerk analytical solution.This study is important in field implementation of infill drilling with well-informed mechanism of fracture interference in the inter-well region.
基金the National Natural Science Foundation of China(No.41272288)。
文摘The infiltration,evaporation and variation of the groundwater table have significant effects on the suction stress of the soils and the supporting earth pressures of the foundation excavation.The distribution of the suction stresses above the ground water table is derived under different fluxes at the ground surface,according to the soil-water characteristic parameters and the effective degree of saturation.In consideration of the cohesive stress formed from the soil suction stress and the relevant anti sliding effect,the calculation model of supporting earth pressures for foundation excavation is established by the variational limit equilibrium method under the steady flow condition.The evolution of the supporting earth pressures is studied in detail for foundation excavation under different fluxes at the ground surface.The effects of the soil-water characteristic parameters,the ground water table and the internal friction angle on the supporting earth pressures are discussed.The results show that the suction stress is reduced because of the infiltration,and thus the supporting earth pressure increases.The larger the air-entry pressures and the pore size are,the smaller the supporting earth pressures are.The higher the ground water table is,the larger the supporting earth pressures are.In order to reduce the construction risk,the effects of the suction stress and the evolution of the potential critical sliding surface should be considered during the calculation of the supporting earth pressures.
基金supported by the National Key R&D Program of China(No.2022YFB3207100)the Hubei Provincial Strategic Scientist Training Plan(No.2022EJD009)the Fundamental Research Funds for the Central Universities(No.2042023kf1041).
文摘In micro-electro-mechanical systems,interface expansion issues are commonly encountered,and due to their small size,they often exist at the micro-or nano-scale.The influence of the micro-structural effect on interface mechanics cannot be ignored.This paper focuses on studying the impact of micro-structural effect on interface crack propagation.Modified couple stress theory(MCST)is used to study the buckling delamination of ultra-thin film-substrate systems.The equivalent elastic modulus(EEM)and equivalent flexural rigidity(EFR)are derived based on MCST.Substituting EEM and EFR into the classical Kirchhoff plate theory,the governing equations of ultra-thin film-substrate system with micro-structural effect can be obtained.The finite element method(FEM)was used to calculate the critical strain energy release rate for crack extension.Differences between the three theoretical approaches of MCST,classical theory(CT),and FEM were compared.The effects of stress ratioσ/σ_(c),initial crack length,film thickness,and micro-structural effect parameters on crack extension were analyzed.The results show that the FEM calculations coincide with the CT calculations.The stress ratioσ/σ_(c),initial crack length,film thickness,and micro-structural effect parameters have significantly influence crack extension.
基金Project supported by National Science and Technology Infrastructure Program“Risk-based key technologies for preventing pipeline accidents”(No.:2011BAK06B01-11).
文摘The current national criteria cannot accurately reflect the multi-crack interference effect of defective gas pipelines,and thus result in conservative assessment results.In order to improve the safety assessment accuracy of defective gas pipelines,we compared the singular element method with the three-dimensional virtual crack closure technology(3D-VCCT),then,selected the 3D-VCCT to build a multi-crack interference model of gas pipelines,and finally,studied the interference effect of subsidiary cracks on the front points of the main crack by analyzing the variation of the interference factors between cracks.The following research results were obtained:first,as the subsidiary crack size changes,the strongest interference effect of a parallel collinear crack occurs near the surface of the crack,exhibiting an enhances effect;second,under certain conditions,the interference effect of a parallel coaxial crack is more significant at the deepest point and the surface point of the crack,and exhibits an enhanced effect at the deepest point and a weakened effect at the surface point;third,when the horizontal spacing between cracks is 6 times greater than the major semi-axis of the main crack,the interference effect of parallel collinear cracks and parallel bias cracks is very little,and multi-cracks can be simplified as a single crack for analysis.The research results not only can provide technical references for fracture analysis of multi-crack gas pipelines,but also can provide theoretical basis for the safety assessment of multi-cracks of the gas pipelines.
