• <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
    陳平過河
    與敵人“坦誠相見”
    慧眼識陳平
    变态另类成人亚洲欧美熟女 | 国产成人影院久久av| 老司机午夜十八禁免费视频| 怎么达到女性高潮| 欧美成人午夜精品| 天天躁狠狠躁夜夜躁狠狠躁| 天天添夜夜摸| 69精品国产乱码久久久| 亚洲国产精品久久男人天堂| 成人国产综合亚洲| 99国产极品粉嫩在线观看| 一区福利在线观看| 狠狠狠狠99中文字幕| e午夜精品久久久久久久| 日日摸夜夜添夜夜添小说| 亚洲精品在线美女| 免费在线观看完整版高清| 亚洲av熟女| www国产在线视频色| 亚洲国产欧美网| 久久久久久人人人人人| 99国产精品免费福利视频| a在线观看视频网站| 多毛熟女@视频| 色老头精品视频在线观看| 日韩成人在线观看一区二区三区| 国产精品爽爽va在线观看网站 | 首页视频小说图片口味搜索| 女人爽到高潮嗷嗷叫在线视频| 欧美+亚洲+日韩+国产| 国产一区二区三区综合在线观看| 亚洲色图 男人天堂 中文字幕| 制服人妻中文乱码| 国产亚洲av高清不卡| 两个人免费观看高清视频| а√天堂www在线а√下载| 在线av久久热| 又黄又粗又硬又大视频| 欧美av亚洲av综合av国产av| 免费看美女性在线毛片视频| 久久这里只有精品19| 丁香六月欧美| 亚洲全国av大片| 啦啦啦 在线观看视频| 国产av又大| 精品国产乱子伦一区二区三区| 精品国产亚洲在线| 亚洲久久久国产精品| 大香蕉久久成人网| 久久人人爽av亚洲精品天堂| 日韩中文字幕欧美一区二区| 黄片小视频在线播放| 大香蕉久久成人网| 免费观看精品视频网站| 一区二区三区高清视频在线| 久久久水蜜桃国产精品网| 大香蕉久久成人网| 国内精品久久久久精免费| 精品少妇一区二区三区视频日本电影| 国产av一区二区精品久久| 一级,二级,三级黄色视频| 午夜成年电影在线免费观看| 一级a爱视频在线免费观看| 两个人免费观看高清视频| 色精品久久人妻99蜜桃| 亚洲天堂国产精品一区在线| 一二三四在线观看免费中文在| 国产主播在线观看一区二区| 18美女黄网站色大片免费观看| 国内久久婷婷六月综合欲色啪| 黄色成人免费大全| 欧美老熟妇乱子伦牲交| 欧美成狂野欧美在线观看| 欧美一级a爱片免费观看看 | 国产亚洲精品久久久久久毛片| 少妇被粗大的猛进出69影院| 色精品久久人妻99蜜桃| 亚洲情色 制服丝袜| 亚洲一区二区三区色噜噜| 亚洲久久久国产精品| 久久久国产成人免费| 啦啦啦韩国在线观看视频| ponron亚洲| 成人亚洲精品一区在线观看| 桃色一区二区三区在线观看| 不卡av一区二区三区| 日韩 欧美 亚洲 中文字幕| 一级毛片高清免费大全| 黑人巨大精品欧美一区二区蜜桃| 免费观看人在逋| 中文字幕色久视频| 手机成人av网站| 午夜福利一区二区在线看| 丝袜人妻中文字幕| 脱女人内裤的视频| 亚洲av成人不卡在线观看播放网| 99国产综合亚洲精品| 成年女人毛片免费观看观看9| 久久久久国产精品人妻aⅴ院| avwww免费| 色精品久久人妻99蜜桃| 欧美一级a爱片免费观看看 | 欧美亚洲日本最大视频资源| 久久久国产精品麻豆| av片东京热男人的天堂| 少妇熟女aⅴ在线视频| 日韩欧美一区视频在线观看| 午夜免费激情av| 一级毛片女人18水好多| 免费看a级黄色片| 两个人视频免费观看高清| 免费久久久久久久精品成人欧美视频| 国产亚洲av高清不卡| 一二三四社区在线视频社区8| 欧美日韩一级在线毛片| 91九色精品人成在线观看| 巨乳人妻的诱惑在线观看| 啪啪无遮挡十八禁网站| 夜夜躁狠狠躁天天躁| 男人舔女人的私密视频| 免费看美女性在线毛片视频| 多毛熟女@视频| 欧美日本亚洲视频在线播放| 日韩视频一区二区在线观看| 久久久久国产一级毛片高清牌| 777久久人妻少妇嫩草av网站| 亚洲成人国产一区在线观看| av福利片在线| 欧美乱色亚洲激情| 日韩欧美一区二区三区在线观看| 国产av又大| 午夜日韩欧美国产| 免费少妇av软件| 亚洲精品久久国产高清桃花| 黄色a级毛片大全视频| 大型黄色视频在线免费观看| 搡老妇女老女人老熟妇| 亚洲天堂国产精品一区在线| 97碰自拍视频| 久久久久久久久久久久大奶| 精品一区二区三区av网在线观看| 波多野结衣高清无吗| 亚洲av第一区精品v没综合| 丰满人妻熟妇乱又伦精品不卡| 国产一级毛片七仙女欲春2 | 亚洲人成伊人成综合网2020| 一边摸一边抽搐一进一出视频| 国产私拍福利视频在线观看| 国产国语露脸激情在线看| 久久精品亚洲精品国产色婷小说| 国产欧美日韩一区二区三| 99在线视频只有这里精品首页| 亚洲av成人一区二区三| 精品一品国产午夜福利视频| 麻豆国产av国片精品| 精品国产国语对白av| av欧美777| 久久午夜亚洲精品久久| 国产成+人综合+亚洲专区| 日韩欧美一区视频在线观看| 可以在线观看的亚洲视频| 一本综合久久免费| 欧美激情 高清一区二区三区| 久99久视频精品免费| 欧美绝顶高潮抽搐喷水| 国产精品,欧美在线| 人人妻人人澡人人看| 少妇熟女aⅴ在线视频| 美女高潮喷水抽搐中文字幕| 最新在线观看一区二区三区| 黑人巨大精品欧美一区二区mp4| 日韩三级视频一区二区三区| 可以免费在线观看a视频的电影网站| 在线观看66精品国产| 啦啦啦 在线观看视频| www国产在线视频色| 亚洲中文字幕日韩| 国产精品自产拍在线观看55亚洲| 成人18禁在线播放| 午夜精品国产一区二区电影| 狠狠狠狠99中文字幕| 国产精品一区二区免费欧美| 午夜影院日韩av| a在线观看视频网站| 99久久99久久久精品蜜桃| avwww免费| 久久香蕉激情| 成年版毛片免费区| 久久久久久久午夜电影| √禁漫天堂资源中文www| 亚洲成人久久性| 免费看美女性在线毛片视频| 精品一区二区三区av网在线观看| 丁香欧美五月| 男人的好看免费观看在线视频 | 国产精品一区二区免费欧美| 麻豆久久精品国产亚洲av| 国产精品影院久久| 成熟少妇高潮喷水视频| 亚洲avbb在线观看| 可以在线观看毛片的网站| 午夜日韩欧美国产| 国产高清视频在线播放一区| 婷婷六月久久综合丁香| 男人的好看免费观看在线视频 | 亚洲熟妇中文字幕五十中出| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩大尺度精品在线看网址 | av福利片在线| 一级毛片女人18水好多| 黄色成人免费大全| 久久精品91蜜桃| 日本欧美视频一区| 欧美激情高清一区二区三区| 欧美日韩一级在线毛片| 亚洲第一电影网av| 欧美最黄视频在线播放免费| 又紧又爽又黄一区二区| 成人国产一区最新在线观看| 一a级毛片在线观看| 国产黄a三级三级三级人| 久久精品91无色码中文字幕| 极品人妻少妇av视频| 怎么达到女性高潮| 岛国在线观看网站| 精品国产一区二区三区四区第35| 久久久久久久精品吃奶| 精品久久久久久久人妻蜜臀av | 久久这里只有精品19| 少妇 在线观看| 黑丝袜美女国产一区| 黑人欧美特级aaaaaa片| 国产xxxxx性猛交| 99香蕉大伊视频| 亚洲午夜精品一区,二区,三区| 美女高潮到喷水免费观看| 亚洲国产中文字幕在线视频| 亚洲伊人色综图| 国产成人影院久久av| 久久精品亚洲精品国产色婷小说| 黄色视频不卡| 国产aⅴ精品一区二区三区波| 涩涩av久久男人的天堂| 欧美色视频一区免费| 岛国在线观看网站| 成在线人永久免费视频| 国产精品一区二区免费欧美| 亚洲国产看品久久| 精品一区二区三区四区五区乱码| 中文亚洲av片在线观看爽| 久久人妻福利社区极品人妻图片| 免费观看精品视频网站| 一进一出好大好爽视频| 性欧美人与动物交配| 老司机靠b影院| 69av精品久久久久久| 成人欧美大片| 丰满的人妻完整版| 欧美黄色淫秽网站| 少妇粗大呻吟视频| 日韩高清综合在线| 一边摸一边抽搐一进一出视频| 亚洲中文av在线| 人妻久久中文字幕网| 女同久久另类99精品国产91| 一夜夜www| 国产1区2区3区精品| 日韩欧美一区视频在线观看| 免费不卡黄色视频| 999久久久国产精品视频| 中国美女看黄片| 成人永久免费在线观看视频| 亚洲无线在线观看| 精品一区二区三区av网在线观看| 婷婷六月久久综合丁香| 色综合欧美亚洲国产小说| 男人舔女人的私密视频| 成人手机av| 超碰成人久久| or卡值多少钱| 国产精品电影一区二区三区| 国产精品,欧美在线| 性欧美人与动物交配| 大型黄色视频在线免费观看| 99在线视频只有这里精品首页| 一级a爱视频在线免费观看| 黄频高清免费视频| 亚洲国产日韩欧美精品在线观看 | 好男人在线观看高清免费视频 | 亚洲成人国产一区在线观看| 校园春色视频在线观看| 99久久久亚洲精品蜜臀av| 国产精品 国内视频| netflix在线观看网站| 亚洲av成人av| 成人18禁高潮啪啪吃奶动态图| 啦啦啦免费观看视频1| 国产国语露脸激情在线看| 不卡一级毛片| 成人手机av| 777久久人妻少妇嫩草av网站| 欧美日韩乱码在线| 精品高清国产在线一区| 欧美日本中文国产一区发布| 欧美久久黑人一区二区| 1024香蕉在线观看| 日韩欧美在线二视频| 亚洲自拍偷在线| 女警被强在线播放| 丝袜在线中文字幕| 国产av又大| 亚洲熟女毛片儿| 两个人视频免费观看高清| 久久精品国产综合久久久| 啦啦啦 在线观看视频| 亚洲情色 制服丝袜| 国产精品久久久久久精品电影 | 亚洲最大成人中文| 国产伦人伦偷精品视频| 美国免费a级毛片| 亚洲中文av在线| 看免费av毛片| 国产精品,欧美在线| 亚洲国产精品久久男人天堂| 日韩 欧美 亚洲 中文字幕| 久久久国产成人精品二区| 亚洲av成人不卡在线观看播放网| 热99re8久久精品国产| 一区二区三区高清视频在线| 麻豆国产av国片精品| 动漫黄色视频在线观看| 午夜免费激情av| 一级黄色大片毛片| 亚洲精品久久国产高清桃花| 精品少妇一区二区三区视频日本电影| 视频在线观看一区二区三区| 亚洲欧美日韩另类电影网站| 1024视频免费在线观看| 精品不卡国产一区二区三区| www.www免费av| 黄色成人免费大全| 免费看十八禁软件| 在线观看免费视频日本深夜| 日韩高清综合在线| 黄网站色视频无遮挡免费观看| 久久久久国内视频| 中文亚洲av片在线观看爽| 中国美女看黄片| 国语自产精品视频在线第100页| 亚洲av熟女| 欧美日韩瑟瑟在线播放| 满18在线观看网站| 久久久久国产精品人妻aⅴ院| 男女之事视频高清在线观看| 久久久久国产精品人妻aⅴ院| 亚洲七黄色美女视频| 超碰成人久久| 亚洲aⅴ乱码一区二区在线播放 | 91在线观看av| 国产精品av久久久久免费| 国产成人精品久久二区二区免费| 91成年电影在线观看| 免费一级毛片在线播放高清视频 | 在线十欧美十亚洲十日本专区| 久久人妻福利社区极品人妻图片| 日韩国内少妇激情av| 国内毛片毛片毛片毛片毛片| 一本大道久久a久久精品| 久久人妻福利社区极品人妻图片| 久久精品亚洲熟妇少妇任你| 大型av网站在线播放| 欧美日韩瑟瑟在线播放| 亚洲五月天丁香| 97人妻天天添夜夜摸| 亚洲国产精品成人综合色| 免费无遮挡裸体视频| 国产精品 国内视频| 国产精品免费视频内射| 中国美女看黄片| 国产精品亚洲一级av第二区| 欧美成狂野欧美在线观看| 欧美亚洲日本最大视频资源| 99香蕉大伊视频| 999久久久国产精品视频| 黄色片一级片一级黄色片| www.自偷自拍.com| 女人被狂操c到高潮| 成人三级黄色视频| 人妻久久中文字幕网| 一区福利在线观看| 成人18禁高潮啪啪吃奶动态图| 99热只有精品国产| 中文字幕av电影在线播放| 国产一区二区激情短视频| 亚洲中文日韩欧美视频| 亚洲九九香蕉| 欧美激情久久久久久爽电影 | 久久婷婷成人综合色麻豆| 欧美人与性动交α欧美精品济南到| 午夜视频精品福利| 很黄的视频免费| 欧美日本亚洲视频在线播放| 亚洲av电影不卡..在线观看| 他把我摸到了高潮在线观看| 九色国产91popny在线| 亚洲精品美女久久久久99蜜臀| 免费在线观看日本一区| 国产一区二区三区在线臀色熟女| 老司机深夜福利视频在线观看| av视频在线观看入口| 九色亚洲精品在线播放| 国产亚洲av嫩草精品影院| 99久久久亚洲精品蜜臀av| 无人区码免费观看不卡| 亚洲成人精品中文字幕电影| 亚洲男人的天堂狠狠| 免费av毛片视频| 91老司机精品| 国产精品爽爽va在线观看网站 | 亚洲中文字幕日韩| 可以在线观看的亚洲视频| 国产av在哪里看| 精品欧美一区二区三区在线| 两个人免费观看高清视频| 九色国产91popny在线| 欧美日本中文国产一区发布| 欧美乱码精品一区二区三区| 男女下面插进去视频免费观看| 亚洲专区国产一区二区| 波多野结衣一区麻豆| 黑人巨大精品欧美一区二区蜜桃| 精品欧美一区二区三区在线| 国产私拍福利视频在线观看| 黄色毛片三级朝国网站| 欧美日韩乱码在线| 亚洲欧美日韩另类电影网站| 午夜精品久久久久久毛片777| 欧美+亚洲+日韩+国产| 女人被狂操c到高潮| 日本 av在线| 欧美久久黑人一区二区| 91成年电影在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 国产1区2区3区精品| 叶爱在线成人免费视频播放| 亚洲精品一卡2卡三卡4卡5卡| 国产色视频综合| 欧美一区二区精品小视频在线| 欧美日韩一级在线毛片| 国产av在哪里看| 色av中文字幕| 一级a爱片免费观看的视频| 精品午夜福利视频在线观看一区| 欧美老熟妇乱子伦牲交| 亚洲中文日韩欧美视频| 满18在线观看网站| 在线观看www视频免费| 免费在线观看黄色视频的| 午夜影院日韩av| 亚洲国产精品久久男人天堂| 老汉色av国产亚洲站长工具| 无限看片的www在线观看| 母亲3免费完整高清在线观看| 精品一区二区三区视频在线观看免费| 亚洲欧美激情综合另类| 午夜福利欧美成人| 91老司机精品| 丝袜美腿诱惑在线| 麻豆成人av在线观看| 男人操女人黄网站| 波多野结衣av一区二区av| 国产精品 欧美亚洲| 亚洲av成人av| 欧美绝顶高潮抽搐喷水| 桃红色精品国产亚洲av| av天堂在线播放| 欧洲精品卡2卡3卡4卡5卡区| 国产精品二区激情视频| 国产精品亚洲美女久久久| 久久天躁狠狠躁夜夜2o2o| 妹子高潮喷水视频| 国产97色在线日韩免费| 久久天堂一区二区三区四区| 男男h啪啪无遮挡| av视频在线观看入口| 亚洲 欧美 日韩 在线 免费| 精品无人区乱码1区二区| 男人舔女人下体高潮全视频| 免费在线观看日本一区| 18禁美女被吸乳视频| 日韩中文字幕欧美一区二区| 欧美色欧美亚洲另类二区 | 日日夜夜操网爽| 精品国产乱码久久久久久男人| 一级a爱视频在线免费观看| 99香蕉大伊视频| 搡老岳熟女国产| 久久热在线av| 亚洲九九香蕉| 亚洲三区欧美一区| 夜夜夜夜夜久久久久| 一区二区日韩欧美中文字幕| 可以在线观看毛片的网站| 亚洲欧美一区二区三区黑人| 亚洲精品国产区一区二| 国产成人欧美在线观看| 岛国视频午夜一区免费看| 丝袜美足系列| 电影成人av| 精品电影一区二区在线| 美女扒开内裤让男人捅视频| АⅤ资源中文在线天堂| 国产精品乱码一区二三区的特点 | 国产男靠女视频免费网站| 丰满人妻熟妇乱又伦精品不卡| 亚洲中文av在线| 曰老女人黄片| 日本一区二区免费在线视频| 久久精品aⅴ一区二区三区四区| 成人欧美大片| 母亲3免费完整高清在线观看| 老熟妇仑乱视频hdxx| 午夜福利成人在线免费观看| 国产一区二区三区综合在线观看| 亚洲最大成人中文| 欧美人与性动交α欧美精品济南到| 国产在线观看jvid| 怎么达到女性高潮| 高清黄色对白视频在线免费看| 999精品在线视频| 中出人妻视频一区二区| 久久精品影院6| 757午夜福利合集在线观看| 欧美色视频一区免费| 国产又爽黄色视频| 国产激情欧美一区二区| 夜夜爽天天搞| www.自偷自拍.com| 亚洲欧洲精品一区二区精品久久久| 国产成人精品在线电影| 一二三四在线观看免费中文在| 国产成人av教育| 精品一品国产午夜福利视频| 波多野结衣巨乳人妻| 老鸭窝网址在线观看| 午夜精品国产一区二区电影| 日本三级黄在线观看| 久久久国产精品麻豆| 变态另类成人亚洲欧美熟女 | av视频免费观看在线观看| 色播在线永久视频| 视频在线观看一区二区三区| 免费看a级黄色片| 两个人看的免费小视频| 久久久精品欧美日韩精品| 亚洲欧美日韩另类电影网站| 操美女的视频在线观看| 国产精品,欧美在线| 成人永久免费在线观看视频| 精品欧美国产一区二区三| 亚洲第一av免费看| 亚洲精品av麻豆狂野| 久久久久久久久久久久大奶| 乱人伦中国视频| 久久午夜综合久久蜜桃| 久久人妻熟女aⅴ| 成人av一区二区三区在线看| 法律面前人人平等表现在哪些方面| 国产三级黄色录像| 一级黄色大片毛片| 成人国语在线视频| 天堂√8在线中文| 国语自产精品视频在线第100页| 非洲黑人性xxxx精品又粗又长| 色综合亚洲欧美另类图片| 美国免费a级毛片| tocl精华| av在线播放免费不卡| 午夜免费鲁丝| 久久久久久久久久久久大奶| 国产av在哪里看| 我的亚洲天堂| 久久久久久久久久久久大奶| 婷婷精品国产亚洲av在线| 人人澡人人妻人| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品中文字幕在线视频| 欧美成人免费av一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品久久成人aⅴ小说| 看黄色毛片网站| 亚洲人成电影观看| 欧美av亚洲av综合av国产av| 日日爽夜夜爽网站| 丝袜人妻中文字幕| 美女大奶头视频| 一区二区三区国产精品乱码| 十八禁人妻一区二区| 亚洲欧美日韩无卡精品| 久久 成人 亚洲| 桃红色精品国产亚洲av| 一a级毛片在线观看| netflix在线观看网站| 欧美色视频一区免费| 黄色女人牲交| 超碰成人久久| 午夜老司机福利片| 男女床上黄色一级片免费看| 亚洲自偷自拍图片 自拍| 一进一出抽搐动态| 免费av毛片视频| www日本在线高清视频| bbb黄色大片| 88av欧美| 亚洲欧美精品综合一区二区三区|