• <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
    陳平過河
    與敵人“坦誠相見”
    慧眼識陳平
    久久精品夜色国产| 欧美3d第一页| 国产aⅴ精品一区二区三区波| 97超碰精品成人国产| 夜夜看夜夜爽夜夜摸| 成人国产麻豆网| 久久久久久大精品| 国产精品久久视频播放| 尤物成人国产欧美一区二区三区| 国产黄片美女视频| 久久这里只有精品中国| 中国美女看黄片| 免费看av在线观看网站| 不卡视频在线观看欧美| 日本色播在线视频| 成人漫画全彩无遮挡| 亚洲成人久久爱视频| 久久人妻av系列| 免费av不卡在线播放| 五月伊人婷婷丁香| 性插视频无遮挡在线免费观看| 99热6这里只有精品| 99热精品在线国产| 此物有八面人人有两片| 亚洲最大成人av| 男插女下体视频免费在线播放| 国产久久久一区二区三区| 黄色配什么色好看| 久久精品国产亚洲av香蕉五月| 免费人成视频x8x8入口观看| 久久久久精品国产欧美久久久| 99久久九九国产精品国产免费| 精品人妻视频免费看| 亚洲性夜色夜夜综合| 51国产日韩欧美| 成人一区二区视频在线观看| 久久精品人妻少妇| a级毛片a级免费在线| 亚洲成人久久爱视频| 天堂av国产一区二区熟女人妻| 成人av一区二区三区在线看| 亚洲熟妇中文字幕五十中出| 一本一本综合久久| 看片在线看免费视频| 天天躁夜夜躁狠狠久久av| АⅤ资源中文在线天堂| 秋霞在线观看毛片| 久久精品国产亚洲av涩爱 | 国产黄a三级三级三级人| 给我免费播放毛片高清在线观看| 国产色爽女视频免费观看| 日本成人三级电影网站| 哪里可以看免费的av片| 乱系列少妇在线播放| 联通29元200g的流量卡| 97碰自拍视频| 最近中文字幕高清免费大全6| 免费看av在线观看网站| 黄色配什么色好看| 久久精品夜夜夜夜夜久久蜜豆| 一个人免费在线观看电影| 久久久国产成人精品二区| 搡老熟女国产l中国老女人| 午夜福利在线在线| 狠狠狠狠99中文字幕| 特大巨黑吊av在线直播| 日本成人三级电影网站| 精品久久国产蜜桃| 99热全是精品| 国产伦精品一区二区三区四那| 麻豆久久精品国产亚洲av| 波多野结衣高清作品| 日本一二三区视频观看| 亚洲av中文字字幕乱码综合| 国产真实伦视频高清在线观看| 免费观看人在逋| 免费看av在线观看网站| 小说图片视频综合网站| 男女边吃奶边做爰视频| 欧美中文日本在线观看视频| 亚洲专区国产一区二区| 国产av不卡久久| 亚洲欧美中文字幕日韩二区| 波多野结衣巨乳人妻| 日韩精品青青久久久久久| 激情 狠狠 欧美| 欧美最新免费一区二区三区| 99久久精品热视频| 国产精品一区www在线观看| 一级毛片aaaaaa免费看小| 91午夜精品亚洲一区二区三区| 成人特级av手机在线观看| 午夜激情福利司机影院| 成人永久免费在线观看视频| www.色视频.com| 亚洲va在线va天堂va国产| 久久久久久久亚洲中文字幕| 神马国产精品三级电影在线观看| 免费搜索国产男女视频| 亚洲欧美日韩高清在线视频| 2021天堂中文幕一二区在线观| 免费大片18禁| 观看免费一级毛片| 国产高潮美女av| 国产一区二区激情短视频| 精品无人区乱码1区二区| 亚洲性久久影院| av国产免费在线观看| 一个人观看的视频www高清免费观看| 别揉我奶头 嗯啊视频| 全区人妻精品视频| 日韩强制内射视频| 精品久久久久久成人av| 三级国产精品欧美在线观看| 午夜亚洲福利在线播放| 亚洲不卡免费看| 成年女人看的毛片在线观看| 偷拍熟女少妇极品色| 最新在线观看一区二区三区| 欧美区成人在线视频| 色噜噜av男人的天堂激情| 女生性感内裤真人,穿戴方法视频| 国产真实乱freesex| 色av中文字幕| 白带黄色成豆腐渣| 国国产精品蜜臀av免费| 99热这里只有精品一区| 日韩欧美 国产精品| 欧美一区二区亚洲| 久久久久精品国产欧美久久久| 国内精品美女久久久久久| 一级毛片电影观看 | 露出奶头的视频| 亚洲一区高清亚洲精品| 白带黄色成豆腐渣| 成人特级av手机在线观看| 男人狂女人下面高潮的视频| 国产午夜福利久久久久久| 国产一区二区激情短视频| 麻豆乱淫一区二区| 亚洲专区国产一区二区| 综合色av麻豆| 国产又黄又爽又无遮挡在线| av天堂在线播放| 欧美zozozo另类| 男女那种视频在线观看| 村上凉子中文字幕在线| 亚洲四区av| 欧美高清成人免费视频www| 美女大奶头视频| 六月丁香七月| 日韩在线高清观看一区二区三区| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲国产欧美人成| 国产精品伦人一区二区| 中文资源天堂在线| 一级黄片播放器| 欧美日韩国产亚洲二区| 我的女老师完整版在线观看| 热99re8久久精品国产| 欧美绝顶高潮抽搐喷水| 亚洲国产欧洲综合997久久,| 三级毛片av免费| 国产精品乱码一区二三区的特点| 亚洲av免费在线观看| 在线免费观看不下载黄p国产| 久久鲁丝午夜福利片| 亚洲av.av天堂| 久99久视频精品免费| 两性午夜刺激爽爽歪歪视频在线观看| 97超视频在线观看视频| 热99re8久久精品国产| 欧美激情久久久久久爽电影| 精品久久久久久久久av| 国产精品久久电影中文字幕| 日韩av不卡免费在线播放| 一区二区三区四区激情视频 | 内射极品少妇av片p| 国产精品av视频在线免费观看| 男人和女人高潮做爰伦理| 久久久久免费精品人妻一区二区| 高清日韩中文字幕在线| 午夜激情福利司机影院| 久久精品91蜜桃| 我要看日韩黄色一级片| 国产乱人偷精品视频| 久久精品人妻少妇| 亚洲五月天丁香| 欧美最黄视频在线播放免费| 久久这里只有精品中国| 国产精品美女特级片免费视频播放器| 一级黄色大片毛片| 国产精品一及| 国产真实乱freesex| 熟妇人妻久久中文字幕3abv| 99在线人妻在线中文字幕| 成年女人毛片免费观看观看9| 熟女电影av网| 亚洲av第一区精品v没综合| 最近最新中文字幕大全电影3| 此物有八面人人有两片| 99久国产av精品国产电影| 欧美日本视频| 日韩中字成人| a级毛色黄片| 国产熟女欧美一区二区| 舔av片在线| 在线播放国产精品三级| 欧美日韩在线观看h| 久久九九热精品免费| 亚洲人成网站在线观看播放| 久久婷婷人人爽人人干人人爱| 国产人妻一区二区三区在| 欧美激情久久久久久爽电影| 成人漫画全彩无遮挡| 国产 一区 欧美 日韩| 日产精品乱码卡一卡2卡三| 日韩欧美精品免费久久| 人妻制服诱惑在线中文字幕| 久久久久久大精品| 精品久久久噜噜| av国产免费在线观看| 亚洲国产色片| 欧美日韩一区二区视频在线观看视频在线 | 亚洲中文字幕一区二区三区有码在线看| 淫秽高清视频在线观看| 国内精品久久久久精免费| 日本色播在线视频| 麻豆乱淫一区二区| 日韩欧美 国产精品| 国产精品久久久久久av不卡| 伊人久久精品亚洲午夜| 综合色av麻豆| 美女被艹到高潮喷水动态| 草草在线视频免费看| 久久精品国产亚洲网站| 亚洲成人久久性| 亚洲aⅴ乱码一区二区在线播放| 美女被艹到高潮喷水动态| 九九久久精品国产亚洲av麻豆| 在线免费观看不下载黄p国产| 99视频精品全部免费 在线| 亚洲性夜色夜夜综合| 99热这里只有是精品50| 九九久久精品国产亚洲av麻豆| 欧美色欧美亚洲另类二区| 国产一级毛片七仙女欲春2| 一夜夜www| a级一级毛片免费在线观看| 日韩欧美一区二区三区在线观看| 欧美极品一区二区三区四区| 少妇人妻一区二区三区视频| 别揉我奶头~嗯~啊~动态视频| 免费无遮挡裸体视频| 日韩强制内射视频| 日日摸夜夜添夜夜添小说| 久久精品国产亚洲av香蕉五月| 99热这里只有是精品在线观看| 欧美三级亚洲精品| 欧美zozozo另类| 亚洲欧美成人精品一区二区| 女生性感内裤真人,穿戴方法视频| 级片在线观看| 久久99热这里只有精品18| a级毛色黄片| 悠悠久久av| 精品久久久久久久末码| 成年版毛片免费区| 69人妻影院| 天堂av国产一区二区熟女人妻| 高清日韩中文字幕在线| 色尼玛亚洲综合影院| 亚洲欧美日韩无卡精品| 在线看三级毛片| 亚洲中文字幕一区二区三区有码在线看| 亚洲国产精品成人综合色| 久久亚洲国产成人精品v| 哪里可以看免费的av片| 国产aⅴ精品一区二区三区波| 黄片wwwwww| 99热全是精品| 亚洲精品日韩在线中文字幕 | 精品国内亚洲2022精品成人| 精品乱码久久久久久99久播| 亚洲国产精品国产精品| 午夜a级毛片| 国产精品一二三区在线看| 寂寞人妻少妇视频99o| 亚洲婷婷狠狠爱综合网| 久久鲁丝午夜福利片| 中文字幕熟女人妻在线| 精品久久久久久久人妻蜜臀av| 黑人高潮一二区| 亚洲精品乱码久久久v下载方式| 亚洲自拍偷在线| 国产精品一区www在线观看| 欧美日韩综合久久久久久| 欧美激情在线99| 在线播放无遮挡| 国产一区二区在线av高清观看| 乱码一卡2卡4卡精品| 精品久久久久久久久亚洲| 国产精品一二三区在线看| 国产爱豆传媒在线观看| 美女cb高潮喷水在线观看| 白带黄色成豆腐渣| 色播亚洲综合网| 国产伦在线观看视频一区| 亚洲18禁久久av| 欧洲精品卡2卡3卡4卡5卡区| 欧美bdsm另类| 欧美成人免费av一区二区三区| 久久人人精品亚洲av| 国产中年淑女户外野战色| 三级国产精品欧美在线观看| 精品一区二区三区视频在线| 国产精品一区二区性色av| 亚洲欧美清纯卡通| 在线a可以看的网站| 亚洲国产精品成人综合色| 精品少妇黑人巨大在线播放 | 亚洲国产欧美人成| 看片在线看免费视频| 久久精品影院6| 乱码一卡2卡4卡精品| 久久婷婷人人爽人人干人人爱| 精品日产1卡2卡| 午夜影院日韩av| 午夜激情欧美在线| 97碰自拍视频| 欧美色欧美亚洲另类二区| 男女啪啪激烈高潮av片| 1024手机看黄色片| 国产精品一区二区性色av| 99久久精品国产国产毛片| 国产精品亚洲一级av第二区| 一本久久中文字幕| 国产精品亚洲一级av第二区| 春色校园在线视频观看| av.在线天堂| 欧美成人一区二区免费高清观看| 女同久久另类99精品国产91| 俺也久久电影网| 国产精品亚洲一级av第二区| 禁无遮挡网站| 国产极品精品免费视频能看的| 18禁裸乳无遮挡免费网站照片| 亚洲av五月六月丁香网| 美女内射精品一级片tv| av福利片在线观看| 春色校园在线视频观看| 人妻夜夜爽99麻豆av| 禁无遮挡网站| 免费在线观看影片大全网站| 麻豆成人午夜福利视频| 日韩av不卡免费在线播放| 日韩 亚洲 欧美在线| 国产精品综合久久久久久久免费| 午夜激情欧美在线| 国产精品福利在线免费观看| 女人被狂操c到高潮| 亚洲成人精品中文字幕电影| 亚洲av不卡在线观看| 久久国产乱子免费精品| 精品久久久久久成人av| a级毛片a级免费在线| 日日摸夜夜添夜夜添小说| 色哟哟·www| 大又大粗又爽又黄少妇毛片口| 久久久久久久久久黄片| 亚洲三级黄色毛片| 97在线视频观看| 少妇的逼水好多| 国产精品爽爽va在线观看网站| 午夜视频国产福利| 男女边吃奶边做爰视频| 日韩成人伦理影院| 日本a在线网址| 日韩精品中文字幕看吧| 97热精品久久久久久| 日本成人三级电影网站| 午夜精品在线福利| 一进一出抽搐gif免费好疼| 1024手机看黄色片| 日韩精品青青久久久久久| 秋霞在线观看毛片| 少妇人妻一区二区三区视频| 欧美高清性xxxxhd video| 婷婷六月久久综合丁香| 十八禁国产超污无遮挡网站| 床上黄色一级片| 欧美潮喷喷水| 露出奶头的视频| 久久亚洲国产成人精品v| 国产精品女同一区二区软件| 少妇熟女aⅴ在线视频| 91在线精品国自产拍蜜月| 成人高潮视频无遮挡免费网站| 日日干狠狠操夜夜爽| 欧美不卡视频在线免费观看| 久久久久国内视频| 亚洲欧美精品自产自拍| 一进一出抽搐动态| 亚洲欧美日韩无卡精品| 别揉我奶头 嗯啊视频| 91麻豆精品激情在线观看国产| 国产aⅴ精品一区二区三区波| 亚洲精品色激情综合| 青春草视频在线免费观看| 成年女人永久免费观看视频| 伦精品一区二区三区| 精品免费久久久久久久清纯| 变态另类成人亚洲欧美熟女| 老司机影院成人| 一区二区三区高清视频在线| 国产日本99.免费观看| 深夜a级毛片| 国产91av在线免费观看| 高清日韩中文字幕在线| av专区在线播放| 舔av片在线| 天天一区二区日本电影三级| 久久久久国内视频| 国产高清视频在线播放一区| 色视频www国产| 黄片wwwwww| 99久久成人亚洲精品观看| 国产成年人精品一区二区| 国产欧美日韩精品一区二区| 激情 狠狠 欧美| 伊人久久精品亚洲午夜| 嫩草影院入口| 亚洲欧美精品自产自拍| 久久精品国产亚洲网站| 一区福利在线观看| 丝袜喷水一区| 亚洲国产欧美人成| 日本熟妇午夜| 久久久久九九精品影院| 男女啪啪激烈高潮av片| 晚上一个人看的免费电影| 日韩欧美一区二区三区在线观看| 联通29元200g的流量卡| 超碰av人人做人人爽久久| 成人高潮视频无遮挡免费网站| 久久天躁狠狠躁夜夜2o2o| 亚洲欧美精品自产自拍| 日本撒尿小便嘘嘘汇集6| 亚洲国产精品sss在线观看| 久久久国产成人精品二区| 久久这里只有精品中国| 亚洲在线观看片| 久久精品综合一区二区三区| 亚洲最大成人中文| 国产精华一区二区三区| 久久热精品热| 欧美中文日本在线观看视频| 最近手机中文字幕大全| 亚洲成人久久性| 麻豆精品久久久久久蜜桃| 少妇人妻一区二区三区视频| 午夜久久久久精精品| 国产黄a三级三级三级人| 91狼人影院| 男女视频在线观看网站免费| 两个人的视频大全免费| 男人狂女人下面高潮的视频| 国产私拍福利视频在线观看| 看黄色毛片网站| 日本爱情动作片www.在线观看 | 最近2019中文字幕mv第一页| 国产 一区精品| 又爽又黄a免费视频| 精品久久久久久久久av| 国产探花极品一区二区| 美女cb高潮喷水在线观看| 亚洲精品亚洲一区二区| 国产精品久久电影中文字幕| 久久精品影院6| 国产精品亚洲美女久久久| 99久久无色码亚洲精品果冻| 色哟哟哟哟哟哟| 国产视频内射| 99久久精品热视频| 亚洲av中文字字幕乱码综合| 丝袜美腿在线中文| 亚洲成人久久性| 又爽又黄无遮挡网站| 99热网站在线观看| 欧美潮喷喷水| 亚洲av一区综合| 老熟妇乱子伦视频在线观看| 亚洲国产色片| 黄片wwwwww| 直男gayav资源| 国产成人福利小说| 天天一区二区日本电影三级| 国内精品宾馆在线| 欧美日本亚洲视频在线播放| 中文在线观看免费www的网站| 亚洲精品国产成人久久av| 老熟妇乱子伦视频在线观看| 亚洲欧美中文字幕日韩二区| www日本黄色视频网| 丝袜喷水一区| 一级毛片久久久久久久久女| 成人特级黄色片久久久久久久| 久久久久性生活片| 日韩制服骚丝袜av| 联通29元200g的流量卡| 男女之事视频高清在线观看| 成人特级黄色片久久久久久久| 日本撒尿小便嘘嘘汇集6| 久久99热这里只有精品18| 天堂av国产一区二区熟女人妻| 丝袜美腿在线中文| 日本免费一区二区三区高清不卡| 老司机福利观看| 国产欧美日韩精品亚洲av| 两性午夜刺激爽爽歪歪视频在线观看| 一本久久中文字幕| 蜜臀久久99精品久久宅男| 99久久成人亚洲精品观看| a级毛片免费高清观看在线播放| 黄色配什么色好看| 久久午夜亚洲精品久久| 国产精品无大码| 国产在线精品亚洲第一网站| 国产精品女同一区二区软件| 亚洲第一区二区三区不卡| 国产日本99.免费观看| 香蕉av资源在线| 看黄色毛片网站| 一个人观看的视频www高清免费观看| 91久久精品电影网| 成熟少妇高潮喷水视频| 天堂√8在线中文| 久久精品国产亚洲av涩爱 | 精品人妻熟女av久视频| 国产精品福利在线免费观看| 国产在线精品亚洲第一网站| 欧美一区二区亚洲| 大香蕉久久网| 国产精品久久久久久av不卡| 成人av一区二区三区在线看| 免费黄网站久久成人精品| 国产一区二区在线av高清观看| 亚洲色图av天堂| 日韩 亚洲 欧美在线| 久久精品国产亚洲av香蕉五月| 午夜日韩欧美国产| 中国美女看黄片| 亚洲真实伦在线观看| 国产在视频线在精品| 国产一区二区三区在线臀色熟女| 国产熟女欧美一区二区| 欧美中文日本在线观看视频| 天堂√8在线中文| 白带黄色成豆腐渣| 99热6这里只有精品| 99久久中文字幕三级久久日本| 国产精品亚洲美女久久久| 在现免费观看毛片| 久久精品久久久久久噜噜老黄 | 天堂网av新在线| 日本爱情动作片www.在线观看 | 蜜臀久久99精品久久宅男| 老师上课跳d突然被开到最大视频| 日韩欧美免费精品| 亚洲欧美精品自产自拍| 校园春色视频在线观看| 综合色av麻豆| av.在线天堂| 精品午夜福利视频在线观看一区| АⅤ资源中文在线天堂| h日本视频在线播放| 少妇人妻精品综合一区二区 | 亚洲精品国产成人久久av| 国产精品亚洲一级av第二区| 久久精品国产99精品国产亚洲性色| 十八禁网站免费在线| 欧美不卡视频在线免费观看| 99热全是精品| 亚洲18禁久久av| 白带黄色成豆腐渣| 麻豆精品久久久久久蜜桃| 亚洲18禁久久av| 亚洲色图av天堂| 久久99热6这里只有精品| 十八禁网站免费在线| 国产白丝娇喘喷水9色精品| 美女 人体艺术 gogo| 97超视频在线观看视频| 又黄又爽又免费观看的视频| 亚洲人与动物交配视频| 丰满人妻一区二区三区视频av| 日韩,欧美,国产一区二区三区 | 在现免费观看毛片| 成人国产麻豆网| 日日摸夜夜添夜夜添av毛片| 最近在线观看免费完整版| 久久亚洲精品不卡| av在线亚洲专区| 亚洲国产精品成人久久小说 | 免费观看的影片在线观看| 97人妻精品一区二区三区麻豆| 国产伦精品一区二区三区四那| 俺也久久电影网| 欧美bdsm另类| 中文资源天堂在线| 在线免费十八禁| 高清日韩中文字幕在线| 久久久久国产网址| 中文字幕精品亚洲无线码一区| 国产精品人妻久久久影院|