• <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).

    色吧在线观看| 十八禁网站网址无遮挡| 不卡视频在线观看欧美| 亚洲av欧美aⅴ国产| 日本vs欧美在线观看视频| 亚洲情色 制服丝袜| 国产黄色视频一区二区在线观看| 日本91视频免费播放| 国产成人精品福利久久| 欧美精品一区二区免费开放| 色视频在线一区二区三区| 亚洲欧美精品综合一区二区三区| h视频一区二区三区| 色94色欧美一区二区| videos熟女内射| 大码成人一级视频| 欧美日韩福利视频一区二区| 成人毛片60女人毛片免费| 国产精品av久久久久免费| 大香蕉久久网| 五月开心婷婷网| 国产一区二区激情短视频 | 国产精品一区二区在线不卡| 欧美久久黑人一区二区| 人妻人人澡人人爽人人| 国产精品国产三级国产专区5o| 天堂俺去俺来也www色官网| 一边摸一边做爽爽视频免费| 亚洲精品成人av观看孕妇| 亚洲一级一片aⅴ在线观看| 日本一区二区免费在线视频| 国产精品免费大片| 亚洲欧洲精品一区二区精品久久久 | 一区二区三区四区激情视频| 亚洲国产最新在线播放| 女性被躁到高潮视频| 侵犯人妻中文字幕一二三四区| 9热在线视频观看99| 老司机靠b影院| 在线观看人妻少妇| 亚洲,欧美精品.| 90打野战视频偷拍视频| 国产欧美日韩综合在线一区二区| 精品少妇内射三级| 少妇精品久久久久久久| 国产女主播在线喷水免费视频网站| xxx大片免费视频| 王馨瑶露胸无遮挡在线观看| 亚洲成人免费av在线播放| 美女高潮到喷水免费观看| 一区在线观看完整版| 亚洲av电影在线观看一区二区三区| 老司机在亚洲福利影院| 久久久精品免费免费高清| 午夜福利网站1000一区二区三区| 亚洲综合色网址| 免费在线观看视频国产中文字幕亚洲 | 国产精品一区二区在线不卡| 欧美 日韩 精品 国产| 亚洲熟女精品中文字幕| 免费高清在线观看日韩| 色吧在线观看| 老司机影院成人| 国产xxxxx性猛交| 亚洲欧美精品自产自拍| 国产精品 欧美亚洲| 777久久人妻少妇嫩草av网站| 男人添女人高潮全过程视频| 在线看a的网站| 最新在线观看一区二区三区 | 精品久久久精品久久久| 97在线人人人人妻| 最近最新中文字幕大全免费视频 | 欧美人与性动交α欧美软件| 丝瓜视频免费看黄片| 成人亚洲精品一区在线观看| 色吧在线观看| 国产片特级美女逼逼视频| 天天躁狠狠躁夜夜躁狠狠躁| 午夜福利网站1000一区二区三区| 欧美久久黑人一区二区| 熟妇人妻不卡中文字幕| 亚洲av日韩在线播放| 亚洲国产欧美在线一区| 久久精品久久久久久久性| 99国产精品免费福利视频| 青草久久国产| 欧美久久黑人一区二区| 亚洲国产成人一精品久久久| 波野结衣二区三区在线| 欧美中文综合在线视频| 啦啦啦啦在线视频资源| 亚洲男人天堂网一区| 两性夫妻黄色片| 人妻人人澡人人爽人人| 在线天堂中文资源库| 国产精品二区激情视频| 五月开心婷婷网| 亚洲国产精品一区二区三区在线| 成人毛片60女人毛片免费| 亚洲欧美成人综合另类久久久| 2018国产大陆天天弄谢| 卡戴珊不雅视频在线播放| 欧美日韩av久久| 黑人巨大精品欧美一区二区蜜桃| 尾随美女入室| 亚洲欧洲国产日韩| 又黄又粗又硬又大视频| 老熟女久久久| 人人妻人人爽人人添夜夜欢视频| 欧美 日韩 精品 国产| 大片免费播放器 马上看| 韩国av在线不卡| 欧美 日韩 精品 国产| 亚洲少妇的诱惑av| 80岁老熟妇乱子伦牲交| 欧美日韩亚洲综合一区二区三区_| 国产黄频视频在线观看| 在线观看三级黄色| 亚洲精品av麻豆狂野| 天天躁日日躁夜夜躁夜夜| 久久久国产欧美日韩av| 大香蕉久久成人网| 一本久久精品| 国产成人一区二区在线| 亚洲欧美成人综合另类久久久| 国产在线一区二区三区精| 欧美日韩av久久| 天堂俺去俺来也www色官网| 亚洲成国产人片在线观看| 精品一区二区免费观看| 少妇人妻久久综合中文| 国产片内射在线| 王馨瑶露胸无遮挡在线观看| 欧美成人午夜精品| 美女午夜性视频免费| 性高湖久久久久久久久免费观看| 亚洲人成电影观看| 亚洲成国产人片在线观看| 老司机在亚洲福利影院| 在线免费观看不下载黄p国产| 亚洲国产欧美在线一区| 亚洲成人手机| 久久精品久久久久久久性| 国产伦理片在线播放av一区| 亚洲精品aⅴ在线观看| 精品一区二区免费观看| 国产成人免费无遮挡视频| 婷婷色av中文字幕| 1024香蕉在线观看| 九九爱精品视频在线观看| 十八禁网站网址无遮挡| 亚洲激情五月婷婷啪啪| 国产精品秋霞免费鲁丝片| 国产精品三级大全| av女优亚洲男人天堂| 秋霞伦理黄片| 欧美日韩亚洲综合一区二区三区_| 极品人妻少妇av视频| 啦啦啦中文免费视频观看日本| 亚洲一区中文字幕在线| 国产爽快片一区二区三区| 国产一卡二卡三卡精品 | 国产精品成人在线| 亚洲精品乱久久久久久| 热re99久久国产66热| 精品国产超薄肉色丝袜足j| 精品视频人人做人人爽| 高清黄色对白视频在线免费看| videos熟女内射| 黄色视频不卡| 日韩大码丰满熟妇| 国产精品 欧美亚洲| 亚洲欧美一区二区三区久久| 国产一级毛片在线| 日韩视频在线欧美| 国产成人精品在线电影| 久久久亚洲精品成人影院| 爱豆传媒免费全集在线观看| 久久国产亚洲av麻豆专区| 99久久精品国产亚洲精品| av在线播放精品| 国产日韩欧美视频二区| 欧美黑人精品巨大| 亚洲av电影在线进入| av免费观看日本| 久久99热这里只频精品6学生| 日韩一区二区三区影片| 国产免费又黄又爽又色| 90打野战视频偷拍视频| 亚洲天堂av无毛| 日本av手机在线免费观看| 不卡av一区二区三区| 女人被躁到高潮嗷嗷叫费观| 亚洲久久久国产精品| 黄片无遮挡物在线观看| 国产人伦9x9x在线观看| 国产精品一二三区在线看| kizo精华| 免费高清在线观看视频在线观看| 亚洲精品国产一区二区精华液| 在现免费观看毛片| 国产亚洲av片在线观看秒播厂| av有码第一页| 69精品国产乱码久久久| 国产视频首页在线观看| 中文字幕精品免费在线观看视频| 国产亚洲av片在线观看秒播厂| 国产成人精品福利久久| 日韩大码丰满熟妇| 欧美97在线视频| 国产精品一区二区在线观看99| 国产黄频视频在线观看| 男女边吃奶边做爰视频| 国产精品熟女久久久久浪| 欧美老熟妇乱子伦牲交| 久久久久久免费高清国产稀缺| 亚洲中文av在线| 日韩精品有码人妻一区| 人人妻人人添人人爽欧美一区卜| 美女扒开内裤让男人捅视频| 国产伦人伦偷精品视频| 天天躁夜夜躁狠狠久久av| 国产精品亚洲av一区麻豆 | 欧美在线黄色| 国产伦人伦偷精品视频| 亚洲一卡2卡3卡4卡5卡精品中文| 久久影院123| 高清欧美精品videossex| 另类精品久久| 久久精品aⅴ一区二区三区四区| 亚洲欧美一区二区三区黑人| 亚洲一级一片aⅴ在线观看| 一本大道久久a久久精品| 国产日韩欧美在线精品| 亚洲av电影在线进入| 日本一区二区免费在线视频| 亚洲国产欧美网| 9191精品国产免费久久| 激情五月婷婷亚洲| 亚洲精品第二区| 80岁老熟妇乱子伦牲交| 久久精品亚洲av国产电影网| 多毛熟女@视频| 一边摸一边做爽爽视频免费| 成人午夜精彩视频在线观看| 亚洲欧美精品自产自拍| 亚洲精品国产色婷婷电影| 水蜜桃什么品种好| 中文字幕人妻丝袜一区二区 | 欧美黄色片欧美黄色片| 1024视频免费在线观看| 免费高清在线观看视频在线观看| a级毛片黄视频| 国产免费又黄又爽又色| 丝袜美腿诱惑在线| 国产精品亚洲av一区麻豆 | 国产精品二区激情视频| 制服诱惑二区| 日日摸夜夜添夜夜爱| www.自偷自拍.com| 国产精品av久久久久免费| 亚洲av福利一区| 国产在视频线精品| 777久久人妻少妇嫩草av网站| 美女福利国产在线| 国产精品久久久久久久久免| 成人免费观看视频高清| 亚洲欧美成人综合另类久久久| 一本久久精品| 国产精品国产av在线观看| 欧美人与性动交α欧美精品济南到| 成人18禁高潮啪啪吃奶动态图| 久久久久久人人人人人| 国产成人精品无人区| 观看美女的网站| 中文字幕制服av| 久久精品人人爽人人爽视色| 日韩熟女老妇一区二区性免费视频| 久久天堂一区二区三区四区| 母亲3免费完整高清在线观看| 国产精品亚洲av一区麻豆 | 精品国产露脸久久av麻豆| 亚洲欧洲日产国产| 97人妻天天添夜夜摸| 久久久精品免费免费高清| 最近最新中文字幕大全免费视频 | 国产免费视频播放在线视频| 一边摸一边抽搐一进一出视频| 午夜免费观看性视频| 欧美日韩亚洲综合一区二区三区_| 美女主播在线视频| 最近中文字幕高清免费大全6| 国产男女内射视频| 交换朋友夫妻互换小说| 人人妻人人澡人人看| 一本大道久久a久久精品| 黑人巨大精品欧美一区二区蜜桃| 成人亚洲精品一区在线观看| 日韩一卡2卡3卡4卡2021年| 亚洲精品日本国产第一区| 国产成人系列免费观看| 久久精品国产亚洲av高清一级| 成年女人毛片免费观看观看9 | 免费观看a级毛片全部| 一本色道久久久久久精品综合| 日本av手机在线免费观看| 自线自在国产av| 国产在线免费精品| 日日撸夜夜添| 日本欧美国产在线视频| 亚洲精品一二三| 最近中文字幕2019免费版| 精品国产一区二区久久| 丰满饥渴人妻一区二区三| 日本黄色日本黄色录像| 亚洲天堂av无毛| xxx大片免费视频| 高清不卡的av网站| 精品视频人人做人人爽| 操出白浆在线播放| av不卡在线播放| 亚洲第一av免费看| 亚洲精品,欧美精品| 日本爱情动作片www.在线观看| 国产深夜福利视频在线观看| 欧美最新免费一区二区三区| www.精华液| 久久久久精品性色| 国产探花极品一区二区| 日日撸夜夜添| 女人爽到高潮嗷嗷叫在线视频| 赤兔流量卡办理| 69精品国产乱码久久久| 女人爽到高潮嗷嗷叫在线视频| 欧美激情 高清一区二区三区| 成年女人毛片免费观看观看9 | 久久ye,这里只有精品| 建设人人有责人人尽责人人享有的| av视频免费观看在线观看| 青草久久国产| 精品久久久精品久久久| www.自偷自拍.com| 深夜精品福利| 中文精品一卡2卡3卡4更新| 久久久精品区二区三区| av又黄又爽大尺度在线免费看| 欧美日韩福利视频一区二区| 女人爽到高潮嗷嗷叫在线视频| 欧美激情 高清一区二区三区| av福利片在线| 91精品伊人久久大香线蕉| 国产野战对白在线观看| 久久国产精品男人的天堂亚洲| 亚洲男人天堂网一区| 精品国产一区二区久久| 亚洲自偷自拍图片 自拍| 精品国产一区二区三区四区第35| 成人手机av| 美女福利国产在线| 久久久亚洲精品成人影院| 久久这里只有精品19| 老司机在亚洲福利影院| 久久精品亚洲熟妇少妇任你| 最近的中文字幕免费完整| 欧美激情 高清一区二区三区| 成人国语在线视频| 婷婷色综合www| 日本91视频免费播放| 色精品久久人妻99蜜桃| 我的亚洲天堂| 日本91视频免费播放| 国产精品成人在线| 婷婷色综合www| 人妻一区二区av| 少妇人妻 视频| 亚洲精品第二区| 亚洲七黄色美女视频| 亚洲成人手机| 亚洲精华国产精华液的使用体验| 大码成人一级视频| 久久婷婷青草| 制服人妻中文乱码| 亚洲精品美女久久av网站| 黄色怎么调成土黄色| 国产毛片在线视频| 国产成人av激情在线播放| 在线精品无人区一区二区三| 国产1区2区3区精品| 成人免费观看视频高清| 欧美乱码精品一区二区三区| 亚洲欧美精品自产自拍| 最近的中文字幕免费完整| 亚洲精品日韩在线中文字幕| 妹子高潮喷水视频| 最黄视频免费看| 男人舔女人的私密视频| 亚洲av男天堂| 久久狼人影院| √禁漫天堂资源中文www| 午夜福利乱码中文字幕| 丰满迷人的少妇在线观看| 一本一本久久a久久精品综合妖精| 亚洲av电影在线观看一区二区三区| 尾随美女入室| 日日爽夜夜爽网站| 人体艺术视频欧美日本| av不卡在线播放| 欧美在线一区亚洲| 少妇 在线观看| 国产精品国产三级专区第一集| 免费黄色在线免费观看| 波多野结衣一区麻豆| 天天躁夜夜躁狠狠躁躁| 成人国语在线视频| 亚洲国产看品久久| 国产成人欧美在线观看 | 日韩视频在线欧美| 伦理电影免费视频| 亚洲精品一二三| 天堂8中文在线网| 日本黄色日本黄色录像| av国产久精品久网站免费入址| 美女脱内裤让男人舔精品视频| xxxhd国产人妻xxx| 免费高清在线观看日韩| 中文乱码字字幕精品一区二区三区| 国产伦人伦偷精品视频| 97人妻天天添夜夜摸| 日韩大码丰满熟妇| 亚洲国产成人一精品久久久| 一区二区av电影网| 巨乳人妻的诱惑在线观看| 国产精品 欧美亚洲| 久久亚洲国产成人精品v| 久久ye,这里只有精品| 久久av网站| 黄色毛片三级朝国网站| 中文字幕人妻丝袜制服| 最近的中文字幕免费完整| 国产精品一国产av| 天天躁夜夜躁狠狠躁躁| 国产97色在线日韩免费| 精品视频人人做人人爽| 在线观看国产h片| 亚洲精品美女久久久久99蜜臀 | 精品久久久精品久久久| 丝袜美腿诱惑在线| 99久国产av精品国产电影| 男女下面插进去视频免费观看| 水蜜桃什么品种好| 成人手机av| 欧美日韩亚洲高清精品| 操美女的视频在线观看| 人妻一区二区av| 中国国产av一级| 精品人妻熟女毛片av久久网站| 国产成人免费无遮挡视频| 欧美日韩亚洲综合一区二区三区_| 免费在线观看视频国产中文字幕亚洲 | 成年人午夜在线观看视频| 高清黄色对白视频在线免费看| 久热爱精品视频在线9| 国产精品一国产av| 新久久久久国产一级毛片| 日韩av免费高清视频| 久久久欧美国产精品| 天美传媒精品一区二区| 夫妻性生交免费视频一级片| 老司机影院成人| 丝袜人妻中文字幕| 纯流量卡能插随身wifi吗| 亚洲国产av影院在线观看| av在线老鸭窝| 国产精品免费大片| 久久av网站| 中文字幕高清在线视频| 最近最新中文字幕大全免费视频 | 亚洲国产成人一精品久久久| 国产精品亚洲av一区麻豆 | 亚洲 欧美一区二区三区| 大香蕉久久网| 欧美人与性动交α欧美软件| 男女边摸边吃奶| 婷婷色av中文字幕| 国产成人91sexporn| 久久国产精品男人的天堂亚洲| 激情五月婷婷亚洲| 亚洲国产日韩一区二区| 天天躁日日躁夜夜躁夜夜| 国产成人精品无人区| 国产日韩欧美视频二区| 又大又黄又爽视频免费| 无遮挡黄片免费观看| www.精华液| 国产精品二区激情视频| 99九九在线精品视频| 成年美女黄网站色视频大全免费| 少妇人妻精品综合一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 国产亚洲av片在线观看秒播厂| 国产激情久久老熟女| 尾随美女入室| 丰满少妇做爰视频| 天天躁夜夜躁狠狠久久av| 晚上一个人看的免费电影| 七月丁香在线播放| 成人午夜精彩视频在线观看| 日日爽夜夜爽网站| 欧美最新免费一区二区三区| 高清欧美精品videossex| 久久天躁狠狠躁夜夜2o2o | 国产一区二区三区综合在线观看| 亚洲国产欧美网| 高清黄色对白视频在线免费看| 热re99久久精品国产66热6| 国产国语露脸激情在线看| 亚洲av成人不卡在线观看播放网 | 视频区图区小说| 午夜影院在线不卡| 国产日韩一区二区三区精品不卡| 岛国毛片在线播放| 色婷婷久久久亚洲欧美| 亚洲美女黄色视频免费看| 乱人伦中国视频| 午夜久久久在线观看| 国产探花极品一区二区| 丰满迷人的少妇在线观看| 欧美日韩综合久久久久久| 九色亚洲精品在线播放| 国产在线一区二区三区精| 久久精品人人爽人人爽视色| 一边摸一边抽搐一进一出视频| 男人操女人黄网站| 午夜影院在线不卡| 久久精品久久精品一区二区三区| 国产福利在线免费观看视频| 免费看av在线观看网站| 亚洲av福利一区| 美女午夜性视频免费| 日韩av免费高清视频| 80岁老熟妇乱子伦牲交| 搡老乐熟女国产| 精品亚洲成a人片在线观看| 国产精品人妻久久久影院| 欧美av亚洲av综合av国产av | 午夜精品国产一区二区电影| av福利片在线| 日韩一区二区三区影片| 久久国产精品男人的天堂亚洲| 国产97色在线日韩免费| 精品亚洲乱码少妇综合久久| 亚洲国产欧美一区二区综合| 亚洲三区欧美一区| 欧美日韩一级在线毛片| 欧美成人精品欧美一级黄| 视频在线观看一区二区三区| 青春草国产在线视频| 91国产中文字幕| 久久久精品区二区三区| 观看av在线不卡| 91成人精品电影| 成年女人毛片免费观看观看9 | 精品一区二区三区av网在线观看 | 91精品三级在线观看| 一二三四中文在线观看免费高清| 亚洲成人国产一区在线观看 | 国产免费一区二区三区四区乱码| 大话2 男鬼变身卡| 一级a爱视频在线免费观看| 国产黄频视频在线观看| 亚洲专区中文字幕在线 | 亚洲成人免费av在线播放| 日韩大片免费观看网站| av在线老鸭窝| 国产成人啪精品午夜网站| 人妻 亚洲 视频| 亚洲一级一片aⅴ在线观看| 十八禁高潮呻吟视频| 精品久久蜜臀av无| 观看美女的网站| 深夜精品福利| 亚洲欧美一区二区三区久久| 国产免费福利视频在线观看| 亚洲四区av| 制服丝袜香蕉在线| 一二三四在线观看免费中文在| 免费高清在线观看视频在线观看| 999久久久国产精品视频| 777米奇影视久久| a级片在线免费高清观看视频| 91精品伊人久久大香线蕉| 亚洲精品aⅴ在线观看| 婷婷色av中文字幕| 久久久精品国产亚洲av高清涩受| 亚洲国产看品久久| 9色porny在线观看| 美女脱内裤让男人舔精品视频| 2021少妇久久久久久久久久久| 国产成人免费观看mmmm| 欧美在线黄色| av一本久久久久| 国产午夜精品一二区理论片| 黄片播放在线免费| 久久久国产精品麻豆| 午夜福利影视在线免费观看| 久久久久久久久久久久大奶| 亚洲精品中文字幕在线视频| 精品久久久久久电影网| 国产精品久久久人人做人人爽| 亚洲色图综合在线观看| 国产成人精品无人区| 国产乱来视频区| 少妇人妻久久综合中文| 男女无遮挡免费网站观看|