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

    Investigation on the interface damage in drilling low-stiあness CFRP/Ti stacks

    2019-10-26 01:15:34BinLUOKaifuZHANGShunuanLIUHuiCHENGRunxiaoWANG
    CHINESE JOURNAL OF AERONAUTICS 2019年9期

    Bin LUO, Kaifu ZHANG, Shunuan LIU, Hui CHENG, Runxiao WANG

    School of Mechanical Engineering, Northwestern Polytechnical University, Xi'an, Shaanxi 710072, China

    KEYWORDS Coupling effects;Drilling;Interface damage;Interlayer gap;Thin walled structures

    Abstract Carbon fiber reinforced plastic and titanium alloy (CFRP/Ti) stacks have been widely used as aerospace structures because of their excellent combination of physical properties.Interface damage caused by interface gaps, significantly different from that of metal/metal stacks, is a common problem in the through-hole drilling of CFRP/Ti stacks with low stiffness. In this study, a force-deformation coupling model was developed to further examine the formation mechanism and the control method of interface damage. Firstly, the coupling model was built considering the interaction between the thrust force and the deformation. To solve this model, a numerical method was proposed in which specific cutting coefficients were calibrated using only the thrust force of rigid stacks. Secondly, drilling experiments were performed with different feed rates and bending stiffness. Experimental results indicate that interface damage mainly includes interlayer chips and surface damage of CFRP layers. The surface damage, which is irreparable, is caused by the rotary extension of metal chips along the interlayer gap. Thirdly, variations of the interface gap were calculated with the coupling model that had been verified by measured thrust forces.The damage area was found to have a linear dependence relation with the interlayer gap. However, in conditions of large gap sizes, the interface damage areas increased with the interlayer gap at high feed rates,while decreasing slightly at low feed rates.This phenomenon was satisfactorily explained by the presented model. Finally, a method was proposed to determine the appropriate pressure exceeding which no interlayer damage will occur. Additional drilling experiments proved the method effective.This study leads to further understanding of the forming mechanism of interlayer damage and of selecting appropriate parameters in drilling low-stiffness composite/metal stacks.?2019 Chinese Society of Aeronautics and Astronautics.Production and hosting by Elsevier Ltd.This is an open access article under the CC BY-NC-ND license(http://creativecommons.org/licenses/by-nc-nd/4.0/).

    1. Introduction

    In modern aircraft manufacturing processes, connecting the thin-wall skin, the stringer, and the frame to the rigid components (e.g. wings or fuselages) is one of the most important procedures.Thousands of holes per assembly are unavoidably drilled for bolting or riveting in the assembly procedure.1Reducing the damage to the holes is crucial to the assembly quality. To achieve higher machining quality and efficiency,the robotic drilling technology has increasingly replaced the traditional manual drilling in the industry.2Components are usually stacked in joint positions and drilled in a single operation of the robotic drilling system. Such one-shot drilling processes not only save time but ensure the proper alignment of the holes on the layers. However, a number of drilling problems emerge in the one-shot drilling of low-stiffness stacked structures, drawing considerable attention from researchers.

    Interfacial burrs are a common problem in the throughhole drilling of metal/metal stacked sheets. Removing these burrs needs non-value-added operations, including disassembly, deburring, and reassembly. Experimental results showed that the height of interfacial burrs was reduced when the drilling position was near the clamp, or when the clamping force was acted on the surface of the stacks.3Since the interlayer gap was mainly affected by the clamping force,the thrust force and the bending stiffness,the gap size was considered to have a significant influence on the formation of interlayer burrs.4A number of simplified models were developed to predict the interlayer gap and analyze its influence.5,6The results indicated that decreasing the gap size through preload pressing was an efficient way to control the interlayer burr height to a satisfying value in the drilling of low-stiffness metal/metal stacks.Non-coaxiality was another problem caused by the interlayer gap in drilling stacked structures.7Experimental work and theoretical analyses were conducted to further understand the interlayer gap formations and non-coaxiality occurrences. It was discovered that increasing the clamping force cannot remarkably reduce the gap size when the force was beyond a threshold.The critical clamping force was proposed in the drilling of aluminum stacks, under which both the interlayer gap and burrs were eliminated.8Gao et al.9claimed that using pre-installed fasteners or decreasing the thrust force could also help reduce interlayer gaps and non-coaxiality.

    Carbon fiber reinforced plastic (CFRP) has been widely used in the modern aerospace industry, and is usually stacked with metal materials. Among the available configurations,CFRP/Ti(Titanium)coupling is the most popular one because of its excellent combination of metallurgical and physical properties.10CFRP/Ti stacks are considered as typical difficult-to-cut materials, whose drilling damage is very different from that of metal/metal stacks. Special issues may arise from severe subsurface damage,rapid tool wear,and excessive interface defects.11Cutting parameters are optimized to improve the drilling quality. Different cutting parameters are applied to different materials.12,13New materials and special geometry are developed to extend the tool life.14,15According to the drilling experiments of low-stiffness CFRP/Al (Aluminium) stacks, non-coaxiality is observed and tends to be more serious than that of the Al/Al stack.16However,interface damage, similar to interfacial burrs of metal/metal stacks, has not yet been studied for the drilling process of composite/metal stacks.

    For any low-stiffness stacks,a high clamping force or a low thrust force helps eliminate interlayer defects caused by the interlayer gaps. A large number of studies developed experimental or theoretical models to determine appropriate clamping forces that match the thrust force and the bending stiffness of the stacks.4-8In effect, the interface gap changes with the dynamic thrust force in drilling stacks. A typical thrust force signal includes several stages according to the cutting materials of drill edges.The variation of thrust forces was ignored in the presented models. In addition, the thrust force is significantly affected by the deformation of work pieces17,18coupled with each other.19Based on this, the modeling of thrust forces and interlayer gaps need to be further studied in the drilling of low-stiffness composite/metal stacks, leading to further understanding of the formation mechanism and the elimination method of interface defects.

    This paper studies the coupling law between the thrust force and the structure deformation in the drilling of low-stiffness CFRP/Ti stacks, making possible the exploration of the formation mechanism of the interface gap and the interface damage. A coupling model was first built considering the interactions between the thrust force and the deformation,followed by the presentation of its numerical solution process. A novel calibration method of cutting coefficients was developed independently through drilling tests of rigid stacks. Secondly,experimental works were carried out to study the interface damage and calibrate the coupling model.Thirdly,the verified model was used to predict the variation of the interface gap.Its influence on the formation of the interlayer damage was also investigated. Finally, a method was proposed to determine an appropriate pressure for eliminating the interlayer damage with a minimum deflection of the stacks.

    2. Force-deformation coupling model in drilling low-stiffness stacks

    2.1. Description of drilling process of low-stiffness stacks

    In a rigid stack drilling process, deformation does not occur under the action of thrust forces. The process is divided into three stages according to tool-work interactions: (a) the composite layer drilling stage,(b)the interface region drilling stage,and (c) the metal layer drilling stage, as shown in Fig. 1 (a).There is no gap or damage at the interlayer region. However,deformation affects the relative feed speed between the cutter and the workpiece in the drilling of low-stiffness stacks.19The stage division of the drilling process should consider both interactions and deformation. Four specific moments are defined in the drilling process.‘‘Contact”refers to the moment when the drill edge touches the second layer;‘‘separation”represents the moment when the interaction force of the layers reduces to zero,where deformation of the second layer exceeds that of the first layer;‘‘entrance”denotes the moment when the drill edge fully engages in the second layer;‘‘springback”refers to the moment when deformation of the second layer reduces.Thus,the drilling process of low-stiffness stacks is divided into five stages based on the four specific moments as shown in Fig. 1(b).

    In stage I,only the first layer is drilled.The deflection of the two layers is the same and increases over time.In stage II,both layers are drilled. The thrust force of the first layer decreases,while that of the second layer increases.The deflection of each layer differs at the end of this stage. In stage III, both layers are still drilled with the deflection of the first layer reducing to zero while that of the second layer increasing.The interlayer gap occurs and increases. A part of the drill edge does not cut the stack because of the gap.In stage IV,only the second layer is drilled. The deflection of the layer as well as the interlayer gap increases to the maximum. In stage V, the second layer is drilled. Both its deflection and the interlayer gap reduce to zero. The region in which stages II and III occur is named interface region where the drill bit cuts both layers.The thrust force in this region equals the sum of the element forces from each layer.

    Fig. 1 Multi-stage drilling process.

    2.2. Calibration of cutting coefficients

    The thrust force in the drilling stack needs to be known before the calculation of the deflection. Since the cutting angles vary along the drill edge, it is difficult to directly model the thrust force. Therefore, the infinitesimal method is proposed, in which the drill edge is divided into small elements. Each element has a constant cutting angle, and is assumed to be orthogonal or oblique cutting. Luis et al. proposed that material removal was the plastic extrusion in the region around the center of the chisel edge, called the indentation zone.20At the remaining portion, material removal was found to be orthogonal cutting with a large negative rake angle.21Since the radius of the indentation zone is much smaller than the chisel length,the orthogonal cutting could be extended to the whole chisel edge.

    The cutting force of the element is proportional to the cutting area.The ratio is defined as the specific cutting coefficient(kc), represented through a power law relationship using the uncut thickness and the rake angle and spindle speed.21Our previous work concluded that the law relationship was effective for the drilling process of the CFRP/Ti stacks.22Since the cutting depth is proportional to the feed rate,and the spindle speed has little effect on the thrust force,22the specific cutting coefficient could be expressed as:

    where γ is the rake angle of the cutting element, and kc1, kc2,and kc3are the coefficients related to the material property.The complexity of the cutting process leads to the lack of methods to obtain coefficients from the material properties.Therefore,the coefficients are always calibrated through experimental data.

    For element M at radius distance r, the elementary forces can be decomposed into the cutting direction (dFc) and the feed direction (dFfd). Thrust force Fthis obtained by summing force dFfdin the axial direction for all the elements, as shown in Fig.2.Then,a simple and effective model of the thrust force can be written as19

    where i is the inclination angle, p the point angle of the drill bit, r the distance of the considered point M on the drill edge from the drill axis, g1and g2are the factors related only to the drill geometry, and r1and r2are the limits of integration, as shown in Fig. 2. Rake angle γ and inclination angle i are calculated using equations in Ref. 22. Since the woven-ply or the cross-ply composite can be assumed as an isotropic material, the cutting force model is effective for the woven CFRP.22

    Since the cutting force in the interface region equals the sum of the forces from both layers,the thrust force in the stack drilling (Fth,stk) can be expressed as

    where Fth,1stis the thrust force in the first layer drilling, Fth,2stis the thrust force in the second layer drilling,h1is the height of the first layer, h2the height of the second layer, vfthe feed speed in mm/s, and hlipis the height of the drill edges which is written as

    Fig. 2 Calculation of thrust force in drilling process.

    where R is the radius of the drill,Rch denotes half of the chisel length, ψ the chisel angle, and ω is half of the web thickness.If the machining time is set to zero at the moment when the drill bits touch the stack, the maximum radius distance of the drill edges involved in the cutting at time t can be calculated by

    Then, the limits of integration in Eq. (3) can be written as

    In previous literature, coefficients kc1, kc2, and kc3need to be calibrated from the thrust force of the entrance stage in the single layer drilling process, indicating the necessity of a series of extra experiments. According to Eq. (3), the thrust force in the interface region equals that from both layers.Since the entrance stage is the inverse process of the exit stage in the well supported case, the thrust force of the second layer in the entrance stage could be separated from the stack drilling force. Based on this, a novel method is proposed to calibrate the specific cutting coefficient directly from the stack drilling force:

    where F1~F3are the thrust forces of the stacks when the integrating ranges are [0, Rch], [Rch, rmid], [Rch, R], respectively,during the drilling in the first layer, and F4~F6are those in the second layer, as shown in Fig. 3. rmidequals (R-Rch)/2.Machining times t1~t6are calculated through Eqs.(5)and(6).

    Fig. 3 Calibration of specific cutting coefficients in drilling stacks.

    2.3. Method for predicting thrust force and interlayer gap

    Since the primary cutting load in the drilling process is applied laterally to the axis of the components, the airframe components are simplified to beams.4,7The mode of deflection is primarily formed through bending. Thin rectangular plates with fixed ends are analyzed in this paper as shown in Fig. 1. It is complicated to obtain an analytical solution for the deformation of this structure. Hence, the finite element method is proposed to calculate the deflection of the stacks. A beam is regarded as an element with a node at each end. Each node undergoes a vertical displacement w and a rotation θ, with two degrees of freedom.The stiffness equation of the beam element can be written in a matrix form:

    or simply as:

    where E is the elastic modulus of the specimen, the area moment of the inertia I=bh3/12, b, h and l are respectively the width,the height,and the length of the stack,and Keis the stiffness matrix relating the nodal displacements matrix qeto the forces matrix Fe.

    The global stiffness matrix K and the global force vector F are assembled. If the number of the beam elements is an odd number n,the size of the assembled is 2(n+1).Boundary conditions are given as qr=0 (r=1, 2, 2n+1, 2n+2). Then,the global matrix element krris set to one, krsand ksrare set to zero (s=1, 2, ..., r-1, r+1, ..., n), and Fris also set to zero. The displacement vector is obtained by F/K. When the thrust force is on the geometric center, the deflection at drill point uthequals w(n+1)/2.Bending stiffness KBis defined as

    where uthis the deflection at the drilling point.

    If the thrust force acts on the first layer, both layers have the same deflection in the fixed end case. Then, the bending stiffness of the stack (KB,stk) equals the sum of the bending stiffness of each layer (KB,1st, KB,2st):

    Thrust force Fth,1stslowly reduces while Fth,2stincreases in stage II. When the deflection of the top layer caused by Fth,1stequals that of the bottom layer caused by Fth,2st,the two layers begin to separate. Limit of integration rin(t) is noted as rspat this moment. Then, rspcould be obtained by solving the following equation

    Before separation, the deflection of the two layers is the same, thus the interlayer gap of the stack is zero in stages I and II. After separation, the first layer resiles and the second layer continues bending. The interlayer gap appears in stages III, IV, and V, which can be expressed as

    where δ(t)is the interlayer gap,uth,1stthe deflection of the first layer, and uth,2stis the deflection of the second layer.

    Bending or release movements change the feed rate f between the drill bit and the stack in Eq. (2). The feed rate of the drill bit is defined as a set value f0.The feed rate is smaller than its set value in the bending stage and larger than that in the release stage. The feed rate in drilling the low-stiffness structure equals the derivation of the deflection subtracted from the set value as follows:19

    where nsis the spindle speed in rev/s.

    In stages I and II, the deflection of each layer is the same,hence their feed rates are the same (f1st=f2st). Therefore,the stack can be regarded as a single layer plate.By combining Eqs. (2), (3), (10), and (14), the coupling model of the thrust force and the deformation in the drilling of low-stiffness stacks can be written as

    Fig. 4 Calculation process of force-deformation coupling model.

    where Fth,stk(0) is the initial condition, and G1st(t) and G2st(t)are noted as

    After separation time tsp, the interlayer gap occurs. Part of the drill edges is not involved in the drilling process. The feed rates of each layer are different(f1st≠f2st).The drilling process of the two layers are independent from each other. Therefore,the thrust force can be expressed as

    where Fth,1st(tsp) and Fth,2st(tsp) are the initial conditions.

    A numerical method is employed to solve Eqs. (15) and(17),since the integration of G1st(t)and G2st(t)are rather complex to solve.Then,the finite element method is used to obtain the deformation of the stack. The solving process of the coupling model is shown in Fig. 4.

    3. Drilling experiments of low-stiffness CFRP/Ti stacks

    3.1. Description of experimental conditions

    Drilling experiments are carried out to verify the proposed model and study interlayer damage caused by interlayer gaps and cutting chips.Since the feed rate is the key factor affecting the thrust force,and bending stiffness is mainly determined by the length of the structures, the experiments are performed with specimens of different lengths and feed rates. The specimens are stacked into two layers, the first layer being CFRP,and the second Ti-6Al-4V.The CFRP/Ti stacks vary in length,leading to varying bending stiffness. The tool used in the experiment is a carbide twist drill with a straight cutting edge.The drilling point is the geometric center of the specimens with fixed ends. The detailed experimental conditions are listed in Table 1. Three feed rates and four lengths are selected. Zeroin length means the specimen is supported by a plate with an 8 mm diameter hole. In this case, the stiffness of the specimen is large enough so that its deflection can be ignored.

    Table 1 Conditions in drilling low-stiffness CFRP/Ti stacks.

    The full factorial method is used to design the experiments.Thus, all levels of variates could be investigated. All the tests are dry drilling operations with a constant spindle speed on a CNC milling machine (XKN713). A dynamometer (Kistler 9257A) and a charge amplifier are used to measure the thrust forces. Both ends of the stacks are fixed to the dynamometer by bolts and fixtures. Force signals are collected and saved by a data recorder. Fig. 5 is the experimental setup for the measurement of the thrust force in the drilling of the CFRP/Ti stack.

    3.2. Experimental results and calibration of the coupling model

    Delamination is known as the main defect in drilling composite laminates. However, it does not appear in stack drilling because of the support of the titanium layer.A camera is used to obtain the damage pictures of the CFRP/Ti stacks. A scale is introduced into the picture according to the hole diameter and the thickness of stacks.As can be seen from Fig.6(a),titanium chips accumulate at the interlayer during the drilling of the second layer of the low-stiffness stacks. The chips scratch and destroy the hole edge of the CFRP layer during the depositing process(Fig.6(b),6(c)),while little damage is found on the surface of the titanium layer (Fig. 6(d)). This is due to the direction of the chip flow, and the higher hardness of the titanium surface. As is shown in Fig. 6, the main interface defects of the low-stiffness CFRP/Ti stacks are the interlayer chips and the surface damage of the CFRP layer. The former can be removed by the reassembly operation,while the latter is irreparable, hence must be avoided.

    Fig. 7 shows the effect of bending stiffness (i.e. specimen length) on the thrust force and the interface damage at f0=0.05 mm/r. As the bending stiffness decreases, the increase of the thrust force slows down, while the maximum thrust force reduces about 16%. The edge damage of the CFRP hole deepens with the decrease of the bending stiffness.

    Three curves of the thrust forces at l=0 mm are used to obtain forces F1~F6shown in Fig.3.Then,the values are substituted into Eq.(7).The logarithm of both sides and the least square method are employed to solve the equation. The calibration results are listed in Table 2.

    Fig. 5 Experimental setup of thrust force measurement.

    Fig. 6 Interlayer chips and induced damage in drilling CFRP/Ti stacks.

    Fig. 7 Effect of bending stiffness.

    Table 2 Calibration results of specific cutting coefficients.

    To quantitatively analyze the influence of experimental conditions on the interface deflects, a damage factor (Id) is identified to evaluate the damage degree as follows:

    where D is the diameter of the surface damage zone and d is the diameter of the hole.

    Fig. 8 Interface damage factor varying with feed rates and specimen lengths.

    Fig. 8 shows the damage factors at various feed rates and specimen lengths.As is seen from the figure,the average value of factor Idincreases when the feed rate rises from 0.025 mm/r to 0.100 mm/rev. Factor Idincreases significantly when the specimen length increases from 80 mm to 100 mm; however,different tendencies emerge when the length increases from 100 mm to 120 mm. Factor Iddecreases when f0=0.025 mm/r, changes slightly in the case of f0=0.05 mm/r, and increases while f0=0.10 mm/r. The reason is explained in Section 4.2.

    4. Results and discussions

    4.1. Verification of force-deformation model

    Fig. 9 Experimental and predicted thrust forces of low-stiffness CFRP/Ti stack.

    Fig. 9 shows the experimental and the predicted thrust forces under the conditions of f0=0.05 mm/r, l=80 mm. Four specific points and five stages are marked on the force curves.As can be seen in the figure, the two results fit well. The peak value of the thrust force appears at point C and point D,respectively, in the drilling of rigid stacks (Fig. 7) and lowstiffness stacks. The variation of the feed rate is the main reason for this difference. Eq. (2) indicates that the thrust force relates to the integration interval[r1,r2]and the feed rate when the material of the stacks and the geometry of the drill bit are the same. Interval [r1, r2] of the second layer increases during stages II and III, remains constant in stage IV, and decreases in stage V.If the feed rate is invariable,the thrust force reaches its maximum value at point C, and then slightly reduces in stage IV, causing local deformation around the drill bit in the feed direction. However, for low-stiffness stacks, the feed rate is smaller than its set value due to the bending movement in stage IV, and is higher due to the resilience movement in stage V(Eq.(14)).Since the stack recovers within a short time,the effect of the interval[r1,r2]reduction is higher than that of the feed rate increase. Thus,the peak value of the thrust force appears at point D.

    Deviations of the experiment from the prediction exist around point D and in stage II. The first deviation is caused by the local deformation in the feed direction which brings forward springback moment D. Solving Eq. (12) shows that radius distance rspis smaller than chisel radius Rch, revealing the occurrence of the separation during the short time when the chisel edge cuts into the second layer. Theoretically, the chisel edge cuts into the second layer immediately after penetrating the first layer. In the zero length case, the decline of the force in stage II indicates the existence of a protrusion on the lower surface of the first layer (Fig. 7) caused by the plastic deformation around the drill bits.Practically,however,the chisel edge touches the second layer after breaking through the protrusion, causing a time delay between the two actions.It brings forward separation moment B (Fig. 9). Therefore,the reason for the second deviation is the small protrusion at the interface.

    According to Fig. 4, the calibrated coefficients in Table 2 are used to calculate the thrust force and the interlayer gap of low-stiffness CFRP/Ti stacks. The predicted thrust forces are compared with the experimental data in Fig. 10. The two groups of forces have the same trend and the same peak magnitude in all the cases.The predicted error of the peak value is within 8%,indicating a good agreement between the predicted and the experimental thrust forces. The force-deformation coupling model is effective in predicting the thrust force.Since deformation has a strong interaction effect with the thrust force,the proposed coupling model is also valid for predicting the deflection of the stack.

    4.2. Interlayer gap and its effect on interlayer damage

    Since the hole is drilled at the geometric center of the CFRP/Ti stack, the height of the interlayer gap is the largest at the drilling location. Fig. 11 shows the maximum height of the gap varying with the machining time. The gap forms during stage III, and reaches the maximum at the moment just before the drill cuts out of the stack instead of the moment when the drill edge fully cuts into the second layer. With the decrease of bending stiffness, both the maximum value and the appearing time of the interlayer gap increase considerably,while the peak value of the thrust force slightly reduces. Therefore, a higher value of the interlayer gap will be obtained when using the thrust force measured under a well supported condition.

    In the drilling process, continuous chips obliquely move upward between the wall of the hole and the drill flutes.When the interlayer gap appears in stages III and IV,the chips of the titanium alloy crash into the hole edge and the low surface of the CFRP layer,subsequently divided into two parts.One part moves along the drill flutes to be removed, and the other extends vortically along the gap scratching the lower surface of the CFRP layer. With the increase of the interlayer gap,the scratch area increases, the chips of the titanium alloy silt more, and the damage around the CFRP surface deepens(see Fig. 7).

    The height of the interlayer gap reaches its maximum value at the‘‘springback”moment(Fig.9).Sections of the interlayer gap at that moment are shown in Fig. 12, when feed rate f0equals 0.05 mm/rev. Since the interface damage is caused by the movement of the metal chips along the interlayer gap,the damage factor should be proportional to the gap size if the titanium chips are sufficient.At a high feed rate,the experimental results agree with our analysis. However, the damage factor decreases at a large gap size when the feed rate is low(Fig. 8). In the case of a large gap size, the metal chips should be strong enough to exert necessary destructive power on the CFRP layer.Considering the influence of both the attenuation of the destructive power and the flow direction of the metal chips, the effective damage zone of the metal chips is similar to an ellipse. Since the destructive power of the metal chip increases with its thickness, the size of the similar ellipse zone increases with the feed rate. Therefore, the similar ellipse is small at a low feed rate. The damage area of the CFRP layer is within the overlap region of the CFRP surface and the similar ellipse,whose location rises in the similar ellipse zone with the increase of the interlayer gap.In the case of a low feed rate and a large gap, the interlayer gap will be higher than the widest location of the similar ellipse, thus the damage factor decreases (Fig. 12). While under the condition of a large gap and a high feed rate, the interlayer gap is still lower than the widest location of the similar ellipse in all cases. The damage time of the chip increases near the hole edge, while decreases in the outer range in the case of a large gap.This is the reason why the damage factors show different tendencies when the length increases from 100 mm to 120 mm (Fig. 8). Nevertheless, the damage of the CFRP surface deepens though its area lessens,because the damage time of the chip is longer at a large interface gap (Fig. 7).

    Fig. 11 Maximum height of interlayer gap varying with drilling parameters.

    Fig. 12 Interlayer gap and schematic damage zone at different lengths.

    Though the damage factors show different tendencies with the interlayer gap, the interface damage is still significantly affected by the interlayer gap in the drilling of low-stiffness stacks.5,6Results of the correlation analysis indicate a linear dependence relation between damage factor Idand the interlayer gap as shown in Fig. 13, hence, reducing the interlayer gap is an effective way to control the interface defect.

    Fig. 13 Linear regression curve for interface damage factor and interlayer gap.

    4.3. Control of interface damage by clamping foot

    Effects of the cutting parameters and the continuous chips on the interlayer damage will be attenuated when the interlayer gap is sufficiently small.The gap could be reduced by selecting a low feed rate.However,this method leads to production efficiency loss. Though increasing bending stiffness could also reduce the interlayer gap, the structure of the fixture needs to be redesigned,which is both time-consuming and expensive in aeronautical manufacturing fields. Another method widely used in robotic drilling systems is to exert pressure upon the surface of the stack through a nosepiece, as shown in Fig. 14. An annular area is a typical contact area between the stack and the nosepiece, whose diameter of the center line(d)is 10 mm.The pressure value should be within an appropriate range, below which the interface defect cannot be eliminated, while beyond this range, the deformation would be too large to result in a bell hole.

    Since diameter d is small compared to the length of the stack, pressure Fp is assumed to be a concentrated force. To ensure the stability of the drilling process, the nosepiece holds still after the pressure rises to the set value. Then, predeflection δ0caused by Fpcan be expressed as

    To eliminate the interlayer gap, the pre-deflection needs to be larger than the deflection of the second layer in the whole drilling process. Inequality Eq. (20) should be met:

    Critical pressure Fp0is defined as the minimum pressure,beyond which no gap exists.

    When l=80 mm and f0=0.05 mm/rev,the maximum gap is 0.56 mm,and the critical pressure is 1142 N.Three pressures below the critical value are selected to investigate the influence of pressures on the drilling process. The clamping force is exerted through bolts, and measured by a round gasket force sensor (Kistler 9133B). Fig. 15 shows the thrust force and the interlayer damage at different pressures 143, 408, and 1081 N. Pre-deflections for each pressure are 0.07, 0.20, and 0.53 mm, respectively. In stage IV, about 3~5 s of Fig. 15(a),the integration interval [r1, r2] identically equals to [0, R].The thrust force varies with the feed rate according to Eq.(2). When the pressure is close to or higher than the critical value,Eq.(14)indicates that the feed rate is a constant,which means the thrust force is also a constant as shown in Line 4 of Fig. 15(a). However, when the pressure is lower than the critical value,the feed rate is firstly smaller than its set value,then increases to the set value. Hence, the thrust force increases as shown in Lines 1, 2, and 3 of Fig. 15(a).

    Fig. 14 Pressure acting on the stack through a clamping foot.

    Fig. 15 Effect of foot pressure.

    As seen from Fig. 15, the quality of the interlayer is improved with the increase of pressure,because both the magnitude and the lifespan of the interlayer gap reduce. A significant reduction in the interface damage is observed under a very small pressure (143N). In addition to the reduction of the gap size,another reason is that the clamping foot prevents the first layer from bending upward, caused by the metal chips at the interface. At pressure 1081N, which is close to the critical value, the interlayer gap still exists. However, no interface defect is observed in this case,because the gap size is too small for the titanium alloy chips to crash into.

    5. Conclusions

    For the drilling process of low-stiffness CFRP/Ti stacks, this paper investigated the force-deformation coupling law, the formation and the control of the interface defects.It is possible to draw the following conclusions from the present work:

    (1) A coupling model has been proposed to describe the interaction between the thrust force and the deformation of the layers in the drilling of low-stiffness stacks, in which specific cutting coefficients are calibrated using only the drilling force of rigid stacks. The predicted and measured thrust forces have the same trend, and the predicted error of the peak value is within 8%. This proves the validness of the coupling model.

    (2) The drilling process of the low-stiffness stack is divided into five stages based on four specific moments. The interlayer gap forms during stage III, reaching its maximum value at the moment just before the drill cuts out of the stack, instead of the moment when the drill edge fully cuts into the second layer.Both the peak value and the lifespan of the interlayer gap increase with the reduction of bending stiffness.

    (3) The surface damage of the CFRP layer and the interlayer chips are the main types of interlayer defects during the drilling of low-stiffness CFRP/Ti stacks. The damage is caused by rotary extensions of the metal chips along the interlayer gap. Damage factor Idis identified to evaluate the irreparable surface damage. Results of the correlation analysis indicate that factor Idhas a linear dependence relation with the interlayer gap(R2=0.938).

    (4) Damage factor Idincreases with the interlayer gap at a high feed rate, and decreases at a low feed rate, when the gap size is large. The effective damage zone of the metal chips is similar to an ellipse whose size increases with the feed rate.At a low feed rate,the damage factor decreases because the maximum interlayer gap would be higher than the widest location of the similar ellipse.

    (5) Imposing pressure on the surface of the stack by a clamping foot proves effective in reducing the interface damage, because the pressure can reduce both the magnitude and the lifespan of the interlayer gap. Then, a method is presented to determine an appropriate pressure for the elimination of the interlayer damage with the minimum deflection of the stacks.

    In future work,a novel finite element method will be developed to analyze detailed deformation of complex structures in the drilling process. Then, both the physical damage and the geometrical error can be modeled for the accurate assembly of complex thin-wall structures.

    Acknowledgments

    This study was co-supported by the National Natural Science Foundation of China (Nos. 51705426 and 51475379) and the Key Research and Development Plan of Shaanxi Province of China (No. 2017GY-101).

    国语对白做爰xxxⅹ性视频网站| av福利片在线| 亚洲一级一片aⅴ在线观看| 免费观看无遮挡的男女| 韩国高清视频一区二区三区| 亚洲三级黄色毛片| 中文字幕人妻丝袜制服| 久久久久视频综合| av不卡在线播放| 晚上一个人看的免费电影| 国产又色又爽无遮挡免| 亚洲精品aⅴ在线观看| 国产97色在线日韩免费| 亚洲精品国产av成人精品| 超碰97精品在线观看| 久久久久精品性色| 欧美激情 高清一区二区三区| 高清欧美精品videossex| 中国国产av一级| 日韩欧美精品免费久久| 国产精品久久久久久精品电影小说| 9色porny在线观看| 国产极品天堂在线| 亚洲,欧美,日韩| 欧美97在线视频| 一边摸一边做爽爽视频免费| 国产片内射在线| 色吧在线观看| 一本—道久久a久久精品蜜桃钙片| 午夜日本视频在线| 成人毛片a级毛片在线播放| 看免费成人av毛片| 国产熟女欧美一区二区| 午夜日本视频在线| 婷婷成人精品国产| 日韩不卡一区二区三区视频在线| 亚洲国产av新网站| 有码 亚洲区| 日本av手机在线免费观看| 最近手机中文字幕大全| 亚洲欧美精品综合一区二区三区 | 美女国产视频在线观看| 亚洲男人天堂网一区| 母亲3免费完整高清在线观看 | 日韩人妻精品一区2区三区| 亚洲男人天堂网一区| 亚洲精品久久成人aⅴ小说| 国产亚洲精品第一综合不卡| 亚洲精品第二区| 一区二区三区四区激情视频| 久久午夜综合久久蜜桃| 国产精品女同一区二区软件| 精品卡一卡二卡四卡免费| 黄色毛片三级朝国网站| 亚洲精品日本国产第一区| 校园人妻丝袜中文字幕| 国产精品99久久99久久久不卡 | 一级黄片播放器| 久久精品国产鲁丝片午夜精品| 国产在线一区二区三区精| 久久毛片免费看一区二区三区| 国产成人精品在线电影| 永久免费av网站大全| 亚洲成人av在线免费| 精品少妇黑人巨大在线播放| 久久鲁丝午夜福利片| 国产成人一区二区在线| 99re6热这里在线精品视频| videosex国产| 日本-黄色视频高清免费观看| 国产精品国产三级专区第一集| 精品一区二区免费观看| av国产精品久久久久影院| av不卡在线播放| 精品一区二区三卡| 美女中出高潮动态图| 高清视频免费观看一区二区| 久久久亚洲精品成人影院| 1024香蕉在线观看| 免费观看在线日韩| 免费人妻精品一区二区三区视频| 久久久久久伊人网av| 亚洲精品久久午夜乱码| 日韩中文字幕欧美一区二区 | 成人国语在线视频| 黄片播放在线免费| 不卡av一区二区三区| 纯流量卡能插随身wifi吗| 午夜福利视频在线观看免费| 日韩制服丝袜自拍偷拍| 欧美成人精品欧美一级黄| 老熟女久久久| 两个人免费观看高清视频| a级毛片在线看网站| 久久久精品94久久精品| 国产精品久久久久久久久免| 国产有黄有色有爽视频| 多毛熟女@视频| 亚洲精品第二区| 人妻一区二区av| 日本av免费视频播放| 精品少妇久久久久久888优播| 亚洲国产成人一精品久久久| 中文字幕人妻丝袜制服| 国产日韩一区二区三区精品不卡| 亚洲伊人色综图| 国产高清不卡午夜福利| 久久久精品94久久精品| 自线自在国产av| 久久女婷五月综合色啪小说| 日本猛色少妇xxxxx猛交久久| 亚洲欧美一区二区三区国产| 欧美激情极品国产一区二区三区| 亚洲精品乱久久久久久| 亚洲欧美日韩另类电影网站| 人妻一区二区av| 欧美日韩精品成人综合77777| 99久久综合免费| 国产野战对白在线观看| 日韩av不卡免费在线播放| 搡老乐熟女国产| 亚洲国产欧美日韩在线播放| 熟女电影av网| 久久热在线av| 精品99又大又爽又粗少妇毛片| 夫妻性生交免费视频一级片| 国产激情久久老熟女| 国产乱来视频区| 国产探花极品一区二区| 十分钟在线观看高清视频www| 9191精品国产免费久久| 啦啦啦啦在线视频资源| 亚洲美女搞黄在线观看| 秋霞在线观看毛片| 婷婷成人精品国产| 午夜免费观看性视频| 久久人妻熟女aⅴ| 亚洲第一青青草原| 国产男女超爽视频在线观看| 国产熟女午夜一区二区三区| 91在线精品国自产拍蜜月| 国产欧美日韩一区二区三区在线| 久久这里有精品视频免费| 国产在线视频一区二区| 成年美女黄网站色视频大全免费| 成人毛片a级毛片在线播放| 大码成人一级视频| 99久久综合免费| 国产爽快片一区二区三区| 国产在线视频一区二区| 高清视频免费观看一区二区| 精品国产一区二区三区久久久樱花| 伊人亚洲综合成人网| 亚洲色图 男人天堂 中文字幕| 亚洲美女视频黄频| 精品国产乱码久久久久久小说| 熟女少妇亚洲综合色aaa.| 久久青草综合色| 色婷婷久久久亚洲欧美| 欧美黄色片欧美黄色片| 最近最新中文字幕免费大全7| 9热在线视频观看99| 日韩在线高清观看一区二区三区| 亚洲国产av影院在线观看| 国产精品熟女久久久久浪| 国产xxxxx性猛交| 免费观看在线日韩| 精品国产超薄肉色丝袜足j| 国产成人欧美| 亚洲国产成人一精品久久久| 午夜影院在线不卡| 久久女婷五月综合色啪小说| 观看av在线不卡| 亚洲情色 制服丝袜| 国产日韩欧美视频二区| 宅男免费午夜| 午夜日本视频在线| 自拍欧美九色日韩亚洲蝌蚪91| 王馨瑶露胸无遮挡在线观看| 黄色视频在线播放观看不卡| 丰满饥渴人妻一区二区三| 国产爽快片一区二区三区| 日韩在线高清观看一区二区三区| 天堂中文最新版在线下载| 国产精品无大码| 一二三四中文在线观看免费高清| 男女无遮挡免费网站观看| 麻豆av在线久日| 最近中文字幕高清免费大全6| 精品人妻一区二区三区麻豆| 大香蕉久久成人网| 亚洲四区av| 高清视频免费观看一区二区| 日本猛色少妇xxxxx猛交久久| 久热久热在线精品观看| av网站免费在线观看视频| 校园人妻丝袜中文字幕| 欧美精品高潮呻吟av久久| 这个男人来自地球电影免费观看 | 精品卡一卡二卡四卡免费| 国产亚洲av片在线观看秒播厂| 国产精品 国内视频| av在线播放精品| 妹子高潮喷水视频| 久久人人97超碰香蕉20202| 国产极品天堂在线| 最新的欧美精品一区二区| 中国三级夫妇交换| 久久国产亚洲av麻豆专区| 啦啦啦啦在线视频资源| 一区二区日韩欧美中文字幕| 久久女婷五月综合色啪小说| 久久av网站| 欧美日本中文国产一区发布| 91午夜精品亚洲一区二区三区| 一级毛片 在线播放| 亚洲欧美日韩另类电影网站| 日本-黄色视频高清免费观看| 两个人看的免费小视频| 一区福利在线观看| 一级毛片 在线播放| 在线观看www视频免费| 天堂中文最新版在线下载| 精品少妇久久久久久888优播| 国产色婷婷99| 极品人妻少妇av视频| 亚洲精品第二区| 国产精品av久久久久免费| 两个人看的免费小视频| 久久久精品免费免费高清| 最新中文字幕久久久久| 黄网站色视频无遮挡免费观看| 老司机影院毛片| 91久久精品国产一区二区三区| www.自偷自拍.com| 亚洲熟女精品中文字幕| 久久人人爽人人片av| 黑丝袜美女国产一区| 国产亚洲最大av| 日日啪夜夜爽| 最黄视频免费看| 九色亚洲精品在线播放| 欧美激情高清一区二区三区 | 王馨瑶露胸无遮挡在线观看| 久久久久国产网址| 在线观看国产h片| 久久久国产欧美日韩av| 搡女人真爽免费视频火全软件| 久久久久国产精品人妻一区二区| 久久99一区二区三区| 一二三四中文在线观看免费高清| 99国产精品免费福利视频| 熟女少妇亚洲综合色aaa.| 丝瓜视频免费看黄片| 18在线观看网站| 性色avwww在线观看| 国产成人一区二区在线| 91久久精品国产一区二区三区| videos熟女内射| 91国产中文字幕| 亚洲一区中文字幕在线| 国产熟女欧美一区二区| videossex国产| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 三级国产精品片| 一级毛片黄色毛片免费观看视频| 国产1区2区3区精品| 国产午夜精品一二区理论片| 熟女少妇亚洲综合色aaa.| 伊人久久国产一区二区| 少妇被粗大猛烈的视频| 国产精品久久久久久精品古装| 最近最新中文字幕大全免费视频 | 黑人欧美特级aaaaaa片| 国产av一区二区精品久久| 国产成人91sexporn| 建设人人有责人人尽责人人享有的| 男女免费视频国产| av电影中文网址| 国产成人精品福利久久| 观看美女的网站| 久久婷婷青草| 又粗又硬又长又爽又黄的视频| 男人舔女人的私密视频| 亚洲天堂av无毛| 侵犯人妻中文字幕一二三四区| av.在线天堂| 国产精品久久久久久av不卡| av有码第一页| av国产精品久久久久影院| 18+在线观看网站| 国产亚洲欧美精品永久| 日日啪夜夜爽| 精品久久久久久电影网| 999久久久国产精品视频| 国产精品偷伦视频观看了| 曰老女人黄片| 国产福利在线免费观看视频| 欧美 亚洲 国产 日韩一| 成年人免费黄色播放视频| 伦理电影免费视频| 久久久欧美国产精品| av又黄又爽大尺度在线免费看| 亚洲伊人色综图| 日韩中字成人| 欧美成人精品欧美一级黄| 只有这里有精品99| 一级黄片播放器| 天天操日日干夜夜撸| 天天躁日日躁夜夜躁夜夜| 久久久久国产网址| 高清不卡的av网站| 亚洲男人天堂网一区| 免费av中文字幕在线| 国产免费又黄又爽又色| 精品久久蜜臀av无| 久久国产精品男人的天堂亚洲| 日本vs欧美在线观看视频| 亚洲国产成人一精品久久久| 少妇被粗大猛烈的视频| 18禁观看日本| 少妇人妻 视频| 女人被躁到高潮嗷嗷叫费观| 最黄视频免费看| 亚洲精品中文字幕在线视频| 日本色播在线视频| 婷婷成人精品国产| 99香蕉大伊视频| 婷婷成人精品国产| 日本午夜av视频| 妹子高潮喷水视频| 日本色播在线视频| 国产欧美日韩一区二区三区在线| 免费在线观看黄色视频的| 国产一级毛片在线| 精品国产乱码久久久久久小说| 亚洲精品中文字幕在线视频| 69精品国产乱码久久久| 女人被躁到高潮嗷嗷叫费观| 最新中文字幕久久久久| 丝瓜视频免费看黄片| 99热国产这里只有精品6| 热re99久久精品国产66热6| 丝袜人妻中文字幕| 国产日韩欧美视频二区| 久久99蜜桃精品久久| 久久久久国产精品人妻一区二区| 丝袜人妻中文字幕| 丰满饥渴人妻一区二区三| 欧美日韩成人在线一区二区| 天堂俺去俺来也www色官网| 男女边摸边吃奶| 丝袜脚勾引网站| 午夜福利网站1000一区二区三区| 国产片特级美女逼逼视频| 免费观看av网站的网址| 99久国产av精品国产电影| 99精国产麻豆久久婷婷| 99久久精品国产国产毛片| 最近中文字幕2019免费版| 99九九在线精品视频| 免费高清在线观看日韩| 婷婷色综合www| 国产熟女欧美一区二区| 热99国产精品久久久久久7| 一级毛片黄色毛片免费观看视频| 大香蕉久久网| 精品午夜福利在线看| 女人高潮潮喷娇喘18禁视频| 少妇的丰满在线观看| 国产亚洲av片在线观看秒播厂| 桃花免费在线播放| 亚洲国产毛片av蜜桃av| 日韩av免费高清视频| 精品一品国产午夜福利视频| 97在线视频观看| 久久久精品免费免费高清| 男女高潮啪啪啪动态图| 国产精品香港三级国产av潘金莲 | 亚洲美女视频黄频| 在线观看人妻少妇| 久久久国产一区二区| 久久精品熟女亚洲av麻豆精品| 久久久久国产一级毛片高清牌| 亚洲国产毛片av蜜桃av| 一区在线观看完整版| 日韩 亚洲 欧美在线| 考比视频在线观看| 欧美日韩成人在线一区二区| 91久久精品国产一区二区三区| 人妻系列 视频| 中文字幕另类日韩欧美亚洲嫩草| 国产亚洲最大av| 欧美成人午夜精品| 国产探花极品一区二区| 18在线观看网站| av在线播放精品| 亚洲国产毛片av蜜桃av| 久久99蜜桃精品久久| 精品卡一卡二卡四卡免费| 国产精品麻豆人妻色哟哟久久| 午夜日本视频在线| 久久综合国产亚洲精品| 自线自在国产av| 成人毛片a级毛片在线播放| 激情五月婷婷亚洲| 老汉色∧v一级毛片| 永久网站在线| 久久99热这里只频精品6学生| 国产又爽黄色视频| 多毛熟女@视频| 国产精品不卡视频一区二区| 91在线精品国自产拍蜜月| 国产精品.久久久| 欧美日韩国产mv在线观看视频| 丝袜美腿诱惑在线| 国产 精品1| 日韩av不卡免费在线播放| 人妻系列 视频| 香蕉国产在线看| 国产成人欧美| 99热网站在线观看| 欧美在线黄色| 青春草亚洲视频在线观看| 国产探花极品一区二区| videossex国产| 午夜福利在线观看免费完整高清在| 美女国产高潮福利片在线看| 亚洲人成电影观看| 18禁观看日本| 国产一区亚洲一区在线观看| 午夜免费男女啪啪视频观看| 久久久久久久精品精品| 亚洲av电影在线进入| 99久久精品国产国产毛片| 中文天堂在线官网| 国产精品国产三级国产专区5o| 亚洲天堂av无毛| av线在线观看网站| 国产亚洲欧美精品永久| 免费在线观看完整版高清| 观看美女的网站| 久久久国产精品麻豆| 咕卡用的链子| 一级片免费观看大全| 男女下面插进去视频免费观看| 男人爽女人下面视频在线观看| 1024视频免费在线观看| 亚洲精品中文字幕在线视频| 又大又黄又爽视频免费| 精品亚洲成国产av| 自拍欧美九色日韩亚洲蝌蚪91| 欧美变态另类bdsm刘玥| 观看美女的网站| 人人妻人人澡人人爽人人夜夜| 国语对白做爰xxxⅹ性视频网站| 成人毛片60女人毛片免费| 国产日韩欧美视频二区| 99香蕉大伊视频| 精品亚洲成a人片在线观看| 一边亲一边摸免费视频| 欧美激情 高清一区二区三区| 中文字幕av电影在线播放| 人体艺术视频欧美日本| 春色校园在线视频观看| 国产成人精品在线电影| 国产熟女午夜一区二区三区| 激情五月婷婷亚洲| 欧美97在线视频| av又黄又爽大尺度在线免费看| 女人高潮潮喷娇喘18禁视频| 国产极品粉嫩免费观看在线| 国产成人欧美| 欧美日韩精品成人综合77777| 飞空精品影院首页| 中文欧美无线码| 久热久热在线精品观看| 如日韩欧美国产精品一区二区三区| 精品国产超薄肉色丝袜足j| 超色免费av| 99热全是精品| 国产精品成人在线| 一边摸一边做爽爽视频免费| 国产极品天堂在线| 亚洲三区欧美一区| 久久久久久久久久久久大奶| 国产在视频线精品| 国产高清不卡午夜福利| xxxhd国产人妻xxx| 欧美黄色片欧美黄色片| 国产精品欧美亚洲77777| 久久人人97超碰香蕉20202| 99香蕉大伊视频| 欧美成人精品欧美一级黄| 青草久久国产| 在线免费观看不下载黄p国产| 高清黄色对白视频在线免费看| 亚洲av免费高清在线观看| 日本-黄色视频高清免费观看| 国产在视频线精品| 国产av码专区亚洲av| videossex国产| 久久久久久久久久人人人人人人| 中文字幕av电影在线播放| 黄色一级大片看看| 免费观看性生交大片5| 99re6热这里在线精品视频| 国产综合精华液| 亚洲在久久综合| 黄色怎么调成土黄色| 久久久久国产网址| 免费观看无遮挡的男女| 久久久国产精品麻豆| 丝袜美足系列| 黑丝袜美女国产一区| 满18在线观看网站| 99九九在线精品视频| 五月开心婷婷网| 亚洲综合精品二区| 亚洲av日韩在线播放| 亚洲欧美日韩另类电影网站| 日韩av在线免费看完整版不卡| 亚洲五月色婷婷综合| 亚洲精品久久久久久婷婷小说| 成人二区视频| 久久久精品国产亚洲av高清涩受| 国产日韩欧美亚洲二区| xxx大片免费视频| 日韩欧美一区视频在线观看| 国产亚洲欧美精品永久| 午夜福利乱码中文字幕| 亚洲欧洲日产国产| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产亚洲午夜精品一区二区久久| 国产精品国产av在线观看| 在线观看美女被高潮喷水网站| 国产成人精品在线电影| 色婷婷久久久亚洲欧美| 国产精品国产三级专区第一集| 天堂俺去俺来也www色官网| 免费播放大片免费观看视频在线观看| 91成人精品电影| 爱豆传媒免费全集在线观看| 久久精品熟女亚洲av麻豆精品| 搡老乐熟女国产| 亚洲国产日韩一区二区| 日日摸夜夜添夜夜爱| 成年人免费黄色播放视频| 91精品国产国语对白视频| 叶爱在线成人免费视频播放| 99久国产av精品国产电影| 1024视频免费在线观看| 成人二区视频| 最新中文字幕久久久久| 久久精品国产鲁丝片午夜精品| 一本色道久久久久久精品综合| 亚洲第一av免费看| 大片电影免费在线观看免费| 在线观看一区二区三区激情| 啦啦啦在线观看免费高清www| 亚洲三级黄色毛片| 男女高潮啪啪啪动态图| 久久毛片免费看一区二区三区| 精品第一国产精品| 男女无遮挡免费网站观看| 成人国产麻豆网| 久久久久久久国产电影| 久久久精品国产亚洲av高清涩受| 天堂俺去俺来也www色官网| 永久免费av网站大全| 侵犯人妻中文字幕一二三四区| 国产黄频视频在线观看| 18禁动态无遮挡网站| 菩萨蛮人人尽说江南好唐韦庄| 国产一区有黄有色的免费视频| 亚洲欧洲日产国产| 国产成人一区二区在线| 91精品国产国语对白视频| 熟女电影av网| 9热在线视频观看99| 免费人妻精品一区二区三区视频| 久热久热在线精品观看| 美女国产视频在线观看| 日韩精品免费视频一区二区三区| 美女福利国产在线| 1024视频免费在线观看| av在线老鸭窝| 高清不卡的av网站| 视频在线观看一区二区三区| 免费黄频网站在线观看国产| 人人妻人人添人人爽欧美一区卜| 欧美精品一区二区大全| 日日啪夜夜爽| 亚洲精品久久午夜乱码| 久久久亚洲精品成人影院| 精品亚洲成国产av| 久久97久久精品| 亚洲av国产av综合av卡| 午夜影院在线不卡| 在线观看免费高清a一片| 只有这里有精品99| 美女中出高潮动态图| 亚洲国产欧美网| av在线观看视频网站免费| 国产色婷婷99| 好男人视频免费观看在线| 亚洲美女黄色视频免费看| 一区福利在线观看| 性色avwww在线观看| 欧美日韩精品成人综合77777| 亚洲人成网站在线观看播放| 中文字幕精品免费在线观看视频| 99热国产这里只有精品6| 久久亚洲国产成人精品v| a级毛片黄视频|