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

    Metal flow hysteresis behavior and shape control strategy during porthole die extrusion of micro heat pipe

    2023-09-05 09:44:38YongdaLIUJieXUXiaoliangWANGDebinSHANBinGUO
    CHINESE JOURNAL OF AERONAUTICS 2023年7期

    Yongda LIU, Jie XU, Xiaoliang WANG, Debin SHAN, Bin GUO

    Key Laboratory of Micro-Systems and Micro-Structures Manufacturing of Ministry of Education, Harbin Institute of Technology, Harbin 150080, China

    School of Materials Science and Engineering, Harbin Institute of Technology, Harbin 150001, China

    KEYWORDS Micro heat pipe;Microforming;Size effect;Porthole die extrusion;Metal flow;Forming integrity;Microstructure

    Abstract The shape control strategy of micro grooves is still unclear and challenging during the porthole die extrusion of grooved micro heat pipe(MHP).Through the simulation and experiment of porthole die extrusion of a MHP profile,the metal flow hysteresis behavior within micro features and the effect of ram speed and extrusion temperature on it and the resulting forming integrity was elucidated.Innovatively,Taguchi design and variance analysis(ANOVA)were introduced to determine their influence magnitude on the metal flow uniformity calculated by simulation results.The main findings are given below.The metal flow hysteresis derives from part feature size effect.The negligible friction-affected area during conventional extrusion severely slows down the metal flow within micro features during the MHP profile extrusion, which is due to the surge in the area ratio of the friction-affected area to the region in which it is located.Neither ram speed nor extrusion temperature can change the distribution of the friction-affected area.However, increasing ram speed multiplies the metal flow hysteresis and severely reduces the forming integrity,whereas extrusion temperature has little effect.Following this strategy,batch extrusion of the profile with microgrooved width of 0.27 ± 0.02 mm was achieved in industrialized conditions.

    1.Introduction

    Heat pipes are heat transfer passive devices that operate on the principle of vapor–liquid phase change and are suitable for aerospace applications and heat dissipation in electronics due to their advantages of high heat transfer performance and compactness.1The capillary wick structure of heat pipes provides the driving force for the liquid in the condenser to flow back to the evaporator, which is crucial for the proper operation of heat pipes.A wick consists mainly of sintered,grooved, and composite types, according to their structural characteristics.2In aerospace, axially grooved heat pipes(AGHPs) made of aluminum alloy and filled with ammonia as working fluid are favored in satellites and their thermal control systems,3,4which is due to their uniform temperature,excellent reliability, process consistency and lightweight.The structure of the AGHPs can be found from this literature of Hoa et al.3, and their operational temperature range of-40 °C to + 80 °C is well suitable for space applications.For example, back in 2003, Baturkin et al.5reported that U-shaped AGHPs have been applied to the elaboration and flight exploitation of thermal control systems for the microsatellites Magion 4, Magion 5 and BIRD.Shao et al.6disclosed that AGHPs helped achieve thermal coupling between the panels on the Chang’e-1 satellite, which solved the problem of high temperatures in the batteries.Regarding the lifetime of AGHPs, Barantsevich and Shabalkin7reported that AGHPs can serve for more than 15 years in the application of cable thermal control for the solar battery drives integrated into the International Space Station.However, with recent advances in space technology towards miniaturization and integration,8thermal control of spacecraft is more challenging than ever.In response to this situation,axially grooved micro heat pipes (MHPs) provide a miniaturized and lightweight thermal management solution.MHPs with crosssectional dimensions of a few millimeters show promising applications in microsatellites and microthermal control systems in aerospace.Baturkin9noted that micro heat pipes with an external circular diameter of 1–6 mm are critical for miniaturized designs of microsatellites.This type of heat pipe indeed exhibits attractive advantages, such as solving the watt-level heat dissipation issues for sensitive components of the charge coupled device (CCD) camera in narrow space in a satellite thermal control system.10

    However,fabrication techniques for MHPs of this type suffer from many problems, and one of the pressing issues is the shape control of the micro grooves during extrusion.The micro grooves of heat pipes are typically less than 1 mm in width and height, resulting in their fabrication falling into the category of microforming.Smaller groove widths of heat pipes represent higher capillary heat transfer limits,11but this also brings great troubles to its hot extrusion, such as the difficulty of fully extruding micro grooves.12According to prior studies,13,14hot extrusion of MHPs has proved to be unique and compelling.For example, complex micro groove and ultra-large extrusion ratio(i.e.,the ratio of cross-sectional area of containers to extrudates) of MHPs can cause many problems,such as extremely poor forming integrity owing to highly uneven metal flow13and severe surface back-end defects caused by rare forward flow behavior of billet skins.14Taking the forming quality of MHPs,which this present paper focuses on, as an example, Liu et al.13reported that in the hot extrusion of an Al MHP profile at an ultra-large extrusion ratio of 205,the micro grooves failed to extrude even at a very small ram speed of 1.4 mm/s,and the primary reason was attributed to the uneven metal flow within the die bearing caused by size effect.Evidently, an in-depth understanding of the metal flow behavior subjected to size effect during hot extrusion is crucial to ensure the forming quality of MHPs.However, to the best of our knowledge, little research has been reported on the hot extrusion technology of complex micro-grooved hollow profiles such as MHPs, and most studies of the size effectaffected metal flow behavior during extrusion and microextrusion used cold extrusion process.To date, the metal flow behavior during cold micro-extrusion,as well as some cases of extrusion process realizations of complex micro parts (mainly solid micro gears), has been widely investigated.

    Compared to conventional extrusion, grain size and workpiece size can cause uneven metal flow during micro-extrusion,which is generally detrimental to the forming quality of extrudates.Parasiz et al.15reported that as the grain size (211 μm)approached the workpiece size (568 μm), the deformation characteristics of the extruded micro pins made of CuZn30 α-brass were dominated by the size and location of specific grains, resulting in an uneven distribution of plastic strain and hardness, which caused the curving tendency.On this basis, Parasiz et al.16illustrated that during micro-extrusion,the penetration of shear deformation occurred for the coarse-grained material, and the imposed strain gradients became steeper with the decrease of workpiece size.Similarly,based on the micro-extrusion of pure coper with various grain sizes,Chan et al.17found that the inhomogeneous deformation occurred in the case with coarse grains, which showed extensive slip bands across the grain boundaries for strain continuity.Interestingly, the size effect is not entirely negative for micro-extrusion process.Experimental results by Ghassemali et al.18revealed that the dead metal zone in the extruded micro pins could be removed by increasing the initial grain size of the workpiece,which was attributed to the effect of the strain gradient on the deformation.Recently, further research on the interactive effect of grain size and crystal structure on the deformation behavior during the progressive microforming process including piercing, two step extrusion, and blanking was investigated by Tang et al.19They revealed that for facecentered cubic (FCC) Cu and body-centered cubic (BCC) Fe,the length of extrudates of the micropart was reduced with the increase of initial grain size owing to the enhanced lateral metal flow in the coarse-grained material, while the dependency of its length on the initial grain size for hexagonal close-packed(HCP)Ti became unnoticeable owing to the subdivision of the original grains by the generated deformation twins.

    Regarding the micro-extrusion process of complex micro parts, micro gears for micro-electro-mechanical systems(MEMS) parts are most usually studied.During the microextrusion of micro gears, the metal flow within micro features(e.g., the teeth of micro-gears) is idiosyncratic compared to other areas, which is similar to the hot extrusion process of MHPs.13Nanthakumar and Rajenthirakumar20recently reported that the hardness of the tooth region of the coldextruded micro gear made of cooper was higher than that of the central region.They believed that the most serious deformation occurred in the tooth region, resulting in the smallest grain size owing to low grain growth rate at room temperature.In the hot micro-extrusion of 7075 aluminum alloy micro gear,Dong et al.21indeed observed that the grain size at the tooth region was smaller than that at the central region, whereas the micro hardness of the tooth region was smaller.Based on the microstructure of the micro gears, it was attributed to the smaller amount of rigid second phase particles at the tooth region.In terms of improving the differences in microstructure distribution between the tooth and central regions of micro gears, Khalilpourazary22determined that fine-grained initial billets prepared by ECAP process were able to reduce the inhomogeneity of the microstructure within micro gears.

    In summary, despite numerous excellent reports on the metal flow behavior subjected to size effect in the microextrusion of micro parts,the reason for the uneven metal flow within micro features is still ambiguous and has not received much attention.More importantly, the influence of extrusion process parameters on the uneven flow behavior is still unclear,which is very important for the formulation of extrusion parameters and seriously restricts the realization of industrial production of micro extrudates and macro extrudates with micro features (e.g., MHPs).

    Following these above studies, the aim of this paper was established to present comprehensive simulated and experimental investigations of the formation mechanism of the uneven metal flow behavior within micro features during hot extrusion of MHPs and the influence magnitude of extrusion process parameters on it and the resulting forming integrity.In this study, the hot extrusion of an aluminum alloy MHP with 0.27 mm wide micro grooves was studied,and its porthole die with unequal length of bearings was designed for use in an industrial extruder.The metal flow behavior at the profile exit of the die during hot extrusion was obtained by HyperXtrude software, and the metal flow uniformity was calculated from the simulated results.Following this, Taguchi design and ANOVA were innovatively introduced to determine the influence magnitude of ram speed and extrusion temperature on the metal flow uniformity.To validate the simulation results, the experimental results including the head profile shape, forming integrity and microstructure of the MHP profile were analyzed.This study provides a crucial basis for the design of extrusion dies and the formulation of extrusion parameters for MHP-like extruded products in practical production using industrial extruders.

    2.Material and method

    2.1.Material and porthole die design

    The as-cast 6063 aluminum alloy was used as the material for the MHP profile, and its chemical composition is given in Table 1.It is most commonly used in the manufacture of MHPs due to good comprehensive properties such as mechanical, welding properties and corrosion resistance.(ED) was depicted to understand the extrusion process.Depending on the structural characteristics of the profile, the design idea was proposed with a relatively small extrusion ratio, short bearing, low weld chamber and local structural strengthening.The specific reasons were listed as follows:(1)Four die orifices were set,which sharply reduced the extrusion ratio from 548 to 137; (2) the overall bearing length was set to 3 mm.The bearing length in aluminum alloy extrusion is generally 2 to 8 mm.Longer bearing can aggravate the uneven metal flow at the profile exit of extrusion dies,whereas shorter bearing can cause periodic cracks on extruded products; (3) to prevent unacceptable deflection deformation of the die core, not only a small weld chamber height of 8 mm,but also a local strengthening structure(shown in Detail A)was designed.In summary,the main die parameters are the extrusion ratio of 137,die core spacing of 28 mm, weld chamber depth of 8 mm, overall bearing length of 3 mm and porthole outer radius of 30 mm.

    Further detailed design of the bearing length was performed to optimize the die structure in Fig.3.This is because numerous studies confirmed that bearing length design could balance the metal flow at the profile exit of extrusion dies.23–27Decreasing the local bearing length results in an increase in the local metal flow velocity, otherwise the opposite.23Generally,designer can slow the metal flow in thick sections of extrudates by increasing the bearing length.24Liu et al.25reported that using an unequal bearing length, the difference of the profile exit velocity (PEV) was decreased from 7.71 to 2.00 mm/s.Therefore, the unequal length of bearing(Scheme 2) was designed to perform the simulated and experimental studies.Considering that a short bearing increases the metal flow velocity at the profile exit, and that, to the best of the author’s knowledge, the bearing length of an aluminum alloy extrusion die is typically not less than 2 mm, otherwise the extrudates may be prone to periodic cracks or other defects.Therefore,a local bearing length of 2 mm in Scheme 2 is set for the micro rib area with difficult metal flow.

    2.2.Steady-state extrusion simulation

    To obtain the velocity and strain distribution at the profile exit of the die,the steady-state extrusion process of the profile was simulated using HyperXtrude based on Arbitrary Lagrange Eulerian (ALE) algorithm under the following extrusion process parameters: billet temperature of 480 °C, die temperature of 460 °C, container temperature of 420 °C, and ram speed of 0.7 mm/s.Both bearing design Schemes in Fig.3 were simulated to incidentally determine the effect of the optimal design of the bearing.

    The preprocessing procedure is that after meshing,entering material properties and applying boundary conditions, the solutions were performed with a nonlinear iteration tolerance of 0.001.Specifically, the billet diameter and length were designed to be ?85 and 200 mm, respectively.The metal flow domain in the container and porthole die was geometrically modeled and discretized by meshing to simulate the steady-

    Fig.1 shows the geometric model and key dimensions of the MHP profile.Obviously, the ultra-large extrusion ratio,complex micro groove and small inner radius are the most prominent characteristics that need to be focused on when designing the porthole extrusion die.A schematic diagram of the extrusion process using a porthole die can be found in the literature by Liu et al.13.

    A porthole die was designed to extrude the MHP profile on an industrial horizontal extruder, and one billet can generally extrude tens of meters of the profile.The overall design of the porthole die is shown in Fig.2, and the extrusion directionstate extrusion process, as shown in Fig.4.The bearing and profile are triangular prism elements, and the others are tetrahedral elements.Notably, the mesh of the micro rib must be locally refined to ensure the simulated accuracy.As shown in the enlarged Detail B view in Fig.4, at least three mesh elements are divided on each edge.For the finite element (FE)model, the minimum element size of the profile, bearing and weld chamber is 0.08 mm, and that of the porthole and billet is 0.5 and 1.5 mm, respectively.The maximum element size of 7 mm comes from the billet.With regard to the simulation parameters and boundary conditions,refer to the study by Liu et al.14.

    Table 1 Chemical composition of 6063 aluminum alloy.

    Fig.1 Key dimensions and geometric model of MHP profile.

    Fig.2 3D model of the overall design of porthole die.

    Fig.3 Bearing design of the porthole die.

    Fig.4 FE model of metal flow domain in hot extrusion of MHP profile.

    Exploratively, this study attempted to simulate using the conventional material constitutive model rather than that obtained from microforming tests (e.g., micro-compression),and discussed the reliability in combination with the experimental results.The specific reasons for this choice are as follows: On the one hand, Chan et al.28found that, using the calibrated flow stress–strain curves obtained from microcompression, the simulated results in a cold micro-double cup extrusion of 6061 aluminum alloy were in good agreement with their experimental results.However,the flow stress–strain curve obtained from micro-compression is clearly not suitable for the extrusion simulation of MHP profiles because it contains not only micro features but also macro features.On the other hand, due to the complex microstructure evolution during porthole die extrusion of MHP profiles,the size effect during the hot extrusion of micro grooves of MHP profiles may be different from that in the case of cold micro-extrusion.

    2.3.Taguchi design for extrusion process parameters

    Combined with the simulated results, Taguchi design and ANOVA were used to analyzed the influence magnitude of ram speed and extrusion temperature on the metal flow uniformity during hot extrusion of the profile.Taguchi design typically requires only a portion of a full factorial experiment,but it ensures balanced comparisons of levels of any one factor.Through analyzing the results of signal-to-noise ratio (referred to as S/N for short) and ANOVA in Taguchi design,Wang et al.29successfully determined the contribution rate and optimal parameter combination of five control factors to the measurement of thermal drift.Regarding the application of Taguchi design in extrusion process,Hsia30combined Taguchi design and FE simulation to determine the better conditions in the micro-extrusion of micro pins, and proved the application potential of this method in microforming process.

    In this study,the investigated extrusion process parameters,including billet,die and container temperatures and ram speed,were selected as the control factors shown in Table 2, and a L16(44) Taguchi design with four factors and four levels was adopted.

    The standard deviation of velocity (SDV) for the PEV was employed to characterize the metal flow uniformity.To reduce unnecessary calculations,the SDV value was calculated by 1/8 of the profile cross section (i.e., the area of interest (AOI) in Fig.1), which was divided into 20 regions.The mathematical expression of SDV is

    where νmis the PEV of the center nodal at the selected region m (mm/s), ν- is the average PEV of all the divided regions (mm/s), and M is the number of the divided regions.Large SDV value represents more uneven metal flow.The SDV values for the 16 sets of extrusion parameters in Table 3 were calculated from the results of 16 steady-state extrusion simulations.

    In Taguchi design,the S/N mainly includes three characteristics, namely larger-is-better, smaller-is-better, and nominalis-best.The SDV value is expected to be minimized, thus the S/N of the smaller-is-better characteristic was selected for the SDV value.For K experimental repetitions with a given sample, the S/N of SDV value is computed according to the following Equation:

    where k is the sequence number of repeated trails, K is the number of repeated trials, and K =1 here.The unit of S/N is dB.The larger the S/N value,the smaller the SDV value,which means more uniform metal flow during extrusion.

    Based on the data of the S/N in Taguchi design, mean values of the S/N for each level and ANOVA for each factor were calculated, and the calculation procedure is detailed in the Appendix A.

    Table 2 Summary of the selected control factors and their levels.

    Table 3 L16(44)Taguchi design with corresponding values of SDV and S/N.

    2.4.Profile extrusion and characterization

    The hot extrusion experiments were performed on a 6.3 MN horizontal industrial extruder (XJ-690UST, Yuanchang Co.,China) with a container diameter of ?90 mm.Before extrusion, the billets were homogenized through holding for 18 hours at 530°C.According to Scheme 2 in Fig.3,the porthole die in Fig.5 was manufactured by the computer numerical control (CNC) machine and electron discharge machine(EDM)method.As depicted in Fig.5,the profile was extruded out of the die as the ram loads along the ED.The profile was cooled by online forced air cooling without straightening and heat treatment.Herein, the billet temperatures of 420 °C,450 °C and 480 °C were selected with other extrusion parameters of die temperature of 460 °C, container temperature of 420 °C and ram speed of 0.7 mm/s.The reasons for only verifying the billet temperature are as follows.Regarding ram speed, Liu et al.13validated that increasing ram speed aggravated the uneven metal flow and degraded the forming integrity of micro grooves during hot extrusion of a MHP profile.In terms of extrusion temperature,the simulated results in this study showed that the effects of billet, die and container temperatures on the metal flow behavior at the profile exit were similar.

    After extrusion, the head profile shape, dimension accuracy,surface micromorphology and microstructure of the profile were characterized to verify the simulated results.Through macroscopic observation of the head profile shape, the metal flow at different positions of the head profile was analyzed when the profile was just extruded out of the die.Subsequently, the forming integrity and dimension accuracy of the profile were measured by a video measuring system (VMS-2515F,Wanhao Co.,China).After that,the surface micromorphology of the micro rib was obtained using a scanning electron microscope (SEM) (Quanta 200FEG, FEI Co., USA).Finally, the electron backscatter diffraction (EBSD) analysis of the cross-section of the homogenized billet and the metallographic observation of the profile extruded at the billet temperature of 480 °C were performed, and the preparation of the electropolished EBSD samples was consistent with the study by Liu et al.13.The metallographic observations were concentrated near the area of interest in Fig.1.Before observations,the polished samples were anodized by electrolytic etching at 20 V in a mixed solution consisting of 5 mL HBF4and 200 mL deionized water.Subsequently, a polarized optical microscope(Axio Imager A2m,Zeiss Co.,Germany)was used to obtain the metallographic images.Notably, the reason for the microstructure characterization is that severe uneven metal flow behavior generally affects the microstructure,so it can not only indirectly verify the simulations, but also serve as a preexploration for the microstructure evolution studies during hot extrusion of the MHP profile.

    3.Results

    3.1.Metal flow hysteresis behavior

    The velocity distribution in Fig.6 shows that the uneven metal flow and high PEV at the profile exit are the two main characteristics during the hot extrusion.We refer to the uneven metal flow phenomenon as metal flow hysteresis behavior, which results from the complex micro-grooved features of the profile.Taking Scheme 1 in Fig.6(a) as an example, the metal flow hysteresis behavior can be categorized into two aspects.On the one hand, the metal flow velocity within the micro rib was considerably lower than that in the base region.On the other hand, significant differences were found in the PEV within the micro rib.Specifically, combining Fig.1 and Fig.6 (a), a large velocity gradient from Area 18 to Area 20 was observed in the area of interest.The closer to the top of the micro rib,the smaller the PEV.As for the high PEV caused by ultra-large extrusion ratio shown in Fig.6, the maximum PEV during extrusion was as high as 154 mm/s, whereas the ram speed was only 0.7 mm/s.An ultra-large extrusion ratio means that the cross-sectional area of the die orifice is much smaller when the billet diameter is constant.However,the volume of the metal flowing out of the die per unit time at a certain ram speed is constant.Thus, the PEV is much greater at an ultra-large extrusion ratio.

    Fig.5 6.3 MN industrial extruder and porthole die.

    Fig.6 Metal flow behavior at the profile exit during hot extrusion of the profile extruded by different dies.

    To deeply understand the metal flow hysteresis behavior,the strain distribution was analyzed.Overall, as can be seen from Fig.6(a), the strain at the thick-walled area (e.g., Point P2)where the metal with a faster flow velocity was easy to flow was the smallest,whereas the strain at the micro rib area(e.g.,Point P1) where the metal with a slower flow velocity was difficult to flow was the largest.The high strain area in Fig.6 indicates that the most severe plastic deformation occurs at the top of the micro rib.It can be determined that the material within the micro ribs undergoes more serious plastic deformation during hot extrusion.Unlike the PEV, the strain distribution at the profile exit was not only related to the shape and wall thickness of the profile,but also to the position at the profile exit.Under the same wall thickness and shape, the strain at the upper left corner of the profile was higher than that at the lower right corner of its axis symmetry owing to the asymmetry of the metal flow in the weld chamber.

    Incidentally, the effect of the optimal bearing design was also investigated.Fig.6 shows that after locally decreasing the bearing length for the micro rib in Scheme 2, the area of low PEV and high strain within the micro rib declined greatly.These results demonstrated that the bearing design of unequal length for micro features can indeed balance the metal flow at the profile exit.

    3.2.Taguchi design result and ANOVA

    The L16(44) Taguchi design and corresponding calculation results are presented in Table 3.

    To intuitively determine the optimal combination of the extrusion process parameters, the mean S/N curves were plotted in Fig.7 for various factors and levels.It was evidently that the optimal parameter combination was A4B3C3D1.Fig.7 indicates that the smaller the ram speed, the greater the S/N value, implying more uniform metal flow at the profile exit during hot extrusion of the MHP profile.Interestingly, unlike ram speed,extrusion temperature had little effect on the metal flow uniformity during extrusion.

    Based on the ANOVA results for the S/N in Table 4, the significance and contribution rate of each factor to the metal flow uniformity at the profile exit were further determined to quantify the influence magnitude of each extrusion process parameter.Regarding ANOVA, a P value less than 0.01 indicates that the association between factor and test index(i.e., S/N) is highly significant, and a P value greater than 0.1 is insignificant.Table 4 shows that the P value of ram speed was less than 0.005,thus the impact of ram speed on the metal flow uniformity during hot extrusion is extremely significant.However, the P values of extrusion temperature (e.g., billet temperature) are greater than 0.10, which confirms their insignificant correlation to the metal flow uniformity during hot extrusion.

    Fig.7 The average S/N for four factors at various levels.

    The contribution rate ηiof each factor can be determined by

    where i is the sequence number of the ith factor in Table 2,niis the number of factors, and Fiis Fisher’s value of the ith factor.During hot extrusion process of the profile, the effect of ram speed on the metal flow uniformity at the profile exit is dominant with a contribution of as high as 99.64 %.However, the influences of extrusion temperature are almost negligible,because the contribution rates of billet,die and container temperatures are only 0.17 %, 0.11 % and 0.09 %, respectively.In other words, extrusion temperature does not need to be taken into account when optimizing the metal flow uniformity at the profile exit during hot extrusion of MHPs.

    3.3.Head profile shape

    The metal flow hysteresis behavior found in the simulation was confirmed through macroscopic observation of the head profile.Fig.8(a) indicates that the metal flow within the micro rib lagged far behind that within the base region, which verified the simulated results in Fig.6.In other words, the PEV within the micro rib region is considerably smaller than that within the base region.Additionally, consistent with the simulated results, there are also significant differences in the metal flow within the micro rib.For example, the metal flow at the top of the micro rib lagged seriously behind that at the bottom.Furthermore, the shape of the head profile extruded ranging from 420 to 480°C in Fig.8 was used to evaluate the influence of billet temperature on the metal flow hysteresis behavior.Through macroscopic observation, no significant difference in the shape of the head profiles extruded at various billet temperatures can be identified.Therefore,the experimental results also show that the billet temperature has no significant effect on the metal flow uniformity during hot extrusion of the profile.

    3.4.Dimension accuracy

    The forming integrity and dimension accuracy of the profile were further investigated.Fig.9 shows that the forming integ-rity of the micro grooves was very good at various billet temperatures, and both the complex micro-grooved feature and the outer fillet of the profile were completely extruded.However,the No.6 micro groove(highlighted by green dashed line)extruded at the billet temperature of 420°C was slightly wider.

    Table 4 ANOVA results for S/N.

    Fig.8 Macroscopic observation of the head profiles extruded at various billet temperatures.

    Fig.9 Cross-sectional images of the profiles extruded at various billet temperatures.

    Fig.10 Dimensions of micro grooves of the profile extruded at various billet temperatures.

    To determine the influence of billet temperature on the dimension accuracy of the profile, the minimum width dimension (i.e., 0.27 mm) of all micro grooves was measured and counted in Fig.10.Under the condition of the billet temperatures of 480 °C and 450 °C, the dimension accuracy of all micro grooves was as high as ± 0.02 mm.Yet, the dimension of the micro groove No.6 extruded at the billet temperature of 420 °C was 0.295 mm in width, and its dimension accuracy decreased slightly.Nevertheless, the dimension accuracy of the micro-grooved width was still better than ± 0.03 mm,which was a good indicator.The slight decline in dimension accuracy of the profile extruded at the billet temperature of 420°C may be due to an increase in uneven metal flow caused by part feature size effect, but also a decrease in metal flow performance.

    3.5.Surface micromorphology

    The SEM images in Fig.11 show the top surface of the micro ribs on both sides of the micro groove No.5 of the profile.

    The low-magnification images show that, consistent with the results in Fig.10, the profiles extruded at various billet temperatures had good consistency in the minimum width of the micro grooves.Under all billet temperature conditions,the surface quality of the micro ribs was highly excellent,showing the characteristics of extrusion streaks parallel to the ED.It is generally believed that in the hot extrusion of aluminum alloys without lubrication, extrusion streaks form due to a metal film adhering to the bearing that has a brushing effect on the surface of the extrudates flowing out of the extrusion die.No micro-cracks were observed, as indicated by the high-magnification images of the edges of the micro ribs near the micro groove No.6.Incidentally, the number and size of the extrusion streaks at the edge of the micro ribs were larger at the higher billet temperature of 480°C.However,it was not statistically significant due to the small observation area.

    Fig.11 SEM images of the profiles extruded at various billet temperatures.

    3.6.Microstructure

    Fig.12 shows the microstructure of the homogenized as-cast billet and extruded profile.According to the literature by Sla′-mova′ et al.31, the metallographic image obtained by anodic coating method can not only show the grain morphology,but also reflect the grain orientation to a certain extent.This is due to the formation of Al2O3film on the sample surface.During the anodic coating of polycrystalline materials, the thickness of the films deposited on grains with different orientations is different.Consequently,they will show different contrast when observed under a polarizing microscope, and the color difference between grains represents the difference of crystal orientation.

    As expected,the metal flow hysteresis behavior had a significant impact on the microstructure within the micro rib.Specifically, under the polarizing microscope, the grain color within the micro rib (e.g., Region 1) was completely different from that within the base region (e.g., Region 2), which indicated that the crystal orientations of the grains were distinct within these two regions.Region 1 consisted mainly of the grains of Orientation a, whereas Region 2 consisted of the grains of Orientation b.Interestingly, the difference in grain size between Region 1 and Region 2 was insignificant.

    4.Discussion

    4.1.Simulation credibility analysis

    In this study, the conventional material constitutive obtained by the non-microforming tests was used to simulate the hot extrusion of the MHP profile.Therefore, necessary simulation credibility analysis was performed.Consistent with the study by Liu et al.13, the metal flow hysteresis behavior during hot extrusion of the MHP profile was found in the simulated results in Fig.6 and well verified by the experimental results in Fig.8.In addition, Fig.12(b) shows that the grain orientation within the micro rib was completely different from that within the base region.However, according to Fig.12(a), the as-cast billet had no obvious preferred crystal orientation,which indirectly proved the existence of the metal flow hysteresis behavior.Therefore, it can be determined that the steadystate extrusion simulation using the conventional material constitutive model is reliable in qualitatively predicting the metal flow behavior during hot extrusion of macro parts with micro features such as grooved MHP profiles.As for the reason, it may be due to the process characteristics of hot extrusion of MHP profiles.Specifically, in microforming, the critical number of grains along a certain direction that falls into the size effect occurrence range depends on the type of materials, such as ‘‘less than 10 grains”mentioned by Hansen32and Kim et al.33.However, unlike cold micro-extrusion, although the grain size of as-cast billets generally used in hot extrusion is large(e.g.,100 μm in this study shown in Fig.12(a)),it may be refined due to dynamic recrystallization when the material reaches the bearing after hot extrusion deformation.For example, Fig.12(b) shows that the number of grains in the width direction of the micro rib can even be up to 10 at most,which indicates that sufficient dynamic recrystallization occurred.It is evident from the reverse inference that considering the grain growth during hot extrusion process,the number of grains in the width direction of the micro rib is likely to be more when the material is about to flow into the bearing.To sum up, this is the reason why the conventional material constitutive model works well in the hot extrusion of micro grooves of the MHP profile.

    Fig.12 Microstructure of billet and extruded profile.

    4.2.Formation mechanism of the metal flow hysteresis behavior

    The formation mechanism of the metal flow hysteresis behavior resulting from size effect during hot extrusion of the MHP profile is discussed based on the simulated results.According to Chan et al.17,the size effect on material deformation behaviors is manifested as grain size,workpiece size,part feature size and interfacial condition.To be precise,in the extrusion of the MHP profile,only the deformation of the micro rib belongs to micro-extrusion,and that of the other parts belong to conventional extrusion.Thus, this metal flow hysteresis behavior is due to the part feature size effect.

    In this study, the friction-affected area was proposed to elaborate the metal flow hysteresis behavior caused by part feature size effect.As shown in Fig.6, the presence of the friction-affected area can be inferred from the velocity distribution near the surface in contact with the die.Specifically,it is subjected to the frictional resistance reverse to the ED when the metal contacting the die (e.g., die core and orifice)is extruded out of the die.Taking Point P3 in Fig.6 as an example, the frictional resistance can reduce the metal flow velocity and increase the deformation intensity of the layer of metal in contact with the die, which leads to the formation of a friction-affected area with a certain thickness.For example, the velocity cloud diagram in Fig.6 (a) shows that the dark blue area near Point P3 is in the friction-affected area where the flow velocity is very low.In the thickness direction of the friction-affected zone, moving away from the die weakens the effect of friction on the metal flow, resulting in an increase in the flow velocity.Moreover,Fig.6(a)indicates that the area of the friction-affected area within the micro rib is close to the order of magnitude of the area of the micro rib.Therefore, the area proportion of the friction-affected area in the micro rib (e.g., the region near Point P1) is far larger than that in the base region(e.g.,the region near Point P2),and the metal is more likely to flow out of the thick-walled area where flow resistance is lower, resulting in slower metal flow within the micro rib.

    Based on the above discussion, and inspired by so-called surface layer models used to account for grain or specimen size effects,34a physical model in Fig.13 was proposed to illustrate the formation mechanism of the metal flow behavior affected by part feature size effect during hot extrusion of the MHP profile.For the conventional extrusion and extrusion of macro parts with micro features(e.g.,MHP profiles),the schematic of the friction-affected area (pink area) within the bearing of the extrusion die was presented.In Fig.13, νEDrepresents the metal flow velocity along the ED.The area ratio of the friction-affected area to the region in which it is located,denoted asξ, is used to illustrate the differences between the extrusion of macro parts with micro features and conventional extrusion, which is expressed as follows,

    where Afis the cross-sectional area of the friction-affected area (mm2), and Amis the cross-sectional area of macro or micro features (e.g., rib or micro rib) (mm2).The larger the ξ value, the more severe the part feature size effect.It is evident that for the conventional extrusion of macro features, the ξ value is much less than 100 %, and the effects of frictionaffected area on the metal flow behavior within the bearing is usually negligible.However, the ξ value increases sharply when the size of macro features of the profile reaches the microforming scale.As depicted in Fig.13, the ξ value gradually approaches 100 % with the decrease of the size of the micro features.It is precisely due to the sharp increase of ξ value that the velocity distribution characteristics depicted in Fig.13 appear during the hot extrusion of the MHP profile,resulting in a significant hysteresis of the metal flow within micro features.Therefore, the effect of friction-affected area,which can be ignored in conventional extrusion, cannot be neglected during the extrusion of macro parts with micro features.

    Fig.13 Schematic of friction-affected area in conventional extrusion and extrusion of macro parts with micro features.

    4.3.Influence magnitude of extrusion parameter on the metal flow hysteresis behavior and forming integrity

    In order to guide the formulation of extrusion process parameters in actual production, their influence on the metal flow hysteresis behavior at the profile exit and resulting forming integrity during hot extrusion of the MHP profile are further discussed.

    Regarding ram speed, as shown in Fig.7, the metal flow uniformity at the profile exit drops seriously as the ram speed increases.Nevertheless,the friction-affected area has not changed.The reason why high ram speed reduces the metal flow uniformity during extrusion subjected to part feature size effect is as follows.According to Eq.(1), the SDV value also increases correspondingly with the increase of ram speed.For example, when νmνmis magnified by a factor of 4, the νand SDV value are evidently also magnified by a factor of 4.As shown in Table 3, when the ram speed was increased from 0.7 mm/s in Case 1 to 2.8 mm/s in Case 4, the SDV value of Case 4 (104.08) was approximately-four times that of Case 1(27.96).Therefore, in essence, the ram speed does not change the distribution of the friction-affected area, but multiplies the inhomogeneity of metal flow at the profile exit.Interestingly, the ANOVA results in Table 4 established that ram speed dominated the metal flow uniformity among all the investigated extrusion process parameters.Table 4 presents that billet temperature has insignificant effects with a contribution rate of approximately 0.17%on the metal flow uniformity at the profile exit, which was verified by the analysis of the head profile shape in Fig.8 and the forming quality in Figs.9–11.On the one hand, Fig.10 shows that the dimension accuracy of the micro-grooved width decreased slightly from ± 0.02 mm to ± 0.03 mm as the billet temperature decreased from 480 to 420 °C.On the other hand, according to Liu et al.13, the increased inhomogeneity of metal flow at the profile exit could deteriorate the surface quality of micro ribs due to the promotion of the initiation and propagation of micro-cracks and even the tearing of the micro ribs.However, as shown in Fig.11, the surface micromorphology of the micro ribs extruded at different billet temperatures was not significantly different,and no micro-cracks were observed.As a result of the discussion, increasing the ram speed multiplies the uneven metal flow and severely reduces the forming integrity of the MHP profile, whereas extrusion temperature has little effect.Nevertheless, neither ram speed nor extrusion temperature can change the distribution of the friction-affected area caused by part feature size effect.

    Fig.14 Successfully extruded MHP profile at Tb = 480 °C.

    Based on the above studies, a strategy for the formulation of extrusion process parameters with adopting a low ram speed and ignoring extrusion temperature can be determined to attenuate the effect of the metal flow hysteresis behavior on the forming integrity during hot extrusion of the MHP profile.In this study,under the combination of unequal bearing design and better extrusion process parameters (i.e., ram speed of 0.7 mm/s, billet temperature of 480 °C, die temperature of 460°C,and container temperature of 420°C),the MHP profile with a micro-grooved width of 0.27 ± 0.02 mm was successfully extruded in large quantities on an industrial extruder.Fig.14 shows the successfully extruded MHP profile that was cut to length.

    5.Conclusions

    (1) During the porthole die extrusion of a micro-grooved MHP profile using coarse-grained(~100 μm)as-cast billets, the conventional material constitutive model obtained by non-microforming tests is reliable in qualitatively predicting the metal flow behavior at the profile exit of the die.

    (2) The metal flow within micro ribs lags far behind that within the base region during the MHP profile extrusion.A physical model was proposed to reveal the metal flow hysteresis behavior caused by part feature size effect.According to this model, the friction-affected area that can be negligible during conventional extrusion seriously slowed the metal flow within micro features (e.g., micro ribs)during the extrusion of macro parts with micro features.The area ratio of the friction-affected area to the region in which is located is denote as ξ.For the extrusion of macro features, the ξ value is much less than 100 %, whereas the ξ value gradually approaches 100%with the decrease of the size of the micro features.

    (3) Neither ram speed nor extrusion temperature can change the distribution of the friction-affected area.However,increasing ram speed multiplies the metal flow hysteresis behavior and severely reduces the forming integrity of MHP profiles, whereas extrusion temperature has little effect.In the actual production of MHP profiles, the ram speed should be as small as possible to prevent multiplying the uneven metal flow, and the effect of extrusion temperature can be ignored.

    (4) The metal flow hysteresis behavior leads to completely distinct grain orientations between the micro rib and base region of the MHP profile.Future research should focus on the microstructure and mechanical properties of the MHP profile.

    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.

    Acknowledgements

    This study was co-supported by the National Natural Science Foundation of China(No.51635005)and the 111 Project(No.B18017).

    Appendix A.A.1.Calculation of mean S/N for each level

    The mean S/N for each level of the control factors can be calculated by

    where S/Nijis the S/N value of the ith factor at the jth level(dB), i is the sequence number of the ith factor in Table 2, j is the sequence number of the jth level in Table 2, and njis the number of levels.The corresponding calculation results are indicated in Table A1,which is used to plot Fig.7 in the paper.

    A.2.Calculation of ANOVA

    ANOVA mainly includes the calculation of sum of deviation squares (referred to as sum of squares), degree of freedom,mean square and Fisher’s value (F value).All calculation results are shown in Table 4.

    The sum of total squares (ST) of S/N can be calculated by

    Table A1 Mean S/N for each level.

    Table A2 kij values for each level.

    where Fiis Fisher’s value of the ith factor (see Tables A1 and A2).

    美女国产高潮福利片在线看| 国产成人免费观看mmmm| 国产精品久久久人人做人人爽| 大型黄色视频在线免费观看| 曰老女人黄片| 老熟妇乱子伦视频在线观看| 99国产精品一区二区三区| 免费观看精品视频网站| 久久久久精品国产欧美久久久| 一进一出抽搐动态| 欧美成狂野欧美在线观看| 午夜福利在线观看吧| 亚洲av第一区精品v没综合| 日韩视频一区二区在线观看| 十八禁网站免费在线| 欧美激情 高清一区二区三区| 精品久久久久久,| 亚洲成人免费av在线播放| 欧美 日韩 精品 国产| 亚洲三区欧美一区| 国产精品一区二区精品视频观看| netflix在线观看网站| 一级毛片高清免费大全| 国产成人精品久久二区二区免费| 国精品久久久久久国模美| 人人妻人人添人人爽欧美一区卜| 黄色a级毛片大全视频| 黄色毛片三级朝国网站| 日日夜夜操网爽| 国产99久久九九免费精品| 啦啦啦视频在线资源免费观看| xxxhd国产人妻xxx| 久久人妻av系列| 欧美大码av| 精品免费久久久久久久清纯 | 三级毛片av免费| 久久久国产一区二区| 亚洲精品美女久久久久99蜜臀| av中文乱码字幕在线| 在线观看一区二区三区激情| tube8黄色片| 少妇被粗大的猛进出69影院| 亚洲久久久国产精品| 黑丝袜美女国产一区| 亚洲av第一区精品v没综合| 男女下面插进去视频免费观看| 亚洲五月天丁香| 亚洲人成电影观看| 18禁裸乳无遮挡免费网站照片 | 老司机影院毛片| 精品久久久久久久毛片微露脸| 亚洲人成77777在线视频| 男女免费视频国产| 高清毛片免费观看视频网站 | 91av网站免费观看| 久久99一区二区三区| 国产成人精品无人区| 久久天躁狠狠躁夜夜2o2o| 国产精华一区二区三区| 久9热在线精品视频| 18禁黄网站禁片午夜丰满| 女警被强在线播放| 在线av久久热| √禁漫天堂资源中文www| 久久影院123| 中文字幕高清在线视频| 中文字幕高清在线视频| 搡老熟女国产l中国老女人| 新久久久久国产一级毛片| 精品人妻1区二区| 亚洲精品久久成人aⅴ小说| 丰满的人妻完整版| 99在线人妻在线中文字幕 | 91国产中文字幕| 国产精品成人在线| 国产熟女午夜一区二区三区| 99精品久久久久人妻精品| 亚洲欧美激情综合另类| 欧美一级毛片孕妇| svipshipincom国产片| 久久人人爽av亚洲精品天堂| 久久青草综合色| 成人永久免费在线观看视频| 少妇 在线观看| 国产精品99久久99久久久不卡| 国产亚洲精品久久久久久毛片 | 黄片大片在线免费观看| 99精品久久久久人妻精品| 欧美成人免费av一区二区三区 | 涩涩av久久男人的天堂| av视频免费观看在线观看| 国产一区二区三区综合在线观看| 男女之事视频高清在线观看| 69av精品久久久久久| 午夜福利在线观看吧| 超色免费av| 欧美在线一区亚洲| 国产精品成人在线| 国产97色在线日韩免费| 夫妻午夜视频| 一级毛片女人18水好多| 色播在线永久视频| 黑人操中国人逼视频| 久久中文看片网| 国产aⅴ精品一区二区三区波| 国产成人av激情在线播放| 精品国产国语对白av| 国产乱人伦免费视频| 亚洲国产欧美日韩在线播放| 91字幕亚洲| 在线观看一区二区三区激情| 国产精品免费一区二区三区在线 | 男女免费视频国产| 18禁裸乳无遮挡免费网站照片 | 99热国产这里只有精品6| 99热网站在线观看| 亚洲精品乱久久久久久| 精品一区二区三区av网在线观看| 久久人人爽av亚洲精品天堂| 日日摸夜夜添夜夜添小说| 欧美日韩瑟瑟在线播放| 国产精品98久久久久久宅男小说| 中文字幕制服av| 一级a爱片免费观看的视频| 国产精品久久久人人做人人爽| 老汉色av国产亚洲站长工具| 久久中文字幕一级| 欧美精品一区二区免费开放| 成年女人毛片免费观看观看9 | avwww免费| 亚洲精品美女久久久久99蜜臀| 黄色视频,在线免费观看| 一级片免费观看大全| 老司机午夜十八禁免费视频| 国产精品欧美亚洲77777| 人人澡人人妻人| 欧美日韩国产mv在线观看视频| 久久国产精品影院| 男人舔女人的私密视频| 国产淫语在线视频| 亚洲中文av在线| 一夜夜www| 久久性视频一级片| 欧美黄色淫秽网站| 亚洲在线自拍视频| 欧美性长视频在线观看| 亚洲成人免费电影在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲熟妇中文字幕五十中出 | 日韩欧美免费精品| 麻豆av在线久日| 国产欧美日韩综合在线一区二区| 国产单亲对白刺激| 黄色毛片三级朝国网站| 手机成人av网站| 十八禁高潮呻吟视频| 欧美在线黄色| 女性生殖器流出的白浆| 18在线观看网站| 欧美丝袜亚洲另类 | x7x7x7水蜜桃| 成人三级做爰电影| 99精品在免费线老司机午夜| 精品久久久精品久久久| 一边摸一边做爽爽视频免费| 午夜影院日韩av| 19禁男女啪啪无遮挡网站| 免费女性裸体啪啪无遮挡网站| 香蕉国产在线看| 欧美乱妇无乱码| 18禁黄网站禁片午夜丰满| av一本久久久久| 王馨瑶露胸无遮挡在线观看| 亚洲欧美精品综合一区二区三区| 亚洲专区中文字幕在线| 久久精品国产清高在天天线| 国产男女超爽视频在线观看| 国产精华一区二区三区| 欧美一级毛片孕妇| 午夜免费鲁丝| 一二三四在线观看免费中文在| 久久久久久人人人人人| 国产主播在线观看一区二区| 一个人免费在线观看的高清视频| 两性夫妻黄色片| 欧洲精品卡2卡3卡4卡5卡区| 超碰成人久久| 乱人伦中国视频| 国产欧美亚洲国产| 男女之事视频高清在线观看| 天天操日日干夜夜撸| 少妇的丰满在线观看| 久久 成人 亚洲| 妹子高潮喷水视频| 18禁国产床啪视频网站| 亚洲五月天丁香| 国产精品亚洲av一区麻豆| 午夜日韩欧美国产| 丰满饥渴人妻一区二区三| 亚洲第一av免费看| 性色av乱码一区二区三区2| 女警被强在线播放| 免费在线观看日本一区| 欧美成人午夜精品| 亚洲片人在线观看| 天堂中文最新版在线下载| 女人精品久久久久毛片| 亚洲国产欧美网| 嫩草影视91久久| 欧美日本中文国产一区发布| 国产av又大| 久久人妻av系列| 国产亚洲av高清不卡| 在线观看一区二区三区激情| 久久中文字幕一级| 两人在一起打扑克的视频| 国产1区2区3区精品| 黄网站色视频无遮挡免费观看| 国产精品国产av在线观看| av网站免费在线观看视频| 少妇的丰满在线观看| 国产免费男女视频| 午夜免费成人在线视频| 亚洲熟妇中文字幕五十中出 | 欧美精品一区二区免费开放| 国产日韩一区二区三区精品不卡| 在线观看免费视频日本深夜| 久久狼人影院| 久久人妻熟女aⅴ| 99久久精品国产亚洲精品| 亚洲成人手机| 午夜视频精品福利| 99国产精品免费福利视频| 欧美成狂野欧美在线观看| 久久ye,这里只有精品| 亚洲精品一卡2卡三卡4卡5卡| 国产在线观看jvid| 丝袜美腿诱惑在线| 精品久久久久久电影网| 久久亚洲真实| 建设人人有责人人尽责人人享有的| 老熟妇仑乱视频hdxx| 日韩 欧美 亚洲 中文字幕| 欧美在线一区亚洲| 欧美日韩成人在线一区二区| 淫妇啪啪啪对白视频| 欧美人与性动交α欧美软件| 久久精品国产亚洲av香蕉五月 | 中国美女看黄片| 三级毛片av免费| 一区二区三区国产精品乱码| 日日爽夜夜爽网站| 97人妻天天添夜夜摸| 亚洲美女黄片视频| 天堂动漫精品| 亚洲九九香蕉| x7x7x7水蜜桃| 两个人看的免费小视频| 色播在线永久视频| 欧美另类亚洲清纯唯美| 久久人人爽av亚洲精品天堂| 国精品久久久久久国模美| 悠悠久久av| 一进一出好大好爽视频| 免费看十八禁软件| 美女 人体艺术 gogo| 成人影院久久| 亚洲成人免费电影在线观看| 免费看a级黄色片| 久久精品国产清高在天天线| 国产精品偷伦视频观看了| www.精华液| 欧美乱色亚洲激情| 999久久久精品免费观看国产| 成人av一区二区三区在线看| 天堂√8在线中文| 亚洲国产精品一区二区三区在线| 麻豆av在线久日| 亚洲情色 制服丝袜| 国产精品国产高清国产av | 飞空精品影院首页| 国产精品亚洲一级av第二区| 亚洲男人天堂网一区| 久久精品亚洲av国产电影网| 人人澡人人妻人| 丝袜在线中文字幕| 男女床上黄色一级片免费看| 日本一区二区免费在线视频| 丁香六月欧美| 免费少妇av软件| 日韩人妻精品一区2区三区| 最新美女视频免费是黄的| 亚洲精品中文字幕在线视频| 天天添夜夜摸| 成人国语在线视频| 国产视频一区二区在线看| 热99久久久久精品小说推荐| 久9热在线精品视频| 黄片播放在线免费| 日韩免费av在线播放| 久久中文字幕一级| 久久久久精品国产欧美久久久| 一边摸一边做爽爽视频免费| 久久精品人人爽人人爽视色| 十分钟在线观看高清视频www| 欧美国产精品va在线观看不卡| 在线视频色国产色| 69av精品久久久久久| 成人18禁高潮啪啪吃奶动态图| 久久久精品区二区三区| 在线免费观看的www视频| 色精品久久人妻99蜜桃| 高潮久久久久久久久久久不卡| 国产欧美日韩一区二区三区在线| 捣出白浆h1v1| 国产精品99久久99久久久不卡| 亚洲全国av大片| 人成视频在线观看免费观看| 国产亚洲欧美98| 国产亚洲欧美98| 水蜜桃什么品种好| 久久国产精品男人的天堂亚洲| 精品人妻在线不人妻| 18禁国产床啪视频网站| 露出奶头的视频| 黄片大片在线免费观看| 在线播放国产精品三级| 久久精品亚洲精品国产色婷小说| 国产精品98久久久久久宅男小说| 91字幕亚洲| 久久香蕉国产精品| av网站免费在线观看视频| 国产国语露脸激情在线看| av电影中文网址| 天天躁夜夜躁狠狠躁躁| 中文字幕人妻丝袜制服| 午夜精品在线福利| 国产精品 欧美亚洲| 又紧又爽又黄一区二区| 91成人精品电影| 女性生殖器流出的白浆| 极品人妻少妇av视频| 午夜免费鲁丝| 女同久久另类99精品国产91| 精品人妻熟女毛片av久久网站| 亚洲黑人精品在线| 女人被狂操c到高潮| 日韩免费高清中文字幕av| 老熟女久久久| 最近最新中文字幕大全免费视频| 在线观看免费午夜福利视频| av天堂久久9| 亚洲av成人不卡在线观看播放网| 欧美日本中文国产一区发布| 亚洲精品乱久久久久久| 精品久久久精品久久久| 色在线成人网| 超色免费av| 多毛熟女@视频| 欧美 日韩 精品 国产| 日韩欧美三级三区| 成人黄色视频免费在线看| 欧美乱码精品一区二区三区| 成人国产一区最新在线观看| 国产91精品成人一区二区三区| 69精品国产乱码久久久| 欧美大码av| 久久精品亚洲精品国产色婷小说| 婷婷成人精品国产| 欧美黑人欧美精品刺激| 国产熟女午夜一区二区三区| 91老司机精品| 欧美精品啪啪一区二区三区| 欧美老熟妇乱子伦牲交| 免费观看人在逋| 精品第一国产精品| 一边摸一边抽搐一进一出视频| 纯流量卡能插随身wifi吗| 人人妻人人澡人人爽人人夜夜| tocl精华| 人妻丰满熟妇av一区二区三区 | 一级片免费观看大全| 久久久久久久久久久久大奶| 国产99久久九九免费精品| 欧美日韩黄片免| 亚洲av日韩精品久久久久久密| 精品久久久久久电影网| 侵犯人妻中文字幕一二三四区| 波多野结衣av一区二区av| 51午夜福利影视在线观看| 午夜精品久久久久久毛片777| svipshipincom国产片| 下体分泌物呈黄色| 日韩欧美一区视频在线观看| 69精品国产乱码久久久| 91大片在线观看| www.自偷自拍.com| 日韩视频一区二区在线观看| 成熟少妇高潮喷水视频| 美女国产高潮福利片在线看| 99在线人妻在线中文字幕 | 中文字幕av电影在线播放| 天天躁夜夜躁狠狠躁躁| 性色av乱码一区二区三区2| 亚洲综合色网址| 一本综合久久免费| 欧美午夜高清在线| 一区二区三区国产精品乱码| 一级a爱视频在线免费观看| 美女高潮喷水抽搐中文字幕| 桃红色精品国产亚洲av| 王馨瑶露胸无遮挡在线观看| 亚洲精品在线观看二区| 人人妻人人澡人人爽人人夜夜| av一本久久久久| 天堂俺去俺来也www色官网| 国产伦人伦偷精品视频| 99热网站在线观看| 一本大道久久a久久精品| 国产又色又爽无遮挡免费看| 欧美亚洲 丝袜 人妻 在线| 国产亚洲精品一区二区www | 伊人久久大香线蕉亚洲五| 母亲3免费完整高清在线观看| 亚洲在线自拍视频| 免费黄频网站在线观看国产| 欧美激情高清一区二区三区| 女人精品久久久久毛片| 在线十欧美十亚洲十日本专区| 国产精品国产高清国产av | 天天影视国产精品| 一二三四社区在线视频社区8| 午夜精品久久久久久毛片777| 亚洲成人手机| 久久天堂一区二区三区四区| 欧美亚洲 丝袜 人妻 在线| 一本大道久久a久久精品| 一a级毛片在线观看| 大型av网站在线播放| 一级毛片女人18水好多| 日韩制服丝袜自拍偷拍| 亚洲五月天丁香| 午夜免费鲁丝| 国产精品香港三级国产av潘金莲| 日韩人妻精品一区2区三区| 精品午夜福利视频在线观看一区| 成人av一区二区三区在线看| 黄色视频,在线免费观看| 怎么达到女性高潮| 嫩草影视91久久| 丝袜人妻中文字幕| 亚洲一区高清亚洲精品| a级毛片黄视频| 国产精品九九99| 大码成人一级视频| av线在线观看网站| 男女免费视频国产| 国精品久久久久久国模美| 丝瓜视频免费看黄片| 国产男女超爽视频在线观看| 亚洲人成伊人成综合网2020| 国产激情久久老熟女| 国产aⅴ精品一区二区三区波| 中文亚洲av片在线观看爽 | 亚洲精品美女久久av网站| 欧美精品av麻豆av| 波多野结衣一区麻豆| 精品国产国语对白av| 久久久水蜜桃国产精品网| 国产精品成人在线| 亚洲中文日韩欧美视频| 韩国精品一区二区三区| 90打野战视频偷拍视频| 精品一区二区三区四区五区乱码| 黑人操中国人逼视频| 国产又色又爽无遮挡免费看| 嫁个100分男人电影在线观看| 91麻豆av在线| 美女视频免费永久观看网站| 国产午夜精品久久久久久| 一级作爱视频免费观看| 午夜福利影视在线免费观看| 黄色女人牲交| 老司机福利观看| 欧美日韩黄片免| 啪啪无遮挡十八禁网站| 国产成人av激情在线播放| 黄色女人牲交| 18禁黄网站禁片午夜丰满| 美女视频免费永久观看网站| 精品免费久久久久久久清纯 | 午夜激情av网站| 国产成人系列免费观看| 亚洲黑人精品在线| 怎么达到女性高潮| 欧美黄色片欧美黄色片| 久久久久视频综合| 国产在线观看jvid| 18禁裸乳无遮挡免费网站照片 | 欧美日韩乱码在线| 人人澡人人妻人| 精品亚洲成a人片在线观看| 女警被强在线播放| 国产精品98久久久久久宅男小说| 老汉色av国产亚洲站长工具| xxx96com| 久久精品国产a三级三级三级| av视频免费观看在线观看| 久久国产精品人妻蜜桃| 欧美最黄视频在线播放免费 | 午夜日韩欧美国产| 国产精品乱码一区二三区的特点 | 天天躁日日躁夜夜躁夜夜| 窝窝影院91人妻| 757午夜福利合集在线观看| 人妻 亚洲 视频| 亚洲一区二区三区欧美精品| 欧美日韩国产mv在线观看视频| 欧美老熟妇乱子伦牲交| 免费av中文字幕在线| 黑人巨大精品欧美一区二区蜜桃| 两性夫妻黄色片| 欧美乱码精品一区二区三区| 国产高清激情床上av| 亚洲精品乱久久久久久| 国产无遮挡羞羞视频在线观看| 黄片大片在线免费观看| 欧美乱码精品一区二区三区| 80岁老熟妇乱子伦牲交| 欧美一级毛片孕妇| 久久久久久人人人人人| 超色免费av| 免费在线观看视频国产中文字幕亚洲| 欧美日韩一级在线毛片| 欧美另类亚洲清纯唯美| 中文亚洲av片在线观看爽 | 91麻豆精品激情在线观看国产 | 亚洲国产欧美网| 深夜精品福利| √禁漫天堂资源中文www| 悠悠久久av| 制服诱惑二区| 美女福利国产在线| 在线视频色国产色| 国产成人一区二区三区免费视频网站| 777久久人妻少妇嫩草av网站| 欧美在线一区亚洲| 黄色毛片三级朝国网站| 国产亚洲欧美在线一区二区| 亚洲国产中文字幕在线视频| 国产亚洲一区二区精品| 成人手机av| 曰老女人黄片| 99国产精品一区二区蜜桃av | 日韩免费高清中文字幕av| 黄色视频,在线免费观看| 在线观看舔阴道视频| 色在线成人网| 亚洲人成电影免费在线| 久热爱精品视频在线9| 久久亚洲精品不卡| 免费在线观看亚洲国产| 一二三四在线观看免费中文在| 亚洲国产欧美一区二区综合| 亚洲成人手机| 国产精品自产拍在线观看55亚洲 | 中文字幕高清在线视频| 亚洲熟妇熟女久久| 高清在线国产一区| 宅男免费午夜| 在线av久久热| 亚洲国产毛片av蜜桃av| 午夜精品久久久久久毛片777| 啦啦啦在线免费观看视频4| 18禁国产床啪视频网站| 亚洲精品乱久久久久久| 成年版毛片免费区| 欧美日韩亚洲综合一区二区三区_| 国产精品欧美亚洲77777| av视频免费观看在线观看| 国产区一区二久久| 国产精品久久久av美女十八| av福利片在线| 777久久人妻少妇嫩草av网站| 亚洲一区二区三区欧美精品| 91精品国产国语对白视频| 91成年电影在线观看| 黑人欧美特级aaaaaa片| 亚洲精品美女久久av网站| 欧美成狂野欧美在线观看| 777久久人妻少妇嫩草av网站| 精品国产乱码久久久久久男人| 一级毛片精品| 国产一卡二卡三卡精品| 久久九九热精品免费| 免费日韩欧美在线观看| 男女下面插进去视频免费观看| 美女扒开内裤让男人捅视频| 丰满的人妻完整版| 精品久久久久久电影网| 免费少妇av软件| 精品久久久久久电影网| 村上凉子中文字幕在线| 欧美日韩视频精品一区| 变态另类成人亚洲欧美熟女 | 午夜福利视频在线观看免费| 免费看十八禁软件| 在线观看舔阴道视频| 一夜夜www| 久久 成人 亚洲| 欧美日韩乱码在线| 国产在视频线精品| 久久久久久久精品吃奶| 在线国产一区二区在线| 18禁美女被吸乳视频|