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

    Effects of orientation on the fatigue crack growth behaviors of the ZK60 magnesium alloy in air and PBS

    2024-04-18 13:44:34JiqiHuZhengLiuZuolingNingHongGo
    Journal of Magnesium and Alloys 2024年1期

    Jiqi Hu ,Zheng Liu ,Zuoling Ning ,Hong Go,b,?

    a School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, China

    b Tianjin Key Laboratory of Chemical Process Safety and Equipment Technology, Tianjin 300072, China

    Abstract Strong anisotropic corrosion and mechanical properties caused by specimen orientations greatly limit the applications of wrought magnesium alloys.To investigate the influences of specimen orientation,the corrosion tests and (corrosion) fatigue crack growth tests were conducted.The rolled and transverse surfaces of the materials show distinct corrosion rate differences in the stable corrosion stage,but the truth is the opposite for the initial stage of corrosion.In air,specimen orientations have a significant influence on the plastic deformation mechanisms near the crack tip,which results in different fatigue fracture surfaces and cracking paths.Compared with R-T specimens,N-T specimens show a slower fatigue crack growth (FCG) rate in air,which can be attributed to crack closure effects and deformation twinning near the crack tip.The corrosion environment will not significantly change the main plastic deformation mechanisms for the same type of specimen.However,the FCG rate in phosphate buffer saline (PBS) is one order of magnitude higher than that in air,which is caused by the combined effects of hydrogen-induced cracking and anodic dissolution.Owing to the similar corrosion rates at crack tips,the specimens with different orientations display close FCG rates in PBS.

    Keywords: Magnesium alloy;Orientation;Corrosion;Fatigue crack growth.

    1.Introduction

    Due to its excellent physical properties such as low density,good machinability,and high specific strength,magnesium(Mg) and its alloys have received increasing attention in numerous fields,including aerospace,automobile,and biomedical [1–3].Especially for biological implantations,Mg alloys possess excellent bio-degradability,comparable mechanical properties with natural bone,and the best bio-compatibilities with human physiology.Over the past two decades,many studies were conducted to explore,understand,and improve their mechanical and corrosion properties.

    Mg and its alloys have a hexagonal close-packed (HCP)crystal structure,and their mechanical properties are different from those of more commonly used cubic materials due to fewer and more difficult to activate slip systems[4].The basal slip (i.e.(0001)<110>) and two non-basal slip modes (the prismatic ({100}<110>) and the pyramidal({101}<110>)) provide deformation parallel to the basal plane only.Furthermore,the pyramidal slip is difficult to activate at room temperature due to the high critical resolved shear stress (CRSS).Therefore,Mg and its alloys must rely on deformation twinning to accommodate c-axis deformation,then accomplish arbitrary homogenous deformation.The most commonly observed deformation twinning mode is the{102}extension twinning,which occurs under c-axis tension[5].The interaction between the above-mentioned deformation mechanisms will lead to a strong mechanical anisotropy in both single crystal and textured polycrystalline Mg.In addition,the texture orientation affects the corrosion behavior of Mg alloys significantly.The basal plane of grains has better corrosion resistance than the non-basal plane due to its higher atomic density and lower surface energy.Compared with the surface with random grain orientations,the surface composed of crystal planes with similar orientations displays higher corrosion resistance due to the suppression of localized micro-galvanic corrosion [6].

    Specimen orientation has a significant influence on deformation mechanisms [7],fatigue behaviors [8],and crack propagation performance [9].Owing to the combined effect of Schmid factor (SF) and CRSS,0° (rolled direction) ~90°(normal direction) orientation specimens display the different yield stresses and deformation modes under tensile tests[10].In air,0° specimens show the highest fatigue strength while the 45° specimens display the lowest fatigue strength and hysteresis with the best symmetry [8].In the Mode-I fatigue crack growth experiments,the orientation of extruded Mg alloy also has a significant influence on the crack growth rate and path [11].As an implant for the human body,Mg alloys not only bear the cyclic load required by body movement but also suffer from corrosion from the human environment.Under the coupling of stress and corrosion,Mg alloys have proved to be susceptible to environmentally assisted cracking in corrosion environments,displaying significantly lower corrosion fatigue resistance in the modified simulated body fluids than in air [12].Solution treatment can increase the corrosion fatigue resistance of as-forged Mg-Zn-Y-Zr alloy by suppressing localized corrosion sites [13].Ultrafine-grains obtained by carefully designed multi-directional compression treatment greatly improve the strength and corrosion resistance of the AZ91 Mg alloy [14].The change of solution pH values also has an obvious influence on the corrosion and crack growth behaviors of the ZK60 Mg alloy [15].Few studies have discussed the effect of orientation on mechanical and chemical coupling failure behaviors of ZK60 Mg alloys.Therefore,revealing the interaction between corrosion and cyclic stress of specimens with different orientations is crucial for Mg alloys’application.

    In the current study,the corrosion and electrochemical corrosion experiments of two oriented surfaces were carried out to analyze their corrosion behavior and provide evidence for the following discussion.The failure behaviors of different specimen orientations were explored through crack growth experiments in air and corrosion environment.Meanwhile,their deformation and failure mechanisms were analyzed by microstructure characterization means.

    2.Materials and methods

    Commercial rolled ZK60 plates of 40 mm thickness are used in the current study,the element composition is listed in Table 1.To distinguish the orientations of rolled plates,the marks of three directions (Rolled Direction (RD),Transverse Direction (TD),Normal Direction (ND)) and three surfaces (Rolled Surface (RS),Transverse Surface (TS),Longitudinal Surface (LS)) of plates are arranged and shown in Fig.1.The metallography of different surfaces was observed by KEYENCE (VHX-900F) optical microscope (OM).The texture of the material was inspected by the FEI (Apreo S LoVac) scanning electron microscope (SEM) equipped with an EDAX (Hikari Super) electron backscattering diffraction(EBSD) detector.

    Table 1 Chemical element composition of ZK60 Mg alloys.

    Fig.1.Schematic diagram of the orientation marks of rolled plate.

    The hydrogen evolution tests were carried out on the RS and TS,and other surfaces were protected with epoxy.Before the tests,the measured surface was also grounded with 1000-grit green silicon carbide papers.Then,the specimens were horizontally immersed for 50 h in 37 °C PBS solution (8 g/L NaCl,0.2 g/L KCl,1.44 g/L Na2HPO4,and 0.24 g/L KH2PO4) with saturated hydrogen,and the evolved hydrogen was collected in a burette above the corroded specimens and regularly recorded.Meanwhile,electrochemical tests were also conducted on the RS and TS using a threeelectrode system with a saturated calomel electrode (SCE)acting as the reference electrode,a platinum electrode acting as the counter electrode,and the specimen acting as the working electrode.Electrochemical impedance spectroscopy(EIS) tests were conducted at different immersion times (2,4,6,8,10,and 12 h).The frequency range and amplitude of EIS were 100 kHz–0.01 Hz and 10 mV,respectively.Following the last EIS test,potentiodynamic polarization tests were conducted with a scan rate of 2 mV/s from -0.5 VSCEto+0.5 VSCErelative to the open circuit potential.

    The mechanical experiments were carried out on IBTC-4000in-situmechanical test system.The dog-bone-shaped specimen was prepared for monotonic mechanical tests,which were controlled by machine displacement with a speed of 0.03 mm/s at room temperature,specimens’ dimensions were shown in Fig.2(a).Compact tension(CT)specimens with two orientations were prepared for fatigue crack growth (FCG)experiments,specimens’ dimensions and orientations were shown in Fig.2(b) and 2(c),respectively.The "X-Y" was used to mark specimens with different orientations,"X" and"Y" represented the normal direction of the notch plane and the direction of crack propagation,respectively.Before the test,the surfaces of the CT specimen were protected by grease,and only the fracture surfaces were exposed to PBS.A triangular wave with a peak value of 235 N,a stress ratio (minimum load over maximum load in a loading cycle) of 0.1,and a frequency of 0.5 Hz was used as the testing load.For the FCG tests in PBS,an online system shown in Fig.3 was built.The CT specimen was vertically set in a corrosion chamber containing PBS,which maintained a constant temperature (37 °C) and a constant flow rate (160 mL/min).The crack length was observed by the camera and recorded by counting screw with an accuracy of 10 μm.

    Fig.2.(a) Dog-bone-shaped specimen,(b) CT specimen,and (c) the orientation of CT specimens.

    Fig.3.Schematic diagram of online fatigue crack growth system.

    The corrosion product layer and corrosion surface of the hydrogen evolution specimens were observed by SEM.After FCG tests,the fracture surfaces and crack growth paths of the specimens were inspected using OM and SEM.In addition,the microstructure near the crack tip in air was observed by SEM and EBSD,before the FCG tests,the CT specimen’s surface was electrochemical polished using AC2 solution.

    3.Results

    3.1. Metallography, texture, and tensile properties

    Fig.4 shows the three-dimensional metallography and texture of rolled plate.According to the metallography in Fig.4(a),the RS is composed of irregular large-grains,equiaxed small-grains,and cluster fine-grains.The TS and LS have a similar grain structure,which is composed of long strip large-grains,equiaxed small-grains,and cluster fine-grains.According to the results of EBSD shown in Fig.4(b,c),the rolled plate has a strong basal texture,that is,the c-axis of most grains are parallel to the ND.

    Fig.4.(a) Three-dimensional metallography of rolled plate,(b) inverse pole figure of the TS,and (c) pole figure of {0001} and {110} crystal planes.

    The uniaxial tensile curves of ZK60 plates in RD,TD,and ND are shown in Fig.5.The mechanical property parameters in different directions are listed in Table 2.Under the tension loading,the {102} tension twin and basal slip are in hard orientation for RD and TD,but ND specimen can activate extension twins easily,thereby showing low yield stress (the stress corresponding to 0.2% plastic strain).The existence of twins hinders the prismatic slip system in the subsequent stage of tension,resulting in greater engineering stress in ND than that in RD and TD.The ND specimen also shows better elongation and ductility due to more available deformation systems.According to the {110} pole figure shown in Fig.4(c),the distribution of grain a-axis is random,resulting in similar mechanical properties of TD and RD.

    Table 2 Uniaxial mechanical property parameters in different directions of ZK60 plates.

    Fig.5.Uniaxial tensile curves of ZK60 plate along RD,TD,and ND.

    3.2. Corrosion and electrochemical corrosion

    Fig.6 shows the hydrogen evolution curves of RS and TS.According to the enlarged figure,the difference in hydrogen evolution volume between the two orientations is small within 4 h.After that,TS specimens directly enter the stable corrosion stage,while RS specimens take a longer time to stabilize.During the stable corrosion stage,the hydrogen evolution rate of the TS specimen is twice as fast as that of the RS specimen.After hydrogen evolution tests,the corrosion morphology of RS and TS is observed.As shown in Fig.7(a,b),the corrosion layers on two types of specimens both consist of cracked block products and fine powdery products.The block products are Mg(OH)2formed by the reaction of magnesium matrix and aqueous solution,and the powdery products are poorly soluble phosphate Mg3(PO4)2formed by magnesium ions and phosphate ions in PBS [15].After immersion for 2 h,RS (Fig.7(c)) and TS (Fig.7(d)) without corrosion products both show typical filiform corrosion characteristics.Due to the different metallographic structures shown in Fig.4(a),their filiform corrosion traces display different geometry.In addition,a few small pitting pits are observed on the path of filiform corrosion.After immersion for 50 h,RS is mainly composed of fine and shallow pitting corrosion pits shown in Fig.7(e),while TS has a lot of large pitting pits shown in Fig.7(f).According to the three-dimensional morphology results shown in Fig.7(g,h),the surface corrosion of the TS specimen is more serious,which is consistent with the results of hydrogen evolution tests.

    Fig.7.The corrosion product layers ((a),(b)),corrosion surfaces for 2 h ((c),(d)) and 50 h ((e) (f)),three-dimensional surface morphology ((g),(h)) after immersion tests.

    Fig.8 shows the electrochemical corrosion results of RS and TS.The Nyquist diagrams of different immersion times shown in Fig.8(a) are fitted using the equivalent circuit in Fig.8(b),and fitting results are listed in Table 3.In the equivalent circuit,Rsrepresents the solution resistance.Doublelayer capacitanceCtin parallel with charge transfer resistanceRtis used to simulate the capacitive reactance arc in the middle-frequency region.Rtis related to the interface charge transferability when the alloys dissolve,the larger value means the slower dissolution rate of the anode i.e.,Mg matrix.The inductive arc in the low-frequency region is simulated by equivalent resistanceRLand equivalent inductanceL[16].The occurrence of inductive arc is concerned with the formation of pitting nucleation,and the inductive response mainly comes from the change of film thickness and adsorption film coverage at pitting corrosion.Fig.8(c,d) summarize the varieties ofRtandRLwith immersion time.The trends ofRtare in good agreement with the results of hydrogen evolution tests shown in Fig.6.In terms of the equivalent pitting resistanceRL,the TS specimen is always lower than the RS specimen,that is,the pitting tendency of the TS specimen is greater.Coupled with the higher anodic dissolution rate of the TS specimen,larger pitting corrosion pits shown in Fig.7(f) are formed on its surface.Fig.8(e) shows the potentiodynamic polarization results of RS and TS specimens after 12 h’ immersion.The Tafel curves of both show a similar shape,which indicates that they have similar electrochemical corrosion behavior.According to the fitting results of the Tafel curves in Table 4,the corrosion potentialVcorrof RS and TS is almost equal.The results of corrosion current densityicorrare in good agreement with the results of hydrogen evolution,in which the corrosion rate of RS is about half of that of TS during the stable corrosion stage.

    Table 3 EIS fitting results of different immersion times for RS and TS specimens.

    Table 4 The fitting results of Tafel curves.

    Fig.8.(a) Nyquist diagram of different immersion time,(b) equivalent circuit,(c) varieties of charge transfer resistance Rt and (d) inductance equivalent resistance RL and with immersion time,and (e) Tafel curves for RS and TS specimens.

    3.3. Fatigue crack growth rate

    Fig.9 shows the FCG rate curves of CT specimens with different orientations in air and PBS.The crack growth rate(da/dN) can be obtained from Eq.(1):

    Fig.9.FCG rate of CT specimens with different orientations in air and PBS.

    whereΔPis the difference between maximum load and minimum load,Bis the thickness of specimens,Wis the width of specimens andα=/W.In the logarithmic coordinate system,da/dNandΔKshow a good linear relationship,which can be described by the following Paris formula:

    whereCandmare material parameters.

    As shown in Fig.9,the corrosion FCG rate of specimens is about one order of magnitude higher than the atmospheric FCG rate.In air,the FCG rate of R-T specimens is higher than that of N-T specimens.In PBS,the difference in FCG rate between them is not obvious.In the current study,no obvious three sub-stages reported in Ref.[9],which is characterized by three distinguishable Paris law slopes during the stable crack growth,are found.Only a single slope commonly occurring in traditional materials is observed,which may be concerned with the relatively continuous crack growth path and limitedΔKrange.

    3.4. Fatigue fracture observation

    Fig.10 shows the overall and partial enlarged detail of fracture for the FCG tests in air and PBS.As shown in Fig.10(a–d),for the same specimen orientation,the fracture surface in air is flatter than that in PBS.Under the same environment,the fracture surface of R-T specimens is significantly smoother than that of N-T specimens.From the partial enlarged detail in Fig.10(e),it can be found that the fracture surface of R-T specimens in air shows typical quasicleavage fracture characteristics,which are composed of many river-like cleavage planes and tear ridges along the crack’s growth direction.In contrast,the tear ridges of N-T specimens shown in Fig.10(f) are disordered,and larger cleavage planes and V-shaped patterns with an included angle of about 122° can be found at local positions.These characteristics indicate that the N-T specimens frequently pass through the long strip large-grains,and similar phenomena have occurred in the Ref.[11].Under the influence of chemical corrosion,the cleavage planes and tear ridges of fracture in PBS become not obvious,which is shown in Fig.10(g,h).In addition,many secondary cracks appear on the corrosion fatigue fracture surface,which is usually related to hydrogen embrittlement [17].Many small-sized secondary cracks along the propagation direction are found on the fracture surface of R-T specimens.While the fracture surface of the N-T specimen appears many large-sized secondary cracks with irregular distribution directions.

    Fig.10.Overall morphology and partially enlarged detail of FCG fracture for ((a),(e)) R-T specimens in air,((b),(f)) N-T specimens in air,((c),(g)) R-T specimens in PBS,((d),(h)) N-T specimens in PBS.

    3.5. Crack propagation path observation

    According to the crack propagation path shown in Fig.11(a–d),the FCG paths of all specimens follow mode-I on the whole.Similar to the features of the aforementioned fractures,the FCG path in air of specimens with the same orientation is straighter than that in PBS.Under the same test environment,the FCG path of R-T specimens is significantly smoother than that of N-T specimens.Especially for the N-T specimens in PBS shown in Fig.11(d),the crack path appears obvious localized deflections,and secondary cracks or bifurcation cracks.

    Fig.11.Overall fatigue crack propagation path: (a) R-T specimen in air,(b) N-T specimen in air,(c) R-T specimen in PBS,and (d) N-T specimen in PBS.

    The microstructure near the propagation path at differentΔK(4 MPa?m1/2,6 MPa?m1/2,9.5 MPa?m1/2) is further examined,and the twins are marked in red as shown in Fig.12.Whether in air or PBS,the main fracture mode is the transgranular fracture,and occasionally intergranular fracture occurs at the position of cluster fine-grains.For the same specimen orientation,the microstructure near the crack is similar in the different environments,which indicates that the corrosion environment won’t significantly change their main plastic deformation mechanisms.As shown in Fig.12(a,c),for the R-T specimens under the smallerΔK,no obvious twins are observed near the crack.UntilΔKreaches 9.5 MPa?m1/2,a few fine twins can be found along the crack edge.Jiang et al.also found similar results for R-T specimens,in which twins are difficult to be activated.Only ifΔKis up to 7 MPa?m1/2only,obvious residual twins can be observed near the crack using OM [9].In PBS,underΔK=9.5 MPa?m1/2,many secondary cracks appear along the edge of the grains and residual twins,which indicates that the coupling of corrosion and stress promotes the tendency of intergranular cracking.For the N-T specimens in Fig.12(b,d),the number and size of twins near the crack increase with theΔK.In addition,the number of twins within the 20 μm area of the crack edge of N-T specimens in PBS (Fig.12(d)) is more than that in air(Fig.12(b)),which will be explained in the following discussion.

    Fig.12.The microstructure along the propagation path at different ΔK: (a) R-T specimens in air,(b) N-T in air,(c) R-T specimens in PBS,(d) N-T specimens in PBS.

    4.Discussion

    4.1. Fatigue crack growth behaviors of different specimen orientations in air

    In the stable FCG stage (Paris region),the crack growth is closely related to the plastic zone at the crack tip.According to the alternating slip model [18],irreversible hardening slip and reversible non-hardening slip successively emanate on the specific system under tension loading.In the following unloading stage,the non-hardening system slip completely returns to the crack tip,while the irreversible slip is retained,which leads to crack propagation along a specific plane.Hardening and non-hardening slip appear alternately on both sides of the crack tip,which leaves regular strips on the fracture surface.To further understand the plastic deformation mechanisms of different specimen orientations in air,the surface morphology and crystal orientation near the crack tip under smallΔKare observed and shown in Fig.13.Deformation traces on the specimens’ surface are marked with three black straight lines,and the corresponding crystal plane traces parallel to the surface deformation traces are marked with a color line.For the R-T specimens shown in Fig.13(a),the plastic zone near the crack tip is very small,and the surface deformation traces are mainly distributed at the edge of the crack in an emission shape.Due to the smallerΔKand hard orientation,no extension twins are directly observed in Fig.13(a).Only by carefully observing the EBSD results(Fig.13(c)) can a small number of fine residual twins be found.The surface deformation traces or bifurcated secondary cracks are mainly parallel to the {101} and {0001} crystal planes,which is similar to the research results of pure magnesium single crystal by Ando et al.[19].At room temperature,this oriented Mg alloy can activate the pyramidal({101}<23>) slips in fatigue crack propagation,while the prismaticslip is rarely active [20].The basalslips can be attributed to the non-ideal crystal orientation of the actual material.It can be seen from Fig.13(a) that the Schmidt factor (SF) at the basalslips can reach 0.29.Ando et al.attributed the basaldislocations found to the interaction of two symmetrical second-order pyramidaldislocations [19].The fatigue cracks mainly propagate along the {110}<100>or {100}<110>,which is consistent with the observation in rolled AZ31Mg alloy of Jiang et al.[9].In addition,microcracks along grain boundaries are also observed at local positions,which indicates that the slip systems available for the current orientation are scarce.

    Fig.13.The surface morphology and crystal orientation near the crack tip under small ΔK: (a) and (c) for R-T specimens;(b) and (d) for N-T specimens.

    For the N-T specimens shown in Fig.13(b),the plastic zone around the crack tip is significantly larger,and the surface deformation traces have a wider distribution range.Since the deformation twinning is in soft orientation,a large number of coarse {102} extension twins (regions painted red) appear at both sides away from the crack path and on the lateral leading edge of the crack tip.Similarly,due to the non-ideal crystal orientation,the local positions also activate the basalslips.Moreover,some surface wrinkled traces parallel to the {102} crystal plane can be found at the edge of the crack.This finding indicates that the detwinning process occurs in these regions,which can be attributed to the crack closure effect.The appearance of a large number of extension twins results in the out-of-plane shrinkage near the crack path due to the rotation of crystals.During the cyclic unloading stage,the separated crack surface will close in advance and squeeze each other,thereby causing the de-twinning behavior at the crack edge.As shown in Fig.13(b),the crack path has obvious signs of squeezing each other.

    The difference in FCG rate between the two specimen orientations in air can be explained from the following two aspects.One is the aforementioned crack closure effect,which can reduce the effective stress intensity factor range of N-T specimens,thus weakening the driving force of crack growth.In addition,from the perspective of energy conservation,the energy input from the outside to the system can be used in three parts: increasing the strain energy of the system,forming the new crack surfaces,and changing the kinetic energy of the system.Among them,the change of kinetic energy can be ignored during the stable cracking stage.Comparing Fig.13(a,b),it is not difficult to find that the deformation of the crack tip of the N-T specimen is larger than that of the R-T specimen,that is,the energy that can be used to form the new crack surface of the N-T specimen is less than that of the R-T specimen,which leads to the lower FCG rate of N-T specimen.

    4.2. Corrosion and FCG behaviors of different specimen orientations in PBS

    In an aqueous solution,the Mg matrix and the secondphase form the corrosion galvanic cell,and the following reactions occur:

    the deposition of Mg(OH)2on the surface (Fig.6(a,b)) can protect the Mg matrix.When the deposition and dissolution of Mg(OH)2reach equilibrium,the corrosion rate of the material enters the stable stage.In the stable corrosion stage,the difference in corrosion rate between the two surfaces is mainly concerned with the crystal orientation.Fig.4(b) shows that the TS is mainly composed of {100} and {110} prismatic planes.Combined with the pole figure in Fig.4(c),it is clear that the RS is mainly composed of {0001} basal planes.Owing to higher binding energy and lower surface energy of the basal planes with higher atomic arrangement density,they often display higher corrosion resistance than non-basal planes.Theoretically,the dissolution rate of prismatic planes is 18–20 times higher than that of basal planes [21].However,as shown in Fig.6,the actual corrosion rate of TS specimens is about twice that of RS specimens due to the influence of microstructure defects (vacancies,dislocations,etc.),secondphase galvanic corrosion,and nonideal texture in polycrystalline materials.It is worth noting that the corrosion rate between RS and TS specimens is no obvious difference in the early stage of the hydrogen evolution test,which may be related to the participation of filiform corrosion.Our previous research [15]results show that ZK60 Mg alloy displays filiform corrosion at the initial stage of immersion in neutral PBS.Because the main factor controlling filiform corrosion is the concentration of Cl-adsorbed on the head of filiform corrosion,which is influenced by the properties of the solution itself (such as pH value,Cl-concentration,etc.),the two specimens show similar hydrogen evolution rate in the early corrosion stage.With the extension of immersion time,pitting corrosion and filiform corrosion develop synchronously,the filiform corrosion area will be covered by shallow pitting pits.After entering the stable corrosion stage,pitting corrosion replaces filiform corrosion as the main corrosion form.For the TS specimens,the rapid entry into the stable corrosion stage makes the pitting corrosion on its surface more serious shown in Fig.7(f).

    According to the EIS results shown in Fig.8(a),no obvious low-frequency capacitive reactance arc related to surface passivation film is found in the early corrosion stage.In addition,at the anodic polarization stage of Tafel curves shown in Fig.8(e),no breakdown potential is found.These results fully indicate that no strong protection passivation film can be formed for the two specimens.The weak corrosion film provides channels for hydrogen atoms generated by Eq.(6) to enter the Mg matrix.Solute hydrogen atoms diffuse between the lattice gaps,and many hydrogen traps with low potential energy near the crack tip,such as dislocations,(twin) grain boundaries,microcracks,and voids,will capture more solute hydrogen atoms.In the deformation process of materials,the diffused hydrogen atoms can significantly promote the nucleation and movement of dislocations,and dislocation will also carry hydrogen atoms to migrate to (twin) grain boundaries,microcracks,voids,and other defect locations faster [17].In PBS,the localized accumulation of dislocations and hydrogen atoms will significantly increase the stress concentration at defect positions,which makes the matrix crack under lower stress and display an obvious embrittlement tendency.The secondary cracks shown in Fig.10(g,h) are favorable evidence of matrix embrittlement.Moreover,the increase in the stress concentration at defect positions also induces more deformation twins for N-T specimens in PBS.During the stage of cyclic tensile unloading,the compression stress caused by the closure of the crack surfaces can be released by the secondary cracks near the crack path and restrain the occurrence of the de-twinning process.The increase in twinning driving force and the suppression to de-twinning make the number of twins within the 20 μm area of the crack edge increase obviously for the N-T specimen in PBS shown in Fig.12(d).

    The crack growth rate of specimens in PBS is greatly higher than that in air,which can be attributed to the following two aspects.Firstly,during the process of FCG in PBS,the crack tip is in the cycle of matrix dissolution,film formation,and rupture.The dissolution of the fresh matrix at the crack tip will accelerate the growth of the crack [22].Secondly,the matrix embrittlement caused by the interaction between hydrogen atoms and cyclic stress can also significantly deteriorate the FCG resistance of the materials.In addition,for the high strength Mg alloys,the crack growth rateof corrosion fatigue can be obtained by the following formula [23]:

    5.Conclusion

    According to the results of corrosion and FCG experiments of different specimen orientations for rolled ZK60 Mg alloys,the following conclusions can be drawn:

    (1) The specimen orientations have a significant influence on the stable corrosion rate (corrosion rate of RS is about half of that of TS),but the truth is the opposite for the initial corrosion rate,which can be attributed to the participation of filiform corrosion in the initial stage of immersion.Compared with RS,TS displays more serious pitting corrosion pits on its surface.

    (2) For the same specimen orientation,the FCG paths and fracture surfaces in air are smoother than those in PBS.For the same testing environment,the FCG paths and fracture surfaces of R-T specimens are greatly smoother than that of N-T specimens.Transgranular cracking is the primary cracking mode for different specimen orientations and testing environments.However,the introduction of a corrosion environment increases the intergranular cracking trend.

    (3) The specimen orientations have a significant influence on the plastic deformation mechanisms near the crack tip in air.For R-T specimens,the pyramidaland basalslips induced by the high-stress concentration of crack tip are the main deformation systems,thereby resulting in the small plastic zone.For N-T specimens,the reversible twinning and detwinning dominate the plastic deformation of the crack tip,forming a large plastic zone.Owing to the crack closure effect and activation of deformation twinning,N-T specimens display a slower crack growth rate than R-T specimens.

    (4) In PBS,the main plastic deformation mechanisms of the two specimen orientations haven’t significantly changed.However,the increase in the stress concentration at the crack tip defect position makes the matrix display an obvious embrittlement tendency.For N-T specimens,the obvious increase in residual twins at the crack edge can be attributed to the comprehensive effects of the enhanced twinning driving force and the suppressed detwinning behaviors.

    (5) Owing to the combined effect of hydrogen-induced cracking and anodic dissolution,the FCG rate of specimens in PBS is greatly higher than that in air.Since the crack tips are always in the early corrosion stage under the influence of cyclic loading.The similar corrosion rates in this stage lead to similar stress corrosion cracking performance of two oriented specimens,which dominates the total corrosion FCG rate.Therefore,specimens with different orientations display close FCG rates in PBS.

    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

    This study was sponsored by the National Natural Science Foundation of China (Nos.52175143 and 51571150)

    男女床上黄色一级片免费看| 欧美av亚洲av综合av国产av| 日韩精品免费视频一区二区三区| 亚洲国产看品久久| 欧美黑人精品巨大| 两个人免费观看高清视频| 国产91精品成人一区二区三区| 亚洲精品粉嫩美女一区| 亚洲狠狠婷婷综合久久图片| 亚洲熟妇中文字幕五十中出| 久久精品国产亚洲av香蕉五月| 亚洲男人天堂网一区| 纯流量卡能插随身wifi吗| 久久人人精品亚洲av| 国产一卡二卡三卡精品| 在线天堂中文资源库| 午夜福利免费观看在线| 久久精品国产亚洲av高清一级| 真人一进一出gif抽搐免费| 自拍欧美九色日韩亚洲蝌蚪91| 国产1区2区3区精品| 亚洲人成网站在线播放欧美日韩| 搡老岳熟女国产| 久久久国产精品麻豆| 男人的好看免费观看在线视频 | 色在线成人网| 国产一区二区三区综合在线观看| 日本 欧美在线| 久久香蕉激情| 久久久国产成人免费| 两人在一起打扑克的视频| 成熟少妇高潮喷水视频| 久久午夜综合久久蜜桃| 久久青草综合色| 一级毛片高清免费大全| 欧美日韩瑟瑟在线播放| 一区二区三区激情视频| 亚洲五月色婷婷综合| 久久国产乱子伦精品免费另类| 国产成人精品在线电影| 欧美日韩一级在线毛片| 级片在线观看| 亚洲精品在线美女| 深夜精品福利| 在线观看日韩欧美| 午夜日韩欧美国产| 最好的美女福利视频网| 成人免费观看视频高清| 桃色一区二区三区在线观看| 欧美国产精品va在线观看不卡| 色在线成人网| 亚洲五月天丁香| 一卡2卡三卡四卡精品乱码亚洲| 欧美 亚洲 国产 日韩一| 涩涩av久久男人的天堂| 一夜夜www| 最新在线观看一区二区三区| 欧美激情久久久久久爽电影 | 麻豆国产av国片精品| 怎么达到女性高潮| 日韩欧美国产一区二区入口| 久久中文字幕人妻熟女| 校园春色视频在线观看| 午夜两性在线视频| 国产一区二区三区综合在线观看| aaaaa片日本免费| 亚洲av片天天在线观看| 丁香欧美五月| 电影成人av| 亚洲国产中文字幕在线视频| 超碰成人久久| 婷婷六月久久综合丁香| 国产色视频综合| 免费观看精品视频网站| 悠悠久久av| 久久精品影院6| 99在线视频只有这里精品首页| 我的亚洲天堂| 此物有八面人人有两片| 欧美日韩乱码在线| 法律面前人人平等表现在哪些方面| 亚洲午夜精品一区,二区,三区| 国产又爽黄色视频| 免费看十八禁软件| 免费一级毛片在线播放高清视频 | 色综合欧美亚洲国产小说| 在线十欧美十亚洲十日本专区| avwww免费| 人成视频在线观看免费观看| 日韩大尺度精品在线看网址 | 国产亚洲欧美精品永久| 国产精品99久久99久久久不卡| 一进一出好大好爽视频| 好男人在线观看高清免费视频 | 老司机靠b影院| 一二三四社区在线视频社区8| 国产欧美日韩一区二区精品| 国产成+人综合+亚洲专区| 在线视频色国产色| 视频在线观看一区二区三区| 午夜久久久久精精品| 亚洲三区欧美一区| 久久久久久久午夜电影| 精品欧美一区二区三区在线| 国产亚洲精品久久久久久毛片| 日韩欧美三级三区| 午夜老司机福利片| 欧美国产精品va在线观看不卡| 日本 av在线| 中国美女看黄片| 9色porny在线观看| 人人妻人人澡欧美一区二区 | 午夜两性在线视频| or卡值多少钱| 精品欧美一区二区三区在线| 亚洲中文字幕日韩| 岛国视频午夜一区免费看| av福利片在线| 色老头精品视频在线观看| 久久精品aⅴ一区二区三区四区| 99riav亚洲国产免费| 久久久久久久久免费视频了| 18禁观看日本| 国产高清有码在线观看视频 | 久久久久九九精品影院| 看片在线看免费视频| 午夜免费成人在线视频| а√天堂www在线а√下载| 亚洲精品国产一区二区精华液| 久久香蕉精品热| 国产精品二区激情视频| 国产亚洲欧美98| 欧美国产精品va在线观看不卡| 日韩欧美一区二区三区在线观看| 精品熟女少妇八av免费久了| 国产亚洲精品av在线| 国产成人系列免费观看| 国产av又大| 精品卡一卡二卡四卡免费| 免费少妇av软件| 日韩大码丰满熟妇| 亚洲午夜理论影院| 宅男免费午夜| 日本vs欧美在线观看视频| 欧美激情久久久久久爽电影 | 黑人操中国人逼视频| 丁香欧美五月| 高清毛片免费观看视频网站| 91在线观看av| 黄片小视频在线播放| 99香蕉大伊视频| 精品久久久精品久久久| 亚洲男人天堂网一区| 无人区码免费观看不卡| 97人妻天天添夜夜摸| 亚洲,欧美精品.| 免费在线观看黄色视频的| 好看av亚洲va欧美ⅴa在| 亚洲国产看品久久| bbb黄色大片| 亚洲自偷自拍图片 自拍| 日韩欧美一区视频在线观看| 69精品国产乱码久久久| 好看av亚洲va欧美ⅴa在| 国产蜜桃级精品一区二区三区| 久久性视频一级片| 欧美性长视频在线观看| 亚洲久久久国产精品| 精品无人区乱码1区二区| 欧美国产精品va在线观看不卡| 国产精品,欧美在线| 99香蕉大伊视频| 老司机深夜福利视频在线观看| 午夜免费观看网址| 桃色一区二区三区在线观看| 国产又爽黄色视频| 男人操女人黄网站| 国产成人免费无遮挡视频| 久久久国产欧美日韩av| 欧美成狂野欧美在线观看| 亚洲av片天天在线观看| 亚洲男人天堂网一区| 亚洲欧美激情综合另类| 国产99白浆流出| 久久久精品欧美日韩精品| 熟女少妇亚洲综合色aaa.| 久久精品91蜜桃| 国产成人精品在线电影| 在线播放国产精品三级| 乱人伦中国视频| 色精品久久人妻99蜜桃| 黑人操中国人逼视频| 少妇被粗大的猛进出69影院| 精品国产一区二区久久| 亚洲一区中文字幕在线| 成人亚洲精品一区在线观看| 亚洲天堂国产精品一区在线| 在线视频色国产色| 香蕉国产在线看| 韩国av一区二区三区四区| 日本 av在线| 精品欧美国产一区二区三| 欧美乱码精品一区二区三区| 1024香蕉在线观看| 97超级碰碰碰精品色视频在线观看| 黄色成人免费大全| 叶爱在线成人免费视频播放| 夜夜躁狠狠躁天天躁| 久久国产乱子伦精品免费另类| 一本久久中文字幕| 欧美 亚洲 国产 日韩一| 日本精品一区二区三区蜜桃| 色av中文字幕| 日韩欧美国产一区二区入口| 午夜福利在线观看吧| 亚洲国产欧美网| 久热这里只有精品99| 18禁黄网站禁片午夜丰满| 亚洲一区中文字幕在线| 成人亚洲精品一区在线观看| 欧美一级a爱片免费观看看 | 男女午夜视频在线观看| 麻豆国产av国片精品| 女性被躁到高潮视频| 国产一区二区在线av高清观看| 婷婷丁香在线五月| 午夜福利高清视频| 黄色成人免费大全| 午夜a级毛片| 久99久视频精品免费| 不卡av一区二区三区| 国产成年人精品一区二区| 最新在线观看一区二区三区| 精品电影一区二区在线| 久久久久九九精品影院| 老熟妇乱子伦视频在线观看| av在线天堂中文字幕| 精品久久久久久久久久免费视频| 国产精品香港三级国产av潘金莲| 久久久精品国产亚洲av高清涩受| 精品一区二区三区视频在线观看免费| 免费观看人在逋| 精品一区二区三区av网在线观看| 国产精品久久久人人做人人爽| 国产成人av教育| 欧美 亚洲 国产 日韩一| АⅤ资源中文在线天堂| 丰满人妻熟妇乱又伦精品不卡| 免费看a级黄色片| 男人的好看免费观看在线视频 | 男男h啪啪无遮挡| 亚洲激情在线av| 狂野欧美激情性xxxx| 老汉色∧v一级毛片| 高清黄色对白视频在线免费看| 啦啦啦韩国在线观看视频| 国产又爽黄色视频| 成人国语在线视频| 伊人久久大香线蕉亚洲五| 久久天躁狠狠躁夜夜2o2o| 搡老熟女国产l中国老女人| 久久精品国产亚洲av香蕉五月| 国产亚洲欧美98| 悠悠久久av| 麻豆av在线久日| 99久久99久久久精品蜜桃| 宅男免费午夜| 久久久久精品国产欧美久久久| 亚洲自拍偷在线| 波多野结衣一区麻豆| 亚洲一区中文字幕在线| 极品教师在线免费播放| 亚洲国产欧美网| www.999成人在线观看| 欧美中文综合在线视频| 亚洲精品在线观看二区| 热99re8久久精品国产| 久久久久久免费高清国产稀缺| 狠狠狠狠99中文字幕| 丁香六月欧美| 国产高清激情床上av| 制服诱惑二区| av电影中文网址| 国产在线观看jvid| 一边摸一边抽搐一进一小说| 高清在线国产一区| 久久香蕉精品热| 国产精品自产拍在线观看55亚洲| 欧美日韩亚洲国产一区二区在线观看| 中文亚洲av片在线观看爽| 88av欧美| 亚洲五月婷婷丁香| 精品第一国产精品| 夜夜看夜夜爽夜夜摸| 老汉色av国产亚洲站长工具| 夜夜夜夜夜久久久久| 欧美精品啪啪一区二区三区| 亚洲国产精品sss在线观看| 男人的好看免费观看在线视频 | 免费高清视频大片| 亚洲aⅴ乱码一区二区在线播放 | 中文字幕另类日韩欧美亚洲嫩草| 99久久精品国产亚洲精品| 免费无遮挡裸体视频| 久久精品人人爽人人爽视色| 在线观看日韩欧美| 午夜福利影视在线免费观看| 国产又色又爽无遮挡免费看| 午夜精品久久久久久毛片777| 国产亚洲精品av在线| 在线永久观看黄色视频| 欧美人与性动交α欧美精品济南到| 91老司机精品| 亚洲专区中文字幕在线| 日韩欧美免费精品| 国产高清videossex| 悠悠久久av| 亚洲一码二码三码区别大吗| 午夜免费鲁丝| 九色亚洲精品在线播放| 在线观看一区二区三区| 老司机深夜福利视频在线观看| 亚洲片人在线观看| 一级作爱视频免费观看| 亚洲av片天天在线观看| 天天一区二区日本电影三级 | 后天国语完整版免费观看| 窝窝影院91人妻| 亚洲国产精品999在线| 老鸭窝网址在线观看| 亚洲精品av麻豆狂野| 99re在线观看精品视频| 人人妻人人澡人人看| 欧美不卡视频在线免费观看 | av天堂久久9| 村上凉子中文字幕在线| 欧美成人一区二区免费高清观看 | 麻豆成人av在线观看| 他把我摸到了高潮在线观看| 两性夫妻黄色片| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美老熟妇乱子伦牲交| 69精品国产乱码久久久| 天天一区二区日本电影三级 | 很黄的视频免费| 黄网站色视频无遮挡免费观看| 亚洲av电影在线进入| 国产精品一区二区在线不卡| 丝袜在线中文字幕| 午夜福利免费观看在线| 国产麻豆成人av免费视频| 又紧又爽又黄一区二区| 午夜福利,免费看| 国产极品粉嫩免费观看在线| 欧美不卡视频在线免费观看 | 黄色 视频免费看| 成人亚洲精品一区在线观看| 亚洲国产中文字幕在线视频| 黄色 视频免费看| 国产精品 国内视频| 一级毛片精品| 国产又爽黄色视频| 视频在线观看一区二区三区| 热re99久久国产66热| 久久狼人影院| 真人一进一出gif抽搐免费| 亚洲va日本ⅴa欧美va伊人久久| 亚洲精品美女久久久久99蜜臀| 成年女人毛片免费观看观看9| 色哟哟哟哟哟哟| 69av精品久久久久久| 精品久久久久久久人妻蜜臀av | 搡老妇女老女人老熟妇| 怎么达到女性高潮| av在线播放免费不卡| 亚洲人成电影免费在线| 亚洲一区中文字幕在线| 露出奶头的视频| 国产人伦9x9x在线观看| 日韩成人在线观看一区二区三区| 大香蕉久久成人网| 一本综合久久免费| 国产av精品麻豆| 国产亚洲精品综合一区在线观看 | 在线免费观看的www视频| 久久久久精品国产欧美久久久| 国产欧美日韩一区二区三| 欧美大码av| 国产私拍福利视频在线观看| 日韩大尺度精品在线看网址 | 天天一区二区日本电影三级 | 黄色 视频免费看| 中文字幕人妻熟女乱码| 成人国产一区最新在线观看| 老司机午夜十八禁免费视频| 搡老岳熟女国产| 亚洲精品国产一区二区精华液| 在线观看免费视频网站a站| 色尼玛亚洲综合影院| 久久婷婷人人爽人人干人人爱 | 亚洲无线在线观看| 午夜福利影视在线免费观看| 国产精品久久久久久人妻精品电影| 欧美不卡视频在线免费观看 | 满18在线观看网站| 侵犯人妻中文字幕一二三四区| 黄片大片在线免费观看| 国产成人欧美| 99久久国产精品久久久| 亚洲精品一区av在线观看| 亚洲成国产人片在线观看| 国产又爽黄色视频| 日韩欧美国产在线观看| 在线天堂中文资源库| 男男h啪啪无遮挡| 丝袜美腿诱惑在线| 欧美日韩精品网址| 久久久久九九精品影院| 纯流量卡能插随身wifi吗| 性少妇av在线| 国产xxxxx性猛交| svipshipincom国产片| 色综合亚洲欧美另类图片| 国产在线精品亚洲第一网站| 19禁男女啪啪无遮挡网站| 国产成人免费无遮挡视频| 国产av精品麻豆| 在线十欧美十亚洲十日本专区| 国产精品 国内视频| 夜夜看夜夜爽夜夜摸| 午夜精品久久久久久毛片777| 日韩欧美国产一区二区入口| 自拍欧美九色日韩亚洲蝌蚪91| avwww免费| 成人18禁高潮啪啪吃奶动态图| 国产一区二区三区综合在线观看| 咕卡用的链子| 国产精品 国内视频| 亚洲成人免费电影在线观看| 18禁黄网站禁片午夜丰满| 亚洲自拍偷在线| 欧洲精品卡2卡3卡4卡5卡区| 欧美午夜高清在线| 脱女人内裤的视频| 91av网站免费观看| 97人妻天天添夜夜摸| 91国产中文字幕| 一进一出抽搐gif免费好疼| 国产成人一区二区三区免费视频网站| 久久国产精品男人的天堂亚洲| 国产成人啪精品午夜网站| 国产精品久久久av美女十八| 亚洲第一欧美日韩一区二区三区| 亚洲激情在线av| 啦啦啦韩国在线观看视频| 大陆偷拍与自拍| 成人国产综合亚洲| 如日韩欧美国产精品一区二区三区| 嫩草影视91久久| 中出人妻视频一区二区| 精品乱码久久久久久99久播| 看片在线看免费视频| 搡老妇女老女人老熟妇| 亚洲中文字幕日韩| 亚洲人成伊人成综合网2020| 午夜a级毛片| 亚洲在线自拍视频| 亚洲九九香蕉| 国产亚洲av高清不卡| 自线自在国产av| 老司机午夜福利在线观看视频| 色综合亚洲欧美另类图片| 久久青草综合色| 亚洲九九香蕉| 精品少妇一区二区三区视频日本电影| 女警被强在线播放| 一级a爱片免费观看的视频| av天堂久久9| 欧美日本中文国产一区发布| 日韩国内少妇激情av| 十八禁人妻一区二区| 91字幕亚洲| 中文字幕色久视频| 岛国视频午夜一区免费看| 亚洲国产欧美网| cao死你这个sao货| 欧美日韩黄片免| x7x7x7水蜜桃| 在线观看免费视频网站a站| 亚洲欧美日韩高清在线视频| 天天添夜夜摸| 777久久人妻少妇嫩草av网站| 波多野结衣一区麻豆| 色在线成人网| 精品第一国产精品| 俄罗斯特黄特色一大片| 国产亚洲av高清不卡| 三级毛片av免费| 日本欧美视频一区| 亚洲在线自拍视频| 免费观看人在逋| 亚洲欧美激情在线| 久久久精品国产亚洲av高清涩受| 久久精品aⅴ一区二区三区四区| 欧美大码av| 国产精品爽爽va在线观看网站 | 午夜视频精品福利| 成年女人毛片免费观看观看9| 亚洲第一青青草原| netflix在线观看网站| 一级,二级,三级黄色视频| 18禁观看日本| √禁漫天堂资源中文www| 免费高清视频大片| 欧美国产精品va在线观看不卡| 久久九九热精品免费| 91成人精品电影| 久久午夜综合久久蜜桃| 91成年电影在线观看| 国产精品一区二区精品视频观看| 黄色片一级片一级黄色片| 一区福利在线观看| 精品久久久久久久人妻蜜臀av | 无限看片的www在线观看| 一本久久中文字幕| 午夜福利18| 级片在线观看| 99国产综合亚洲精品| 亚洲人成网站在线播放欧美日韩| 性少妇av在线| 黄片小视频在线播放| 亚洲国产欧美网| 久久精品成人免费网站| 国产男靠女视频免费网站| 中文字幕最新亚洲高清| 欧美日韩亚洲综合一区二区三区_| 中文亚洲av片在线观看爽| 夜夜夜夜夜久久久久| 国产精品,欧美在线| 香蕉久久夜色| 高清毛片免费观看视频网站| 午夜福利18| 91老司机精品| 最近最新中文字幕大全免费视频| 国产高清视频在线播放一区| 精品久久久久久成人av| 啦啦啦 在线观看视频| 亚洲色图av天堂| 精品熟女少妇八av免费久了| 欧美黑人精品巨大| 真人做人爱边吃奶动态| 国产精品久久视频播放| 免费看a级黄色片| 国产在线精品亚洲第一网站| 日本a在线网址| 亚洲精品中文字幕在线视频| 亚洲av电影不卡..在线观看| 变态另类丝袜制服| 国产精品98久久久久久宅男小说| 性少妇av在线| 国产成年人精品一区二区| 欧美精品啪啪一区二区三区| 国产精品美女特级片免费视频播放器 | 人人妻人人澡人人看| 国产精品久久久人人做人人爽| 国产亚洲精品av在线| 国产亚洲av嫩草精品影院| 久久国产精品男人的天堂亚洲| 一级毛片精品| 1024视频免费在线观看| 久久人人精品亚洲av| 色尼玛亚洲综合影院| 国产精品自产拍在线观看55亚洲| 在线观看www视频免费| 精品午夜福利视频在线观看一区| 国产一区二区在线av高清观看| 国产精品1区2区在线观看.| 99久久精品国产亚洲精品| 人人妻人人爽人人添夜夜欢视频| 757午夜福利合集在线观看| 一进一出抽搐动态| 91麻豆精品激情在线观看国产| 无限看片的www在线观看| 欧美人与性动交α欧美精品济南到| 在线观看免费视频网站a站| 亚洲av成人一区二区三| 给我免费播放毛片高清在线观看| 一a级毛片在线观看| 免费一级毛片在线播放高清视频 | 嫩草影视91久久| 国产亚洲欧美精品永久| 国产黄a三级三级三级人| 亚洲av五月六月丁香网| 性少妇av在线| 国产欧美日韩精品亚洲av| 国产精品国产高清国产av| x7x7x7水蜜桃| 亚洲成av片中文字幕在线观看| 757午夜福利合集在线观看| 午夜老司机福利片| 亚洲男人天堂网一区| 欧美不卡视频在线免费观看 | 99精品久久久久人妻精品| 国产精品久久久人人做人人爽| 中文字幕精品免费在线观看视频| 日韩有码中文字幕| 亚洲第一青青草原| 亚洲午夜精品一区,二区,三区| 日本撒尿小便嘘嘘汇集6| 国产欧美日韩精品亚洲av| 亚洲专区国产一区二区| 国产成人欧美在线观看| 黄色丝袜av网址大全| 91av网站免费观看| 99久久久亚洲精品蜜臀av|