• <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影院在线观看| 岛国在线观看网站| 久久精品久久久久久噜噜老黄| 欧美黄色片欧美黄色片| 国产精品国产三级国产专区5o| 搡老乐熟女国产| 久久精品成人免费网站| 精品视频人人做人人爽| 黄色a级毛片大全视频| 亚洲五月婷婷丁香| 日韩精品免费视频一区二区三区| 建设人人有责人人尽责人人享有的| 久久ye,这里只有精品| 久久精品国产亚洲av高清一级| 老司机在亚洲福利影院| bbb黄色大片| 老汉色av国产亚洲站长工具| 操出白浆在线播放| av视频免费观看在线观看| 水蜜桃什么品种好| 高清av免费在线| 天天操日日干夜夜撸| 国产精品麻豆人妻色哟哟久久| 成年人午夜在线观看视频| 91大片在线观看| 国产欧美日韩一区二区精品| 19禁男女啪啪无遮挡网站| 男人舔女人的私密视频| 真人做人爱边吃奶动态| tocl精华| 国产精品一区二区在线观看99| 亚洲国产毛片av蜜桃av| 精品少妇一区二区三区视频日本电影| 久久国产亚洲av麻豆专区| 国产精品免费大片| 性高湖久久久久久久久免费观看| 在线 av 中文字幕| 久久国产精品人妻蜜桃| 美女主播在线视频| 大陆偷拍与自拍| 久热这里只有精品99| 成年人午夜在线观看视频| 人妻久久中文字幕网| 1024视频免费在线观看| 宅男免费午夜| 亚洲全国av大片| 欧美精品啪啪一区二区三区 | 欧美亚洲 丝袜 人妻 在线| 久久久久精品人妻al黑| 一个人免费看片子| 国产精品一区二区精品视频观看| 国产精品亚洲av一区麻豆| 又大又爽又粗| 国产在线免费精品| 在线观看免费视频网站a站| 黄色 视频免费看| 黑人巨大精品欧美一区二区mp4| 国产1区2区3区精品| 十分钟在线观看高清视频www| 汤姆久久久久久久影院中文字幕| 欧美日本中文国产一区发布| 精品亚洲成国产av| 国精品久久久久久国模美| 精品熟女少妇八av免费久了| 久久精品国产亚洲av高清一级| 日韩三级视频一区二区三区| 一区二区日韩欧美中文字幕| 热re99久久精品国产66热6| 久久久久久久精品精品| 亚洲欧美成人综合另类久久久| 肉色欧美久久久久久久蜜桃| 91成人精品电影| 啦啦啦啦在线视频资源| 手机成人av网站| 777久久人妻少妇嫩草av网站| av网站在线播放免费| 国产成人精品久久二区二区91| 国产三级黄色录像| 欧美黄色片欧美黄色片| 午夜影院在线不卡| 午夜精品国产一区二区电影| 亚洲一区中文字幕在线| 色播在线永久视频| 精品国产国语对白av| 青春草亚洲视频在线观看| 久久久精品国产亚洲av高清涩受| 亚洲一卡2卡3卡4卡5卡精品中文| 99久久人妻综合| 一级黄色大片毛片| netflix在线观看网站| av又黄又爽大尺度在线免费看| 超色免费av| 成在线人永久免费视频| 日韩欧美一区视频在线观看| 热re99久久精品国产66热6| 女人爽到高潮嗷嗷叫在线视频| av在线播放精品| 一区二区三区乱码不卡18| www日本在线高清视频| 久热爱精品视频在线9| 精品亚洲乱码少妇综合久久| 国产亚洲精品第一综合不卡| 欧美老熟妇乱子伦牲交| 91av网站免费观看| 后天国语完整版免费观看| 午夜激情久久久久久久| 色播在线永久视频| 蜜桃国产av成人99| 亚洲第一av免费看| www.av在线官网国产| 1024视频免费在线观看| 欧美一级毛片孕妇| 老司机福利观看| 一级毛片精品| 国产成人系列免费观看| 1024视频免费在线观看| 亚洲中文日韩欧美视频| 亚洲成人免费av在线播放| 中文字幕人妻丝袜制服| 久久久久久久久免费视频了| 97精品久久久久久久久久精品| 少妇粗大呻吟视频| 日韩免费高清中文字幕av| 老熟妇乱子伦视频在线观看 | 久久人人爽人人片av| av在线app专区| 亚洲精品美女久久久久99蜜臀| 亚洲av日韩在线播放| 久久这里只有精品19| 久久人人爽av亚洲精品天堂| 丝袜脚勾引网站| 极品少妇高潮喷水抽搐| 性高湖久久久久久久久免费观看| 大香蕉久久成人网| 成人三级做爰电影| 王馨瑶露胸无遮挡在线观看| 18禁观看日本| 久久久久精品人妻al黑| 一区二区三区精品91| 亚洲 欧美一区二区三区| 亚洲少妇的诱惑av| 这个男人来自地球电影免费观看| 91精品伊人久久大香线蕉| 91成人精品电影| 中亚洲国语对白在线视频| 精品一区二区三区av网在线观看 | 涩涩av久久男人的天堂| 中亚洲国语对白在线视频| 久久国产精品人妻蜜桃| 国产亚洲av高清不卡| 一本色道久久久久久精品综合| 伊人亚洲综合成人网| 亚洲欧美精品自产自拍| 大型av网站在线播放| 嫁个100分男人电影在线观看| 亚洲综合色网址| 精品国产国语对白av| 久久久久久久大尺度免费视频| 人妻一区二区av| 男女边摸边吃奶| 中文字幕色久视频| 午夜免费鲁丝| 午夜免费鲁丝| 精品少妇一区二区三区视频日本电影| 青春草视频在线免费观看| 国产有黄有色有爽视频| 热99re8久久精品国产| 我的亚洲天堂| 日韩中文字幕欧美一区二区| 国产精品久久久久久精品电影小说| 十八禁人妻一区二区| 丰满迷人的少妇在线观看| 婷婷成人精品国产| 国产不卡av网站在线观看| 亚洲av日韩在线播放| 一级片'在线观看视频| 又紧又爽又黄一区二区| 搡老岳熟女国产| 黑丝袜美女国产一区| 亚洲国产欧美一区二区综合| 操出白浆在线播放| 国产精品一区二区免费欧美 | 999久久久国产精品视频| 1024香蕉在线观看| 免费高清在线观看日韩| 叶爱在线成人免费视频播放| 黄色片一级片一级黄色片| 日日夜夜操网爽| 精品国产超薄肉色丝袜足j| 1024香蕉在线观看| 免费高清在线观看视频在线观看| 美女扒开内裤让男人捅视频| 曰老女人黄片| 日韩精品免费视频一区二区三区| 久久99一区二区三区| 日韩中文字幕视频在线看片| 黑人操中国人逼视频| 亚洲中文字幕日韩| 两个人免费观看高清视频| 亚洲欧美一区二区三区黑人| 亚洲欧洲精品一区二区精品久久久| 中文字幕制服av| 国产精品久久久久久人妻精品电影 | 成年美女黄网站色视频大全免费| 国精品久久久久久国模美| 日韩 欧美 亚洲 中文字幕| 99精品久久久久人妻精品| 色老头精品视频在线观看| 深夜精品福利| 国产男女内射视频| 9热在线视频观看99| 亚洲欧美激情在线| 午夜免费观看性视频| 亚洲一区二区三区欧美精品| 国产成人av激情在线播放| 乱人伦中国视频| 久久久久国内视频| 岛国在线观看网站| 丝袜美腿诱惑在线| 日日爽夜夜爽网站| 久久中文看片网| 在线av久久热| 女人被躁到高潮嗷嗷叫费观| 欧美亚洲 丝袜 人妻 在线| 一区二区三区乱码不卡18| 久久av网站| 欧美日韩福利视频一区二区| 精品少妇久久久久久888优播| 色婷婷久久久亚洲欧美| 激情视频va一区二区三区| www.av在线官网国产| 久久av网站| 亚洲精品美女久久av网站| 老司机午夜十八禁免费视频| 欧美亚洲日本最大视频资源| 中亚洲国语对白在线视频| 欧美日韩国产mv在线观看视频| 精品乱码久久久久久99久播| 高清av免费在线| 亚洲国产成人一精品久久久| 久久九九热精品免费| 我要看黄色一级片免费的| 高清视频免费观看一区二区| 欧美另类亚洲清纯唯美| 无遮挡黄片免费观看| 欧美成狂野欧美在线观看| 一边摸一边抽搐一进一出视频| 国产免费现黄频在线看| 我的亚洲天堂| 日韩欧美一区二区三区在线观看 | 国产精品.久久久| 久久久久久免费高清国产稀缺| 精品久久久久久久毛片微露脸 | 国产真人三级小视频在线观看| 操出白浆在线播放| 丝袜在线中文字幕| a在线观看视频网站| 午夜日韩欧美国产| av在线app专区| 精品福利永久在线观看| 女人爽到高潮嗷嗷叫在线视频| 亚洲成av片中文字幕在线观看| 亚洲欧美成人综合另类久久久| 黑人欧美特级aaaaaa片| 咕卡用的链子| 天天躁夜夜躁狠狠躁躁| 国产成人精品在线电影| 亚洲av男天堂| 老熟女久久久| 黄色片一级片一级黄色片| 亚洲三区欧美一区| 最近最新免费中文字幕在线| 美女国产高潮福利片在线看| 国产免费福利视频在线观看| 午夜91福利影院| 人人妻人人爽人人添夜夜欢视频| 91麻豆av在线| 男女午夜视频在线观看| 99久久国产精品久久久| 欧美 日韩 精品 国产| 电影成人av| 亚洲欧美清纯卡通| 丝袜美腿诱惑在线| 免费少妇av软件| 亚洲精华国产精华精| 夜夜骑夜夜射夜夜干| 国产一区二区激情短视频 | 大香蕉久久网| xxxhd国产人妻xxx| 亚洲久久久国产精品| 日本猛色少妇xxxxx猛交久久| 91成人精品电影| bbb黄色大片| 欧美日韩黄片免| 久久天躁狠狠躁夜夜2o2o| 女人爽到高潮嗷嗷叫在线视频| 国产欧美亚洲国产| 日日夜夜操网爽| 高清视频免费观看一区二区| 成人18禁高潮啪啪吃奶动态图| 国产精品香港三级国产av潘金莲| 午夜激情久久久久久久| 精品国产一区二区三区四区第35| 各种免费的搞黄视频| 国产精品影院久久| 精品国产一区二区久久| 老熟妇仑乱视频hdxx| 久久九九热精品免费| 久久久久国产精品人妻一区二区| 99国产精品一区二区蜜桃av | 91精品伊人久久大香线蕉| 男女床上黄色一级片免费看| 淫妇啪啪啪对白视频 | 国产精品熟女久久久久浪| 女人久久www免费人成看片| 色视频在线一区二区三区| 蜜桃在线观看..| 高清av免费在线| 亚洲欧美色中文字幕在线| 99久久精品国产亚洲精品| 好男人电影高清在线观看| 久久久国产成人免费| 啦啦啦啦在线视频资源| 后天国语完整版免费观看| 亚洲熟女毛片儿| 久久天躁狠狠躁夜夜2o2o| 国产老妇伦熟女老妇高清| 91av网站免费观看| 国产精品熟女久久久久浪| 香蕉丝袜av| 咕卡用的链子| 国产亚洲精品久久久久5区| 国产成人a∨麻豆精品| 十八禁网站免费在线| 久久久久久久久免费视频了| 一二三四社区在线视频社区8| 正在播放国产对白刺激| 国产一卡二卡三卡精品| 国产精品欧美亚洲77777| 精品一区二区三卡| 欧美中文综合在线视频| 亚洲精品久久成人aⅴ小说| 各种免费的搞黄视频| 国产亚洲av片在线观看秒播厂| 日韩一卡2卡3卡4卡2021年| 纯流量卡能插随身wifi吗| 欧美人与性动交α欧美精品济南到| 秋霞在线观看毛片| 久久性视频一级片| 亚洲精品成人av观看孕妇| 精品福利永久在线观看| 97精品久久久久久久久久精品| 美女大奶头黄色视频| 老司机亚洲免费影院| 免费观看人在逋| 午夜91福利影院| 亚洲性夜色夜夜综合| 亚洲欧美一区二区三区久久| 中国国产av一级| 久久狼人影院| 欧美另类一区| a级片在线免费高清观看视频| 伦理电影免费视频| 18禁国产床啪视频网站| 国产精品免费大片| 女人被躁到高潮嗷嗷叫费观| 国产1区2区3区精品| 精品亚洲乱码少妇综合久久| 一级,二级,三级黄色视频| 巨乳人妻的诱惑在线观看| 肉色欧美久久久久久久蜜桃| 亚洲精品久久午夜乱码| 国产精品久久久av美女十八| 免费在线观看黄色视频的| 国产成人免费无遮挡视频| 精品亚洲乱码少妇综合久久| cao死你这个sao货| 高清av免费在线| 欧美精品亚洲一区二区| 国产精品香港三级国产av潘金莲| 美女大奶头黄色视频| 欧美人与性动交α欧美软件| 人人妻人人爽人人添夜夜欢视频| 免费久久久久久久精品成人欧美视频| 丝瓜视频免费看黄片| 亚洲精华国产精华精| 搡老熟女国产l中国老女人| 午夜福利一区二区在线看| 国产97色在线日韩免费| 秋霞在线观看毛片| 国产男女超爽视频在线观看| 亚洲伊人久久精品综合| 另类精品久久| 亚洲欧美一区二区三区黑人| 精品高清国产在线一区| 精品乱码久久久久久99久播| 国产黄频视频在线观看| 中文字幕最新亚洲高清| 超碰97精品在线观看| 国产亚洲欧美在线一区二区| 亚洲精品国产色婷婷电影| 午夜福利视频在线观看免费| 欧美精品一区二区免费开放| 欧美 亚洲 国产 日韩一| 午夜精品国产一区二区电影| 国产又色又爽无遮挡免| 免费观看av网站的网址| 免费在线观看黄色视频的| 一本色道久久久久久精品综合| 国产高清videossex| 国产在视频线精品| 国产高清视频在线播放一区 | 丝袜美腿诱惑在线| 免费观看a级毛片全部| www.精华液| 色播在线永久视频| 热re99久久精品国产66热6| 欧美在线黄色| 久久精品人人爽人人爽视色| 超色免费av| 亚洲精品国产一区二区精华液| 91精品伊人久久大香线蕉| 日韩中文字幕欧美一区二区| 美女午夜性视频免费| 一个人免费在线观看的高清视频 | 丝袜美足系列| 亚洲精品成人av观看孕妇| 少妇精品久久久久久久| 19禁男女啪啪无遮挡网站| 久久久久国内视频| av网站在线播放免费| 午夜91福利影院| 女人被躁到高潮嗷嗷叫费观| kizo精华| 母亲3免费完整高清在线观看| 十分钟在线观看高清视频www| 国产av国产精品国产| 97精品久久久久久久久久精品| av在线老鸭窝| 成人免费观看视频高清| 亚洲熟女精品中文字幕| 国产人伦9x9x在线观看| 9热在线视频观看99| 欧美av亚洲av综合av国产av| 午夜福利免费观看在线| 香蕉国产在线看| 91大片在线观看| av免费在线观看网站| 日韩视频在线欧美| 免费女性裸体啪啪无遮挡网站| 丰满饥渴人妻一区二区三| 极品少妇高潮喷水抽搐| 国产精品.久久久| 99香蕉大伊视频| 日本av手机在线免费观看| 亚洲熟女毛片儿| kizo精华| 久久久久国内视频| 高清在线国产一区| 亚洲专区国产一区二区| 国产xxxxx性猛交| 欧美xxⅹ黑人| av又黄又爽大尺度在线免费看| 午夜福利影视在线免费观看| 免费观看人在逋| 精品久久久久久电影网| 脱女人内裤的视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品美女久久久久99蜜臀| 91国产中文字幕| 我的亚洲天堂| 每晚都被弄得嗷嗷叫到高潮| 日韩,欧美,国产一区二区三区| 久久精品亚洲av国产电影网| 丰满迷人的少妇在线观看| 天堂俺去俺来也www色官网| 久久亚洲国产成人精品v| 国产在线视频一区二区| 国产欧美日韩一区二区精品| 啦啦啦中文免费视频观看日本| 少妇粗大呻吟视频| 国产一区二区激情短视频 | 欧美日韩av久久| 少妇被粗大的猛进出69影院| 久久免费观看电影| 伦理电影免费视频| 天堂俺去俺来也www色官网| 人人妻人人爽人人添夜夜欢视频| 国产色视频综合| 日日夜夜操网爽| 国产精品影院久久| 女人精品久久久久毛片| 黑人操中国人逼视频| 波多野结衣av一区二区av| 51午夜福利影视在线观看| av网站在线播放免费| 少妇 在线观看| 国产91精品成人一区二区三区 | 人人妻人人澡人人爽人人夜夜| 色婷婷av一区二区三区视频| 欧美 亚洲 国产 日韩一| 午夜福利乱码中文字幕| 十八禁人妻一区二区| 亚洲精品久久午夜乱码| 自线自在国产av| 91av网站免费观看| 欧美中文综合在线视频| 激情视频va一区二区三区| 国产三级黄色录像| 欧美黄色淫秽网站| 国产精品99久久99久久久不卡| 久久亚洲国产成人精品v| 久久精品久久久久久噜噜老黄| 亚洲欧美一区二区三区久久| 免费少妇av软件| 亚洲精品国产一区二区精华液| 侵犯人妻中文字幕一二三四区| 美女脱内裤让男人舔精品视频| 黄色毛片三级朝国网站| 一二三四社区在线视频社区8| 国产高清videossex| 天天影视国产精品| 极品人妻少妇av视频| 午夜视频精品福利| 丝袜人妻中文字幕| a级毛片黄视频| 国产精品av久久久久免费| 一级毛片女人18水好多| 亚洲国产日韩一区二区| 熟女少妇亚洲综合色aaa.| 亚洲精品中文字幕在线视频| 日韩欧美免费精品| 丰满迷人的少妇在线观看| 日本五十路高清| 一级片免费观看大全| 一级毛片女人18水好多| 午夜影院在线不卡| 如日韩欧美国产精品一区二区三区| 男男h啪啪无遮挡| a级毛片在线看网站| 久久人人爽av亚洲精品天堂| 97在线人人人人妻| 18在线观看网站| 王馨瑶露胸无遮挡在线观看| 考比视频在线观看| 99热国产这里只有精品6| 曰老女人黄片| 亚洲中文av在线| 国产男女内射视频| 国产淫语在线视频| 亚洲欧洲日产国产| 午夜免费观看性视频| 午夜久久久在线观看| 欧美精品高潮呻吟av久久| 99久久精品国产亚洲精品| 亚洲 国产 在线| 99香蕉大伊视频| 老汉色∧v一级毛片| 亚洲五月婷婷丁香| 97人妻天天添夜夜摸| 纵有疾风起免费观看全集完整版| 亚洲专区字幕在线| 国产成人精品在线电影| 亚洲国产精品一区三区| 天天影视国产精品| 99久久国产精品久久久| 麻豆乱淫一区二区| 亚洲性夜色夜夜综合| 久久人人爽av亚洲精品天堂| 巨乳人妻的诱惑在线观看| 日本av手机在线免费观看| 丝瓜视频免费看黄片| 中国美女看黄片| 无限看片的www在线观看| 精品国产一区二区三区四区第35| 欧美精品一区二区免费开放| 国产成人欧美| 亚洲一码二码三码区别大吗| 法律面前人人平等表现在哪些方面 | 秋霞在线观看毛片| 叶爱在线成人免费视频播放| 亚洲 欧美一区二区三区| 国产免费av片在线观看野外av| 十分钟在线观看高清视频www| 成人av一区二区三区在线看 | 午夜福利视频精品| 久久av网站| 亚洲avbb在线观看| 久久久国产欧美日韩av| tocl精华| 一边摸一边做爽爽视频免费| 国产免费福利视频在线观看| netflix在线观看网站| 国产在线视频一区二区| 搡老乐熟女国产| 天天影视国产精品| 亚洲av男天堂| 成年av动漫网址| 波多野结衣一区麻豆| 国产亚洲av高清不卡| 人成视频在线观看免费观看| 国产亚洲精品一区二区www | 夜夜夜夜夜久久久久| 国产无遮挡羞羞视频在线观看| av超薄肉色丝袜交足视频| 热99国产精品久久久久久7| 亚洲中文av在线| 久久99热这里只频精品6学生| 黄片小视频在线播放| √禁漫天堂资源中文www| 欧美在线黄色| 各种免费的搞黄视频| 精品人妻一区二区三区麻豆| 热99re8久久精品国产|