During the development of hot dry rock,the research on thermal fatigue damage caused by thermal shock of cold and heat cycles is the basis that ensures the long-term utilization of geothermal resources,but there are n...During the development of hot dry rock,the research on thermal fatigue damage caused by thermal shock of cold and heat cycles is the basis that ensures the long-term utilization of geothermal resources,but there are not enough relevant studies at present.Based on this,the thermal damage tests of granite at different temperatures(250,350,450°C)and quenching cycles(1,5,10,15 cycles)were carried out.Preliminary reveals the damage mechanism and heat transfer law of the quenching cycle effect on hot dry rock.The results show that with the increase of temperature and cycles,the uneven thermal expansion of minerals and the thermal shock caused by quenching promote the crack development of granite,resulting in the decrease of P-wave velocity,thermal conductivity and uniaxial compressive strength of granite.Meanwhile,the COMSOL was used to simulate the heat transfer of hot dry rock under different heat treatment conditions.It concluded that the increase in the number of quenching cycles reduced the heat transfer capacity of the granite,especially more than 10 quenching cycles,which also reflects that the thermal fatigue damage leads to a longer time for the temperature recovery of the hot dry rock mass.In addition,the three-dimensional nonlinear fitting relationship among thermal conductivity,temperature and cycle number was established for the first time,which can better reveal the change rule of thermal conductivity after quenching thermal fatigue effect of hot dry rock.The research results provide theoretical support for hot dry rock reservoir reconstruction and production efficiency evaluation.展开更多
The former studies indicate that loading rates significantly affect dynamic behavior of brittle materials,for instance,the dynamic compressive and tensile strength increase with loading rates.However,there still are m...The former studies indicate that loading rates significantly affect dynamic behavior of brittle materials,for instance,the dynamic compressive and tensile strength increase with loading rates.However,there still are many unknown or partially unknown aspects.For example,whether loading rates have effect on crack dynamic propagating behavior(propagation toughness,velocity and arrest,etc).To further explore the effect of loading rates on crack dynamic responses,a large-size single-cleavage trapezoidal open(SCTO)specimen was proposed,and impacting tests using the SCTO specimen under drop plate impact were conducted.Crack propagation gauges(CPGs)were employed in measuring impact loads,crack propagation time and velocities.In order to verify the testing result,the corresponding numerical model was established using explicit dynamic software AUTODYN,and the simulation result is basically consistent with the experimental results.The ABAQUS software was used to calculate the dynamic SIFs.The universal function was calculated by fractal method.The experimental-numerical method was employed in determining initiation toughness and propagation toughness.The results indicate that crack propagating velocities,dynamic fracture toughness and energy release rates increase with loading rates;crack delayed initiation time decreases with loading rates.展开更多
The prevention and the reduction of the rockfall are the common measures of the prevention and the reduction of disasters.When the rock-shed resists the impact of the rockfall,the force that acts on the structure cons...The prevention and the reduction of the rockfall are the common measures of the prevention and the reduction of disasters.When the rock-shed resists the impact of the rockfall,the force that acts on the structure consists of the cushion dead load and the impact-induced load,of which the dynamic process of the propagation of the impactinduced load is complex.Therefore,we conducted a numerical study to investigate the impact of the rockfall.Considering the highly discrete characteristic of the sand,we developed a numerical model on the basis of the discrete element method(DEM).The numerical model,which simulation results were validated by the results of real-scale experiments,was used to investigate the dynamic response of the impact force of the rockfall and the transmission of the impact force under the different magnitude of the falling height and the different thickness of the sand cushion.The results of our study indicated that the cushion thickness had little effect on the impact of the rockfall,and the dense sand cushion generated higher impact force than did the loose sand cushion.Although the high thickness enhanced the buffer performance of the sand cushion,the additional force induced by the dead load of sand cushion was significant.Therefore,to determine the appropriate thickness of the sand cushion,we suggested designers consider the buffer performance and the dead load of the sand cushion.The analysis presented in this paper provided a practical estimation of the impact-induced force of the thick sand cushion.展开更多
Natural gas hydrate(NGH)holds great promise as a source of clean energy.It is critical for acquiring the largest possible in situ NGH core for NGH eigen features and resource assessment.However,the existing NGH coring...Natural gas hydrate(NGH)holds great promise as a source of clean energy.It is critical for acquiring the largest possible in situ NGH core for NGH eigen features and resource assessment.However,the existing NGH coring technology has limitations,such as temperature increments,limited coring diameters,low coring rates,and complex coring structures.Therefore,this study designs and proposes an NGH freezing coring(NGHFC)method and verifies the freezing and coring capacities of the NGHFC method in laboratories and experimental wells.Results suggest that NGHFC shows good freezing and heat-retention properties.A freezing core heat transfer model is developed.According to the actual air temperature and operating time,the optimum initial temperature of the cold source can be determined using this model.The average coring rate of NGHFC can reach 77.86%.The research results will provide a new idea of coring gas hydrates.展开更多
In-situ pressure-preserved coring(IPP-Coring)is considered to be the most reliable and efficient method for the identification of the scale of oil and gas resources.During IPP-Coring,because the rotation behavior of t...In-situ pressure-preserved coring(IPP-Coring)is considered to be the most reliable and efficient method for the identification of the scale of oil and gas resources.During IPP-Coring,because the rotation behavior of the pressure controller valve cover in different medium environments is unclear,interference between the valve cover and inner pipe may occur and negatively affect the IPP-Coring success rate.To address this issue,we conducted a series of indoor experiments employing a high-speed camera to gain greater insights into the valve cover rotation behavior in different medium environments,e.g.,air,water,and simulated drilling fluids.The results indicated that the variation in the valve cover rotation angle in the air and fluid environments can be described by a one-phase exponential decay function with a constant time parameter and by biphasic dose response function,respectively.The rotation behavior in the fluid environments exhibited distinct elastic and gravitational acceleration zones.In the fluid environments,the density clearly impacted the valve cover closing time and rotation behavior,whereas the effect of viscosity was very slight.This can be attributed to the negligible influence of the fluid viscosity on the drag coefficient found in this study;meanwhile,the density can increase the buoyancy and the time period during which the valve cover experienced a high drag coefficient.Considering these results,control schemes for the valve cover rotation behavior during IPP-Coring were proposed for different layers and geological conditions in which the different drilling fluids should be used,e.g.,the use of a high-density valve cover in high-pore pressure layers.展开更多
Cross-laminated bamboo(CLB)have a high strength to weight ratio and stable bidirectional mechanical properties.Inspired by the investigation on cross-laminated timber(CLT)rocking walls,CLB rocking walls with conventio...Cross-laminated bamboo(CLB)have a high strength to weight ratio and stable bidirectional mechanical properties.Inspired by the investigation on cross-laminated timber(CLT)rocking walls,CLB rocking walls with conventional friction dampers(CFDs)are studied in this paper.To investigate the mechanical properties of the CLB rocking wall,seven tests are conducted under a cyclic loading scheme,and different test parameters,including the existence of the CFDs,the moment ratio,and the loading times,are discussed.The test results show a bilinear behavior of the CLB rocking wall.The small residual displacements of the CLB rocking wall demonstrate an idealized self-centering capacity.The cumulative energy dissipation curves indicate that the energy dissipation capacity of the CLB rocking wall can be greatly improved with CFDs.The limit states of the CLB rocking wall under a lateral force are proposed based on the strains,stress,and damage level of the CLB material and posttensioned rebar.In addition,an analytical model of the CLB rocking wall is developed based on the proposed limit states of the CLB rocking wall to evaluate the hysteretic response of the CLB rocking wall,and the model is validated by the experimental data.The comparison results show the potential value of the analytical model for engineering design.展开更多
In-situ upgrading by heating is feasible for low-maturity shale oil,where the pore space dynamically evolves.We characterize this response for a heated substrate concurrently imaged by SEM.We systematically follow the...In-situ upgrading by heating is feasible for low-maturity shale oil,where the pore space dynamically evolves.We characterize this response for a heated substrate concurrently imaged by SEM.We systematically follow the evolution of pore quantity,size(length,width and cross-sectional area),orientation,shape(aspect ratio,roundness and solidity)and their anisotropy—interpreted by machine learning.Results indicate that heating generates new pores in both organic matter and inorganic minerals.However,the newly formed pores are smaller than the original pores and thus reduce average lengths and widths of the bedding-parallel pore system.Conversely,the average pore lengths and widths are increased in the bedding-perpendicular direction.Besides,heating increases the cross-sectional area of pores in low-maturity oil shales,where this growth tendency fluctuates at<300℃ but becomes steady at>300℃.In addition,the orientation and shape of the newly-formed heating-induced pores follow the habit of the original pores and follow the initial probability distributions of pore orientation and shape.Herein,limited anisotropy is detected in pore direction and shape,indicating similar modes of evolution both bedding-parallel and bedding-normal.We propose a straightforward but robust model to describe evolution of pore system in low-maturity oil shales during heating.展开更多
The sealing performance of contact interfaces plays the most important role in the design and operation of the in-situ pressure-preserved coring system.To meet the demand of ultra-high pressure-retained coring for oil...The sealing performance of contact interfaces plays the most important role in the design and operation of the in-situ pressure-preserved coring system.To meet the demand of ultra-high pressure-retained coring for oil and gas exploration in deep reservoirs,a quantitative analysis of the contact mechanical behavior of the pressure controller was performed.Based on the micro-contact theory of rough surfaces,a three-dimensional numerical model of the rough contact interface between the valve cover and the valve seat was constructed,and the micro-contact behavior of the metal contact surfaces was comprehensively studied.The results show that the actual contact area of the valve interface increases with the increase of surface roughness before the critical contact point,but decreases after that.Compared with the real contact model with double rough surfaces,although the simplified hard-contact model with a single rough surface can reflect the micro-contact behavior of the rough surface to a certain extent,it cannot truly reveal the microchannel morphology between the sealing interfaces under pressure.Therefore,the realistic double-rough-surface model should be recommended to evaluate the sealing performance of coring tools,particularly for high pressure conditions.The material properties of valves have a significant effect on the contact characteristics of rough surfaces,which suggested that the actual contact area decreases with the increase of the elastic modulus of the contact material under the same loading conditions.The knowledge of this work could help to enhance the seal design of pressure controllers for in-situ pressure-preserved coring.展开更多
Texture is inevitably introduced during the manufacturing of most NiTi shape memory alloys(SMAs),and the textured nanocrystalline NiTi has been extensively employed in engineering.However,the effect of texture,and the...Texture is inevitably introduced during the manufacturing of most NiTi shape memory alloys(SMAs),and the textured nanocrystalline NiTi has been extensively employed in engineering.However,the effect of texture,and the joint effect of grain size(GS)and texture on the functional properties of NiTi SMAs and the corresponding microscopic mechanisms have not been clearly understood yet.In this work,based on the phase field method,the effect of texture on the GS-dependent functional properties of NiTi SMAs,including super-elasticity(SE),one-way shape memory effect(OWSME),and stress-assisted two-way shape memory effect(SATWSME),is investigated,and the corresponding microscopic mechanisms are revealed.Moreover,the samples with discrete geometrical gradients and/or texture gradients are designed to achieve graded functional properties.The simulation results indicate that the dependence of functional properties on texture is due to the effect of crystallographic orientation on martensite transformation and reorientation,which can lead to different inelastic strains.In the designed samples with texture gradients,the stress–strain responses of sheets with various textures are different,allowing for the coordination of overall deformation of the sample by combining such sheets,with varying inelastic deformation degrees.Thus,the overall response of the sample differs from that without texture gradient,leading to the achievement of graded functional properties.The simulation results and new findings in this work contribute to a deeper understanding of the effects of texture,GS,and their interaction on the functional properties of SMAs,and provide valuable reference for the design and development of SMA-based devices with desired functional properties.展开更多
The inelastic deformations of shape memory alloys(SMAs)always show poor controllability due to the avalanche-like martensite transformation,and the effective control for the deformation of precision de-vices has been ...The inelastic deformations of shape memory alloys(SMAs)always show poor controllability due to the avalanche-like martensite transformation,and the effective control for the deformation of precision de-vices has been not yet mature.In this work,the phase field method was used to investigate the shape memory effects(SMEs)of NiTi SMAs undergoing grain size(GS)engineering,to obtain tunable one-way and stress-assisted two-way SMEs(OWSME and SATWSME).The OWSME and SATWSME of the systems with various gradient-nanograin structures and bimodal grain structure,as well as that with geometric gradients were simulated.The simulated results indicate that due to the GS dependences of martensite transformation and reorientation,the occurrence and expansion of martensite reorientation,martensite transformation and its reverse can be efficaciously controlled via the GS engineering.When combining the GS engineering and geometric gradient design,since the effects of GS and stress gradient can be su-perimposed or competing,and the responses of martensite reorientation,martensite transformation and its reverse to this are different,the OWSME and SATWSME of the geometrically graded systems with various nanograin structures can exhibit different improvements in controllability.In short,the reorienta-tion hardening modulus during OWSME is increased and the transformation temperature window during SATWSME is widened by GS engineering,indicating the improved controllability of SMEs.The optimal GS engineering schemes revealed in this work provide the basic reference and guidance for designing tun-able SMEs and producing NiTi-based driving devices catering to desired functional performance in various engineering fields.展开更多
In-situ thermal upgrading is used to tune the pore system in low-maturity oil shales. We introduce fractal dimension(D), form factor(ff) and stochastic entropy(H) to quantify the heating-induced evolution of pore morp...In-situ thermal upgrading is used to tune the pore system in low-maturity oil shales. We introduce fractal dimension(D), form factor(ff) and stochastic entropy(H) to quantify the heating-induced evolution of pore morphological complexity and azimuthal disorder and develop a model to estimate the impact on seepage capacity via permeability. Experiments are conducted under recreated in-situ temperatures and consider anisotropic properties—both parallel and perpendicular to bedding. Results indicate that azimuthal distribution of pores in the bedding-parallel direction are dispersed, while those in the bedding-perpendicular direction are concentrated. D values indicate that higher temperatures reduce the uniformity of the pore size distribution(PSD) in the bedding-parallel direction but narrow the PSD in the bedding-perpendicular direction. The greater ff(> 0.7) values in the bedding-parallel direction account for a large proportion, while the dominated in the bedding-perpendicular direction locates within 0.2-0.7, for all temperatures. The H value of the bedding-parallel sample remains stable at ~0.925 during heating, but gradually increases from 0.808 at 25℃ to 0.879 at 500℃ for the beddingperpendicular sample. Congruent with a mechanistic model, the permeability at 500℃ is elevated~1.83 times(bedding-parallel) and ~6.08 times(bedding-perpendicular) relative to that at 25℃—confirming the effectiveness of thermal treatment in potentially enhancing production from low-maturity oil shales.展开更多
This paper studied the rock dynamic fracture propagation under impact loads elaborately with a determination method proposed to calculate crack propagation dynamic stress intensity factor(DSIF).By utilizing the split-...This paper studied the rock dynamic fracture propagation under impact loads elaborately with a determination method proposed to calculate crack propagation dynamic stress intensity factor(DSIF).By utilizing the split-Hopkinson pressure bar,the impact experiments with an improved single cleavage semi-circle(ISCSC)specimen were conducted to illuminate the dynamic crack propagation behaviour.Meanwhile,the fracture characteristics and crack propagation velocity were obtained by the crack propagation gauges.Coordinating experiments with a numerical approach,the crack propagation dynamic stress intensity factors were calculated by an experimental—numerical method with fractal theory.Then,a finite difference model was developed based on the tensile fracture softening damage criterion.With the analysis of numerical and experimental results,the crack propagation behaviour and mechanism of crack arrest were discussed sophisticatedly.The results demonstrate that the novel ISCSC specimen shows a definite advantage in determining crack propagation and arrest DSIF.Additionally,the crack arrest DSIF is larger than the average propagation DSIF with a sharp increase.Meanwhile,the numerical simulation results which agree well with the actual crack propagation illustrate that the crack arrest should be dominated by the compressive stress perpendicular to the crack path,and there were several arrest pauses existing in the transitory crack arrest process.展开更多
In this work,a three-dimensional crystal-plasticity-based phase-field model considering three kinds of inelastic deformation mechanisms,i.e.,martensitic transformation,dislocation slip in austenite,and dislocation sli...In this work,a three-dimensional crystal-plasticity-based phase-field model considering three kinds of inelastic deformation mechanisms,i.e.,martensitic transformation,dislocation slip in austenite,and dislocation slip in martensite,is established to simulate the stress-assisted two-way shape memory effect(SATWSME)of NiTi single crystals and its cyclic degradation.The simulation results show that the ability of the SATWSME of NiTi single crystal increases as increasing the constant stress in the range discussed in this work(10–100 MPa),which is due to the increase of reoriented martensite formed in the cooling process due to the enhanced variant-selection capability of increased constant stress.The martensitic transformation and its reverse in the cyclic process reflecting the SATWSME show more and more obvious localization characteristics,resulting in the accumulation of significantly heterogeneous plastic deformation(mainly caused by the dislocation slip in austenite),which leads to the cyclic degradation of SATWSME.The simulation results and the conclusions drawn from this work are helpful for further understanding the mechanism of functional cyclic degradation of NiTi alloys.展开更多
It is well known that coarse-grained super-elastic NiTi shape memory alloys(SMAs)exhibit localized rather than homogeneous martensite transformation(MT),which,however,can be strongly influenced by either internal size...It is well known that coarse-grained super-elastic NiTi shape memory alloys(SMAs)exhibit localized rather than homogeneous martensite transformation(MT),which,however,can be strongly influenced by either internal size(grain size,GS)or the external size(geometric size).The coupled effect of GS and geometric size on the functional properties has not been clearly understood yet.In this work,the super-elasticity,one-way,and stress-assisted two-way shape memory effects of the polycrystalline NiTi SMAs with different aspect ratios(length/width for the gauge section)and different GSs are investigated based on the phase field method.The coupled effect of the aspect ratio and GS on the functional properties is adequately revealed.The simulated results indicate that when the aspect ratio is lower than about 4:1,the stress biaxiality and stress heterogeneity in the gauge section of the sample become more and more obvious with decreasing the aspect ratio,which can significantly influence the microstructure evolution in the process involving external stress.Therefore,the corresponding functional property is strongly dependent on the aspect ratio.With decreasing the GS and the aspect ratio(to be lower than 4:1),both the aspect ratio and GS can affect the MT or martensite reorientation in each grain and the interaction among grains.Thus,due to the strong internal constraint(i.e.,the constraint of grain boundary)and the external constraint(i.e.,the constraint of geometric boundary),the capabilities of the functional properties of NiTi SMAs are gradually weakened and highly dependent on these two factors.展开更多
Aiming at obtaining an accurate porosity of gas shale,various approaches are attempted.Therein,nuclear magnetic resonance(NMR),being treated as a kind of new-developed technique,possesses the representative significan...Aiming at obtaining an accurate porosity of gas shale,various approaches are attempted.Therein,nuclear magnetic resonance(NMR),being treated as a kind of new-developed technique,possesses the representative significance.However,as a booming technique,the reliability of NMR-based porosity of shale is not exactly defined.Depending on NMR device,this work measured the porosity of shale experiencing different water soaking time,accordingly,judging the reliability of NMR-based porosity.Results indicate the NMR outcomes vary with the water soaking time,making a doubt about the objectivity of NMRbased porosity in reflecting the real shale porosity.Furthermore,some supplementary means were adopted to verify the water soaking-induced variation in the pore system of shale sample,which intensities the suspicion if the NMR-based porosity is reliable or not.To sum up,this work considers that the NMR-based porosity of shale is not reliable enough when water is used as the probe.Besides,this work also offers some suggestions on how to enhance the reliability of NMR-based porosity of shale sample.Basically,this work selects a fresh perspective to analyze the NMR approach in determining shale porosity,which is hopefully helpful in promoting the development of NMR technique in the shale-related field.展开更多
Tunnels in fractured rock masses are typically damaged by dynamic disturbances from various directions.To investigate the influence of blasting load directions on the stability of a tunnel with a precrack nearby,blast...Tunnels in fractured rock masses are typically damaged by dynamic disturbances from various directions.To investigate the influence of blasting load directions on the stability of a tunnel with a precrack nearby,blasting tests were conducted on the physical models of an external crack around a tunnel(ECT)in this study.Failure modes of the tunnels were analysed based on stress wave theory.The Riedel-Hiermaier-Thoma(RHT)material model was employed to perform the numerical simulations on ECT models.Stress distribution around the tunnels and final failure patterns of the tunnels were characterised.The results show that,under blasting loads,the pre-crack propagates and then new cracks initiates on the incident side of the tunnel.These cracks extend towards each other and eventually coalesce.Blasting load directions significantly influence the ultimate failure mode of the tunnel in the fractured rock masses.The new cracks on the shadow side of the tunnel appear at different positions when the blasting stress waves come from various directions.The results are meaningful to the analysis of tunnel stability and optimisation of the tunnel support scheme.展开更多
Dynamic crack propagation in brittle materials plays an important role in unJerstanding the fracture mechanism.Numerical simulations on crack propagation in the polymethyle methacrylate(PMMA)under explosive loads were...Dynamic crack propagation in brittle materials plays an important role in unJerstanding the fracture mechanism.Numerical simulations on crack propagation in the polymethyle methacrylate(PMMA)under explosive loads were carried out through Autodyn 3D code.Although the Johnson-Holmquist(JH-2)constitutive model has been widely treated as an acceptable scheme,it cannot improve the description in the post-failure response of material and reduce the mesh dependency,so the crack softening failure model are introduced in our present work based on two failure criterions.First,material parameters of JH-2 model,failure criterion,and crack softening failure model are determined from available data and calculations.The circular and rectangle thin plates are modeled to explore the fracture mechanisms for single-borehole and dual-borehole explosions.The simulation results well reproduced the entire dynamic evolutionary process of the crushed and fractured zones,crack initiation,propagation and arrest as well as secondary propagation,which successfully proved the relibilities of the combination of JH-2 constitutive model,failure criterion,and crack softening failure model and corresponding material parameters.For dual-borehole explosion,the crack linkage is well performed when borehole distance L is 20 and 30 cm;cracks are failed to link with each other at L=40 cm although main cracks have arrived into the opposite fractured zone.More importantly,it is found that the crack linkage mainly depends on L,which has an important effect on linkage style and its location.Despite a larger L makes more difficult in crack linkage,this difficulty lies not in the shortage of crack length,but in losing control of directional fracture of cracks between two boreholes.展开更多
Rock structures are often subjected to dynamic loads,such as blasts,impacts and earthquakes,and their loading rates differ largely.To investigate the effect of loading rates on the dynamic behavior of crack propagatio...Rock structures are often subjected to dynamic loads,such as blasts,impacts and earthquakes,and their loading rates differ largely.To investigate the effect of loading rates on the dynamic behavior of crack propagation,impact tests were conducted on large single-cleavage semicircle compression(LSCSC)specimens using a drop weight impact test system.Five types of rock materials were selected to prepare the LSCSC specimens,and crack propagation gauges were mounted along the crack propagation paths to measure crack initiation time and propagation speeds.Finite element models were established by using ABAQUS code,and the dynamic stress intensity factors(SIFs)were calculated.The curves of dynamic SIFs versus time were obtained,and the initiation toughness was determined by using these curves and the initiation time measured in the impact tests.The results show that loading rate has a significant effect on crack propagation behavior,and both the crack propagation speed and initiation toughness increase with the loading rate,whereas the delayed fracture time decreases with the increase in loading rate.展开更多
Dynamic fractures occur frequently in geophysical processes and engineering applications.It is thus essential to study crack and failure behaviors,such as crack time-to-initiation,crack growth rate and arrest period u...Dynamic fractures occur frequently in geophysical processes and engineering applications.It is thus essential to study crack and failure behaviors,such as crack time-to-initiation,crack growth rate and arrest period under dynamic loading.In this study,impact experiments were implemented by utilizing the single-flawed tunnel specimens under drop-hammer impacts.Four brittle materials,i.e.,green sandstone,red sandstone,black sandstone and polymethyl methacrylate,were selected to make single-flawed tunnel specimens.Strain gauges and crack extension gauges were employed to measure the crack extension parameters.The properties of crack growth rate,crack time-to-initiation and arrest period of these four brittle materials were discussed and analyzed.The corresponding numerical simulation was performed by using the commercial software AUTODYN.The numerical results of crack growth rate and crack time-to-initiation agreed with the impact test results.The commercial software ABAQUS was applied to compute the dynamic stress intensity factors.The results show that both the dynamic initiation fracture toughness and the crack growth rate increase with the elastic moduli of these four types of brittle materials under the same loading conditions,whereas the crack time-to-initiation decreases with the increase in elastic moduli of the brittle materials under the same loading conditions.展开更多
Stress waves affect the stress field at the crack tip and dominate the dynamic crack propagation.Therefore,evaluating the influence of blasting stress waves on the crack propagation behavior and the mechanical charact...Stress waves affect the stress field at the crack tip and dominate the dynamic crack propagation.Therefore,evaluating the influence of blasting stress waves on the crack propagation behavior and the mechanical characteristics of crack propagation is of great significance for engineering blasting.In this study,ANSYS/LS-DYNA was used for blasting numerical simulation,in which the propagation characteristics of blasting stress waves and stress field distribution at the crack tip were closely observed.Moreover,ABAQUS was applied for simulating the crack propagation path and calculating dynamic stress intensity factors(DSIFs).The universal function was calculated by the fractalmethod.The results show that:the compressive wave causes the crack to close and the reflected tensile wave drives the crack to initiate and propagate,and failure mode is mainly tensile failure.The crack propagation velocity varies with time,which increases at first and then decreases,and the crack arrest occurs due to the attenuation of stress waves and dissipation of the blasting energy.In addition,crack arrest toughness is smaller than the crack initiation toughness,applied pressure waveforms(such as the peak pressure,duration,waveforms,wavelengths and loading rates)have a great influence on DSIFs.It is conducive to our deep understanding or the study of blasting stress waves dominated fracture,suggesting a broad reference for the further development of rock blasting in engineering practice.展开更多
基金financially supported by National Natural Science Foundation of China(Nos.U2013603,52004167,52078477,51827901,and U1965203)Sichuan International Technological innovation Cooperation Project(No.2018HH0159)。
文摘During the development of hot dry rock,the research on thermal fatigue damage caused by thermal shock of cold and heat cycles is the basis that ensures the long-term utilization of geothermal resources,but there are not enough relevant studies at present.Based on this,the thermal damage tests of granite at different temperatures(250,350,450°C)and quenching cycles(1,5,10,15 cycles)were carried out.Preliminary reveals the damage mechanism and heat transfer law of the quenching cycle effect on hot dry rock.The results show that with the increase of temperature and cycles,the uneven thermal expansion of minerals and the thermal shock caused by quenching promote the crack development of granite,resulting in the decrease of P-wave velocity,thermal conductivity and uniaxial compressive strength of granite.Meanwhile,the COMSOL was used to simulate the heat transfer of hot dry rock under different heat treatment conditions.It concluded that the increase in the number of quenching cycles reduced the heat transfer capacity of the granite,especially more than 10 quenching cycles,which also reflects that the thermal fatigue damage leads to a longer time for the temperature recovery of the hot dry rock mass.In addition,the three-dimensional nonlinear fitting relationship among thermal conductivity,temperature and cycle number was established for the first time,which can better reveal the change rule of thermal conductivity after quenching thermal fatigue effect of hot dry rock.The research results provide theoretical support for hot dry rock reservoir reconstruction and production efficiency evaluation.
基金Projects(11672194,U19A2098)supported by the National Natural Science Foundation of ChinaProject(2018SCU12047)supported by Fundamental Research Funds for the Central Universities,ChinaProject(2018JZ0036)supported by the Project of Science and Technology of Sichuan Province,China。
文摘The former studies indicate that loading rates significantly affect dynamic behavior of brittle materials,for instance,the dynamic compressive and tensile strength increase with loading rates.However,there still are many unknown or partially unknown aspects.For example,whether loading rates have effect on crack dynamic propagating behavior(propagation toughness,velocity and arrest,etc).To further explore the effect of loading rates on crack dynamic responses,a large-size single-cleavage trapezoidal open(SCTO)specimen was proposed,and impacting tests using the SCTO specimen under drop plate impact were conducted.Crack propagation gauges(CPGs)were employed in measuring impact loads,crack propagation time and velocities.In order to verify the testing result,the corresponding numerical model was established using explicit dynamic software AUTODYN,and the simulation result is basically consistent with the experimental results.The ABAQUS software was used to calculate the dynamic SIFs.The universal function was calculated by fractal method.The experimental-numerical method was employed in determining initiation toughness and propagation toughness.The results indicate that crack propagating velocities,dynamic fracture toughness and energy release rates increase with loading rates;crack delayed initiation time decreases with loading rates.
基金Sichuan Transportation Science and Technology Project(Grant Nos.2020-MS3-101/2020-B-01 and 2019-ZL-12 and 2018-B-03)the Science and Technology Department of Sichuan Province(Nos.2021YFH0048 and 2021YFH0118).
文摘The prevention and the reduction of the rockfall are the common measures of the prevention and the reduction of disasters.When the rock-shed resists the impact of the rockfall,the force that acts on the structure consists of the cushion dead load and the impact-induced load,of which the dynamic process of the propagation of the impactinduced load is complex.Therefore,we conducted a numerical study to investigate the impact of the rockfall.Considering the highly discrete characteristic of the sand,we developed a numerical model on the basis of the discrete element method(DEM).The numerical model,which simulation results were validated by the results of real-scale experiments,was used to investigate the dynamic response of the impact force of the rockfall and the transmission of the impact force under the different magnitude of the falling height and the different thickness of the sand cushion.The results of our study indicated that the cushion thickness had little effect on the impact of the rockfall,and the dense sand cushion generated higher impact force than did the loose sand cushion.Although the high thickness enhanced the buffer performance of the sand cushion,the additional force induced by the dead load of sand cushion was significant.Therefore,to determine the appropriate thickness of the sand cushion,we suggested designers consider the buffer performance and the dead load of the sand cushion.The analysis presented in this paper provided a practical estimation of the impact-induced force of the thick sand cushion.
基金supported by the Program for Guangdong Introducing Innovative and Entrepreneurial Teams(No.2019ZT08G315)National Natural Science Foundation of China No.51827901 and U2013603Shenzhen Basic Research Project(JCYJ20190808153416970)
文摘Natural gas hydrate(NGH)holds great promise as a source of clean energy.It is critical for acquiring the largest possible in situ NGH core for NGH eigen features and resource assessment.However,the existing NGH coring technology has limitations,such as temperature increments,limited coring diameters,low coring rates,and complex coring structures.Therefore,this study designs and proposes an NGH freezing coring(NGHFC)method and verifies the freezing and coring capacities of the NGHFC method in laboratories and experimental wells.Results suggest that NGHFC shows good freezing and heat-retention properties.A freezing core heat transfer model is developed.According to the actual air temperature and operating time,the optimum initial temperature of the cold source can be determined using this model.The average coring rate of NGHFC can reach 77.86%.The research results will provide a new idea of coring gas hydrates.
基金The authors are grateful for the financial support from the National Natural Science Foundation of China(No.51827901&No.52274133)the Program for Guangdong Introducing Innovative and Enterpreneurial Teams(No.2019ZT08G315)the Shenzhen National Science Fund for Distinguished Young Scholars(RCJC20210706091948015).
文摘In-situ pressure-preserved coring(IPP-Coring)is considered to be the most reliable and efficient method for the identification of the scale of oil and gas resources.During IPP-Coring,because the rotation behavior of the pressure controller valve cover in different medium environments is unclear,interference between the valve cover and inner pipe may occur and negatively affect the IPP-Coring success rate.To address this issue,we conducted a series of indoor experiments employing a high-speed camera to gain greater insights into the valve cover rotation behavior in different medium environments,e.g.,air,water,and simulated drilling fluids.The results indicated that the variation in the valve cover rotation angle in the air and fluid environments can be described by a one-phase exponential decay function with a constant time parameter and by biphasic dose response function,respectively.The rotation behavior in the fluid environments exhibited distinct elastic and gravitational acceleration zones.In the fluid environments,the density clearly impacted the valve cover closing time and rotation behavior,whereas the effect of viscosity was very slight.This can be attributed to the negligible influence of the fluid viscosity on the drag coefficient found in this study;meanwhile,the density can increase the buoyancy and the time period during which the valve cover experienced a high drag coefficient.Considering these results,control schemes for the valve cover rotation behavior during IPP-Coring were proposed for different layers and geological conditions in which the different drilling fluids should be used,e.g.,the use of a high-density valve cover in high-pore pressure layers.
基金would like to extend their sincere gratitude for the financial support from the Integrated Key Precast Components and New Wood-bamboo Composite Structure Foundation of China(2017YFC0703502)the National Natural Science Foundation of China(51978152)the Fundamental Research Funds for the Central Universities(YJ202061).
文摘Cross-laminated bamboo(CLB)have a high strength to weight ratio and stable bidirectional mechanical properties.Inspired by the investigation on cross-laminated timber(CLT)rocking walls,CLB rocking walls with conventional friction dampers(CFDs)are studied in this paper.To investigate the mechanical properties of the CLB rocking wall,seven tests are conducted under a cyclic loading scheme,and different test parameters,including the existence of the CFDs,the moment ratio,and the loading times,are discussed.The test results show a bilinear behavior of the CLB rocking wall.The small residual displacements of the CLB rocking wall demonstrate an idealized self-centering capacity.The cumulative energy dissipation curves indicate that the energy dissipation capacity of the CLB rocking wall can be greatly improved with CFDs.The limit states of the CLB rocking wall under a lateral force are proposed based on the strains,stress,and damage level of the CLB material and posttensioned rebar.In addition,an analytical model of the CLB rocking wall is developed based on the proposed limit states of the CLB rocking wall to evaluate the hysteretic response of the CLB rocking wall,and the model is validated by the experimental data.The comparison results show the potential value of the analytical model for engineering design.
基金financially supported by the National Key Research and Development Program of China(Grant No.2022YFE0129800)the National Natural Science Foundation of China(Grant No.42202204)。
文摘In-situ upgrading by heating is feasible for low-maturity shale oil,where the pore space dynamically evolves.We characterize this response for a heated substrate concurrently imaged by SEM.We systematically follow the evolution of pore quantity,size(length,width and cross-sectional area),orientation,shape(aspect ratio,roundness and solidity)and their anisotropy—interpreted by machine learning.Results indicate that heating generates new pores in both organic matter and inorganic minerals.However,the newly formed pores are smaller than the original pores and thus reduce average lengths and widths of the bedding-parallel pore system.Conversely,the average pore lengths and widths are increased in the bedding-perpendicular direction.Besides,heating increases the cross-sectional area of pores in low-maturity oil shales,where this growth tendency fluctuates at<300℃ but becomes steady at>300℃.In addition,the orientation and shape of the newly-formed heating-induced pores follow the habit of the original pores and follow the initial probability distributions of pore orientation and shape.Herein,limited anisotropy is detected in pore direction and shape,indicating similar modes of evolution both bedding-parallel and bedding-normal.We propose a straightforward but robust model to describe evolution of pore system in low-maturity oil shales during heating.
基金supported by the Program for Guangdong Introducing Innovative and Entrepreneurial Teams(No.2019ZT08G315)Shenzhen Basic Research Program(General Program,No.JCYJ20190808153416970)National Natural Science Foundation of China No.51827901
文摘The sealing performance of contact interfaces plays the most important role in the design and operation of the in-situ pressure-preserved coring system.To meet the demand of ultra-high pressure-retained coring for oil and gas exploration in deep reservoirs,a quantitative analysis of the contact mechanical behavior of the pressure controller was performed.Based on the micro-contact theory of rough surfaces,a three-dimensional numerical model of the rough contact interface between the valve cover and the valve seat was constructed,and the micro-contact behavior of the metal contact surfaces was comprehensively studied.The results show that the actual contact area of the valve interface increases with the increase of surface roughness before the critical contact point,but decreases after that.Compared with the real contact model with double rough surfaces,although the simplified hard-contact model with a single rough surface can reflect the micro-contact behavior of the rough surface to a certain extent,it cannot truly reveal the microchannel morphology between the sealing interfaces under pressure.Therefore,the realistic double-rough-surface model should be recommended to evaluate the sealing performance of coring tools,particularly for high pressure conditions.The material properties of valves have a significant effect on the contact characteristics of rough surfaces,which suggested that the actual contact area decreases with the increase of the elastic modulus of the contact material under the same loading conditions.The knowledge of this work could help to enhance the seal design of pressure controllers for in-situ pressure-preserved coring.
基金The National Natural Science Foundation of China(12202294 and 12022208)the Project funded by China Postdoctoral Science Foundation(2022M712243)the Fundamental Research Funds for the Central Universities(2023SCU12098)are acknowledged.
文摘Texture is inevitably introduced during the manufacturing of most NiTi shape memory alloys(SMAs),and the textured nanocrystalline NiTi has been extensively employed in engineering.However,the effect of texture,and the joint effect of grain size(GS)and texture on the functional properties of NiTi SMAs and the corresponding microscopic mechanisms have not been clearly understood yet.In this work,based on the phase field method,the effect of texture on the GS-dependent functional properties of NiTi SMAs,including super-elasticity(SE),one-way shape memory effect(OWSME),and stress-assisted two-way shape memory effect(SATWSME),is investigated,and the corresponding microscopic mechanisms are revealed.Moreover,the samples with discrete geometrical gradients and/or texture gradients are designed to achieve graded functional properties.The simulation results indicate that the dependence of functional properties on texture is due to the effect of crystallographic orientation on martensite transformation and reorientation,which can lead to different inelastic strains.In the designed samples with texture gradients,the stress–strain responses of sheets with various textures are different,allowing for the coordination of overall deformation of the sample by combining such sheets,with varying inelastic deformation degrees.Thus,the overall response of the sample differs from that without texture gradient,leading to the achievement of graded functional properties.The simulation results and new findings in this work contribute to a deeper understanding of the effects of texture,GS,and their interaction on the functional properties of SMAs,and provide valuable reference for the design and development of SMA-based devices with desired functional properties.
基金The National Natural Science Foundation of China(12022208)the Project funded by China Postdoctoral Science Foundation(2022M712243)the Fundamental Research Funds for the Cen-tral Universities are acknowledged.
文摘The inelastic deformations of shape memory alloys(SMAs)always show poor controllability due to the avalanche-like martensite transformation,and the effective control for the deformation of precision de-vices has been not yet mature.In this work,the phase field method was used to investigate the shape memory effects(SMEs)of NiTi SMAs undergoing grain size(GS)engineering,to obtain tunable one-way and stress-assisted two-way SMEs(OWSME and SATWSME).The OWSME and SATWSME of the systems with various gradient-nanograin structures and bimodal grain structure,as well as that with geometric gradients were simulated.The simulated results indicate that due to the GS dependences of martensite transformation and reorientation,the occurrence and expansion of martensite reorientation,martensite transformation and its reverse can be efficaciously controlled via the GS engineering.When combining the GS engineering and geometric gradient design,since the effects of GS and stress gradient can be su-perimposed or competing,and the responses of martensite reorientation,martensite transformation and its reverse to this are different,the OWSME and SATWSME of the geometrically graded systems with various nanograin structures can exhibit different improvements in controllability.In short,the reorienta-tion hardening modulus during OWSME is increased and the transformation temperature window during SATWSME is widened by GS engineering,indicating the improved controllability of SMEs.The optimal GS engineering schemes revealed in this work provide the basic reference and guidance for designing tun-able SMEs and producing NiTi-based driving devices catering to desired functional performance in various engineering fields.
基金financially supported by the National Key Research and Development Program of China (Grant No. 2022YFE0129800)the National Natural Science Foundation of China (Grant No. 42202204)support from the G. Albert Shoemaker endowment。
文摘In-situ thermal upgrading is used to tune the pore system in low-maturity oil shales. We introduce fractal dimension(D), form factor(ff) and stochastic entropy(H) to quantify the heating-induced evolution of pore morphological complexity and azimuthal disorder and develop a model to estimate the impact on seepage capacity via permeability. Experiments are conducted under recreated in-situ temperatures and consider anisotropic properties—both parallel and perpendicular to bedding. Results indicate that azimuthal distribution of pores in the bedding-parallel direction are dispersed, while those in the bedding-perpendicular direction are concentrated. D values indicate that higher temperatures reduce the uniformity of the pore size distribution(PSD) in the bedding-parallel direction but narrow the PSD in the bedding-perpendicular direction. The greater ff(> 0.7) values in the bedding-parallel direction account for a large proportion, while the dominated in the bedding-perpendicular direction locates within 0.2-0.7, for all temperatures. The H value of the bedding-parallel sample remains stable at ~0.925 during heating, but gradually increases from 0.808 at 25℃ to 0.879 at 500℃ for the beddingperpendicular sample. Congruent with a mechanistic model, the permeability at 500℃ is elevated~1.83 times(bedding-parallel) and ~6.08 times(bedding-perpendicular) relative to that at 25℃—confirming the effectiveness of thermal treatment in potentially enhancing production from low-maturity oil shales.
基金the National Natural Science Foundation of China(Grant Nurmbers 11702181,11672194)the Sichuan Science and Technology Program(Grant Numbers 2019YFG0047).
文摘This paper studied the rock dynamic fracture propagation under impact loads elaborately with a determination method proposed to calculate crack propagation dynamic stress intensity factor(DSIF).By utilizing the split-Hopkinson pressure bar,the impact experiments with an improved single cleavage semi-circle(ISCSC)specimen were conducted to illuminate the dynamic crack propagation behaviour.Meanwhile,the fracture characteristics and crack propagation velocity were obtained by the crack propagation gauges.Coordinating experiments with a numerical approach,the crack propagation dynamic stress intensity factors were calculated by an experimental—numerical method with fractal theory.Then,a finite difference model was developed based on the tensile fracture softening damage criterion.With the analysis of numerical and experimental results,the crack propagation behaviour and mechanism of crack arrest were discussed sophisticatedly.The results demonstrate that the novel ISCSC specimen shows a definite advantage in determining crack propagation and arrest DSIF.Additionally,the crack arrest DSIF is larger than the average propagation DSIF with a sharp increase.Meanwhile,the numerical simulation results which agree well with the actual crack propagation illustrate that the crack arrest should be dominated by the compressive stress perpendicular to the crack path,and there were several arrest pauses existing in the transitory crack arrest process.
基金The National Natural Science Foundation of China(12202294,12022208,and 12322203)the Project funded by China Postdoctoral Science Foundation(2022M712243)the Fundamental Research Funds for the Central Universities(2023SCU12098)are acknowledged.
文摘In this work,a three-dimensional crystal-plasticity-based phase-field model considering three kinds of inelastic deformation mechanisms,i.e.,martensitic transformation,dislocation slip in austenite,and dislocation slip in martensite,is established to simulate the stress-assisted two-way shape memory effect(SATWSME)of NiTi single crystals and its cyclic degradation.The simulation results show that the ability of the SATWSME of NiTi single crystal increases as increasing the constant stress in the range discussed in this work(10–100 MPa),which is due to the increase of reoriented martensite formed in the cooling process due to the enhanced variant-selection capability of increased constant stress.The martensitic transformation and its reverse in the cyclic process reflecting the SATWSME show more and more obvious localization characteristics,resulting in the accumulation of significantly heterogeneous plastic deformation(mainly caused by the dislocation slip in austenite),which leads to the cyclic degradation of SATWSME.The simulation results and the conclusions drawn from this work are helpful for further understanding the mechanism of functional cyclic degradation of NiTi alloys.
基金supported by the National Natural Science Foundation of China (Grant Nos.12202294 and 12022208)the Project funded by China Postdoctoral Science Foundation (Grant No.2022M712243)the Fundamental Research Funds for the Central Universities (Grant No.2023SCU12098).
文摘It is well known that coarse-grained super-elastic NiTi shape memory alloys(SMAs)exhibit localized rather than homogeneous martensite transformation(MT),which,however,can be strongly influenced by either internal size(grain size,GS)or the external size(geometric size).The coupled effect of GS and geometric size on the functional properties has not been clearly understood yet.In this work,the super-elasticity,one-way,and stress-assisted two-way shape memory effects of the polycrystalline NiTi SMAs with different aspect ratios(length/width for the gauge section)and different GSs are investigated based on the phase field method.The coupled effect of the aspect ratio and GS on the functional properties is adequately revealed.The simulated results indicate that when the aspect ratio is lower than about 4:1,the stress biaxiality and stress heterogeneity in the gauge section of the sample become more and more obvious with decreasing the aspect ratio,which can significantly influence the microstructure evolution in the process involving external stress.Therefore,the corresponding functional property is strongly dependent on the aspect ratio.With decreasing the GS and the aspect ratio(to be lower than 4:1),both the aspect ratio and GS can affect the MT or martensite reorientation in each grain and the interaction among grains.Thus,due to the strong internal constraint(i.e.,the constraint of grain boundary)and the external constraint(i.e.,the constraint of geometric boundary),the capabilities of the functional properties of NiTi SMAs are gradually weakened and highly dependent on these two factors.
基金financially supported by the Science and Technology Department of Sichuan Province(Grant Nos.2021YFH0048 and 2021YFH0118)the Project funded by China Postdoctoral Science Foundation(Grant No.2020M683253)
文摘Aiming at obtaining an accurate porosity of gas shale,various approaches are attempted.Therein,nuclear magnetic resonance(NMR),being treated as a kind of new-developed technique,possesses the representative significance.However,as a booming technique,the reliability of NMR-based porosity of shale is not exactly defined.Depending on NMR device,this work measured the porosity of shale experiencing different water soaking time,accordingly,judging the reliability of NMR-based porosity.Results indicate the NMR outcomes vary with the water soaking time,making a doubt about the objectivity of NMRbased porosity in reflecting the real shale porosity.Furthermore,some supplementary means were adopted to verify the water soaking-induced variation in the pore system of shale sample,which intensities the suspicion if the NMR-based porosity is reliable or not.To sum up,this work considers that the NMR-based porosity of shale is not reliable enough when water is used as the probe.Besides,this work also offers some suggestions on how to enhance the reliability of NMR-based porosity of shale sample.Basically,this work selects a fresh perspective to analyze the NMR approach in determining shale porosity,which is hopefully helpful in promoting the development of NMR technique in the shale-related field.
基金funded by the National Natural Science Foundation of China(Grant No.U19A2098)the open fund of MOE Key Laboratory of Deep Underground Science and Engineering(Grant No.DESEYU202101)the Sichuan Science and Technology Program(Grant No.2021YJ0511)。
文摘Tunnels in fractured rock masses are typically damaged by dynamic disturbances from various directions.To investigate the influence of blasting load directions on the stability of a tunnel with a precrack nearby,blasting tests were conducted on the physical models of an external crack around a tunnel(ECT)in this study.Failure modes of the tunnels were analysed based on stress wave theory.The Riedel-Hiermaier-Thoma(RHT)material model was employed to perform the numerical simulations on ECT models.Stress distribution around the tunnels and final failure patterns of the tunnels were characterised.The results show that,under blasting loads,the pre-crack propagates and then new cracks initiates on the incident side of the tunnel.These cracks extend towards each other and eventually coalesce.Blasting load directions significantly influence the ultimate failure mode of the tunnel in the fractured rock masses.The new cracks on the shadow side of the tunnel appear at different positions when the blasting stress waves come from various directions.The results are meaningful to the analysis of tunnel stability and optimisation of the tunnel support scheme.
基金the Key Research Foundation of Education Bureau of Sichuan Province(Grant No.17zd1122)the Natural Science Foundation of Southwest University of Science and Technology(Grant No.19zx7168)+2 种基金the Shock and Vibration of Engineering Materials and Struc-tures Key Laboratory of Sichuan Province(Grant No.I9kfgk07)the Sci-ence and Technology Department of Sichuan Province(Grant No.2021 YJ0525)National Natural Science Foundation of China(Grant No.11802255).
文摘Dynamic crack propagation in brittle materials plays an important role in unJerstanding the fracture mechanism.Numerical simulations on crack propagation in the polymethyle methacrylate(PMMA)under explosive loads were carried out through Autodyn 3D code.Although the Johnson-Holmquist(JH-2)constitutive model has been widely treated as an acceptable scheme,it cannot improve the description in the post-failure response of material and reduce the mesh dependency,so the crack softening failure model are introduced in our present work based on two failure criterions.First,material parameters of JH-2 model,failure criterion,and crack softening failure model are determined from available data and calculations.The circular and rectangle thin plates are modeled to explore the fracture mechanisms for single-borehole and dual-borehole explosions.The simulation results well reproduced the entire dynamic evolutionary process of the crushed and fractured zones,crack initiation,propagation and arrest as well as secondary propagation,which successfully proved the relibilities of the combination of JH-2 constitutive model,failure criterion,and crack softening failure model and corresponding material parameters.For dual-borehole explosion,the crack linkage is well performed when borehole distance L is 20 and 30 cm;cracks are failed to link with each other at L=40 cm although main cracks have arrived into the opposite fractured zone.More importantly,it is found that the crack linkage mainly depends on L,which has an important effect on linkage style and its location.Despite a larger L makes more difficult in crack linkage,this difficulty lies not in the shortage of crack length,but in losing control of directional fracture of cracks between two boreholes.
基金the National Natural Science Foundation of China(1167219411702181)+1 种基金by Sichuan Administration of Work Safety(aj20170515161307)the project of Science and Technology of Sichuan province(2018JZ0036).
文摘Rock structures are often subjected to dynamic loads,such as blasts,impacts and earthquakes,and their loading rates differ largely.To investigate the effect of loading rates on the dynamic behavior of crack propagation,impact tests were conducted on large single-cleavage semicircle compression(LSCSC)specimens using a drop weight impact test system.Five types of rock materials were selected to prepare the LSCSC specimens,and crack propagation gauges were mounted along the crack propagation paths to measure crack initiation time and propagation speeds.Finite element models were established by using ABAQUS code,and the dynamic stress intensity factors(SIFs)were calculated.The curves of dynamic SIFs versus time were obtained,and the initiation toughness was determined by using these curves and the initiation time measured in the impact tests.The results show that loading rate has a significant effect on crack propagation behavior,and both the crack propagation speed and initiation toughness increase with the loading rate,whereas the delayed fracture time decreases with the increase in loading rate.
基金This work was financially supported by the National Natural Science Foundation of China(11672194)Sichuan Administration of Work Safety(aj20170515161307)the Project of Science and Technology of Sichuan Province(2018JZ0036).
文摘Dynamic fractures occur frequently in geophysical processes and engineering applications.It is thus essential to study crack and failure behaviors,such as crack time-to-initiation,crack growth rate and arrest period under dynamic loading.In this study,impact experiments were implemented by utilizing the single-flawed tunnel specimens under drop-hammer impacts.Four brittle materials,i.e.,green sandstone,red sandstone,black sandstone and polymethyl methacrylate,were selected to make single-flawed tunnel specimens.Strain gauges and crack extension gauges were employed to measure the crack extension parameters.The properties of crack growth rate,crack time-to-initiation and arrest period of these four brittle materials were discussed and analyzed.The corresponding numerical simulation was performed by using the commercial software AUTODYN.The numerical results of crack growth rate and crack time-to-initiation agreed with the impact test results.The commercial software ABAQUS was applied to compute the dynamic stress intensity factors.The results show that both the dynamic initiation fracture toughness and the crack growth rate increase with the elastic moduli of these four types of brittle materials under the same loading conditions,whereas the crack time-to-initiation decreases with the increase in elastic moduli of the brittle materials under the same loading conditions.
基金This researchwas supported by the National Natural Science Foundation of China(No.52227805)the Fundamental Research Funds for Central Universities(No.2022JCCXLJ01).Awards were granted to the author Liyun Yang.
文摘Stress waves affect the stress field at the crack tip and dominate the dynamic crack propagation.Therefore,evaluating the influence of blasting stress waves on the crack propagation behavior and the mechanical characteristics of crack propagation is of great significance for engineering blasting.In this study,ANSYS/LS-DYNA was used for blasting numerical simulation,in which the propagation characteristics of blasting stress waves and stress field distribution at the crack tip were closely observed.Moreover,ABAQUS was applied for simulating the crack propagation path and calculating dynamic stress intensity factors(DSIFs).The universal function was calculated by the fractalmethod.The results show that:the compressive wave causes the crack to close and the reflected tensile wave drives the crack to initiate and propagate,and failure mode is mainly tensile failure.The crack propagation velocity varies with time,which increases at first and then decreases,and the crack arrest occurs due to the attenuation of stress waves and dissipation of the blasting energy.In addition,crack arrest toughness is smaller than the crack initiation toughness,applied pressure waveforms(such as the peak pressure,duration,waveforms,wavelengths and loading rates)have a great influence on DSIFs.It is conducive to our deep understanding or the study of blasting stress waves dominated fracture,suggesting a broad reference for the further development of rock blasting in engineering practice.
