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

    In-situ analysis of slip transfer and heterogeneous deformation in tension of Mg-5.4Gd-1.8Y-1.5Zn alloy

    2021-01-04 04:55:36ShaosongJiangYongJiaXiaojunWang
    Journal of Magnesium and Alloys 2020年4期

    Shaosong Jiang,Yong Jia,Xiaojun Wang

    School of Materials Science and Engineering,Harbin Institute of Technology,Harbin,150001,China

    Received 25 July 2019;received in revised form 27 November 2019;accepted 13 January 2020 Available online 2 June 2020

    Abstract Slip transfer is influential in determining damage nucleation of polycrystalline material.The interactions between dislocations and grain boundaries(GBs)was investigated using in-situ tension test in a multi-directionally forged Mg-5.4Gd-1.8Y-1.5Zn(wt%)alloy.It was found that strain accommodation of individual grains by means of slip occurred more easily than slip transfer when several slip systems were operable.The basal-basal slip transfer occurred when the GB misorientation was smaller than 34.2°,whereas basal-pyramidal type took place when the crystallographic misorientation was larger than 48.8°.The product of Luster-Morris mfactor and the sum of the Schmid factors of the two correlated slip systems indicated that the threshold for basal-basal slip transfer may exist,however,basal-pyramidal slip transfer shows no such threshold and is more complicated.These results presented here demonstrated that besides the geometrical alignment,the deformation details(such as the number of operable slip systems)and stress state in each individual grain must be considered.

    Keywords:Slip transfer;Mg-RE alloy;Heterogeneous deformation;In-situ tensile test.

    1.Introduction

    Heterogeneous deformation of polycrystals has been an intense field of research[1,2].During deformation,the mechanical integrity and continuity of GBs are maintained,and abundant studies have proven that GBs act as effective barriers for dislocation motion[3-7].It is also possible that plastic deformation takes place across GBs through the traverse of dislocations[1,8-11].High-angle or low-angle GBs,which facilitates or hinders the mobility of dislocations,will determine the mechanical behavior of polycrystalline metals[12,13].The strain field,accommodated by the GBs,will influence the boundary energy and,consequently,cohesion and nucleation of the crack[8,14-16].Thus,a thorough understanding of the localized deformation at or near GBs will be necessary to establish an accurate model capable of predicting the damage nucleation.

    Full-field homogenization in crystal plasticity modeling,using crystal plasticity finite element method(CPFEM)or spectral methods by fast Fourier transformation(FFT),has the potential to accurately predict the mechanical properties of polycrystalline metals[17,18].The principal advantage of CPFEM or spectral solver would be their great flexibility to handle the complicated internal and/or external boundary condition.However,simulations generally show partial agreements with experimental measures and exhibit rather shaper orientation gradients near GBs[19-21].This is because one consequence of the crystalline anisotropy is the spatially heterogeneous boundary conditions,resulting from stress or strain concentrations that normally precedes the onset of fracture.Dislocation pile-ups and geometrically necessary dislocations at GBs,which accommodates the strain gradients,is generally ignored in crystal plasticity modeling.Therefore,to predict effective materials properties,strain field interactions between dislocations and GBs must be quantitatively understood and imported into the crystal plasticity modeling.

    Fig.1.Geometric relationship between the incoming slip system and outgoing slip system across GBs for a bicrystal

    To introduce the grain boundaries properties into the mesoscale model,a clear understanding of the details of the microstructural hierarchy and stress state near GBs is necessary.Slip transfer process across GBs has a double-sided effect on mechanical properties as this phenomenon promotes or prevents crack initiation.M.Jiménez et al.[22]found that the short cracks propagated along the easy slip transfer path across the GBs in Ni-based superalloys.Britton et al.[11]measured the stress field at the head of the blocked slip band using high-angular resolution electron backscatter diffraction(HREBSD)and found that blocked slip bands with or without stress concentration are instructive for studying fracture initiation and propagation mechanism.Wang et al.[23]reported that slip transfer enables mechanical twinning to occur in an adjacent grain despite the fact that the twins have a very low Schmid factor relative to the global(uniaxial tensile)stress tensor in commercial purity titanium.

    In the past research,multiple criteria have been proposed to assess the ease of slip transfer in polycrystalline metals,including geometric parameter and stress criteria.The geometric relationship of the slip transfer is indicated in Fig.1.Themparameter,from Luster and Morris[24],considers the coplanarity of incoming and outgoing slip system and is promising to predict the slip transfer:

    whereψrepresents the angle between a pair of slip planes normal andκis the angle between the incoming and outgoing slip directions.Themfactor lies in the range of 0~1,where 1 means that dislocations could traverse GBs without resistance and 0 represents GB is impenetrable.Nonetheless,in-situstraining and diffraction-contrast electron tomography in TEM experiments show that the interactions between lattice dislocations and GBs have a rather complicated process and that the emitted slip system is governed by the magnitude of the residual boundary dislocation[25].Atomic and molecular dynamics simulations demonstrate that many other factors,such as local resolved stress,boundary structures,and the boundary strain energy density,also determine the slip transmission process[26-28].

    While most of the studies focus on the nucleation of dislocation at the head of blocked dislocations,there is very limited work investigating the threshold criterion for slip transmission.Hémery et al.[29]researched slip transfer between primaryαnodules in Ti-6Al-4V and confirmed thatmparameter is a convincing indicator of the probability of slip transfer.Bieler et al.[30]investigated the slip transfer cases within near-cube oriented grains in pure Al and found that,besides the geometric alignments between slip systems,the deformation details in each grain also play an important role.The purpose of this work is to extract the threshold for slip transfer and analyze the heterogeneous deformation near GBs in a Mg-5.4Gd-1.8Y-1.5Zn alloy.Such work would contribute to crystal plasticity modeling used to address the deformation kinematics.

    2.Experimental procedure

    The as-cast magnesium alloy with a nominal composition of 5.4 wt% Gd,1.8 wt% Y,and 1.5 wt% Zn was provided from Shenyang Institute of Metal Research with coarse grain diameter about 220μm.To refine the grain size,the alloy was forged.Before forging the ingots were homogenized at 500°C for 12 h followed by rapid cooling in the water.As general hot-deformation process,like rolling and extrusion,produces strong basal texture and the multi-directional forging results in weak texture,the alloy was multi-directionally forged at 450°C by six passes.Thus,the effect of the texture on the initiation of the slip systems is excluded from the alloy.After forging,the wrought alloy was heated up to 210°C for 24 h to achieve peak aging.

    The dog-bone tensile specimens having a thickness of 0.8 mm,gauge width of 1.8 mm,and gauge length of 7 mm,were cut from the aged alloy by electrical discharge machining.Specimens surface were ground using silicon carbide paper with various grit sizes and then polished using colloidal silica suspension.To obtain a stress-free surface for electron backscatter diffraction(EBSD)characterization,the specimen surface was electro-polished at a voltage of 10 V and the temperature of-20 °C using an electrolyte solution of 10 vol.%perchloric acid with 90 vol.% alcohol for 100 s.The grain orientation of the tensile specimen was measured using EBSD before and afterin-situtension in a field-emission gun scanning electron microscope(FEG-SEM),a Nova NanoSEM 430 instrument(FEI Company,Hillsboro,OR),equipped with a Nordlys EBSD detector(Oxford Instruments,Abingdon,Oxfordshire).Images were made using an acceleration voltage of 20 kV and a current of 6.1 nA.Automated beam scanning on the sample surface covered an area of 0.8 mm2with a step size between points of 0.9μm.Afterin-situtension,the number of nonindexed or incorrectly indexed electron backscatter patterns(EBSPs)increased,and therefore a medium level of zero-solution extrapolation was performed by Channel 5 software.

    Table 1Possible slip systems for the alloy;the same coloring scheme represents the same family of crystallographic planes and is used throughout this study.

    A tension test was carried out at room temperature within the NanoSEM 430 instrument chamber,using a voltage of 10 kV,which is sensitive to surface topography.The dog-bone specimen was stressed at a strain rate of 1×10?3s?1and was mounted horizontally on the 2 kN tensile stage(Deben,London,UK)driven by a dual threaded leadscrew to allow real-time observation of the gauge section.The displacement was stopped at different strains up to 22% to monitor slip bands.

    The crystallographic orientation data was imported into a MATLAB toolbox[31]to calculate slip lines on the sample surface and themfactor from Luster-Morris.The MATLAB toolbox also plots pairs of cell lattice orientations with highmvalues around GBs,corresponding to the appropriate slip system and slip traces.The computed slip traces,according to the cross product of slip plane normal and sample surface normal(i.e.T=nsp×e,whereTis the slip bands vector,nsprepresents the slip plane normal vector in sample coordinate system,andeis[001]),is compared to the SEM micrograph of slip lines.Theoretical slip lines that deviate from experimental observation by more than 5° are not considered in the following analysis.The possible slip systems initiated during tension are listed in Table 1 and are colored according to three different hexagonal close-packed(HCP)slip planes.The basal,prismatic,and pyramidal slip traces are colored blue,red,and gold respectively throughout this study.

    3.Results

    3.1.Microstructure evolution

    The inverse pole figures(IPFs)in the tensile direction(X direction)and the corresponding pole figures(PFs)of the sample at the same area before tension and after~22% strain are shown in Fig.2.The arrow indicates the same grain in Fig.2a and Fig.2c.It should be mentioned that the rumpled surface after straining affects the EBSP data acquisition,but does not influence the precision of the orientation data.The constituent grains in Fig.2a have a wide distribution of sizes,from 10μm to 100μm.The orientation of each grain is displayed in different colors in IPFs,and as shown in Fig.2c,the fraction of<10-10>fiber texture colored in blue increases with strain.After straining,the grains have developed a sharper texture(as indicated by the circular inset in Fig.2d)where four prismatic plane poles are distinct in the deformed pole figure,indicating that prismatic slip has rotated the grains to this more stable orientation[32].In addition,there is an increase in the strength of the clustering of<0001>poles perpendicular to X-axis in Fig.2d(from 5.6 to 6.5 times random).The poles of the<0001>texture component in PF that are parallel to the X direction(tensile direction),which are illustrated by the square inset in Fig.2b,rotate nearly 90° toward the Y direction in Fig.2d,suggesting the formation of tensile twins(i.e.<-1011>{10-12}tensile twins rotate<2-1-10>axis by 86.2°).

    3.2.Slip transfer analysis

    3.2.1.Analysis strategy

    The portion of the tensile specimen surface after straining is shown in Fig.3,and it appears that heterogeneous deformation occurs for most of the grains,as evidenced by welldefined boundary ledges and rough topography.Slip lines are visible on the sample surface and are kinked or curved due to the surface rumpling,which would influence the interpretation of slip traces.It seems that a single slip system dominates the deformation within each grain and the operated slip system number varies from grain to grain.

    Since it is impossible to monitor the slip lines evolution in real-time within the SEM.The identification of the incoming or outgoing slip assumes that if the basal-basal slip transfer occurs,the one with the higher Schmid factor is identified as the incoming slip.If the basal-pyramidal slip occurs,the basal slip was considered as the incoming slip.This is because at room temperature the CRSS of basal slip is much smaller than the pyramidal slip with typical CRSSnon-basal/CRSSbasalratios ranging between 20 and 100.The evaluation of themfactor associated with the pair of slip systems at GB 1 is given using the method illustrated in Table 2.The coloring scheme of slip systems is identical to slip planes in Table 1.The slip that commenced on a specific slip system is determined by the critical resolved shear stress(CRSS)and the Schmid factor,and thus,the slip systems are ranked according to the Schmid factor in the table ofmfactor.It should be noted that as the local stress tensors may vary substantially from grain to grain,so the ranking of the Schmid factor is only approximate.The leftmost column of Table 2 notes the slip system number among all potential slip systems(3{0001},3{10-10},and 6{11-22})for grain 565,and the magnitude of the Schmid factor for these slip systems are arranged in decreasing order,which is listed in the second column.Similarly,the slip system number in the top row of Table 2 for grain 568 are sorted by corresponding Schmid factor.For a bicrystal,twelve slip systems within each grain results in a 12×12 array with 144mfactors,and eachmvalue involves the geometry of two slip systems.Slip systems with Schmid factor larger than 0.25 are considered in Table 2,which produces 24mvalues;mfactors associated with the observed basal and pyramidal slip systems in Fig.3a are shown in bold.

    Fig.2.EBSD inverse pole figure(IPF)in the tensile direction for the tensile specimen(a)before straining and(c)after~22% strain.The corresponding pole figures of(b)undeformed and(d)strained sample are adjacent.The fraction of the<10-10>texture component colored in blue in Fig.2c increases with strain,indicating the occurrence of the prismatic slip.

    Table 2mvalue related to the geometry of two slip systems at GB 1;the slip systems are ranked according to the Schmid factor.The mfactor associated with slip lines near GB 1 are highlighted in bold.

    Table 2mvalue related to the geometry of two slip systems at GB 1;the slip systems are ranked according to the Schmid factor.The mfactor associated with slip lines near GB 1 are highlighted in bold.

    ?

    The magnitude of the Schmid factors for grain 568 ranges between 0.297 and 0.48,while grain 565 has only four slip systems with Schmid factor greater than 0.25.Although basal slip dominates the deformation at room temperature,nonbasal and cross-slip are also observed in polycrystalline magnesium alloys.Slip systems having lower Schmid factors are also possible when the CRSS reaches some critical value.One way of identifying a meaningfulmfactor is to check the upper left portion of Table 2,as this part constitutes slip systems interactions that are likely to act.This method developed by Bieler is useful in understanding the potential for slip transfer[30].

    In the following analysis,strain field interactions between slip bands and GBs are categorized as follows:(a)slip transfer is not observed,where sharp topography develops at GBs region;(b)slip transfer does happen without ledges or with faint topography.

    3.2.2.Slip transfer does not happen

    As mentioned above,the surface topography of grain 565 and 568 are shown in the scanning electron micrograph of Fig.3a.Due to the curved topography,slip bands of grain 568 follows the“C-type”path,and the curvature of slip lines increases as the bands approach grain boundaries.By virtue of the orientation data,it can be known that for grain 568,the basal plane is nearly parallel to the tensile stress axis,and for this extreme orientation the slip deformation forms as steps on the sample surface.The basal plane of grain 565 is oriented nearly 45°to the applied stress direction,and,therefore,dense slip lines that are parallel to one another are observed.The slip trace calculations based upon the orientation data indicate that basal slip commences in grain 565,while slip operates on the pyramidal plane in grain 568.Slip lines on the sample surface are closely aligned with the computed slip traces,as illustrated in Fig.3a.In grain 565,slip occurs along<2-1-10>-type directions on the(0001)plane(slip system 1 or 2),as indicated by the blue dashed line;while the(2-1-12)[-2113]slip system,represented by the gold line in Fig.3,commences in grain 568(slip system 10).

    Fig.3.(a)Surface topography of the tensile specimen after straining,where it appears that a single slip system is dominant for most of the grains.The operated slip systems in grain 565 and 568 are represented by the blue and gold lines;The number in the circle denotes the slip system in Table 1.(b)Visualization of the bold mvalue in Table 2 associated with the geometry of two slip systems in the hexagonal unit cell:the blue and gold lines represent the basal and pyramidal slip traces,respectively.Because of the distinct GB ledge,slip transfer across GB 1 is not observed.

    The possible slip traversing parameters for grain 565 and 568 are listed in Table 2.The highestmfactor,0.91,accompanies with second highly stressed slip system in grain 565,but the other slip system(slip system 12 ranked 5th)is not observed in the SEM photograph for grain 568.The computed slip traces at boundary 1 result in themfactor of 0.7 or 0.81,which are associated with the geometry of two slip systems in the hexagonal unit cell as visualized in Fig.3b.Themvalues of 0.7 or 0.81,resulting from the most highly stressed slip system in grain 565(slip systems 2 and 1),correspond to slip system 10 in grain 568,which has the 3rdhighest value of the Schmid factor.However,the well-defined boundary ledge at boundary 1 indicates heterogeneous strain between grain 565 and 568,such that slip transfer seems unlikely.

    As shown in Fig.4,it appears that dislocations traverse GB 2 during deformation;rather,the slip lines on either side of the GB 2 are misaligned with each other and the existence of a boundary ledge indicates heterogeneous strain between grain 456 and 457.Two sets of parallel yet intersecting sets of slip lines can be observed in grain 456;slip trace calculations based on orientation data predict that(-2112)[2-1-13]pyramidal slip and basal slip operate in these two grains,as indicated in Fig.4a.It is noteworthy that grain 456 is oriented in an unfavorably orientation-that is,the basal plane is parallel to the applied stress and in this condition,the Schmid factor of basal slip system approaches zero,as listed in the last row of Table 3.Still,slip occurs within{0001}plane,which means that stress tensors in the alloy are relatively complex(thus,as defined in section 3.2.1,the ranking of the Schmid factor is only approximate).In addition to the distinct boundary ledge,given that the poor geometric alignment between two operated slip system with themvalue of 0.6,slip transfer across boundary 2 is highly unlikely.

    Table 3Schmid factor values for grain 456 and 457 and the mfactor for boundary 2.The mfactor 0.6,which is typeset in bold,corresponds to the operated slip system in grain 456 and 457.

    Table 3Schmid factor values for grain 456 and 457 and the mfactor for boundary 2.The mfactor 0.6,which is typeset in bold,corresponds to the operated slip system in grain 456 and 457.

    ?

    Table 4mfactor between potential slip systems in grain 335 and 334 and corresponding Schmid factors.Two sets of operated slip systems produce mfactor 0.77 or 0.85.

    Table 4mfactor between potential slip systems in grain 335 and 334 and corresponding Schmid factors.Two sets of operated slip systems produce mfactor 0.77 or 0.85.

    ?

    3.2.3.Slip transfer occurs convincingly

    Fig.5 demonstrates that slip lines traverse boundary 3 with misorientation angle of 76° as the GB is not effective in prohibiting dislocation motion.The GB seems transparent,and in this case no boundary ledge forms,indicating homogeneous strain in both grains.It seems that slip lines in grain 335 could be from either the slip system 10 or 7;the upper right portion of slip bands is more consistent with the slip system 10,whereas the lower parts are nearly aligned with slip system 7.Also,the slip systems 10 and 7 in grain 335 have similarly high Schmid factors;however,the slip system 10 is highly misaligned with the initiated slip system 3 in grain 334(m=0.34),as underlined in Table 4.The computed slip system traces in grain 335 disagree with observations because the topography of the sample surface is sloped toward the lower left corner in Fig.5a,such that the inclined slip lines for slip system 7 in grain 335 look more like slip system 10.The geometrical configurations between observed slip systems 7 and 3 in grain 335 and 334 results inmfactor,0.77,which is shown in bold in Table 4,and both of the slip systems correspond to the highest Schmid factor.Considering the lower value of the Schmid factor for basal slip system 2,which has 7thhighest Schmid factor in grain 334,and complex stress tensors during tension,another m?factor,0.85,is also possible and is illustrated in Fig.5c.

    Fig.4.(a)Scanning electron micrograph showing slip lines at boundary 2 with misorientation 68.2° between grain 456 and 457.It appears that slip transfer occurs at boundary 2,however distinct boundary ledge indicates heterogeneous deformation between two grains.(b)The slip trace calculations,in comparison to slip lines in Fig.4a,demonstrate that pyramidal and basal slip occur in two grains as indicated by the dashed lines,both of which correspond to highest Schmid factor.The poor alignment between two slip systems at boundary 2 and the boundary ledge show that slip transfer across boundary 2 is prohibited in this case.

    Fig.5.(a)Scanning electron micrograph overlaid with IPF map shows continuous slip or dislocation motion through boundary 3 with misorientation 76°between grain 334 and 335.(b)The mfactor is obtained by geometrical relationships between observed slip systems at boundary 3 with values 0.77 or 0.85.

    Slip transfer initiated from basal to basal plane in a bicrystal is also observed in this study,where the blue slip system 3 traverses the boundary 4 with themfactor of 0.99,as illustrated in Fig.6a.Grain 455 and 459 having an angle of misorientation 3.7°is oriented unfavorably,in which the basal plane is almost perpendicular to the applied stress direction,such that the CRSS has small value on{0001}plane.As shown in Fig.6b,the magnitude of the Schmid factor for the operated basal slip system in two grains has a lower ranking.Actually,slip on a pyramidal plane is most favorable relative to the tensile axis,but pyramidal slip system is not observed in these grains,illustrating the significance of different CRSS values for different slip systems.

    Fig.6.(a)Scanning electron micrograph overlaid with IPF map showing slip transfer takes place at boundary 4 with misorientation 3.7°,in comparison with GB 3 illustrated above.(b)The basal slip system operates in grain 455 and 459,resulting in the mfactor of 0.99.

    3.2.4.Statistics from 94 observations

    As mentioned previously,slip transfer across GBs seems independent of the crystallographic misorientation.In this section,statistical results in conjunction with 94 grain boundaries are used to identify the relationships between the slip transfer events and geometrical configuration at the common boundary,as shown in Fig.7,where the relationships are plotted asmvalue versus grain boundary misorientation.Slip transfer from basal to prismatic plane is not observed in this study,and the observed slip transfer events from basal to basal and basal to pyramidal planes are represented by blue and gold triangles respectively.In addition,the maximummvalue in each case is displayed in solid squares in a similar way,indicating that the observed cases are sometimes different from the maximum m?value.The×symbols denote the cases when slip transfer is prohibited at the boundary,and the relevant maximummparameter is also noted in the figure.Statistical results in Fig.7 indicate that a highmvalue,larger than 0.9,is associated with low-angle grain boundaries(more specifically,the degree of misorientation<23.8°)and highangle grain boundaries with angles of misorientation greater than 48.8°-that is,themvalues follow the downward trend when the angle of GBs misorientation ranges between 6.5°and 48.8°.Several slip transfer parameters have values less than the maximummfactor.It can be clearly seen that slip transfer by basal to basal plane occurs when the grain boundary misorientation is smaller than 34.2°,and basal-pyramidal type takes place when the angle of misorientation is greater than 48.8°.It’s not clear why slip transfer is not observed when the geometrical relationships between incoming and outgoing slip systems have good coplanarity,but stress state adjacent to the boundary area certainly plays an important role in determining the slip transfer process.Furthermore,Fig.7 suggests that the threshold value of the slip transfer parameter,m,should larger than 0.77 in this alloy.

    Fig.7.(a)Plot of factor versus grain boundary misorientation for the alloy.Basal to basal slip is observed when the GB misorientation is smaller than 34.2°,and basal-pyramidal type occurs when the angle of misorientation is larger than 48.8°.

    3.3.Fracture mechanisms

    Heterogeneous deformation of polycrystals not only varies from grain to grain,but significant orientation gradients arise within each grain,as demonstrated in Fig.8b.The higher degree of misorientation is seen in the brighter yellow color.The kernel average misorientation(KAM)map,which is a local misorientation map,indicates that a higher level of misorientation and therefore strain develops around GBs compared to the interior of grains.High strain or high misorientation are generally seen at grain triple junctions during tension.Heterogeneous deformation within grains is expected to induce localized stress concentrations and promote local cracking.As indicated in Fig.8a,a crack propagates along the GB,in which grain 76 and 92 are oriented in hard and soft orientations respectively.In fact,small cracks tend to form at some points with both hard and soft orientation more often than regions that are both hard[1,33],where the extent and magnitude of dislocation motion in two grains give rise to the incompatible shape change.

    Fig.8.(a)Scanning electron micrograph showing crack propagation,in which a crack propagates between grain 76 and 92;(b)kernel average misorientation(KAM)map of grains analyzed in Fig.8a.The higher degree of misorientation and therefore strain is seen in the bright yellow.High strain or high misorientation are observed around the GBs and at grain triple junctions.

    4.Discussion

    Slip transfer in the rare earth magnesium alloy has been studied using anin-situtensile test.Basal to basal slip with the sametype Burgers vector is identified;whereas basal to pyramidal slip involvesand-type slip direction when the crystallographic misorientation is large.Previous work by Wang et al.[33]indicates that basal-prismatic slip transfer is mostly seen in an extruded Mg-10Gd-3Y-0.5Zr(wt.%)alloy with strong basal texture when the GBs have misorientation angles larger than 75°.In contrast to Wang’s work,observed slip transfer is most often basal-pyramidal type in this study when the angles of the GBs misorientation are greater than 48.8°.The maximum texture strength value of 5.6,relative to that from a random distribution in Fig.2b,indicates that the initial material has a weak texture.The difference in slip transfer type illustrates the significance of crystallographic orientation and texture on slip activity for these two alloys.

    Since the slip system with higher Schmid factor is most likely to activate and therefore enables slip transfer.The effect of the Schmid factor is evaluated in Fig.9,which plots the product ofmfactor and the sum of the Schmid factors versus the angles of the GB misorientation.Considering the different CRSS for basal and pyramidal slip,slip transfer types of basal-basal and basal-pyramidal are displayed separately in Fig.9a and b.For the data points representing slip transfer cases,values ofmtimes the sum of the two Schmid factors,SF,are expected to concentrate in the upper region of the figure under the assumption of uniaxial stress state.However,there is significant variation and scatter in the calculated results.The dashed line in Fig.9a indicates that a threshold for basal to basal slip may exist,below which instances of no slip transfer are most often observed and the slip transfer events occasionally arise.Basal to pyramidal slip has a wide range of Schmid factors on slip systems as shown in Fig.9b,and significant scatter in the calculated results ofmtimesSF indicates a complex threshold criteria for basal to pyramidal slip.The reason for the considerable scatter in the calculated data points is probably because the local stress tensors differ from uniaxial tension,which affects the Schmid factors.

    In agreement with previous studies,dislocations traversing GBs can take place in HCP polycrystals when pairs of slip systems have relatively high values ofmfactor.However,in spite of the fact that there is favorable geometrical continuity between slip systems from one grain into the other,instances of no slip transfer are increasingly observed when crystallographic misorientation larger than 6.3°,as shown in Fig.7.The slip transfer process is a function of the orientation of slip components upon which the stresses are taken to act.Coincident with the slip activity is the introduction of a complex stress state adjacent to the GBs.Crystal plasticity modeling provides that during deformation,the stress tensor within individual grain deviate from the uniaxial stress state[17],and this further affects the driving force for dislocation motion.The effect of the stress tensor,which can be predicted using CPFEM and/or measured by high angular resolution EBSD,needs to be examined in further work.

    As indicated in Fig.10a,sometimes the continuity of slip lines on the sample surface are clear(i.e.slip bands intersect at the common point across GB),but the presence of boundary ledge indicates heterogeneous strain on either side.Fig.10a demonstrates that more than one slip system operates in grain 204:basal and pyramidal slip system are both active in this grain.Although there is a highmfactor(Fig.10b),the magnitude of lattice distortion from the primary slip system in grain 204 causes incompatible deformation among two grains,such that slip transfer enabling a homogeneous strain is not significant in this case.The aligned slip traces with the presence of boundary ledge is also reported in a near-cube oriented Aluminum alloy[30].Additionally,cross slip process by which screw dislocations move from(-2112)plane to(11-22)plane along the<2-1-13>type direction is observed in grain 204,and thus wavy slip lines are observed on grains surface.Investigations by transmission electron microscopy(TEM)reveal that the addition of rare earth element dramatically changes the slip activity in magnesium alloys[34,35].Cross-slip is often reported in FCC metals;for magnesium alloy density functional theory and first-principle calculations suggest that the stacking fault energies(SFEs)decrease with the addition of rare earth element,and the stacking fault acts as nucleation source fordislocation[36-38].Therefore,pyramidal slip can operate and cross slip would be possible in an HCP crystal structure.

    Fig.9.The product of factor and the sum of the Schmid factors versus the angles of the GBs misorientation for(a)basal to basal slip and(b)basal to pyramidal slip respectively.The dashed line in Fig.9a suggests that the threshold criterion for basal to basal slip may exist.

    Fig.10.(a)It appears that basal to pyramidal slip takes place from grain 201 to 204;however,the boundary opens apart,indicating heterogeneous strain on either side.Meanwhile,more than one slip system operates in grain 204 and cross-slip occurs from pyramidal type plane to another.(b)mvalue of 0.78 for the observed basal-pyramidal slip pair.

    The strain state on the sample surface is analyzed based upon the crystallographic orientation data.As might be expected,poor alignment between slip systems across GBs may make dislocations traversing GBs more difficult during tension;rather the blocked dislocations introduce geometrically necessary dislocations(GNDs),which causes lattice distortion near GBs.Fig.11 gives an enlarged portion of strain maps on the deformed specimen surface,as indicated by the rectangle insets in Fig.2c.The local lattice curvature measured by a KAM map is used to examine the relationship with the angle of GBs misorientation.Small orientation changes adjacent to GBs region can be found in the Fig.11b,which means that there are some accumulations of GNDs along some boundaries.The grain reference orientation deviation(GROD)angle map in Fig.11c visualizes the long-range orientation gradients within the individual grains and demonstrates that nonuniform deformation develops within each grain.In addition,the GB misorientation map of Fig.11d indicates that there is no relationship between the density of GNDs and the high or low angle GBs.

    Fig.11.(a)SEM images showing the surface topography of the tensile specimen after straining;(b)KAM map,(c)GROD angle map,and(d)GB misorientation map of Fig.11a.

    The fact that slip transfer cases on the sample surface are highly aligned may suggest that due to the lack of constraint of adjacent grains the free surface facilitates the nucleation of dislocation at the GBs.Although only specimen surface is examined byin-situtension test,the maximum stress within structural parts occurs at surface positions for many common loading situations.Consequently,most of the time the failure of engineering materials is originating from the cracks on the component surface.Therefore,the deformation behavior on surface grains still has a significant influence on its failure mechanism[1].

    5.Conclusion

    A combination ofin-situtension test and EBSD orientation data is used to investigate the slip transfer process and heterogeneous deformation of a forged Mg-5.4Gd-1.8Y-1.5Zn alloy.The following conclusions can be drawn:

    1 Slip transfer process of basal-basal and basal-pyramidal type is observed for the material with weak texture.Basalbasal slip occurs when the GBs misorientation smaller is than 34.2°,whereas basal-pyramidal type takes place when the crystallographic misorientation is larger than 48.8°.

    2 As a non-sufficient but necessary condition for slip transfer cases,the Luster-Morris factor,m,should larger than 0.77 for this alloy.

    3 Strain accommodation of individual grains by means of slip occurred more easily than slip transfer when several slip systems were operable in this alloy.

    4 The product ofmfactor and the sum of the Schmid factors demonstrates that the threshold for basal-basal slip may exist,however,there is significant scatter in the calculated results ofmtimesSF for basal to pyramidal slip,indicating complex stress tensor within the grains to initiate pyramidal slip.

    5 Besides the geometrical criterion,the stress state in each grain plays an important role in determining the slip transfer process.The effect of the stress tensor,which can be predicted using crystal plasticity modeling,needs to be examined in further work.

    Acknowledgements

    The authors acknowledge the funding support from National Natural Science Foundation of China via grant 51775135.

    Data availability

    The raw/processed data required to reproduce these findings cannot be shared at this time as the data also forms part of an ongoing study.

    能在线免费看毛片的网站| 亚洲国产欧美在线一区| 黄片wwwwww| 色吧在线观看| 我的女老师完整版在线观看| 欧美精品国产亚洲| 狂野欧美激情性xxxx在线观看| 日本猛色少妇xxxxx猛交久久| 色哟哟·www| 国国产精品蜜臀av免费| 亚洲乱码一区二区免费版| 午夜福利高清视频| 成人鲁丝片一二三区免费| 简卡轻食公司| 欧美xxxx黑人xx丫x性爽| 男女国产视频网站| 3wmmmm亚洲av在线观看| 久久久久久久久久久免费av| 99久久人妻综合| 国产精品美女特级片免费视频播放器| 免费看a级黄色片| 成人亚洲欧美一区二区av| 日本猛色少妇xxxxx猛交久久| 99视频精品全部免费 在线| 国产精品女同一区二区软件| 成年女人看的毛片在线观看| 亚洲内射少妇av| 亚洲欧美日韩东京热| 久久久久九九精品影院| 国产一区二区在线观看日韩| 欧美xxxx性猛交bbbb| 久久久久国产网址| 亚洲欧美精品专区久久| 岛国在线免费视频观看| 国产高清三级在线| 伊人久久精品亚洲午夜| 99在线视频只有这里精品首页| 噜噜噜噜噜久久久久久91| 亚洲国产欧洲综合997久久,| 26uuu在线亚洲综合色| 最新中文字幕久久久久| 亚洲成人精品中文字幕电影| 国产精品.久久久| 黄色一级大片看看| 中文亚洲av片在线观看爽| 国产精品国产高清国产av| 看黄色毛片网站| 亚洲无线观看免费| 亚洲,欧美,日韩| 亚洲成人av在线免费| 亚洲精品日韩在线中文字幕| 男女啪啪激烈高潮av片| 亚洲欧美精品自产自拍| 日韩欧美 国产精品| 久久久午夜欧美精品| 日本与韩国留学比较| 国产精品三级大全| 韩国av在线不卡| 天堂av国产一区二区熟女人妻| 欧美三级亚洲精品| 国产国拍精品亚洲av在线观看| 一区二区三区四区激情视频| 午夜福利视频1000在线观看| 国内精品宾馆在线| 亚洲人与动物交配视频| 国产亚洲av嫩草精品影院| 日本欧美国产在线视频| 国产美女午夜福利| 少妇熟女欧美另类| 欧美3d第一页| 特大巨黑吊av在线直播| 精品人妻一区二区三区麻豆| 欧美潮喷喷水| 久久精品熟女亚洲av麻豆精品 | 中文字幕av在线有码专区| 99热这里只有是精品在线观看| 精品久久久久久久末码| av又黄又爽大尺度在线免费看 | 国产一区二区在线观看日韩| 91久久精品电影网| 欧美不卡视频在线免费观看| 如何舔出高潮| 亚洲av免费高清在线观看| 舔av片在线| 欧美bdsm另类| 久久热精品热| 波野结衣二区三区在线| 国产精品人妻久久久久久| 级片在线观看| 久久精品久久久久久久性| 毛片女人毛片| 国模一区二区三区四区视频| 少妇的逼好多水| 亚洲av福利一区| 舔av片在线| 久久国内精品自在自线图片| 18禁动态无遮挡网站| 天堂网av新在线| 亚洲国产精品合色在线| 三级国产精品欧美在线观看| 少妇丰满av| 精品久久国产蜜桃| 91久久精品国产一区二区三区| 综合色av麻豆| 插阴视频在线观看视频| 建设人人有责人人尽责人人享有的 | 国产午夜精品一二区理论片| 国内精品一区二区在线观看| 国产亚洲精品av在线| .国产精品久久| 国产精品福利在线免费观看| 男人舔女人下体高潮全视频| 久久久精品94久久精品| 一本久久精品| 特大巨黑吊av在线直播| 欧美变态另类bdsm刘玥| 国产av码专区亚洲av| 亚洲自拍偷在线| 99久久精品热视频| 午夜免费激情av| 赤兔流量卡办理| 免费无遮挡裸体视频| 中文字幕精品亚洲无线码一区| 国模一区二区三区四区视频| 久久精品影院6| 国产真实伦视频高清在线观看| 欧美日韩一区二区视频在线观看视频在线 | 波多野结衣高清无吗| eeuss影院久久| 最近视频中文字幕2019在线8| 国内精品宾馆在线| videos熟女内射| 天天一区二区日本电影三级| 婷婷色麻豆天堂久久 | 精品久久久久久成人av| 搡女人真爽免费视频火全软件| 亚洲精品亚洲一区二区| 国产精品福利在线免费观看| 26uuu在线亚洲综合色| 日韩制服骚丝袜av| 久久婷婷人人爽人人干人人爱| 久久久亚洲精品成人影院| 日本-黄色视频高清免费观看| 少妇高潮的动态图| 久久亚洲国产成人精品v| 中文字幕精品亚洲无线码一区| 午夜激情欧美在线| 亚洲一区高清亚洲精品| 久久精品国产亚洲av天美| 成人鲁丝片一二三区免费| 91精品国产九色| 别揉我奶头 嗯啊视频| 精品久久久噜噜| 别揉我奶头 嗯啊视频| 国产精品国产三级国产专区5o | 99热这里只有精品一区| 亚洲一区高清亚洲精品| 国产麻豆成人av免费视频| 男人和女人高潮做爰伦理| 97热精品久久久久久| 欧美激情在线99| 欧美一区二区亚洲| 亚洲av中文av极速乱| 免费一级毛片在线播放高清视频| 日韩高清综合在线| 亚洲成人精品中文字幕电影| 国产色爽女视频免费观看| 亚洲内射少妇av| 永久网站在线| 国产乱人偷精品视频| av在线天堂中文字幕| 男的添女的下面高潮视频| 国产精品久久电影中文字幕| 欧美性猛交黑人性爽| 欧美日韩综合久久久久久| 亚洲欧美一区二区三区国产| 我的女老师完整版在线观看| 精品酒店卫生间| 欧美丝袜亚洲另类| 男插女下体视频免费在线播放| 超碰av人人做人人爽久久| 国产精品日韩av在线免费观看| 日本wwww免费看| 亚洲自拍偷在线| 自拍偷自拍亚洲精品老妇| 天堂中文最新版在线下载 | 亚洲av熟女| 国产精品爽爽va在线观看网站| 亚洲av免费在线观看| 中文欧美无线码| 22中文网久久字幕| 在线a可以看的网站| 国产探花在线观看一区二区| 亚洲av免费高清在线观看| 国产精品精品国产色婷婷| 爱豆传媒免费全集在线观看| 国产精品国产三级专区第一集| 亚洲精品久久久久久婷婷小说 | 亚洲精品影视一区二区三区av| 中文精品一卡2卡3卡4更新| 日本wwww免费看| 在线观看一区二区三区| 91av网一区二区| 婷婷色麻豆天堂久久 | 精品午夜福利在线看| 99视频精品全部免费 在线| 国产精品乱码一区二三区的特点| 久久人妻av系列| 极品教师在线视频| 国产精品一区www在线观看| 亚洲国产最新在线播放| 久久99热这里只频精品6学生 | 国国产精品蜜臀av免费| 黑人高潮一二区| 我的女老师完整版在线观看| 亚洲伊人久久精品综合 | av又黄又爽大尺度在线免费看 | 我要看日韩黄色一级片| 久久久精品欧美日韩精品| 国产一区二区在线av高清观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产黄色小视频在线观看| 熟女人妻精品中文字幕| 乱码一卡2卡4卡精品| 99久久九九国产精品国产免费| 波多野结衣高清无吗| 亚洲欧洲日产国产| 国内精品宾馆在线| 中文在线观看免费www的网站| 亚洲人成网站在线观看播放| 黄色日韩在线| 久久久午夜欧美精品| 国产私拍福利视频在线观看| 狂野欧美激情性xxxx在线观看| 免费观看的影片在线观看| 亚洲精品国产av成人精品| 99久国产av精品国产电影| 亚洲欧美中文字幕日韩二区| 午夜激情福利司机影院| 国产精品久久久久久精品电影| 男人舔奶头视频| 2022亚洲国产成人精品| 只有这里有精品99| 国产午夜精品久久久久久一区二区三区| 村上凉子中文字幕在线| www.色视频.com| 精品久久国产蜜桃| 美女内射精品一级片tv| 成人性生交大片免费视频hd| 成人鲁丝片一二三区免费| 亚洲精品乱码久久久v下载方式| 国产免费又黄又爽又色| 国产亚洲精品av在线| 熟妇人妻久久中文字幕3abv| 精品不卡国产一区二区三区| ponron亚洲| 深夜a级毛片| 七月丁香在线播放| 男女国产视频网站| 一边摸一边抽搐一进一小说| 国产精品国产三级国产专区5o | ponron亚洲| 久久亚洲精品不卡| 国产国拍精品亚洲av在线观看| 久久精品久久精品一区二区三区| 亚洲综合精品二区| 国产69精品久久久久777片| 国产av一区在线观看免费| 国产精品不卡视频一区二区| 欧美日韩一区二区视频在线观看视频在线 | 亚洲在线自拍视频| 97超碰精品成人国产| 99久久九九国产精品国产免费| 纵有疾风起免费观看全集完整版 | 啦啦啦韩国在线观看视频| 女人久久www免费人成看片 | 一本久久精品| 女人被狂操c到高潮| 级片在线观看| 看片在线看免费视频| 国产精品人妻久久久久久| 内射极品少妇av片p| 波多野结衣巨乳人妻| 日韩国内少妇激情av| 日韩av不卡免费在线播放| 九九久久精品国产亚洲av麻豆| 欧美zozozo另类| 成人亚洲精品av一区二区| 亚洲久久久久久中文字幕| 白带黄色成豆腐渣| 久久99热这里只有精品18| 禁无遮挡网站| 麻豆一二三区av精品| 91aial.com中文字幕在线观看| 舔av片在线| 老司机福利观看| 伊人久久精品亚洲午夜| 久久国产乱子免费精品| 毛片一级片免费看久久久久| 成人二区视频| 在线播放无遮挡| 哪个播放器可以免费观看大片| 欧美激情在线99| 在线播放国产精品三级| 成人高潮视频无遮挡免费网站| 人人妻人人澡欧美一区二区| 国产精品爽爽va在线观看网站| 欧美色视频一区免费| 狂野欧美白嫩少妇大欣赏| 亚洲美女搞黄在线观看| 久久久欧美国产精品| 大话2 男鬼变身卡| 国产私拍福利视频在线观看| 国产免费视频播放在线视频 | 午夜亚洲福利在线播放| 秋霞在线观看毛片| 国产淫语在线视频| 99视频精品全部免费 在线| 看非洲黑人一级黄片| 欧美性感艳星| 午夜福利视频1000在线观看| 久久精品国产鲁丝片午夜精品| 高清视频免费观看一区二区 | 搞女人的毛片| 欧美激情久久久久久爽电影| 日韩大片免费观看网站 | 国产精品美女特级片免费视频播放器| 成人亚洲精品av一区二区| 两性午夜刺激爽爽歪歪视频在线观看| 午夜福利在线观看免费完整高清在| 精品人妻视频免费看| 免费av毛片视频| 日韩av在线大香蕉| 特大巨黑吊av在线直播| 亚洲性久久影院| 波多野结衣高清无吗| 色哟哟·www| 日韩亚洲欧美综合| 建设人人有责人人尽责人人享有的 | 免费观看的影片在线观看| 免费黄色在线免费观看| 最近手机中文字幕大全| 国产乱来视频区| 中文字幕人妻熟人妻熟丝袜美| 久久久欧美国产精品| 亚洲中文字幕日韩| 日韩大片免费观看网站 | 亚洲国产成人一精品久久久| 免费av观看视频| 免费观看在线日韩| 综合色av麻豆| 美女高潮的动态| 欧美日本视频| 中文字幕av成人在线电影| 成年女人永久免费观看视频| 久久韩国三级中文字幕| 非洲黑人性xxxx精品又粗又长| 国产欧美日韩精品一区二区| 国产免费一级a男人的天堂| 久久精品影院6| 1024手机看黄色片| 国产精品电影一区二区三区| 国产爱豆传媒在线观看| 男女边吃奶边做爰视频| 精品无人区乱码1区二区| 91狼人影院| 日韩欧美精品v在线| 我要看日韩黄色一级片| 成人鲁丝片一二三区免费| 三级男女做爰猛烈吃奶摸视频| 黑人高潮一二区| 成年免费大片在线观看| 中国美白少妇内射xxxbb| 亚洲自拍偷在线| 日日摸夜夜添夜夜爱| 高清毛片免费看| 插逼视频在线观看| 国产高清视频在线观看网站| 熟女人妻精品中文字幕| 日韩在线高清观看一区二区三区| 亚洲乱码一区二区免费版| 两个人视频免费观看高清| 99久久人妻综合| 午夜福利在线观看免费完整高清在| 又粗又硬又长又爽又黄的视频| 国产精品美女特级片免费视频播放器| 国产欧美日韩精品一区二区| 亚洲av中文av极速乱| 纵有疾风起免费观看全集完整版 | 级片在线观看| h日本视频在线播放| 日日撸夜夜添| 最近手机中文字幕大全| 黄色欧美视频在线观看| 91久久精品电影网| 噜噜噜噜噜久久久久久91| 久热久热在线精品观看| 亚洲最大成人中文| 能在线免费观看的黄片| 国产探花在线观看一区二区| 久久久亚洲精品成人影院| 中文字幕熟女人妻在线| 又爽又黄无遮挡网站| 九草在线视频观看| 色综合亚洲欧美另类图片| 日本一二三区视频观看| 国产伦精品一区二区三区四那| 欧美又色又爽又黄视频| 日韩成人av中文字幕在线观看| 中文字幕亚洲精品专区| h日本视频在线播放| 97超视频在线观看视频| 国产精品久久久久久精品电影小说 | 美女国产视频在线观看| 亚洲国产最新在线播放| 国产欧美日韩精品一区二区| 黄片无遮挡物在线观看| 久久99热这里只频精品6学生 | 天堂av国产一区二区熟女人妻| 日产精品乱码卡一卡2卡三| 91在线精品国自产拍蜜月| 精品人妻一区二区三区麻豆| 欧美精品一区二区大全| 国产精品日韩av在线免费观看| 亚洲欧美中文字幕日韩二区| h日本视频在线播放| 舔av片在线| 亚洲美女搞黄在线观看| 国产精品.久久久| 国产私拍福利视频在线观看| 好男人在线观看高清免费视频| 国产伦精品一区二区三区四那| 日本免费一区二区三区高清不卡| 成人国产麻豆网| or卡值多少钱| 国产精华一区二区三区| 久热久热在线精品观看| 国产精品电影一区二区三区| 国产精品.久久久| 久久欧美精品欧美久久欧美| 美女国产视频在线观看| 欧美性感艳星| 中文乱码字字幕精品一区二区三区 | 日韩视频在线欧美| 看十八女毛片水多多多| 欧美日本亚洲视频在线播放| 婷婷色av中文字幕| 人妻制服诱惑在线中文字幕| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品久久久久久精品电影小说 | 青春草国产在线视频| 免费观看在线日韩| 国产精品久久视频播放| 熟妇人妻久久中文字幕3abv| 欧美性猛交黑人性爽| 国产在视频线在精品| 国产精品无大码| 午夜激情福利司机影院| 亚洲欧美精品自产自拍| 好男人视频免费观看在线| 日韩欧美国产在线观看| 国产成人精品久久久久久| 色综合站精品国产| 99久国产av精品| 精品一区二区三区人妻视频| 狂野欧美激情性xxxx在线观看| 亚洲av不卡在线观看| 亚洲av日韩在线播放| 午夜亚洲福利在线播放| 老女人水多毛片| 91午夜精品亚洲一区二区三区| 天美传媒精品一区二区| 日韩欧美国产在线观看| 久久欧美精品欧美久久欧美| 欧美激情在线99| 国产淫语在线视频| 五月伊人婷婷丁香| 久久国内精品自在自线图片| 亚洲欧美成人精品一区二区| 精华霜和精华液先用哪个| 好男人在线观看高清免费视频| 伦精品一区二区三区| 久久久久久久久久黄片| 久久久久性生活片| 女人久久www免费人成看片 | 久久精品国产99精品国产亚洲性色| 亚洲经典国产精华液单| 久久久久久久午夜电影| 国产高清三级在线| 成人三级黄色视频| 成人漫画全彩无遮挡| 午夜激情欧美在线| 精品国产露脸久久av麻豆 | 欧美日韩在线观看h| 国产麻豆成人av免费视频| 中文天堂在线官网| 熟女人妻精品中文字幕| 男的添女的下面高潮视频| 99在线视频只有这里精品首页| 大话2 男鬼变身卡| 国产精品熟女久久久久浪| 两个人的视频大全免费| 国产精品一区www在线观看| 日本黄大片高清| 日本欧美国产在线视频| av在线蜜桃| 成人性生交大片免费视频hd| 午夜福利高清视频| 久久精品国产亚洲网站| 亚洲精品,欧美精品| 啦啦啦啦在线视频资源| 久久国产乱子免费精品| 欧美成人免费av一区二区三区| 国产精华一区二区三区| 国产精品一区二区性色av| 亚洲欧美中文字幕日韩二区| 在线观看av片永久免费下载| 能在线免费观看的黄片| 日韩精品青青久久久久久| 欧美xxxx性猛交bbbb| 2022亚洲国产成人精品| 亚洲一区高清亚洲精品| 哪个播放器可以免费观看大片| 特大巨黑吊av在线直播| 久久亚洲精品不卡| 国产av一区在线观看免费| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 一区二区三区高清视频在线| 国产精品一区www在线观看| 日韩精品青青久久久久久| 桃色一区二区三区在线观看| 久久久午夜欧美精品| 亚洲一区高清亚洲精品| 少妇被粗大猛烈的视频| 精品久久久久久成人av| 国产精品久久久久久精品电影| 国产大屁股一区二区在线视频| 亚洲精品色激情综合| 成人综合一区亚洲| 成人特级av手机在线观看| av在线亚洲专区| 成人午夜高清在线视频| 99久久九九国产精品国产免费| 亚洲婷婷狠狠爱综合网| 国产精品久久久久久久电影| 1024手机看黄色片| 亚洲精品亚洲一区二区| 国产69精品久久久久777片| 少妇丰满av| 直男gayav资源| 欧美精品国产亚洲| 人人妻人人看人人澡| 最近视频中文字幕2019在线8| 色5月婷婷丁香| 女人十人毛片免费观看3o分钟| 免费av不卡在线播放| 1024手机看黄色片| 久久久午夜欧美精品| 91精品国产九色| 九九久久精品国产亚洲av麻豆| 国产精品久久久久久久久免| 搡女人真爽免费视频火全软件| 国产大屁股一区二区在线视频| 爱豆传媒免费全集在线观看| 国产色婷婷99| 久久久久久大精品| 麻豆成人av视频| 久久精品久久精品一区二区三区| 久久久久免费精品人妻一区二区| 日本爱情动作片www.在线观看| 99久久九九国产精品国产免费| 六月丁香七月| 国内精品美女久久久久久| 少妇人妻精品综合一区二区| 国产精品国产三级专区第一集| 亚洲欧美精品综合久久99| 精品久久久久久电影网 | 在线免费观看的www视频| 又黄又爽又刺激的免费视频.| 亚洲内射少妇av| 久久亚洲国产成人精品v| 久久精品久久精品一区二区三区| 老司机影院毛片| 国产老妇伦熟女老妇高清| 波多野结衣高清无吗| 91精品国产九色| 菩萨蛮人人尽说江南好唐韦庄 | 日本黄色视频三级网站网址| 男人和女人高潮做爰伦理| 精品一区二区三区人妻视频| 久久久国产成人免费| 亚洲精品乱码久久久久久按摩| 久久久色成人| 伦理电影大哥的女人| videossex国产| 寂寞人妻少妇视频99o| 日韩一区二区三区影片| 少妇人妻一区二区三区视频| 中文字幕亚洲精品专区| 亚州av有码| 国产日韩欧美在线精品| 亚洲熟妇中文字幕五十中出| 日本三级黄在线观看| 国内揄拍国产精品人妻在线| 国产一区有黄有色的免费视频 | 蜜臀久久99精品久久宅男| 亚洲精品456在线播放app| 99久久精品热视频| 国产成人精品一,二区| 噜噜噜噜噜久久久久久91| 欧美激情久久久久久爽电影| 欧美xxxx黑人xx丫x性爽| 纵有疾风起免费观看全集完整版 | 欧美性感艳星| 国产成人freesex在线|