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

    A plastic strain energy method exploration between machined surface integrity evolution and torsion fatigue behaviour of low alloy steel

    2022-11-13 07:30:12YongWANGXiinWANGZhiingLIUShuyoLIUSitoWANGHongtoCHENCiSONGYifnBAIPiWANGYngLIU
    CHINESE JOURNAL OF AERONAUTICS 2022年10期

    Yong WANG, Xiin WANG, Zhiing LIU,*, Shuyo LIU, Sito WANG,Hongto CHEN, Ci SONG, Yifn BAI, Pi WANG, Yng LIU

    a School of Mechanical Engineering, Beijing Institute of Technology, Beijing 100081, China

    b Key Laboratory of Fundamental Science for Advanced Machining, Beijing Institute of Technology, Beijing 100081, China

    c Beijing North Vehicle Group Corporation, Beijing 100072, China

    KEYWORDS Cyclic hysteresis loop;Multistage machining processes;Strain energy;Surface integrity evolution;Torsional fatigue

    Abstract To explore the evolution mechanism of multistage machining processes and torsional fatigue behaviour based on strain energy for the first time and provide process optimization of axis parts of low alloy steel for service performance,four multistage machining processes were applied to the 45CrNiMoVA steel,including the Rough Turning process(RT),RT+the Finish Turning process(FRT),FRT+the Grinding process(GFRT)and RT+the Finish Turning process on dry cutting condition (FRT0). The result showed that the FRT process’s average low-cycle torsional fatigue life increased by 50% when it evolved from the RT process. The lower surface roughness of Ra 1.3 μm caused the total strain energy to increase by 163.8 Pa mm/mm instead of the unchanged strain energy density, and the crack feature evolved from some specific bulges to flat shear plane characteristics. When the GFRT process evolved from the FRT process, its average fatigue life increased by 1.45 times, compared with the RT process. Plastic strain amplitude decreased by 21%, and the strain energy density decreased by 4% due to more considerable compressive residual stress(-249 MPa).Plastic deformation layer depth had a consistent tendency with surface roughness. In this paper, surface integrity evolutions on cyclic characteristics and fatigue behaviour have also been explained. A fatigue life prediction model based on the energy method for machined surface integrity is proposed.

    1. Introduction

    As low alloy steel, 45CrNiMoVA steel is widely used in highspeed rail, automobile and other torsion bars due to high specific strength, high hardenability,etc1,35.Numerous factors affect the fatigue life of torsion bars in a natural service environment, such as material and metallurgy, part geometry,machining process, surface condition, etc. After machining,the free surface of the torsion bar bears a large load, which is affected by the external environment,such as thermal stress,oxidation, wear, etc. Especially in the actual structure and components, the crack initiates from the concentrated stress region due to cyclic fatigue load. The crack initiation and propagation are attributed to the surface integrity of the machining process2.Thus,it is vital to study the effects of surface integrity of the machining process on the torsional fatigue of shaft materials. At present, the research of the low cyclic fatigue life of low alloy steel by many scholars mainly focuses on the following aspects: (A) cyclic stress-strain curves3; (B)cyclic hardening-softening properties4;(C)relationship of cyclic responses with dislocation structure evolution5; (D) cracks initiation and propagation behaviour6.Although several studies have been implemented on cyclic stress-strain response,most of them focus on the cyclic deformation behaviour and dislocation structure features of average stress or heat treatment under uniaxial tensile load.However,due to the complex shape of the torsion bar or the harsh environment, the stress born by the torsion bar is often more complex than the simple uniaxial tension-compression fatigue load, and the unreasonable machined surface layer will make the parts failure occur prematurely, which brings significant risk to the application of parts. At the same time, few studies focus on the evolution of machined surface layers induced in multistage machined processes on the low cycle fatigue properties of low alloy steel.Therefore,a comprehensive exploration of the cyclic hysteresis loop behaviour and fatigue mechanism of torsion bar induced by the multi-functional process of machining will be beneficial for preventing structural failure and life loss, and providing a theoretical basis for the optimization design process of the machining process from the aspect of strain energy method.

    This paper aims at surface integrity evolution and fatigue evaluation after rough turning, finish turning and grinding mode for 45CrNiMoVA steel under pure torsional fatigue loadings.Much effort is taken to correlate the cyclic hysteresis loop behaviour and stress-strain response with the surface integrity evolution in every link of machined processes at different cyclic stages. Moreover, it is attempted to associate the surface integrity evolution with the fatigue life prediction model based on the energy method.

    2.Preparation of specimens with multistage machining processes

    2.1. Material

    The common shaft parts material(45CrNiMoVA steel)is prepared in this test,and chemical compositions are shown as follows: 0.8wt%-1.1 wt%Cr, 1.3wt%-1.8 wt%Ni,0.42wt%-0.49 wt%C, 0.5wt%-0.8 wt%Mn, 0.17wt%-0.37 wt%Si, 0.2wt%-0.3 wt%Mo,0.1wt%-0.2 wt%V.Fig.1 illustrates grain orientation and size distributions of 45CrNiMoVA steel with the grain sizes of 25 μm by EBSD, and Fig. 1(d) show the final machined shape of torsional specimens for the actual service conditions of materials based on the ISO 1352:2011.The material is annealed and mechanical properties (ASTM A370) are shown in Table 1.

    2.2. Multistage machining processes

    Table 1 Mechanical properties of 45CrNiMoVA steel.

    Four groups of multistage processes were applied to assess the effect of surface integrity induced in every link of the machining process on low-cycle torsional fatigue of 45CrNiMoVA steel. i.e., the rough turning process (RT), RT+ the finish turning process (FRT), FRT+ the grinding process (GFRT)and RT+ the finish turning process on dry cutting condition(FRT0). Multistage machining processes arrangements were carried out in Fig. 2. The turning processes were prepared in a CNC lathing-milling centre (Mazak nexus 200-II ML,Japan), and the turning tools VNMG160408 VP15TF (Mitsubishi, Japan) with a tool corner radius of 0.8 mm were prepared for the RT process.On the other hand,the turning tools VNMG160404-MS VP05RT (Mitsubishi, Japan) with a tool corner radius of 0.4 mm were prepared for the FRT and FRT0 processes.

    Multistage machining processes parameters are shown in Table 2. The RT group were operated for the first group,and the machining process parameters under wet condition are n=1200 r/min;ap=0.5 mm;f=0.35 mm/r,respectively.It should be noted that the final size of the four groups machining process is 12.5 mm, thereby reducing the influence of specimen size on the fatigue performance. As shown in Fig. 2, FRT specimens’ final machining process parameters are n = 2037 r/min, ap= 0.25 mm, f = 0.12 mm/r, respectively,and the machining allowance(ap)of 0.25 mm is retained for passing RT parameters. For GRFT specimens, the final speed of the grinding wheel nf= 25 m/s, ns= 150 r/min,ap= 0.01 mm, f = 200 mm/min, respectively, and the finish turning allowance(ap1+ap2)passing RT process and grinding allowance (ap2) passing FRT parameters are 0.5 mm and 0.25 mm, respectively. One group of finish turning under dry cutting conditions (FRT0) is operated to study the cyclic deformation behaviour of dry cutting under pure torsional fatigue load. FRT0 has the same machining process parameters as the FRT group, and the other three groups are under the wet cutting process.

    2.3. Surface integrity measurements

    A 3D Laser Scanning Microscope (VK-X100, Keyence) was used to measure surface morphology and roughness.Consider-ing that the crack initiation passing the fatigue torsion tests is always located in the arc transition region near the clamping end, measurements were carried out every 120° in the circumferential region at the arc transition region, and the final roughness value was averaged. An X-ray stress analyzer (X-350A,Astor)was used to analyze the circumferential and axial residual stresses on the surface and the conventional X-ray diffraction method. The LXRD instrument was prepared to measure different circumferential areas at position A three times using the sin2ψ method,and the final residual stress value was averaged.

    Table 2 Multistage machining processes applied on 45CrNiMoVA steel.

    A field-emission scanning electron microscope (ERA-8900FE, Kanagawa) was used to get the grain properties of the machined surface layer through electron backscatter diffraction analysis. The OIM Data Collection 7.0 software was prepared to analyze the OIM data. Moreover, all specimens were obtained in cross-section, which was perpendicular to the axis. Considering that the surface layer induced by machining processes is easily removed by electropolishing,the specimens were obtained by mechanical polishing and vibration polishing instead of electrolytic polishing7.

    2.4. Testing condition

    Low cyclic fatigue behaviours of specimens were completed on MTS-250kN809 closed-loop servo-hydraulic testing machine controlled by a microcomputer.The total strain-controlled test was set to a strain ratio of 0, a sine wave with a frequency of 0.2 Hz.The epsilon circumferential and axial extensometer was used to measure the strain changes in each torsional cycle by controlling the shear strain at room temperature, where an extensometer of 25 mm was located on the gauge and the maximum torsion angle is 3°,i.e.the strain amplitude is 0.013 mm/mm.

    3. Experimental results

    3.1. Machined surface layer evolution

    3.1.1. Surface topography and roughness

    Surface topography evolution of 45CrNiMoVA steel in multistage machining processes is shown in Figs.3-7.A series of 3D surface topography of height distribution is shown on the right side, and optical images are on the left. In the rough turning process,considerable cutting depth and feed rate remove most machining allowances,and make the surface show great peaks and valleys(Fig.3(b))with scratches(Fig.3(a))mainly caused by the tool nose radius and feed rate. The theoretical surface roughness Rmaxfollows f2/8r, and the measured Ryof 24.2 μm is larger than the ideal Rmaxof 19.1 μm in the sampling length.

    Under the finish process parameters, the introduction of dry cutting condition causes chatter marks where high plastic flow (Fig. 4(a)) on the surface and the peaks in the troughs become irregular(Fig.4(b)).Under dry conditions,the friction in the machined process increases, and scratches (Fig. 4(a))appears. The temperature in the cutting zone increases, which results in the increment of high plastic flow and the formation of chatter marks. At the same time, mechanical adhesion(Fig. 4(a)) can be seen, and the measured value Ryis 43.2 μm.

    When the finish turning process under wet conditions was introduced,the low feed rate improves the geometric accuracy and smoothen the rough surface.Comparing Fig.5(b)to Fig.3(b), small tool nose radius and slow feed rate make the width significantly narrower even if doubling the image the distance of peaks and valleys.Both wet and dry cutting conditions feature high plastic flow, as shown in Fig. 4(a), Fig. 5(a). The ideal surface roughness Rmaxis 4.5 μm,and the measured value Ryis 8.9 μm.

    The grinding process is taken as the last process, and the grinding surface is composed of lots of numerous minor scratches between grits and the workpiece.Long and thin grits scratches show minimal spacing(Fig.6(b)).After grinding,the surface morphology is relatively smooth,thus reducing the surface roughness Ra. However, Ry(23.6 μm) does not decrease.Although the small grits spacing removes the great peaks and valleys in the finish process,it also makes the surface produce narrower and sharper peaks and valleys(Fig.6(a)).Consequently,most of the surface roughness parameters do not get smoother(Fig.6(c)),even if they are inferior to that in the finish process under wet cutting conditions, except for Rsm.Roughness values of Ra, Ry, Rzand Rsmin the multistage machining processes are shown in Fig. 7.

    3.1.2. Residual stress change

    Compressive residual stress (CRS) reduces the equivalent external load by offsetting the moderate applied load and improves fatigue crack initiation life through the reduction of the initial microcrack growth factor. Residual stress evolution in different processes is illustrated in Fig. 8. Firstly, the residual stresses on the surface of the original round bar through wire electro-discharge machining were measured prior to the rough turning process, when the residual stresses in circumferential or axial directions are 427 MPa and 325 MPa,respectively. When the rough turning process was introduced,the residual tensile stress decreases because compressive residual stress induced by large plastic deformation offsets the moderate tensile stress.When the finish turning processes under dry or wet cutting conditions were introduced,the plastic deformation increases, and the tensile residual stress in the axial direction decreases. However, the residual stresses in the circumferential direction almost remain unchanged and are within the error band ranges. Compared with the axial direction, the smaller curvature radius in the circumferential direction is more challenging to produce a considerable large plastic strain8. After the last grinding process, the residual stress transforms from tension to compression compared with the finish turning process. The axial residual stresses of-312 MPa are much higher than that(-249 MPa)in a circumferential direction, consistent with 300 M steel under wet cutting conditions reported by Dang et al.9. It indicates that the compressive residual stress between grits and the machined surface is greater than the tensile residual stress induced by the thermal effect,and the compressive residual stress appears.Balart et al.10proposed a critical transition temperature for the onset of residual stresses transforming from tension to compression on the surface of steel materials, and the axial CRS is greater than the circumferential CRS.

    After fatigue tests,each process on the surface tends to present a specific compressive residual stress (Fig. 8). This can be attributed to the fact that the asymmetric cyclic load has a specific pre-twisting effect11during the fatigue test. Plastic deformation caused compressive residual stress that occurred on the surface when shear stress was unloaded after fatigue.Simultaneously, the cyclic loading process results in plastic deformation between the grains with more considerable residual stress tending to be uniform12,13,and compressive residual stress becoming smaller.

    3.1.3. Grain size and grain boundary evolution

    The machined layer’s grain size and grain boundary have an important influence on fatigue crack initiation. The EBSD results of the RT process are shown in Fig. 9, where Fig. 9(a), Fig. 9(b) and Fig. 9(c) present the crystal orientation,image quelity (IQ) and grain boundary distributions of different machined processes, respectively. The first line in yellow(Fig. 9(a), Fig. 9(b)) shows the surface, while the area above the first line in yellow presents noise characteristic in crystal orientation and turns black after discarding. The region with an IQ smaller than 46.4 was regarded as the surface. On the layer presented between the first and second line in Fig. 9(a),due to this layer has an IQ smaller than 77.8,the crystal orientation is presented randomly. The crystal grains of this layer presented a smaller grain size due to hard working14. Thus,the image quality (IQ) on layer I is poor. It belongs to the sub-grain layer (named layer II) indicated by the red arrow,which is caused by the increased geometrically necessary(GN) dislocation density due to hard working7. Layers I and II are collectively referred to as the grain refinement layers.Many high-angle grain boundaries were presented on the two layers due to the formation of many refined grains,as indicated by the red boundary in Fig.9(c).Layer III occurs due to the plastic deformation in a grain because of the change of dislocation orientation induced by multistage machining processes. As shown with the black-oriented arrows in Fig. 9(a),there are many plastic deformation traces, and the grain boundaries mainly range from 5° to 15° (as shown with the black circles in Fig.9(c)).Layers I,II and Layer III are collectively referred to as plastic deformation layers.

    Quantifying the depths of the grain refinement layer and the plastic deformation layer is of vital importance to the study of the influencing factors of fatigue life. KAM refers to the local dislocation parameter related to geometrically necessary(GN)dislocation density15, where a new refined grain satisfies grain tolerance angle 15°, as shown in Fig. 10(a). Fig. 10(d) shows that there is average grain size 10.9 μm for RT process at the depth of 11 μm. Considering grain size (25 μm) in the matrix(Fig.10(e)),it can be regarded as refined grain before grain size stabilizes at 25 μm. While GROD represents the local misorientation parameter related to the plastic deformation16, large plastic deformation still makes the dislocation angle between points change in grain. It corresponds to higher-angle grain boundaries of disguise more than 5°in a grain. The plastic deformation depth layer was measured using the same method as that used for refined grained layer depth. Fig. 10(g) shows the distribution of GROD at the subsurface of the rough turning process.

    In the rough turning process, a large machining allowance produces not only a large grain refinement layer of 16 μm(Fig. 10(a), Fig. 10(e)) but also a deep plastic deformation layer of 39 μm(Fig.10(b),Fig.10(g))at the subsurface.When the finish turning process under wet conditions was introduced, it removed the sizeable plastic deformation layer of the RT process.In the meantime,it produced a shallow plastic deformation layer of 20.1 μm and a deep grain refined layer of 8 μm.However,the finish turning process under dry conditions introduced the deep plastic deformation layer of 46.4 μm,probably attributed to the chatter between the tool and workpiece under dry conditions. After grinding, it has the smallest grain refinement layer of 4.8 μm and plastic deformation layer of 12.2 μm (as indicated in Fig. 10(c)). The depth of the grain refinement layer and plastic deformation layer exhibit the same evolution trend.

    3.1.4. Microhardness evolution

    The microhardness evolution on the surface before and after fatigue tests is shown in Fig.11(a).The microhardness in multistage processes after fatigue tests is markedly sturdier than that before fatigue tests,which is attributed to the increment in dislocation density induced by reciprocating dislocation motion during fatigue tests.As shown in Fig.11(b),the microhardness of the machined surface layer decreases with the increasing of depth, and it is attributed to work hardening on metal surface of the specimens under the combined action of plastic deformation and cutting heat. Due to the existence of work hardening,the plastic deformation of the surface layer goes deep into a certain distance below the surface,resulting in elastic deformation in the middle layer to balance the plastic deformation in the subsurface layer, and the innermost layer is the undeformed matrix. Therefore, the microhardness of specimen has a certain distribution law along the layer depth.As shown in Fig.11(b),microhardness depth is within 100 μm.It should be noted that the greater the surface microhardness,the deeper the impacted layer. However, the finish turning under dry conditions and grinding result in lower microhardness on the surface (Fig. 11(a)). Machining chatter or sliding friction of abrasive particles caused the enhancement of surface thermal softening effect, which offsets the strengthening effect of plastic deformation on the surface.

    3.2. Effects of surface integrity evolution on cyclic stress-strain response characteristics

    3.2.1. Monotonic torsional characteristics Fig. 12 shows the monotonic mechanical properties in multistage processes, where the torsional properties are presented in Table 3.Two same monotonic mechanical experiments were carried out in FRT and GFRT processes to test the repeatability of the results.As illustrated in Fig.12,the torsional characteristic curves approach each other. The torsional coefficient and exponent fit each other best in teams from the Ramberg-Osgood type relation17, as shown:

    When the finish processes were introduced, it can be seen from Table 3 that the evolutions of all the shear property parameters of the FRT0 specimens change little compared with those of the RT specimens, and the trend curves are almost coincident. However, FRT specimens have smaller shear yield strength, which the refined grain layer may cause.

    The literature18indicated similar conclusions for yield strength on low alloy steel. The grain refinement layer on the surface makes yield strength decrease. However, the ultimate strength changes little. The grinding process shows the highest yield strength, shear modulus, strength coefficient and strain hardening exponent than the other specimens. This can be attributed to the plastic deformation of the GRFT specimens. The grinding process produces a smaller depth of plastic deformation layer, indicating that the resistance damage is minimal.More importantly, grinding generates the most considerable compressive residual stress, offsets the external load to complete the torsional fracture.

    3.2.2. Shear cyclic characteristics

    Mechanical characteristics of cyclic torsion in multistage machined processes are presented in Fig. 13. It is evident that from Fig. 13(a), the cyclic maximum shear stress responses on the turning processes are mainly characterized by similar cyclic softening rates (as indicated by the red, blue and black curves of Fig. 13(a)). However, the cyclic maximum shear stress response in the grinding process is characterized by a slower cyclic hardening rate, and the increment of τminexceeds the decrement of τmaxfrom Fig. 13(b). When the FRT0 process evolves from the RT process, τmax, τminchange little, this is probably due to the same plastic deformation depth causedby solid flutter in the FRT0 process.In the FRT process,there is a minimum shear stress τmindue to the refined grain.

    Table 3 Hysteresis loop characteristics of machined specimens.

    Grinding is the final machining process,where the considerable compressive residual stress leads to τminpresent cyclic hardening (Fig. 13(b)). Several researcher19,20showed that the cyclic hardening of the material was ascribed to microstructural barriers or plastic incompatibilities. However,τmaxshows cyclic softening.The cyclic softening was attributed to the rearrangement and changes in dislocations microstructure21,22. Saitova et al.23showed that high dislocation density in the undeformed material tends to reconfigure and decrease due to cyclic plastic deformation.

    Considering that elastic deformation is reversible, only plastic strain can result in irreversible fatigue damage. Therefore,it is more reasonable to use plastic strain energy as a damage parameter.When the coordinate origin is set at the bottom point of the hysteresis loop,as shown in Fig.14(a),it satisfies:

    The summaries of the first cycle and 0.5Nfcycle in multistage processes are listed in Table 3. By comparing the τ-γ behaviours of the grinding and turning processes at 0.5Nfcycle. It is observed that the turning processes resulted in a similar stable hysteresis loop,and the grinding has the smallest plastic strain amplitude Δγp. When life evolves from the first cycle to the 0.5Nfcycle,the plastic strain amplitude Δγpof the turning processes gets more extensive (as shown in Table 3).However, the grinding process shows the opposite result,

    where the plastic strain amplitude Δγpdecreases from 0.0056 to 0.0052 MPa·mm/mm, and the reverse cyclic hardening (as shown in the green curve in Fig. 13(b)) dominates the decrement of plastic strain amplitude.

    The comparison of ΔWpin multistage processes in Fig. 15 indicates that the all plastic area has decreased.The reciprocating cyclic motion annihilates the lattice friction,foreign atoms,etc.,and eventually reduces the friction energy24,25.The plastic properties mainly included friction energy and the stored energy. Besides, the decrement of friction energy causes the total plastic strain energy density to decrease.As per Giordana et al.26, the friction stress of ferrite/martensitic steel EUROFER 97 decreased with the decrement of the free dislocation density in the sub-grain, resulting in cyclic softening room temperature.

    However, the grinding process has the most significant decrement compared with the turning process. This can be attributed to stored energy produced in the grinding specimens in the machining process. The compressive residual stress resulted in the increment of the reversed back stress by the microstructural barriers or strain incompatibilities in the material27. The generation of the compressive residual stress resulted in the fact that the total plastic deformation ΔWpdecreases significantly at 0.5 Nffatigue life (Table 3). When the decrement of plastic deformation is dominant for Eq. (5),the total plastic strain energy density decreases significantly.

    3.2.3. Dislocation morphology characteristics Fig. 16 shows the dislocations distribution passing fatigue experiments of different machining processes. As shown in Fig. 16(a), there are apparent high-density dislocation tangles at two-phase areas between ferrite and cementite. In addition,in some ferrite grains with more significant plastic deformation, well-developed dislocation walls and dislocation cells are presented (Fig. 16(c)) compared with the GFRT process.This phenomenon corresponds with the increment of plastic energy (as indicated in the red line of Fig. 15).

    In the FRT0 process, the dislocations are arranged in the form of plane slip bands in ferrite grains(Fig.16(e)).Moreover,dislocation slip meets many carbide grains. The detour phenomenon occurs, which indicated that the plastic deformation degree of ferrite grain becomes lower, as shown in Fig.16(f).

    TEM observations of failed specimens in the FRT process are presented in Fig. 17, and dislocation configurations occur with dislocation tangles in two adjoining boundaries (Fig. 17(c)). Specifically, the ill-developed dislocation loop is showed in the ferrite grains in Fig. 17(a). It indicates that the degree of plastic deformation in ferrite becomes lower, attributed to the decrement of plastic energy density (as indicated in the black line of Fig. 15). However, in dislocation configurations by dislocation tangles at the grain boundary becomes unclear.There is a lower dislocation density. It tends to show transgranular behaviour from dislocation configurations at the grain boundary.Consequently,the GFRT process has the lowest plastic energy density compared with other processes.

    3.3. Fatigue behaviour evolution

    3.3.1. Fatigue life Machined surface layer evolution ultimately affects the torsional fatigue performance, and the fatigue life is shown in Table 4.It can be seen that the RT process has the shortest life of 5704 cycles due to the rougher surface.When the FRT process under wet cutting conditions was introduced, the average fatigue life increases by 50%, indicating that the surface morphology and the plastic deformation layer strongly affect the fatigue life. However, the fatigue life for the FRT0 process under dry conditions increases by 20% compared with the RT process. As the final process, the Grinding process has the highest average fatigue life of 13,950 cycles increasing by 1.45 times compared with the RT process. Aboulkhair NT et al.28suggested that the machined process does not remarkably improve fatigue life at high stress levels regardless of heat treatment. The literature18shows that specimens by laser shock peening have a 50% higher fatigue strength compared with the machined specimens. However, a factor of nearly ten times is obtained at a low-stress level.

    3.3.2. Fatigue fracture analysis

    The behaviour of 45CrNiMoVA steel is ductile, and shows a shear failure mechanism. The fracture of torsional fatigue occurs on the maximum shear plane and is perpendicular to the axial direction of the specimen. The surface morphology of fatigue fracture in multistage machined processes was presented from Figs.18-20.Fatigue crack initiates on the surface,and then propagates to the inside with obvious microplastic deformation trace (Fig. 18(e) and Fig. 20(b)). Eventually, the fracture occurs in the centre of Fig.18(c)and Fig.20(b).Interestingly, it is observed that many fatigue fracture surfaces are located at the arc transition near the clamping end.

    Table 4 Effect of surface integrity evolution on fatigue life.

    Fatigue fracture in the RT process is presented in Fig.18.It is shown from Fig. 18(d) that there are apparent bulges with extremely thin tearing slices due to shear slip. In the crack propagation area located at 1/2 R of the surface (Fig. 18(b)),the obvious microplastic deformation trace is perpendicular to the radial direction,and Fig.18(e)shows the magnified view Fig.18(b).The microplastic deformation trace in the whole circumferential region presents a circular distribution. Fig. 18(c)illustrates the transient breaking area located in the centre of the specimen, and Fig. 18(f) shows the magnified view of Fig. 18(c). It can be observed that the transient breaking area presents a circular shape, and obvious cleavage steps are presenting the typical characteristic of the final fracture.

    Fig. 19 shows the fatigue fracture of the FRT and FRT0 processes. Fig. 19(a), Fig. 19(b) show the two crack initiation sites on the surface of the FRT0 under dry cutting conditions,and Fig.19(d),Fig.19(e)illustrate respectively the corresponding magnified views. It can be seen that the surface crack sites of the FRT0 process also present bulges which are similar to that of the RT process, and the tear line with a certain depth can be seen at the bottom of the bulges (presented by the red line in Fig. 19(d), Fig. 19(e)).

    However, when SEM observes the entire circumferential surface of the FRT process, no noticeable bulges are presenting on the surface of the FRT process under wet cutting conditions (Fig. 19(c), Fig. 19(f)). In addition, a deeper shear line is observed at the flat shear plane,and the crack initiation area is smoother. In the FRT0 and RT processes, microplastic deformation trace in the propagation area, and the cleavage steps in the transient breaking area are similar to the fracture morphology of the FRT process,respectively,so it is not illustrated any further in this series.

    Fig. 20 shows the fatigue fracture of the grinding process.Fig.20(a)illustrates the fatigue crack initiation site on the surface,and Fig.20(d)is the magnified view of the crack initiation site.The crack initiation region is flatter,and there is no bulge similar to that of the RT and FRT0 process, where it presents more obvious transgranular characteristics.A noticeable shear microplastic deformation trace around the crack site, is perpendicular to the radial direction (Fig. 20(b)). However, the transient breaking region transfers from the regular circular area (dotted line in Fig. 18(c)) of the turning process to the irregular ring area (presented by the yellow line in Fig. 20(c))of the grinding process. Moreover, apparent cleavage steps are presenting a typical characteristic of final fracture. The irregular ring reduces the crack propagation rate and contributes to the torsional fatigue life improvement. At the same time, the conical dimples can also be seen in the lower right corner of the instantaneous breaking region.

    3.4. Effects of surface integrity evolutions on torsional fatigue behaviour

    The axial or circumferential residual stress shows the same trend with fatigue life when it evolves from the FRT process to the GFRT process. When the residual stress evolves from the dry condition to the wet condition,the circumferential residual stress which is almost unchanged,improves fatigue life.This is probably attributed to the decrement of surface morphology or plastic deformation layer.

    The microhardness in multistage machining processes can’t strongly reflect the fatigue performance because the microhardness variation range (0-30 HV) in Eq. (8) changes little.

    3.5. Fatigue life prediction model based on energy method considering machined surface integrity

    In section 3.4, the relationship between surface integrity evolution and fatigue properties is discussed from the aspects of geometric-mechanical-metallurgical characteristics. Considering that surface roughness Rahas the same evolution tend with the plastic deformation layer depth and makes the measurement more direct, the revised model considering the machined surface integrity is proposed.

    From the perspective of the physical meaning of the energy method, the increment of surface roughness contributes to the fatigue fracture occurred prematurely, which is mainly reflected by the decrease of the total energy.Thus,the geometric parameter of BRais added. On the one hand, considerable compressive residual stress reduces the plastic strain energy density in one cycle,while,changing the total plastic properties due to the rearrangement and changes in dislocations microstructure. It should be noted that the plastic energy density caused by residual stress has been considered in the Morrow model. Thus the mechanical characteristics of Cτ2resis added. The parameters calculated by the Morrow model and the proposed revised model are listed in Table 5, respectively.

    Fig.24(a)shows the fatigue prediction of the revised model considering surface integrity evolution. It is evident that the prediction life of the revised model is closer to the experimental data than the Morrow model. Fig. 24(b) compares the results of the fatigue life prediction model of the Morrow energy model and the revised model, where the life dispersion band represents the deviation of the predicted life from the experimental life.The prediction results of the Morrow energy model are in the 30% dispersion band, while the modified model results are in the 10% dispersion band, indicating that the model considering surface integrity parameters has higher prediction accuracy. The Morrow model’s minimum and average prediction precisions are 67.6%and 85.4%,respectively,while the minimum and average prediction precisions of the revised model considering surface integrity are 89.0%and 93.0%.The precision of the revised energy model is higher than that of the traditional Morrow energy model. Fatigue life predictionmodel based on energy method considering machined surface integrity shows a more explicit physical meaning.

    Table 5 Prediction model parameters.

    4. Conclusions

    Low-cycle torsional fatigue experiments were carried out on 45CrNiMoVA steel bars to explore the strain energy relationship between machined surface layer evolution and torsional fatigue behaviour of high strength steel. An energy-based model for the prediction of the torsional fatigue life considering the machined surface integrity was proposed, as follows:

    (1) Plastic deformation layer variation with multistage machining processes substantially impacts the lowcycle torsional fatigue life of 45CrNiMoVA steel bars,and it is consistent with a surface roughness Rain geometric characteristics. However, microhardness variation in metallurgical characteristics fails to affect the torsional fatigue life of 45CrNiMoVA steel. When the shallower plastic deformation layer and the smoother surface of the FRT processes evolve from the RT processes, the fatigue crack initiation feature evolves from some specific bulges to no raised shear plane characteristic. There is an obvious indication of high density dislocation tangles at two-phase areas between ferrite and cementite. The plastic strain energy density of one cycle change little at 0.5Nffatigue life, and the average torsional fatigue life of the FRT process increases by 50%.

    (2) Residual stress variation in mechanical characteristics significantly impacts the torsional fatigue life of 45CrNi-MoVA steel bars by changing the hysteresis loop. Cycle characteristics transform from softening to hardening when more significant compressive residual stress of the GFRT processes evolves from the tensile residual stress of the FRT process. Plastic strain amplitude decreases by 21%, and plastic strain energy density per cycle reduces 4% at 0.5Nffatigue life, and the average low-cycle torsional fatigue life of the GFRT process increases by 63.2%.

    (3) An energy-based model for the predicting torsional fatigue life of 45CrNiMoVA steel bars considering the machined surface integrity was proposed.Surface geometry characteristics were modelled using surface roughness Raand residual stress τres, considering that the plastic deformation layer depth in metallurgical characteristics is consistent with a surface roughness Rain geometric characteristics. It has a significant physical significance that surface topography contributes to the fatigue fracture that occurred prematurely, and is mainly reflected in the decrement of total energy. In addition, residual stress mainly contributes to the decrement of plastic energy density in one cycle.The mechanical characteristics are added to the total energy considering the rearrangement and changes in dislocations microstructure.

    (4) In a revised model of the prediction fatigue based on the energy method,the dispersion band of the revised model results with experimental life transfers from 30%to 10%compared with the conventional Morrow energy model.This model can be applied to precisely predict the fatigue life of 45CrNiMoVA steel when variation in surface geometry characteristics is generated in different machined processes. The model result indicates that the influence of the compressive residual stress on fatigue life is mainly attributed tothe plastic strain density of one cycle due to a smaller fair value.

    Declaration of Competing Interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgement

    This work supported by the National Natural Science Foundation of China (No. 52075042).

    国产精品 国内视频| 国产免费又黄又爽又色| 黄色毛片三级朝国网站| 久久久精品国产亚洲av高清涩受| 丝袜在线中文字幕| 99久久综合免费| 在线观看三级黄色| 亚洲色图 男人天堂 中文字幕| 亚洲欧美色中文字幕在线| 97人妻天天添夜夜摸| 国产成人啪精品午夜网站| 日韩大码丰满熟妇| 国产精品久久久久久人妻精品电影 | 最近手机中文字幕大全| 亚洲免费av在线视频| 看十八女毛片水多多多| 女人高潮潮喷娇喘18禁视频| 日韩 亚洲 欧美在线| 成人影院久久| 制服诱惑二区| 国产成人精品在线电影| kizo精华| 99久久人妻综合| 18禁观看日本| 亚洲图色成人| 9热在线视频观看99| 中文字幕高清在线视频| svipshipincom国产片| 欧美日韩视频精品一区| 国产精品av久久久久免费| 久久久精品免费免费高清| 丰满饥渴人妻一区二区三| av不卡在线播放| 亚洲,一卡二卡三卡| 国产成人免费观看mmmm| 久久久精品国产亚洲av高清涩受| 国产成人系列免费观看| 免费黄网站久久成人精品| av电影中文网址| 看免费av毛片| 久久久久国产一级毛片高清牌| 一级毛片黄色毛片免费观看视频| 精品一区二区三卡| 国产一区二区 视频在线| 国产一区二区激情短视频 | 国产免费又黄又爽又色| 国产精品免费视频内射| 丰满少妇做爰视频| 欧美少妇被猛烈插入视频| 99久久精品国产亚洲精品| 亚洲av中文av极速乱| 如日韩欧美国产精品一区二区三区| 天天添夜夜摸| 日韩熟女老妇一区二区性免费视频| 十八禁网站网址无遮挡| 亚洲国产欧美日韩在线播放| 中文乱码字字幕精品一区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 免费少妇av软件| 激情视频va一区二区三区| 日韩欧美精品免费久久| 在线观看国产h片| 日本vs欧美在线观看视频| 久久99精品国语久久久| 亚洲成人av在线免费| 美女高潮到喷水免费观看| 精品一区二区三卡| 亚洲精品国产区一区二| 一本大道久久a久久精品| 亚洲精品av麻豆狂野| 国产亚洲午夜精品一区二区久久| 香蕉丝袜av| 黄色一级大片看看| 久久精品aⅴ一区二区三区四区| 99热全是精品| netflix在线观看网站| 人人澡人人妻人| 又粗又硬又长又爽又黄的视频| 交换朋友夫妻互换小说| 午夜福利免费观看在线| 国产av一区二区精品久久| 亚洲少妇的诱惑av| 国产精品久久久久久精品电影小说| 国产精品免费大片| 国产麻豆69| 久久精品亚洲av国产电影网| 亚洲一码二码三码区别大吗| 制服丝袜香蕉在线| 大话2 男鬼变身卡| 老司机影院毛片| 精品国产一区二区三区久久久樱花| 免费在线观看完整版高清| 波多野结衣av一区二区av| 老熟女久久久| 国产野战对白在线观看| 欧美精品高潮呻吟av久久| 亚洲伊人色综图| 日韩 亚洲 欧美在线| 亚洲三区欧美一区| 精品国产超薄肉色丝袜足j| 曰老女人黄片| 久久久久网色| 啦啦啦在线免费观看视频4| 午夜免费鲁丝| av天堂久久9| 两个人看的免费小视频| 亚洲情色 制服丝袜| 欧美精品av麻豆av| 欧美国产精品一级二级三级| 美女扒开内裤让男人捅视频| 最新在线观看一区二区三区 | 超碰97精品在线观看| 亚洲欧美中文字幕日韩二区| 精品久久久精品久久久| 一本大道久久a久久精品| 亚洲自偷自拍图片 自拍| 人妻人人澡人人爽人人| 精品国产国语对白av| 大片免费播放器 马上看| 午夜福利影视在线免费观看| 欧美日韩av久久| a 毛片基地| 国产片特级美女逼逼视频| 精品少妇久久久久久888优播| 国产精品亚洲av一区麻豆 | 午夜91福利影院| 欧美久久黑人一区二区| 国产精品一区二区在线不卡| 男女国产视频网站| 日本欧美视频一区| av在线老鸭窝| kizo精华| 久久鲁丝午夜福利片| 啦啦啦 在线观看视频| 99国产精品免费福利视频| a级片在线免费高清观看视频| 97在线人人人人妻| 国产亚洲av片在线观看秒播厂| 丰满饥渴人妻一区二区三| 狠狠精品人妻久久久久久综合| 亚洲av综合色区一区| 亚洲成人一二三区av| 久久久久国产精品人妻一区二区| 最近手机中文字幕大全| 欧美乱码精品一区二区三区| 日本爱情动作片www.在线观看| 男女国产视频网站| 成年女人毛片免费观看观看9 | 日韩中文字幕视频在线看片| 十分钟在线观看高清视频www| 成人黄色视频免费在线看| 两性夫妻黄色片| 国产精品无大码| 亚洲一级一片aⅴ在线观看| 又大又黄又爽视频免费| 日本色播在线视频| 国产精品成人在线| 国产黄频视频在线观看| 国产精品无大码| 大香蕉久久成人网| 老鸭窝网址在线观看| 亚洲精品视频女| 美女大奶头黄色视频| 亚洲精品aⅴ在线观看| 狠狠婷婷综合久久久久久88av| 亚洲人成77777在线视频| 不卡视频在线观看欧美| 十八禁网站网址无遮挡| 在线天堂中文资源库| 精品国产一区二区三区久久久樱花| 国产毛片在线视频| 岛国毛片在线播放| 中文字幕另类日韩欧美亚洲嫩草| 纵有疾风起免费观看全集完整版| 国产成人精品久久二区二区91 | 国产 一区精品| 最黄视频免费看| 中文精品一卡2卡3卡4更新| 国产精品三级大全| 久久午夜综合久久蜜桃| av卡一久久| 又黄又粗又硬又大视频| 老司机影院成人| 国产日韩欧美视频二区| 日韩中文字幕欧美一区二区 | 热re99久久国产66热| 亚洲成人手机| 最近手机中文字幕大全| 国产精品国产三级国产专区5o| 99久国产av精品国产电影| 午夜福利影视在线免费观看| 国产精品嫩草影院av在线观看| 99久久99久久久精品蜜桃| 一区二区三区乱码不卡18| 亚洲精品日本国产第一区| 侵犯人妻中文字幕一二三四区| 国产在视频线精品| 欧美日韩av久久| 美女视频免费永久观看网站| 九九爱精品视频在线观看| 亚洲精品一二三| 国产精品 国内视频| 日韩av免费高清视频| 日韩 亚洲 欧美在线| 久久毛片免费看一区二区三区| 午夜免费男女啪啪视频观看| 国产午夜精品一二区理论片| 人人妻人人澡人人爽人人夜夜| 波野结衣二区三区在线| 国产极品天堂在线| 精品国产一区二区久久| 国产人伦9x9x在线观看| 99久久99久久久精品蜜桃| 丝袜脚勾引网站| 悠悠久久av| 99久久人妻综合| 国产不卡av网站在线观看| 国产免费福利视频在线观看| 亚洲av男天堂| 久久国产亚洲av麻豆专区| a级毛片黄视频| 18禁观看日本| 欧美黑人欧美精品刺激| 亚洲七黄色美女视频| 最新的欧美精品一区二区| 女人精品久久久久毛片| 亚洲精品久久午夜乱码| av不卡在线播放| 超碰97精品在线观看| 天堂俺去俺来也www色官网| 99久久人妻综合| 久热这里只有精品99| 免费在线观看完整版高清| 亚洲av电影在线进入| 黄频高清免费视频| 中文欧美无线码| 久久久久精品久久久久真实原创| 最新的欧美精品一区二区| 又大又黄又爽视频免费| 大片电影免费在线观看免费| 99久久综合免费| 少妇 在线观看| 国产一区有黄有色的免费视频| 国产淫语在线视频| 亚洲在久久综合| 国产av码专区亚洲av| 亚洲人成电影观看| 亚洲av福利一区| 欧美日韩国产mv在线观看视频| 色婷婷久久久亚洲欧美| 999精品在线视频| 午夜影院在线不卡| 捣出白浆h1v1| 母亲3免费完整高清在线观看| 黄片无遮挡物在线观看| 久久这里只有精品19| 成人18禁高潮啪啪吃奶动态图| 我的亚洲天堂| 亚洲精品日本国产第一区| 性色av一级| 在线观看一区二区三区激情| 亚洲久久久国产精品| 精品福利永久在线观看| 别揉我奶头~嗯~啊~动态视频 | 黄色毛片三级朝国网站| 免费黄色在线免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品国产亚洲av涩爱| 人成视频在线观看免费观看| 日本色播在线视频| 久久久久人妻精品一区果冻| 亚洲精品国产av成人精品| 女的被弄到高潮叫床怎么办| 国产免费又黄又爽又色| 亚洲精品美女久久av网站| 人成视频在线观看免费观看| 亚洲一卡2卡3卡4卡5卡精品中文| 观看av在线不卡| 国产男女内射视频| 好男人视频免费观看在线| a级毛片黄视频| 韩国精品一区二区三区| 1024香蕉在线观看| 免费观看人在逋| 欧美人与性动交α欧美精品济南到| 欧美日韩视频精品一区| 国产精品欧美亚洲77777| 不卡av一区二区三区| 三上悠亚av全集在线观看| 高清视频免费观看一区二区| 美女大奶头黄色视频| 啦啦啦在线观看免费高清www| 欧美精品av麻豆av| 成年人午夜在线观看视频| av线在线观看网站| 人人妻人人爽人人添夜夜欢视频| 免费观看av网站的网址| 国产欧美日韩综合在线一区二区| 男人操女人黄网站| 日韩制服骚丝袜av| 啦啦啦中文免费视频观看日本| 欧美国产精品va在线观看不卡| 日韩制服丝袜自拍偷拍| 90打野战视频偷拍视频| 看非洲黑人一级黄片| 在线观看国产h片| 97精品久久久久久久久久精品| 亚洲国产精品国产精品| 欧美人与善性xxx| 久久免费观看电影| 亚洲 欧美一区二区三区| 少妇 在线观看| 国产淫语在线视频| 香蕉丝袜av| 最近最新中文字幕免费大全7| 多毛熟女@视频| 国产 精品1| 国产av一区二区精品久久| 亚洲av日韩在线播放| 最黄视频免费看| 久久影院123| 欧美日韩精品网址| 国产日韩欧美亚洲二区| 免费高清在线观看日韩| av网站在线播放免费| 亚洲精品第二区| 日韩一区二区视频免费看| 国产免费福利视频在线观看| 18禁裸乳无遮挡动漫免费视频| 亚洲精品国产色婷婷电影| 国产精品久久久久久精品电影小说| 天堂俺去俺来也www色官网| 另类亚洲欧美激情| 麻豆av在线久日| 国产男人的电影天堂91| 一边摸一边抽搐一进一出视频| 99热全是精品| 亚洲精华国产精华液的使用体验| 国产成人精品在线电影| 99香蕉大伊视频| 啦啦啦视频在线资源免费观看| 一边摸一边抽搐一进一出视频| 欧美日韩亚洲高清精品| 精品少妇黑人巨大在线播放| 久久久精品区二区三区| 中文字幕制服av| 国产野战对白在线观看| 19禁男女啪啪无遮挡网站| 久久女婷五月综合色啪小说| 日本wwww免费看| 王馨瑶露胸无遮挡在线观看| 乱人伦中国视频| 性高湖久久久久久久久免费观看| 亚洲精品久久久久久婷婷小说| 中文字幕色久视频| 又黄又粗又硬又大视频| 搡老岳熟女国产| 久久人妻熟女aⅴ| 一本色道久久久久久精品综合| 日日爽夜夜爽网站| 女的被弄到高潮叫床怎么办| 亚洲精品美女久久久久99蜜臀 | 男女无遮挡免费网站观看| 在线看a的网站| 男女边吃奶边做爰视频| 韩国av在线不卡| 黄色一级大片看看| 亚洲美女视频黄频| 久久精品亚洲熟妇少妇任你| 午夜福利乱码中文字幕| 国产极品天堂在线| 国产高清国产精品国产三级| 蜜桃国产av成人99| 久久鲁丝午夜福利片| 日本wwww免费看| 一区二区三区乱码不卡18| 国产午夜精品一二区理论片| 亚洲色图 男人天堂 中文字幕| 亚洲国产精品一区三区| 黄色 视频免费看| 免费女性裸体啪啪无遮挡网站| e午夜精品久久久久久久| 老司机在亚洲福利影院| 熟女av电影| 操出白浆在线播放| 国产精品久久久久成人av| 十八禁网站网址无遮挡| 国产精品免费大片| 亚洲精品在线美女| 色94色欧美一区二区| 国产精品亚洲av一区麻豆 | 91成人精品电影| 亚洲精品aⅴ在线观看| 久久久久人妻精品一区果冻| 亚洲成人一二三区av| 黑丝袜美女国产一区| 久久综合国产亚洲精品| av线在线观看网站| 成人毛片60女人毛片免费| 青青草视频在线视频观看| 午夜精品国产一区二区电影| 国产av码专区亚洲av| 国产精品欧美亚洲77777| 如何舔出高潮| 最近中文字幕2019免费版| 国产97色在线日韩免费| av女优亚洲男人天堂| 赤兔流量卡办理| 男女边摸边吃奶| 狂野欧美激情性xxxx| a级毛片黄视频| 国产成人一区二区在线| 两性夫妻黄色片| 国产免费一区二区三区四区乱码| 美女国产高潮福利片在线看| 欧美激情 高清一区二区三区| 只有这里有精品99| 欧美久久黑人一区二区| 亚洲精品一区蜜桃| 欧美在线一区亚洲| 少妇猛男粗大的猛烈进出视频| 高清在线视频一区二区三区| 亚洲激情五月婷婷啪啪| 午夜日韩欧美国产| 国产毛片在线视频| 可以免费在线观看a视频的电影网站 | 中文字幕最新亚洲高清| 九色亚洲精品在线播放| 亚洲,欧美,日韩| 美女主播在线视频| 午夜日韩欧美国产| 日本欧美视频一区| 国产精品国产三级专区第一集| 亚洲欧洲日产国产| 黄网站色视频无遮挡免费观看| 精品一品国产午夜福利视频| 麻豆av在线久日| 日本爱情动作片www.在线观看| 男人爽女人下面视频在线观看| 在线亚洲精品国产二区图片欧美| 亚洲视频免费观看视频| 亚洲熟女精品中文字幕| av在线播放精品| 久久精品国产a三级三级三级| 精品酒店卫生间| 色婷婷av一区二区三区视频| 夜夜骑夜夜射夜夜干| 久热这里只有精品99| 亚洲熟女精品中文字幕| 亚洲五月色婷婷综合| 18禁动态无遮挡网站| 国产一区二区激情短视频 | 精品国产国语对白av| 国产色婷婷99| 十八禁人妻一区二区| 女人爽到高潮嗷嗷叫在线视频| 欧美成人精品欧美一级黄| 亚洲av成人精品一二三区| 国语对白做爰xxxⅹ性视频网站| 亚洲美女搞黄在线观看| 国产精品av久久久久免费| 国产精品三级大全| 啦啦啦啦在线视频资源| 国产免费现黄频在线看| 在线观看三级黄色| 麻豆精品久久久久久蜜桃| 亚洲中文av在线| 国产精品人妻久久久影院| 少妇人妻精品综合一区二区| 在线天堂中文资源库| 成人国产av品久久久| 亚洲精品国产一区二区精华液| 超碰成人久久| 国产不卡av网站在线观看| 国产国语露脸激情在线看| 国产精品免费视频内射| 国产成人精品福利久久| 成人三级做爰电影| 丝袜在线中文字幕| 亚洲欧美一区二区三区黑人| 亚洲国产精品一区三区| 日韩制服丝袜自拍偷拍| 日韩大码丰满熟妇| 王馨瑶露胸无遮挡在线观看| 久久久欧美国产精品| 国产麻豆69| 99热全是精品| 欧美日韩综合久久久久久| 国产亚洲午夜精品一区二区久久| 一区二区三区激情视频| 免费女性裸体啪啪无遮挡网站| 黑人猛操日本美女一级片| 飞空精品影院首页| 天堂中文最新版在线下载| 我的亚洲天堂| 人妻 亚洲 视频| av片东京热男人的天堂| 亚洲熟女毛片儿| 狠狠婷婷综合久久久久久88av| 国产亚洲欧美精品永久| 婷婷色av中文字幕| 天堂中文最新版在线下载| 日本av手机在线免费观看| 国产亚洲av片在线观看秒播厂| 亚洲国产精品一区二区三区在线| 欧美国产精品一级二级三级| 亚洲专区中文字幕在线 | 国产激情久久老熟女| 999久久久国产精品视频| 久久99一区二区三区| 91精品伊人久久大香线蕉| 久久毛片免费看一区二区三区| 午夜av观看不卡| 五月天丁香电影| 国产av码专区亚洲av| 精品一区二区三区av网在线观看 | 亚洲精品国产色婷婷电影| 极品少妇高潮喷水抽搐| 欧美精品一区二区免费开放| 久久精品亚洲熟妇少妇任你| 午夜福利在线免费观看网站| 精品第一国产精品| 国产日韩欧美在线精品| 婷婷色综合大香蕉| 国产 精品1| 国产精品av久久久久免费| 国产成人av激情在线播放| 成年动漫av网址| 亚洲免费av在线视频| 国产成人欧美| 美女大奶头黄色视频| 91精品伊人久久大香线蕉| 亚洲,欧美精品.| 国产精品三级大全| 久久精品熟女亚洲av麻豆精品| 欧美日韩成人在线一区二区| 女性被躁到高潮视频| 亚洲成色77777| 制服人妻中文乱码| 男人添女人高潮全过程视频| 日韩大片免费观看网站| 18禁观看日本| 欧美人与性动交α欧美精品济南到| 国产1区2区3区精品| 另类亚洲欧美激情| av又黄又爽大尺度在线免费看| 亚洲三区欧美一区| 纵有疾风起免费观看全集完整版| xxxhd国产人妻xxx| 久久毛片免费看一区二区三区| 日本av免费视频播放| 一边摸一边抽搐一进一出视频| 亚洲久久久国产精品| 午夜激情久久久久久久| 熟妇人妻不卡中文字幕| 欧美精品一区二区免费开放| 亚洲精品成人av观看孕妇| 久久99一区二区三区| 欧美乱码精品一区二区三区| www.av在线官网国产| 国产片内射在线| avwww免费| xxx大片免费视频| 老熟女久久久| 日本一区二区免费在线视频| 婷婷色av中文字幕| 欧美日韩成人在线一区二区| 免费黄色在线免费观看| 欧美精品一区二区大全| 亚洲四区av| 香蕉丝袜av| 一区二区三区乱码不卡18| 成人免费观看视频高清| 亚洲色图综合在线观看| 国产精品二区激情视频| 亚洲av电影在线观看一区二区三区| 亚洲精品国产av蜜桃| 美女主播在线视频| 三上悠亚av全集在线观看| 欧美日韩一级在线毛片| 晚上一个人看的免费电影| 日本av免费视频播放| 亚洲成人一二三区av| 在线观看免费午夜福利视频| 少妇精品久久久久久久| 大码成人一级视频| 日韩视频在线欧美| av网站在线播放免费| 日韩视频在线欧美| www.熟女人妻精品国产| 人妻 亚洲 视频| 亚洲av电影在线观看一区二区三区| 无遮挡黄片免费观看| 精品久久蜜臀av无| 搡老岳熟女国产| av网站在线播放免费| 久久婷婷青草| 男女高潮啪啪啪动态图| 亚洲精华国产精华液的使用体验| 熟女少妇亚洲综合色aaa.| 欧美黑人精品巨大| www.精华液| 欧美老熟妇乱子伦牲交| 国产成人av激情在线播放| 亚洲,欧美,日韩| 日本爱情动作片www.在线观看| 国产精品久久久人人做人人爽| 亚洲欧洲日产国产| 亚洲精品久久成人aⅴ小说| 美女扒开内裤让男人捅视频| 国产午夜精品一二区理论片| 黄片播放在线免费| 老鸭窝网址在线观看| 国产无遮挡羞羞视频在线观看| 久久久亚洲精品成人影院|