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

    The slip activity during the transition from elastic to plastic tensile deformation of the Mg-Al-Mn sheet

    2021-10-28 10:09:48DobroDrozdenkoFeketeKnpekBohlenChmel
    Journal of Magnesium and Alloys 2021年3期

    P.Dobroˇn,D.Drozdenko,K.Fekete,M.Knpek,J.Bohlen,F.Chmelík

    a Charles University,Faculty of Mathematics and Physics,Department of Physics of Materials,Ke Karlovu 5,12116,Prague 2,Czech Republic

    b The Czech Academy of Sciences,Nuclear Physics Institute,Hlavní130,250 68,ˇRe?,Czech Republic

    c Helmholtz-Zentrum Geesthacht,MagIC-Magnesium Innovation Centre,Max-Planck-Stra?e 1,D21502 Geesthacht,Germany

    Abstract The deformation behavior of the Mg-Al-Mn sheet was investigated during tensile loading along the rolling(RD)and transversal direction(TD)with special attention to the early stage of deformation.The activity of dislocation slip systems during the transition from elastic to plastic deformation was revealed by the acoustic emission(AE)technique.The parametrization and statistical AE analysis using the adaptive sequential k-mean(ASK)clustering provided necessary information about the individual deformation mechanisms and their evolution.The AE findings were supported by microstructural analyses,including in-situ secondary electron(SE)imaging and Schmid factor estimation for the activity of particular dislocation slip systems with respect to the loading direction.It was found that basalslip is the dominating mechanism up to the stress of~80 MPa in both loading directions with an absolute dominance during the RD-loading,while during the TD-loading,the contribution of prismaticslip to the deformation at stresses above 50 MPa was determined.Below the yielding in both loading directions,the predominance of prismaticover pyramidalslip was found at the stress in the range of 80–110 MPa and the opposite tendency occurred at stresses between 110 and 140 MPa.? 2021 Chongqing University.Publishing services provided by Elsevier B.V.on behalf of KeAi Communications Co.Ltd.This is an open access article under the CC BY-NC-ND license(http://creativecommons.org/licenses/by-nc-nd/4.0/)Peer review under responsibility of Chongqing University

    Keywords:Magnesium sheet;Yield asymmetry;Acoustic emission;Clustering.

    1.Introduction

    In general,Mg sheets exhibit anisotropic mechanical properties and low formability,which limit their use in technical applications.These negative effects can be effectively reduced by using appropriate alloying elements affecting the grain size and the initial texture.The strengthening effect of Al and the grain refinement effect of Mn in Mg alloys have been reported,e.g.,in[1–6].Therefore,Mg-Al-Mn sheets with good strength and formability even at room temperature(RT)are promising candidates for their application,for example,in the automotive industry.

    The macroscopic deformation behavior of Mg sheets has been extensively studied and several works have focused on revealing the activity of individual deformation mechanisms.Particularly,the activity ofanddislocations has been experimentally documented using electron microscopy(e.g.,[7]),diffraction methods[8,9],and theoretically estimated using molecular dynamics simulations[10]or by applying dislocation and elastoplastic polycrystal models to deformation and texture data[11–13].

    In-situobservations of dislocation activity using scanning and transmission electron microscopy(SEM,TEM)are becoming more popular nowadays[14–17],despite providing information from a limited volume of samples.On the other hand,diffraction methods operate within higher penetration depths.Both the electron microscopy and diffraction methods,however,require a certain time for data collection and,therefore,can be rather classified as semiin-situmethods.In contrast,the acoustic emission(AE)monitoring during deformation provides information on local dynamic changes from the entire volume of the material in real-time.The development of the AE technique in the last decade enables the storage and processing of data with a very high time resolution.Moreover,the application of modern statistical analyses to the AE data,particularly clustering techniques,significantly contributes to the identification and understanding of the processes that occur during deformation of the material[18].The combination of the AE technique with other methods provides complementary information about deformation mechanisms.For example,the AE technique is sensitive to collective dislocation movement and twin nucleation[19,20],and the neutron diffraction(ND)method refers to dislocation density and changes in twin volume fraction(twin growth).Using those techniques,it was shown that during tension of cast(texturefree)pure Mg[21],the extension twins were nucleated during the entire test(thin twins),while during compression,twins were formed preferentially around the yield point(YP)and then only their massive growth/thickening was observed.Thus,to fulfill the von-Mises criterion for the compatibility of plastic deformation[22],a significant twinning activity occurs during both tensile and compressive loading of cast Mg.In wrought Mg alloys,the basal texture strongly affects the activity of individual deformation mechanisms with respect to the loading direction(and deformation temperature),which is also reflected in the AE response[12,23–26].Particularly,during the compressive loading of wrought Mg alloys along the extrusion/rolling direction,the extension twins are preferably activated around the YP,providing the necessary compatibility of plastic deformation along the c-axis,as described above for cast Mg alloys.On the contrary,during the tensile loading of wrought Mg alloys,dislocation glide is the dominating deformation mechanism and the contribution of basal and non-basal slip systems to the deformation is still intensively discussed.Previous AE results obtained during the dislocation slip-dominated deformation in wrought Mg alloys,e.g.,[23–25,27],indicate that the AE technique can be advantageously used to study the collective dislocation movement,particularly at the early stage of deformation.

    According to generally accepted dislocation-based AE models,a release of the AE energy during deformation originates from(i)relaxation of the elastic stress field in the lattice caused by the passage of dislocations[28],(ii)annihilation of dislocations[29],and/or(iii)the radiation of wave energy from accelerating or decelerating dislocations[30–32].Based on many dislocation models,it has been assumed by Vinogradov et al.[33]that the AE energy is proportional to the density of mobile dislocationsρmand the squared mean velocity of dislocations〈v〉2.

    To track the changes in the dynamics of mobile dislocations in the material(also concerning variations in the mean free path),the AE count rate can be effectively used.This AE parameter is defined as a number of times the AE signal exceeds a pre-set threshold per time unit[34].It follows from the above statements that the AE energy and AE count rate provide different information about dynamic processes in the material.The AE energy depends on the mobile dislocation density and the AE count rate is proportional to the rate of change in the density of mobile dislocations.The mobility of dislocations is very closely related to the mean free path〈L〉of dislocations and depends on the dislocation densityρ.

    The adaptive sequential k-mean(ASK)clustering analysis,developed by Pomponi and Vinogradov[35],operates with the power spectral density(PSD)function G(f)which is calculated by the Welch’s method.The basic AE parameters,such as the AE energyEand the median frequencyfmof the PSD function,are defined in the working frequency rangefminandfmaxas follows:

    This statistical method has been successfully applied to AE data received during the deformation of cast or wrought Mg alloys.Clusters created by this algorithm exhibit specific features and,taking into account the microscopic and diffraction observations,can be assigned to individual deformation mechanisms[36–38].

    The aim of the present study is to obtain comprehensive information on the deformation behavior of the Mg-Al-Mn sheet in an earlier stage of deformation;particularly,on the transition between elastic to plastic deformation with respect to the loading direction.The activity of individual slip systems at the early deformation stage crucially influences the yielding behavior of the Mg sheet and the AE monitoring with the statistical analysis of the AE data enables us to track this activity in real-time at a very high time resolution.

    2.Experimental procedure

    The AM61(Mg+5.9 wt.% Zn+0.6 wt.% Mn)alloy was gravity cast at 730 °C into a rectangular steel crucible from which a slab with a thickness of 20 mm was machined.Before rolling,homogenization annealing at 430 °C for 20 h was carried out on the slab.After rolling at 430 °C(14 passes with intermediate re-heating for 10 min between passes),a final sheet with a thickness of 1.8 mm was heat-treated at 430 °C for 1 h.More information on the sheet production is given in[2].

    Dog bone-shaped tensile samples with gauge sizes of 10 mm×1.9 mm×30 mm were machined from the sheet with the loading direction along(i)the rolling direction(RD)and(ii)the transversal direction(TD).Deformation tests were performed using the Instron 5882(Norwood,MA,USA)and the Zwick Z50(ZwickRoell,Ulm-Einsingen,Germany)universal testing machine at RT and a constant crosshead speed giving an initial strain rate of 10?3s?1.The possible influence of the equipment setting and environment(background noise)on the AE response was also elaborated by measuring AE during tensile tests performed on two testing machines located in different laboratories,showing practically no differences in the AE results.

    AE was monitored during the tensile tests using the Physical Acoustics Corporation(PAC)AE system(MISTRAS Group,Princeton Jct,NJ,USA)consisting of a PC-integrated PCI-2 AE measurement board,a voltage preamplifier(type 2/4/6),and a broadband AE sensor(MICRO-30S).The AE sensor was attached directly to the sample with a wooden clip and a good acoustic contact between the sample and the sensor was achieved using the vacuum grease Apiezon M(Cryoandmore,Neuss,Germany).The preamplifier was set to 40 dB to ensure a high signal-to-noise ratio,and the AE data were recorded and stored with a sampling frequency of 1 MHz.Threshold-level signal detection(26 dB)was used to determine the AE count rate.

    The ASK analysis was performed on the continuously recorded AE data during the tensile tests.The AE data from the RD-and TD-loading were merged into a single file(RDTD or TD-RD)to ensure the same cluster identification for both loading directions.The ASK algorithm was applied to both merged files to eliminate any influence of the arrangement on the AE results.By comparing the ASK results obtained from both merged files,no differences were observed in the identified clusters.

    The AE data were firstly divided into successive frames containing 8192 samples representing a time interval of 8.192 ms,and subsequently,the PSD function was calculated for each frame.Thereafter,the PSD function was normalized by dividing this function by the AE energy in the working frequency range(Eq.3).The features of the PSD function in the first frame define the initial cluster.As a reference point for the analysis,the initial cluster was associated with background noise as the AE data acquisition began before the start of particular tensile tests.The successive frames were analyzed one by one and conditions for the formation of a new cluster were based on the k-means algorithm.More details on the algorithm can be found in[35].After the clustering procedure,each cluster is assigned to a dominant AE source mechanism.As mentioned above,the first cluster represents the background noise.In general,the assignment of other clusters to individual deformation mechanisms operating in Mg alloys is based on the characteristic features of the PSD function(energy,frequency distribution,etc.)for each deformation mechanism,and is supported by other experimental techniques(texture measurement,microscopy,etc.)or theoretical modeling.

    The global characterization of the texture of the Mg sheet was performed using X-ray diffraction(XRD).The Malvern Panalytical X-ray diffractometer X’Pert(Malvern,United Kingdom)setup using CuKαradiation was employed to measure six pole figures on the polished sample in a reflection geometry to a sample tilt of 70°.

    Fig.1.Initial microstructure and texture of the conventionally rolled and heat-treated AM61 sheet.

    For the microstructure investigations,conventional metallographic procedures including grinding with SiC paper and polishing with diamond pastes down to 1μm particle size were used.Finally,the surface was polished by means of the LectroPol-5 electrolytic polisher(Struers Inc.,Cleveland,OH,USA)with Struers-AC2 electrolyte at-40 °C.The electron backscatter diffraction(EBSD)measurement of the sample in the initial condition was performed in the RD-TD plane using the Zeiss Auriga Compact SEM(Jena,Germany)equipped with the FIB system and EDAX EDX/EBSD detector.The EBSD measurements were carried out at a working distance of 9 mm with 10 kV acceleration voltage.To reveal the activity of the deformation mechanisms,in-situsecondary electron(SE)imaging of the microstructure was performed at defined stresses during tensile loading along RD and TD using the MTEST Quattro(Materials Testing system ADMET)deformation stage installed inside a chamber of the Zeiss Auriga Compact SEM.The stress level selection(75,125,and 175 MPa)was based on the previously recorded AE response.

    3.Experimental results

    Fig.1 shows the microstructure(SEM)and texture(XRD)of the AM61 sheet after conventional rolling and heat treatment at 430 °C for 1 h.The investigated sheet exhibits a relatively homogeneous microstructure with an average grain size of 13±2μm and a strong basal texture with a slightly broader spread of basal planes towards RD than TD.The precipitates are distributed in the matrix and,occasionally,some particle strings were found in the microstructure.

    Table 1Tensile mechanical properties of the rolled AM61 sheet.TYS–tensile yield strength,UTS–ultimate tensile strength.

    The EBSD maps and Schmid factor(SF)analysis for theanddislocation slip systems from the perspective of RD-and TD-loading of the alloy in the initial state are presented in Fig.2.The majority of the grains in the EBSD map are oriented with the basal<0001>planes nearly parallel to the RD-TD plane of the sheet(i.e.with the c-axis perpendicular to both,RD and TD)and there are only a few grains with〈100〉or〈110〉planes oriented parallel to RD and/or TD.The SF maps were evaluated in order to assess the activation of the typical slip systems:(0001)〈110〉–basal,{100}〈120〉–prismatic,{101}〈110〉–pyramidal(first-order)and{112}〈11〉–pyramidal(second-order).From the SF maps and distributions for the basalslip,it is clearly seen that this slip system is more favorable(indicated by a higher fraction of area with a SF in the range of 0.4-0.5)for tensile loading along RD than TD.The SF distribution for prismatic and pyramidal slip systems is plotted in the range between 0.3 and 0.5 to highlight the difference between RD-and TD-loading behavior.In contrast to the basalglide(preferred rather for RD loading),the activation of the prismaticglide is expected to be more significant during the tensile loading along TD compared to the RD-loading(higher fraction of area with a SF of 0.4-0.5).The SF maps for the prismaticslip system indicate a relatively weak probability for its activation(i.e.,SF values are relatively low)and do not show a significant difference between the RD-and TD-loading(not shown here).The SF distribution of pyramidalglide does not show any orientation preference and it is rather comparable for both loading directions.Consequently,among all possibleslip systems,the highest activity of basalglide and the lowest one of pyramidalglide are predicted for the early stage of deformation.However,differences in the critical resolved shear stress(CRSS)for the activation of an individual slip system need to be taken into account.Furthermore,a slightly higher activity of the pyramidalslip system is expected during loading along TD compared to RD.

    Tensile deformation curves and the corresponding AE activity,represented by the AE count rate,AE energy,and median frequency of the AE signal,are presented in Fig.3.The deformation curves have a similar shape for both RD and TD-loading directions,with a higher yield strength for the TD-loading.A similar elongation of about 25%(with the initial strain rate of 10?3s?1it corresponds to 250 s)was achieved in both loading directions(Table 1).Monitoring of the AE activity during the deformation test revealed that a higher AE count rate and AE energy appear in the sample subjected to the RD-loading.Both AE parameters exhibit a typical evolution with their maximum at yielding,which is followed by a relatively distinct decrease.However,a change in the AE count rate is observed already before yielding and this effect is especially pronounced in the TD-loading.The median frequency of the AE signal is starting with its decrease at the beginning of the deformation,followed by an increase,which is strongly pronounced for the TD-loading,and further progressing with a gradual increase at later stages of deformation.In the case of the RD-loading,an additional drop around the time of 200 s can be seen on the median frequency–time curve,Fig.3a.To increase the visibility of the AE parameters at the beginning of the test(i.e.,the transition between elastic and plastic region),they were plotted against the tensile stress(Fig.4).The AE count rate increases rapidly from the beginning of the deformation test and then,in the stress range between 20 and 75 MPa,it increases either slightly(RD-loading,Fig.4a)or remains almost constant(TD-loading,Fig.4b).Thereafter,in the case of the RD-loading,a massive increase in the AE count rate can be observed immediately,while during the TD-loading,the AE count rate first increases slightly up to the tensile stress of~115 MPa and then increases rapidly up to the YP.After reaching its maxima at the YP,the AE count rate gradually decreases.During tensile loading up to the stress of~115 MPa,both loading directions exhibit similarities in the median frequency development.The median frequency is almost constant up to the stress level of~75 MPa and then drops sharply with a more pronounced effect in the case of the TD-loading.The following rise in the median frequency development is,therefore,only small for the RD-loading and more pronounced in the case of the TD-loading,Fig.4b.

    To reveal a dislocation glide and possible twinning activity during the transition from elastic to plastic tensile deformation,deformation tests performed inside the SEM chamber were interrupted at 75,125 and 175 MPa(black circles on the deformation curves,Fig.4),and the corresponding microstructures were investigated together with the initial state using SE imaging(Figs.5 and 6).The stress levels were chosen with respect to changes in the AE response in order to correlate these changes to the development of the microstructure.For the RD-loading,the first basal slip traces are observed in the microstructure subjected to the stress level of 125 MPa(Fig.5).In the case of the TD-loading,the first evidence of plastic deformation was found on the surface of the sample deformed up to 175 MPa(Fig.6).Moreover,the twinning activity was not observed in any of the loading directions in the stress range between 0 and 175 MPa.

    The ASK analysis was used to distinguish the changes in the activity of individual deformation mechanisms and to contribute to the understanding of the atypical AE response in the elastoplastic deformation region.For easier comparison of the cluster activity revealed by the ASK analysis and the above presented AE results(Figs.3–4),the dependence of the AE energy on the median frequency is presented for individual clusters in Fig.7.The background noise was measured before the tensile test and its energy is significantly lower than the energy released during the test.By comparing the results for the RD-and TD-loading,it is clearly seen that the released AE energy is an order of magnitude higher for the RD-loading and features a much narrower median frequency distribution in this loading direction.

    Fig.2.Initial microstructure of the investigated alloy:orientation map and SF analysis for the activity of basal and non-basalandslip systems for the RD-and TD-loading.Note that different format and scale(0-0.5 for basal,0.3-0.5 for prism and pyramidal,respectively)in SF plots are used for better visualization of the data.

    Fig.3.Tensile stress-time curves and AE parameters for(a)RD-and(b)TD-loading.

    Fig.4.Tensile stress-time curves and AE parameters plotted in the time interval of 0-120 s for(a)RD-and(b)TD-loading.The stress level of 75,125 and 175 MPa used for the SEM analysis is marked by circles on the deformation curves.Note that for the sake of clarity the scale of the AE energy in Fig.3b is 10×higher than that in Fig.4b.

    Fig.5.In-situ SE imaging of the microstructure during tensile loading along RD.Black arrows indicate the first visible traces of basalslip at image taken at stress of 125 MPa.

    Fig.6.In-situ SE imaging of the microstructure during tensile loading along TD.Black arrows indicate the first visible slip traces at image taken at stress of 175 MPa.

    Fig.7.Dependence of energy on median frequency corresponding to individual AE clusters for(a)RD-and(b)TD-loading.

    Fig.8.Dependence of a relative AE cluster activity on tensile stress for(a)RD-and(b)TD-loading.

    The relative activity of individual clusters(background noise is not included)as a function of the applied tensile stress is shown in Fig.8.Almost immediately after the beginning of the tensile test,cluster 1(red color code)becomes dominant in both loading directions.The subsequent transition from cluster 1 to other clusters is realized at different stresses for the RD-and TD-loading and this transition differs in its course.In the RD-loading,the transition from cluster 1 is realized abruptly due to a concurrent activation of the other two clusters(clusters 2 and 3-green and blue color code,respectively,in Fig.8a)at a tensile stress of~80 MPa,while during the TD-loading,the increment of the cluster 2 is a gradual process(Fig.8b)which starts to occur at a stress of~50 MPa.

    In the case of the RD-loading,cluster 2 shows higher activity compared to cluster 3 in the stress range of 80–110 MPa,then loses its dominance between stresses of 110 and 140 MPa and finally becomes dominant again up to stresses beyond the YP.

    A different scenario was observed in the TD-loading,where cluster 1 loses its dominance gradually at the expense of a higher activity of cluster 2.Cluster 2 reaches its maximum activity around the tensile stress of 90 MPa,at which cluster 3 begins to be active.Then,both cluster 2 and cluster 3 show the competitive behavior with increasing tensile stress beyond the YP with higher activity of cluster 3 in the stress range of 110–145 MPa and at stresses above 230 MPa.The absolute dominance of cluster 2 at the expense of cluster 3 was observed in the stress range of 160-210 MPa.

    4.Discussion

    The homogeneous microstructure with a strong initial basal texture having a slightly wider spread of basal planes towards RD than TD(see Fig.1),is commonly observed in rolled Mg alloys[39–42].The homogenization temperature of 430°C was used to dissolve the intermetallic Mg17Al12phase and,therefore,the precipitates in the microstructure are assumed to be Al-Mn particles[2].With respect to the texture(Fig.1)andin-situSE imaging(Figs.5-6),dislocation glide has to be a preferential deformation mechanism during the tensile loading of the AM61 sheet.Twinning,as a competitive deformation mechanism in Mg alloys,can play an important role during the tensile loading of Mg sheets if the basal poles are extensively tilted from the normal plane of the sheet,as observed,for example in the ZE10 sheet[27].During loading in both directions up to 175 MPa(Figs.5-6),no evidence of twinning activity was found and,therefore,the SF analysis was performed only for theandslip systems(i.e.excluding twinning)in order to estimate their possible activity at the beginning of the deformation test,see Fig.2.

    A slight variation in the texture(distribution of the(0001)intensity along RD and TD,Fig.1)with a more favorable orientation of grains for the basal glide along the RD-loading(higher fraction of grains with SF in the range of 0.4–0.5,see Fig.2)results in a planar anisotropy of mechanical properties with higher tensile yield strength for the TD-loading(Table 1).

    The AE activity,observed almost immediately at the beginning of the test in both loading directions,indicates a long mean free path〈L〉due to a relatively low dislocation density in the initial microstructure caused by the solution heat treatment applied after rolling.

    According to Eq.1,the AE energy is proportional to the number of mobile dislocationsρmand a higher cumulative release of energy during tension along RD than TD is,therefore,caused by a more favorable orientation of the grains for easier basal glide in this loading direction(Fig.2).A local increase in the released AE energy in the stress range of 50–75 MPa for the RD-loading can be attributed to a higher collective dislocation movement compared to the TD-loading and can be explained(in terms of the SF analysis)by a higher number of well-oriented grains for basalslip in case of the RDloading.With the progress of the tensile stress,an increasing number of grains involved in the deformation contribute to the total AE energy in both loading directions.The following massive multiplication of dislocations leads to the maximal release of AE energy at the macroscopic YP.Beyond this point,increasing dislocation densityρtogether with increasing number of immobile dislocations results in the decreasing mean free path〈L〉of dislocations and,therefore,reduces the amount of released AE energy.

    On the other hand,the AE count rate is sensitive to changes in the dynamics of mobile dislocations.First,collectively moving dislocations produce detectable AE and a higher AE count rate in the RD-loading can be explained by a higher number of well-oriented grains for basal slip,which is typically the first active slip system in Mg alloys.The practically constant AE count rate(Fig.4)in the stress range of 25–75 MPa correlates with the fact that in crystalline materials showing no alternation of deformation mechanisms due to a change in the strain rate,a proportionality between the AE count rateand strain ratecan be described according to Fisher and Lally[43]as:

    Considering a constant strain rate during the entire tensile loading,the basalglide is the dominant deformation mechanism in the above-mentioned stress range.A similar,but not so pronounced,flat AE response for the Mg alloys before the YP was also reported for the AZ31 sheet[23].

    The following increase in the AE count rate at 75 MPa indicates an activation of further deformation mode,which is supposed to be the prismaticslip considering the CRSS required for its activation and the orientation of the grains with respect to the loading direction(SF analysis in Fig.2).The values of CRSS for the prismatic and pyramidal slip systems in Mg alloys are estimated to be around 70 MPa and 110 MPa,respectively[44,45].Increasing applied stresses give rise to the activation of pyramidalandslip resulting in the massive increase in the AE count rate above the stress of~115 MPa.The SF analysis(Fig.2)indicates that pyramidalslip is supposed to be more active during the TD-loading rather than in RD-loading,which can explain a more rapid increase in the AE count rate in the TD-loading.A continuously increasing dislocation density,especially of immobile dislocations after the YP,reduces the mean free path〈L〉of moving dislocations and mitigates the AE count rate.

    The activity of individual slip systems is also reflected in the median frequencyfmevolution.The constantfmin the stress range of 25–75 MPa(Fig 4.)can be explained by using the phenomenological relation forfm,presented for example in[33,46],

    and assuming that(i)the low density of mobile dislocationsρmat the beginning of the tensile loading of the heat-treated(i.e.stress-free)sheet should not significantly affect the mean free path〈L〉of moving dislocations and(ii)the mean velocity of dislocations〈v〉can be changed only slightly due to a constant strain rate.The subsequent decrease in the median frequencyfmat tensile stresses between 75 and 115 MPa can be attributed to a massive increase inρm,which subsequently reduces〈v〉.It was shown in[33]that during tensile loading of pure metals with a face-centred cubic(fcc)lattice(Al,Cu,Ni,Ag),fmdecreases from the beginning of the tests with increasing loading up to yielding and after this point it constantly increases,showing only a different slope for each investigated metal.In those pure fcc metals,deformation is realized via several equivalent slip systems(having the same Burgers vector and CRSS)and the minimum in the median frequencyfmoccurs at the YP.However,in hcp metals,the slip systems are not equivalent and require different conditions for their activation.Thus,the increase infmbefore the YP in the investigated Mg sheet can be associated with a rather significant transition between slip systems.A more pronounced drop in the case of the TD-loading is in good agreement with the SF analysis presented in Fig.2 indicating higher SF values for the activation of prismaticslip during the TD-loading compared to those observed during the RD-loading.Moreover,the activity of pyramidalslip is assumed with respect to CRSS for its activation[44,45].Furthermore,a higher activity of pyramidalslip during the TD-loading is expected due to a more significant increase infm,which is in good agreement with the SF analysis(Fig.2)showing a more preferable condition for the activation of pyramidalslip in the case of the TDcompared to RD-loading.

    The above-described activity of slip systems is strongly supported by results of the ASK analysis.All three clusters exhibit a similar shape(Fig.7)in the investigated time interval of 0–120 s and,consideringin-situSE imaging(Figs.5 and 6),the SF analysis(Fig.2)and previous studies on Mg alloys[36-38],they can be related to the dislocation slip activity.Therefore,the dominant mechanism during tensile loading along both RD and TD is dislocation glide.Cluster 1 dominates at lower tensile stresses in both loading directions and,considering the CRSS required for the activation of individual slip systems in Mg alloys,can be assigned to basalslip(Fig.8).The assignment of cluster 2 to prismatic and cluster 3 to the pyramidal slip system is based on the same idea taking into account the CRSS values for both slip systems[44,45].The dominance of theslip system up to the stress of~80 MPa during the RD-loading(Fig.8a)can be explained by a lower contribution of other slip systems to the deformation(SF analysis,Fig.2)and the principle of the ASK algorithm(revealing only the most dominant process in a given time interval).The subsequent activation of the prismatic and pyramidal slip systems represents a competitive process with a higher contribution of the prismatic slip up to a tensile stress of 110 MPa.The less favorable orientation of grains for the basal slip in the TD-loading(Figs.1 and 2)leads to the activation of prismaticslip(cluster 2,Fig.8b)at tensile stresses above 50 MPa and this slip system gradually takes control of the deformation.Interestingly,for both loading directions,the maximal activity of prismatic slip(cluster 2)is observed at the stress of 90 MPa and dominance of the pyramidalslip system(cluster 3)is seen in the stress range of 110–140 MPa.

    The uniqueness of the AE technique is manifested by its comparison with thein-situSE imaging,where the first slip traces on the surface of the sample,as the evidence of slip activity,were observed at stresses only above 125 MPa in the RD-loading.To date,the ASK analysis has been successfully used in the cases where deformation mechanisms show significant differences in their AE characteristics,such as slip vs.twinning or kinking[18,37,47,48].The presented results demonstrate that by applying the ASK algorithm on AE data it is also possible to reveal important real-time information about the activity of different dislocation slip systems from the entire sample volume leading,in turn,to the detailed determination of their contribution in the anisotropic tensile deformation behavior of the Mg-based sheets.

    5.Conclusions

    The dislocation activity in the rolled Mg-Al-Mn alloy sheet was investigated in tension along rolling(RD)and transversal directions(TD),complemented by the acoustic emission(AE)technique using the parametrization of the AE signal and its clustering analysis.The AE results were supported by the microstructure observation,particularly by the estimation of the Schmid factor(SF)for the activation of individual slip systems.The following conclusions can be drawn:

    ?Higher AE energy and count rate in the RD-loading than those observed in the TD-loading is caused by a higher number of favorably oriented grains for dislocation glide,i.e.,the higher number of AE sources.

    ?The constant AE count rate and the median frequency in the stress range of 25–75 MPa are related to a relatively low initial dislocation density in the heat-treated(i.e.stress-free)sheet for both loading directions.The subsequent increase in the AE count rate and decrease in the median frequency is explained in terms of an increasing number of moving dislocations influencing the mean velocity of dislocations.

    ?Increase in the median frequency before the macroscopic yield point(YP)indicates the transition between individual dislocation slip systems for both loading directions.Thein-situsecondary electron(SE)imaging excludes the twinning activity in the microstructure.

    ?The ASK(clustering)analysis supported byin-situSE imaging and the SF analysis reveals that basalslip dominates during tensile loading of the Mg sheet up to the stress of~80 MPa.It operates as a single deformation mode in the RD-loading,while in the TD-loading,at stresses above 50 MPa,it gradually loses its dominance at the expense of a higher activity of the prismaticslip system.

    ?For both loading directions,prismaticslip predominates in the stress range of 80–110 MPa,then it loses its dominance at the expense of pyramidalslip and finally,at stresses above 140 MPa,shows the absolute dominance in the RD-loading or exhibits higher activity compared to pyramidalslip in the TD-loading.

    Declaration of Competing Interest

    The authors declare that they have no conflict of interest.

    Acknowledgements

    This work was financially supported by the Czech Science Foundation(grant No.19-22604S),and the Operational Program Research,Development and Education,The Ministry of Education,Youth,and Sports(OP RDE,MEYS)(grant No.CZ.02.1.01/0.0/0.0/16_013/0001794).

    国产探花在线观看一区二区| 亚洲男人的天堂狠狠| 亚洲国产色片| 别揉我奶头 嗯啊视频| 亚洲欧美日韩高清在线视频| 熟女电影av网| 婷婷色综合大香蕉| 成人特级av手机在线观看| 俄罗斯特黄特色一大片| 久久九九热精品免费| x7x7x7水蜜桃| 天堂影院成人在线观看| 亚洲无线观看免费| 国产一级毛片七仙女欲春2| 97碰自拍视频| 深夜精品福利| а√天堂www在线а√下载| 久久久久九九精品影院| 欧美色欧美亚洲另类二区| 国产av一区在线观看免费| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 女人被狂操c到高潮| 99热这里只有是精品50| av福利片在线观看| 深夜a级毛片| 久久久久久久久大av| 偷拍熟女少妇极品色| 高清毛片免费观看视频网站| 如何舔出高潮| 尤物成人国产欧美一区二区三区| 日本熟妇午夜| 午夜福利在线在线| 亚洲精品粉嫩美女一区| 成年人黄色毛片网站| 欧美黄色淫秽网站| 免费观看精品视频网站| 夜夜看夜夜爽夜夜摸| 久久久久亚洲av毛片大全| 久久久久久久久大av| 欧美成人一区二区免费高清观看| 亚洲在线观看片| 九九在线视频观看精品| 国产免费av片在线观看野外av| 国产精品一区二区三区四区免费观看 | 国产在线男女| 99热这里只有是精品50| 国产精品乱码一区二三区的特点| 1000部很黄的大片| 欧美黑人欧美精品刺激| 18禁裸乳无遮挡免费网站照片| 成年免费大片在线观看| 十八禁国产超污无遮挡网站| 欧美精品啪啪一区二区三区| 午夜影院日韩av| 精品久久久久久久久av| 一个人免费在线观看的高清视频| 午夜福利视频1000在线观看| av黄色大香蕉| av福利片在线观看| 国产精品av视频在线免费观看| 午夜福利免费观看在线| 亚洲成人久久爱视频| 国内精品久久久久久久电影| 中文字幕av在线有码专区| 又紧又爽又黄一区二区| 久久精品国产亚洲av香蕉五月| 中文字幕熟女人妻在线| 国产美女午夜福利| 久久久久久久亚洲中文字幕 | 日本 欧美在线| 久久国产乱子免费精品| 97碰自拍视频| 免费在线观看日本一区| 老女人水多毛片| 日韩欧美在线二视频| 亚洲第一电影网av| 丰满人妻熟妇乱又伦精品不卡| 黄色丝袜av网址大全| 中文字幕精品亚洲无线码一区| 69人妻影院| 两个人视频免费观看高清| 国产高清有码在线观看视频| 99久久精品热视频| 免费观看人在逋| 精品人妻熟女av久视频| 成人美女网站在线观看视频| 国产一级毛片七仙女欲春2| 最新中文字幕久久久久| 最近最新免费中文字幕在线| 内地一区二区视频在线| 国产精品国产高清国产av| av天堂在线播放| АⅤ资源中文在线天堂| 麻豆成人午夜福利视频| 亚洲,欧美,日韩| 午夜老司机福利剧场| 亚洲美女黄片视频| 久久欧美精品欧美久久欧美| 哪里可以看免费的av片| 久久午夜亚洲精品久久| 中文在线观看免费www的网站| 色综合婷婷激情| 日韩人妻高清精品专区| 淫秽高清视频在线观看| 国产精品美女特级片免费视频播放器| 久久久久免费精品人妻一区二区| 亚洲精品成人久久久久久| av在线天堂中文字幕| 精品熟女少妇八av免费久了| 国产伦一二天堂av在线观看| 午夜福利欧美成人| 久久九九热精品免费| 欧美一级a爱片免费观看看| 12—13女人毛片做爰片一| 99在线人妻在线中文字幕| 午夜精品久久久久久毛片777| 丁香六月欧美| 成人三级黄色视频| 午夜老司机福利剧场| 精品一区二区三区av网在线观看| 国产爱豆传媒在线观看| 欧美日韩乱码在线| 日韩欧美在线二视频| 久久国产精品影院| 国产精品三级大全| 久久草成人影院| 欧美精品啪啪一区二区三区| 黄片小视频在线播放| 99在线视频只有这里精品首页| 美女被艹到高潮喷水动态| 大型黄色视频在线免费观看| 亚洲七黄色美女视频| 亚洲男人的天堂狠狠| 小说图片视频综合网站| 国产一区二区亚洲精品在线观看| 亚洲第一欧美日韩一区二区三区| 午夜两性在线视频| 国产精品爽爽va在线观看网站| 欧美丝袜亚洲另类 | 国产aⅴ精品一区二区三区波| 国产午夜福利久久久久久| 有码 亚洲区| 亚洲欧美激情综合另类| 一个人看的www免费观看视频| 欧美午夜高清在线| 久久6这里有精品| 免费人成视频x8x8入口观看| 免费观看精品视频网站| 首页视频小说图片口味搜索| 亚洲av成人精品一区久久| 成年版毛片免费区| 久久精品夜夜夜夜夜久久蜜豆| 琪琪午夜伦伦电影理论片6080| 亚洲av日韩精品久久久久久密| 国产欧美日韩精品亚洲av| a级毛片免费高清观看在线播放| 十八禁网站免费在线| 国产私拍福利视频在线观看| 亚洲片人在线观看| 久久精品人妻少妇| 我的老师免费观看完整版| 女同久久另类99精品国产91| 一本综合久久免费| 麻豆一二三区av精品| 亚洲国产精品久久男人天堂| 热99re8久久精品国产| 亚洲国产色片| 午夜精品在线福利| 18美女黄网站色大片免费观看| 日本免费a在线| 少妇丰满av| 自拍偷自拍亚洲精品老妇| 精品人妻偷拍中文字幕| 国产老妇女一区| 搞女人的毛片| 欧美高清成人免费视频www| 成人欧美大片| 欧美最黄视频在线播放免费| 91九色精品人成在线观看| 午夜福利高清视频| 国产免费av片在线观看野外av| 国产成人福利小说| 亚洲人成伊人成综合网2020| 亚洲第一欧美日韩一区二区三区| 午夜日韩欧美国产| 一进一出抽搐gif免费好疼| 我的女老师完整版在线观看| 久久天躁狠狠躁夜夜2o2o| 成年女人毛片免费观看观看9| 麻豆一二三区av精品| 在线观看av片永久免费下载| 亚洲美女视频黄频| 偷拍熟女少妇极品色| 精品人妻熟女av久视频| 国产高清三级在线| 伊人久久精品亚洲午夜| 日本黄大片高清| 欧美黄色片欧美黄色片| 狠狠狠狠99中文字幕| 一进一出抽搐动态| 九色国产91popny在线| 国产白丝娇喘喷水9色精品| 欧美在线黄色| 听说在线观看完整版免费高清| 精品日产1卡2卡| 又黄又爽又免费观看的视频| 午夜影院日韩av| 欧美一区二区亚洲| 麻豆成人午夜福利视频| 啦啦啦观看免费观看视频高清| 变态另类成人亚洲欧美熟女| 久久中文看片网| 免费人成视频x8x8入口观看| 婷婷丁香在线五月| 黄色配什么色好看| 搡老岳熟女国产| 国产伦人伦偷精品视频| 精品午夜福利视频在线观看一区| 一个人看的www免费观看视频| 一级a爱片免费观看的视频| 黄片小视频在线播放| 久久性视频一级片| 最近最新中文字幕大全电影3| 国产免费男女视频| 动漫黄色视频在线观看| 国产黄a三级三级三级人| 熟女人妻精品中文字幕| 99久久99久久久精品蜜桃| 我的老师免费观看完整版| 久久久久久国产a免费观看| 免费在线观看成人毛片| 欧美激情国产日韩精品一区| 久久国产乱子免费精品| 床上黄色一级片| 别揉我奶头~嗯~啊~动态视频| 午夜免费男女啪啪视频观看 | 性色avwww在线观看| 一进一出抽搐动态| 久久国产精品人妻蜜桃| 在线观看66精品国产| 3wmmmm亚洲av在线观看| 日本成人三级电影网站| 国产高潮美女av| 亚洲国产欧美人成| 日韩高清综合在线| 欧美3d第一页| 国产精品一区二区免费欧美| av在线蜜桃| 日本精品一区二区三区蜜桃| 美女cb高潮喷水在线观看| 能在线免费观看的黄片| 人妻丰满熟妇av一区二区三区| 欧美最新免费一区二区三区 | 精品久久久久久久久av| 禁无遮挡网站| 韩国av一区二区三区四区| 日韩精品青青久久久久久| 亚洲av第一区精品v没综合| 日韩欧美国产在线观看| ponron亚洲| 国产精品自产拍在线观看55亚洲| 成年女人毛片免费观看观看9| 国产乱人伦免费视频| 免费av毛片视频| 国产成人福利小说| 国产伦人伦偷精品视频| 老司机午夜十八禁免费视频| 真实男女啪啪啪动态图| 午夜福利在线观看免费完整高清在 | 人妻制服诱惑在线中文字幕| www日本黄色视频网| 亚洲黑人精品在线| 高清在线国产一区| 99视频精品全部免费 在线| 国产高清视频在线观看网站| 一本一本综合久久| 看免费av毛片| 波多野结衣高清无吗| 免费大片18禁| 精品久久久久久久人妻蜜臀av| 色综合亚洲欧美另类图片| 成年免费大片在线观看| 亚洲精品在线观看二区| 日韩有码中文字幕| 三级国产精品欧美在线观看| 日韩高清综合在线| 精品不卡国产一区二区三区| 黄色视频,在线免费观看| 性欧美人与动物交配| 人妻久久中文字幕网| 欧美成人性av电影在线观看| 日韩 亚洲 欧美在线| 国产亚洲av嫩草精品影院| 两人在一起打扑克的视频| netflix在线观看网站| 又粗又爽又猛毛片免费看| 久久久久久大精品| 99精品在免费线老司机午夜| 男女床上黄色一级片免费看| 精品久久久久久久久久免费视频| 亚洲av美国av| 国产大屁股一区二区在线视频| 在线观看一区二区三区| 国产v大片淫在线免费观看| 亚洲精品亚洲一区二区| 日韩有码中文字幕| 亚洲专区中文字幕在线| 男女床上黄色一级片免费看| 欧美色欧美亚洲另类二区| 一区二区三区激情视频| 禁无遮挡网站| 国产一级毛片七仙女欲春2| 欧美成人a在线观看| 亚洲va日本ⅴa欧美va伊人久久| 听说在线观看完整版免费高清| 国产激情偷乱视频一区二区| 国产成人福利小说| 99久久精品一区二区三区| 精品乱码久久久久久99久播| 中文在线观看免费www的网站| 久久亚洲精品不卡| 亚洲av不卡在线观看| 自拍偷自拍亚洲精品老妇| 午夜激情福利司机影院| 又粗又爽又猛毛片免费看| 午夜精品久久久久久毛片777| 久久精品91蜜桃| 亚洲五月天丁香| 亚洲va日本ⅴa欧美va伊人久久| 极品教师在线免费播放| 亚洲,欧美,日韩| 欧美性猛交╳xxx乱大交人| 波多野结衣巨乳人妻| 亚洲内射少妇av| 久久久久久久亚洲中文字幕 | 精品99又大又爽又粗少妇毛片 | 51午夜福利影视在线观看| 欧美高清成人免费视频www| 麻豆国产97在线/欧美| 亚洲av成人精品一区久久| 国产精品,欧美在线| 日本三级黄在线观看| 国产一区二区在线观看日韩| 色综合欧美亚洲国产小说| 91av网一区二区| 国产精品亚洲av一区麻豆| 欧美又色又爽又黄视频| 亚洲欧美激情综合另类| 亚洲成人精品中文字幕电影| 特大巨黑吊av在线直播| 国产一区二区在线观看日韩| 欧美高清性xxxxhd video| av欧美777| 精品日产1卡2卡| 18禁黄网站禁片免费观看直播| 99热这里只有是精品在线观看 | 亚洲男人的天堂狠狠| 91av网一区二区| 蜜桃亚洲精品一区二区三区| 中文字幕人成人乱码亚洲影| 精品久久久久久成人av| 色综合亚洲欧美另类图片| 精品久久久久久久久久久久久| 久久精品夜夜夜夜夜久久蜜豆| 亚洲avbb在线观看| 亚洲综合色惰| av在线老鸭窝| 国产欧美日韩精品亚洲av| 脱女人内裤的视频| 国产一区二区激情短视频| 脱女人内裤的视频| 一个人免费在线观看的高清视频| 国产精品不卡视频一区二区 | 一级黄色大片毛片| 亚洲精华国产精华精| 久久久久免费精品人妻一区二区| 亚洲av美国av| 午夜精品久久久久久毛片777| 国产美女午夜福利| 能在线免费观看的黄片| 亚洲人成网站高清观看| 日韩中字成人| 男插女下体视频免费在线播放| 怎么达到女性高潮| av国产免费在线观看| 1024手机看黄色片| 搡老熟女国产l中国老女人| 亚洲片人在线观看| 免费在线观看亚洲国产| 久久香蕉精品热| 精品免费久久久久久久清纯| 免费人成在线观看视频色| 国产蜜桃级精品一区二区三区| 麻豆国产av国片精品| 日本一本二区三区精品| 男女下面进入的视频免费午夜| 在线看三级毛片| 精品一区二区三区人妻视频| 久久久久久久精品吃奶| 99热精品在线国产| 色综合亚洲欧美另类图片| 欧美中文日本在线观看视频| 身体一侧抽搐| 欧美激情在线99| 日韩欧美在线乱码| 我的女老师完整版在线观看| 精品午夜福利视频在线观看一区| 免费在线观看影片大全网站| 久久久国产成人免费| 久久精品人妻少妇| 老司机午夜福利在线观看视频| 精华霜和精华液先用哪个| 亚洲av不卡在线观看| 99久久精品热视频| 中文字幕久久专区| 亚洲av电影在线进入| 欧美激情国产日韩精品一区| 精品久久久久久久末码| 国产精品一区二区性色av| 99精品久久久久人妻精品| 久久热精品热| 欧美丝袜亚洲另类 | 国产av不卡久久| 国产激情偷乱视频一区二区| 天美传媒精品一区二区| 日本成人三级电影网站| 成人性生交大片免费视频hd| 亚洲自拍偷在线| 亚洲第一欧美日韩一区二区三区| 一进一出好大好爽视频| 国产色婷婷99| or卡值多少钱| 在线十欧美十亚洲十日本专区| 色精品久久人妻99蜜桃| 精品国产亚洲在线| 国内精品久久久久久久电影| 少妇被粗大猛烈的视频| 久久久久久久亚洲中文字幕 | 婷婷亚洲欧美| 欧美成狂野欧美在线观看| 久久精品人妻少妇| 精品日产1卡2卡| 国产精品一及| 亚洲中文日韩欧美视频| 欧美激情在线99| 国产亚洲精品久久久久久毛片| 久久性视频一级片| av天堂在线播放| 午夜福利成人在线免费观看| 亚洲av免费在线观看| 中文字幕高清在线视频| 日本a在线网址| 国产老妇女一区| 亚洲人成电影免费在线| 美女大奶头视频| 欧美极品一区二区三区四区| 日本精品一区二区三区蜜桃| 日日摸夜夜添夜夜添av毛片 | 国产av一区在线观看免费| 亚洲av成人av| 国产精品99久久久久久久久| 少妇的逼好多水| 窝窝影院91人妻| 嫩草影院精品99| 亚洲欧美日韩东京热| 国产精品久久视频播放| 在线观看av片永久免费下载| 精品日产1卡2卡| 欧美午夜高清在线| 国产精品电影一区二区三区| 最近视频中文字幕2019在线8| 国产精品乱码一区二三区的特点| 日本一二三区视频观看| 又粗又爽又猛毛片免费看| 特大巨黑吊av在线直播| www.熟女人妻精品国产| 欧美黑人欧美精品刺激| 中文资源天堂在线| 欧美黑人巨大hd| 免费在线观看影片大全网站| 成年女人毛片免费观看观看9| 亚洲激情在线av| 精品久久久久久久人妻蜜臀av| 国产高清有码在线观看视频| 亚洲中文字幕一区二区三区有码在线看| 国内精品久久久久久久电影| 久久精品国产亚洲av香蕉五月| 观看免费一级毛片| 极品教师在线免费播放| 欧美区成人在线视频| 亚洲人与动物交配视频| 亚洲精品成人久久久久久| 男女那种视频在线观看| 亚洲熟妇中文字幕五十中出| 欧美绝顶高潮抽搐喷水| 成年人黄色毛片网站| av在线蜜桃| 久久精品国产99精品国产亚洲性色| 久久婷婷人人爽人人干人人爱| 成人一区二区视频在线观看| 天美传媒精品一区二区| 国产极品精品免费视频能看的| 在线免费观看不下载黄p国产 | 乱人视频在线观看| 婷婷精品国产亚洲av| 亚洲最大成人手机在线| 国产真实伦视频高清在线观看 | 精品一区二区三区视频在线观看免费| x7x7x7水蜜桃| 国产欧美日韩一区二区三| 免费观看精品视频网站| 国内毛片毛片毛片毛片毛片| 一区二区三区激情视频| 性欧美人与动物交配| 99热只有精品国产| 久久久久久久久久成人| 欧美黑人巨大hd| 91字幕亚洲| 亚洲av不卡在线观看| 欧美精品国产亚洲| 亚洲国产精品合色在线| 久久性视频一级片| 久久国产乱子免费精品| 中文字幕人妻熟人妻熟丝袜美| 变态另类丝袜制服| 伦理电影大哥的女人| 国产欧美日韩一区二区精品| 亚洲精品日韩av片在线观看| 人妻丰满熟妇av一区二区三区| 51午夜福利影视在线观看| 不卡一级毛片| 成人鲁丝片一二三区免费| 男人舔女人下体高潮全视频| 色av中文字幕| 香蕉av资源在线| 午夜免费男女啪啪视频观看 | 久久精品国产自在天天线| 免费在线观看日本一区| 日日干狠狠操夜夜爽| 国内精品美女久久久久久| 婷婷丁香在线五月| 色哟哟·www| 亚洲av二区三区四区| 日韩有码中文字幕| 黄色女人牲交| 久久久久久久久中文| 欧美黑人巨大hd| 国产精品久久久久久亚洲av鲁大| 婷婷丁香在线五月| 国产一区二区激情短视频| 国产精品98久久久久久宅男小说| 国产中年淑女户外野战色| 国产私拍福利视频在线观看| 熟妇人妻久久中文字幕3abv| 久久久国产成人免费| 午夜a级毛片| 人人妻人人看人人澡| 一个人观看的视频www高清免费观看| 国产高清视频在线播放一区| 久久精品国产亚洲av天美| 五月玫瑰六月丁香| 亚洲av五月六月丁香网| 亚洲成av人片免费观看| 亚洲欧美日韩无卡精品| 夜夜躁狠狠躁天天躁| 丝袜美腿在线中文| 久久热精品热| 中文字幕人妻熟人妻熟丝袜美| 深爱激情五月婷婷| 国产极品精品免费视频能看的| 老司机福利观看| a级毛片免费高清观看在线播放| 少妇高潮的动态图| 国产人妻一区二区三区在| 精品国产三级普通话版| 91在线观看av| 欧美成狂野欧美在线观看| 一a级毛片在线观看| 久久久久久久久中文| 日本一二三区视频观看| 午夜福利欧美成人| 国产精品人妻久久久久久| 国产精华一区二区三区| 91av网一区二区| 精品久久久久久久人妻蜜臀av| 日韩高清综合在线| 亚洲成人中文字幕在线播放| 脱女人内裤的视频| 中国美女看黄片| 国产野战对白在线观看| 夜夜躁狠狠躁天天躁| 老鸭窝网址在线观看| 村上凉子中文字幕在线| 人妻丰满熟妇av一区二区三区| 毛片一级片免费看久久久久 | 久久中文看片网| 人妻制服诱惑在线中文字幕| 美女 人体艺术 gogo| 亚洲精品色激情综合| 人妻制服诱惑在线中文字幕| 在线免费观看不下载黄p国产 | 久99久视频精品免费| 国产精品av视频在线免费观看| 欧美日韩国产亚洲二区| 亚洲人成电影免费在线| 我的女老师完整版在线观看| 757午夜福利合集在线观看| 一级黄片播放器| 美女高潮喷水抽搐中文字幕| 国产白丝娇喘喷水9色精品| 最新在线观看一区二区三区| 有码 亚洲区| 一区二区三区激情视频| 免费黄网站久久成人精品 | 在线a可以看的网站|