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

    Hole surface strengthening mechanism and riveting fatigue life of CFRP/aluminum stacks in robotic rotary ultrasonic drilling

    2023-11-10 02:16:36SongDONGWenheLIAOKanZHENGFengXUELianjunSUN
    CHINESE JOURNAL OF AERONAUTICS 2023年10期

    Song DONG, Wenhe LIAO, Kan ZHENG, Feng XUE, Lianjun SUN

    School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China

    KEYWORDS

    Abstract Carbon fiber reinforced plastic(CFRP)and aluminum stacks are widely used in aviation industry due to light weight and high performance.Millions of rivet holes need to be drilled on body materials, and more than 80% of fatigue cracks occur at the connection holes, so the damage and residual stress of hole surface have crucial effect on the riveting fatigue life of CFRP/aluminum stacks and the flight performance.Recently, robotic rotary ultrasonic drilling (RRUD) technology is a promising method to machine the stacks.However, the hole surface strengthening mechanism in RRUD and the service performance of the riveting joint are not verified.Thus,in this paper,the hole surface strengthening mechanism of RRUD for CFRP/aluminum stacks is investigated,a theoretical residual stress model is established, and the fatigue life experiment of riveted joints is conducted.Firstly, analysis on residual stress in RRUD is carried out with consideration of strengthening force and cutting temperature.Residual stress model is established based on the calculation of elastic stress, plastic stress and stress release.Validation experiment results show that ultrasonic vibration changes residual stress from tensile stress to compressive stress.At the same time,comparative damage analysis of CFRP hole exit and hole surface in robotic conventional drilling(RCD)and RRUD is presented.Finally,fatigue strength experiments of riveted joints are conducted for performance verification.Experimental results indicate that fatigue life of single-hole riveted joints is increased by 68% with ultrasonic vibration, and four-hole riveted joint arranged according to aerospace design standards is increased by more than 86%.

    1.Introduction

    In order to meet the demands of lightweight and structural strength,CFRP/aluminum stacks have been widely used in aircraft manufacturing including the F35 Lightning II,the Airbus A380 and the Boeing 787 Dreamliner.1–6During airframe manufacturing and assembly process, millions of rivet holes need to be drilled on body stacks materials.7–9According to the statistics, 70% of aircraft body fatigue failure originates from structural connections, and more than 80% of fatigue cracks occur at the connection holes.10So the quality of hole,especially the damage and residual stress of hole surfaces, is crucially important in the riveting fatigue life of CFRP/aluminum stacks and even the flight performance.Compared with manual hole drilling, robotic drilling has the advantages of high flexibility, high quality consistency and high phase accuracy.11–136-DOF serial robot equipped with drilling end effector offers operation flexibility, which is regarded as advanced option to accomplish aircraft component drilling.14–17.

    In recent years, rotary ultrasonic drilling (RUD) technology, as a promising machining method, has been introduced to machine difficult-to-cut materials.18–22Compared with other machining methods,RUD enables many potentials,such as lower cutting force, less tool wear and smaller burr size.Yarar and Karabay23carried out RUD experiments on aluminum AA6061 with different heat treatments, and studied the effects of cutting parameters on machining.The results showed that ultrasonic vibration cutting can reduce cutting force and tool wear, and obtain better surface roughness at the same time.Li et al.24conducted an investigation on RUD of titanium alloys.Experimental results indicated that the thrust force, torque, cutting temperature, surface roughness, hole diameter deviation and burr height could be decreased by vibration assistance.Huang et al.25designed the RUD experiments of CFRP and systematically studied the tool wear behavior in machining.Compared with conventional machining, the introduction of ultrasonic vibration could reduce the average width of tool flank wear effectively,with a maximum reduction of 13.0%.Chang and Bone26compared burr size of aluminum which was machined by conventional drilling and RUD, and presented that RUD is superior in both burr height and width.Thus,Dong et al.27proposed a new method which combined RUD and robotic drilling.Their investigation realized the suppression of lateral chatter and the reduction of exit burr height by adopting ultrasonic vibration during robotic drilling, and the experimental results indicated that the stability lobe of robotic drilling was improved significantly and the burr height of aluminum plate was reduced by more than 40%.28.

    However, the research of robotic ultrasonic drilling is mainly based on the problems of drilling quality, machining efficiency and so on, without considering the impact of this method on assembly performance.Investigation on damage and residual stress of hole surface during RRUD is still deficient so far.In recent years, high-performance manufacturing has become the international research frontier.Damage and residual stress have crucial effect on the riveting fatigue life of CFRP/aluminum stacks.Previous studies focused on the suppression effect of ultrasonic vibration on the robotic drilling chatter and exit burr, but the influence of ultrasonic field energy on the residual stress formation, drilling damage and fatigue performance has not been clarified.Therefore,in order to improve the investigation on riveting fatigue life,this paper firstly presents a study on the hole surface strengthening mechanism and riveting fatigue life of CFRP /aluminum stacks in RRUD.It is organized in five sections.Following this introduction section, the procedure of developing the theoretical model for residual stress is shown in Section 2.Hole surface damage of CFRP drilling is analyzed in Section 3.In Section 4,fatigue life experiments of riveted joints are conducted to verify the effect of RRUD.Conclusions are drawn in Section 5.

    2.Residual stress analysis of RRUD

    2.1.Strengthening force

    2.1.1.Strengthening force from average drilling force

    In RRUD, the factors that influence the thrust force generation can be divided into two groups.One group consists of certain factors,including mechanical properties of material,uncut chip thickness and tool geometric sizes.The other includes uncertain factors, including the stability of RRUD system(Fig.1) and the cutting temperature.

    From previous study,28thrust force Fthcan be expressed as

    where Fth0represents the certain component of thrust force;K is the correction parameter denoting the uncertain component;n is the spindle speed; vfis the feed rate; A is the amplitude of ultrasonic vibration;is the thrust force on cutting edge;is the thrust force on chisel edge; d is the diameter of drill bit; d’ is the diameter of chisel edge; τAland σAlare the shear and yield stress of aluminum;ψ is the chisel angle;ω is the half chisel thickness; r is the distance to the tool center; E is the elastic modulus of aluminum; vAlis the Poisson ratio of aluminum; kchis the coefficient of chisel edge; γfis the angle between feeding direction and the force on chisel edge; γwis the wedge angle.

    In Eq.(1), acis the axial uncut chip thickness and can be calculated as

    where T is the period of rotational angle, T=; F is the frequency of ultrasonic vibration.

    As shown in Fig.2, aluminum hole surface is subjected to high-frequency extrusion by twist drill in RRUD.This compression is similar to the effect of strengthening process,which can change the internal stress on surface.Extrusion force perpendicular to the hole surface shown in Fig.2(a) is defined as the strengthening force.

    Fig.1 Structure and dynamic model of RRUD system.

    During RRUD, total strengthening force is generated by dynamic impact force and average cutting force, and it can be expressed as

    where Fqis the total strengthening force; Fqdand Fqsare the strengthening force formed by dynamic impact force and average cutting force, respectively.According to the force analysis in Fig.2(b),the total force FRcan be decomposed into Faand Fqin the axial and radial directions, and into Ffand Fnin the tangential and vertically tangential directions.By the meaning of angle relationship in the diagram,Fqsgenerated by the average cutting force can be calculated as

    In Eq (2), Fasis the average thrust force on cutting edge,and csand dsare the upper and lower limits of the axial force integral in the fillet part of the cutting tool,respectively.βcand γcare the reference angles shown in Fig.2(b),which can be calculated as

    where μ is the friction coefficient between tool and workpiece,u is the thickness of extrusion layer,and rqis the fillet radius of the twist drill.

    2.1.2.Strengthening force generated by vibration impact

    In the formation of hole surface strengthening force,the ultrasonic vibration effect is manifested by the introduction of instantaneous dynamic impact.Strengthening force by dynamic impact Fqdshould be analyzed with kinematics of RRUD.The displacement equation,axial velocity and acceleration at any point of cutting edge in RRUD are expressed as

    According to Newton’s second law of motion,the dynamic force of cutting tool fillet under high frequency vibration can be expressed as

    In Eq.(8),m=4/3ρπ (rq)3,which is the mass at tool fillet,and ρ is the density of tool material; Fadmaxis the maximum dynamic impact force.And the strengthening force by Fqddynamic impact can be obtained from the geometric relationship shown in Fig.2 (b):

    2.2.Cutting temperature

    During drilling,the work done by axial force is the fundamental cause of cutting temperature.29,30The linear velocity of the tool rotation is perpendicular to the axial force direction, so the linear velocity has no effect on the axial force work and can be ignored in this model.The axial velocity based on kinematics analysis of RRUD can be expressed as

    The power of the work done by the axial force generating heat can be expressed as

    Therefore, the heat flux q can be calculated by

    where Fthis the thrust force of aluminum in RRUD; Vais the axial velocity; η1is the proportionality factor of heat transfer.In this study, the cutting temperature is defined as the maximum temperature in drilling area.It is assumed that the heat distribution ratio of the workpiece and the tool is proportional to their thermal conductivity.Thus,η1=λ1/(λ1+λ2),where λ1and λ2are the thermal conductivity of aluminum and carbide tool.P is the thermal work equivalent, P = 4.1840 J/cal; S is the heat source area which can be expressed as

    According to the temperature field theory of continuous infinite heat source,31the maximum temperature of RRUD can be expressed as

    where q is the heat flux; λ1is the thermal conductivity of aluminum workpiece; β1is the thermal diffusivity of the workpiece, β1= λ1/cρAl,c is the specific heat capacity of the workpiece material, and ρAlis the density of the workpiece material; t1is the temperature rise time.

    Substitute Eq.(12)into Eq.(14),and the expression of temperature prediction model for RRUD of aluminum is shown as

    Substitute Eq.(1) and Eq.(13) into Eq.(15), and the relationship among cutting parameters, vibration parameters and cutting temperature can be expressed as

    2.3.Residual stress model in RRUD

    2.3.1.Elastic stress

    In RRUD,the contact between the tool fillet and the surface of hole surface is a special form of Hertz contact between two spheres,that is,one sphere is a flat plate.Deformation characteristics of Hertz contact are shown in Fig.3.Radial displacement of two elastic spheres when they impact and extrude against each other can be calculated as

    where U1(l) and U2(l) are the displacements of ball 1 in ball 2 and ball 2 in ball 1,respectively;ξ is the total displacement on the connecting line, ξ = d1+ d2; leis the equivalent contact radius, and can be expressed as

    Fig.3 Schematic diagram of deformation characteristics of two spheres in Hertz contact.

    where l1and l2are the radius of ball 1 and ball 2,respectively.

    When ball 2 is replaced by a flat plate, the radius l2is infinite.Assuming that the carbide tool is rigid, the displacement U2(l)of ball 2 on ball 1 at this point is close to 0.Therefore,the tool fillet and the extrusion displacement of hole surface can be expressed as

    According to Hertz contact theory,32the pressure distribution in contact area and the relationship between displacement and pressure can be expressed as follows:

    where P(l)is the contact pressure in the direction of l;Pmis the maximum pressure at the center-connecting line;lmis the maximum contact radius;Eeis the equivalent elastic modulus,and can be obtained by

    where E1and E2are the elastic modulus of the two contacts,and ν1and ν2are the Poisson ratio of two contacts.

    According to the research results of Teimouri and Amini,33the relationship among the maximum pressure Pm, the maximum contact radius lmand the hole surface strengthening force Fqcan be expressed as

    Based on the elastic stress theory of Hertz contact,the elastic stress expression along the hole surface direction is

    where σea, σetand σeqare the elastic stresses in the axial, tangential and hole surface thickness directions of the tool,respectively;νAlis the elastic modulus of aluminum material;s is the thickness value from the hole surface.Von Mises yield rule is used to analyze the internal stress of the material, and equivalent von Mises elastic stress and strain are described as follows:

    2.3.2.Plastic stress

    Based on the theory of linear relationship between elastic strain and elastic–plastic strain,34elastic–plastic strain can be expressed as

    where c1is the ratio of plastic deformation to elastic deformation.It can be expressed as

    According to the Hertz contact principle, the elastic deformation diameter decan be obtained by

    At the same time, the plastic deformation diameter can be calculated by

    where Sbis the ultimate stress of aluminum material.

    Introduction of ultrasonic vibration makes tool and hole surface press repeatedly at high frequency.And the hole surface material shows an extremely high strain rate, which can be calculated as

    The vibration frequency of ultrasonic machining is usually above 20000 Hz.According to Eq.(32),the material strain rate caused by ultrasonic impact is above 6.28×104s-1.When the strain rate is above 104,Johnson-Cook(J-C)model is the most suitable method for calculating the flow stress on the workpiece surface.35

    where A, B, C, m and n are the material physical constants of the J-C model.A is the initial yield stress of the workpiece material; B is the hardening constant; C is the strain rate constant;m is the thermal soft index;n is the hardening index;θ is the cutting temperature of RRUD;Tris room temperature;Tmis the melting point of aluminum.

    2.3.3.Stress release and equalization

    As shown in Fig.4, stress–strain curve (red curve) of metallic materials can be simplified to multi-segment lines.The ultimate elastic stress based on isotropic hardening theory can be expressed as

    In order to facilitate the calculation of residual stress, this study makes the following assumptions: (1) Deformation of hole surface is small; (2) Unloading process is elastic until the reverse yielding begins; (3) Static pressure does not cause plastic deformation.Based on the above assumptions, preresidual stress can be calculated by elastic stress and plastic deviatoric stress:

    For isotropic materials, the above expression can be expressed as36

    According to the material model shown in Fig.4, reverse yielding and hardening occur when σe≥2σp.At this point,the pre-residual stress can be calculated as

    Based on the findings of Teimouri33and Li et al.34, the residual stress type applied on the surface is a complete plane stress when the whole hole surface is subjected to multiple impacts.That is, the above stress is only applied to the depth of material (in the direction of hole surface thickness), and there is no change in the axial and tangential direction of tool.This indicates that the strain value in the plane is zero,and the strain-dependent nonzero value only occurs in the direction of the hole surface thickness.Therefore, the only possible preresidual stress and strain boundary conditions are

    From Eq.(37)and Eq.(38), it is clear that the pre-residual stresses and strains do not satisfy the above equilibrium boundary conditions, and stress relaxation must be generated.The corresponding relaxationandcan be obtained according to Hooke’s law:

    Therefore,the surface residual stress model of robot rotary ultrasonic drilling hole surface can be expressed as

    From the above analysis on residual stress theoretical model of hole surface, it can be seen that the ultrasonic vibration enhances the extrusion effect of tool on hole surface.Under the extrusion effect on hole surface,plastic deformation occurs in the machining area.After the unloading of external load, the material in the plastic deformation area cannot be recovered.The grain dislocation and slip are generated inside the workpiece.In addition, by Johnson-Cook continuity model, the plastic stress is greatly affected by the cutting temperature during dry drilling, and the uneven temperature field distribution is formed in the surface and interior of the hole surface material.The volume expansion trend of surface material is restricted by internal material, which generates residual stress on the surface.

    2.4.Validation experiment of residual stress model in RRUD

    2.4.1.Experimental setup

    Drilling experiments are performed on the robotic machining system shown in Fig.5.Laser vibrometer (Polytec OFV-5000)is used to calibrate the ultrasonic amplitude of tool tip.The stacks are composed of CFRP (CCF300) and aluminum(Al7075).Material physical constants of J-C model for aluminum are shown in Table 1.And the thermal conductivity λ1is 173 W/(m?K), heat capacity c is 0.96 J/(g?K), and room temperature is 20°C.Mechanical properties of CFRP and aluminum materials are shown in Table 2 and Table 3 respectively.The tool used in the experiment is a carbide twist drill with a diameter of 4.5 mm.Tool geometric parameters are shown in Table 4.Both single factor experiments and orthogonal experiments are carried out to verify the accuracy of residual stress model.Experimental parameters and levels are shown in Table 5 and Table 6.Residual stress on aluminum hole surface is measured by a X-ray stress tester (X-350A).Laser microscope(KEYENCE VK-X150)and field emission scanning electron microscope(SEM Quant 250FEG)are used to observe the hole surface damage and outlet damage of CFRP.

    2.4.2.Experimental results and analysis

    Measurement results of residual stress under RRUD and RCD are shown in Fig.6.It can be seen that residual stresses at the aluminum hole surface in RRUD are negative, and the intakeof ultrasonic vibration energy causes residual stresses to change from tensile stresses (positive values) to compressive stresses (negative values).

    Property Value EAl (MPa) 7.17 × 104 τAl (MPa) 150 σAl (MPa) 572 vAl 0.32 Density (g/cm3) 2.81

    Parameter Value Parameter Value d(mm) 4.5 2w(mm) 0.81 p(deg) 59 γw (deg) 30 ψ(deg) 54 γf (deg) 15 β(deg) 30

    Fig.5 Robotic rotary ultrasonic drilling system.

    Factor Unit Level A B C D E F G H Spindle speed r/min 1500 2000 2500 3000 3500 4000 4500 5000 Feed rate mm/s 0.8 1.0 1.1 1.3 1.4 1.6 1.7 1.9 Ultrasonic amplitude μm 4 6 8 10 12 14

    Factor Unit Level A B C D Spindle speed r/min 1500 2500 3500 4500 Feed rate mm/s 0.8 1.1 1.4 1.7 Ultrasonic amplitude μm 8 10 12 14

    Fig.6(a) shows the variation of residual stress with the spindle speed.It can be found that the residual stresses tend to decrease with the increase of spindle speed for both machining methods, and the residual stresses decrease more when the speed increases from 1500 r/min to 3500 r/min for RCD.Therefore, RCD at low spindle speed has a greater effect on the residual stresses.In contrast, the curve change is relatively flat in RRUD, and the residual stresses are relatively weakly influenced by the spindle speed.

    The effect of feed rate on the residual stress is shown in Fig.6(b).Residual stresses on the hole surface increase when the feed rate raises, and this trend is observed in both RCD and RRUD.This is due to the fact that an increase in feed rate indicates increment in feed per revolution and material removal per unit time, eventually causing an increase in the average cutting force.The raise in average cutting force contributes to the increase in residual stress.

    Fig.6(c) shows the effect of ultrasonic amplitude on the residual stress during RRUD.Residual stress tends to enhance significantly with the increasing amplitude.When the ultrasonic amplitude raises from 4 μm to 14 μm, the residual stress value increases from 39.4 MPa to 107.6 MPa,with an augmentation of 173.1%.This is because there is a positive correlation between the ultrasonic amplitude and the instantaneous vibration impact force.

    Substituting the RRUD parameters and material properties into theoretical model, the predictive values of residual stress can be obtained.The relative error of residual stress prediction is shown in Fig.7.The relative error of theoretical prediction ranges from 1.6% to 12.1%.Only one group has a prediction error of more than 10%, and rest of groups are within the acceptable range, with an average relative error of only 6.9%.Therefore,the results of experiments verified the validity and correctness of residual stress theoretical model.The model has a high prediction accuracy and can be used for prediction of residual stresses in RRUD.

    3.Drilling damage of CFRP in RRUD

    3.1.Comparison of CFRP hole surface damage

    As shown in Fig.8(a), severe defects of hole surface are observed in RCD, with a large depth of pits.The deepest pit is magnified(×1000)using SEM as shown in Fig.9(a).A large number of microscopic defects such as exposed fibers, fiber pull-out and micro-cracks exist in the pit.These defects limit the load-bearing capacity of hole surface greatly.In contrast,the surface of CFRP is relatively flat in RRUD.The same SEM magnification (×250) observation of hole surface is shown in Fig.8(b), which has only mild scratching, so the quality of hole surface is significantly improved.From Fig.8(d),we can see that, at spindle speed of 5000 r/min in RRUD,CFRP hole surface defects are minimal and the hole surface quality is the best.

    Fig.6 Comparison curve of residual stresses at hole surface by two processing methods.

    Fig.7 Prediction accuracy of theoretical model for residual stresses.

    Effects of ultrasonic vibration and feed speed on the defects are shown in Fig.10.Overall analysis of Fig.10(a)-(d) shows that RRUD decreases hole surface defects when the feed speed is varied.CFRP defects under both RCD and RRUD tend to deteriorate as the feed rate increases.SEM images (×1000) of the deepest position of the pit are shown in Fig.11(a)and(b),respectively.The pit of the defect in RCD is relatively deeper.Fiber fracture disorder, fiber voids and crack clustering in the pit are obvious,which will have a great impact on the riveting strength of stacks inevitably.In RRUD, the CFRP hole surface also shows certain defects.Fiber exposure phenomenon is more significant, but depth of the pit and the defect range are significantly smaller than those in RCD.These indicate that the effect of ultrasonic vibration can still inhibit the formation of defects in the robotic drilling at higher feed rate.

    Effect of ultrasonic amplitude on hole surface defects in RRUD is shown in Fig.12.Surface defects improve gradually as the ultrasonic amplitude increases from 8 μm to 12 μm.When ultrasonic amplitude is 8 μm, hole surface defects show a shallow scratch in depth.When the amplitude is raised to 10 μm, the area of scratch damage is rapidly reduced, and flat surface can be obtained when it is increased to 12 μm.However, when the ultrasonic amplitude is 14 μm, the surface defects of the hole surface appear again.The SEM magnification observation is shown in Fig.13.Defects are in the form of matrix spalling and fiber exposure, and deep pits disappear.

    3.2.Exit damage in RRUD

    CFRP drilling exit is prone to defect in the form of edge chipping and delamination, as shown in Fig.14.As can be seen in Fig.14(c), the depth of the edge crack is large, reaching more than 135 μm.Effect of ultrasonic vibration and spindle speed on the exit defect of robotic drilling is shown in Fig.15.A comprehensive comparison shows that RRUD decreases edge chipping and delamination, and obtains a relatively complete edge.In addition, with the gradual increase of the spindle speed, exit defects show a trend of improvement by both machining methods.This is because the thrust force of drilling decreases as the feed per tooth reduces.As a result, the tool force on the uncut layer at the exit is reduced and the drilling exit defect is improved.

    4.Fatigue life of riveted joints in RRUD

    4.1.Optimization of RRUD parameters

    In this section,optimization of RRUD parameters for CFRP/aluminums stacks is carried out.First,objective function is set according to the material removal efficiency and machining quality requirements.Then, the constraint conditions are determined according to the performance requirements of end-effector and ultrasonic machining device.The volume of material removed per unit time is usually used as the main indicator to characterize the cutting efficiency, which can be expressed as

    where Q is the material removal rate, which is the material removal volume in unit time; π is a constant; R is the drilling radius;n is the spindle speed;fris the feed per revolution.Penalty function method of particle swarm algorithm is chosen to transform this optimization mathematical model into an unconstrained problem.With the penalty function method,incorporating the constraints into the objective function yields:

    Fig.8 Effect of ultrasonic vibration and spindle speed on the shape of defects (Feed rate: 1.2 mm/s).

    Fig.10 Effect of ultrasonic vibration and feed rate on shape of defects.

    Particle swarm algorithm is used to solve this optimization mathematical model.The initial particle population is generated in a random way.The number of particles is set as 200 where each particle has four dimensions.The acceleration constants are all defined as 2 and the maximum number of iterations is set as 300.The optimization results are shown in Table 7.The frequency of the ultrasonic machining system is determined by the system structure and the design frequency is 20 kHz.

    4.2.Fatigue life experiments

    Fig.11 SEM inspection diagram of local amplification of CFRP hole surface defects (spindle speed: 3500 r/min, feed rate: 1.9 mm/s).

    Fig.12 Effect of ultrasonic vibration amplitude on hole surface defects.

    Fig.13 SEM inspection map of local amplification of CFRP hole surface defects (Ultrasonic amplitude: 14 μm).

    Single- and four-hole riveted joints of fatigue test workpiece are shown in Fig.16(a).The rivets are made from aluminum 2117-T4 with 5 mm diameter bar and 12 mm length, and are drilled using 5.1 mm diameter tool in accordance with the standard riveting technology.In order to avoid the influence of riveting process on fatigue life, the riveting process is completed with the assistance of aviation industry Jiangxi Hongdu Aviation Industry Group Co., Ltd., and the riveting process ensures the consistency of riveting force of each workpiece.Furthermore,the addition of the balance piece makes the load direction perpendicular to the rivet axial direction, which makes the working condition closer to the actual flight state of the aircraft.After riveting, the clamping of experimental parts and the site of fatigue test are shown in Fig.17.A universal fatigue test machine(LFV100kN)is selected for the fatigue test.The maximum static load is 100 kN, the maximum dynamic load is 50 kN,and the maximum frequency is 250 Hz.

    Fig.14 Typical shape of exit defect of CFRP drilling.

    Fig.15 Effect of ultrasonic vibration and spindle speed on exit defects.

    Tensile fatigue tests of single-hole and four-hole riveted stacks are all carried out at room temperature.80% of static tensile ultimate strength is selected as the maximum load stress.The stress ratio and tensile frequency are set as 0.1and 50 Hz respectively.When the displacement between plates of stacks exceeds 5% of the hole diameter, the joint is judged to fail.Each group of fatigue tests is repeated 20 times and any large deviation is averaged out to avoid chance error effect.

    Parameter Value Number of iterations 300 Spindle speed (r/min) 4489.99 Feed rate (mm/s) 1.63 Ultrasonic amplitude (μm) 10.43

    Fig.16 Structure of fatigue test workpiece.

    Fig.17 Clamping of test pieces and fatigue test site.

    4.2.1.Static tensile test results

    The static tensile experiment of riveted joints is firstly carried out to lay the foundation for the tensile fatigue test.Ten groups of riveted joints are tested with and without ultrasonic robotic drilling, and those groups with abnormal results are removed.The load–displacement curves for single-hole riveted joints with both RRUD and RCD are shown in Fig.18.

    At the initial stage, when the displacement is less than 0.05 mm,there is an approximately linear relationship between load and displacement.In addition to the elastic deformation of the assembly under tensile load, friction also exists at this point to resist shear load.This stage is mainly due to the gap fit between the nail hole and the rivet, and tensile load is mainly used to overcome the friction between the plates.There is no significant effect of process methods on the load–displacement curve at this stage, as the rivet does not exert an extrusion force on the hole surface.When the load exceeds the elastic limit σe, the curve enters the stable load-bearing phase.Compared with RCD, better quality of hole surface and residual compressive stress are obtained in RRUD processing, which inhibits crack initiation and expansion effectively, resulting in a higher load carrying capacity.The final stage of the curve is the damage and destruction stage of the material.By comparing the ultimate load σb, the static load limit of the riveted joint is increased by more than 12% via the action of ultrasonic vibration.

    4.2.2.Fatigue test results

    Average fatigue life of single-hole riveted CFRP/aluminum stacks processed by RCD is 731,237 cycles, compared with 1,231,275 cycles processed by RRUD.Fatigue life of singlehole riveted joints drilled by RRUD is increased by more than 68%.Average fatigue life of four-hole riveted stacks machined by RRUD and RCD are 713,627 cycles and 381,787 cycles respectively.The introduction of ultrasonic vibration makes an increase of more than 86% for the four-hole riveted joint.In addition, the fatigue life of the four-hole riveted joint is lower than that of single-hole riveted joint for both RRUD and RCD.In RCD, the fatigue life of the four-hole riveted joint is reduced by 47.8% compared with the single-hole riveted joint.The reduction in RRUD is 42.0%.As the number of holes raises,the likelihood of hole surface damage increases.Crack sprouting and fatigue damage in any riveted hole have an effect on the fatigue life of the overall joint.Therefore,fatigue life experiments effectively verify the effect of ultrasonic vibration on the drilling quality improvement of CFRP/aluminum stacks.

    Fig.18 Load-displacement curves of riveted joints.

    4.3.Analysis of fatigue damage morphology

    Laser microscopy is used to analyze the damage to the aluminum surface of the laminated joint after 380,000 cycles of tensile fatigue, and the results are shown in Fig.19.After tensile fatigue, the hole surfaces of the riveted joints show increased damage and cracking on the aluminum surface under repeated compression of the rivets.By comparing the results of surface inspection between RCD and RRUD,ultrasonic vibration effectively inhibits the expansion of damage and crack initiation on the aluminum hole surface.The degree of damage is significantly less than that of RCD.

    The reasons for the above phenomenon can be specifically analyzed from the following aspects: From the analysis of the reinforcement effect in RRUD,residual stresses in the hole surface are transformed from tensile stresses to compressive stresses, which effectively increase the load-bearing capacity of the hole surface and extend the fatigue life of stacks.In addition, from the analysis of the stability of RRUD, it is found that the machining state is stable and the hole diameter deviation is significantly improved.27Meantime, from the analysis of the exit burr of aluminum in RRUD, the effect of ultrasonic vibration on the reduction of the burr height is significant.The improvement of hole diameter and exit burr size has an important role for improvement of fatigue life.Combined with the inhibition of microscopic cracking by residual compressive stress,the superposition of the two effects effectively extends the fatigue life of the riveted joint.

    Fig.19 Surface profile of aluminum hole after 380,000 cycles of tensile fatigue.

    Fig.20 Surface morphology of CFRP hole after 380,000 cycles of tensile fatigue (×1000).

    CFRP hole surface damage after 380,000 cycles of tensile fatigue is observed as shown in Fig.20.The surface of CFRP hole exhibits severe damage under repeated compression by the rivet.By comparing the CFRP hole surface morphology of RRUD and RCD, it can be seen that the hole surface of RCD shows deep pits, cracks and fiber breaks after 380,000 cycles of tensile fatigue and the riveted joint reaches failure criteria.At the same number of fatigue cycles,CFRP hole surface with RRUD only shows micro-cracks and shallow pits, which is a significant improvement.

    5.Conclusions

    This study proposes the strengthening mechanism of residual stress and surface damage during RRUD of CFRP/aluminum stacks.Meantime, tensile fatigue experiment for riveted joints of CFRP/aluminum stacks is carried out to verify the improvement effect on the service performance.The main conclusions are drawn as follows:

    (1) A theoretical model of residual stresses with consideration of the combined effects of drilling temperature and thrust force is established.Residual compressive stresses are effectively formed by RRUD.Effect of ultrasonic amplitude on the residual compressive stresses is significant.When the ultrasonic amplitude is increased from 4 μm to 14 μm,the increase of the residual stress reaches 173.1%.

    (2) CFRP damages of hole surface and hole exit are significantly improved under the action of ultrasonic vibration.As the feed rate increases, the CFRP hole surface defects tend to deteriorate for both RCD and RRUD.As the spindle speed increases, the hole surface quality gradually improves.

    (3) Compared with RCD,the hole surface of RRUD shows a higher static load carrying capacity.The fatigue life of single-hole riveted joints is increased by 68%, and the fatigue life of four-hole riveted joints arranged according to aerospace design standards is increased by more than 86%.

    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 the Project on the Technological Leading Talent Teams Led by Frontiers Science Center for Complex Equipment System Dynamics (No.FSCCESD220401) and the Jiangsu Funding Program for Excellent Postdoctoral Talent (No.2022ZB264).

    日韩免费av在线播放| 国产91精品成人一区二区三区| 天堂√8在线中文| 99国产综合亚洲精品| 午夜福利视频1000在线观看 | 国产成+人综合+亚洲专区| 久久国产精品影院| 亚洲国产精品sss在线观看| 亚洲熟妇中文字幕五十中出| 国产精品久久久久久精品电影 | 久久精品国产亚洲av高清一级| 久久久久九九精品影院| 女人精品久久久久毛片| 久久国产亚洲av麻豆专区| 国产精品免费视频内射| 亚洲色图综合在线观看| 在线观看66精品国产| 久久国产精品男人的天堂亚洲| 欧美日本中文国产一区发布| 亚洲人成电影观看| 欧美另类亚洲清纯唯美| 国产精品电影一区二区三区| 夜夜看夜夜爽夜夜摸| 97人妻精品一区二区三区麻豆 | 天天躁狠狠躁夜夜躁狠狠躁| 免费高清在线观看日韩| 国产av又大| 国产亚洲欧美精品永久| 韩国av一区二区三区四区| 黄色成人免费大全| 亚洲精品久久成人aⅴ小说| 视频区欧美日本亚洲| 十八禁人妻一区二区| 久久草成人影院| 麻豆一二三区av精品| 亚洲人成网站在线播放欧美日韩| 久久精品国产亚洲av香蕉五月| av在线天堂中文字幕| 国产色视频综合| 国产单亲对白刺激| 精品久久久精品久久久| 亚洲国产欧美日韩在线播放| 一区二区日韩欧美中文字幕| 欧美成人性av电影在线观看| 亚洲第一欧美日韩一区二区三区| 色老头精品视频在线观看| 嫩草影视91久久| 日韩欧美三级三区| 这个男人来自地球电影免费观看| 高清在线国产一区| 夜夜看夜夜爽夜夜摸| 午夜福利一区二区在线看| 国产亚洲精品久久久久5区| 一进一出抽搐动态| 国产精品永久免费网站| 亚洲色图av天堂| 午夜影院日韩av| 国产精品九九99| 禁无遮挡网站| 亚洲情色 制服丝袜| 在线观看免费视频日本深夜| 国产91精品成人一区二区三区| 丝袜人妻中文字幕| 国产aⅴ精品一区二区三区波| 又黄又爽又免费观看的视频| av天堂在线播放| 午夜福利18| av片东京热男人的天堂| 日韩大码丰满熟妇| 桃红色精品国产亚洲av| 日本a在线网址| 久久久久久久久久久久大奶| 国产三级黄色录像| 亚洲无线在线观看| 51午夜福利影视在线观看| 国产一区二区三区综合在线观看| 午夜福利,免费看| 久久婷婷人人爽人人干人人爱 | 最近最新免费中文字幕在线| 日日干狠狠操夜夜爽| 美女高潮喷水抽搐中文字幕| 精品不卡国产一区二区三区| 日本免费一区二区三区高清不卡 | 一个人观看的视频www高清免费观看 | 国产欧美日韩综合在线一区二区| 久久久久久亚洲精品国产蜜桃av| 国产高清videossex| 亚洲一区二区三区不卡视频| 制服丝袜大香蕉在线| 亚洲中文av在线| 成人特级黄色片久久久久久久| 中文亚洲av片在线观看爽| 亚洲专区中文字幕在线| 嫩草影视91久久| 国产蜜桃级精品一区二区三区| 亚洲人成77777在线视频| 亚洲人成电影观看| 18禁黄网站禁片午夜丰满| 国产欧美日韩一区二区精品| 狂野欧美激情性xxxx| 我的亚洲天堂| 可以免费在线观看a视频的电影网站| 国产视频一区二区在线看| 欧美日韩瑟瑟在线播放| 亚洲人成网站在线播放欧美日韩| 涩涩av久久男人的天堂| 中文字幕人成人乱码亚洲影| 久热这里只有精品99| 一a级毛片在线观看| 老汉色av国产亚洲站长工具| 欧美日韩乱码在线| 亚洲精品国产精品久久久不卡| 亚洲中文字幕日韩| 天天一区二区日本电影三级 | 制服诱惑二区| 欧美成人性av电影在线观看| av网站免费在线观看视频| 狠狠狠狠99中文字幕| 久久久精品国产亚洲av高清涩受| 国产成人一区二区三区免费视频网站| 少妇被粗大的猛进出69影院| 国内毛片毛片毛片毛片毛片| 午夜免费激情av| www.www免费av| 自拍欧美九色日韩亚洲蝌蚪91| 国产成人av激情在线播放| 欧美绝顶高潮抽搐喷水| 久久精品亚洲熟妇少妇任你| 亚洲国产日韩欧美精品在线观看 | av免费在线观看网站| 男女午夜视频在线观看| 波多野结衣av一区二区av| 变态另类成人亚洲欧美熟女 | 亚洲第一青青草原| 91精品三级在线观看| 午夜福利免费观看在线| 少妇裸体淫交视频免费看高清 | 最好的美女福利视频网| 亚洲国产高清在线一区二区三 | 免费一级毛片在线播放高清视频 | 亚洲熟妇熟女久久| 国内精品久久久久久久电影| 亚洲专区字幕在线| 国产xxxxx性猛交| 两个人视频免费观看高清| 夜夜夜夜夜久久久久| 桃色一区二区三区在线观看| 不卡一级毛片| 日本五十路高清| 欧美激情 高清一区二区三区| 午夜福利18| 十八禁网站免费在线| 亚洲黑人精品在线| 国产成人一区二区三区免费视频网站| 九色国产91popny在线| 韩国av一区二区三区四区| 亚洲五月色婷婷综合| av天堂在线播放| 欧美黄色片欧美黄色片| 久久性视频一级片| 国产成人精品久久二区二区免费| 午夜视频精品福利| 91成年电影在线观看| 热99re8久久精品国产| 日日爽夜夜爽网站| 男人舔女人的私密视频| 欧美激情高清一区二区三区| e午夜精品久久久久久久| 久久久久久久久中文| 亚洲av成人不卡在线观看播放网| 久久人妻熟女aⅴ| 大码成人一级视频| 国产精品99久久99久久久不卡| 他把我摸到了高潮在线观看| 给我免费播放毛片高清在线观看| 在线国产一区二区在线| 欧美精品啪啪一区二区三区| 国产91精品成人一区二区三区| 两性夫妻黄色片| 国产99久久九九免费精品| 波多野结衣高清无吗| 俄罗斯特黄特色一大片| 成人三级做爰电影| 欧美日韩黄片免| 久久中文字幕一级| 久久伊人香网站| 日韩国内少妇激情av| 日韩精品免费视频一区二区三区| 好男人在线观看高清免费视频 | 伦理电影免费视频| 久久人人97超碰香蕉20202| 久久久国产精品麻豆| 免费不卡黄色视频| 熟妇人妻久久中文字幕3abv| 一进一出好大好爽视频| 两性夫妻黄色片| 午夜日韩欧美国产| 欧美黑人精品巨大| 日本一区二区免费在线视频| 欧美绝顶高潮抽搐喷水| 男女床上黄色一级片免费看| 日本一区二区免费在线视频| 天天一区二区日本电影三级 | 国产区一区二久久| 久热爱精品视频在线9| 欧美激情 高清一区二区三区| 久久香蕉激情| 色婷婷久久久亚洲欧美| av视频免费观看在线观看| 午夜精品在线福利| 超碰成人久久| 日韩欧美三级三区| 操美女的视频在线观看| 麻豆成人av在线观看| 久久中文看片网| 久99久视频精品免费| 亚洲最大成人中文| 99热只有精品国产| 欧美 亚洲 国产 日韩一| 国产又爽黄色视频| 亚洲色图 男人天堂 中文字幕| 9191精品国产免费久久| 久久婷婷人人爽人人干人人爱 | 欧美激情久久久久久爽电影 | 最新美女视频免费是黄的| 18禁国产床啪视频网站| 午夜福利视频1000在线观看 | 欧美黑人欧美精品刺激| 极品人妻少妇av视频| 成人永久免费在线观看视频| 人人妻人人爽人人添夜夜欢视频| 最新在线观看一区二区三区| 日日夜夜操网爽| 午夜影院日韩av| 老鸭窝网址在线观看| 国产精品久久电影中文字幕| 天天躁夜夜躁狠狠躁躁| 男女下面插进去视频免费观看| 国产aⅴ精品一区二区三区波| 日韩免费av在线播放| 久久人妻熟女aⅴ| 9色porny在线观看| 日韩 欧美 亚洲 中文字幕| 如日韩欧美国产精品一区二区三区| 成人欧美大片| 美女大奶头视频| 18禁裸乳无遮挡免费网站照片 | 50天的宝宝边吃奶边哭怎么回事| 国产亚洲精品第一综合不卡| 十八禁人妻一区二区| 男男h啪啪无遮挡| 在线观看一区二区三区| 热99re8久久精品国产| 欧美日韩精品网址| 天堂动漫精品| 欧美日韩亚洲国产一区二区在线观看| 99久久久亚洲精品蜜臀av| 琪琪午夜伦伦电影理论片6080| 最新在线观看一区二区三区| 成人特级黄色片久久久久久久| 免费高清在线观看日韩| 91大片在线观看| 亚洲av成人不卡在线观看播放网| 丝袜人妻中文字幕| 一本综合久久免费| 免费高清视频大片| www.熟女人妻精品国产| 日本五十路高清| 一级,二级,三级黄色视频| 高潮久久久久久久久久久不卡| 美女大奶头视频| 69精品国产乱码久久久| 免费无遮挡裸体视频| 久久精品人人爽人人爽视色| 天堂影院成人在线观看| 在线观看免费视频网站a站| 欧美日韩福利视频一区二区| 岛国在线观看网站| 亚洲av电影在线进入| 亚洲第一电影网av| 一级a爱片免费观看的视频| cao死你这个sao货| 免费一级毛片在线播放高清视频 | 美女扒开内裤让男人捅视频| 最好的美女福利视频网| av福利片在线| 一级a爱视频在线免费观看| 久久久精品欧美日韩精品| 91大片在线观看| 亚洲欧美日韩无卡精品| 无人区码免费观看不卡| 日本黄色视频三级网站网址| 久久久久久久久久久久大奶| 熟妇人妻久久中文字幕3abv| 女性生殖器流出的白浆| 伦理电影免费视频| 精品久久久久久,| 国产色视频综合| 欧美 亚洲 国产 日韩一| 99国产精品免费福利视频| 色播亚洲综合网| 成人精品一区二区免费| 宅男免费午夜| 久久久国产成人精品二区| 亚洲电影在线观看av| 国产av精品麻豆| 天堂影院成人在线观看| 97人妻精品一区二区三区麻豆 | 男女之事视频高清在线观看| 亚洲欧美精品综合久久99| 亚洲一区二区三区色噜噜| 国产高清有码在线观看视频 | 国产一区在线观看成人免费| 成人手机av| 99国产精品一区二区三区| 免费在线观看完整版高清| 国产单亲对白刺激| 久久久久久人人人人人| 久久久久九九精品影院| 国产免费男女视频| 男女做爰动态图高潮gif福利片 | 国产成人av激情在线播放| 午夜福利欧美成人| 日本撒尿小便嘘嘘汇集6| 国产一级毛片七仙女欲春2 | 精品一区二区三区四区五区乱码| 亚洲一区二区三区色噜噜| 亚洲人成电影观看| 精品国产亚洲在线| 久久九九热精品免费| 午夜激情av网站| 国产成人精品无人区| 亚洲一区二区三区不卡视频| or卡值多少钱| 国产亚洲精品一区二区www| 国产成人精品在线电影| 黑人操中国人逼视频| 国产精品日韩av在线免费观看 | 国产三级黄色录像| 欧美久久黑人一区二区| 亚洲av美国av| 日韩欧美在线二视频| 午夜福利18| 免费不卡黄色视频| 日本免费一区二区三区高清不卡 | 久久国产精品影院| 亚洲av第一区精品v没综合| 色播在线永久视频| 亚洲精品国产精品久久久不卡| av天堂久久9| 久久国产精品男人的天堂亚洲| 99精品久久久久人妻精品| 成人国产综合亚洲| 18禁美女被吸乳视频| 日日摸夜夜添夜夜添小说| 亚洲男人天堂网一区| 精品久久久久久成人av| 在线视频色国产色| 性欧美人与动物交配| 亚洲 国产 在线| 丝袜人妻中文字幕| 搡老熟女国产l中国老女人| 19禁男女啪啪无遮挡网站| 国产aⅴ精品一区二区三区波| avwww免费| 夜夜夜夜夜久久久久| 亚洲全国av大片| 91国产中文字幕| 老司机在亚洲福利影院| 色播亚洲综合网| 亚洲午夜理论影院| 精品乱码久久久久久99久播| 亚洲国产精品久久男人天堂| 97人妻精品一区二区三区麻豆 | 老司机午夜十八禁免费视频| 男人舔女人的私密视频| 亚洲人成电影免费在线| 亚洲成人久久性| 久久九九热精品免费| 一边摸一边做爽爽视频免费| 丁香六月欧美| 国产1区2区3区精品| 99国产精品免费福利视频| 成人18禁高潮啪啪吃奶动态图| 日韩精品免费视频一区二区三区| 欧美一级毛片孕妇| 国产又色又爽无遮挡免费看| 国产精品亚洲av一区麻豆| 91成人精品电影| 日本vs欧美在线观看视频| 亚洲全国av大片| 中文字幕人妻熟女乱码| 午夜精品国产一区二区电影| 中亚洲国语对白在线视频| 久热爱精品视频在线9| av网站免费在线观看视频| 亚洲自偷自拍图片 自拍| 黄色视频,在线免费观看| 久久 成人 亚洲| 国产黄a三级三级三级人| 欧美黄色淫秽网站| 欧美亚洲日本最大视频资源| 69精品国产乱码久久久| 不卡一级毛片| 欧美色视频一区免费| 妹子高潮喷水视频| 91国产中文字幕| 黄网站色视频无遮挡免费观看| 日韩高清综合在线| 黑人巨大精品欧美一区二区mp4| av超薄肉色丝袜交足视频| 手机成人av网站| 狠狠狠狠99中文字幕| 欧美 亚洲 国产 日韩一| 黄色视频,在线免费观看| 国产1区2区3区精品| 国产又色又爽无遮挡免费看| 国产精品久久电影中文字幕| 在线观看一区二区三区| 91麻豆av在线| 亚洲国产毛片av蜜桃av| 成年女人毛片免费观看观看9| 亚洲成人精品中文字幕电影| 97超级碰碰碰精品色视频在线观看| 日韩精品青青久久久久久| 国产区一区二久久| 美女大奶头视频| 午夜久久久在线观看| 国产主播在线观看一区二区| 久久久久精品国产欧美久久久| 国产免费av片在线观看野外av| 日本三级黄在线观看| 一区二区日韩欧美中文字幕| svipshipincom国产片| av视频免费观看在线观看| 成人18禁在线播放| 日本 av在线| 老司机午夜福利在线观看视频| 美女免费视频网站| 法律面前人人平等表现在哪些方面| 亚洲国产精品999在线| 国产国语露脸激情在线看| 国产精品综合久久久久久久免费 | 一本综合久久免费| 99国产综合亚洲精品| 国产成+人综合+亚洲专区| 国产伦人伦偷精品视频| 一本综合久久免费| 久久香蕉国产精品| 中文字幕高清在线视频| 日韩视频一区二区在线观看| 一本综合久久免费| 亚洲av第一区精品v没综合| 免费一级毛片在线播放高清视频 | 在线播放国产精品三级| 日本欧美视频一区| 欧美日韩一级在线毛片| 亚洲欧美精品综合一区二区三区| 欧美 亚洲 国产 日韩一| av在线播放免费不卡| 一区在线观看完整版| 亚洲国产精品合色在线| 一二三四社区在线视频社区8| 精品午夜福利视频在线观看一区| 老熟妇仑乱视频hdxx| 欧美在线一区亚洲| 看黄色毛片网站| 日韩 欧美 亚洲 中文字幕| 久9热在线精品视频| 欧美日韩瑟瑟在线播放| 亚洲男人天堂网一区| 夜夜夜夜夜久久久久| 久久婷婷成人综合色麻豆| 深夜精品福利| 久久国产精品人妻蜜桃| 久久精品aⅴ一区二区三区四区| 欧美日韩一级在线毛片| 久久精品亚洲熟妇少妇任你| 一a级毛片在线观看| 免费在线观看黄色视频的| 国产精品一区二区精品视频观看| 国产一区二区三区在线臀色熟女| 亚洲国产精品久久男人天堂| 亚洲少妇的诱惑av| 久久久国产成人免费| 国产精品99久久99久久久不卡| www.999成人在线观看| 午夜成年电影在线免费观看| 精品人妻1区二区| 亚洲色图 男人天堂 中文字幕| 老司机在亚洲福利影院| 俄罗斯特黄特色一大片| 黑人欧美特级aaaaaa片| 日本欧美视频一区| 国产国语露脸激情在线看| 日本黄色视频三级网站网址| 免费人成视频x8x8入口观看| 中文字幕久久专区| 深夜精品福利| 搡老熟女国产l中国老女人| 亚洲精华国产精华精| 欧美亚洲日本最大视频资源| 十分钟在线观看高清视频www| 久久婷婷人人爽人人干人人爱 | 免费av毛片视频| 日日爽夜夜爽网站| 满18在线观看网站| 色哟哟哟哟哟哟| 欧美中文日本在线观看视频| 亚洲第一青青草原| 国产亚洲精品av在线| 一级片免费观看大全| 国产成人欧美| 精品欧美一区二区三区在线| 国产精品,欧美在线| 亚洲av片天天在线观看| 亚洲精华国产精华精| 久久人人97超碰香蕉20202| 精品久久久久久成人av| 欧美激情极品国产一区二区三区| 女性被躁到高潮视频| 国产欧美日韩一区二区精品| 丁香欧美五月| 午夜免费鲁丝| 亚洲av第一区精品v没综合| 欧美日韩黄片免| 少妇的丰满在线观看| 国产成年人精品一区二区| 成人永久免费在线观看视频| 男人的好看免费观看在线视频 | 日本免费一区二区三区高清不卡 | 欧美日本中文国产一区发布| 亚洲最大成人中文| 69精品国产乱码久久久| 国产精品香港三级国产av潘金莲| 窝窝影院91人妻| 欧美成人性av电影在线观看| 悠悠久久av| 精品久久蜜臀av无| 88av欧美| 国产成+人综合+亚洲专区| 禁无遮挡网站| 国产成人系列免费观看| 日韩大尺度精品在线看网址 | 动漫黄色视频在线观看| 19禁男女啪啪无遮挡网站| 女人爽到高潮嗷嗷叫在线视频| 丝袜人妻中文字幕| 亚洲美女黄片视频| 亚洲在线自拍视频| 满18在线观看网站| 纯流量卡能插随身wifi吗| 欧美中文综合在线视频| 久久久久久国产a免费观看| 精品福利观看| 老司机午夜十八禁免费视频| 国产精品,欧美在线| 熟妇人妻久久中文字幕3abv| 男人操女人黄网站| 午夜a级毛片| 国产成人精品久久二区二区免费| 免费看十八禁软件| 亚洲自拍偷在线| 国产aⅴ精品一区二区三区波| 91av网站免费观看| 悠悠久久av| 亚洲国产精品久久男人天堂| 嫩草影院精品99| 国产亚洲精品久久久久久毛片| 最新在线观看一区二区三区| 中文字幕人妻丝袜一区二区| 色老头精品视频在线观看| 国产成人av教育| 亚洲成人国产一区在线观看| 又黄又粗又硬又大视频| 别揉我奶头~嗯~啊~动态视频| 久久精品aⅴ一区二区三区四区| 99久久99久久久精品蜜桃| 禁无遮挡网站| 国产不卡一卡二| 色播在线永久视频| 成人手机av| avwww免费| 91在线观看av| 男人操女人黄网站| 精品人妻在线不人妻| 女人爽到高潮嗷嗷叫在线视频| 后天国语完整版免费观看| 午夜a级毛片| 韩国av一区二区三区四区| 久久午夜亚洲精品久久| 精品久久久久久,| 可以在线观看毛片的网站| 视频在线观看一区二区三区| 午夜福利18| 日日爽夜夜爽网站| 99riav亚洲国产免费| 国产又爽黄色视频| 国产成人精品久久二区二区免费| 国产精华一区二区三区| 男人舔女人下体高潮全视频| 精品高清国产在线一区| 久久九九热精品免费| 日本 av在线| 一本综合久久免费| av福利片在线| 国产精品免费一区二区三区在线| 成人永久免费在线观看视频| 在线观看免费午夜福利视频| 18禁国产床啪视频网站| 亚洲国产精品合色在线| 日韩一卡2卡3卡4卡2021年| 亚洲男人天堂网一区| 97人妻天天添夜夜摸| 久久久国产成人精品二区|