• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Effects of thawing-induced softening on fracture behaviors of frozen rock

    2024-03-25 11:06:26TingWngHilingJiQingSunXinjunTnLiyunTng

    Ting Wng,Hiling Ji,Qing Sun,Xinjun Tn,Liyun Tng

    a College of Geology and Environment, Xi’an University of Science and Technology, Xi’an, 710054, China

    b College of Architecture and Civil Engineering, Xi’an University of Science and Technology, Xi’an, 710054, China

    c State Key Laboratory of Geomechanics and Geotechnical Engineering,Institute of Rock and Soil Mechanics,Chinese Academy of Sciences,Wuhan,430064,China

    Keywords: Frozen sandstone Different thawing temperature Fracture toughness Microcrack propagation process Unfrozen water content

    ABSTRACT Due to the presence of ice and unfrozen water in pores of frozen rock,the rock fracture behaviors are susceptible to temperature.In this study,the potential thawing-induced softening effects on the fracture behaviors of frozen rock is evaluated by testing the tension fracture toughness (KIC) of frozen rock at different temperatures(i.e.-20 °C,-15 °C,-12 °C,-10 °C,-8 °C,-6 °C,-4 °C,-2 °C,and 0 °C).Acoustic emission (AE) and digital image correlation (DIC) methods are utilized to analyze the microcrack propagation during fracturing.The melting of pore ice is measured using nuclear magnetic resonance(NMR)method.The results indicate that:(1)The KIC of frozen rock decreases moderately between-20 °C and -4 °C,and rapidly between -4 °C and 0 °C.(2) At-20 °C to-4 °C,the fracturing process,deduced from the DIC results at the notch tip,exhibits three stages: elastic deformation,microcrack propagation and microcrack coalescence.However,at -4 °C-0 °C,only the latter two stages are observed.(3) At-4 °C-0 °C,the AE activities during fracturing are less than that at -20 °C to -4 °C,while more small events are reported.(4) The NMR results demonstrate a reverse variation trend in pore ice content with increasing temperature,that is,a moderate decrease is followed by a sharp decrease and-4 °C is exactly the critical temperature.Next,we interpret the thawing-induced softening effect by linking the evolution in microscopic structure of frozen rock with its macroscopic fracture behaviors as follow:from-20 °C to-4 °C,the thickening of the unfrozen water film diminishes the cementation strength between ice and rock skeleton,leading to the decrease in fracture parameters.From-4 °C to 0 °C,the cementation effect of ice almost vanishes,and the filling effect of pore ice is reduced significantly,which facilitates microcrack propagation and thus the easier fracture of frozen rocks.

    1.Introduction

    Increasing global warming now leads to numerous large-scale,high-potential rock slope failures,rockfalls,and other geological hazards in alpine regions(Gruber and Haeberli,2007;Gobiet et al.,2014),and it brings serious threats to infrastructure,human life,and property safety(Phillips,2006;Huggel et al.,2008).Numerous evidences indicate that melting of the permafrost in high mountains is a major factor for destabilization of rock slope(Huggel et al.,2010;Ravanel et al.,2010;Etzelmüller et al.,2022).Mountain permafrost is generally characterized by thick frozen bedrock,and the mechanical properties of which primarily determine the stability of rock slopes in high mountains (Harris et al.,2009;Phillips et al.,2017).Fracture toughness(KC),a common mechanical parameter of frozen rock,can be considered as a criterion of fracture initiation and propagation to estimate the instability of rock slopes (Fischer et al.,2010;Krautblatter et al.,2013).Due to the coexistence of unfrozen water and ice in pores of frozen rock,its fracture behaviors are susceptible to temperature change.Thus,investigating the fracture behaviors of frozen rock during thawing is crucial for engineering activities in mountain permafrost regions.

    In fracture mechanics,there are basically three modes of rock fracturing (Whittaker et al.,1992): tension (mode I),sliding shear(mode II),and tearing shear mode (mode III).As a result of the comparatively low tensile strength of rocks,tension fracture is the most common in failure analysis of rock slopes (Stead and Wolter,2015;Chen et al.,2021).In addition,tension fracture is encountered in other areas of rock engineering(He et al.,2010;Zhou et al.,2018),including hydraulic fracturing (Pan et al.,2021),shale gas extraction (Li et al.,2019),and deep nuclear waste (Peng et al.,2020).Therefore,numerous researches have been conducted on the tension fracture properties of rocks,including test methods of tension fracture toughnessKIC(Mu?oz-Ibá?ez et al.,2021;Kramarov et al.,2020),various influence factors (Cai et al.,2020;Deng et al.,2021)and fracturing process(Lin et al.,2020;Han et al.,2022).

    Previous researches have provided a basic understanding of the tension fracture process of rock.The International Society of Rock Mechanics and Rock Engineering (ISRM) proposed four standard ways to testKIC,namely,the short rod (SR),chevron bend (CB),cracked chevron notched Brazilian disc(CCNBD)and semi-circular bend(SCB)methods(Kuruppu et al.,2014).Studies on the influence factors ofKICshow that in addition to lithology and water content,temperature also has a significant influence.At high or low temperatures,the fracture toughness of rocks differs greatly from that at room temperature(Dwivedi et al.,2000;Feng et al.,2017;Miao et al.,2020).For the tension fracture,due to the singularity of cracks,a rather extensive plastic zone,termed as the fracture process zone(FPZ),is generated by microcracks around the cracks’tips.Variation of the FPZ directly influences the fracture propagation process(Pan et al.,2021).Methods such as digital image correlation(DIC) and acoustic emission (AE) (Wong and Guo,2019;Lin et al.,2020) have been frequently employed to determine the range of FPZ and investigate its evolution(Wei et al.,2015;Wu et al.,2011;Lin and Labuz,2013).

    While previous researches on frozen rock are mainly focused on the strength,deformability and failure features(Jia et al.,2020a)or the effects of freeze-thaw cycles on its mechanical properties(Jamshidi et al.,2017;Jia et al.,2020b;Zhu et al.,2021;Jamshidi,2021;Wang et al.,2022;Ren et al.,2022),little attention is paid to the fracture behaviors of frozen rock.Heating of frozen rock leads to melting of the pore ice and hence alters the pore structure of the rock,which will have a substantial effect on the fracture mechanical properties of frozen rock (Wang et al.,2022).Investigating the fracture behaviors of frozen rock becomes challenging.However,a proper evaluation of the effects of thawing on the fracture features of frozen rock is rarely reported.For this,a series of semi-circular bending (SCB) tests of frozen rocks at different thawing temperatures(T=-20°C,-15°C,-12°C,-10°C,-8°C,-6°C,-4°C,-2°C and 0°C) were conducted.AE and DIC techniques were then utilized to monitor the AE signals and surface strain of the sample during microcrack propagation,and the impact of thawing was also analyzed.The nuclear magnetic resonance(NMR)method was used to examine the change in unfrozen water content of rock.Combined with the change of unfrozen water content and microcrack propagation process,the mechanism of how thawing modifies the fracture behaviors of frozen rock was discussed.

    2.Methodology

    2.1.Test method of KIC

    The SCB samples are used in this study to evaluate theKICof frozen rock due to the easy preparation and testing procedure(Kuruppu et al.,2014).The size of the SCB sample is determined as shown in Fig.1a,and theKICis calculated according to Eq.(1)-(2)in Kuruppu et al.(2014).

    Fig.1.(a) Size of SCB sample;and (b) SCB rock samples.

    2.2.Digital image correlation (DIC) method

    DIC can be used to accurately measure optical surface deformation.The device consists of a light,a charge coupled device(CCD) camera,and a computing apparatus with relevant software.The correlation calculation approach is used to calculate the deformation information by analyzing the two images before and after deformation.The deformation analysis process of DIC is as follows:

    (1) A pre-deformation image is selected as the reference image and the region of interest is determined (see Fig.2);

    Fig.2.The DIC method to calculate deformation.

    (2) By specifying the subset size and distance,the region of interest is meshed to obtain uniformly spaced virtual grids;

    (3) A group of pixels on each point of the virtual grid is defined as a subset.It is assumed that each subset is in a rigid motion,ignoring the deformation of the subset itself.Then,correlation calculations are performed on each subset using a predefined correlation function and certain search strategies.In the deformed image,the region that has the strongest correlation with the reference subset is determined as the target subset,which marks the position of the reference subset after deformation,and the displacement between them is obtained accordingly(Fig.2) (Pan and Li,2011);and

    (4) The deformation information of the whole field can be obtained by computing all subsets.A zero-mean normalized sum of squared difference is utilized as the similarity measure of the reference and the target subsets in this study due to its robustness to scale and offset variations of illumination as well as its simplicity for optimization (Pan et al.,2010).

    2.3.NMR method

    If a rock sample is in the permanent magnet field of the NMR equipment,the hydrogen (H) protons of pore water generate magnetization vectors along the direction of the magnetic field.While a pulsed electromagnetic field is emitted,these H protons flip to a new equilibrium position,and then back to its original equilibrium position as the alternating electromagnetic field is withdrawn.This process is called relaxation.During the relaxation,the H protons will cut the magnetic flux line and produce an electrical signal,and the electrical signals continues to decay,which is called free induction decay (FID).The NMR signals of FID are positively correlated with the number of water molecules (Zhao et al.,2017).While,the NMR signals of H protons in ice can hardly be detected in low-intensity magnetic field.Therefore,the NMR signals of FID can be used to determine the liquid water (i.e.unfrozen water) content in frozen rock and soils (Tice et al.,1982;Watanabe and Wake,2009):

    The ice content can be calculated by

    whereSrwis the relative unfrozen water content(%),FIDfrozenis the NMR signals of FID at different thawing temperatures,FIDunfrozenis the FID value at 20°C,andSriis the relative ice content (%).

    2.4.Experimental design

    2.4.1.Sample preparation

    (1)Petrographic studies.The sample used in this study is yellow feldspar quartz sandstone bearing a massive structure without obvious bedding layers.In addition,it is composed of 90% quartz,5% feldspar sand and 5% Kaolinite.Fig.3 exhibits the single polarized and orthogonal polarized microscopic thin-section images of the sandstone.

    Fig.3.(a) The single polarized;and (b) orthogonal polarized microscopic thin-section images.

    (2)SCB sample preparation.Twenty two cylindrical samples(diameter×height:100 mm×45 mm)were drilled from the same sandstone block.They were oven-dried at 105°C and placed in a vacuum-saturated cylinder to pump for 4 h under-0.1 MPa.Water was then injected into the saturated cylinder to completely immerse the sample,the water-saturated mass of samples was measured after 72 h,and their porosity was determined.The Pwave test of each sample was also conducted.Among the 22 cylindrical samples,18 samples with similar porosity and P-wave velocity were selected and sliced into 36 semicircle SCB samples(Fig.1b).The dry density of samples was (2.1 ± 0.1) g/cm3,the porosity of samples was(19.18 ±0.2)%,and the P-wave velocity of samples was (2.37 ±0.6) km/s.

    Thirty six saturated SCB samples were wrapped with plastic films and placed in a moisturizing vat(at 20°C,relative humidity of 50%) for a week to ensure uniform water distribution inside the samples.The samples were then frozen for 48 h in an environmental chamber at -20°C.In order to conduct the DIC test,the sample’s surface should be treated before the test to show random speckle patterns.Therefore,the surface of specimens was first painted white;once the white paint dried,black dots were dotted with a marker to create randomly scattered speckle patterns.

    2.4.2.Tests of KIC

    TheKICtests at different temperatures were carried out via a hydraulic servo press,which is equipped with a high-precision environmental chamber (Fig.4a).Loading was controlled by the axial displacement at a rate of 0.12 mm/min,which satisfies the ISRM’s recommendation that the loading rate of the SCB test should not exceed 0.2 mm/min.The chamber has a temperature control accuracy of ±0.1°C.

    Fig.4.(a) The press machine which is fitted out with a high-precision environment chamber,and the AE and DIC strain capture systems;and (b) NMR test system.

    The lower temperatures limit in cold regions (except in polar regions)is usually above-20°C(Ma et al.,2022),thus it is rational to set -20°C as the starting point of thawing for frozen rock.TheKICof 36 SCB saturated samples at nine thawing temperatures(T=-20°C,-15°C,-12°C,-10°C,-8°C,-6°C,-4°C,-2°C,and 0°C)were tested.The testing steps are as follows:a frozen sample is taken out from the-20°C chamber,and four AE probes are fixed on the sample according to the positions given in Fig.1a,and it is then placed back in the chamber equipped to the press,as shown in Fig.4a.Afterwards,the temperature of the chamber was set as a preset test temperature and kept for 4 h,and theKICtest was conducted.

    The DS5-8B AE instrument (Beijing Softland Times Technology Co.,Ltd.)was employed to record AE signals during sample loading.The AE system can collect data at a rate of 3 Mhz (Fig.4a).The preamplifier and trigger threshold of AE acquisition system are both set as 40 dB.

    The variation of surface strain during loading was analyzed using the MatchID-2d/Stereo full field strain measurement system(Fig.4a).A photon FASTCAM Mini UX100 camera was vertically fixed at a tripod 0.5 m away from the press machine.The frame rate is 50 fps(frames per second),and the light source is LED cold light.

    2.4.3.Measurement of the unfrozen water content in frozen rock during thawing

    The MacroMR12-150H NMR instrument (Suzhou Niumai Analyzer Co.,Ltd.)is used in this study(Fig.4b).A low temperature control system is equipped with NMR receiving coil,and it can provide a uniform subzero temperature field.The parameters for NMR test is the same as that of Wang et al.(2021).The following are the steps for measuring the unfrozen water content of frozen rock during thawing:(1)A sample was firstly placed in the receiving coil and the temperature is set as-20°C;(2)After 4 h,the NMR test is conducted;(3)The temperatureTis then set to-15°C,-12°C,-10°C,-8°C,-6°C,-5°C,-4°C,-3°C,-2°C,-1°C,0°C and 20°C,successively,and (4) the NMR signal is measured at each temperature after 4 h.

    3.Results

    3.1.Changes in KIC and fracture energy during thawing

    TheKICrepresents the capacity of a rock to resist unstable propagation of cracks.The fracture energy refers to as the required energy per unit area of crack propagation,which is dissipated during crack propagation.The two parameters are used to describe the initiation and propagation of cracks for assessing the fracture characteristics of rocks.

    TheKICof frozen rock at various thawing temperatures is displayed in Fig.5a.It drops remarkably due to thawing and exhibits a two-stage mode of (a) slow decrease,KICdeclined by 19.5% asTdrops from-20°C to-4°C;and(b)rapid decrease,KICdeclined by 74%as-4°C ≤T<0°C.Therefore,-4°C is a critical temperature that marks the starting of rapid decrease inKICof frozen rock during thawing.

    Fig.5.(a) Changes of KIC and (b) fracture energy during thawing.

    During the SCB sample fracture,the fracture energy can be computed by the work done by the press system on the sample during loading,that is,the area is enclosed by the axial loaddisplacement curve before fracture (Peng et al.,2020).Based on this,the fracture energy of samples at various temperatures is displayed in Fig.5b.It shows that the fracture energy falls dramatically as the temperature increases and also exhibits a twostage mode of being essentially unchanged from-20°C to-12°C and rapidly reduces from-12°C to 0°C.

    3.2.Influence of thawing on fracture process of frozen rock

    Fig.6 illustrates the evolution of horizontal strain (HS) at the notch tip of SCB samples during fracture at various temperatures.It shows that the HS increases with the rise of axial stress,and thawing plays a dominant role.As -20°C ≥T>-4°C (Fig.6a-f),the HS vs.displacement curve exhibits a three-stage mode of (I)changeless stage;(II) moderate increase stage;and (III) rapid increase stage.At-4°C ≤T<0°C(Fig.6g-i),it exhibits a two-stage mode of(II)moderate increase stage,and(III)rapid increase stage.The mode I stage is not observed.

    Fig.6.(a-i)Variations in HS at notch tip during fracture at different temperatures of-20 °C;-15 °C;-12 °C;-10 °C;-8 °C;-6 °C;-4 °C;-2 °C;and 0 °C;and(j)HS of notch tip under peak load at different temperatures.

    (1) At -20°C ≥T>-4°C (Fig.6a-f),in the changeless stage(mode I),the HS is essentially unaffected by increase in the displacement and horizontal strain field (HSF) is basically evenly distributed(HSF at pointAon the HS vs.displacement curve).This suggests that the sample is undergoing a vertically elastic compression stage and the horizontal strain localization does not appear.In the moderate increase stage(mode II),the HS increases linearly as the displacement increases,and the HS around the notch tip is much higher than that in other regions(pointB),and this local region is the FPZ.The presence of FPZ implies the occurrence of strain localization due to expansion of the microcracks.In the rapid increase stage (mode III),the HS grows rapidly and the FPZ expands substantially(pointC).It demonstrates that the FPZ gradually develops,and as the microcrack coalesces,the precut macrocracks grow.Therefore,when the temperature is between-20°C and-4°C,the three stages of HS variation actually correspond to three stages in microcrack evolution at the notch tip of frozen rock: elastic deformation stage,microcrack propagation stage,microcrack coalescence stage.

    (2) At -4°C ≤T<0°C (Fig.6g-i),in mode II stage,the HS increases linearly and the FPZ(pointB)occurs,indicating that microcracks have propagated.At the mode III stage,the HS grows rapidly and the FPZ range expands significantly(pointC),which is the precursor of microcrack propagation.The absence of the changeless stage suggests that microcracks at the notch tip of frozen rock are more likely to propagate at higher temperatures.

    (3) Aiming to understand the effect of thawing on microcrack propagation in frozen rock,the HS at the notch tip of peak load (pointC) is measured at various temperatures,as depicted in Fig.6j.The HS at the notch tip grows approximately linearly with temperature increasing from -20°C to-4°C,while the increase in HS accelerates as temperature surpasses-4°C.This result indicates that thawing induces a significant softening effect on the fracture behaviors of frozen rock,and confirms again that -4°C is a critical temperature for the softening effect.

    According to the above analyses,pointAindicates initiation of the microcrack propagation,pointBdenotes the beginning of microcrack coalescence,and pointCsignifies the occurrence of macrocrack propagation.Therefore,loads at pointA,pointBand pointCcan be marked as the microcrack propagation load,microcrack coalescence load and macrocrack propagation load,respectively.The three loads are extracted and their variation trend with temperature is illustrated in Fig.7.It demonstrates that the thawing leads to a gradual decrease in microcrack propagation load,microcrack coalescence load,and macrocrack propagation load,and they all change in two stages.At-20°C to-4°C,the load at pointAis decreased to 0 kN,and the loads at pointBand pointCare moderately decreased.From-4°C to 0°C,the loads at pointsBandCfall rapidly,while the load at pointAremains constant at 0 kN.

    Fig.7.The variations in microcrack propagation load,microcrack coalescence load,and macrocrack propagation load at different thawing temperatures.

    As shown in Fig.6(pointC),the range of FPZ varies at different temperatures.The FPZ at the peak load depicts the result of microcrack propagation and coalescence prior to macro-fracture.Therefore,many scholars determine the range of FPZ by the variation of the HSF,or displacement field (Wu et al.,2011;Lin et al.,2019).The method of determining the FPZ length by displacement field is shown in Fig.8a and b.

    Fig.8.(a) Incremental horizontal displacement contours at peak load;(b)opening displacements along the notch of the sample at -8 °C;and (c) FPZ length at different thawing temperatures.

    Fig.8a shows the incremental horizontal displacement contours on the sample surface as computed by DIC.To reduce the calculation error,horizontal displacement values on two reference lines along the notch length 2 mm to the side of the notch are extracted(the lines of L1and L2in Fig.8a),and shown in Fig.8b.The distance between the notch tip and the confluence point of two horizontal displacement lines delineates the length of FPZ:Lp(Lin and Labuz,2013).Fig.8c depicts the variation of FPZ length at the peak load with temperatures.Due to thawing,the FPZ length increases and demonstrates a two-stage mode of (a) moderate increase and (b)rapid increase.Unsurprisingly,the increase in FPZ length begins to accelerate at -4°C.

    3.3.The influence of thawing on AE signals during fracturing of frozen rock

    3.3.1.Changes in AE counts and cumulative AE counts

    The changes in AE signals during fracturing of frozen rock under various testing temperatures are shown in Fig.9.At the beginning of loading,the AE counts are weak;with expansion of the microcrack,the AE counts rise.The cumulative AE counts grow in stages during microcrack propagation and exhibit a strong correlation with the change of HS.However,variation of the AE counts is sensitive to the temperature.

    At-20°C ≥T>-4°C,the AE counts gradually intensify during microcracks expansion.The cumulative AE counts rise in three stages (Fig.9a-f).In the elastic deformation stage,microcracks basically do not expand,and the cumulative AE counts rise slowly.In the microcrack propagation stage,increase in the cumulative AE counts is greater;while in the microcrack coalescence stage,the cumulative AE counts rise rapidly.However,as the testing temperature rises over-4°C,the AE counts are much lower than that at-20°C to-4°C(Fig.9g-h).Moreover,the cumulative AE counts grow in two stages: a slow increase during the microcrack propagation stage and a rapid increase during the microcrack coalescence stage.When the ice in the sample is completely melted (0°C),the AE counts during microcrack expansion are nearly zero,and the intensity of the AE signal is extremely low(Fig.9i).It demonstrates that at-4°C-0°C,the expansion of a large number of microcracks releases less elastic energy.

    3.3.2.Change of AE b value

    To quantify the magnitude of seismicity,the G-R formula was used to express the statistical correlation between earthquake magnitude and frequency (Gutenberg and Richter,1944):

    whereMis the magnitude,Nis the number of events for which the magnitude is higher thanM,andaandbare the constants.

    The parametersbandahave specific physical significance.Thebvalue represents the ratio of small-magnitude events to largemagnitude ones.The destruction of a rock sample is comparable to a natural earthquake,with the exception of differences in magnitude and radiation frequency.Consequently,the relationship between the amplitude of an earthquake and its frequency can be utilized to analyze AE events during rock failure.Eq.(3)is expressed as follows (Colombo et al.,2003;Goebel et al.,2013):

    whereAdBis the AE event amplitude,andN′is the hit number of AE events whose amplitude exceedsAdB.Thebvalue is used to analyze the variation in AE events and rock failure.This demonstrates that during rock failure,when microcrack propagation is dominant,higherb-values can be noticed;whereas macrocrack propagation occurs,bvalues fall (Colombo et al.,2003;Nejati and Ghazvinian,2013).

    Fig.10 shows variation of the AEbvalue during thawing.The AEbvalue exhibits two stages with temperatures:basically unchanged and rapidly increasing.At-20°C to-4°C,the AEbvalues remain essentially unchanged,indicating that thawing has slight impact on the amplitude and number of AE events during the loading process.At-4°C-0°C,AEbvalues increase rapidly,indicating that the large amplitude AE events decrease while the small amplitude events rise.In other words,the microcrack propagation is more significant at-4°C-0°C.This is consistent with the results displayed in Fig.9.

    Fig.10.AE b-value at different thawing temperatures.

    3.4.Variation of pore ice during thawing of frozen rock

    The fracture behaviors of frozen rock are susceptible to temperature change,which should be attributed to thawing of the pore ice.Based on the principles in Section 2.3,the variation of pore ice content during thawing can be determined by subtracting the unfrozen water content from overall water content.

    The variation of ice and unfrozen water content is shown in Fig.11.It shows that at-20°C to-4°C,the unfrozen water content gradually rises and the ice content decreases,indicating that a small amount of ice is melted.At-4°C-0°C,the unfrozen water content rapidly grows,and the ice content quickly decreases,showing that a lot of ice is melted.At 0°C,the ice content is zero,which suggests that the ice in the pores of frozen rock is completely melted.Coincidentally,the critical temperature for the variation of ice content is also -4°C,which may hold the interpretation of thawing-induced softening effects on the fracture behaviors of frozen rock.

    Fig.11.The relative unfrozen water content and ice content at varying thawing temperatures.

    4.Discussion

    Thawing can induce a significant softening effect on the fracture behaviors of frozen rock,which is evidenced by the decrease inKIC,fracture energy (Fig.5),and the increases in horizontal strain(Fig.6)and AE signals(Fig.9)with rising temperature from-20°C to 0°C.Moreover,-4°C is recognized as a critical temperature for the softening effect,where the softening effect intensifies (see Figs.5-8,10).Coincidentally,the NMR results demonstrate a reverse variation trend in pore ice content with increase of temperature,that is,a moderate decrease is followed by a sharp decrease and-4°C is exactly the critical point(Fig.11).Accordingly,we interpreted the thawing-induced softening effect by linking evolution in the microscopic structure of frozen rock with its macroscopic fracture behaviors.

    When water-saturated rocks were frozen at -20°C,a vast majority of pore water was froze but the rest remained unfrozen in the form of a thin film between ice and mineral grains(Jia et al.,2019;Wang et al.,2022).Pore ice and unfrozen water film has two strengthening effects on the mineral skeleton (Wang et al.,2021):

    (1) The filling effect of pore ice.A stress concentration effect will occur at a crack tip when the crack is under tension,leading to an increased stress.The stress concentration effect in the crack is positively correlated with its length-diameter ratio(Hertzberg,1996).The presence of ice drastically reduces the effective length of the crack (Fig.12,at -20°C),hence relieves the stress concentration effect at the crack tip.The filling effect depends on the pore ice content.

    Fig.12.Influence mechanisms of thawing on fracture of frozen rock.

    (2) The cementation effect of the unfrozen water film.The ice is cemented to the rock skeleton through the unfrozen water film (Fig.12,at -20°C),hence increasing the resistance to tension failure of the crack.The cementation effect is negatively related to the thickness of unfrozen water film,which depends on temperature (Guerin et al.,2016;Chen et al.,2018;Wang et al.,2021).The melting of ice within the cracks diminishes the two strengthening effects,resulting in a degradation of the fracture mechanical properties of frozen rock.

    The melting of crack ice initiated from its outer boundary,i.e.the unfrozen water film.When the temperature increases from-20°C to-4°C,a small quantity of ice is melted,resulting in thickening of the unfrozen water film (Fig.12,at -4°C).So,the cementation strength between ice and rock skeleton decreases,making microcracks easier to propagate (Figs.6a-g and 7).This process also is responsible for the increase of FPZ length(see Fig.8),and the crack is more likely to fracture.However,as the above two strengthening effects are still significant,high-amplitude AE events dominate during the microcrack propagation process,and the AE signals are relatively strong(Figs.9a-g and 10).The weakening of the fracture mechanical parameters of frozen rocks during thawing is primarily determined by the thickness of the unfrozen water film within this temperature interval.

    As temperature is greater than-4°C,much of the ice is melted(Fig.12,at 0°C),and the unfrozen water content rises rapidly(Fig.11b).As a result,the cementation effect of ice almost vanishes and the filling effect of pore ice is reduced significantly,facilitating the propagation of microcracks(Figs.6g-i and 7).In addition,due to the weakening of the two strengthening effects,small amplitude AE events dominate the microcrack expansion process,and the amplitude of AE signals decreases significantly (Figs.9g-i and 10).After the specimen is completely melted,both strengthening effects disappear,and the AE signals are very small (Fig.9i).

    5.Conclusions

    A series of SCB tests were conducted on frozen rocks at various thawing temperatures,and the microcrack propagation process was monitored by AE and DIC methods.By linking the evolution in microscopic structure of frozen rock with its macroscopic fracture behaviors,the influence mechanisms of thawing were interpreted.

    (1) As the thawing temperature rises,the variation inKICof frozen rock exhibits a two-stage mode of slow decrease stage,KICdeclined by 19.5%from-20°C to-4°C;and rapid decrease,KICdeclined by 74%from-4°C to 0°C.The fracture energy varies in two stages: (a) essentially unchanged(-20°C~-12°C);and(b)rapidly decreased(-12°C-0°C).

    (2) During the loading process,when the thawing temperature is at -20°C to-4°C,the HS at the notch tip varies in three stages:(I)changeless stage;(II)moderate increase stage;and(III) rapid increase stage,which correspond to the three stages: elastic deformation stage,microcrack propagation stage,and the microcrack coalescence stage.However,the change of HS at the notch tip shows a two-stage mode of(II)moderate increase stage and (III) rapid increase stage at -4°C-0°C,which corresponds to the microcrack propagation stage and microcrack coalescence stage.

    (3) The microcrack propagation load,microcrack coalescence load,and macrocrack propagation load all decrease in two stages.At-20°C to-4°C,the microcrack propagation load is reduced to 0 kN,while the microcrack coalescence load and macrocrack propagation load are moderately decreased.At-4°C-0°C,the microcrack coalescence load and macrocrack propagation load fall rapidly,while the microcrack propagation load remains constant at 0 kN.Moreover,when the temperature increases,the FPZ length grows significantly.

    (4) During fracturing of frozen rock,the AE counts change in three stages at -20°C to -4°C and in two stages at 0°C to-4°C,which is consistent with the variation of HS at the notch tip.The AE counts at 0°C to-4°C are much fewer than that at -20°C to -4°C.In addition,at -4°C-0°C,the microcrack propagation at the notch tip is more prominent than that at -20°C to-4°C.

    (5) The weakening mechanism of the fracture mechanical parameters of frozen rocks differs at various thawing stages.At-20°C to-4°C,the thickening of the unfrozen water film diminishes the cementation strength between ice and rock skeleton,so the weakening of fracture mechanical parameters is primarily determined by the thickness of the unfrozen water film.However,at-4°C-0°C,the cementation effect of ice almost vanishes,and the filling effect of pore ice is reduced significantly,which facilitates microcrack propagation.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgments

    We acknowledge the funding support from the National Natural Science Foundation of China (Grant No.42271148).

    九色亚洲精品在线播放| 不卡视频在线观看欧美| 亚洲欧美一区二区三区黑人 | 国产精品熟女久久久久浪| 999久久久国产精品视频| kizo精华| 久久免费观看电影| 欧美中文综合在线视频| 在线免费观看不下载黄p国产| tube8黄色片| 免费在线观看完整版高清| 在线精品无人区一区二区三| 纵有疾风起免费观看全集完整版| 黄片无遮挡物在线观看| 人人妻人人添人人爽欧美一区卜| 三级国产精品片| 日韩精品免费视频一区二区三区| 国产女主播在线喷水免费视频网站| 十分钟在线观看高清视频www| 啦啦啦视频在线资源免费观看| 另类精品久久| 欧美黄色片欧美黄色片| 日本av免费视频播放| 国产亚洲最大av| 2021少妇久久久久久久久久久| 国产深夜福利视频在线观看| 精品国产露脸久久av麻豆| 日韩精品有码人妻一区| av国产久精品久网站免费入址| 亚洲国产成人一精品久久久| 精品酒店卫生间| 亚洲熟女精品中文字幕| 啦啦啦中文免费视频观看日本| 激情视频va一区二区三区| 少妇被粗大的猛进出69影院| 精品福利永久在线观看| 男女边吃奶边做爰视频| 亚洲精品在线美女| 久久久久精品人妻al黑| 国产又爽黄色视频| 香蕉丝袜av| 国产成人91sexporn| av卡一久久| 狠狠精品人妻久久久久久综合| 精品少妇黑人巨大在线播放| 日产精品乱码卡一卡2卡三| 黑人猛操日本美女一级片| 久久久久精品性色| 免费看不卡的av| 男女啪啪激烈高潮av片| 18禁裸乳无遮挡动漫免费视频| 一本久久精品| 国产乱人偷精品视频| 亚洲成人av在线免费| 久久久精品国产亚洲av高清涩受| 午夜福利视频精品| 黑人巨大精品欧美一区二区蜜桃| 久久女婷五月综合色啪小说| av又黄又爽大尺度在线免费看| 天天躁日日躁夜夜躁夜夜| 日韩制服丝袜自拍偷拍| 一级毛片黄色毛片免费观看视频| 欧美日韩亚洲国产一区二区在线观看 | 99国产综合亚洲精品| 亚洲欧美成人综合另类久久久| 久久国产亚洲av麻豆专区| 18禁动态无遮挡网站| 精品酒店卫生间| 国产精品秋霞免费鲁丝片| 久久午夜福利片| 十八禁高潮呻吟视频| 亚洲成人av在线免费| 18禁裸乳无遮挡动漫免费视频| 亚洲四区av| 亚洲国产看品久久| 2022亚洲国产成人精品| 国产av国产精品国产| 久久久久久久国产电影| av.在线天堂| 国产男人的电影天堂91| 久久人人爽人人片av| 亚洲视频免费观看视频| 国产极品天堂在线| 看免费av毛片| 在线观看人妻少妇| 男女边摸边吃奶| 亚洲av电影在线观看一区二区三区| 一级黄片播放器| 亚洲欧美成人精品一区二区| 26uuu在线亚洲综合色| 久久国产亚洲av麻豆专区| 99久久人妻综合| 亚洲成人av在线免费| 亚洲av.av天堂| 欧美人与性动交α欧美软件| 18+在线观看网站| 最新的欧美精品一区二区| 中国国产av一级| 妹子高潮喷水视频| 亚洲中文av在线| 久久这里只有精品19| 午夜91福利影院| 老女人水多毛片| 久久人妻熟女aⅴ| av网站在线播放免费| 一级a爱视频在线免费观看| 宅男免费午夜| 黄片无遮挡物在线观看| 国语对白做爰xxxⅹ性视频网站| 精品亚洲成国产av| 日韩av在线免费看完整版不卡| 在线观看国产h片| 啦啦啦啦在线视频资源| 亚洲av电影在线进入| 亚洲av.av天堂| 国产精品嫩草影院av在线观看| 精品视频人人做人人爽| 国产熟女午夜一区二区三区| 我要看黄色一级片免费的| 一区福利在线观看| 最黄视频免费看| 大片免费播放器 马上看| 国产精品久久久av美女十八| av在线观看视频网站免费| 久久久久久久精品精品| 亚洲少妇的诱惑av| 亚洲精品视频女| 天堂8中文在线网| 午夜精品国产一区二区电影| 狠狠婷婷综合久久久久久88av| 欧美日韩亚洲国产一区二区在线观看 | 免费在线观看黄色视频的| 精品久久久精品久久久| 国产不卡av网站在线观看| 波多野结衣一区麻豆| 精品卡一卡二卡四卡免费| 国产一区有黄有色的免费视频| 美女大奶头黄色视频| av天堂久久9| 午夜精品国产一区二区电影| 欧美精品人与动牲交sv欧美| av在线播放精品| av在线app专区| 亚洲精品aⅴ在线观看| 制服人妻中文乱码| 高清在线视频一区二区三区| av网站免费在线观看视频| 青青草视频在线视频观看| 韩国精品一区二区三区| 亚洲av中文av极速乱| 人人妻人人添人人爽欧美一区卜| 国产色婷婷99| 丝袜人妻中文字幕| 久久精品国产综合久久久| 美女午夜性视频免费| freevideosex欧美| 午夜91福利影院| h视频一区二区三区| 黄片无遮挡物在线观看| 欧美国产精品va在线观看不卡| 1024香蕉在线观看| av免费在线看不卡| 亚洲精品中文字幕在线视频| 国产有黄有色有爽视频| 久久久久久久久免费视频了| 天天躁狠狠躁夜夜躁狠狠躁| 国产一区二区三区av在线| 一级毛片我不卡| 欧美变态另类bdsm刘玥| 亚洲伊人久久精品综合| 国产精品国产三级国产专区5o| 久久久久久久久久久免费av| 男女国产视频网站| 波多野结衣av一区二区av| 18在线观看网站| 毛片一级片免费看久久久久| 久久毛片免费看一区二区三区| 精品国产露脸久久av麻豆| 黑人猛操日本美女一级片| 亚洲国产欧美日韩在线播放| 亚洲av欧美aⅴ国产| 综合色丁香网| 久热这里只有精品99| 成人黄色视频免费在线看| 99九九在线精品视频| 精品国产一区二区三区四区第35| 99香蕉大伊视频| 免费日韩欧美在线观看| 99国产精品免费福利视频| 一级毛片 在线播放| 一级黄片播放器| 亚洲人成77777在线视频| av网站在线播放免费| 国产精品一二三区在线看| 免费高清在线观看日韩| 综合色丁香网| 午夜福利一区二区在线看| 少妇被粗大猛烈的视频| 少妇 在线观看| 久久久国产精品麻豆| 高清av免费在线| 国产精品久久久久久精品古装| 丝袜美腿诱惑在线| 精品少妇久久久久久888优播| 精品人妻一区二区三区麻豆| 国产熟女午夜一区二区三区| 精品少妇一区二区三区视频日本电影 | 国产成人午夜福利电影在线观看| 桃花免费在线播放| 麻豆精品久久久久久蜜桃| 国产精品av久久久久免费| 国产片内射在线| 女的被弄到高潮叫床怎么办| av在线观看视频网站免费| 午夜福利影视在线免费观看| 免费大片黄手机在线观看| 在线观看三级黄色| 十八禁高潮呻吟视频| 亚洲精品视频女| 咕卡用的链子| 国产成人一区二区在线| 免费不卡的大黄色大毛片视频在线观看| 日韩一区二区三区影片| 亚洲一码二码三码区别大吗| 精品少妇久久久久久888优播| 免费少妇av软件| av免费观看日本| a 毛片基地| www.av在线官网国产| 亚洲精品国产一区二区精华液| 国产一级毛片在线| 人成视频在线观看免费观看| 日韩欧美精品免费久久| 亚洲精品美女久久久久99蜜臀 | 国产97色在线日韩免费| 欧美成人精品欧美一级黄| 纵有疾风起免费观看全集完整版| 国产免费视频播放在线视频| 欧美黄色片欧美黄色片| 丰满乱子伦码专区| 十八禁高潮呻吟视频| 成人手机av| 欧美 亚洲 国产 日韩一| 亚洲精品美女久久av网站| 午夜久久久在线观看| 伊人久久大香线蕉亚洲五| 亚洲精品视频女| 亚洲精品国产av成人精品| 91精品三级在线观看| 各种免费的搞黄视频| 国产欧美日韩一区二区三区在线| 亚洲精品视频女| 午夜福利影视在线免费观看| 国产成人精品在线电影| 高清不卡的av网站| 亚洲精品美女久久久久99蜜臀 | 中文字幕色久视频| 久久人人97超碰香蕉20202| 亚洲精品久久久久久婷婷小说| 美国免费a级毛片| 亚洲精品美女久久久久99蜜臀 | 国产精品人妻久久久影院| 国产精品久久久久久久久免| 亚洲人成电影观看| 免费黄色在线免费观看| 久久久久久久久久久久大奶| 色婷婷久久久亚洲欧美| 国产精品av久久久久免费| 看免费av毛片| a级片在线免费高清观看视频| 中文字幕色久视频| 免费黄频网站在线观看国产| 久久久精品国产亚洲av高清涩受| 久久久久久久国产电影| 日韩中字成人| 啦啦啦啦在线视频资源| 亚洲人成电影观看| 秋霞在线观看毛片| 久久久久久久久久久久大奶| 亚洲国产色片| 亚洲一码二码三码区别大吗| 国语对白做爰xxxⅹ性视频网站| 日本黄色日本黄色录像| 国产精品女同一区二区软件| 久久久欧美国产精品| 亚洲成国产人片在线观看| 18禁观看日本| 日本色播在线视频| 亚洲精品成人av观看孕妇| 日本免费在线观看一区| 日韩欧美精品免费久久| 国产一区二区激情短视频 | 国产av精品麻豆| 亚洲,一卡二卡三卡| 1024香蕉在线观看| 国产人伦9x9x在线观看 | 天美传媒精品一区二区| 久久亚洲国产成人精品v| 啦啦啦在线观看免费高清www| 黑人巨大精品欧美一区二区蜜桃| 亚洲av成人精品一二三区| 超碰97精品在线观看| 亚洲av国产av综合av卡| 九色亚洲精品在线播放| 天天操日日干夜夜撸| 99九九在线精品视频| 在线天堂中文资源库| 黑人欧美特级aaaaaa片| 宅男免费午夜| 精品第一国产精品| 亚洲精华国产精华液的使用体验| av一本久久久久| 天天影视国产精品| 国产成人精品婷婷| 免费少妇av软件| 最近最新中文字幕免费大全7| 在线亚洲精品国产二区图片欧美| av福利片在线| 如何舔出高潮| 国产一区二区三区av在线| 久久午夜福利片| 日韩在线高清观看一区二区三区| 国产高清国产精品国产三级| 97在线人人人人妻| 欧美日韩成人在线一区二区| 久久久精品区二区三区| 午夜av观看不卡| 亚洲五月色婷婷综合| 久久99一区二区三区| 久久国内精品自在自线图片| 国产成人精品久久久久久| 午夜91福利影院| 色94色欧美一区二区| av在线播放精品| 91精品三级在线观看| 黑人欧美特级aaaaaa片| 精品少妇黑人巨大在线播放| 亚洲 欧美一区二区三区| 国产毛片在线视频| 久久99热这里只频精品6学生| 亚洲视频免费观看视频| 9色porny在线观看| 久久精品国产亚洲av涩爱| 中文字幕最新亚洲高清| 国产在线一区二区三区精| 久久久久精品久久久久真实原创| 熟妇人妻不卡中文字幕| 亚洲欧洲国产日韩| 女人高潮潮喷娇喘18禁视频| 一本—道久久a久久精品蜜桃钙片| 久久精品国产亚洲av天美| 欧美人与性动交α欧美精品济南到 | 亚洲精品一二三| 亚洲精品国产色婷婷电影| 秋霞伦理黄片| av一本久久久久| 赤兔流量卡办理| 久久久久国产一级毛片高清牌| tube8黄色片| 免费不卡的大黄色大毛片视频在线观看| 国产日韩欧美亚洲二区| 一本色道久久久久久精品综合| 国产野战对白在线观看| 少妇的丰满在线观看| 精品亚洲成a人片在线观看| 极品人妻少妇av视频| 国产在线一区二区三区精| 日韩制服骚丝袜av| 天天躁狠狠躁夜夜躁狠狠躁| 高清av免费在线| 一区福利在线观看| 欧美另类一区| 宅男免费午夜| 最近中文字幕高清免费大全6| 国产又色又爽无遮挡免| 久久午夜综合久久蜜桃| 国产毛片在线视频| 中文字幕av电影在线播放| 精品酒店卫生间| 国产成人a∨麻豆精品| 亚洲av电影在线进入| 午夜91福利影院| 国产熟女欧美一区二区| 在现免费观看毛片| 欧美精品人与动牲交sv欧美| 国产乱来视频区| 18在线观看网站| 丰满饥渴人妻一区二区三| 香蕉国产在线看| 亚洲成人一二三区av| 日本午夜av视频| 18禁动态无遮挡网站| 视频区图区小说| 亚洲图色成人| 最近最新中文字幕免费大全7| 一级爰片在线观看| 日本猛色少妇xxxxx猛交久久| 两个人免费观看高清视频| 国产精品 欧美亚洲| 国产成人午夜福利电影在线观看| 人人妻人人添人人爽欧美一区卜| 国产免费现黄频在线看| 2021少妇久久久久久久久久久| 国产成人精品无人区| 在线 av 中文字幕| 2022亚洲国产成人精品| 欧美 日韩 精品 国产| 人人妻人人澡人人看| 亚洲av在线观看美女高潮| av免费在线看不卡| 精品一品国产午夜福利视频| 亚洲精品国产色婷婷电影| 亚洲,欧美精品.| 久久人人爽av亚洲精品天堂| 女性生殖器流出的白浆| 人妻人人澡人人爽人人| 日本-黄色视频高清免费观看| 国产精品熟女久久久久浪| av网站免费在线观看视频| 99热全是精品| 国产熟女午夜一区二区三区| 亚洲少妇的诱惑av| 日本猛色少妇xxxxx猛交久久| 国产又色又爽无遮挡免| av国产精品久久久久影院| 黄色一级大片看看| 我要看黄色一级片免费的| 欧美精品一区二区大全| 啦啦啦啦在线视频资源| 亚洲国产av新网站| 亚洲第一av免费看| av一本久久久久| 日韩欧美一区视频在线观看| 精品第一国产精品| 国产在线免费精品| 嫩草影院入口| 高清黄色对白视频在线免费看| 97精品久久久久久久久久精品| 成人手机av| 狠狠精品人妻久久久久久综合| 亚洲综合色惰| 久热这里只有精品99| 青青草视频在线视频观看| 成人免费观看视频高清| 老司机影院成人| a级毛片在线看网站| 热re99久久精品国产66热6| 五月伊人婷婷丁香| 国产日韩欧美亚洲二区| 欧美少妇被猛烈插入视频| 国产伦理片在线播放av一区| 热re99久久国产66热| 丝袜美腿诱惑在线| 国产成人精品福利久久| 色播在线永久视频| 亚洲av日韩在线播放| 亚洲欧美精品综合一区二区三区 | 777久久人妻少妇嫩草av网站| 免费人妻精品一区二区三区视频| 久久人人97超碰香蕉20202| 亚洲成国产人片在线观看| 日韩欧美一区视频在线观看| 91在线精品国自产拍蜜月| 午夜福利在线免费观看网站| 亚洲一区中文字幕在线| 国产精品一国产av| 日日啪夜夜爽| 一级毛片电影观看| av有码第一页| 亚洲一区中文字幕在线| 香蕉国产在线看| 国产乱人偷精品视频| 日本猛色少妇xxxxx猛交久久| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 人妻 亚洲 视频| 黄色怎么调成土黄色| 人人妻人人澡人人爽人人夜夜| 亚洲四区av| 欧美亚洲 丝袜 人妻 在线| 久久久久久人妻| 亚洲精品乱久久久久久| 免费看不卡的av| 少妇猛男粗大的猛烈进出视频| 精品国产露脸久久av麻豆| 最近2019中文字幕mv第一页| 美女主播在线视频| 波多野结衣av一区二区av| 国产亚洲最大av| 国产精品.久久久| 久久久久网色| 哪个播放器可以免费观看大片| 国产日韩欧美亚洲二区| 国产人伦9x9x在线观看 | 久久久久精品性色| 69精品国产乱码久久久| tube8黄色片| 成人漫画全彩无遮挡| 伦理电影大哥的女人| 在线 av 中文字幕| 国产熟女午夜一区二区三区| 午夜影院在线不卡| 天美传媒精品一区二区| 亚洲欧美一区二区三区黑人 | 一级毛片 在线播放| 国产精品国产av在线观看| 国产av国产精品国产| 寂寞人妻少妇视频99o| 麻豆乱淫一区二区| av片东京热男人的天堂| 少妇的逼水好多| 中文字幕av电影在线播放| 满18在线观看网站| 国产精品久久久久久精品古装| 日韩大片免费观看网站| 国产成人一区二区在线| 久久国内精品自在自线图片| 如日韩欧美国产精品一区二区三区| 国产野战对白在线观看| 久久热在线av| 99国产综合亚洲精品| 99热全是精品| 国产97色在线日韩免费| 久久av网站| 国产免费一区二区三区四区乱码| 国产黄色免费在线视频| 老熟女久久久| 成人影院久久| 精品国产乱码久久久久久男人| 黄网站色视频无遮挡免费观看| 日韩不卡一区二区三区视频在线| 日韩一区二区视频免费看| 亚洲,一卡二卡三卡| 纵有疾风起免费观看全集完整版| 久久久久久人人人人人| 国产成人精品婷婷| 久久久久国产网址| 中文乱码字字幕精品一区二区三区| 国产熟女午夜一区二区三区| 亚洲三级黄色毛片| 国产成人免费无遮挡视频| 老鸭窝网址在线观看| 如何舔出高潮| 国产精品三级大全| 午夜福利在线免费观看网站| 国产日韩欧美视频二区| 最近最新中文字幕免费大全7| 午夜福利,免费看| 欧美人与性动交α欧美软件| 亚洲天堂av无毛| 久久久久久久亚洲中文字幕| 久久人人爽av亚洲精品天堂| 久久精品国产亚洲av高清一级| 欧美 日韩 精品 国产| 久久久国产一区二区| 伦理电影大哥的女人| 久久久久国产精品人妻一区二区| 日韩制服骚丝袜av| 尾随美女入室| 国产精品嫩草影院av在线观看| 中文字幕人妻熟女乱码| 两个人免费观看高清视频| 亚洲国产欧美网| 免费人妻精品一区二区三区视频| 精品国产露脸久久av麻豆| 韩国高清视频一区二区三区| 爱豆传媒免费全集在线观看| 日韩,欧美,国产一区二区三区| 下体分泌物呈黄色| 精品久久久久久电影网| 欧美黄色片欧美黄色片| 女性被躁到高潮视频| 又大又黄又爽视频免费| 久久毛片免费看一区二区三区| 美国免费a级毛片| 欧美日韩视频精品一区| 一区福利在线观看| av国产久精品久网站免费入址| 91国产中文字幕| 精品一区二区三卡| 久久国产亚洲av麻豆专区| 在线看a的网站| 中文字幕亚洲精品专区| 男女边吃奶边做爰视频| 人体艺术视频欧美日本| 久热久热在线精品观看| 免费av中文字幕在线| 久久国产精品大桥未久av| 亚洲精品在线美女| 免费久久久久久久精品成人欧美视频| 亚洲经典国产精华液单| 一边摸一边做爽爽视频免费| 久久久久国产一级毛片高清牌| 女人久久www免费人成看片| 亚洲欧美精品综合一区二区三区 | 999久久久国产精品视频| 欧美精品av麻豆av| 狠狠婷婷综合久久久久久88av| 亚洲国产成人一精品久久久| 男女国产视频网站| 免费女性裸体啪啪无遮挡网站| 纯流量卡能插随身wifi吗| 在线观看三级黄色| 久久久久久久久久久免费av| 一区二区三区激情视频| 日韩av免费高清视频| 波野结衣二区三区在线| 婷婷色av中文字幕| 欧美另类一区| 午夜福利在线观看免费完整高清在| 欧美成人精品欧美一级黄| 国产视频首页在线观看| 寂寞人妻少妇视频99o| 久久精品国产亚洲av天美| 观看美女的网站| 麻豆av在线久日|