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

    Microstructure and surface texture driven improvement in in-vitro response of laser surface processed AZ31B magnesium alloy

    2021-10-30 12:49:22TsoChngWuSmeehnJoshiYeeHsienHoMngeshPntwneSubhsisSinhNrendrDhotre
    Journal of Magnesium and Alloys 2021年4期

    Tso-Chng Wu ,Smeehn S.Joshi,b ,Yee-Hsien Ho,b ,Mngesh V.Pntwne,b ,Subhsis Sinh,c ,Nrendr B.Dhotre,b,*

    aDepartment of Materials Science and Engineering,University of North Texas,3940 N Elm St,Denton,TX 76207,USA

    b Center for Agile and Adaptive Additive Manufacturing,University of North Texas,3940 N Elm St,Denton,TX 76207,USA

    c Department of Metallurgical Engineering,Indian Institute of Technology (BHU),Varanasi,Uttar Pradesh 221005,India

    Abstract The present work explored effects of laser surface melting on microstructure and surface topography evolution in AZ31B magnesium alloy.Thermokinetic effects experienced by the material during laser surface melting were simulated using a multiphysics finit element model.Microstructure and phase evolution were examined using scanning electron microscopy,X-ray diffraction,and electron back scatter diffraction.Surface topography was evaluated using white light interferometry.The interaction of surface melted samples with simulated body flui was monitored by contact angle measurements and immersion studies up to 7 days.Laser surface melting led to formation of a refine microstructure with predominantly basal crystallographic texture.Concurrently,the amount of β phase (Mg17Al12) increased with an increase in the laser fluence β phase preferentially decorated the cell boundaries.In terms of topography,the surface became progressively rougher with an increase in laser fluence As a result,upon immersion in simulated body fluid the laser surface melted samples showed an improved wettability,corrosion resistance,and precipitation of mineral having composition closer to the hydroxyapatite bone mineral compared to the untreated sample.

    Keywords: Laser surface melting;Magnesium alloy;Laser surface engineering;Biomineralization.

    1.Introduction

    Biodegradable magnesium (Mg) alloys have gained considerable interest for consumable implant application due to their appropriate properties,such as low density (1.74 g/cm3),elastic modulus similar to natural bone,and biocompatibility[1-6].Recently,Mg alloys have been developed for several implant devices such as cardiovascular stents,bone fixation and porous scaffolds [7,8].Although,biodegradable Mg alloys are potentially attractive for the orthopedic field there are number of issues which limit the development of Mg based biomaterials.The primary issue with Mg based materials is their propensity to accelerated corrosion,especially in chloride ion containing physiological environment [9-15],resulting in the deterioration of mechanical properties of implant before the host bone tissue has been sufficientl healed.Furthermore,the evolved hydrogen bubbles during corrosion reaction can form gas pockets in the vicinity of the implant,leading to separation of tissue layer [16].Moreover,the released Mg ions dissolved into the surrounding blood plasma can induce localized alkalization near the implant leading to damage of the cells and the regenerated tissue [17].These phenomena become a challenge and demand quite an effort to resolve the associated drawbacks in order to increase the success rate of implantation.Hence,it is critical to reduce and control the degradation rate of Mg based implant in physiological environment.

    It has been reported that the corrosion resistance of engineering alloys,including Mg alloys,is significantl influence by the microstructural features,such as grain size,crystallographic texture (orientation),and chemical composition [18].Besides the improved corrosion properties,a favorable surface (physical texture/roughness) to enhance cell interaction in physiological environment is also essential for a successful implant [19].In order to efficientl reduce degradation rate and improve interfacial biocompatibility of Mg based implants,various surface modification have been developed.The methods of surface modificatio include chemical modification such as plasma electrolytic oxidation [20] as well as microarc oxidation [21],and physical modification such as heat treatment [22],friction stir processing [23-25],and laser processing [26-29].Among these methods,laser offers a rapid and reproducible processing solution for the implant material.Laser provides a non-contact and clean high energy density source by focusing high power concentrated into a small diameter (few microns) spot on material surface.Laser processing can treat site specifi surfaces to have a spatial control on the surface properties without disturbing the bulk material.In addition,due to use of lasers such as Nd-YAG laser with an optical fibe based beam delivery,laser processing is highly compatible with automation and remote processing in intended manufacturing environment and can treat the contoured surfaces further enhancing these advantages.Laser beam generates a small-scale localized interaction zone with the material (few hundredμm depth) which leads to a self-quenching effect (104-105K/s) occurring due to rapid heat extraction from the bulk of the substrate producing a refine microstructure [30,31].These thermophysical effects can be employed for various end effects such as material heat treatments via phase transitions,surface melting (LSM)for producing physical texture,and material removal for machining [4,31-33].By controlling various parameters,such as power,scanning rate,and laser track overlap;the process can be optimized to produce a microstructure and suitable surface roughness on the implants to improve the biocompatibility and biointegration characteristics [27,34].Several studies revealed the refine grain structure within Mg alloys surfaces after laser surface melting (LSM) [27,28,35-37].Liu et al.used 10 kW CO2laser to carry out surface melting of AM60B Mg alloy and reported the homogeneous distribution of secondary-phase with refineα-Mg grains[35].In previous studies from the present research group [27,28],the optimal laser fluenc was investigated to produce enriched Mg17Al12βphase at theα-Mg grain/cell boundaries in AZ31B Mg alloy for the enhanced corrosion resistance.Rakesh et al.studied[38] the positive correlation between surface roughness and wettability of laser processed Mg-Zn-Dy alloy.These find ings from the literature indicated that laser process provides multiple advantages to produce Mg based implant with suffi cient corrosion resistance and biocompatibility.

    Even though the previous studies pointed towards potential of LSM for Mg based implant materials,there were gaps such as understanding about evolution of crystallographic texture and exact nature of mineral deposited upon immersion in simulated body flui (SBF).As a continuation of the investigation and development of Mg based bioimplant materials by the present research group [24,25,27,28,39],the current work focused on probing into the thermokinetic effects and the influenc of laser processing on microstructure evolution in AZ31B Mg alloy using an integrated experimental and computation approach.The thermal profile experienced by the material were simulated using a multiphysics computational model.The phase and microstructure evolution were probed in the laser surface melted zone using analytical techniques to get an idea about microstructure and phase evolution,development of crystallographic texture,and surface roughness in the LSM samples.LSM samples were subjected toin-vitroevaluation in the SBF using contact angle measurements and immersion studies.Post immersed samples were observed using optical and electron microscopy.Furthermore,to get an idea about precipitated mineral compounds upon immersion,X-ray photo electron spectroscopy (XPS) was utilized.The present study shed a light on combined effects of various aspects associated with the laser processing to deepen the understanding about utilizing laser as an effective processing tool to produce improved Mg based bioimplant materials.

    2.Experimental procedures

    2.1.Materials and sample preparation

    Commercially available AZ31B Mg alloy (composition:3.0 wt% Al,1.0 wt% Zn,0.5 wt% Mn,and Mg balance) was used for this study.Prior to laser surface modification the AZ31B Mg sheets were cut using slow speed diamond saw into the rectangular blocks of dimension 50 mm × 50 mm ×7 mm (width × length × thickness) (Fig.1).All the coupons were then ground using silicon carbide (SiC) abrasive papers from 60 grit to 400 grit to remove the oxide layer and contamination.The coupons were then cleaned by rinsing with water and methanol.

    2.2.Laser surface melting (LSM)

    A 3 kW diode pumped Ytterbium fibe laser in continuous mode (IPG Photonics,model# YLS-3000,wavelength of 1064 nm) was employed during this study to conduct LSM.The laser beam had a Gaussian profil (TEM00mode) with a diameter of 0.6 mm on the sample surface.The focal length of the laser was 120 mm.Such a long laser focal length ensured minimal/no deviation in spot size which remained fully focused on the sample surface.The laser scanning speed of 500 mm/s and laser track fil spacing (distance between centers of consecutive neighboring track) of 0.15 mm were kept constant during all the laser processing experiments while the laser power was varied in the range of 250-750 W (Table 1).These parameters were chosen based on several trials and optimization studies previously published by the present research group [27,28,39].Ar(g)was used as the shielding gas during all the laser processing experiments.The processed samples were then be used for microstructure characterization,contact angle measurements and immersion studies in SBF.

    Fig.1.Schematic of laser surface melting (LSM) process.

    Table 1 LSM parameters used to process AZ31B Mg alloy.

    2.3.Phase identificatio

    Phases in untreated and the LSM samples were identifie using X-ray diffraction (XRD) technique.A Rigaku Ultima diffractometer with Cu Kαradiation (wavelength -0.15418 nm),operating at 40 kV and 40 mA was employed to carry out XRD analysis.Samples were scanned in the 2θrange of 20?to 90?using a step size of 0.025 and a scan rate of 1.5?/mm.Phase identificatio was carried out by comparing the XRD spectra with standard International Center for Diffraction Data (ICDD) file obtained from the Joint Committee of Powder Diffraction Standards (JCPDS).

    2.4.Microstructure characterization

    Microstructure examination was conducted in the crosssection perpendicular to the laser track and laser treated surface.The cross-sections of samples were cut into smaller rectangular blocks (5 mm × 3 mm × 7 mm) using a slow speed diamond wafer cutting machine with oil based lubricant to achieve stress free cut.The cut samples were mounted in epoxy molds.These mounted samples were ground using series of (120-1200 grit) SiC abrasive papers,followed by fin surface polishing in 1-0.03 μm alumina (Al2O3) water-based abrasive suspension.A smooth mirror finis on the polished surface was obtained by employing 0.02 μm colloidal silica(SiO2) suspension for 40 min.The finishe polished samples were then cleaned by ultrasonic setup in succession with deionized water and methanol as the media.Prior to the microstructural investigation,the prepared samples were etched by the acetic picric solution (5 ml acetic acid+6 g picric acid+100 ml ethyl alcohol+10 mL H2O) to denote the microstructural features.The etched coupons were characterized using scanning electron microscope (SEM) equipped with Energy-Dispersive X-ray Spectroscopy (EDS) (FEI Quanta 200 environmental SEM).The crystallographic orientation of laser treated samples on the surface and in the cross-section were observed and analyzed by using FEI Nova NanoSEM 230 equipped with a Hikari super electron backscatter diffraction (EBSD) detector (operated at 20 kV).The EBSD scanning step size was 0.5 μm.TSL OIM Analysis 8.0 software was used to obtain microstructure maps and pole figure from the EBSD data.

    2.5.Numerical modeling

    The main purpose of implementing an integrated experimental and computational approach during the current work was to realize the microstructure evolution in LSM AZ31B alloy and its effect on the response of the LSM AZ31B Mg alloy in SBF environment.Laser is a rapid process characterized by μs-ms laser beam-material interaction times and small laser spot sizes in the range of μm-mm.Complex physical phenomena such as plasma and plume formation occur during these time frames in a small region making it difficul to monitor the laser beam-material interaction and carry out experimental measurements of temperature evolution experienced by the material.Material rapidly undergoes solid state heating,melting,vaporization,and sever cooling within the short time frames (μs-ms) which leads to a refine and novel microstructure.In light of this,authors modeled the LSM process using a multiphysics process model developed on a COMSOLT Mplatform to gather better understanding about the rapid time-temperature cycles experienced by AZ31B material while undergoing LSM treatment.The model mainly dwells on thermokinetics involved during laser-material interaction while considering the complex hydrodynamics associated with the melt pool.The development of the model detailing the governing equations,assumptions,boundary conditions,meshing,and model configuratio is adopted from the earlier reports by the authors [27,40].

    2.6.Topographic characterization

    The surface roughness of untreated and LSM AZ31B Mg samples were measured with 10X objective using white light by the Rtec Universal Tribometer (Rtec Instrument,San Jose,CA,USA).The detected surface area was 1.3 mm× 1.6 mm.The surface profil data was fi as per the standard protocol according to the International Organization for Standardization [41].Advanced data processing was conducted by a modular program software,Gwyddion 2.50.As a result,sample surface roughness was presented in form of the arithmetic mean height (Ra,μm),the root mean square height (Rz,μm),as well as the topographic maps.

    2.7.Wettability

    The surface wettability of untreated Mg and LSM samples was evaluated by a static liquid sessile drop method using CAM-Plus? contact angle goniometer (Cheminstruments Inc.Fairfiled OH),equipped with fibe optic light source.Freshly prepared SBF was used as testing medium,and the preparation process has been described in the previous publication [39].A liquid droplet of volume 3 μL (drop diameter of~1 mm) was placed on the thoroughly cleaned sample by a hypodermic syringe.The liquid droplet was placed on the surface for approximately 10 s to stabilize before the reading was taken.The measurement was carried out at room temperature (25.0?C) and a minimum of ten contact angle readings were taken on each sample at various locations to minimize errors in the measurement.In order to eliminate the error associated with the arbitrary tangential alignment,the contact angle was measured based on the patented half angle method(US Patent 5268733).

    2.8.Biomineralization evaluation

    To examine the biomineralization behavior,the untreated and LSM AZ31B Mg samples were cut into coupon of 5 mm× 5 mm and mounted in epoxy.The cut coupons were immersed in SBF maintained at a constant temperature of 36.5?C with the aid of constant temperature water bath equipped with a temperature sensor and digital temperature read out.The samples were placed vertically in a 250 ml beaker with SBF.The volume of SBF for the immersion study was determined by Vs=Sa/10,where Vsis the volume of SBF (mL) and Sais the surface area of the coupon (mm2).The periods of immersion time were setup as 1 day,3 days,5 days,and 7 days.In order to maintain a pH value of the SBF solution at 7.4,flui was refreshed every 24 h.After immersion,the samples were removed from SBF solution and rinsed with distilled water and then dried at room temperature.The morphological and elemental characterization of the surface mineralization was conducted by SEM and EDS,and the chemical compound analysis was conducted using Xray photoelectron spectroscopy (XPS) (make:VersaProbeTM5000).The monochromatic X-ray beam source at 1486.6 eV,49.3 W,and 200 μm beam diameter was operated in the vacuum (5 x 10-6Pa) to perform scans on the mineralized surface for the analysis.

    2.9.Corrosion behavior

    The corrosion behavior of untreated and LSM AZ31B Mg samples was studied using weight loss measurement.It is worth mentioning here that the electrochemical corrosion behavior was investigated in detail in previous report published by the current research group[28].The experimental setup for weight loss measurement was the same as that used in biomineralization evaluation described above.However,in this case,post-experimental sample cleaning process was different from that in the biomineralization testing.The immersed samples were firs cleaned with boiling chromic acid (180 g/l) to remove mineralization products and surface corrosion layer.Then rinsing with adequate methyl alcohol was performed to wash out the residual acid followed by air drying at room temperature.The weight loss is deduced as per the following Eq.(1) [42,43]:

    According to the measured weight loss (W),the corrosion rate can be further derived as expressed in Eq.(2):

    WhereApresents the original surface area of the sample(cm2),Tis immersion time (h),andDis the density of the sample (g/cm3).

    3.Results and discussion

    3.1.Phase identificatio

    Fig.2.X-ray diffraction spectra for untreated and LSM AZ31B Mg samples.

    Phase analysis conducted using XRD revealed the peaks corresponding to solid solution and intermetallic phases in untreated and LSM AZ31B Mg Samples (Fig.2).Theα-Mg phase with hexagonal close packed (HCP) structure was dominant phase whereas minor peaks corresponding to an intermetallicβphase (Mg17Al12) were also observed in both untreated and LSM AZ31B (Fig.2).Theβphase peaks with higher intensity were observed among the laser processed samples.With an increase in the laser fluence there was a progressive increase in the intensity of peaks associated withβphase as reported in previous papers by the authors[27,28].Additionally,the peaks corresponding to theα-Mg phase were shifted toward higher angles by~0.2?for LSM samples.Furthermore,an increase in peak width was observed for the laser treated samples suggesting a refinemen of the microstructure[44].The change in the intensity ofβphase was likely to be an outcome of changing volume fraction ofβphase,which,in turn,can be due to changing composition and thermokinetics of the melt pool with varying laser fluence This increase in volume fraction ofβ-Mg phase was presented in previous publications by the authors along with the detailed methodology of the calculations [27,28,45].The calculated volume fraction ofβphase did increase with an increase in the laser fluenc (volume percent range<0.5% to 3.3%).To further probe into these initial observations about microstructure refinemen and phase evolution,it was necessary to understand the thermokinetic effects experienced by the samples during LSM.Thus,a transient temperature distribution was computationally predicted and analyzed.

    The transient temperature computationally probed at the surface of the melt pool corresponding to the different laser fluenc varying from 1.06 to 3.18 J/mm2(250 to 750 W) at a constant scanning rate (500 mm/s) is presented in Fig.3.The thermokinetic parameters such as heating rate,cooling rate,and peak temperature obtained from these temperature-time plots are listed in Table 2.Due to increased laser fluence the heating rate can be seen to increase (0.748 to 1.952 ×105K/s)with laser fluenc (Table 2).The peak temperature of the melt pool surface was reached at the end of residence timetr=D/V=1.2 ms,where D is the laser beam diameter on the sample surface and V is the laser beam scanning velocity.As the laser fluenc increased,the peak temperature increased from 1359 K (at 1.06 J/mm2) to 2955 K (at 3.18 J/mm2).Vaporization temperatures of Mg,Al,and Zn which are the key constituent elements of AZ31B alloy,are 1364 K,2743 K,and 1180 K respectively[46].The processing condition corresponding to laser fluenc of 1.06 J/mm2led to surface melting of the alloy and reaching the temperature for very short duration of time (<0.5 ms) above the vaporization temperature Zn (Fig.3).Thus,the composition likely remained unaffected as the material experienced rapid thermokinetics driven melt pool dynamics.At the laser fluenc of 2.12 J/mm2,the peak temperature of 2192 K surpassed vaporization temperature of both Zn and Mg (Fig.3).Finally,for the laser fluenc of 3.18 J/mm2,the temperature exceeded vaporization points of Zn,Mg,and Al (Fig.3).

    Fig.3.Computationally predicted variation in time-temperature profile for all LSM AZ31B Mg samples.Vaporization temperature (Tv) for Zn,Mg,and Al is also shown.In addition,the region of mushy zone between solidus(Tsolidus) and liquidus (Tliquidus) is shown.

    Table 2 Computationally derived thermokinetic parameters.

    Fig.4.SEM BSE mode micrographs in cross sectional view (XZ plane) of (a) Untreated AZ31B Mg,(b)-(c) low magnificatio views corresponding to 3.18 J/mm2 revealing a cellular morphology and (d) high magnificatio view corresponding to 3.18 J/mm2 sample showing underlying equiaxed grain structure as a result of channeling contrast.

    To analyze the peak shifts observed in XRD spectra(Fig.2),the finding from the model were coupled with the experimental observations to better realize the phase evolution during LSM of AZ31B Mg alloy.The lattice parameters forα-Mg andβphase were calculated from the XRD spectra of LSM and untreated AZ31B Mg alloy.The values of lattice parameters were reported in Fig.2.As a firs observation,the value of lattice parameter ofβphase remained identical in all the cases.βphase is an intermetallic compound with a fi ed composition at room temperature[47].However,the lattice parameters ofα-Mg appeared to increase for 1.06 J/mm2LSM sample compared to the untreated AZ31B Mg alloy(Fig.2).The lattice parameters then steadily decreased at the laser fluenc of 2.12 J/mm2and became comparable to the untreated AZ31B Mg alloy for the sample processed with laser fluenc of 3.18 J/mm2(Fig.2).Such a variation in lattice parameter can be attributed to factors such as solute concentration of alloying elements in the matrix and micro residual stresses which can change the inter planar spacing [48].For LSM sample processed with 1.07 J/mm2,as the melting occurred with temperature exceeding vaporization of Zn for a short time (Fig.3),it is possible that density driven segregation of heavy elements such as Zn would occur towards the depth of the melt pool of liquid AZ31B Mg alloy.Simultaneously,minor quantities of Zn would vaporize from the surface layer.This would leaveα-Mg matrix lean of Zn in the surface region.It has been reported that the lattice parameters of Mg increase with a decrease in solute contents of Zn and Al [49].With an increase in laser fluenc to 2.12 J/mm2the AZ31B Mg alloy experienced temperature well above vaporization of Zn and Mg(Fig.3).This would have possibly led to higher lattice parameters forα-Mg compared to the untreated AZ31B Mg alloy[49]as shown in Fig.2.Lastly,for the LSM AZ31B Mg sample corresponding to 3.18 J/mm2laser flu ence,the temperature exceeded vaporization of Zn,Mg,and Al(Fig.3).The computed lattice parameters for this condition were similar to the untreated AZ31B Mg alloy (Fig.2).It is likely that effect of localized variation in chemistry occurring due to redistribution and evaporation of constituent elements driven by thermokinetics and complex physical phenomena such as Marangoni convection during LSM can affect the lattice parameters ofα-Mg phase,which is reflecte as the shift(~0.2?) in the peaks associated withαphase in the XRD spectra of the LSM samples (Fig.2).

    The cooling rate was mainly computationally determined in the solidificatio range (839 K-905 K) which dictates the size of the microstructural features.The highest cooling rate of 5.283 ×104K/s was computed at the laser flu ence of 1.06 J/mm2,which decreased up to 4.481×104K/s at 3.18 J/mm2.The increased size of the melt pool with increased laser fluenc lowers the rate of heat extraction by the surrounding solid mass,causing reduction in cooling rate[32,40].Such a rapid thermokinetics,in addition to the unique combination of other auxiliary parameters such as thermal gradient,and solidificatio rate up to 500 mm/s,produced distinct microstructural morphology which was further probed using SEM and EBSD.

    The SEM micrographs of untreated AZ31B mg in back scatter electron (BSE) mode mainly revealed the coarse equiaxedα-Mg grains with an average grain size of 10±1.2 μm (Fig.4(a)).On the other hand,the LSM AZ31B Mg viewed in XZ plane (Fig.4 (b)-(d)) mainly displayed the occurrence of solidificatio under cellular mode throughout the melt pool.The cellular structure was highly refine compared to the microstructure of the untreated sample (Fig.4(a)).Interestingly,even with the presence of cellular structure,the underlying grain structure was clearly visible as a result of channeling contrast when observed carefully at a higher magnificatio (Fig.4 (c)) for the LSM AZ31B Mg samples.Moreover,the EDS elemental point probe indicated the segregation of Al along the intercellular boundaries as reported in previous papers by the authors [27,28].The Al concentration in this intercellular region was in the range of 3.5 to 4 wt%,which decreased to~1 to 2 wt% in the intracellular region.Thus,βphase is likely to occur in the intercellular or aluminum rich region.As the cooling rate decreased with increased laser fluence the more aluminum likely to diffused in the intercellular region thereby increasing theβphase formation (Fig.2).Additional insights into grain morphology and crystallographic texture were gained using EBSD.

    Fig.5.EBSD results showing (a) inverse pole figur map for top surface (b) 0001 pole figur for the top surface (c) image quality map for the top surface,(d) inverse pole figur map for cross-section,(e) 0001 pole figur for cross section,and (f) image quality map for the cross-section.

    The EBSD analysis performed on top surface and in the cross-section uncovered the nature of grain structure and the orientation effects in these planes (XY and XZ planes respectively)(Fig.5).The top surface grain structure had an average grain size of 2.5 μm as indicated by the inverse pole figur(IPF) maps from the XY plane is nearly 5 times smaller than grain size (~10 μm) of untreated AZ31B Mg material (Fig.5(a)).Corresponding 0001 pole figur texture plot and image quality (IQ) map are presented in Fig.5 (b)-(c) respectively.On the same lines,the IPF maps from the cross-sectional XZ plane with corresponding 0001 pole figur texture plot and image quality (IQ) map are presented in Fig.5 (d) and (f) respectively.Although,the average grain size of LSM AZ31B Mg in this section (XZ and XY plane) is smaller than that in untreated region (Fig.5 (a) and (d)),the grains are relatively smaller in the top region compared to those in the bottom region near the interface with the underlying base microstructure.This is likely as the thermokinetic parameters vary within the melt pool,the solidificatio rate near the bottom of the melt pool is lowest (1-10 mm/s) and increases rapidly towards the top and attains its maximum solidificatio rate(500 mm/s)near the surface.The laser processing usually results in high thermal gradients (105-107K/m) mainly directed towards laser heat source (along Z direction) [40].This produces strong crystallographic texture in LSM region.In the bottom region of the LSM (Fig.5 (d)-(e)),with competitive epitaxial growth accompanied by the highest thermal gradient along Z direction,the grains grow in the<0001>direction,which is the highest thermal conductivity direction for HCP structure of AZ31B Mg.Consequently,prismatic plane direction appears to lie in the cross sectional XZ plane.This is further confirme by 0001 pole figure (Fig.5 (d)-(e)) that clearly show a fibe texture with peak intensity away from the basal in the cross-section (XZ plane) and corresponding higher intensity of basal poles in the top section (XY plane).Along the same lines,majority of grains on XY plane were aligned close to basal plane (0001) (Fig.5 (a)-(b)).Further evaluations indicated that 40% of the grains in the XY planes had (0001) texture.Therefore,the results demonstrated that the laser surface modificatio did not only evolve a distinct cellular grain morphology but it also influence crystallographic texture in AZ31B Mg alloy.Such a distinct crystallographic texture with unique microstructural morphology is likely to have an effect on various properties such as corrosion resistance and biomineralization,which is discussed in the subsequent sections.Apart from crystallographic texture,the surface topography/texture also influence the response of an implant material to a physiological environment.In light of this,as one of the preceding step to thein-vitroevaluation,surface roughness and topography of the samples were examined.

    3.2.Topographic characterization

    The three-dimensional surface morphologies of untreated AZ31B coupon and three AZ31B Mg samples processed with various laser process conditions are presented in Fig.6.The smoothest surface corresponded to the lowest laser fluenc sample(Fig.6(a)).The roughness information of samples was expressed in area root mean square height (Rz)and area arithmetic mean height

    Fig.6.Three-dimensional surface morphology of (a) untreated AZ31B Mg,(b) 1.06 J/mm2, (c) 2.12 J/mm2, and (d) 3.18 J/mm2 LSM AZ31B Mg sample.

    (Ra).The overall surface roughness of the samples increased with an increase in the laser fluenc (Fig.6).The higher laser fluenc incident on the surface results in deeper melting depths.In addition,the temperature profil via computational modeling further indicated the higher vaporization occurrence with higher laser power treatment (Fig.3).The peak temperatures of 2057 K and 2728 K for 3.18 J/mm2and 2.12 J/mm2laser fluence respectively were higher than the evaporation temperature (Tv) of Mg (1373 K) in Fig.3.Therefore,a portion of material will be evaporated from the surface thereby affecting both the roughness and chemical composition within the surface/subsurface region via physical phenomena such as flui fl w and recoil pressure effects[33,50].The maximum temperature of 1285 K for the lowest laser fluenc processing (1.06 J/mm2) was lower than the Tv,hence generating lowest roughness profil among the three different laser processing conditions employed in the present work.All the other conditions experienced temperature higher than Tv.Evolution of such a rough surface topography is expected to add to the surface energy.High energy surfaces are deemed suitable for the bioactivity of the implant as they can have better affinit to the plasma.In order to evaluate the combined effect of crystallographic texture,phase evolution,and surface texture;the samples were subjected toin-vitroevaluation as explained in the following subsections.

    3.3.Wettability in SBF

    The contact angle in SBF steadily decreased with an increase in the laser fluenc (Fig.7).Moreover,there was an inverse correlation between surface roughness (Ra) and SBF contact angle (Fig.7).These observations indicated that surface became hydrophilic for higher laser fluence Such as modifie surface would be suitable for bio-integration and cell attachment [27,51,52].It is worth mentioning here that the surface wetting behavior is dictated by surface energy of the sample under consideration.Surface energy also bares an inverse relationship with contact angle wherein a hydrophilic(low contact angle) sample possesses a relatively higher surface energy [27].Roughness,chemical composition,and microstructure (grain size,crystallographic texture,and phase make up) are some of the key factors which influenc the surface energy and thus the wetting behavior of a material[25,27,29].Based on the discussion presented in previous sections,for all the LSM AZ31B Mg samples the grain size(~1.5 μm) and crystallographic orientation (predominantly basal plane on the top surface) were nearly same due to partially narrow range of thermokinetic conditions experienced by these samples under processing conditions employed in the present work.However,there were two key differentiating factors among the LSM samples:i) surface roughness and ii) fraction ofβphase both of which increased as a function of laser fluenc likely due to composition variation (Al enrichment) experienced at higher laser fluences

    Fig.7.SBF contact angle and roughness parameter as functions of laser processing condition for both untreated and LSM AZ31B Mg samples.

    Fig.8.Surface morphology of untreated and LSM AZ31B Mg samples after various immersion periods in SBF.

    For a biomaterial,SBF contact angle (wettability) is related to cell attachment and growth.Reports indicated that the cell adhesion and growth is not favorable for both highly hydrophobic and hydrophilic surface [52,53].Arima et al.reported the optimal contact angle for cell adhesion is in the range 40?-60?[51].Webb et al.pointed out the hydrophilic surfaces with contact angle 20?-40?promoted the highest levels of cell attachment [54].It is clear that cell growth and adhesion are favored with higher surface wettability.To examine the influenc of wetting behavior,samples were subjected to immersion studies in SBF and were observed at regular intervals as presented in following subsections.

    3.4.Biomineralization

    Examination of nature of mineral precipitated on a material upon immersion in SBF is a good indicator of biointegration tendency of an implant material.As a firs level of observations,surface morphology of samples was examined optically and differences were noted as a function of immersion time in SBF (Fig.8).Untreated AZ31B Mg after 1-day immersion displayed several localized pits due to corrosion.Furthermore,the precipitated products did not have a uniform distribution onto the surface.In addition,the severity of pitting increased with time leading to separation of the precipitates from the surface dissolving back into the surrounding SBF solution(Fig.8).

    On the contrary,LSM AZ31B Mg samples experienced much lesser pitting compared to the untreated AZ31B Mg sample during all the time periods of immersion (Fig.8).In addition,the precipitates appeared to remain adhered to the surface with greater stability.While comparing all laser processed samples together the samples processed with higher laser fluence appeared to mineralize better (Fig.8).It has been reported that a refine microstructure apparently stabilizes the layers of deposits [55].Furthermore,the increased fraction of grain boundaries would also contribute to the higher surface energy as indicated in the previous publication by current research group [27].

    The nature of mineral deposited is also a critical parameter which would facilitate cell attachment and growth on an implant surface.The mineralized surfaces were observed using SEM+EDS and XPS techniques.Observations of SEM+EDS suggested predominant presence of Mg and O on the as mineralized surface of untreated AZ31B Mg sample (Fig.9 (a)).In contrast,the peaks of Ca and P are strong for LSM AZ31B Mg samples indicating a shift in composition of the mineralized products (Fig.9(b-d)).As observed in previous studies,the Mg rich deposit is associated with a plate like morphology,and the Ca and P rich deposit with a globular morphology [27,29].

    It is worth mentioning here that SEM+EDS is a semi quantitative technique and cannot accurately determine the exact nature of the precipitated mineral.This was also a gap in the previously published work related to laser treated Mg based implant materials by the current group [27,29] where in SEM+EDS was used to semi quantitatively evaluate the mineralized surfaces.To fil in this gap and expand the understanding about mineralized products further,XPS technique was utilized to obtain the exact nature of the mineral products (Fig.10).The Ca element split into two peaks of 2P1/2and 2P3/2attributed to the spin orbital splitting,and attributed to compound Ca10(PO4)6(OH)2(calcium hydroxyapatite,HA) [56] (Fig.10(a)).Even though,this peak splitting of Ca was observed for all the samples,the two peaks for higher laser fluence were fully resolved (Fig.10(a)).These peaks matched well with Ca2+from HA mineral as shown by fitte results in Fig.10(b).However,peaks shifted toward lower binding energy as observed for both untreated and 1.06 J/mm2laser fluenc samples,suggesting precipitation of multiple types of calcium apatite (eg.octacalcium phosphate and dicalcium phosphate) [57,58].

    Fig.9.SEM surface morphologies of the (a) untreated AZ31B Mg,(b) 1.06 J/mm2, (c) 2.12 J/mm2, and (d) 3.18 J/mm2 after 7 days immersion in SBF along with corresponding EDS spectra.

    Fig.10.XPS analysis of apatite precipitated on surface of both untreated AZ31B Mg and all LSM samples after 7-days immersion in SBF:(a) Ca 2P (b)Curve fittin of Ca 2P showing matching with calcium phosphate Ca bonds,(c) O 1s peaks (d) curve fittin and deconvolution of O 1S showing calcium phosphate O-2 bonds for high laser fluenc samples and OH- bond for untreated AZ31B Mg sample indicating hydroxide formation,(e) P 2P peaks (f) curve fittin P 2P peaks showing P bonds matching with calcium phosphate for LSM AZ31B Mg samples,and (g) Ca/P atomic ratios for various samples.

    Oxygen appeared as twin peak for untreated AZ31B Mg samples revealing the presence of a shoulder peak(Fig.10(c)).On the other hand,such a shoulder peak was absent for LSM samples AZ31B samples (Fig.10(c)).The broader peak for untreated sample was deconvoluted to resolve the nature of oxygen bonding (Fig.10(d)).Observation of the resolved peak for untreated AZ31B Mg sample indicated presence of both OH-bonding state of oxygen and the presence of O2-found in HA mineral (Fig.10(d)) [59].The precipitates in LSM AZ31B Mg samples predominantly possessed O2-bonding state (as an example for 3.18 J/mm2in Fig.10(d)).

    The bonding state of element P in precipitates of the untreated AZ31B Mg sample is associated with the other apatite other than HA mineral (Fig.10(c)).On the contrary,with presence of a small shoulder peak with the bonding phases of P,both HA-like and non-HA like phases appeared to be present in the precipitates of the AZ31B Mg samples treated with low laser fluenc of 1.06 J/mm2(Fig.10(b)).AZ31B Mg sample processed with 3.18 J/mm2laser fluenc upon deconvolution of the P peak revealed the combinations predominantly consisting ofalong with presence(Fig.10(d)) [59].This is an indication that the mineralized surface of laser melted samples primarily consisted of HA along with minor amounts of other types apatite.

    As a concluding examination of XPS results,Ca/P peak intensity ratio was calculated from the spectra.The AZ31B Mg samples treated at 2.12 J/mm2and 3.18 J/mm2were associated with a much closer value to stoichiometric HA mineral(Ca/P:1.67) confirmin the individual elemental peak analyses.On these lines,the untreated AZ31B Mg sample had a much lower value of Ca/P ratio illustrating formation of nonstoichiometric HA like minerals and hydroxides upon immersion in SBF solution.It is worth noting here that apart from HA there are 7 different types of apatite minerals with Ca/P ratio in the range of 1-2 [57,58].

    Apart from observations of the mineralized surface and understanding the chemical nature of the deposited products,measuring the weight loss as a function of time is an important indicator of the implant performance in the physiological environment.The mineralized samples were cleaned using chromic acid and their weight loss was measured to calculate the corrosion rate.In case of untreated AZ31B Mg sample,a continuous weight loss was due to the hydration of surface via formation of a hydroxide/oxide layer that in turn simultaneously continued to dissolve in surrounding SBF (Fig.11).It has been reported that the Mg hydroxide is especially unstable in chloride ion solution such as SBF [10].On the other hand,the corrosion resistance of laser processed AZ31B Mg samples was enhanced with immersion period as indicated by a negative slope (Fig.11).This behavior can be attributed to the formation of apatite layer on the surface (as identifie by XPS analysis in Fig.10).Apatite layer is thought to protect the underlying Mg surface from corrosion due to its extremely low solubility in SBF (-log(Ksp)25?C)=116.8) [60].

    Fig.11.Corrosion rate as a function of immersion time for untreated AZ31B Mg and LSM AZ31B Mg samples in SBF as measured by weight loss technique.

    The results indicate that laser surface modificatio facilitate precipitation of apatite,possibly due to higher energy state of the fin grain boundary region acting as active nucleation site.Combining with the improvement in wettability (lower SBF contact angle) represents the accelerated mineralization rate for higher laser fluenc treated samples.Additionally,several reports [61,62] indicated the closed pack crystallographic planes showed higher biocompatibility for HCP biomaterials.The results suggest that the substantial increase in the volume fraction (40%) of the crystallographic orientation of (0001)basal plane of HCP structure of Mg in laser treated surface in this work was highly beneficia to biocompatibility (Fig.5).Furthermore,the basal texture has also been reported to be beneficia for corrosion resistance [63].According to the equation proposed by Song et al.[63],the electrochemical dissolution rate of basal plane (0001) was 18~20 times slower thanplanes of AZ31B Mg alloy.

    where n is the number of electrons involved in the electrochemical reaction;k is a reaction constant;F,R,T and E are Faraday constant,gas constant,absolute temperature and the electrode potential,respectively.αis a transit coefficient Q is an activation energy for a metallic ion to escape from the metal lattice and dissolve into the solution.

    Due to the high value of Q for a densely packed atomic plane,the dissolution rate of basal plane was slower compared to any other loosely packed atomic plane.In the present work,the~40% volume fraction of (0001) plane in 3.18 J/mm2LSM AZ31B Mg sample led to reduced corrosion rate of LSM samples as a result of enhanced growth of mineral apatite.

    4.Conclusions

    The LSM process with various laser fluenc (1.06,2.12,and 3.18 J/mm2) led to a reduction in grain size from 10 μm for untreated sample to 2.5 μm for LSM AZ31B Mg samples.The grain refinemen was accompanied by formation of a cellular microstructure withβphase decorating the cell boundaries.EBSD analysis uncovered presence of 40%area fraction of the basal(0001)crystallographic texture on the surface and sub surface regions.Correspondingly,in the cross-section,the grain structure was associated predominantly prismatic crystallographic texture.From topography point of view,the surface roughness (Ra) of LSM AZ31B Mg samples steadily increased as a function of laser fluenc from 1.34 μm for 1.06 J/mm2to 3.95 μm for 3.18 J/mm2laser fluence The microstructure and surface texture significantl influence thein-vitroevaluation of the LSM AZ31B Mg samples.The contact angle in SBF steadily reduced from 65?for untreated sample to 43?for 3.18 J/mm2LSM AZ31B Mg sample.Furthermore,the Ca/P ratio in the mineral phase(1.683-1.691)in LSM AZ31B Mg samples was close to the stoichiometric HA mineral (1.67),whereas,the untreated sample was associated with a Ca/P ratio of 1.526.Corrosion resistance measured via weight loss revealed a steady increase in corrosion rate from 4 to 8 mm/year as a function of time for the untreated AZ31B Mg sample.For LSM AZ31B Mg samples,the corrosion rate steadily decreased from 2.5 to less that 0.5 mm/year as a function of time.Therefore,the combined effects of crystallographic texture,grain refinement phase evolution,and surface texture induced by LSM process led to a significan improvement inin-vitroresponse of AZ31B Mg alloy demonstrated by a low contact angle in SBF,formation of HA like mineral upon 7 day immersion in SBF,and a concurrent reduction in corrosion rate.Present work helped to shed light on effects if LSM process and further reinforced its suitability for improving performance of Mg based consumable implant materials.

    Declaration of Competing Interest

    None.

    Acknowledgment

    The authors acknowledge the Materials Research Facility(MRF) at the University of North Texas for access to microscopy and phase analysis facilities.Authors thank Dr.Sundeep Mukherjee for the access to the optical profilomete.

    欧美国产日韩亚洲一区| 99精品欧美一区二区三区四区| 99国产综合亚洲精品| 麻豆国产av国片精品| www.熟女人妻精品国产| 午夜激情av网站| 波多野结衣高清无吗| 亚洲精品美女久久久久99蜜臀| 亚洲一码二码三码区别大吗| 亚洲成av人片免费观看| 在线免费观看的www视频| 在线观看66精品国产| 一二三四在线观看免费中文在| 亚洲欧美精品综合一区二区三区| 久久国产亚洲av麻豆专区| 亚洲欧美日韩高清在线视频| 高潮久久久久久久久久久不卡| 级片在线观看| 欧洲精品卡2卡3卡4卡5卡区| 久久久久精品国产欧美久久久| 久久久久国产一级毛片高清牌| 欧美一级毛片孕妇| 欧美色欧美亚洲另类二区 | 久久久久久久精品吃奶| 午夜福利在线观看吧| www.www免费av| 久久精品成人免费网站| 免费看十八禁软件| 美国免费a级毛片| 女性被躁到高潮视频| 成年版毛片免费区| 黄色 视频免费看| 天天躁狠狠躁夜夜躁狠狠躁| 伊人久久大香线蕉亚洲五| 婷婷六月久久综合丁香| 久久香蕉国产精品| 亚洲七黄色美女视频| 免费人成视频x8x8入口观看| 又大又爽又粗| 性少妇av在线| 在线播放国产精品三级| 中文字幕av电影在线播放| 精品国产国语对白av| 色老头精品视频在线观看| 老熟妇仑乱视频hdxx| 日韩 欧美 亚洲 中文字幕| 黄色视频,在线免费观看| 无人区码免费观看不卡| av电影中文网址| av中文乱码字幕在线| 亚洲性夜色夜夜综合| 非洲黑人性xxxx精品又粗又长| 日日摸夜夜添夜夜添小说| 每晚都被弄得嗷嗷叫到高潮| 天天添夜夜摸| 精品日产1卡2卡| 不卡av一区二区三区| 国产xxxxx性猛交| 最近最新中文字幕大全免费视频| www.熟女人妻精品国产| 成人av一区二区三区在线看| 久久精品影院6| a在线观看视频网站| xxx96com| ponron亚洲| av有码第一页| 老司机午夜福利在线观看视频| 超碰成人久久| 两个人免费观看高清视频| 夜夜爽天天搞| 国产在线精品亚洲第一网站| 免费在线观看影片大全网站| 天天添夜夜摸| svipshipincom国产片| 黄色成人免费大全| 国产91精品成人一区二区三区| 这个男人来自地球电影免费观看| 国产99久久九九免费精品| 午夜视频精品福利| 亚洲av电影在线进入| aaaaa片日本免费| 久久精品国产综合久久久| 亚洲精品一区av在线观看| 日韩欧美免费精品| 伦理电影免费视频| 亚洲熟女毛片儿| 亚洲av美国av| 日本vs欧美在线观看视频| 久久国产亚洲av麻豆专区| 午夜两性在线视频| 中文亚洲av片在线观看爽| 91成人精品电影| 亚洲av日韩精品久久久久久密| 夜夜爽天天搞| a级毛片在线看网站| 热99re8久久精品国产| 操美女的视频在线观看| 淫妇啪啪啪对白视频| 日韩精品中文字幕看吧| 国产午夜福利久久久久久| 97人妻精品一区二区三区麻豆 | 无限看片的www在线观看| 亚洲一区中文字幕在线| 日本黄色视频三级网站网址| 91成人精品电影| 深夜精品福利| 国产成人精品在线电影| 母亲3免费完整高清在线观看| 欧美人与性动交α欧美精品济南到| 午夜老司机福利片| 国产av一区在线观看免费| 国产激情欧美一区二区| 9热在线视频观看99| 日本免费一区二区三区高清不卡 | 精品久久蜜臀av无| 色综合婷婷激情| 午夜免费成人在线视频| 黄色视频不卡| 久9热在线精品视频| 久久久久久久久中文| 91字幕亚洲| 免费高清在线观看日韩| 国产精品永久免费网站| 男人舔女人的私密视频| 免费搜索国产男女视频| 欧美一级a爱片免费观看看 | 伊人久久大香线蕉亚洲五| 俄罗斯特黄特色一大片| 色尼玛亚洲综合影院| 妹子高潮喷水视频| 国产精品爽爽va在线观看网站 | 国产精品 欧美亚洲| 亚洲成人免费电影在线观看| 亚洲熟妇中文字幕五十中出| 免费看十八禁软件| 脱女人内裤的视频| 99久久久亚洲精品蜜臀av| 丰满的人妻完整版| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成+人综合+亚洲专区| av欧美777| 亚洲五月婷婷丁香| 黑人操中国人逼视频| 美女免费视频网站| 夜夜爽天天搞| 亚洲国产欧美网| 精品久久久久久,| 久久亚洲真实| 99精品在免费线老司机午夜| av在线播放免费不卡| 91国产中文字幕| 校园春色视频在线观看| 老司机在亚洲福利影院| 欧美日韩中文字幕国产精品一区二区三区 | 夜夜躁狠狠躁天天躁| 欧美一区二区精品小视频在线| 亚洲av电影在线进入| 欧美另类亚洲清纯唯美| 久久婷婷成人综合色麻豆| av在线天堂中文字幕| 精品久久久久久成人av| 999久久久国产精品视频| 国产欧美日韩一区二区三| 欧美日韩中文字幕国产精品一区二区三区 | 国产在线精品亚洲第一网站| 亚洲自拍偷在线| 亚洲五月天丁香| 国产精品电影一区二区三区| 婷婷六月久久综合丁香| 亚洲成a人片在线一区二区| 欧美日韩精品网址| 黄色视频不卡| 丰满人妻熟妇乱又伦精品不卡| 亚洲中文字幕一区二区三区有码在线看 | 亚洲情色 制服丝袜| 男人舔女人的私密视频| 十八禁人妻一区二区| 亚洲av片天天在线观看| 美女高潮喷水抽搐中文字幕| 中出人妻视频一区二区| 精品国产亚洲在线| 69av精品久久久久久| 日韩大尺度精品在线看网址 | 国产aⅴ精品一区二区三区波| ponron亚洲| 亚洲av熟女| 精品久久久久久久久久免费视频| 丁香欧美五月| 国产精品美女特级片免费视频播放器 | 国产一区二区激情短视频| 成年女人毛片免费观看观看9| 亚洲熟女毛片儿| 成在线人永久免费视频| 黄片播放在线免费| 欧美黑人精品巨大| 18禁黄网站禁片午夜丰满| 国产激情欧美一区二区| 久久久久亚洲av毛片大全| 最近最新免费中文字幕在线| 精品欧美一区二区三区在线| 午夜激情av网站| 琪琪午夜伦伦电影理论片6080| 精品午夜福利视频在线观看一区| 亚洲中文字幕一区二区三区有码在线看 | 99精品在免费线老司机午夜| 99re在线观看精品视频| 在线观看免费午夜福利视频| 国产视频一区二区在线看| 淫秽高清视频在线观看| 国产在线精品亚洲第一网站| 欧美av亚洲av综合av国产av| 大香蕉久久成人网| 好看av亚洲va欧美ⅴa在| 91麻豆精品激情在线观看国产| 美女 人体艺术 gogo| 欧美日韩精品网址| 91成年电影在线观看| 伦理电影免费视频| 国产一区在线观看成人免费| 69精品国产乱码久久久| 亚洲国产精品合色在线| 极品教师在线免费播放| 免费观看精品视频网站| 亚洲成人国产一区在线观看| 欧美另类亚洲清纯唯美| 亚洲成人精品中文字幕电影| 国产成人一区二区三区免费视频网站| 亚洲国产高清在线一区二区三 | 久久精品成人免费网站| 久久久久国产精品人妻aⅴ院| 啦啦啦 在线观看视频| 亚洲avbb在线观看| 不卡一级毛片| 久久精品人人爽人人爽视色| 免费女性裸体啪啪无遮挡网站| 久久久国产成人精品二区| 成人国产综合亚洲| 天天躁夜夜躁狠狠躁躁| 亚洲熟妇中文字幕五十中出| 亚洲五月婷婷丁香| 国产99久久九九免费精品| 在线观看66精品国产| 亚洲国产看品久久| 亚洲精品中文字幕一二三四区| 久久九九热精品免费| 亚洲国产精品久久男人天堂| 亚洲精品一卡2卡三卡4卡5卡| 精品久久久久久,| а√天堂www在线а√下载| 中文字幕色久视频| 久久精品亚洲精品国产色婷小说| 午夜福利高清视频| 欧美日本中文国产一区发布| 欧美色视频一区免费| 精品一区二区三区av网在线观看| 国产一区二区三区视频了| www.999成人在线观看| tocl精华| 亚洲 欧美一区二区三区| 国产精品久久久久久亚洲av鲁大| 成在线人永久免费视频| 亚洲欧美日韩高清在线视频| 日韩欧美一区视频在线观看| 可以免费在线观看a视频的电影网站| 亚洲国产欧美日韩在线播放| 午夜免费激情av| 超碰成人久久| 中文字幕另类日韩欧美亚洲嫩草| 在线观看舔阴道视频| 久久精品影院6| 国产av一区在线观看免费| 美女大奶头视频| 日本一区二区免费在线视频| 国产成人欧美| 制服诱惑二区| 夜夜看夜夜爽夜夜摸| 最近最新中文字幕大全免费视频| 欧美乱色亚洲激情| 欧美午夜高清在线| 日韩欧美在线二视频| 欧美一区二区精品小视频在线| 亚洲美女黄片视频| 亚洲专区字幕在线| 国产午夜福利久久久久久| 少妇被粗大的猛进出69影院| 国产伦人伦偷精品视频| 欧美黄色淫秽网站| 国产1区2区3区精品| 999久久久国产精品视频| 国产午夜精品久久久久久| 午夜精品久久久久久毛片777| 丁香六月欧美| 久久久国产成人精品二区| 美女国产高潮福利片在线看| 男人舔女人下体高潮全视频| 色播亚洲综合网| av天堂久久9| 亚洲精品中文字幕在线视频| 久久久国产精品麻豆| 久久久久久亚洲精品国产蜜桃av| 色老头精品视频在线观看| 免费av毛片视频| tocl精华| 欧美乱色亚洲激情| 丝袜美腿诱惑在线| netflix在线观看网站| 亚洲精品粉嫩美女一区| 国产激情欧美一区二区| 亚洲熟妇熟女久久| 99精品在免费线老司机午夜| 变态另类丝袜制服| 国产精品电影一区二区三区| 久久久水蜜桃国产精品网| 伦理电影免费视频| 在线观看一区二区三区| 老司机福利观看| 亚洲少妇的诱惑av| 亚洲成a人片在线一区二区| 好男人电影高清在线观看| 后天国语完整版免费观看| 亚洲五月色婷婷综合| 两个人视频免费观看高清| 9191精品国产免费久久| 久久久久久大精品| 色综合婷婷激情| 欧美激情极品国产一区二区三区| 亚洲精品av麻豆狂野| 男女床上黄色一级片免费看| 怎么达到女性高潮| 天堂√8在线中文| 制服人妻中文乱码| 美女高潮喷水抽搐中文字幕| 两个人看的免费小视频| 黄片小视频在线播放| 多毛熟女@视频| 午夜福利18| 国产精品秋霞免费鲁丝片| 国产成人欧美| 日韩有码中文字幕| 日韩欧美国产一区二区入口| 国内精品久久久久久久电影| 亚洲精品在线观看二区| 免费不卡黄色视频| 男女午夜视频在线观看| 少妇的丰满在线观看| 免费看a级黄色片| 非洲黑人性xxxx精品又粗又长| 天天躁夜夜躁狠狠躁躁| 极品人妻少妇av视频| ponron亚洲| av福利片在线| 神马国产精品三级电影在线观看 | 国产精品爽爽va在线观看网站 | 9热在线视频观看99| 欧美成人性av电影在线观看| 9热在线视频观看99| 大型av网站在线播放| www日本在线高清视频| 精品免费久久久久久久清纯| 亚洲国产精品合色在线| 亚洲精华国产精华精| 午夜激情av网站| 久久亚洲真实| 欧美成人午夜精品| 亚洲色图 男人天堂 中文字幕| 夜夜看夜夜爽夜夜摸| 一区二区三区精品91| 一区二区三区激情视频| 欧美成狂野欧美在线观看| 亚洲欧美激情在线| 一本久久中文字幕| 免费在线观看影片大全网站| 无人区码免费观看不卡| 亚洲精品一区av在线观看| 免费不卡黄色视频| 久久国产精品男人的天堂亚洲| 中文字幕色久视频| 日本 av在线| 法律面前人人平等表现在哪些方面| 精品一品国产午夜福利视频| 国产精品一区二区在线不卡| 久久天堂一区二区三区四区| 午夜免费鲁丝| 午夜福利免费观看在线| 搡老熟女国产l中国老女人| 亚洲性夜色夜夜综合| 91九色精品人成在线观看| 亚洲九九香蕉| 亚洲一区中文字幕在线| 这个男人来自地球电影免费观看| 身体一侧抽搐| bbb黄色大片| 日韩欧美国产一区二区入口| 亚洲一区高清亚洲精品| 一a级毛片在线观看| 亚洲男人的天堂狠狠| 国产主播在线观看一区二区| 黄色丝袜av网址大全| 国产蜜桃级精品一区二区三区| 国产精品一区二区在线不卡| 亚洲熟妇熟女久久| 99久久99久久久精品蜜桃| 久久中文字幕一级| 亚洲av五月六月丁香网| 这个男人来自地球电影免费观看| 999久久久精品免费观看国产| tocl精华| 久久人人97超碰香蕉20202| 国产麻豆成人av免费视频| 日韩欧美国产一区二区入口| 亚洲免费av在线视频| 搡老岳熟女国产| 午夜久久久在线观看| 色综合亚洲欧美另类图片| 制服人妻中文乱码| 国产一区二区在线av高清观看| 久久久久久久久中文| 国内毛片毛片毛片毛片毛片| 成年版毛片免费区| 成年人黄色毛片网站| 国产成人欧美| 桃色一区二区三区在线观看| 精品一品国产午夜福利视频| 两个人看的免费小视频| videosex国产| 99国产精品一区二区三区| www.999成人在线观看| 亚洲欧美激情综合另类| 久久久精品国产亚洲av高清涩受| 久久久久国产精品人妻aⅴ院| 久久久国产精品麻豆| 最新美女视频免费是黄的| 亚洲片人在线观看| 国产97色在线日韩免费| 成人国产一区最新在线观看| a在线观看视频网站| 91在线观看av| 中文字幕人妻丝袜一区二区| 桃色一区二区三区在线观看| 咕卡用的链子| 亚洲成人国产一区在线观看| 十分钟在线观看高清视频www| 妹子高潮喷水视频| 久久精品成人免费网站| 国产精品国产高清国产av| 黄色视频,在线免费观看| 国产精品电影一区二区三区| 人人妻人人爽人人添夜夜欢视频| 亚洲国产精品久久男人天堂| 啦啦啦免费观看视频1| 桃色一区二区三区在线观看| 日韩免费av在线播放| 亚洲五月婷婷丁香| 亚洲精品一卡2卡三卡4卡5卡| 18禁美女被吸乳视频| 91九色精品人成在线观看| 欧美乱码精品一区二区三区| 精品高清国产在线一区| 欧美午夜高清在线| 色尼玛亚洲综合影院| 男人舔女人下体高潮全视频| 欧美丝袜亚洲另类 | 国产亚洲精品综合一区在线观看 | 欧洲精品卡2卡3卡4卡5卡区| 国产又爽黄色视频| 久久久久久亚洲精品国产蜜桃av| www.自偷自拍.com| 在线观看免费午夜福利视频| 身体一侧抽搐| 午夜精品在线福利| 亚洲五月色婷婷综合| 在线观看午夜福利视频| 亚洲 国产 在线| 精品欧美一区二区三区在线| 日韩欧美三级三区| 精品久久久精品久久久| 1024视频免费在线观看| 欧美+亚洲+日韩+国产| 欧美成人免费av一区二区三区| 免费av毛片视频| 日日摸夜夜添夜夜添小说| 久久人妻熟女aⅴ| 精品不卡国产一区二区三区| 久99久视频精品免费| 97人妻精品一区二区三区麻豆 | 色播在线永久视频| 老鸭窝网址在线观看| 妹子高潮喷水视频| 欧美黄色淫秽网站| 国产精品免费一区二区三区在线| 9191精品国产免费久久| 中国美女看黄片| 天天躁夜夜躁狠狠躁躁| 免费久久久久久久精品成人欧美视频| 天堂√8在线中文| 无限看片的www在线观看| 丝袜人妻中文字幕| 久久欧美精品欧美久久欧美| 日韩欧美在线二视频| 黄色视频不卡| 国产精品国产高清国产av| 国产av又大| 亚洲自拍偷在线| 欧美大码av| 91成年电影在线观看| 成年版毛片免费区| 国产极品粉嫩免费观看在线| 在线观看免费视频网站a站| 脱女人内裤的视频| 国产成人精品久久二区二区91| 国产一区二区三区视频了| 国产精品亚洲av一区麻豆| 国产精品,欧美在线| 精品卡一卡二卡四卡免费| a级毛片在线看网站| 999久久久国产精品视频| 亚洲av电影不卡..在线观看| 亚洲男人天堂网一区| 美国免费a级毛片| 十八禁人妻一区二区| 人人妻,人人澡人人爽秒播| 成人18禁在线播放| 91成年电影在线观看| 一级a爱视频在线免费观看| 国产极品粉嫩免费观看在线| 欧美黑人精品巨大| 亚洲中文日韩欧美视频| 18美女黄网站色大片免费观看| 午夜福利高清视频| 午夜成年电影在线免费观看| 色综合站精品国产| 国产精品98久久久久久宅男小说| 午夜老司机福利片| 欧美午夜高清在线| 香蕉久久夜色| 国产精华一区二区三区| 成人三级黄色视频| 精品国内亚洲2022精品成人| 亚洲国产精品成人综合色| 天天躁夜夜躁狠狠躁躁| 丰满人妻熟妇乱又伦精品不卡| 国产成人精品在线电影| 国产视频一区二区在线看| 桃色一区二区三区在线观看| 亚洲av美国av| 在线十欧美十亚洲十日本专区| 一区在线观看完整版| av在线播放免费不卡| 国产1区2区3区精品| 亚洲五月天丁香| 久久久久亚洲av毛片大全| 天堂影院成人在线观看| 国产精品av久久久久免费| 动漫黄色视频在线观看| 亚洲第一欧美日韩一区二区三区| 99在线人妻在线中文字幕| 老熟妇乱子伦视频在线观看| 18禁黄网站禁片午夜丰满| 午夜福利一区二区在线看| 男男h啪啪无遮挡| 又大又爽又粗| 亚洲国产精品sss在线观看| 亚洲在线自拍视频| 男女下面插进去视频免费观看| 国产欧美日韩一区二区三区在线| 欧美国产日韩亚洲一区| 久久中文字幕人妻熟女| 999精品在线视频| 免费久久久久久久精品成人欧美视频| videosex国产| 亚洲成人免费电影在线观看| 成年女人毛片免费观看观看9| 久久久水蜜桃国产精品网| 国语自产精品视频在线第100页| 久久久国产成人免费| 精品人妻在线不人妻| 韩国av一区二区三区四区| 国产精华一区二区三区| 午夜福利欧美成人| 在线国产一区二区在线| 国产亚洲av嫩草精品影院| 精品国产一区二区三区四区第35| 成人精品一区二区免费| 午夜免费激情av| 国产xxxxx性猛交| 男男h啪啪无遮挡| 一本久久中文字幕| 巨乳人妻的诱惑在线观看| 日日摸夜夜添夜夜添小说| 免费看十八禁软件| 亚洲国产毛片av蜜桃av| 丁香六月欧美| 亚洲人成电影观看| 国产一区在线观看成人免费| 久久久久久久久中文| 国产色视频综合| av欧美777| 咕卡用的链子| 久久香蕉激情| 欧美午夜高清在线| 欧美国产精品va在线观看不卡| 一级作爱视频免费观看| 国产精品久久久久久精品电影 | 午夜福利欧美成人| 久久久久久久久中文| 给我免费播放毛片高清在线观看| 在线观看一区二区三区| 咕卡用的链子| 欧美日韩一级在线毛片| 亚洲成人国产一区在线观看| 欧美激情久久久久久爽电影 | 色综合亚洲欧美另类图片| 免费在线观看视频国产中文字幕亚洲| 99香蕉大伊视频| 亚洲视频免费观看视频|