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

    Sodium dodecyl sulfate (SDS) as an effective corrosion inhibitor for Mg-8Li-3Al alloy in aqueous NaCl: A combined experimental and theoretical investigation

    2023-03-24 09:24:34HonggunSongZhidongXuLhouriBenouZhengYinHongyuGunHongYnLuoChoZhiHuXudongWng
    Journal of Magnesium and Alloys 2023年1期

    Honggun Song ,Zhidong Xu ,Lhouri Benou ,Zheng Yin ,Hongyu Gun ,Hong Yn ,Luo Cho,Zhi Hu,?,Xudong Wng,?

    a Department of Mechanical and Electrical Engineering,Nanchang University,Nanchang,330031,China

    b Université Paris-Saclay,UVSQ,LISV,78124,Vélizy -Villacoublay,France

    c Institute for Advanced Study,Nanchang University,Nanchang,330031,China

    Abstract The corrosion inhibition behavior of Mg-8Li-3Al alloy in NaCl solution with sodium dodecyl sulfate (SDS) was investigated by hydrogen analysis,scanning electron microscopy (SEM),electrochemical test,scanning Kelvin probe force microscopy (SKPFM) and computational methods.Results showed that the corrosion resistance of Mg-8Li-3Al alloy in NaCl solution was effectively improved with SDS.The SEM and SKPFM results confirme a dense,200 nm-thick SDS-adsorbed layer had formed on the alloy surface.The separation energy ΔEgap and adsorption energy Eads of SDS on the Mg surface were calculated by density functional theory and molecular dynamics simulations,respectively.And the corrosion inhibition mechanism was hypothesized and described.

    Keywords: Mg-Li alloy;Corrosion;Sodium dodecyl sulfate;Scanning kelvin probe force microscopy;Density functional theory;Molecular dynamics.

    1.Introduction

    Mg-Li alloy is one of the lightest structural metal alloys and has a wide range of application prospects in light industrial applications,aerospace,and the automotive industries[1–6];however,the widespread use of Mg-Li alloys is still limited due to their relatively high corrosion rates[7,8].Many coating technologies have been developed to improve the corrosion resistance of traditional Mg-Li alloys,including corrosion inhibitors,physical vapor deposition,conversion coating,anodizing,plasma electrolytic oxidation,and electroplating[9–15].Among them,corrosion inhibitors are usually used to improve corrosion resistance due to their environmental friendliness and low costs.

    Corrosion inhibitors have shown adsorption ability on the surface of Mg alloys,leading to improved corrosion resistance of the alloys [16–17].Hu et al.[18] synthesized an environmentally-friendly tetraphenylporphyrin and studied its inhibition effect on AZ91D magnesium alloys in 0.05 wt.% NaCl solution.They claimed that the inhibition corrosion efficien y of tetraphenylporphyrin reached 90%,because the inhibitor chelated with Mg ions to form a fil able to retard Mg dissolution.Sodium dodecyl sulfate(SDS)was used as an effective corrosion inhibitor for AZ91 Mg alloys in a study by Lu et al.,which found that the generated H2from the alloy immersed in 3.5 wt.% NaCl for 48 h with 0.05 M SDS addition was reduced from 9.1 to 0.3 mL/cm2[19].The corresponding inhibition mechanism of SDS on Mg alloy was shown to occur through chemical adsorption of SDS on AlMn cathodic intermetallics to inhibit micro-galvanic corrosion,which then led kinetically inhibited cathodic hydrogen evolution and anodic dissolution of Mg [20].However,there are only a few reports regarding the characterization of the adsorbed layer of corrosion inhibitors and its influenc on micro-galvanic corrosion of alloys,and the theoretical adsorption of inhibitors.

    Theoretical calculations for the adsorption capacity of the inhibitor on the alloy surface have become an important judgment basis for corrosion inhibition evaluation [21,22].The work of Lu et al.[19] fi the Langmuir adsorption isotherm according to the surface coverage of Mg alloy under various concentrations of SDS,and inferred that the standard free adsorption energy of SDS on AZ91 surface is ?27.07 kJ/mol.The inhibition effect of methionine and tryptophane on the corrosion of LA141 Mg-Li alloy sheets in 3.5 wt.% NaCl solution was investigated by molecular dynamics (MD) simulation,and the results showed that the binding energy between the Mg (0 0 0 1) surface and methionine molecule is?12.18 kcal/mol,while the binding energy between the Mg(0 0 0 1) surface and tryptophane molecule is ?34.37 kcal/mol.The combined effect of the physical and chemical absorption on the LA141 surface was considered as the main reason of the suppressed Mg dissolution [23,24].

    SDS with the chemical formula C12H25NaSO4and 0.0082 M critical micelle concentration in pure water at 25 °C,as a kind of anionic surfactant,had been used as an effective corrosion inhibitor for Mg alloys,which consists of a hydrophobic 12-carbon chain and a polar sulfate head group [25].During the process of inhibition for Mg alloy,the polar sulfate head group and hydrophobic 12-carbon chain contacted with the matrix and the solution,respectively.The hydrophobic 12-carbon chain tail,which was perpendicular to the matrix surface,made the water contact angle increased from 73.45°to 118.54°with the increasing of immersion time,thus built up a protective corrosion layer [19].However,the inhibition mechanism of SDS for Mg alloy is not precisely clear.Therefore,in this work,the corrosion inhibition behaviors and mechanism of Mg-8Li-3Al alloy in NaCl aqueous with the addition of SDS were investigated by a combined experimental and molecular simulation approaches (DFT and MD),with specifi focus on the characterization of the SDSadsorbed layer and the elaboration of the corrosion inhibition mechanism of SDS on Mg-8Li-3Al surface.

    2.Experimental

    Using pure magnesium (99.9 wt.%),pure lithium(99.9 wt.%),and pure aluminum (99.9 wt.%) as raw materials,the as-extruded Mg-8Li-3Al sample was prepared.The mixed alloys were held in a protective argon environment at 740 °C for 600 s for complete dissolution.The melt was then poured into a preheated 200 °C,Φ45×120 mm iron mold.Subsequently,the homogenized samples were held at 250 °C for 1 h.Under an extrusion speed of 2.5 mm/s at 250 °C,the sample was processed into cuboids.The actual composition of the tested sample was 7.89 wt.% Li,2.78 wt.% Al and the balance Mg,which was measured by inductively coupled plasma-atomic emission spectroscopy.

    Prior to corrosion testing,all tested alloys were mechanically ground using 600#to 3000#grit SiC paper and 0.5–5 μm Al2O3suspensions.The various concentrations of corrosion inhibitor SDS (0.02–0.06 M,from XiLong Scientifi Co.,Ltd) were mixed into a 3.5 wt.% NaCl solution during the process of immersion testing and electrochemical measurements of the tested alloy.Each tested alloy was soaked in a 3.5 wt.% NaCl solution+(0,0.02,0.04,0.06 M) SDS inhibitor for 1 d The corrosion products on the metal surface are removed by pickling solution.The components of the pickling solution are CrO3and AgNO3(20:1).Five identical samples for each group of experiments were prepared,and the average values were counted.The corrosion rate was evaluated by measuring the hydrogen evolution volume of the specimens,and the inhibition efficien y was calculated using Eq.(1):

    The electrochemical measurements of studied alloys were obtained with a volume of 400 mL aqueous solution(3.5 wt.%NaCl+0,0.02,0.04,0.06 M SDS inhibitor) at 25 °C on an electrochemical workstation (Princeton P4000).The electrochemical measurements on the specimens were performed with three electrodes (the working electrode,reference electrode,and counter electrode).Before the electrochemical measurements,all tested alloys were immersed for 600 s in aqueous solution (3.5 wt.% NaCl+0,0.02,0.04,0.06 M SDS inhibitor) to achieve open circuit potential.For the measurement of the polarization curve,the scanning rate and range were 1 mVSCEs?1and ±350 mVSCE,respectively.In the measurement of electrochemical impedance spectroscopy,the frequency ranged from 105Hz to 10?2Hz,and the disturbance amplitude was 5 mV.The impedance spectrum was analyzed and fitte by Zview software.

    The surface morphology and composition of the studied alloys were characterized by scanning electron microscopy(SEM,TESCAN VEGA 3)with energy dispersive X ray spectroscopy (EDXS).The Volta potentials maps of the studied alloys were measured by the scanning Kelvin probe force microscope (SKPFM,Agilent 5500) at room temperature.A 20 nm thick silicon tip with a Pt/Ir-coating was operated at 75 kHz mechanical drive frequency and 1 V AC modulation at 10 kHz for simultaneous acquisition of topography and surface potential signal.The surface topographies of the studied alloys were characterized in tapping mode,and then the tip was lifted to 100 nm above the alloy surface to record the potential distribution.

    Fig.1.(a) Hydrogen evolution volume,(b) pH change during immersion for 24 h in 3.5 wt.% NaCl solution with and without the addition of SDS.

    Furthermore,Material Studio software was used to carry out all calculations.As it is known that SDS has been widely used as anionic surfactant [26,27],the adsorption of the SDS anion,DS–,is concentrated in this part.The DS–was firs generated and optimized under DFT approximation by Dmol3module,using the GGA-BLYP exchange-correlation functional and the COSMO implicit solvent model.The convergence thresholds for energy,maximum force,and maximum atomic displacement were set equal to 1×10?5Ha,2×10?3Ha?1,and 5×10?3,respectively,with a Fermi smearing of 5×10?3Ha.The adsorption of DS–on the Mg (0 0 0 1) surface,which is the most frequently discussed surface for adsorption research,was then performed by MD methodology with the Forcite module.Herein,the calculation was carried out in a simulation box (dimension 55.6× 32.1×83.4) with periodic boundary conditions.From the top to the bottom,the box consisted of a vacuum layer with a thickness of 30,a solvent layer (with 1800 H2O and one DS–) and an Mg (0 0 0 1) slab.The uppermost part of water and all Mg atoms were fi ed,while other atoms were released during the calculation.The insertion of the vacuum layer and the fixatio of the upper water were meant to eliminate the associated error due to periodic boundary conditions,ensuring that the upper molecules were not influence by the Mg surface.The CVFF (Consistent Valence Force Field) was adopted,and the whole system was calculated under the Andersen thermostat and a canonical (NVT)ensemble at 298.0 K,with a time step of 1.0 fs.The system was simulated for 60 ps to ensure that both the temperature and energy reach equilibrium.Besides,electrostatic interactions were calculated by the Ewald method and van der Waals interactions were calculated by the atom-based method with a cutoff distance of 15.5.Copyright line should read: ?2020 National Engineering Research Center for Magnesium Alloys of China,Chongqing University.

    3.Results

    3.1.Hydrogen evolution measurement

    Fig.1 shows the H2evolution volume and pH of Mg-8Li-3Al alloys in pure NaCl solution and SDS-containing NaCl solution over 24 h.As shown in Fig.1a,the Mg-8Li-3Al specimen immersed in NaCl solution shows a fast H2evolution rate,but the H2evolution rate of samples significantl decreases when immersed in SDS-containing NaCl solution.Moreover,the higher the content of SDS in the solution,the smaller the volume of released H2,which indicates that the effica y of the corrosion inhibitor is affected by the concentration of the corrosion inhibitor.After the specimen immersion,the pH of the NaCl corrosive medium significantl increased,from 6.7 to 10.9,as depicted in Fig.1b.However,the pH of the NaCl+SDS corrosive medium has a little change with the increasing of immersion time.Specially,the pH value of the NaCl+0.06 M SDS corrosive medium is basically unchanged at about 8.9,which agrees with the previous results of the hydrogen evolution experiment.These results suggest that the corrosion behavior of Mg-8Li-3Al alloy is signifi cantly suppressed by the addition of SDS.

    3.2.Surface and cross-sectional morphology study

    Fig.2 shows the scanning electron microscopy (SEM) images of the surface and cross-section of the Mg-8Li-3Al alloy with and without the addition of SDS.To obtain the morphology after corrosion occurred,the Mg-8Li-3Al alloys were immersed in pure NaCl solution and 0.06 M SDS containing NaCl solution for 24 h.As shown in Fig.2a,the morphology of the Mg-8Li-3Al alloy revealed large corrosion pits with ~500 μm lengths and many corrosion products after the sample was immersed in pure NaCl solution.However,only several fin corrosion products were formed on the sample surface that had been immersed in NaCl+0.06 M SDS solution for 24 h (Fig.2b).As exhibited in Fig.2c,when immersed in the pure NaCl solution,larger corrosion pits occurred in the cross-sectional image of the Mg-8Li-3Al alloy.It is worth noting that severe localized corrosion pits were found to extend into the interior of the Mg-8Li-3Al alloy,as shown by the arrow in Fig.2c.Additionally,a uniform corrosion fil was observed on the cross-section SEM micrograph of Mg-8Li-3Al alloy,which had been immersed in NaCl+0.06 M SDS solution for 24 h,as depicted in Fig.2d.

    Fig.2.SEM images of the Mg-8Li-3Al alloy surface (a,b) and cross-section (c,d): (a,c) after immersion for 24 h in blank NaCl solution and (b,d) NaCl during immersion for 24 h with the addition of 0.06 M SDS.

    3.3.Electrochemical polarization test

    Fig.3 shows the potentiodynamic curves of the tested samples after immersion in aqueous solution(3.5 wt.%NaCl+0,0.02,0.04,0.06 M SDS inhibitor).Table 1 lists the fit ting data of polarization curves of the tested alloys.According to Fig.3 and Table 1,the Mg-8Li-3Al alloy with 0.02 M SDS inhibitor (?1.416 VSCE) showed a more noble corrosion potential (Ecorr) than the matrix alloy (?1.448 VSCE).As listed in Table 1,the corrosion current densities(Icorr) of Mg-8Li-3Al alloys with SDS inhibitor are smaller than that of the matrix alloy.Notably,the Mg-8Li-3Al alloy+0.06 M SDS inhibitor has the smallest corrosion current density (Icorr=61.6 mA/cm2).Additionally,the inhibition efficien y (η) of the samples increases with increasing SDS inhibitor concentration.Notably,the inhibition efficien y(η) reaches its maximum value of 89.1% for the sample in 3.5 wt.% NaCl+0.06 M SDS inhibitor.These results suggested that the aqueous solution (3.5 wt.% NaCl+0,0.02,0.04,0.06 M SDS inhibitor) can suppress the corrosion rates of the tested samples.

    Table 1 Electrochemical parameters obtained from the potentiodynamic polarization curves of the tested samples and corresponding inhibition efficien y.

    Table 2 Electrochemical parameters of studied alloys immersed for 1 h.

    Table 3 Electrochemical parameters of studied alloys immersed for 24 h.

    Fig.3.Potentiodynamic curves of the Mg-8Li-3Al alloys after immersion in aqueous solution (3.5 wt.% NaCl+0,0.02,0.04,0.06 M SDS inhibitor).

    3.4.Electrochemical impedance spectroscopy (EIS)

    Fig.4 shows the Nyquist plots and Bode plots,performed in aqueous solution (3.5 wt.% NaCl+0,0.02,0.04,0.06 M SDS inhibitor) at different times (1 h,24 h).As shown in Fig.4a,after 1 h of immersion,the Mg-8Li-3Al alloys were characterized by a capacitive loop and an inductive loop.Additionally,the diameter of the capacitive loops increases with increasing SDS concentration,and the Mg-8Li-3Al alloy with 0.06 M SDS inhibitor shows the largest capacitive loop radius.As exhibited in Fig.4b,the impedance moduli |Z| of the modifie alloys were larger than that of the tested alloy.Similarly,Fig.4c,d show the Nyquist plots and Bode plots of the tested alloys immersed for 24 h.All tested samples immersed for 24 h consist of a capacitive loop at high–frequency spectra,as shown in Fig.4c.It is apparent that the diameters of the capacitive loops in the SDS inhibitor were bigger than that in the pure 3.5 wt.% NaCl solution.Additionally,the corresponding Bode curves display that the tested alloy with 0.06 M SDS inhibitor has the highest value of impedance modulus |Z|,as depicted in Fig.4(d).

    Fig.5 shows the corresponding equivalent circuit of the tested alloys.RS,Rct,andCPEdlare the solution resistance,charge transfer resistance,and double layer capacitance,respectively.Additionally,RLand L represent inductance resistance and inductance,respectively.Polarization resistance(Rp) is calculated as follows [28]:

    Fig.4.EIS plots for the tested alloys after immersion in aqueous solution (3.5 wt.% NaCl+0,0.02,0.04,0.06 M SDS inhibitor): (a) Nyquist plots for 1 h,(b) Bode plots for 1 h,(c) Nyquist plots for 24 h,(d) Bode plots for 24 h.

    Fig.5.The model with the fitte equivalent circuits of tested alloys: (a) samples immersed for 1 h,(a) samples immersed for 24 h.

    The fittin results by Zview 2.8 software are shown in Tables 2 and 3.As listed in Table 2,theRpvalues of the tested samples after immersion in aqueous solution (3.5 wt.%NaCl+0.02,0.04,0.06 M SDS inhibitor) for 1 h are larger than that of the unmodifie alloy.Additionally,theRpvalues of the Mg-8Li-3Al alloy with 0.06 M SDS inhibitor increased to 358.6Ω?cm2,about 2.2 times that of the unmodifie alloy(162.5Ω?cm2).Similarly,theRctvalues of the tested samples after 24 h immersion increased with increasing of SDS concentration,and the Mg-8Li-3Al alloy with 0.06 M SDS inhibitor showed the largest value of 1296Ω?cm2.

    Table 4 The calculated ELUMO,EHOMO,ΔEgap(=ELUMO–EHOMO) and LUMO–HOMO gaps for the Mg-DS– interaction.

    3.5.Langmuir adsorption isotherms

    In order to investigate the adsorption and inhibition mechanism of SDS,the Langmuir adsorption isotherm plots for the SDS at different concentrations is evaluated using hydrogen evolution method,and the results are shown in Fig.6.The plot of Cinh/θagainst inhibitor concentration Cinhyield a straight line with correlation coefficien of 0.997 and slop of 1.17(in Fig.6),confirmin that the adsorption of SDS on Mg-8Li-3Al alloy surface fi the Langmuir adsorption isotherm[29]:

    Fig.6.Langmuir adsorption isotherm plots for the SDS at different concentrations using hydrogen evolution method.

    Fig.7.Main elemental distribution (b) O,(c) Mg,(d) Al and (e) S) on Mg-8Li-3Al surface after immersed for 24 h with the addition of 0.06 M SDS.

    Where R is the universal gas equilibrium constant,T is absolute temperature (298 K) and 55.5 is the molar concentration of water in solution.

    3.6.The corrosion layer and SPKFM study

    3.6.1.The corrosion layer

    To study the surface corrosion influenc of 0.06 M SDS on the Mg-8Li-3Al alloy,Fig.7 shows a BSE-SEM image of the sample and the corresponding EDXS maps.As shown in Fig.7,the surface morphology of the tested alloy revealed that the main elemental distribution is O,Mg,Al and S after the sample had been immersed for 24 h with the addition of 0.06 M SDS.Notably,the bulk material on the sample surface was mainly enriched with elemental S but not elemental Mg,suggesting that SDS inhibitor can completely cover the surface of the alloy in a dense manner,as exhibited in Fig.6a,c,e.

    Fig.8.SKPFM analysis of the Mg-8Li-3Al alloy without the addition of SDS: (a) surface topographic image;(b) 3D topographic map;(c) surface potential image;(d) 3D potential map;(e) potential profile

    3.6.2.Scanning kelvin probe force microscope (SKPFM)

    SKPFMhas become a recent choice to study the weak force of atoms on the surface of the tested alloy,reflectin the surface morphology and surface potential of the sample [31–35].Coy et al.[36–37] used SKPFM to evaluate the role of the phases in Mg alloys on potential galvanic interactions,and the results showed that the smaller Volta potential differences of Al2Gd/α–Mg (250 mV) reduce micro-galvanic corrosion;however,the moniliform AlCuMg phase with the higher Volta potential difference(680 mV)relatively toα–Mg dramatically accelerates micro-galvanic corrosion.Moreover,the SKPFM technique is also able to detect if a metallic surface is active or passive comparing their Volta potential values for corrosion inhibition mechanism [38–39].Vega et al.[38] study the corrosion of a bare metal exposed to a droplet of electrolyte using SKP technique.In the droplet test,the potential value of a mixed surfmer (PMB_S2) coating was below ?0.4 V since the beginning of the measurement and remained below this value during the droplet test.In contrast,the PMB_S1 system showed a positive potential value (around 0.3 V) and it remained in the range of 0.3 ?0.4 V during the entire measurement.According to the SKP potential values and the visual inspection,the metal/coating interface with phosphatization (PMB_S1) did not undergo corrosion compared to unreacted interface (PMB_S2).In this work,the topographic map,surface potential and potential profil of the Mg-8Li-3Al alloy are given and compared with the sample immersed in SDS (Figs.8–10).

    Fig.9.SKPFM analysis of the Mg-8Li-3Al alloy with the addition of SDS: (a) surface topographic image;(b) 3D topographic map;(c) surface potential image;(d) 3D potential map;(e) potential profile

    Fig.8 shows the topographic map,surface potential and potential profil of the tested alloy without the addition of SDS.As shown in Fig.8a,theα-Mg phase in the sample is surrounded by the grayβ-Li phase.It can be inferred that the particle with ~3 μm diameter and ~0.5 μm height is the Al-Li phase,as exhibited in the 3D topographic map of Fig.8b.The Al-Li particles (135 mV) exhibit cathodic behavior compared to theα-Mg phase (25 mV),as shown in Fig.8c–e.However,with respect to the magnesium matrix,theβ-Li inclusion exhibits lower potential differences of 60 mV.

    Fig.9 shows the topographic map,surface potential and potential profil of the tested alloy with the addition of SDS.As shown in Fig.9a,b,a dense and fla covering layer is observed on the surface of the sample with the addition of SDS.The surface potential map and potential profil showed that the potential value of the SDS layer is only about 10 mV(Fig.9c–e).These results suggest that the secondary phases on the surface of the specimens were all covered by SDS inhibitor,thus inhibiting to micro-galvanic corrosion.

    Fig.10 presents the SKPFM images of the SDS layer with rift.As shown in Fig.10a,b,e,it can be seen that the thickness of the SDS layer is about 200 nm on the surface of the tested alloy.Additionally,a rift with about 10 μm length is found in the surface topographic map image(Fig.10b),which indicated that this area is not completely covered by the SDS layer.According to the surface potential images(Fig.10c,d,f),the surface potential near the rift is about 70 mV.Therefore,the micro-galvanic corrosion in the rift of the SDS layer probably occurred when the sample was immersed in solution.

    Fig.10.SKPFM analysis of the SDS layer with rift: (a) surface topographic image;(b) 3D topographic map;(c) surface potential image;(d) 3D potential map;(e) topographic profile (f) potential profile

    Fig.11.(a) Optimized anion structure,(b) HOMO,and (c) LUMO of DS–.

    3.7.Quantum chemical calculations and molecular dynamic simulation

    3.7.1.Quantum chemical calculations

    The relationship between the electronic properties of the surfactant structure and the corrosion inhibition behaviors on the metallic surface are frequently investigated by frontier molecular orbital theory [22,40–45].Fig.11a–c show optimized geometric structures of the DS anion,its highest occupied molecular orbital (HOMO),and its lowest unoccupied molecular orbital (LUMO),respectively.It was found that the distribution of HOMO and LUMO mainly gather around the sulfonic acid group rather than the alkyl tail,indicating preference for the interactions with metallic surfaces [29].In addition,this interaction can be characterized by donoracceptor interactions [46].This means that a higher energy of HOMO (EHOMO) describes the tendency to donate electron(s) to the suitable acceptor molecule with low energy and empty molecular orbital.In contrast,a lower energy of LUMO(ELUMO) represents the ability to accept electron(s).The calculatedELUMO,EHOMO,the separation energy ΔEgap(=ELUMO–EHOMO),and theLUMO–HOMOgaps for the Mg-DS–interaction are given in Table 4.

    It can be demonstrated that asEHOMOof Mg is higher than that of DS–,Mg acts as a Lewis base while DS–plays the role of Lewis acids.Therefore,the HOMO orbital of Mg tends to neutralize the LUMO orbital of DS–,which enhances the adsorption of DS–on the Mg surface.Further,the corrosion conducted by ion or molecule,such as Cl?or H2O,can be inhibited by the coated surfactant film This interaction provides ionic characteristics since the value ofELUMO_DS-–EHOMO_Mgis positive.As a reference,a strong covalent band can be generated only if this value is approximately zero [46,47].ΔEgapis an important parameter to express the reactivity of inhibitor molecules towards the adsorption on metallic surfaces,where a smaller value of ΔEgapimplies higher inhibitor efficien y[44].A magnitude of 3.8238 eV describes a high adsorption reactivity [45].

    3.7.2.MD consideration

    MD calculations were further implemented to provide additional proof regarding adsorption behavior.The initial state of the simulation box is shown in Fig.12a.The sulfonic acid head is deliberately initialized in the solvent farther away from the Mg surface while the alkyl tail is placed closer to the surface.The top view and main view of the equilibrium system after MD simulation are illustrated in Fig.12b and c,respectively.It can be clearly observed that the DS–adsorbs on the surface tightly,especially the sulfonic acid head.As such,it can be deduced that the surfactant DS–“grows”gradually on the metal surface so that a hydrophobic fil is generated,resulting in corrosion inhibition.The adsorption energyEadscan be calculated by Eq.(5) [24]:

    WhereEtotalis the total energy of the simulated system,is the total energy of the Mg crystal with all H2O molecules,denotes the total energy of DS–and H2O,andis equal to the total energy of H2O.Generally,the higher the absolute value ofEads,the stronger the adsorption that can be obtained by the inhibitor molecule.The calculated value in this work is ?29.0565 kJ/mol,which shows that the DS–spontaneously adsorbed on the Mg surface [46].Overall,both the DFT and MD studies are in good agreement with experimental results.

    4.Discussions

    Micro-galvanic corrosion,as a common mode of corrosion failure for alloys,can induce or even accelerate the occurrence of stress corrosion,crevice corrosion,hydrogen embrittlement and other corrosion processes [36,47].Organic inhibitors had been reported to be adsorbed on the cathodic AlMn intermetallics of Mg alloys and inhibit the cathodic reaction,leading to a reduction in the cathodic current of hydrogen evolution [19].However,in this work,some coarse Al-Li cathode phases were still exposed on the Mg alloy surface after immersed for 24 h with the addition of 0.06 M SDS,while the surface ofα-Mg matrix was covered by SDS adsorption layer with the thickness of ~200 nm.Moreover,the potential value of the SDS adsorption layer is about 10 mV,which is far less than the potential difference between the cathodic phase andα-Mg.The following is a discussion on the influenc mechanism of corrosion inhibitor SDS on micro-galvanic corrosion,and the schematic diagram is shown in Fig.13.

    Fig.12.3D adsorption configuratio of DS– on the Mg (0 0 0 1) surface of (a) initial state;(b) top and (c) main views of equilibrium state.

    Fig.13.Schematic diagram of the corrosion and corrosion inhibition mechanism of the Mg-8Li-3Al alloy: (a) and (b) without the addition of SDS;and (c)and (d) with the addition of SDS.

    In the matrix alloy,micro-galvanic couples (Al-Li/α-Mg,β-Li/α-Mg) were observed in Fig.13a,because different secondary phases have different potentials.It can be seen that hydrogen is produced on the surface of the tested alloy and there is electron exchange between the different phases.Due to the large potential difference,the Al-Li particle is the main cathode phase.As the corrosion process continued,the corrosion around the main cathode phase was serious,and the main cathode phase fell off from the alloy.As the corrosion time increased,corrosion spread along the secondary cathode(β-Li) towards the alloy interior,thus forming large pits,as shown in Fig.13b.

    When adding SDS corrosion inhibitors,SDS is spontaneously adsorbed on the alloy surface.The polar sulfate head group of SDS adsorbed with theα-Mg matrix,because the adsorption energyEadsof SDS on the Mg (0 0 0 1) surface is?29.0565 kJ/mol.While the hydrophobic 12-carbon chain tail of SDS was floate on theα-Mg matrix surface,thus form a dense and uniform SDS adsorbed layer with a thickness of about 200 nm,as shown in Figs.10b and 13c.The principal cathode Al-Li and secondary cathodeβ-Li phases were covered by the SDS-adsorbed layer,leading to a reduced potential difference between the cathode phase and anodeα-Mg and inhibited micro-galvanic corrosion with the addition of SDS(in Fig.13c).Although some coarse Al-Li cathode phases were still exposed on the alloy surface with SDS adsorption layer (as shown in Figs.8 and 13c),the surrounding area of Al-Li particles was adequately covered by the SDS adsorption layer and leading to less micro-galvanic corrosion.Therefore,it can be inferred that the formation of the dense,uniform SDS adsorbed layer effectively achieves isolation between the cathode phases,anode phase and electrolyte solution.

    As shown in Figs.10 and 13d,when the rift was formed in the SDS-adsorbed layer,the cathode phases and anodeα-Mg in the rift would soak in electrolyte solution,thus causing the formation of micro-galvanic corrosion.With increased of the corrosion time,corrosion pits more easily formed on the alloy surface.Meanwhile,a few severe corrosion sites occurred at the rift and continually expanded into the alloy.Hence,it can be concluded that the formation of the rift may be the main reason as to why the highest corrosion inhibition efficien y of the SDS-adsorbed layer on reached 89.1%.

    5.Conclusion

    The inhibition behaviors and mechanism of SDS for Mg-8Li-3Al alloy in aqueous NaCl solution have been investigated experimentally and theoretically,and the results may be summarized as follows:

    (1) The combined results of immersion testing,potentiodynamic curves and electrochemical impedance spectroscopy revealed that the Mg-8Li-3Al alloy in the NaCl+0.06 M SDS corrosive medium presented the best inhibition efficien y of 89.1%.Additionally,a uniform corrosion layer was observed on the crosssectional SEM micrograph of the Mg-8Li-3Al alloy with 0.06 M SDS addition,whereas severe localized corrosion was found to extend into the interior of the Mg-8Li-3Al alloy.

    (2) A 200-nm-thick SDS adsorption layer was observed by SEM and SKPFM,which evenly covered the surface of the Mg-8Li-3Al alloy.Meanwhile,the potential value of the SDS layer was only about 10 mV,suggesting that the potential difference between the secondary phases and Mg matrix dramatically decreased,thus inhibiting the occurrence of surface micro-galvanic corrosion.

    (3) The calculatedELUMO,EHOMO,separation energy ΔEgap,as well as theLUMO–HOMOgaps for the Mg-DS–interaction have been given,and the MD calculation indicated that DS–spontaneously adsorbed on the Mg surface due to the adsorption energyEadsvalue(?29.0565 kJ/mol).

    (4) The formation of the dense,uniform SDS-adsorbed layer effectively achieved isolation between the cathode phases,anode phase and electrolyte solution.Also,the cathode Al-Li and anodeα-Mg in the rift on the SDS adsorbed layer caused to formation of micro-galvanic corrosion and continually expanded into the alloy.

    Acknowledgments

    The authors would like to acknowledge the financia support by the National Natural Science Foundation of China(51961026)and the Interdisciplinary Innovation Fund of Nanchang University (Project No.2019-9166-27060003).

    a级毛片黄视频| 中文字幕制服av| 亚洲少妇的诱惑av| 女人久久www免费人成看片| 男女边摸边吃奶| 久久久久国产一级毛片高清牌| 精品国产国语对白av| 老汉色∧v一级毛片| 国产不卡av网站在线观看| 亚洲一码二码三码区别大吗| 国产精品免费大片| 国产片特级美女逼逼视频| 免费黄频网站在线观看国产| 国产av精品麻豆| 99热国产这里只有精品6| 我的亚洲天堂| 欧美日韩亚洲高清精品| 看免费成人av毛片| 日本色播在线视频| 高清不卡的av网站| 亚洲精品一区蜜桃| 美女脱内裤让男人舔精品视频| av福利片在线| 一区二区三区四区激情视频| 丝瓜视频免费看黄片| 亚洲国产日韩一区二区| 亚洲色图 男人天堂 中文字幕| 国产视频首页在线观看| 青春草国产在线视频| 久久久精品免费免费高清| 国产免费福利视频在线观看| 国产精品人妻久久久影院| 午夜福利影视在线免费观看| 欧美在线黄色| 欧美日本中文国产一区发布| 亚洲国产欧美一区二区综合| 亚洲欧洲日产国产| 久久久久国产精品人妻一区二区| 青春草视频在线免费观看| 精品免费久久久久久久清纯 | 亚洲精品久久午夜乱码| 自线自在国产av| 国产精品.久久久| 免费观看a级毛片全部| 两个人看的免费小视频| 亚洲五月色婷婷综合| 久久精品aⅴ一区二区三区四区| 最近最新中文字幕大全免费视频 | 热99国产精品久久久久久7| 欧美另类一区| 一本大道久久a久久精品| 免费黄频网站在线观看国产| 国产熟女欧美一区二区| 2018国产大陆天天弄谢| 侵犯人妻中文字幕一二三四区| 国产精品麻豆人妻色哟哟久久| 欧美变态另类bdsm刘玥| 国产xxxxx性猛交| 人人妻人人添人人爽欧美一区卜| 青青草视频在线视频观看| 国语对白做爰xxxⅹ性视频网站| 日日啪夜夜爽| 黄频高清免费视频| 伦理电影免费视频| 妹子高潮喷水视频| 看十八女毛片水多多多| 各种免费的搞黄视频| 王馨瑶露胸无遮挡在线观看| 中文字幕亚洲精品专区| 捣出白浆h1v1| 亚洲伊人久久精品综合| 女性被躁到高潮视频| 欧美在线黄色| 亚洲国产欧美一区二区综合| 午夜免费鲁丝| 新久久久久国产一级毛片| 99国产综合亚洲精品| 色婷婷av一区二区三区视频| 在现免费观看毛片| 尾随美女入室| 精品国产一区二区久久| 最近手机中文字幕大全| 欧美黑人精品巨大| 婷婷色麻豆天堂久久| 欧美精品av麻豆av| 另类精品久久| 久久久国产精品麻豆| 亚洲国产av新网站| 五月天丁香电影| 精品视频人人做人人爽| 国产黄色免费在线视频| 校园人妻丝袜中文字幕| 男女午夜视频在线观看| 母亲3免费完整高清在线观看| 男女边吃奶边做爰视频| 亚洲国产欧美日韩在线播放| 国产熟女欧美一区二区| 少妇人妻 视频| 亚洲免费av在线视频| 国产日韩欧美在线精品| 国产乱人偷精品视频| 在线观看一区二区三区激情| 亚洲欧美成人精品一区二区| 精品国产一区二区三区久久久樱花| 一级,二级,三级黄色视频| 青草久久国产| 十分钟在线观看高清视频www| 青春草亚洲视频在线观看| 9191精品国产免费久久| 99热全是精品| 亚洲精品一区蜜桃| 亚洲人成网站在线观看播放| 国产精品一区二区在线观看99| 成年美女黄网站色视频大全免费| 丰满乱子伦码专区| 亚洲精品国产av蜜桃| 波野结衣二区三区在线| 欧美日韩亚洲国产一区二区在线观看 | 国产欧美日韩综合在线一区二区| 热99国产精品久久久久久7| 午夜免费观看性视频| 18在线观看网站| 午夜免费鲁丝| 亚洲国产成人一精品久久久| 捣出白浆h1v1| 国产男人的电影天堂91| 婷婷色av中文字幕| videosex国产| 多毛熟女@视频| 性色av一级| www.av在线官网国产| 日本wwww免费看| 亚洲精品国产av成人精品| 国产淫语在线视频| 别揉我奶头~嗯~啊~动态视频 | 在线观看国产h片| 亚洲欧美激情在线| 成人午夜精彩视频在线观看| 天堂中文最新版在线下载| 在线观看免费日韩欧美大片| 久久久久久久大尺度免费视频| 少妇 在线观看| 中文字幕亚洲精品专区| 晚上一个人看的免费电影| 亚洲精品久久成人aⅴ小说| 久久久亚洲精品成人影院| 男女边吃奶边做爰视频| 国语对白做爰xxxⅹ性视频网站| 婷婷成人精品国产| 国产精品香港三级国产av潘金莲 | 老汉色av国产亚洲站长工具| 国产成人精品久久久久久| 69精品国产乱码久久久| 男女边吃奶边做爰视频| 日韩欧美精品免费久久| 自线自在国产av| 捣出白浆h1v1| 欧美精品一区二区免费开放| 亚洲欧洲日产国产| 一区二区三区乱码不卡18| 亚洲成人一二三区av| 热99久久久久精品小说推荐| 满18在线观看网站| 大码成人一级视频| 久久精品国产综合久久久| 日韩制服丝袜自拍偷拍| av有码第一页| 少妇被粗大猛烈的视频| 久久青草综合色| 欧美日韩亚洲高清精品| 毛片一级片免费看久久久久| 99久久综合免费| 久久久精品区二区三区| 秋霞伦理黄片| 男女之事视频高清在线观看 | 精品国产超薄肉色丝袜足j| 丝袜喷水一区| 性高湖久久久久久久久免费观看| 免费高清在线观看视频在线观看| 美女扒开内裤让男人捅视频| 亚洲成人手机| 亚洲图色成人| 高清视频免费观看一区二区| 老鸭窝网址在线观看| 日韩中文字幕视频在线看片| 大片免费播放器 马上看| 999精品在线视频| 精品一区二区三卡| 国产免费福利视频在线观看| 国产成人免费观看mmmm| 亚洲精品久久午夜乱码| 丰满饥渴人妻一区二区三| 一本一本久久a久久精品综合妖精| 国产福利在线免费观看视频| 中文字幕色久视频| 最近手机中文字幕大全| av线在线观看网站| 亚洲国产日韩一区二区| 十八禁人妻一区二区| 激情五月婷婷亚洲| 超碰成人久久| tube8黄色片| 欧美日韩亚洲综合一区二区三区_| 成人国产麻豆网| 欧美日韩视频高清一区二区三区二| 免费在线观看视频国产中文字幕亚洲 | 蜜桃在线观看..| 中文天堂在线官网| 18禁裸乳无遮挡动漫免费视频| 国精品久久久久久国模美| a级毛片黄视频| 亚洲第一av免费看| 女人久久www免费人成看片| 日本午夜av视频| 日本av手机在线免费观看| 少妇 在线观看| 欧美亚洲日本最大视频资源| 深夜精品福利| 中文欧美无线码| 亚洲av电影在线观看一区二区三区| 亚洲欧美中文字幕日韩二区| 欧美在线黄色| 亚洲精品aⅴ在线观看| av网站在线播放免费| 最近手机中文字幕大全| 亚洲国产精品国产精品| 宅男免费午夜| 超色免费av| 久久ye,这里只有精品| 久久精品国产综合久久久| 男人舔女人的私密视频| 啦啦啦 在线观看视频| 国产精品成人在线| 国产视频首页在线观看| 亚洲欧美中文字幕日韩二区| 巨乳人妻的诱惑在线观看| 国产av精品麻豆| 成人国产av品久久久| 亚洲精品国产av成人精品| 国产成人一区二区在线| 啦啦啦在线观看免费高清www| 国产成人欧美| 十八禁人妻一区二区| 日韩免费高清中文字幕av| 99re6热这里在线精品视频| 国产男人的电影天堂91| 精品亚洲成国产av| 男人操女人黄网站| 麻豆精品久久久久久蜜桃| 男女之事视频高清在线观看 | 国产精品久久久久成人av| 国产精品二区激情视频| 成人国产av品久久久| 欧美xxⅹ黑人| 日韩欧美精品免费久久| 性高湖久久久久久久久免费观看| 久久99热这里只频精品6学生| 日韩电影二区| 成年女人毛片免费观看观看9 | 亚洲欧洲精品一区二区精品久久久 | 咕卡用的链子| 丝瓜视频免费看黄片| 亚洲国产日韩一区二区| 视频区图区小说| 好男人视频免费观看在线| 免费在线观看完整版高清| av片东京热男人的天堂| 亚洲图色成人| 如日韩欧美国产精品一区二区三区| 美女国产高潮福利片在线看| 国产成人精品久久久久久| 毛片一级片免费看久久久久| 19禁男女啪啪无遮挡网站| 黑人猛操日本美女一级片| 国产亚洲av片在线观看秒播厂| 欧美日韩一级在线毛片| 亚洲av综合色区一区| 中文天堂在线官网| 精品久久久精品久久久| 久久女婷五月综合色啪小说| 人人妻人人添人人爽欧美一区卜| 欧美黄色片欧美黄色片| 久久久久人妻精品一区果冻| 久久婷婷青草| 亚洲综合色网址| av线在线观看网站| 美女主播在线视频| 天堂俺去俺来也www色官网| 亚洲精品aⅴ在线观看| 一区二区三区乱码不卡18| 欧美激情极品国产一区二区三区| 91老司机精品| 国产精品国产av在线观看| 国产精品免费视频内射| 在线天堂中文资源库| 丰满乱子伦码专区| 欧美变态另类bdsm刘玥| 夫妻午夜视频| 久热这里只有精品99| 成年女人毛片免费观看观看9 | 狠狠精品人妻久久久久久综合| 久久久精品免费免费高清| 久久久久久人人人人人| 亚洲男人天堂网一区| 国产一区二区三区av在线| 久久久久久久精品精品| 新久久久久国产一级毛片| 国产精品久久久久久人妻精品电影 | 欧美人与性动交α欧美精品济南到| 欧美乱码精品一区二区三区| 精品国产一区二区三区四区第35| 香蕉国产在线看| 久久精品国产综合久久久| 午夜福利一区二区在线看| av免费观看日本| 老司机深夜福利视频在线观看 | 夜夜骑夜夜射夜夜干| 另类亚洲欧美激情| 黑人猛操日本美女一级片| 国产xxxxx性猛交| 成人国语在线视频| 欧美av亚洲av综合av国产av | 久久精品aⅴ一区二区三区四区| 久久久久视频综合| 国产又爽黄色视频| 在线观看三级黄色| 日韩 欧美 亚洲 中文字幕| 亚洲成国产人片在线观看| 国产乱人偷精品视频| 免费在线观看黄色视频的| 婷婷色av中文字幕| 国产熟女午夜一区二区三区| 女人被躁到高潮嗷嗷叫费观| 婷婷色综合大香蕉| 国产欧美日韩一区二区三区在线| 别揉我奶头~嗯~啊~动态视频 | 国产精品免费视频内射| 亚洲精品自拍成人| 麻豆乱淫一区二区| 久久99精品国语久久久| 另类精品久久| 美女主播在线视频| 自线自在国产av| 久久精品久久久久久噜噜老黄| 国产精品久久久久久久久免| 在线观看免费午夜福利视频| 日韩制服丝袜自拍偷拍| 色婷婷av一区二区三区视频| 亚洲成人免费av在线播放| 99久久精品国产亚洲精品| 99久久人妻综合| 精品国产一区二区久久| av在线app专区| 日韩制服骚丝袜av| 97在线人人人人妻| 少妇的丰满在线观看| 免费黄频网站在线观看国产| 777米奇影视久久| 午夜影院在线不卡| 美女扒开内裤让男人捅视频| 另类亚洲欧美激情| 欧美久久黑人一区二区| 欧美日本中文国产一区发布| 美女福利国产在线| 人体艺术视频欧美日本| 我要看黄色一级片免费的| 免费观看人在逋| 一级毛片我不卡| 色网站视频免费| 女人爽到高潮嗷嗷叫在线视频| 亚洲av电影在线进入| 中文字幕高清在线视频| 国产熟女欧美一区二区| 国产老妇伦熟女老妇高清| 国产成人精品久久久久久| 热99国产精品久久久久久7| 欧美xxⅹ黑人| 国产 精品1| 久久久久精品人妻al黑| 久久韩国三级中文字幕| 国产有黄有色有爽视频| 国产精品熟女久久久久浪| 精品免费久久久久久久清纯 | h视频一区二区三区| 国产激情久久老熟女| avwww免费| 国产精品三级大全| 久久天堂一区二区三区四区| 成人漫画全彩无遮挡| 熟女av电影| 亚洲国产精品一区三区| 亚洲精品美女久久av网站| 高清不卡的av网站| 久久久久久久久久久免费av| 美国免费a级毛片| 狠狠婷婷综合久久久久久88av| 免费在线观看视频国产中文字幕亚洲 | 黑人巨大精品欧美一区二区蜜桃| 七月丁香在线播放| 国产日韩欧美在线精品| 激情五月婷婷亚洲| 亚洲在久久综合| 日韩成人av中文字幕在线观看| 久久久久精品人妻al黑| 王馨瑶露胸无遮挡在线观看| 国产有黄有色有爽视频| 亚洲欧美中文字幕日韩二区| 男女边摸边吃奶| 韩国高清视频一区二区三区| 精品国产乱码久久久久久小说| 黑人巨大精品欧美一区二区蜜桃| 爱豆传媒免费全集在线观看| 精品一品国产午夜福利视频| 亚洲国产精品一区三区| 伦理电影大哥的女人| 国产成人a∨麻豆精品| 亚洲av日韩精品久久久久久密 | 国产 精品1| 日韩大码丰满熟妇| 欧美成人精品欧美一级黄| 亚洲欧洲国产日韩| 国产精品久久久av美女十八| √禁漫天堂资源中文www| 亚洲精品一二三| xxxhd国产人妻xxx| 人体艺术视频欧美日本| 黄色视频在线播放观看不卡| 青草久久国产| 国产深夜福利视频在线观看| 在线观看三级黄色| 国产一卡二卡三卡精品 | 久久毛片免费看一区二区三区| 国产人伦9x9x在线观看| 欧美激情极品国产一区二区三区| 欧美乱码精品一区二区三区| 男的添女的下面高潮视频| 在线免费观看不下载黄p国产| 欧美在线一区亚洲| 久久影院123| 精品第一国产精品| 在线观看一区二区三区激情| 最近最新中文字幕免费大全7| 九色亚洲精品在线播放| 欧美激情极品国产一区二区三区| 在线观看www视频免费| 人人妻人人澡人人爽人人夜夜| 青青草视频在线视频观看| 美女中出高潮动态图| 女人精品久久久久毛片| 中文字幕制服av| 国产精品成人在线| 久久久久国产精品人妻一区二区| 精品免费久久久久久久清纯 | 欧美最新免费一区二区三区| 免费日韩欧美在线观看| 免费高清在线观看视频在线观看| 亚洲专区中文字幕在线 | 欧美亚洲日本最大视频资源| 一级,二级,三级黄色视频| 中文天堂在线官网| 狂野欧美激情性xxxx| 亚洲成人一二三区av| 国产熟女欧美一区二区| 国产免费视频播放在线视频| 国产免费福利视频在线观看| 天堂俺去俺来也www色官网| 又大又黄又爽视频免费| 亚洲国产毛片av蜜桃av| 久久久久精品国产欧美久久久 | 老司机深夜福利视频在线观看 | 精品国产乱码久久久久久男人| 久久精品亚洲熟妇少妇任你| 亚洲av中文av极速乱| 丝袜在线中文字幕| 久久久久视频综合| 久久亚洲国产成人精品v| 一二三四中文在线观看免费高清| 在线看a的网站| 亚洲熟女毛片儿| 人人妻人人澡人人看| 巨乳人妻的诱惑在线观看| 亚洲av电影在线观看一区二区三区| 国产精品.久久久| 不卡av一区二区三区| 国产精品国产三级专区第一集| 国产爽快片一区二区三区| av又黄又爽大尺度在线免费看| 精品国产一区二区三区四区第35| 两个人看的免费小视频| 大陆偷拍与自拍| av女优亚洲男人天堂| 91aial.com中文字幕在线观看| 女人高潮潮喷娇喘18禁视频| www.熟女人妻精品国产| 女性生殖器流出的白浆| 曰老女人黄片| 哪个播放器可以免费观看大片| 亚洲精品一区蜜桃| 精品一区在线观看国产| 人人妻人人澡人人爽人人夜夜| 国产精品久久久久久精品电影小说| 如日韩欧美国产精品一区二区三区| 新久久久久国产一级毛片| www.熟女人妻精品国产| 国产在线免费精品| 国产片特级美女逼逼视频| 亚洲色图 男人天堂 中文字幕| 亚洲成人免费av在线播放| 亚洲精品国产色婷婷电影| 成人免费观看视频高清| 九色亚洲精品在线播放| 少妇人妻 视频| 在线观看免费视频网站a站| xxxhd国产人妻xxx| 男女免费视频国产| 免费看av在线观看网站| 日韩大片免费观看网站| 亚洲男人天堂网一区| 亚洲综合精品二区| 国产精品偷伦视频观看了| 女人高潮潮喷娇喘18禁视频| 久久午夜综合久久蜜桃| videosex国产| 一二三四在线观看免费中文在| 日韩精品有码人妻一区| 一区二区日韩欧美中文字幕| 久久国产亚洲av麻豆专区| 在线看a的网站| 飞空精品影院首页| 777米奇影视久久| 老汉色av国产亚洲站长工具| av不卡在线播放| 制服诱惑二区| 亚洲第一av免费看| 啦啦啦中文免费视频观看日本| 在线看a的网站| 女的被弄到高潮叫床怎么办| 色综合欧美亚洲国产小说| 精品人妻熟女毛片av久久网站| 热99国产精品久久久久久7| 成人黄色视频免费在线看| 国产xxxxx性猛交| 看免费av毛片| 成人毛片60女人毛片免费| 国产欧美日韩一区二区三区在线| 国产成人啪精品午夜网站| 亚洲精品美女久久av网站| a级毛片黄视频| 久久久久精品久久久久真实原创| 考比视频在线观看| 丰满少妇做爰视频| 黄片无遮挡物在线观看| videos熟女内射| 黑人欧美特级aaaaaa片| 在线天堂最新版资源| 亚洲成人一二三区av| 色综合欧美亚洲国产小说| 19禁男女啪啪无遮挡网站| 国产 精品1| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品一二三| 一本大道久久a久久精品| 国产 一区精品| 巨乳人妻的诱惑在线观看| 亚洲精品av麻豆狂野| av网站免费在线观看视频| 亚洲精品国产av蜜桃| 男女国产视频网站| 香蕉丝袜av| 国产精品久久久久久久久免| 亚洲国产中文字幕在线视频| 欧美最新免费一区二区三区| 激情视频va一区二区三区| 亚洲av福利一区| 欧美激情极品国产一区二区三区| 一二三四在线观看免费中文在| 99热网站在线观看| 亚洲在久久综合| 丝袜脚勾引网站| 性少妇av在线| 在线精品无人区一区二区三| 五月天丁香电影| 日韩精品有码人妻一区| 亚洲精品一二三| 热99国产精品久久久久久7| 日本一区二区免费在线视频| 久久天堂一区二区三区四区| 2021少妇久久久久久久久久久| 在线观看免费午夜福利视频| 热99久久久久精品小说推荐| 最近手机中文字幕大全| 99香蕉大伊视频| 久久毛片免费看一区二区三区| 高清欧美精品videossex| 亚洲熟女精品中文字幕| 欧美黑人精品巨大| 国产精品秋霞免费鲁丝片| 亚洲人成网站在线观看播放| 最近手机中文字幕大全| 美女主播在线视频| 国产精品久久久人人做人人爽| 无限看片的www在线观看| 欧美日韩一级在线毛片| 亚洲视频免费观看视频| av不卡在线播放| 亚洲av成人精品一二三区| a级片在线免费高清观看视频| 热99国产精品久久久久久7| 一区二区av电影网| 中文字幕人妻丝袜一区二区 | 91aial.com中文字幕在线观看| 一区福利在线观看| 亚洲精品国产av成人精品| 哪个播放器可以免费观看大片| 赤兔流量卡办理|