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

    Chatter stability of robotic rotary ultrasonic countersinking

    2023-11-10 02:16:30ZhenwenSUNWenheLIAOKnZHENGSongDONGPeiLEILinjunSUN
    CHINESE JOURNAL OF AERONAUTICS 2023年10期

    Zhenwen SUN, Wenhe LIAO,*, Kn ZHENG, Song DONG, Pei LEI,Linjun SUN

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

    b Chengdu Aircraft Industrial (Group) Co., Ltd., Chengdu 610091, China

    KEYWORDS

    Abstract In recent years,industrial robots have received extensive attention in manufacturing field due to their high flexibility and great workspace.However, the weak stiffness of industrial robots makes it extremely easy to arouse chatter,which affects machining quality inevitably and generates noise pollution in severe cases.Compared with drilling, the chatter mechanism of robotic countersinking is more complex.The external excitation changes with cutting width and depth in countersinking.This characteristic results in time-varying and nonlinearity of robotic countersinking dynamics.Thus,it is urgent to propose a new method of chatter suppression and provide an accurate stability analysis model.As a new special machining technology, rotary ultrasonic machining has been proved to improve robotic drilling and milling stability effectively.Based on this, robotic rotary ultrasonic countersinking(RRUC)is proposed to improve the robotic countersinking stability in this paper.A three-dimensional stability domain method of RRUC is established.First, the countersinking process was divided into ρ parts.The dynamic model of every unit was constructed based on ultrasonic function angle(γ)and dynamic chip area.Then,the stability region of RRUC is obtained based on the semi-discrete method(SDM).Compared with the robotic conventional countersinking (RCC), RRUC improves the stability by 27%.Finally, the correctness and effectiveness of the stability region model are proved by robotic ultrasonic countersinking experiments.

    1.Introduction

    In modern aerospace industry,riveting plays an important role to connect skin and structure parts of the aircraft.1–3Materials used in modern aerospace industry include composite materials, titanium alloys and aluminum alloys, which are difficult to be machined.4–6According to statistics,the number of rivets on C919 passenger aircraft is more than one million.Drilling and countersinking are the basis of riveting.7–10For such a huge number of drilling and countersinking,it is urgent to find a higher-precision more efficient and economic machining method.In recent years, industrial robots have been widely used in manufacturing field due to their high flexibility, great machining accuracy and economic performance.11–14However, the weak stiffness of industrial robots has become an important factor restricting the further development.15–17During countersinking,the axial force is positively correlated with the cutting width and depth.18The time-varying makes the chatter of industrial robots easier to occur.Serious chatter will not only affect countersinking depth, but also cause irreversible damage to industrial robots and tools.19–22Therefore,the chatter suppression of industrial robots has become a hot topic.23–26

    Nowadays, several methods are proposed to suppress the chatter of industrial robots.The stability lobe diagram is used to predict robotic stability.27–29The zero order approximation(ZOA), the semi-discretization method (SDM) and the fulldiscretization method(FDM)are usually used to obtain stability lobe diagram.30–33Roukema and Altintas34presented the dynamic chip thickness of drilling.The ZOA was used to establish the drilling stability model.Cen et al.35proposed that the angle between the main stiffness direction of the robot and the milling force vector direction can be controlled by using the conservative congruence transformation (CCT) stiffness model.It has been proved to avoid the vibration mode coupling chatter of the robotic milling.Meanwhile,Cordes et al.36established the dynamic cutting layer thickness and dynamic milling force model of robot milling.The stability lobe diagram of robot milling is obtained by the ZOA and SDM.Xin et al.37analyzed the tool-workpiece separation model,time-varying process damping model and axial cutting depth modulation model of robotic milling.The time delay coefficient was introduced to explain the improvement of stability region when there was a low-frequency vibration.Shi et al.38analyzed the machining stability of the 5-DOF parallel robot for helical milling.It is found that the spindle speed and the machining pose of the parallel robot play an important role in the stability.Ozturk et al.39put forward a model to predict the cutting thickness of helical milling.According to system dynamic model, the time delay between the spindle speed and the tool path is obtained.And the helical milling stability is solved in discrete-time and frequency domains.It has been proved that strengthening the structural rigidity of the robot and installing various dampers can reduce the vibration during robot machining effectively.40Guo et al.41presented that robotic boring chatter was a forced vibration with displacement feedback.Drigalski et al.42showed an end effector with vibration damping angle.It was used to suppress chatter vibration in robotic drilling.The vibration damping angle can improve the stability of industrial robot efficiently by increasing the frictional force.Chen and Zhao43designed a new eddy current damper which is connected to the spindle.The eddy current damper suppressed the chatter of industrial robots by changing the dynamic characteristics of tool tips.The results indicated that the peak value of the frequency response function at the tool tip decreases by 22.1%and 12.4%in the Z direction, respectively.Yuan et al.44carried out a semi-active MREs shock absorber with adjustable stiffness.The device was fixed on the robot spindle to absorb the vibration energy of mode coupling flutter during machining.In recent years,ultrasonic machining is used to remove difficult-to-machine materials.45–49It has been proved to expand the stability lobe diagram efficiently as a new processing method.50,51Sun et al.52,53established the stability model of rotary ultrasonic machining (RUM).It avoided the robotic milling chatter successfully.The results showed that the stability lobe diagram of milling has increased by 133%when ultrasonic amplitude was 10 μm.At the same time, longitudinal-torsional ultrasonic milling was also proposed.With the intake of torsional ultrasonic, there was a separation between tool and workpiece,which expanded the stability lobe diagram by 124.42%.Dong et al.54proposed the robotic ultrasonic drilling (RUD) firstly.And the stability model of RUD was established.The drilling stability region was extended successfully.In this paper, ultrasonic vibration is used in industrial robots to suppress the chatter.

    However, the above research shows that the existing study on robotic chatter stability is mainly focused on drilling and milling.There are few reports on countersinking stability.With the increase of feeding time,the width and depth of countersinking also increase.During countersinking, the force is positively correlated with both the width and depth.Since the excitation of countersinking is dynamic instead of constant, the current research method of stability analysis is not suitable for the robotic countersinking.Therefore, it is urgent to propose a new method for RRUC stability analysis.In Section 2,the dynamic force model of differential unit in countersinking is presented.In Section 3, the mathematical model of SDM is detailed.In Section 4, the experiments are carried out to verify the correctness of the stability lobe diagrams.In Section 5, the reasons for stability region extension of RRUC and effect of ultrasonic parameters on RRUC stability are discussed.

    2.Dynamic model of robotic rotary ultrasonic countersinking

    2.1.Structural dynamics

    In RRUC, the force on each flute is composed of tangential force Ft, radial force Frand axial force Fa, as shown in Fig.1.Since ultrasonic vibration is along the axial direction and tangential direction of flute, a three-degree-of-freedom(3-DOF) equivalent model is established to solve RRUC stability region diagram.The dynamic countersinking system can be formulated as follows:

    Fig.1 Elemental forces analysis of a two-fluted countersinking bit.

    where F is the dynamic force and the vector U represents the deflection of the countersinking in the Cartesian coordinates.Meanwhile, F and U can be expressed as

    The matrices M,C and K contain the lumped mass,damping and stiffness characteristics of the RRUC system, respectively.The parameter bdrepresents the diameter of drill.

    2.2.Motion analysis of RRUC

    Compared with robotic conventional countersinking (RCC),RRUC has different trajectories which lead to a tendency of separation between tool and workpiece.Moreover, the tool kinematic characteristics are changed.The motion of countersinking bit consists of axial feed,circumferential rotation,axial ultrasonic vibration and torsional ultrasonic vibration.The RRUC system and the motion trajectory of cutting edge are shown in Fig.2.The kinematics equations can be written as

    Fig.2 RRUC system and motion trajectory.

    where R is the horizontal distance of countersinking bit from the central axis(mm);the parameter ω and t represent angular velocity (rad/s) and time (s), respectively; FUdenotes the frequency of ultrasonic vibration (Hz); AU-Land AU-Trepresent amplitudes of longitudinal vibration(μm)and torsional vibration (μm), respectively; vfis the feed rate (mm/s).

    With the intake of ultrasonic energy, longitudinal and torsional velocity are introduced on the cutter.According to Eq.(3), the longitudinal vibration velocity (vU-L) and torsional vibration velocity (vU-T) are written as

    As shown in Fig.3, orthogonal coordinate system OXYZ represents the tool coordinate with origin(point O)at the center of the maximum cutting radius.Z-axis is defined as the feed direction.To analyze the effect of RRUC chatter suppression better,the orthogonal coordinate system OUXUYUZUis established.OUis a point on the cutter where the cutting radius is bξ.XU-axis represents the rotating direction.And the connection of vector OUO is set as YU-axis.ZU-axis is the same as Zaxis in the orthogonal coordinate system OXYZ.

    In RCC, cutting speed is the compound movement of rotary speed (vr) and feed speed (vf).As the feeding speed is far less than rotary speed, the cutting direction is following the XU-axis.However,the intake of ultrasonic energy changes the cutting direction by introducing axial vibration velocity(vU-L) and torsional vibration velocity (vU-T).Axial vibration velocity(vU-L)and torsional vibration velocity(vU-T)are along the ZU-axial direction and rotating direction,respectively.The cutting direction (vC) changes in the shaded area periodically.To describe periodic change of cutting direction better, the ultrasonic function angle (γ) is defined.γ represents the angle between negative direction of Z-axis and the cutting direction(vC).η is the ratio of torsional vibration amplitude to longitudinal vibration amplitude.The value of γ can be expressed as follows:

    Fig.3 Intake of ultrasonic function angle.

    As shown in Eq.(5),when vU-L>|vf|,γ is an obtuse angle;when vU-L<| vf|, γ is a sharp angle.

    2.3.Dynamic chip area in countersinking

    Different from drilling, the chip area which is composed of chip width and thickness changes with time.The countersinking is divided into ρ parts to simplify the model.For the ξ(1 < ξ < ρ) part, as shown in Fig.4, the chip width bξ(mm)is represented by feeding speed and countersinking time.Meanwhile, it can be turned into the relationship of countersinking depth.According to kinematic relationship, the parameter bξcan be expressed as follows:

    where dt is the processing time under every infinitesimal(s);the parameter H and f represent the final depth of countersinking(mm)and the axial feeding speed(mm/s),respectively;κtis the half tip angle of the countersinking bit.

    As for the chip thickness of ξ part, it can be divided into static chip and dynamic thickness.For the static thickness, it is the result of feed per revolution and the axial ultrasonic vibration.The feed per revolution and the axial ultrasonic vibration in this work can be expressed as

    The dynamic thickness is caused by the regenerative chatter in three directions.The regenerative displacements of the orthogonal coordinate system OXYZ are du,which is the combination of dx,dy and dz.Based on the work of Altintas,34the dynamic thickness on each lip is:

    Fig.4 Analysis of chip area for ξ part.

    2.4.Dynamic countersinking force

    The linear function is assumed to describe the force on each flute.There is a linear relationship between the tangential force and actual width of chip and dynamic chip thickness.Meanwhile,the radial and axial force are proportional to the tangential force.The countersinking force on each lip can be expressed as

    As illustrated in Fig.3, the directions of tangential and axial force are fluctuated with the period change of ultrasonic vibration.The forces projection on three axis of the orthogonal coordinate system OUXUYUZUcan be expressed as

    The cutting force in the orthogonal coordinate system OXYZ is given by the geometric transformation:

    Since the static chip thickness hshas no contribution to stability, it could be neglected.The dynamic force in orthogonal coordinate system OXYZ is given by the dynamic thickness.Substituting Eq.(9) and Eq.(10) into Eq.(11), we rewrite the result in matrix form as

    All the elements of α(t) and β(t) read as

    3.Stability analysis by SDM

    At present, SDM was proposed as a time-domain method.It plays an important role to obtain stability lobe diagrams efficiently.According to the work of Insperger,31the ndimensional delay differential equation can be written as

    where τ is the time delay and T is the time period.In SDM,the length of Δt represents the construction of the time interval division [ti, ti+1], which means T = kΔt.And k is an integer that can be considered as an approximation parameter regarding the time period.

    In this paper,the third-order SDM is obtained based on the second-order system.In the time interval division [ti, ti+1], it can be expressed as

    where x(t-τ)i,y(t-τ)iand z(t-τ)iare the items of time delay;J,E and S are written by three matrixes:

    where ξnx, ξnyand ξnzare the relative damping; ωnx, ωnyand ωnzare the angular natural frequency; mtx, mtyand mtzare the modal mass.

    Substitute Eqs.(18)–(20) into Eq.(17), and it can be expressed as

    where I represents unit matrix.The parameter m can be written as

    As for ui,it is the function of ti,which can be represented as ui= u(ti).And u(ti) can be written as

    Similarly,ui+1is determined by the initial condition,and it can be written as

    where Pi= exp(AiΔt) and

    A 3 m + 6 dimensional state vector can be expressed as

    From above analysis, zi+1can be calculated by the transform vector Diand zi.Diis expressed as

    By coupling Eq.(27)for i=0,1,...,k-1,the 3 m+6 dimensional transition matrix Φ can be determined as

    For a given spindle speed and depth of countersinking,Φ is a judgement index.If all Φ eigenvalues in modulus are less than 1,the system is stable.Otherwise,the system is at the loss of stability.As mentioned in Section 2,the width and depth of countersinking change in every differential unit.It is necessary to guarantee the stability in every unit.Then the stability boundary curve is determined by the judgement of every unit.For the machining centers with good stiffness,only the fundamental frequency and corresponding damping ratio are considered in stability analysis.However, the industrial robot has a weaker stiffness compared to machining centers.In order to make stability lobe diagrams more accurate, multi-order modes of robotic countersinking system should be considered.In this paper,a main purpose is to figure out the effect of ultrasonic energy on robotic countersinking stability.It is assumed that the modes are independent.And the final lobe diagram is obtained by taking the lowest envelop of the stability lobe diagram.Based on the analysis model in Section 2 and Section 3,the chatter stability lobe diagrams of each mode are presented in Fig.5 and Fig.6.Fig.5 represents the predicted stability curves of RRUC (AU-L= 5 μm, AU-T= 3 μm).And Fig.6 is achieved without ultrasonic vibration.It can be found that the curve of mode 1 is not shown in both Fig.5 and Fig.6.The mode 1 of robotic system has no impact on stability lobes at high spindle speeds(more than 2000 r/min).The final stability regions are shown in Fig.7 based on Fig.5 and Fig.6.

    Fig.5 RRUC stability region(AU-L=5 μm,AU-T=3 μm)with different modes (modes 2–5).

    Fig.6 RCC stability region (AU-L = 0 μm, AU-T = 0 μm) with different modes (modes 2–5).

    Fig.7 Comparison of stability region with or without ultrasonic amplitude.

    As shown in Fig.7,the blue points are the verified machining parameters.Circles represent stable countersinking and crosses mean unstable countersinking.The feed rate per revolution and ultrasonic vibration frequency are 0.02 mm/r and 40 kHz, respectively.It demonstrates that the stability region has a certain improvement with ultrasonic vibration.Compared with the high spindle speed,the improvement of stability boundary curve is greater under the low spindle speed.There is a maximum increase in stability when the spindle speed is 2300 r/min.Under RRUC and RCC, the critical countersinking depth is 1.70 mm and 1.33 mm, respectively.The RRUC stability is improved by 27 %.

    4.Experimental verification

    In order to verify the accuracy of stability lobe diagram prediction, countersinking experiments are carried out on robotic machining system shown in Fig.8.The RRUC system consists of robotic countersinking system and miniaturized rotary ultrasonic system.Robotic countersinking system includes a KUKA KR500 robot and an end effector.The miniaturized ultrasonic machining system is mounted on end effector.It consists of a non-contact power transmission, a transducer and a horn.During countersinking,the transducer can convert the low-frequency voltage to high-frequency mechanical vibration.Through the horn, the vibration form is changed and amplitude is amplified.Finally, the mechanical vibration is transmitted to the drilling-countersinking tool.In this paper,longitudinal-torsional ultrasonic vibration is used and the ultrasonic vibration frequency is 40 kHz.The longitudinal vibration amplitude is 5 μm and the torsional vibration amplitude is 3 μm.

    The influence of the ultrasonic vibration on the robotic countersinking stability is investigated by countersinking carbon fiber reinforced plastic (T700) plate.The plate with size of 280 mm-260 mm-6 mm is fixed via pressing plate on both sides.The drill diameter of drilling-countersinking tool is 4.1 mm, the countersink angle is 100 degrees, and the number of lips is 2.The three-axis accelerometer PCB-356A15 is connected with the spindle in end effector system.The acceleration signal in machining is time-varying recorded.The sampling frequency is 10240 Hz.The spindle of end effector is Renaud 430–120380-HSK32E.The parameter is presented in Table 1.According to the method mentioned in Ref.,51the results of countersinking force coefficients kt, krand kaare shown in Table 2.The modal parameters in robotic countersinking posture are obtained from the frequency response functions(FRFs) measured by standard modal hammering test, as shown in Table 3.As mentioned in Section 3, modes 1–5 are selected to predict the stability lobes.

    As shown in Fig.7, twenty points are selected to verify the accuracy of stability region.The stability of each machining point is marked in Fig.7, where circles and crosses are used to indicate whether the system is stable or not, respectively.The red points are in RRUC and the black ones are not in ultrasonic state.In countersinking, chatter is more likely to occur in feeding direction.Therefore, the acceleration signal in feeding direction is selected to determine whether the system is stable.The point 17 in Fig.7 is taken as an example.In RCC and RRUC, the time domain signals and the frequency spectrum after fast Fourier transform(FFT)are drawn respectively in Fig.9.Firstly, from the perspective of time-domain signal,the acceleration value of RRUC is significantly smaller than that of RCC, which means that its fluctuation is smaller and tends to be stable in countersinking.Secondly, from the spectrum analysis,the spindle speed of point 17 is 4000 r/min.The frequency of each order in RRUC spectrum is in rotational frequency and frequency doubling, and there is no large deviation.In RCC, there is a large deviation between 889 Hz and 957 Hz.And the relevant doubling of rotational frequency are 867 Hz and 933 Hz.It indicates that the energy gathers at the frequency doubling of the offset rotational frequency.Meantime, the peak value is much larger than RRUC in the above two order frequencies, which further proves the existence of energy concentration.In summary, RRUC is stable at point 17 while RCC is unstable.

    Fig.8 Experimental setup for stability identification.

    Table 1 Spindle parameter of end effector.

    Table 2 Values of countersinking coefficients.

    Table 3 Robotic modes in X/Y/Z directions (KUKA KR500).

    Through the above method, the stable state of the twenty points is shown in Fig.7.For RRUC, the prediction of the machining state is consistent with the stable lobe diagram.However, the prediction of the state at point 2 is wrong for RCC.It is mainly because the machining parameter is located at the edge of the stable lobe diagram.The accuracy of the prediction cannot be fully guaranteed at this time.In addition to point 2,the stability results of its machining parameters are the same as the theoretical prediction.

    5.Discussion

    5.1.Reasons for stability region extension of RRUC

    Compared with RCC, the stability region of RRUC has a significant improvement, especially at low spindle speed.There are two main reasons.Firstly, the intake of ultrasonic energy changes the direction of external excitation in robotic countersinking.As shown in Fig.3, each force is no longer coaxial with the axis of the original coordinate system OXYZ.There is an ultrasonic function angle (γ) between each force and the original coordinate system.As the force is vector, tangential force, radial force and axial force are mapped to the original coordinate system.And then the vector is added and subtracted.From Fig.3, it can be found that the mapping of tangential force and axial force offsets each other, reducing the size of the external excitation effectively.Thus it makes the stability region of RRUC have a certain improvement.

    Secondly,since the form of ultrasonic vibration in this theory is the combination of longitudinal vibration and torsional vibration, torsional vibration and spindle rotation affect each other in a vibration period.In general,the two have the opposite velocity direction at the same time in half of an ultrasonic period, and the other half has the same velocity direction.When the velocity of torsional vibration is opposite to the spindle rotation velocity, its value is greater than that of the spindle rotation.There will be a separation between the tool and workpiece.This phenomenon explains the enhancing effect of ultrasonic vibration at low spindle speed compared to high spindle speed.Under different countersinking depth,the critical separation spindle speed is also different, and the calculation formula can be given as

    When the torsional vibration speed is the same as the spindle speed,the cutting speed is faster,which is equivalent to the increase of the spindle speed,so that the material is easier to be removed.It reduces the external excitation generated in material removal and improves the stability of RRUC.

    5.2.Effect of ultrasonic parameters on RRUC stability

    According to the above analysis, ultrasonic vibration mainly improves the stability region of RRUC by influencing the ultrasonic function angle.This section mainly explores the influence of ultrasonic parameters on the ultrasonic function angle.According to Eq.(5), the parameters which affect the ultrasonic functional angle mainly include ultrasonic frequency FU,longitudinal vibration amplitude AU-Land torsional vibration amplitude AU-T.The variation range is from γminto γmax,which are obtained when the ultrasonic vibration takes the extreme value, respectively, and can be expressed as

    To figure out the influence of ultrasonic parameters on stability region,the difference between γmaxand γminis defined as Δγ, which can be expressed as

    Fig.11 Effect of different η on Δγ (spindle speed is 4000 r/min,FU = 40 kHz).

    As vfis much lower than vrand ultrasonic velocity, it is ignored in computation to simplify computation.Δγ represents the range of ultrasonic functional angle.The larger the Δγ,the larger the shadow area in Fig.3.There is a positive correlation between Δγ and stable region.

    Firstly,the influence of ultrasonic frequency and longitudinal vibration amplitude on stability region is explored.As shown in Fig.10, when the spindle speed is 4000 r/min and η is 0.6, different ultrasonic frequency is taken to explore the variation of ultrasonic function angle Δγ under different longitudinal vibration amplitude.With the increase of longitudinal vibration amplitude and ultrasonic frequency, Δγ also has a great improvement.However, when the ultrasonic frequency is large enough, Δγ will increase slowly with the increase of longitudinal ultrasonic amplitude.When the ultrasonic frequency is 10 kHz,there is almost a linear relationship between Δγ and the longitudinal vibration amplitude.And when the ultrasonic frequency is 40 kHz, there is a slow-growth parabolic form, which is no longer a linear relationship.It can be found that both ultrasonic frequency and longitudinal vibration amplitude have a great effect on Δγ.In the case of the same spindle speed and torsional vibration amplitude to longitudinal vibration amplitude, the influence of ultrasonic frequency is greater than the longitudinal vibration ultrasonic amplitude.

    When the amplitude of longitudinal vibration is constant,the amplitude of torsional vibration is determined by η.Thus,the effect of torsional vibration amplitude on Δγ is converted to the effect of η.The comparison of Δγ under different η is shown in Fig.11.With the increase of torsional vibration amplitude, Δγ also has an improvement.When η is 0.6, Δγ improves by 10.4 % compared to the case when η is 0.1.

    5.3.Effect of modal parameters on stability region

    In RRUC, the stability region can be affected by all of modes selected in Table 3.However, different ultrasonic parameters have different effects on each mode.In order to explore the sensitivity of each mode to ultrasonic energy,the stability lobe diagrams of each mode under different ultrasonic parameters are drawn respectively.

    Fig.12 Stability region of each mode (modes 2–5) under different ultrasonic energy (η = 0.6).

    As shown in Fig.12, the increase of ultrasonic frequency and amplitude expands the stability lobe diagram in different modes.The area between stability domain curve and horizontal axis is defined to evaluate the improvement of stability lobe diagram.With the increase of ultrasonic energy, the extension radio of stability lobe diagram in different modes are 15.4%,13.4%,17.3%and 18%.According to the results,the stability lobe diagram of mode 5 has the greatest improvement.And the improvement rate of mode 4 is slightly smaller than that of mode 5.As for mode 2 and mode 3, they are relatively insensitive to the increase of ultrasonic parameters.

    6.Conclusions

    In this paper, a stability analysis method is proposed for robotic rotary ultrasonic countersinking (RRUC) systems.Based on the robotic ultrasonic countersinking carbon fiber reinforced plastic (CFRP-T700) plates, the feasibility of this method is confirmed.According to the prediction model and validation results, three conclusions can be drawn as follows:

    (1) This study establishes the RRUC stability model on the basis of the countersinking dynamic analysis.Meanwhile,the experimental results agree with RRUC stability boundary curve well.Compared with robotic conventional countersinking(RCC),the RRUC stability is improved by 27%.

    (2) The intake of ultrasonic energy improves the RRUC stability region by changing the tool kinematics.Among them, the reduction of dynamic countersinking force is the direct reason.RRUC stability region enlarges with the growth of the ultrasonic parameters such as the ultrasonic vibration frequency and amplitude.

    (3) The stability lobe diagram of RRUC is affected by multi-order modes of robotic countersinking system.With the increase of ultrasonic energy, the extension radio of stability region in various modes are different.The improvement effect in higher order mode is larger than that in lower order mode.

    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 Project on the Technological Leading Talent Teams Led by Frontiers Science Center for Complex Equipment System Dynamics (No.FSCCESD220401)and the National Natural Science Foundation of China (No.52075265).

    免费av毛片视频| 成年人黄色毛片网站| 女人被狂操c到高潮| 丁香六月欧美| 国产欧美日韩一区二区三区在线| 19禁男女啪啪无遮挡网站| 又大又爽又粗| 大码成人一级视频| 色婷婷久久久亚洲欧美| 国产又色又爽无遮挡免费看| 日韩av在线大香蕉| 免费久久久久久久精品成人欧美视频| 欧美日韩av久久| 久久国产乱子伦精品免费另类| 五月开心婷婷网| 在线观看免费日韩欧美大片| 国产人伦9x9x在线观看| 欧美日韩精品网址| 在线永久观看黄色视频| 在线观看免费日韩欧美大片| 夜夜躁狠狠躁天天躁| 久久国产精品影院| 国产免费av片在线观看野外av| 久久精品国产亚洲av高清一级| 日本vs欧美在线观看视频| 麻豆一二三区av精品| 999精品在线视频| 久久天堂一区二区三区四区| 亚洲精华国产精华精| 亚洲欧美精品综合久久99| 精品久久久精品久久久| 国产成人系列免费观看| 久久久久久久精品吃奶| www.熟女人妻精品国产| 亚洲专区国产一区二区| 黄片小视频在线播放| 一级,二级,三级黄色视频| 欧美丝袜亚洲另类 | 国产免费男女视频| 村上凉子中文字幕在线| 欧美在线一区亚洲| 99久久综合精品五月天人人| 亚洲va日本ⅴa欧美va伊人久久| 新久久久久国产一级毛片| 色综合站精品国产| 亚洲免费av在线视频| 亚洲av成人一区二区三| 新久久久久国产一级毛片| 日韩视频一区二区在线观看| 久久精品国产99精品国产亚洲性色 | 国产国语露脸激情在线看| 黄片播放在线免费| 女人被躁到高潮嗷嗷叫费观| 国产三级黄色录像| 亚洲免费av在线视频| 国产伦一二天堂av在线观看| 国产xxxxx性猛交| 免费女性裸体啪啪无遮挡网站| 操美女的视频在线观看| 女性被躁到高潮视频| 一级a爱片免费观看的视频| 国产伦人伦偷精品视频| 夫妻午夜视频| 又紧又爽又黄一区二区| 极品教师在线免费播放| aaaaa片日本免费| 亚洲欧美一区二区三区久久| 久久人人精品亚洲av| 国产精品av久久久久免费| 免费女性裸体啪啪无遮挡网站| 国产极品粉嫩免费观看在线| 制服人妻中文乱码| 亚洲va日本ⅴa欧美va伊人久久| 美女大奶头视频| 亚洲精品中文字幕一二三四区| 少妇的丰满在线观看| 亚洲自偷自拍图片 自拍| 巨乳人妻的诱惑在线观看| 国产成人精品久久二区二区91| 色综合欧美亚洲国产小说| 女性被躁到高潮视频| 久久亚洲精品不卡| 久久午夜亚洲精品久久| 欧美最黄视频在线播放免费 | 1024视频免费在线观看| 国产日韩一区二区三区精品不卡| 精品电影一区二区在线| 中文字幕另类日韩欧美亚洲嫩草| 黑人猛操日本美女一级片| 久热爱精品视频在线9| 免费av毛片视频| 男人的好看免费观看在线视频 | 欧美黑人欧美精品刺激| 80岁老熟妇乱子伦牲交| 巨乳人妻的诱惑在线观看| 国产精品一区二区三区四区久久 | 亚洲精品av麻豆狂野| 国产精品一区二区精品视频观看| www.自偷自拍.com| 在线永久观看黄色视频| 一个人免费在线观看的高清视频| 视频在线观看一区二区三区| 女人精品久久久久毛片| 亚洲在线自拍视频| 亚洲午夜精品一区,二区,三区| av天堂在线播放| 长腿黑丝高跟| 高清黄色对白视频在线免费看| 精品午夜福利视频在线观看一区| 久久香蕉精品热| 国产精品一区二区三区四区久久 | 国产精品自产拍在线观看55亚洲| www.精华液| 精品少妇一区二区三区视频日本电影| 五月开心婷婷网| 国产区一区二久久| 在线观看免费视频日本深夜| 国产成年人精品一区二区 | 在线观看一区二区三区激情| 国产精品偷伦视频观看了| xxxhd国产人妻xxx| 国产av在哪里看| 免费女性裸体啪啪无遮挡网站| 午夜成年电影在线免费观看| 天天躁夜夜躁狠狠躁躁| a级毛片在线看网站| 成人18禁高潮啪啪吃奶动态图| 丝袜在线中文字幕| 国产亚洲精品久久久久5区| 免费久久久久久久精品成人欧美视频| 女人被躁到高潮嗷嗷叫费观| 69av精品久久久久久| 黄色丝袜av网址大全| 欧美日韩黄片免| 精品久久久久久久毛片微露脸| 欧美中文日本在线观看视频| 亚洲精品一卡2卡三卡4卡5卡| 国产一卡二卡三卡精品| 国产区一区二久久| 激情视频va一区二区三区| 女性生殖器流出的白浆| 多毛熟女@视频| 999精品在线视频| 黑丝袜美女国产一区| 成人三级做爰电影| 欧美日韩亚洲高清精品| cao死你这个sao货| 窝窝影院91人妻| 亚洲精品一区av在线观看| а√天堂www在线а√下载| 黑人猛操日本美女一级片| 国产成人欧美| 满18在线观看网站| 脱女人内裤的视频| 久久亚洲精品不卡| 中文字幕人妻熟女乱码| 亚洲一区中文字幕在线| 亚洲中文av在线| 少妇的丰满在线观看| 精品少妇一区二区三区视频日本电影| 丝袜美腿诱惑在线| 亚洲在线自拍视频| 亚洲自拍偷在线| 一级,二级,三级黄色视频| 可以在线观看毛片的网站| 国产男靠女视频免费网站| cao死你这个sao货| 精品国内亚洲2022精品成人| avwww免费| 香蕉国产在线看| 久久精品国产综合久久久| 在线观看一区二区三区| 女警被强在线播放| 欧美精品亚洲一区二区| 亚洲精品美女久久久久99蜜臀| 日本vs欧美在线观看视频| 国产精品98久久久久久宅男小说| 国产三级黄色录像| 欧美大码av| 伊人久久大香线蕉亚洲五| 在线看a的网站| 黑人欧美特级aaaaaa片| 一级a爱片免费观看的视频| 大香蕉久久成人网| 久久精品亚洲熟妇少妇任你| 亚洲,欧美精品.| 日本wwww免费看| 国产精品一区二区精品视频观看| 天天添夜夜摸| 精品一区二区三区av网在线观看| 精品一区二区三卡| 一区在线观看完整版| 欧美人与性动交α欧美软件| 99在线人妻在线中文字幕| 可以免费在线观看a视频的电影网站| 亚洲精品成人av观看孕妇| 琪琪午夜伦伦电影理论片6080| 91麻豆av在线| 99国产极品粉嫩在线观看| 老鸭窝网址在线观看| 欧美日韩亚洲国产一区二区在线观看| 12—13女人毛片做爰片一| 精品国产乱码久久久久久男人| 国产片内射在线| 男人舔女人的私密视频| 最近最新免费中文字幕在线| 亚洲一区二区三区不卡视频| 日本黄色日本黄色录像| 国内久久婷婷六月综合欲色啪| 岛国在线观看网站| 女同久久另类99精品国产91| 免费久久久久久久精品成人欧美视频| 午夜久久久在线观看| 精品第一国产精品| 精品一区二区三区视频在线观看免费 | 精品国产国语对白av| 午夜精品国产一区二区电影| 啦啦啦 在线观看视频| 国产熟女午夜一区二区三区| 制服人妻中文乱码| 亚洲中文字幕日韩| 1024香蕉在线观看| 桃红色精品国产亚洲av| 国产精品一区二区免费欧美| 国产亚洲精品久久久久久毛片| 1024香蕉在线观看| 国产欧美日韩一区二区精品| 高潮久久久久久久久久久不卡| 国产一区二区三区视频了| 女生性感内裤真人,穿戴方法视频| 人人妻人人澡人人看| 在线观看舔阴道视频| av欧美777| 电影成人av| 国产精品永久免费网站| 久久精品国产99精品国产亚洲性色 | 欧美激情久久久久久爽电影 | 青草久久国产| 91在线观看av| 中文字幕人妻丝袜一区二区| 久久伊人香网站| 悠悠久久av| 人成视频在线观看免费观看| 在线国产一区二区在线| 欧美激情 高清一区二区三区| 最新在线观看一区二区三区| 999久久久国产精品视频| 色精品久久人妻99蜜桃| 在线观看舔阴道视频| 人人妻人人爽人人添夜夜欢视频| 国产av一区在线观看免费| 国产三级在线视频| 国产麻豆69| 亚洲三区欧美一区| 一级片免费观看大全| 国产熟女xx| 亚洲激情在线av| 性欧美人与动物交配| 欧美亚洲日本最大视频资源| 日韩高清综合在线| 欧美日韩av久久| 嫩草影院精品99| 夜夜躁狠狠躁天天躁| 久久亚洲精品不卡| 美女扒开内裤让男人捅视频| 亚洲第一av免费看| 国产成人欧美在线观看| 一进一出抽搐gif免费好疼 | 日韩免费av在线播放| 这个男人来自地球电影免费观看| 曰老女人黄片| 精品国产乱码久久久久久男人| 亚洲 欧美一区二区三区| 一级a爱片免费观看的视频| 99久久国产精品久久久| 每晚都被弄得嗷嗷叫到高潮| 真人一进一出gif抽搐免费| 9热在线视频观看99| 精品久久久久久久毛片微露脸| 我的亚洲天堂| 欧美午夜高清在线| svipshipincom国产片| 国产成人系列免费观看| 国产精品亚洲av一区麻豆| 亚洲全国av大片| 午夜免费激情av| 欧美日韩瑟瑟在线播放| 狠狠狠狠99中文字幕| 精品高清国产在线一区| 美女高潮到喷水免费观看| 国产精品一区二区在线不卡| 亚洲精品久久成人aⅴ小说| 在线av久久热| 亚洲少妇的诱惑av| 美女 人体艺术 gogo| 久久久久久免费高清国产稀缺| 免费人成视频x8x8入口观看| 岛国在线观看网站| 免费观看人在逋| 看黄色毛片网站| 一二三四社区在线视频社区8| 高清在线国产一区| 美女午夜性视频免费| 久久这里只有精品19| 国产一区二区三区视频了| 在线观看午夜福利视频| 激情在线观看视频在线高清| 亚洲国产精品一区二区三区在线| 757午夜福利合集在线观看| 俄罗斯特黄特色一大片| 午夜91福利影院| 黄色视频不卡| 免费在线观看日本一区| 在线天堂中文资源库| 级片在线观看| 宅男免费午夜| 国产精品一区二区免费欧美| 欧美黑人欧美精品刺激| 精品国产乱子伦一区二区三区| 日韩欧美一区二区三区在线观看| 丰满人妻熟妇乱又伦精品不卡| 免费少妇av软件| 两人在一起打扑克的视频| 黑人巨大精品欧美一区二区蜜桃| 国产麻豆69| 午夜精品国产一区二区电影| 视频区图区小说| 亚洲一区二区三区色噜噜 | 日韩精品中文字幕看吧| 99国产精品一区二区蜜桃av| 97人妻天天添夜夜摸| 亚洲七黄色美女视频| 欧美人与性动交α欧美软件| а√天堂www在线а√下载| 国产伦一二天堂av在线观看| 免费女性裸体啪啪无遮挡网站| 9191精品国产免费久久| 久久中文字幕一级| 久久精品亚洲熟妇少妇任你| 久久久久久久精品吃奶| 国产精品99久久99久久久不卡| 亚洲av第一区精品v没综合| 亚洲一区二区三区不卡视频| 男女之事视频高清在线观看| 99国产精品一区二区三区| 99国产综合亚洲精品| 国产成人精品久久二区二区免费| 五月开心婷婷网| 日日爽夜夜爽网站| 淫秽高清视频在线观看| 国产亚洲精品综合一区在线观看 | 国产精品一区二区精品视频观看| 中国美女看黄片| 97人妻天天添夜夜摸| 啪啪无遮挡十八禁网站| 侵犯人妻中文字幕一二三四区| 夜夜爽天天搞| 水蜜桃什么品种好| 18禁裸乳无遮挡免费网站照片 | 国产乱人伦免费视频| 亚洲av电影在线进入| 色综合婷婷激情| 午夜福利免费观看在线| 亚洲一区二区三区不卡视频| 狠狠狠狠99中文字幕| www日本在线高清视频| 中文字幕人妻丝袜制服| 老司机深夜福利视频在线观看| 欧美日韩乱码在线| a级毛片黄视频| 亚洲黑人精品在线| 亚洲av五月六月丁香网| 日韩高清综合在线| 国产成人欧美| 精品国产一区二区久久| 国产三级在线视频| 成人三级做爰电影| 丰满的人妻完整版| 国产在线精品亚洲第一网站| 黄色视频不卡| 久久精品国产综合久久久| 国产av精品麻豆| 亚洲专区字幕在线| 久久中文看片网| 国产区一区二久久| 桃红色精品国产亚洲av| 美女福利国产在线| 久热这里只有精品99| 99久久国产精品久久久| 久久久久国产精品人妻aⅴ院| 国产精品九九99| 婷婷丁香在线五月| 老司机亚洲免费影院| 美女高潮喷水抽搐中文字幕| aaaaa片日本免费| 国产精品免费视频内射| 免费高清在线观看日韩| 色尼玛亚洲综合影院| 欧美+亚洲+日韩+国产| 麻豆av在线久日| 美国免费a级毛片| 亚洲精品美女久久av网站| 国产精品永久免费网站| 色在线成人网| 99久久国产精品久久久| 午夜视频精品福利| 最近最新免费中文字幕在线| 热re99久久国产66热| 免费在线观看完整版高清| 欧美日韩亚洲国产一区二区在线观看| 乱人伦中国视频| 欧美成狂野欧美在线观看| 97人妻天天添夜夜摸| 极品教师在线免费播放| 中文字幕av电影在线播放| 欧美老熟妇乱子伦牲交| 亚洲第一欧美日韩一区二区三区| 超碰97精品在线观看| 麻豆av在线久日| 久久精品亚洲熟妇少妇任你| 精品久久久久久久毛片微露脸| 日日摸夜夜添夜夜添小说| 老汉色av国产亚洲站长工具| 丰满迷人的少妇在线观看| 亚洲精品一二三| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品一区av在线观看| 欧美成人性av电影在线观看| 午夜亚洲福利在线播放| 黑人猛操日本美女一级片| 悠悠久久av| 男女之事视频高清在线观看| 亚洲色图综合在线观看| 国产免费av片在线观看野外av| 亚洲第一欧美日韩一区二区三区| 国产精品影院久久| 国产精品成人在线| av免费在线观看网站| 中文字幕人妻丝袜制服| 99国产精品99久久久久| 长腿黑丝高跟| 国产人伦9x9x在线观看| 巨乳人妻的诱惑在线观看| 精品一区二区三区四区五区乱码| 成人av一区二区三区在线看| 久久久久久大精品| 91大片在线观看| 亚洲精品美女久久av网站| 亚洲av片天天在线观看| 在线看a的网站| 国产亚洲av高清不卡| 久久久国产成人免费| 伊人久久大香线蕉亚洲五| 桃色一区二区三区在线观看| 80岁老熟妇乱子伦牲交| av国产精品久久久久影院| 极品教师在线免费播放| 亚洲精品一卡2卡三卡4卡5卡| 中文欧美无线码| 黑人巨大精品欧美一区二区蜜桃| 久久久久国产一级毛片高清牌| 久久精品aⅴ一区二区三区四区| 国产精品一区二区在线不卡| 久久久久久大精品| 亚洲av五月六月丁香网| 一区在线观看完整版| 亚洲少妇的诱惑av| 国产精品永久免费网站| 女生性感内裤真人,穿戴方法视频| 一本综合久久免费| 国产精品久久视频播放| 18禁观看日本| 一级片'在线观看视频| 欧美激情 高清一区二区三区| 在线观看免费高清a一片| 18禁美女被吸乳视频| 国产视频一区二区在线看| 麻豆国产av国片精品| 国产亚洲精品第一综合不卡| 久久久久精品国产欧美久久久| 十八禁网站免费在线| 91在线观看av| 中国美女看黄片| 三上悠亚av全集在线观看| 国产成人欧美在线观看| 在线免费观看的www视频| 在线观看免费高清a一片| 满18在线观看网站| 欧美亚洲日本最大视频资源| 性欧美人与动物交配| 久久天堂一区二区三区四区| 色综合婷婷激情| 免费搜索国产男女视频| 国产精品综合久久久久久久免费 | 99久久精品国产亚洲精品| 国产不卡一卡二| 夜夜躁狠狠躁天天躁| 久久中文字幕一级| 亚洲成人免费电影在线观看| 色婷婷av一区二区三区视频| 免费av毛片视频| 看黄色毛片网站| 母亲3免费完整高清在线观看| 少妇粗大呻吟视频| av网站在线播放免费| 日本精品一区二区三区蜜桃| 搡老熟女国产l中国老女人| 老鸭窝网址在线观看| 纯流量卡能插随身wifi吗| 少妇的丰满在线观看| 日韩欧美三级三区| 91九色精品人成在线观看| 久久精品成人免费网站| www.自偷自拍.com| 高清在线国产一区| 亚洲av电影在线进入| 一二三四在线观看免费中文在| 国产欧美日韩综合在线一区二区| 在线观看免费视频网站a站| 精品国产乱子伦一区二区三区| 少妇裸体淫交视频免费看高清 | 国产精品久久久久久人妻精品电影| 国产成人欧美| 一级片'在线观看视频| 99在线人妻在线中文字幕| 国产精品乱码一区二三区的特点 | 免费观看精品视频网站| 欧美成人性av电影在线观看| 悠悠久久av| av网站在线播放免费| 日韩三级视频一区二区三区| 曰老女人黄片| bbb黄色大片| 两个人免费观看高清视频| 看片在线看免费视频| 老鸭窝网址在线观看| 伊人久久大香线蕉亚洲五| 黑人猛操日本美女一级片| 久久精品人人爽人人爽视色| 一级黄色大片毛片| 久久这里只有精品19| 另类亚洲欧美激情| 91九色精品人成在线观看| 男人的好看免费观看在线视频 | 麻豆国产av国片精品| 欧美av亚洲av综合av国产av| 日韩大尺度精品在线看网址 | 欧美另类亚洲清纯唯美| av有码第一页| 久热爱精品视频在线9| aaaaa片日本免费| a级片在线免费高清观看视频| 久久精品国产综合久久久| 高清毛片免费观看视频网站 | 黄色片一级片一级黄色片| 亚洲精品久久成人aⅴ小说| 一区二区三区国产精品乱码| 免费看十八禁软件| 村上凉子中文字幕在线| 精品一区二区三区av网在线观看| 狠狠狠狠99中文字幕| ponron亚洲| 国内久久婷婷六月综合欲色啪| 精品久久久久久久毛片微露脸| 亚洲熟女毛片儿| www国产在线视频色| 日韩三级视频一区二区三区| 777久久人妻少妇嫩草av网站| 18禁黄网站禁片午夜丰满| 国产精品国产av在线观看| 国产亚洲精品一区二区www| 久久精品亚洲av国产电影网| 巨乳人妻的诱惑在线观看| 日日摸夜夜添夜夜添小说| 中亚洲国语对白在线视频| 乱人伦中国视频| xxx96com| 在线看a的网站| 欧美+亚洲+日韩+国产| 国产精品一区二区在线不卡| 男女之事视频高清在线观看| 纯流量卡能插随身wifi吗| 高清在线国产一区| 91av网站免费观看| 麻豆久久精品国产亚洲av | 国产一区二区在线av高清观看| 国产精品一区二区精品视频观看| 亚洲成a人片在线一区二区| 人人妻人人爽人人添夜夜欢视频| 日韩精品中文字幕看吧| 精品一区二区三区四区五区乱码| 亚洲av美国av| 欧美不卡视频在线免费观看 | 一进一出好大好爽视频| 高清av免费在线| 亚洲第一av免费看| 国产人伦9x9x在线观看| 国产三级黄色录像| 国产亚洲精品综合一区在线观看 | 美女 人体艺术 gogo| 久久国产乱子伦精品免费另类| 久久久国产一区二区| 国产高清激情床上av| 搡老乐熟女国产| 天堂√8在线中文| 国产亚洲精品一区二区www| 啪啪无遮挡十八禁网站| 不卡一级毛片| 亚洲一区高清亚洲精品| 亚洲精品久久成人aⅴ小说| 黄色a级毛片大全视频| 黄色视频,在线免费观看| 一级黄色大片毛片| 三级毛片av免费| 巨乳人妻的诱惑在线观看|