• <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麻豆| 综合色丁香网| 少妇熟女欧美另类| 日韩欧美在线乱码| 日韩av在线免费看完整版不卡| 久久久久久久久久黄片| 午夜福利视频1000在线观看| 久久久a久久爽久久v久久| 国产白丝娇喘喷水9色精品| av播播在线观看一区| 国产精品三级大全| 国产真实乱freesex| 国产精品麻豆人妻色哟哟久久 | 两个人的视频大全免费| 久久精品影院6| 男人的好看免费观看在线视频| 国产精品伦人一区二区| 免费看光身美女| 久久99蜜桃精品久久| 亚洲精品国产成人久久av| 99在线视频只有这里精品首页| 午夜久久久久精精品| 成年版毛片免费区| 色综合亚洲欧美另类图片| 女人十人毛片免费观看3o分钟| 啦啦啦啦在线视频资源| 精品少妇黑人巨大在线播放 | 内地一区二区视频在线| 一个人看的www免费观看视频| www日本黄色视频网| 免费观看人在逋| 人妻系列 视频| 精品国内亚洲2022精品成人| 中文字幕av成人在线电影| 99在线人妻在线中文字幕| 日本一二三区视频观看| 国产色婷婷99| 日韩一本色道免费dvd| 高清av免费在线| 美女cb高潮喷水在线观看| 尤物成人国产欧美一区二区三区| 夜夜看夜夜爽夜夜摸| 日韩制服骚丝袜av| 免费观看的影片在线观看| 高清在线视频一区二区三区 | 亚洲av一区综合| 国产精品国产三级国产专区5o | 99久久无色码亚洲精品果冻| 日韩欧美国产在线观看| 身体一侧抽搐| 综合色av麻豆| 极品教师在线视频| 免费av不卡在线播放| 精品国产三级普通话版| 日韩av在线大香蕉| 日日摸夜夜添夜夜爱| 国产午夜精品论理片| 97超碰精品成人国产| 亚洲经典国产精华液单| 日韩视频在线欧美| av线在线观看网站| 日韩亚洲欧美综合| 久久草成人影院| 国产精品99久久久久久久久| 91av网一区二区| 日韩成人av中文字幕在线观看| 日本黄色片子视频| 伦理电影大哥的女人| 26uuu在线亚洲综合色| 欧美激情在线99| 精品99又大又爽又粗少妇毛片| 黄色一级大片看看| 我的女老师完整版在线观看| 国产美女午夜福利| 少妇熟女aⅴ在线视频| 亚洲国产精品sss在线观看| 91在线精品国自产拍蜜月| 91精品国产九色| 黄色配什么色好看| 在线观看一区二区三区| 日本欧美国产在线视频| 大又大粗又爽又黄少妇毛片口| 两个人的视频大全免费| 又粗又硬又长又爽又黄的视频| 九九久久精品国产亚洲av麻豆| 国产91av在线免费观看| 亚洲精品456在线播放app| 日韩亚洲欧美综合| 免费电影在线观看免费观看| 久久99精品国语久久久| 国产午夜精品久久久久久一区二区三区| 成人漫画全彩无遮挡| 成人二区视频| 国产精品.久久久| 久久这里只有精品中国| 欧美色视频一区免费| 小蜜桃在线观看免费完整版高清| 哪个播放器可以免费观看大片| 最近中文字幕高清免费大全6| 国产单亲对白刺激| 日本免费在线观看一区| 一级二级三级毛片免费看| 啦啦啦啦在线视频资源| 美女脱内裤让男人舔精品视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日韩欧美精品v在线| 国产黄色小视频在线观看| 麻豆av噜噜一区二区三区| 九草在线视频观看| 国产成人精品婷婷| 亚洲人与动物交配视频| 国产精品.久久久| 99久久中文字幕三级久久日本| 久久精品久久精品一区二区三区| АⅤ资源中文在线天堂| 国产极品精品免费视频能看的| 全区人妻精品视频| 亚洲第一区二区三区不卡| 久久人妻av系列| 蜜桃亚洲精品一区二区三区| 乱码一卡2卡4卡精品| 波多野结衣巨乳人妻| 大话2 男鬼变身卡| 精品免费久久久久久久清纯| 内地一区二区视频在线| 国产乱人视频| 色播亚洲综合网| 又爽又黄a免费视频| 天堂中文最新版在线下载 | 亚洲四区av| 人妻夜夜爽99麻豆av| 91狼人影院| 国产精品国产av在线观看| 亚洲第一区二区三区不卡| 搡女人真爽免费视频火全软件| 日日爽夜夜爽网站| 纯流量卡能插随身wifi吗| 在线观看人妻少妇| 少妇的逼水好多| 亚洲,欧美,日韩| 欧美精品一区二区大全| 在线观看免费日韩欧美大片| 97在线视频观看| 久久99热这里只频精品6学生| 少妇人妻精品综合一区二区| 久久这里有精品视频免费| 婷婷色麻豆天堂久久| 多毛熟女@视频| av卡一久久| 久久人人爽人人片av| 久久精品久久久久久久性| 精品一区二区三区视频在线| 人妻 亚洲 视频| 黄色怎么调成土黄色| 一本久久精品| 黄网站色视频无遮挡免费观看| 丝袜喷水一区| 亚洲国产日韩一区二区| 色哟哟·www| 在线观看www视频免费| 欧美最新免费一区二区三区| 两个人看的免费小视频| 国产亚洲av片在线观看秒播厂| 中文天堂在线官网| 在线观看免费高清a一片| 涩涩av久久男人的天堂| 91aial.com中文字幕在线观看| 亚洲成人手机| 色5月婷婷丁香| 国产成人91sexporn| 九色亚洲精品在线播放| 最近中文字幕2019免费版| 国产男女内射视频| 黑人巨大精品欧美一区二区蜜桃 | 人妻人人澡人人爽人人| 免费看光身美女| 色吧在线观看| 晚上一个人看的免费电影| 香蕉丝袜av| 国产亚洲最大av| 国产爽快片一区二区三区| 永久免费av网站大全| 最近中文字幕高清免费大全6| 在线观看人妻少妇| 久久精品国产自在天天线| 美女国产高潮福利片在线看| 性高湖久久久久久久久免费观看| 久久久国产一区二区| 飞空精品影院首页| 亚洲av中文av极速乱| 久久 成人 亚洲| 国产色婷婷99| 国产成人91sexporn| 亚洲精品成人av观看孕妇| 少妇高潮的动态图| 精品第一国产精品| 国产色婷婷99| 极品人妻少妇av视频| 欧美精品亚洲一区二区| 国产精品.久久久| 午夜久久久在线观看| 丝袜在线中文字幕| 天天操日日干夜夜撸| 久久人人爽av亚洲精品天堂| 一边摸一边做爽爽视频免费| 99热全是精品| 久久午夜综合久久蜜桃| 精品卡一卡二卡四卡免费| 亚洲欧美日韩另类电影网站| 日韩制服丝袜自拍偷拍| 人妻 亚洲 视频| 国产免费福利视频在线观看| 亚洲一级一片aⅴ在线观看| 欧美日韩av久久| 欧美激情极品国产一区二区三区 | 国产免费又黄又爽又色| 亚洲国产精品专区欧美| 欧美成人精品欧美一级黄| 国产成人精品久久久久久| 午夜免费鲁丝| 97超碰精品成人国产| 日日撸夜夜添| 免费黄网站久久成人精品| 中文天堂在线官网| videossex国产| 欧美激情 高清一区二区三区| 国产成人免费观看mmmm| 国产爽快片一区二区三区| 1024视频免费在线观看| 性色avwww在线观看| 熟女av电影| 黑人高潮一二区| 亚洲精品456在线播放app| 中文天堂在线官网| 99久国产av精品国产电影| 免费大片黄手机在线观看| 日本猛色少妇xxxxx猛交久久| 黑人巨大精品欧美一区二区蜜桃 | 久久精品国产鲁丝片午夜精品| 汤姆久久久久久久影院中文字幕| 下体分泌物呈黄色| 大陆偷拍与自拍| 成人18禁高潮啪啪吃奶动态图| 又黄又爽又刺激的免费视频.| 观看av在线不卡| 丝袜脚勾引网站| 大陆偷拍与自拍| 黄色怎么调成土黄色| 久久精品aⅴ一区二区三区四区 | 少妇的逼好多水| 高清毛片免费看| 汤姆久久久久久久影院中文字幕| 日韩 亚洲 欧美在线| 欧美最新免费一区二区三区| 一边亲一边摸免费视频| 秋霞在线观看毛片| h视频一区二区三区| 亚洲中文av在线| 少妇人妻 视频| 最近中文字幕高清免费大全6| 在线亚洲精品国产二区图片欧美| 免费在线观看完整版高清| 日本黄大片高清| 亚洲av福利一区| 好男人视频免费观看在线| 亚洲精品456在线播放app| 成人综合一区亚洲| 亚洲av欧美aⅴ国产| 国产精品一国产av| 久久精品久久久久久噜噜老黄| 啦啦啦中文免费视频观看日本| 少妇 在线观看| 国产成人一区二区在线| 国产麻豆69| av女优亚洲男人天堂| 久久久a久久爽久久v久久| 麻豆乱淫一区二区| 人妻 亚洲 视频| 国产熟女欧美一区二区| 国产成人免费观看mmmm| xxxhd国产人妻xxx| 欧美成人午夜免费资源| 精品人妻熟女毛片av久久网站| 中文字幕另类日韩欧美亚洲嫩草| 国产一区二区在线观看av| 精品久久久精品久久久| 亚洲精品中文字幕在线视频| 欧美激情极品国产一区二区三区 | 老司机影院成人| 女人久久www免费人成看片| 精品第一国产精品| 亚洲,欧美精品.| 欧美国产精品一级二级三级| 午夜福利乱码中文字幕| 亚洲内射少妇av| 精品人妻一区二区三区麻豆| 色婷婷av一区二区三区视频| 国产深夜福利视频在线观看| 亚洲av成人精品一二三区| 国产69精品久久久久777片| 久久免费观看电影| 国产 精品1| 狠狠婷婷综合久久久久久88av| 久久 成人 亚洲| 亚洲成av片中文字幕在线观看 | 精品少妇内射三级| 99re6热这里在线精品视频| 97在线视频观看| 男人添女人高潮全过程视频| 高清av免费在线| 国产精品久久久av美女十八| 国产精品国产三级国产av玫瑰| 久久精品国产亚洲av天美| 婷婷色综合大香蕉| 亚洲国产欧美日韩在线播放| 日韩中字成人| 内地一区二区视频在线| 草草在线视频免费看| 蜜桃国产av成人99| 观看av在线不卡| 国产无遮挡羞羞视频在线观看| 男女啪啪激烈高潮av片| 一级黄片播放器| 亚洲精品aⅴ在线观看| 亚洲精品国产av成人精品| 免费高清在线观看日韩| 精品久久久精品久久久| 成人免费观看视频高清| 亚洲欧洲日产国产| 男人爽女人下面视频在线观看| av天堂久久9| 国产乱来视频区| 尾随美女入室| 免费在线观看黄色视频的| 又黄又粗又硬又大视频| 国产亚洲精品第一综合不卡 | 午夜久久久在线观看| 色婷婷久久久亚洲欧美| 久久久精品免费免费高清| 亚洲美女视频黄频| 精品一区二区三区视频在线| 亚洲国产日韩一区二区| 国精品久久久久久国模美| 高清视频免费观看一区二区| 久久热在线av| 久久久久精品人妻al黑| 国产1区2区3区精品| 捣出白浆h1v1| 美女福利国产在线| 最近中文字幕2019免费版| 午夜91福利影院| 国产成人一区二区在线| 精品久久蜜臀av无| 国产成人a∨麻豆精品| 啦啦啦视频在线资源免费观看| 午夜精品国产一区二区电影| 男女下面插进去视频免费观看 | 亚洲国产精品一区三区| 超碰97精品在线观看| 国产一区二区激情短视频 | 天天躁夜夜躁狠狠躁躁| 国产一区二区在线观看av| 亚洲图色成人| 亚洲精品乱久久久久久| 久久精品熟女亚洲av麻豆精品| 亚洲av福利一区| 精品一区二区三卡| 欧美少妇被猛烈插入视频| 精品一品国产午夜福利视频| 欧美精品av麻豆av| 99国产精品免费福利视频| 久久午夜综合久久蜜桃| 精品少妇黑人巨大在线播放| 美女福利国产在线| 高清视频免费观看一区二区| 国产国语露脸激情在线看| 亚洲欧洲国产日韩| 嫩草影院入口| 久久人人97超碰香蕉20202| 国产福利在线免费观看视频| 国产综合精华液| 91成人精品电影| 欧美日韩成人在线一区二区| 亚洲,一卡二卡三卡| 国产深夜福利视频在线观看| 满18在线观看网站| 黄色毛片三级朝国网站| 亚洲国产看品久久| 国产精品秋霞免费鲁丝片| 最近手机中文字幕大全| 肉色欧美久久久久久久蜜桃| 9热在线视频观看99| 日韩人妻精品一区2区三区| 国产精品一区www在线观看| 一级片免费观看大全| 久久精品夜色国产| 亚洲熟女精品中文字幕| 亚洲av在线观看美女高潮| 中国三级夫妇交换| 久久精品熟女亚洲av麻豆精品| 欧美另类一区| 亚洲,欧美,日韩| 亚洲欧美日韩卡通动漫| 久久国产精品大桥未久av| 高清欧美精品videossex| 欧美xxⅹ黑人| 97精品久久久久久久久久精品| 欧美日韩成人在线一区二区| 中文字幕制服av| 亚洲国产精品专区欧美| 美女中出高潮动态图| 午夜久久久在线观看| 亚洲欧美成人精品一区二区| 亚洲精品av麻豆狂野| 久久国产精品大桥未久av| 国产亚洲精品久久久com| 大话2 男鬼变身卡| 欧美成人午夜精品| 美女主播在线视频| 精品亚洲成a人片在线观看| 亚洲av福利一区| 赤兔流量卡办理| 丰满迷人的少妇在线观看| 久久鲁丝午夜福利片| 好男人视频免费观看在线| 免费女性裸体啪啪无遮挡网站| 国产精品久久久久久精品电影小说| 午夜日本视频在线| 中文字幕制服av| 国产成人免费观看mmmm| 日韩欧美一区视频在线观看| 黑人欧美特级aaaaaa片| 成人亚洲精品一区在线观看| www日本在线高清视频| 1024视频免费在线观看| 久久久久久人人人人人| 99热网站在线观看| 久久精品aⅴ一区二区三区四区 | 成人影院久久| 女人久久www免费人成看片| 亚洲综合精品二区| 亚洲精品视频女| 美女视频免费永久观看网站| 国产成人免费观看mmmm| 日产精品乱码卡一卡2卡三| 日韩一区二区三区影片| 视频区图区小说| 在线观看人妻少妇| 国产片内射在线| 日本黄色日本黄色录像| 亚洲欧美清纯卡通| 这个男人来自地球电影免费观看 | 午夜激情久久久久久久| 亚洲国产欧美在线一区| 三级国产精品片| 成年人免费黄色播放视频| 高清av免费在线| 亚洲国产欧美日韩在线播放| 国产精品 国内视频| 国产亚洲精品久久久com| 少妇高潮的动态图| 国产精品一区二区在线不卡| 青春草视频在线免费观看| 美女国产视频在线观看| 青春草国产在线视频| 高清黄色对白视频在线免费看| 亚洲高清免费不卡视频| 亚洲第一区二区三区不卡| 久久久久精品性色| 成人手机av| 99久久人妻综合| 大片电影免费在线观看免费| 国产毛片在线视频| 国语对白做爰xxxⅹ性视频网站| 成人毛片a级毛片在线播放| 国产精品国产三级国产专区5o| 又大又黄又爽视频免费| 国产淫语在线视频| av在线播放精品| 精品熟女少妇av免费看| 韩国av在线不卡| 亚洲欧美成人综合另类久久久| 成人综合一区亚洲| xxxhd国产人妻xxx| 天堂中文最新版在线下载| 麻豆乱淫一区二区| 只有这里有精品99| 国产乱来视频区| 亚洲欧美成人精品一区二区| 一级,二级,三级黄色视频| 午夜福利影视在线免费观看| 久久女婷五月综合色啪小说| 一二三四中文在线观看免费高清| a级片在线免费高清观看视频| 97在线视频观看| 青春草视频在线免费观看| 最近最新中文字幕免费大全7| 十八禁高潮呻吟视频| 色婷婷av一区二区三区视频| 亚洲伊人色综图| 人人妻人人澡人人看| 国产爽快片一区二区三区| 久久精品国产自在天天线| 纯流量卡能插随身wifi吗| 中文精品一卡2卡3卡4更新| 欧美老熟妇乱子伦牲交| 丰满少妇做爰视频| 国产精品国产三级专区第一集| 久久精品国产a三级三级三级| √禁漫天堂资源中文www| 国产一区二区激情短视频 | 午夜精品国产一区二区电影| 2018国产大陆天天弄谢| 九草在线视频观看| 亚洲国产精品成人久久小说| 免费观看a级毛片全部| 天天影视国产精品| 91在线精品国自产拍蜜月| 久久毛片免费看一区二区三区| 最新的欧美精品一区二区| 精品久久蜜臀av无| 男人添女人高潮全过程视频| 久久人人爽人人爽人人片va| 国产成人精品一,二区| 国产男人的电影天堂91| 毛片一级片免费看久久久久| 成人综合一区亚洲| 精品午夜福利在线看| 韩国精品一区二区三区 | 91精品国产国语对白视频| 十八禁网站网址无遮挡| 在线观看美女被高潮喷水网站| 国产亚洲精品第一综合不卡 | 日韩大片免费观看网站| 99久久人妻综合| 极品人妻少妇av视频| 亚洲av在线观看美女高潮| 亚洲国产av影院在线观看| 18禁观看日本| 亚洲国产av影院在线观看| 新久久久久国产一级毛片| 国产一区二区在线观看av| 99久久人妻综合| 国产成人91sexporn| 又黄又爽又刺激的免费视频.| 又粗又硬又长又爽又黄的视频| 日韩一区二区视频免费看| 多毛熟女@视频| 久久av网站| 亚洲精品av麻豆狂野| videosex国产| 9热在线视频观看99| 少妇被粗大猛烈的视频| 成年动漫av网址| 中文字幕av电影在线播放| 22中文网久久字幕| 熟妇人妻不卡中文字幕| 亚洲欧洲精品一区二区精品久久久 | 免费大片黄手机在线观看| 国产精品久久久久久av不卡| 国产成人精品在线电影| 又黄又粗又硬又大视频| 欧美日韩综合久久久久久| 美国免费a级毛片| 寂寞人妻少妇视频99o| 黄色视频在线播放观看不卡| 大香蕉久久成人网| 午夜影院在线不卡| 一二三四中文在线观看免费高清| 1024视频免费在线观看| 久久久国产欧美日韩av| av在线播放精品| 男女边摸边吃奶| 乱码一卡2卡4卡精品| 亚洲精品久久午夜乱码| 伊人亚洲综合成人网| 黑人巨大精品欧美一区二区蜜桃 | 女人久久www免费人成看片| 少妇 在线观看| 一本—道久久a久久精品蜜桃钙片| 精品人妻熟女毛片av久久网站| 欧美日韩av久久| 亚洲综合色惰| 欧美+日韩+精品| 下体分泌物呈黄色| 美女国产视频在线观看| 精品人妻熟女毛片av久久网站| 国产成人aa在线观看| 99久久中文字幕三级久久日本| 欧美日本中文国产一区发布| 狂野欧美激情性bbbbbb| 亚洲av男天堂| 日韩,欧美,国产一区二区三区| 欧美xxxx性猛交bbbb| 亚洲精品乱久久久久久| av有码第一页| 天堂俺去俺来也www色官网| 国产国拍精品亚洲av在线观看| 亚洲人成77777在线视频| 超碰97精品在线观看| 男女边吃奶边做爰视频| 你懂的网址亚洲精品在线观看| 国产欧美日韩一区二区三区在线| 国产精品一区二区在线不卡| 寂寞人妻少妇视频99o| av线在线观看网站| 日韩成人av中文字幕在线观看| a级毛片在线看网站| 国产精品麻豆人妻色哟哟久久| 久久韩国三级中文字幕| 两性夫妻黄色片 | 亚洲国产欧美在线一区|