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

    Numerical research on effect of overlap ratio on thermal-stress behaviors of the high-speed laser cladding coating*

    2021-01-21 02:15:38XiaoxiQiao喬小溪TonglingXia夏同領(lǐng)andPingChen陳平
    Chinese Physics B 2021年1期
    關(guān)鍵詞:陳平小溪

    Xiaoxi Qiao(喬小溪), Tongling Xia(夏同領(lǐng)), and Ping Chen(陳平)

    School of Mechanical Engineering,University of Science and Technology Beijing,Beijing 100083,China

    Keywords: high-speed laser cladding,overlap ratio,thermal-stress evolution,residual stress and deformation,numerical simulation

    1. Introduction

    Laser cladding is a kind of surface modification technology,which is widely used in the fields of aerospace,automobile, petroleum, and so on. It utilizes a high-energy-density laser beam to melt the alloy power and form a cladding layer metallurgically bonded to the metallic substrate. The cladding coating has the advantages of low dilution ratio and good bond condition, which can significantly improve the anti-corrosion and wear resistance of the substrate surface.[1,2]However,due to the rapid heating and cooling phenomenon of the molten pool during laser cladding, and the difference in thermal expansion properties between the coating and the substrate, the cladding coating has obvious local plastic deformation and residual stress,which will easily lead to cracks and affect the fatigue life of the mechanical parts.[3–6]Therefore, it is crucial to study the thermal-stress evolution, residual stress, and deformation of laser cladding coating.

    The laser cladding process involves complex physical phenomena,such as metallurgy,phase transformation,molten pool expansion, and solidification shrinkage.[7–11]For the multi-track laser cladding coating, its temperature and stress evolution processes are even more complicated. Due to the limitation of experimental methods,[12]the numerical method is widely adopted for its advantages in obtaining the thermalstress evolution and the three-dimensional distribution of residual stress and deformation, which can help to reveal the generation mechanism of residual stress and deformation, as well as the high-risk area caused by the residual stress concentration.

    Residual stress is influenced by the scan speed, power,overlap ratio, etc.[13–15]So far, researchers have done many numerical simulations on the thermal-mechanical evolution during the laser cladding process. The models of powder flow and molten pool flow fields are improved to study the complex physical phenomena of laser cladding.[16–20]Farahmand et al.[15]established a numerical model to study the multitrack laser cladding process with an overlap ratio of 25%.Results show that the residual stress is high at the interface of the cladding coating and substrate, where it has a high crack tendency. Sun et al.[21]developed a numerical model with an overlap ratio of 33% to study Ni60A coatings, and found that the longitudinal residual stress(LRS)increased the crack susceptibility. Chew et al.[22]developed a 3D numerical model with an overlap ratio of 50% to investigate the residual stress for the laser cladding of AISI 4340 steel powder onto a similar substrate material.Results indicated that the deposited cladding layer could reduce the tensile residual stress of the previous track cladding coating,and the resulting local stress peak may damage the fatigue performance. Ghorashi et al.[23]also analyzed the residual stress of the multi-track laser cladding process numerically and captured the stress relaxation and hardening phenomena.

    In practical industrial applications, due to the limitation of the width and thickness of single-track laser cladding coating, multi-track laser cladding coating is used to cover the whole work surface and to improve its corrosion resistance and wear resistance. In the multi-track laser cladding process, the overlap ratio is an important parameter, which will directly affect the quality of multi-track cladding coating.[24]Meanwhile, the corresponding thermal-stress evolution during the laser cladding process,the residual stress distribution,and residual deformation will be more complex;and the complex interactions between each track cladding layer, as well as between the cladding coating and the substrate cannot be neglected. Prasad et al.[25]found experimentally that increasing the overlap ratio can improve the wear resistance of the coating. Zhao et al.[26]studied the effects of overlap ratio on the mechanical property of YCF101 laser cladding coating, and results showed that overlap mode could improve the re-melting area, control the cracks and pores in the cladding coating,and refine the microstructure.

    In the process of laser cladding, scanning speed is also an important parameter, which will affect the coating quality and cladding efficiency. Since 2013,Fraunhofer proposed the high-speed laser cladding technique first,[27,28]and it has gained wide attention. The high-speed laser cladding coating has the advantages of low dilution ratio,small surface roughness, few pores and inclusions, and small thickness, which is more suitable for industrial application.[29]However, during the process of high-speed laser cladding, the heating speed and solidification rate of the molten pool are much higher than those of the traditional cladding coating;thus,the”rapid heating and rapid cooling” phenomenon is more severe than that of the traditional cladding.

    However,research on the influence of the overlap ratio on the thermal-mechanical behaviors of high-speed laser cladding coating is still lacking,and the mechanisms of heat accumulation and thermal deformation need further study.[30]In this paper, the thermal-stress evolution mechanism during the highspeed multi-track laser cladding process,as well as the effect rules of the overlap ratio on the thermal-stress evolution and the residual stress and deformation,are investigated.Section 2 introduces the experimental equipment and the laser cladding parameters. Section 3 introduces the numerical details, including control equation, material parameters, and boundary conditions, as well as the developed geometric model of the high-speed laser cladding coatings with overlap ratios of 10%,30%, and 50% based on the experimental results. Section 4 presents the results and discussion about the influence rules of the overlap ratio, mainly focusing on the thermal-stress evolution mechanism of the high-speed laser cladding,as well as the residual stress and deformation.

    2. Experiment detail

    An RF-J1500 laser cladding equipment is used to prepare the cladding coating,as shown in Fig.1(b). The device adopts a coaxial powder feeding mode,and the schematic diagram of the laser cladding coating process is shown in Fig.1(a).Argon is used as the carrier gas and protective gas. In our study,the cladding parameters used in the numerical simulations are the same as those in experiments. The laser power, scan speed,and feeding rate are 1200 W,25 mm/s,and 120 mg/s,respectively. And the overlap ratios used are 10%, 30%, and 50%.The powder and substrate materials are Ni60A nickel-based self-fluxing alloy and 316L stainless steel,respectively.

    Fig. 1. (a) Schematic diagram of laser cladding. (b) High-speed laser cladding equipment.

    3. Numerical simulation detail

    3.1. Geometric model and meshing

    The cross-section geometry of the high-speed laser cladding coating is different from that of the traditional lowspeed cladding coating. In the simulation, the geometric model of a single-track cladding coating is established based on the experimental result, as shown in Fig. 2(a). Since there is no dilution zone before laser cladding, the geometry model does not need to consider the shape of the dilution zone. The size of the substrate used in the study is 30mm×30mm×4mm, as shown in Fig. 2(b). In order to reduce the discontinuity of stress distribution caused by geometric end effects,[31]the length of the cladding coating in this study is 25 mm, which is slightly shorter than the length of the substrate. The laser scanning direction is defined as the Y-direction (i.e., longitudinal direction), and the vertical direction of laser scanning is defined as the X-direction (i.e.,transverse direction).

    Three overlap ratios of 10%, 30%, and 50% are analyzed in this study, and the corresponding geometric model and mesh model are shown in Figs. 2(e)–2(g). And the track number used in the cladding is five. For the overlap ratio of 10%, its cross-section geometry has no obvious difference with that of the single-track coating due to the small overlapping area between each track coating. Then the geometric model is established, as shown in Fig. 2(e). When the overlap ratios are 30% or 50%, the interaction between each track coating is obvious. According to the coating thicknesses of 0.25 mm and 0.35 mm obtained experimentally,the geometric models established are shown in Figs.2(f)and 2(g).[22]

    The mesh model adopted in the study is shown in Fig. 2(c). In order to improve the accuracy of numerical simulation, the meshes of the cladding coating and its adjacent substrate are encrypted (see Figs. 2(c), 2(e), 2(f), and 2(g)). Further, the grid independence is verified, and the results are shown in Fig. 2(d). According to change rules of the maximum temperature, the maximum stress and strain,and the calculating time varying with mesh number, Grid 3 is finally adopted, which can meet the calculation accuracy requirement. The minimum element size of grid 3 is 0.05mm×0.05mm×0.3mm,the relative error of molten pool temperature is less than 0.5%, and the relative error of maximum residual stress is less than 1.7%.

    Fig. 2. Geometric and mesh models of highspeed laser cladding coating with overlap ratios of 10%, 30%, and 50%: (a) single-track coating geometry model,(b)the geometry model,and(c)the mesh model of the substrate and the multi-tracks coating,(d)mesh independence verification,and the geometric model and mesh model of the coating for overlap ratios of(e)10%,(f)30%,and(g)50%.

    3.2. Material parameters

    The substrate material used is 316L stainless steel, and its melting range is 1400?C–1460?C. The powder material is nickel-based self-fluxing alloy Ni60A,which has excellent corrosion resistance and wear resistance,and its melting range is 960–1040?C. To simplify the analysis, the solidification points of 316L and Ni60A are set to be 1400?C and 960?C,respectively,and the freezing point of 316L and Ni60A contact interface is assumed to be 1180?C.[14]That is to say,when the node temperatures of the substrate material,coating material,and the interface of the two materials are reduced to 1400?C,960?C,and 1180?C,respectively,the node can be regarded as solidified. Specific parameters, such as thermal conductivity,specific heat capacity, young’s modulus, yield strength, and density of the substrate and powder materials,can be found in the literatures.[21,32]

    3.3. Laser heat source model

    Gaussian heat source is used in our simulation, and the model equation is[18,33]

    3.4. Heat conduction differential equation

    In laser cladding,the heat conduction obeys Fourier law,and its partial differential equation is[3,18]

    where ρ is the material density, cp(T) and λ(T) are the specific heat capacity and thermal conductivity varying with temperature T.

    3.5. Boundary conditions

    In the simulation,the initial thermal boundary conditions are

    When t =3600 s, the substrate and cladding coating are regarded to be completely cooled down to the room temperature. In the laser cladding process, the energy loss mainly includes thermal conduction, thermal radiation, and thermal convection.Thus the relationship among them is defined as[15]

    where n is the unit normal vector, ?T/?n is the derivative of temperature in n direction, and hcis the thermal convection coefficient, which is 80 W/m2·K[23]under the condition of forced convection between the carrier and protective gas.σ =5.67×10-8W/m2·K4is Stefan–Boltzmann constant. ε(0 <ε <1) is the emission coefficient, due to the relatively high radiation coefficient of the molten pool, ε =0.8 is used in this study.[34]Γ is the area affected by the laser beam. T is the surface temperature,and T0=20?C is the ambient temperature.

    3.6. Thermal stress coupling analysis method

    The indirect coupling method is adopted for the thermal analysis and stress analysis. Firstly, the temperature field of laser cladding is calculated, and then each step results of the transient temperature field are used as the input of the stress field model. The calculation of the mechanical stress is based on the thermo-elastoplastic model. The substrate material used is austenitic stainless steel,and most of the nickel-based corrosion-resistant alloys have austenitic structure, so the effect of phase transformation needs not to be considered.[35]In this study, elastic strain, thermal strain, and plastic strain are considered,and the total strain is[36,37]

    In the simulation, the element activation method[38]is used. It means that the properties of the coating material are close to zero before the formation of the molten pool. When the temperature exceeds the material melting point,its thermal and mechanical properties are activated. In the laser cladding process, the laser beam not only interacted with the coating powder but also with the substrate,[39]as shown in Fig. 3.In other words, when the heat source passes, both the active cladding material and a part of the substrate surface absorb the laser energy at the same time. Figure 3(a)shows the area where the substrate directly interacts with the laser, and the corresponding thickness of the substrate interacting with laser energy is 0.05 mm. Figure 3(b) shows the energy absorption region.

    Fig.3. (a) The interaction area of laser and substrate. (b) The laser action area during the cladding process.

    4. Results and discussion

    4.1. Validation of heat source efficiency and temperature field

    Due to the energy loss during the photo-thermal conversion process of the laser equipment,the heat source efficiency is calculated by comparing the experimental and numerical simulation results. The results show that when the energy utilization rate is 0.4,the simulation result is consistent with the experimental result in terms of appearance and size of the coating,as shown in Fig.4(b).Further,the solidification rates at AJ points are extracted,as shown in Fig.4(a). The results show that with solidification distance increasing from point A to point J,the solidification rate first decreases sharply and then gradually tends to be stable, which is consistent with the results obtained in the reference[14](see Fig.4(a)). In summary,the results indicate that the temperature field model adopted in this study is accurate.

    Fig.4.Heat source check:(a)the heat source and the corresponding solidification ratio at the crosssection of the coating,(b)the appearance and size of the cladding coating obtained by experiment and simulation.

    4.2. Temperature evolution

    The cloud images of the temperature field and the temperature rate for the overlap ratio of 10%at t=1 s are extracted,as shown in Fig.5.During the laser cladding process,there are hightemperature gradients around the molten pool. Especially in the front of the molten pool,the highest temperature gradient is more than 1200?C/mm (960?C/0.8 mm). Meanwhile,the maximum heating rate of the molten pool is also located in the front of the molten pool,reaching 23180.7?C/s. The main reason is that, due to the high scanning speed, both powder material and the substrate material in the front of the molten pool are heated rapidly to the melting point temperature and melt,but the heat cannot be transferred out immediately.

    Because the heat flux density at the center of the laser heat source is the highest,thus the highest temperature point,along with the critical point of temperature rising and cooling of the molten pool, are just in the center of the laser heat source,and the corresponding highest temperature of molten pool is 1879.2?C.After the laser source passing through,the molten pool begins to cool. The higher energy loss and lower energy input make the cooling rate of the latter part of the molten pool large. The maximum cooling rate is-8909.5?C/s,located at the geometric center of the molten pool,as Fig.5(b)shows.

    The results show that, due to the high scanning speed,the maximum heating rate and the maximum cooling rate obtained in this study are much higher than that of the traditional lowspeed laser cladding. Therefore, the ‘rapid heating and rapid cooling’ phenomenon in the high-speed laser cladding process is more pronounced. However,it does not mean it will result in high residual stress and high residual deformation. In contrast, high scanning speed can greatly reduce the heat absorbed per unit length by the material, which is conducive to reducing the dilution ratio and local plastic deformation of the material near the molten pool.

    The temperature evolution process of multi-track cladding coating with a 10% overlap ratio is extracted further, as shown in Fig. 6. It can be seen that there is apparent heat accumulation generated during the laser cladding process. And the heat accumulation can reheat the previous track cladding coating and preheat the substrate of the following track cladding coating. When the heat is accumulated to a certain extent, the preheat effect will be significantly enhanced, which will influence the residual stress and deformation.[8,40,41]As Fig.6(d)shows,the temperature of the fifth track cladding substrate reaches 310?C during the fourth track cladding process. It will reduce the temperature difference between the cladding coating and the substrate,which is conducive to reduce the residual stress at the position.

    Fig.5. The cloud images of(a)temperature and(b)temperature rate for the overlap ratio of 10%at t=1 s.

    Fig.6. The temperature evolution process of the five-track laser cladding coating with an overlap rate of 10%: (a)t=0.5 s;(b)t=1.5 s;(c)t=2.5 s;(d)t=3.5 s;(e)t=4.5 s;(f)t=10 s in the cooling process.

    During the cooling process after t=5 s,the highest temperature position of the coating is located at the center of the fourth track coating due to heat accumulation, which will result in the highest residual stress and increase the crack susceptibility in this position. Increasing the overlap ratio will gradually shift the highest temperature position in the cooling process to the center of the cladding coating. When the overlap ratio is 30%,the position is in the middle of the third and fourth coating. When the overlap ratio is 50%, the highest temperature position is located in the third track cladding coating.

    At the same time,the scanning speed is also an important factor affecting the highest temperature position in the cooling process. The simulation results in the Refs.[15,23]show the heat accumulation always tends to the last track coating for the laser cladding with a scanning speed of 3 mm/s and an overlap ratio of 50%. While in our study, a high scanning speed of 25 mm/s is used. High scanning speed puts the middle of the substrate in poor heat dissipation condition,and then the final heat accumulate appears in the region.

    Further,the temperature evolution processes of the points on the top center of each track cladding coating are monitored.The positions of the five points, P1, P2, P3, P4, and P5, are shown in Fig.7(a). The corresponding temperature evolution curves for the overlap ratios of 10%,30%,and 50%are shown in Figs.7(b),7(c),and 7(d),respectively.

    Fig.7. (a)The positions of the five monitored points and their corresponding temperature evolution curves for the overlap ratio of(b)10%,(c)30%,and(d)50%.

    Figure 7 shows that the laser heat source will reheat the previous track coating,and there is an apparent second temperature peak caused by the reheat effect during the laser cladding process.For example,the second peak temperature of point P1 rises again to 406.7?C during the second track cladding process,as Fig.7(b)shows. Due to the heat accumulation effect,the second peak temperatures of each monitored point increase with the proceeding of the laser cladding. And the laser heat source can preheat the substrate of the next track coating;thus,the maximum temperature of each monitored point increases with the proceeding of the laser cladding. From Fig. 7(b),we can see that the maximum temperature of the point P1 is 1879.2?C for the overlap ratio of 10%, while the maximum temperature of the point P2 rises to 1902.3?C.

    Meanwhile, the maximum temperature and the second peak temperature of each monitored point also increase with the overlap ratio increasing. For example,with the overlap ratio increasing from 10% to 50%, the maximum temperature of the point P5 changes from 1957.3?C to 2028.1?C,as well as the second peak temperature of the point P1 changes from 406.7?C to 967.1?C.When the overlap ratio increases to 50%,the second peak temperature of the point P1 reaches 967.1?C,higher than the melting point of Ni60A(960?C),which means that more than half of each track cladding coating will be remelted. This“re-melting”phenomenon[42]can strengthen the metallurgical bonding between the adjacent track coating and improve the cladding coating performance.[25]

    4.3. Stress and deformation evolution

    Firstly,the stress and strain distribution of the coating and substrate at t =1 s is analyzed, including Von Mises stress(VMS),transverse stress,and longitudinal stress,as shown in Fig. 8. In the laser cladding process, powder materials and a part of the substrate are melted, forming a high-temperature molten pool. Because the VMS and thermal stress of the melted metal are very tiny, then it is convenient to judge roughly the solid-liquid interface of the molten pool with VMS and thermal strain clouds, as Fig. 8(a) shows. The materials adjacent to the molten pool also are heated to high temperatures and will have a large thermal expansion (Fig. 8(b)),but its expansion will be limited by the surrounding lowtemperature materials (Fig. 8(e)), thus resulting in compressive stress and large thermal strain in the area. After the laser heat source passing through, the molten pool material begins to cool down and solidify. Due to the solidification shrinkage of the cladding coating, it generates compressive stress in the transverse direction (see Fig. 8(c), and tensile stress in the longitudinal direction(Fig.8(d)). The complex stress distribution will pull the substrate upward to bend and deform(see Fig.8(e)),which damages the original appearance of mechanical parts and even causes cracks and other defects, thus affecting their corrosion and wear resistance and fatigue life.

    Fig.8. The cloud images of(a)Von Mises stress,(b)longitudinal section thermal strain,(c)longitudinal section transverse stress,and(d)longitudinal section longitudinal stress,as well as(e)the stress and deformation mechanism of the laser cladding process.

    The transverse compressive stress is beneficial to suppress the generation of cracks. In contrast,the excessive longitudinal tensile stress will lead to the generation and propagation of cracks.[40]Especially under the cyclic load working condition, cracks will be easy to occur and propagate at the interface of the coating and substrate. Meanwhile, the stress state of the cladding coating is also closely related to the overlapping ratio, which will make the stress distribution at the joint of coating and substrate more complex. Further,the longitudinal stress evolution process at path 1 for different overlap ratios are given, where path 1 is in the middle of the interface between the coating and substrate and perpendicular to the cladding track,as Fig.9 shows.

    The results show that due to the reheat effect, the deposition of each new track coating reduces the longitudinal stress of the previous track coating,which is defined as‘stress release’.[23]For example, when the overlap ratio is 10%, the stress value decreases from 377 MPa at t=3 s to 199 MPa at t =4 s during the fourth track cladding. Meanwhile, a local stress peak will form at path 1 after each track coating is deposited. And due to the preheat effect, the longitudinal stress peak value decreases with the tracks increasing. That is because the preheat effect can reduce the temperature difference between the cladding coating and the substrate of the next track. For the coating with an overlap ratio of 10%, its local stress peak value changes from 377 MPa at t =3 s to 307 MPa at t=4 s,as Fig.9(b)shows.

    There are stress release and the decreased stress peak phenomena for all the three overlap ratios, which are closely related to the reheat and preheat effects. Moreover,with the increase of the overlap ratio,the lapped area between each track coating increases,and part of the high-temperature surface of the previous track cladding coating will become the substrate of the next track cladding layer. Due to the heat treatment effect between each track coating, the average residual stress for the multi-track coating is relatively lower than that of the single-track laser cladding coating.

    Stress and deformation are closely related to each other.Further,the longitudinal deformation evolution for the overlap ratio of 10%is obtained,as shown in Fig.10. In order to observe the deformation process clearly,the deformation value of the cloud image is magnified by 50 times,and the middle point of one end of the upper surface of the substrate is fixed as the clamping condition of the numerical model(see Fig.10(a)).

    From Figs.9 and 10, we can see that after the first track cladding coating is finished (t =1 s), the substrate presents a slight upward “convex” trend, as shown in Fig. 9(b). In addition, there is large tensile stress inside the coating due to its solidification shrinkage, and compressive stress inside the substrate surrounding the coating due to its thermal expansion limited by the surrounding low-temperature material, as shown in Fig.10(b). At t=2 s,there is an apparent stress relaxation phenomenon of the first track coating due to the reheat effect. During the subsequent cladding process,the stress distribution and stress relaxation phenomena of each track coating are similar to that of the first track coating. However, the tensile stress area increases,and the compressive stress amplitude and area of the surrounding substrate also gradually increase.Meanwhile,the substrate”concave”deformation gradually becomes more serious.

    When the fifth track coatings are finished at t =5 s, the substrate shows an apparent“concave”deformation,as shown in Fig.10(f). Finally,when the substrate is fully cooled down at t=3600 s,it shows longitudinal tensile stresses in the coating area,and the surrounding substrate exists large longitudinal compressive stresses, as shown in Fig. 9(b). The reason is that the tensile stress, caused by the solidification shrinkage of the cladding coatings and the nearby high-temperature materials,leads to the upward bending of the substrate,which will compress the surrounding substrate,thus resulting in compressive stress in the region. At the same time,the transverse shrinkage of the cladding coating will also cause the substrate to bend upward, then finally the substrate is in the“concave”state.

    Fig.9.The location diagram of(a)path 1 at the interface between coating and substrate,and the corresponding longitudinal stress evolution on path 1 for the overlap ratios of(b)10%,(c)30%,and(d)50%.

    Fig.10. For the coating with an overlap ratio of 50%,its stress and deformation evolution varying with time: (a)t=0.1 s,(b)t=1 s,(c)t=2 s,(d)t=3 s,(e)t=4 s,and(f)t=5 s.

    The stress evolution mechanism and the substrate deformation trend for the three overlap ratios are consistent.However, the final longitudinal stress distribution and amplitude of the cladding coating vary with the overlap ratio. At t =3600 s, the maximum longitudinal residual stress (LRS)tends to locate at the middle of the substrate with the overlap ratio increasing. For example, the maximum LRS appears near the fourth cladding coating for the overlap ratio of 10%,while the maximum LRS appears near the third cladding layer for the overlap ratios of 30% and 50%, as Figs. 9(c)and 9(d) show. It is closely related to the heat accumulation during the laser cladding. The heat accumulation region will be the last one to finish the solidification shrinkage process, which leads to the largest residual tensile stress in the region.

    It is worth noting that there is an apparent residual stress mutation at the overlap joint of each track coating for the overlap ratio of 10% (Fig. 9(b)), which may cause geometric defects. In order to analyze the reason, the temperature, stress,and strain evolution of node 3(N3)and middle node(MN)on path 1 are extracted,as Fig.11 shows.

    The difference in stress and strain curves of points N3 and MN is mainly related to their temperature evolution. During the laser cladding process, the highest historical temperature of point N3 is 1185.9?C (larger than the melting point of 1180?C), and it will melt when the temperature exceeds the melting point. While the historical maximum temperature of point MN is 1089.9?C, which means it is in a solid phase throughout(see Fig.12(c)). When the laser heat source passes through the two points, the material at point N3 will melt, its stress reduces to zero instantly, and then the tensile stress will be produced during the subsequent cooling process. However, for the point MN, the compressive stress will be generated due to the restriction of surrounding materials on its high-temperature thermal expansion, and the compressive stress reaches the maximum value when the heat source passes through. Their different stress evolution behaviors finally result in the different residual stress of the two points.

    Fig. 11. The position of (a) the monitoring points N3 and MN, and their corresponding(b)temperature,stress,and strain evolution curves.

    Fig.12. (a)The location diagram,(b)longitudinal residual stress cloud image,and(c)maximum temperature history cloud image of the L-L section for the coating with an overlap ratio of 10%.

    Further, the residual longitudinal stress and maximum temperature history cloud images of the L-L cross-section including N3 and MN points are obtained,as shown in Fig.12.The results show that due to the low overlap ratio,the input energy at the lap joint is insufficient, which makes the lap joint in a compressive stress state,while the rest region of this section is in a tensile stress state. The sudden stress change will lead to geometric defects at the lap joint.[13,43,44]Therefore,the overlap ratio of multi-track cladding coating should not be too low,otherwise unnecessary defects and abnormal residual stress may appear at the lap joint(see Fig.12(b)),resulting in cracks and reducing the fatigue life of the cladding coating.

    4.4. Residual stress and deformation

    Fig.13. The residual stress distributions of VMS,TRS,and LRS for the overlap ratios of 10%,30%,and 50%.

    Fig. 14. The (a) VMS (b) TRS and (c) LRS on path 1 for coating with overlap ratios of 10%,30%,and 50%.

    After the complex temperature, stress, and strain evolution process, the substrate and cladding coating will produce not only residual stresses but also large residual deformations. In order to analyze the residual deformation of the cladding substrate,the deformation of the fixed point(given in Fig. 10(a)) is set to zero. The residual deformation cloud of the cladding substrate with the overlap ratio of 10%is shown in Fig.15(a),and its maximum deformation is about 0.46 mm.At the same time,the residual deformations at the four corners of the cladding substrate for the three different overlap ratios are obtained, as shown in Fig. 15(b). Results show that with the increase of the overlap ratio, the deformations of the four corners decrease, and the deformation trends for all the three overlap ratios are consistent,that is,concave in the center and warping in four corners of the substrate.

    Fig.15. (a)Residual deformation cloud image of the coating with an overlap ratio of 10%and(b)effect of overlap ratios on the residual deformation at the four corners of the cladding substrate.

    Fig.16. The cloud images of(a)maximum temperature history and(b)longitudinal residual stress of the crosssection including path 1 for the coatings of 10%,30%,and 50%.

    High overlap ratio enhances the heat accumulation depth,and makes the tensile residual stress extend to the deeper position of the substrate,as shown in Fig.16,which will balance part of the LRS of the upper surface caused by the coating shrinkage,and decrease the upward warping trend of the substrate,but increase the final VMS and LRS.Therefore,according to the thickness of the substrate and the allowable deformation,the overlap ratio can be increased appropriately to improve the cladding coating quality and to reduce residual deformation and the possibility of cracks. However,a high overlap ratio will increase the heat amount per unit volume, thus enlarging the heat-affected zone and hardening depth of the substrate,and increasing the average residual stress.

    5. Conclusion

    In the paper, the geometric model of high-speed laser cladding coating is established based on the corresponding laser cladding experiment. And a thermal-stress sequential coupling numerical model is developed to reveal the evolution mechanism of temperature, stress, and strain during the laser cladding process and study the influence of the overlap ratio on their evolution rules.

    The temperature,strain,and stress evolution processes of the multi-track cladding coating are extremely complex. During the laser cladding process,with the increase of the overlap ratio,the heat accumulation tends to the center of the cladding coating, which is related to the scanning speed and heat dissipation conditions. Thus, the highest LRS also inclines to locate at the center of the laser cladding coating. When the overlap ratio is 10%, the largest LRS appears near the fourth track cladding coating.For the overlap ratios of 30%and 50%,the maximum residual LRS appears near the third cladding coating.

    For the coating with a small overlap ratio,there are apparent residual stress mutations at the lap joint, which will tend to form cracks. High overlap ratios will slow down stress relaxation, enhance the preheating treatment strength, and reduce the residual deformation. When the overlap ratio is 50%,more than half of each track coating will be re-melted,which will enhance the metallurgical bond strength and reduce the geometric defects. However,increasing the overlap ratio will enlarge the heat-affected zone and hardening depth of the substrate,and also will increase the average residual stress.

    猜你喜歡
    陳平小溪
    與敵人“坦誠相見”\t
    張負的慧眼
    張負的慧眼
    春天的小溪
    民族音樂(2019年6期)2020-01-15 06:50:20
    小溪
    小溪上有橋了
    小學生導刊(2018年1期)2018-03-15 08:02:41
    小溪啊小溪
    草原歌聲(2018年3期)2018-01-24 07:06:52
    陳平過河
    與敵人“坦誠相見”
    慧眼識陳平
    视频在线观看一区二区三区| 欧美日韩亚洲综合一区二区三区_| 亚洲中文字幕日韩| 亚洲精品久久国产高清桃花| 精品乱码久久久久久99久播| 美女大奶头视频| av欧美777| 亚洲国产欧洲综合997久久, | 日韩视频一区二区在线观看| 人成视频在线观看免费观看| 18禁裸乳无遮挡免费网站照片 | 精品高清国产在线一区| 黄色成人免费大全| 亚洲五月色婷婷综合| 神马国产精品三级电影在线观看 | 国产亚洲精品第一综合不卡| 亚洲成人精品中文字幕电影| av中文乱码字幕在线| 亚洲午夜精品一区,二区,三区| 91老司机精品| 淫妇啪啪啪对白视频| www.熟女人妻精品国产| 久久久久精品国产欧美久久久| 欧美一级a爱片免费观看看 | 精品久久久久久久人妻蜜臀av| 免费观看精品视频网站| 中文字幕高清在线视频| 激情在线观看视频在线高清| 丝袜在线中文字幕| 麻豆久久精品国产亚洲av| 日韩精品青青久久久久久| 精华霜和精华液先用哪个| 2021天堂中文幕一二区在线观 | 欧美成狂野欧美在线观看| 黄片大片在线免费观看| 一进一出抽搐动态| 欧美色欧美亚洲另类二区| 亚洲第一欧美日韩一区二区三区| 亚洲欧美精品综合久久99| 少妇的丰满在线观看| 制服人妻中文乱码| 欧美激情久久久久久爽电影| 国产男靠女视频免费网站| 亚洲国产欧洲综合997久久, | 不卡av一区二区三区| 一级黄色大片毛片| 日韩欧美 国产精品| 国产高清激情床上av| 亚洲国产精品999在线| 亚洲av成人av| 一夜夜www| 欧美乱妇无乱码| 欧美中文综合在线视频| 看免费av毛片| 午夜激情av网站| 香蕉久久夜色| 日本三级黄在线观看| 啦啦啦韩国在线观看视频| 免费在线观看视频国产中文字幕亚洲| 国产精品 欧美亚洲| 日本在线视频免费播放| 久久九九热精品免费| 非洲黑人性xxxx精品又粗又长| 欧美人与性动交α欧美精品济南到| 男女之事视频高清在线观看| 国产亚洲av嫩草精品影院| 正在播放国产对白刺激| 黄色 视频免费看| 日本精品一区二区三区蜜桃| 久久午夜亚洲精品久久| 国产又黄又爽又无遮挡在线| 一级黄色大片毛片| 19禁男女啪啪无遮挡网站| 别揉我奶头~嗯~啊~动态视频| 国产99久久九九免费精品| 国产一级毛片七仙女欲春2 | 91在线观看av| 国产又爽黄色视频| 免费人成视频x8x8入口观看| 成人18禁高潮啪啪吃奶动态图| 国产精品98久久久久久宅男小说| 午夜a级毛片| 久久亚洲真实| 亚洲第一电影网av| 白带黄色成豆腐渣| 两个人免费观看高清视频| 国产成人精品久久二区二区91| 狂野欧美激情性xxxx| 国产精品乱码一区二三区的特点| 久久亚洲精品不卡| 制服诱惑二区| 国产精品一区二区免费欧美| 国产伦在线观看视频一区| 久久精品91蜜桃| 伊人久久大香线蕉亚洲五| 午夜免费成人在线视频| 国产一卡二卡三卡精品| 国产真实乱freesex| 成人一区二区视频在线观看| 女人高潮潮喷娇喘18禁视频| 国产高清有码在线观看视频 | 亚洲天堂国产精品一区在线| 亚洲熟女毛片儿| 久久精品aⅴ一区二区三区四区| 久久精品国产亚洲av香蕉五月| 午夜福利成人在线免费观看| 久久久精品国产亚洲av高清涩受| 亚洲av电影不卡..在线观看| 窝窝影院91人妻| 欧美日本视频| 国产成人欧美| 在线永久观看黄色视频| av在线天堂中文字幕| tocl精华| 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产欧美网| 精品国产超薄肉色丝袜足j| av在线播放免费不卡| 国产一区二区激情短视频| 天天添夜夜摸| av福利片在线| 看免费av毛片| 国产一级毛片七仙女欲春2 | 99在线视频只有这里精品首页| 欧美成人免费av一区二区三区| 麻豆一二三区av精品| 国产亚洲av高清不卡| 97碰自拍视频| 久久国产精品男人的天堂亚洲| 国产精品香港三级国产av潘金莲| 久久中文看片网| 亚洲精品中文字幕一二三四区| 无人区码免费观看不卡| 国产精品 欧美亚洲| av电影中文网址| 久久久久久人人人人人| 久久中文看片网| 一级片免费观看大全| 国产v大片淫在线免费观看| 日本一本二区三区精品| 成年版毛片免费区| www日本黄色视频网| 一个人免费在线观看的高清视频| 俺也久久电影网| 亚洲国产中文字幕在线视频| 看黄色毛片网站| 亚洲熟女毛片儿| 久久久久久人人人人人| 亚洲精品美女久久久久99蜜臀| 午夜亚洲福利在线播放| 国产熟女午夜一区二区三区| 高清在线国产一区| 无限看片的www在线观看| 国产黄片美女视频| 嫩草影视91久久| 久久久久久久午夜电影| 黄片大片在线免费观看| 亚洲国产中文字幕在线视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品永久免费网站| 韩国av一区二区三区四区| 欧美亚洲日本最大视频资源| 12—13女人毛片做爰片一| 美女午夜性视频免费| 少妇被粗大的猛进出69影院| 在线永久观看黄色视频| 免费女性裸体啪啪无遮挡网站| 别揉我奶头~嗯~啊~动态视频| 亚洲精品国产一区二区精华液| 黄色丝袜av网址大全| 啦啦啦观看免费观看视频高清| 搡老熟女国产l中国老女人| 波多野结衣av一区二区av| av视频在线观看入口| 法律面前人人平等表现在哪些方面| 亚洲国产精品合色在线| 亚洲全国av大片| 午夜免费观看网址| 搡老妇女老女人老熟妇| 中文字幕精品亚洲无线码一区 | 成熟少妇高潮喷水视频| 国产精品1区2区在线观看.| 亚洲熟妇中文字幕五十中出| 99热6这里只有精品| 成人永久免费在线观看视频| 国产高清激情床上av| 国产熟女午夜一区二区三区| 一本大道久久a久久精品| 色综合亚洲欧美另类图片| 在线观看日韩欧美| 美女国产高潮福利片在线看| 国产精品九九99| 婷婷丁香在线五月| 99久久久亚洲精品蜜臀av| 欧美性猛交╳xxx乱大交人| 亚洲欧美精品综合久久99| 国产精品一区二区精品视频观看| АⅤ资源中文在线天堂| 欧美另类亚洲清纯唯美| 一区二区三区高清视频在线| 日韩欧美免费精品| 亚洲天堂国产精品一区在线| 自线自在国产av| 日本免费a在线| 日韩欧美国产在线观看| 亚洲自拍偷在线| 久久精品91无色码中文字幕| 在线观看舔阴道视频| 久久国产乱子伦精品免费另类| 18美女黄网站色大片免费观看| 精品国产乱子伦一区二区三区| 久久国产乱子伦精品免费另类| 黄片大片在线免费观看| 色老头精品视频在线观看| 亚洲久久久国产精品| 女人被狂操c到高潮| 日本黄色视频三级网站网址| а√天堂www在线а√下载| xxx96com| avwww免费| av福利片在线| 侵犯人妻中文字幕一二三四区| 亚洲无线在线观看| xxx96com| 草草在线视频免费看| 三级毛片av免费| 777久久人妻少妇嫩草av网站| av有码第一页| 国内久久婷婷六月综合欲色啪| 久久久久久人人人人人| 男人的好看免费观看在线视频 | 亚洲自拍偷在线| 一区二区三区精品91| 欧美乱色亚洲激情| 久久久国产精品麻豆| 法律面前人人平等表现在哪些方面| 亚洲专区中文字幕在线| 成人午夜高清在线视频 | 欧美日韩乱码在线| 热re99久久国产66热| 一进一出好大好爽视频| 婷婷丁香在线五月| 女人高潮潮喷娇喘18禁视频| 99精品久久久久人妻精品| 国产精品国产高清国产av| 美女高潮到喷水免费观看| svipshipincom国产片| 大型黄色视频在线免费观看| 国产精华一区二区三区| 777久久人妻少妇嫩草av网站| 亚洲无线在线观看| 国产99久久九九免费精品| 亚洲精品久久成人aⅴ小说| bbb黄色大片| 黄色成人免费大全| 亚洲av熟女| 久久 成人 亚洲| 男女做爰动态图高潮gif福利片| 国内精品久久久久久久电影| 久久婷婷成人综合色麻豆| 人成视频在线观看免费观看| 男男h啪啪无遮挡| 亚洲av中文字字幕乱码综合 | 久久草成人影院| 热re99久久国产66热| 好男人电影高清在线观看| 国产三级黄色录像| a级毛片在线看网站| 欧美国产日韩亚洲一区| 曰老女人黄片| av欧美777| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲中文字幕日韩| 午夜福利在线观看吧| 亚洲av成人一区二区三| 国产精品九九99| 欧美大码av| 国产真人三级小视频在线观看| 中文字幕人妻熟女乱码| 动漫黄色视频在线观看| 1024视频免费在线观看| 国产黄片美女视频| 天堂√8在线中文| 久久精品91无色码中文字幕| 欧美激情高清一区二区三区| 亚洲天堂国产精品一区在线| 一进一出抽搐gif免费好疼| 香蕉丝袜av| 国产高清videossex| 人人妻,人人澡人人爽秒播| 欧美成人一区二区免费高清观看 | 亚洲欧美一区二区三区黑人| 亚洲欧洲精品一区二区精品久久久| ponron亚洲| 51午夜福利影视在线观看| 久久人妻av系列| 亚洲精品一区av在线观看| 亚洲国产欧美一区二区综合| 99国产精品99久久久久| 国产男靠女视频免费网站| 欧美性猛交╳xxx乱大交人| 午夜久久久在线观看| 一区二区三区激情视频| 一个人免费在线观看的高清视频| 欧美国产日韩亚洲一区| 免费在线观看影片大全网站| 欧美日韩中文字幕国产精品一区二区三区| 一个人观看的视频www高清免费观看 | 亚洲男人的天堂狠狠| 欧美不卡视频在线免费观看 | 午夜免费激情av| 曰老女人黄片| 人成视频在线观看免费观看| 又黄又爽又免费观看的视频| av视频在线观看入口| 好看av亚洲va欧美ⅴa在| 国产亚洲欧美在线一区二区| 国产99久久九九免费精品| 久久精品国产清高在天天线| 国产精品永久免费网站| 少妇被粗大的猛进出69影院| 老司机福利观看| 色在线成人网| 国产日本99.免费观看| 免费在线观看影片大全网站| 久久久国产成人免费| 欧美三级亚洲精品| 嫁个100分男人电影在线观看| 国语自产精品视频在线第100页| 精品国产乱码久久久久久男人| 亚洲av片天天在线观看| 欧美性猛交黑人性爽| 亚洲 欧美一区二区三区| 美女高潮喷水抽搐中文字幕| 两个人看的免费小视频| 欧美激情 高清一区二区三区| 亚洲第一电影网av| 香蕉av资源在线| 校园春色视频在线观看| 午夜福利欧美成人| 人妻久久中文字幕网| 2021天堂中文幕一二区在线观 | 美女午夜性视频免费| 成人国语在线视频| 无遮挡黄片免费观看| 美女 人体艺术 gogo| 国产爱豆传媒在线观看 | 国产成人欧美在线观看| 美女国产高潮福利片在线看| 亚洲成av片中文字幕在线观看| 国产97色在线日韩免费| 夜夜看夜夜爽夜夜摸| 欧美激情极品国产一区二区三区| 国产人伦9x9x在线观看| 99国产极品粉嫩在线观看| 国产成人影院久久av| 韩国av一区二区三区四区| 长腿黑丝高跟| 久久午夜综合久久蜜桃| 国产精品香港三级国产av潘金莲| 国产亚洲av嫩草精品影院| 久久精品亚洲精品国产色婷小说| 91麻豆精品激情在线观看国产| 精品一区二区三区四区五区乱码| 国产亚洲精品久久久久久毛片| av视频在线观看入口| 日日爽夜夜爽网站| 欧美在线黄色| 在线观看日韩欧美| 久热爱精品视频在线9| av在线天堂中文字幕| 欧美最黄视频在线播放免费| www日本黄色视频网| 别揉我奶头~嗯~啊~动态视频| 一二三四在线观看免费中文在| av在线天堂中文字幕| 国产精品电影一区二区三区| 免费电影在线观看免费观看| 久久久国产成人精品二区| 每晚都被弄得嗷嗷叫到高潮| 一夜夜www| 精品久久久久久久人妻蜜臀av| 国产精品,欧美在线| 国产97色在线日韩免费| 丝袜在线中文字幕| 亚洲久久久国产精品| 亚洲精品久久国产高清桃花| 亚洲精品国产一区二区精华液| 99国产精品一区二区三区| 人人妻,人人澡人人爽秒播| 国产三级黄色录像| 久9热在线精品视频| 我的亚洲天堂| 精品国产美女av久久久久小说| 午夜亚洲福利在线播放| 美女 人体艺术 gogo| 国产成人av激情在线播放| 老熟妇乱子伦视频在线观看| 久久久久亚洲av毛片大全| 99热只有精品国产| 无人区码免费观看不卡| 欧美三级亚洲精品| a级毛片在线看网站| 黄色毛片三级朝国网站| 国产一区在线观看成人免费| 青草久久国产| 亚洲五月色婷婷综合| 免费观看人在逋| 午夜久久久久精精品| 成人av一区二区三区在线看| 亚洲av片天天在线观看| 国产精品久久久久久精品电影 | 精品少妇一区二区三区视频日本电影| 狠狠狠狠99中文字幕| 成人永久免费在线观看视频| 亚洲精品美女久久久久99蜜臀| 国产国语露脸激情在线看| 好看av亚洲va欧美ⅴa在| 精华霜和精华液先用哪个| 国产成人一区二区三区免费视频网站| 国产亚洲精品综合一区在线观看 | 91麻豆精品激情在线观看国产| 老司机福利观看| 久久欧美精品欧美久久欧美| 大香蕉久久成人网| 国产亚洲精品av在线| 丁香欧美五月| 在线看三级毛片| 老熟妇仑乱视频hdxx| 1024视频免费在线观看| 天天一区二区日本电影三级| 黑人巨大精品欧美一区二区mp4| 最好的美女福利视频网| 国产黄a三级三级三级人| 这个男人来自地球电影免费观看| 婷婷丁香在线五月| 好看av亚洲va欧美ⅴa在| 亚洲全国av大片| 长腿黑丝高跟| 久久性视频一级片| 亚洲欧洲精品一区二区精品久久久| 757午夜福利合集在线观看| 2021天堂中文幕一二区在线观 | 亚洲欧美激情综合另类| 一级毛片女人18水好多| 国产亚洲精品久久久久久毛片| 亚洲片人在线观看| а√天堂www在线а√下载| 欧美黄色淫秽网站| 亚洲av成人一区二区三| 日韩免费av在线播放| 99在线人妻在线中文字幕| 给我免费播放毛片高清在线观看| √禁漫天堂资源中文www| 不卡av一区二区三区| 亚洲自拍偷在线| 亚洲九九香蕉| 免费在线观看亚洲国产| 久久 成人 亚洲| 熟妇人妻久久中文字幕3abv| 国产成人欧美| 欧美黄色片欧美黄色片| 一本久久中文字幕| 黄色片一级片一级黄色片| 欧美国产日韩亚洲一区| 国产熟女xx| 老司机在亚洲福利影院| av欧美777| 老鸭窝网址在线观看| √禁漫天堂资源中文www| 亚洲av五月六月丁香网| 一a级毛片在线观看| 欧美zozozo另类| 熟女少妇亚洲综合色aaa.| 亚洲av五月六月丁香网| 久久精品aⅴ一区二区三区四区| 久久亚洲真实| 天堂影院成人在线观看| 欧美三级亚洲精品| 哪里可以看免费的av片| 国产国语露脸激情在线看| 国产亚洲精品综合一区在线观看 | 欧美三级亚洲精品| 久久久水蜜桃国产精品网| 亚洲人成电影免费在线| 亚洲av中文字字幕乱码综合 | 午夜福利在线观看吧| 久久热在线av| 亚洲精品国产区一区二| 国产精品九九99| 欧美中文综合在线视频| 精华霜和精华液先用哪个| xxxwww97欧美| 男女午夜视频在线观看| 两人在一起打扑克的视频| 久久香蕉精品热| 国产精品二区激情视频| 日本一区二区免费在线视频| 99久久99久久久精品蜜桃| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲黑人精品在线| 欧美久久黑人一区二区| 亚洲专区中文字幕在线| 欧美日韩亚洲综合一区二区三区_| 色综合欧美亚洲国产小说| 精品少妇一区二区三区视频日本电影| 欧美黑人欧美精品刺激| 黄色片一级片一级黄色片| 一本大道久久a久久精品| 国产主播在线观看一区二区| 青草久久国产| 亚洲av片天天在线观看| www日本黄色视频网| 十八禁人妻一区二区| 国产激情偷乱视频一区二区| 村上凉子中文字幕在线| 最近最新免费中文字幕在线| xxx96com| 婷婷精品国产亚洲av在线| 大型av网站在线播放| 50天的宝宝边吃奶边哭怎么回事| 午夜福利在线在线| 一个人观看的视频www高清免费观看 | 久久中文字幕人妻熟女| 法律面前人人平等表现在哪些方面| 999久久久精品免费观看国产| 日本免费a在线| 看片在线看免费视频| 成人免费观看视频高清| 成熟少妇高潮喷水视频| 不卡一级毛片| 日本免费a在线| 九色国产91popny在线| 午夜福利在线在线| 色精品久久人妻99蜜桃| 免费一级毛片在线播放高清视频| 天天躁狠狠躁夜夜躁狠狠躁| 高清毛片免费观看视频网站| 午夜日韩欧美国产| 午夜福利一区二区在线看| 国产真实乱freesex| 国产精品国产高清国产av| 免费高清视频大片| 久久久国产成人精品二区| 动漫黄色视频在线观看| 听说在线观看完整版免费高清| 丰满人妻熟妇乱又伦精品不卡| 伊人久久大香线蕉亚洲五| 亚洲国产欧美一区二区综合| 波多野结衣巨乳人妻| 婷婷六月久久综合丁香| 三级毛片av免费| 精品不卡国产一区二区三区| 黄色视频,在线免费观看| 久久午夜综合久久蜜桃| 国产成人系列免费观看| 熟女少妇亚洲综合色aaa.| 久久国产乱子伦精品免费另类| 久久九九热精品免费| 免费在线观看视频国产中文字幕亚洲| 久久久久久久久久黄片| 可以免费在线观看a视频的电影网站| 亚洲全国av大片| 久久这里只有精品19| 人人妻人人澡欧美一区二区| 在线av久久热| 天天躁狠狠躁夜夜躁狠狠躁| 婷婷丁香在线五月| 免费av毛片视频| 国产成人av激情在线播放| 国产精品二区激情视频| 99国产精品99久久久久| 国产精品精品国产色婷婷| 美女国产高潮福利片在线看| 大香蕉久久成人网| 午夜福利高清视频| 美女 人体艺术 gogo| 18禁黄网站禁片免费观看直播| 国产人伦9x9x在线观看| 长腿黑丝高跟| 9191精品国产免费久久| 一夜夜www| 一级a爱视频在线免费观看| 国产亚洲av嫩草精品影院| 丝袜在线中文字幕| 99热6这里只有精品| 日韩欧美国产一区二区入口| 无人区码免费观看不卡| 精品电影一区二区在线| 在线观看免费午夜福利视频| 日韩欧美 国产精品| 国产精品美女特级片免费视频播放器 | 亚洲自偷自拍图片 自拍| 亚洲国产精品成人综合色| 手机成人av网站| 午夜福利在线观看吧| 亚洲精品中文字幕在线视频| 麻豆成人午夜福利视频| 男人舔女人的私密视频| 国产97色在线日韩免费| 亚洲精品久久成人aⅴ小说| 欧美久久黑人一区二区| 色综合站精品国产| 国产伦在线观看视频一区| 啪啪无遮挡十八禁网站| 欧美日韩黄片免| 欧美色欧美亚洲另类二区| 久久精品aⅴ一区二区三区四区| 中文资源天堂在线| 欧美成狂野欧美在线观看| 日本 av在线| 女人爽到高潮嗷嗷叫在线视频| 麻豆国产av国片精品| 久久精品aⅴ一区二区三区四区| 久久久久免费精品人妻一区二区 |