• <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一区二区三区在线看| 男女床上黄色一级片免费看| 欧美黑人巨大hd| 亚洲欧美日韩高清专用| 中文在线观看免费www的网站| 亚洲欧美日韩东京热| 老熟妇仑乱视频hdxx| 热99在线观看视频| 我要搜黄色片| 一本一本综合久久| 亚洲不卡免费看| 中文字幕av在线有码专区| 国产视频一区二区在线看| 村上凉子中文字幕在线| 18禁国产床啪视频网站| 岛国在线免费视频观看| 日韩亚洲欧美综合| 亚洲国产精品久久男人天堂| 欧美乱码精品一区二区三区| 国产精品久久久久久久电影 | 亚洲狠狠婷婷综合久久图片| 俺也久久电影网| 国产v大片淫在线免费观看| 国产亚洲精品综合一区在线观看| 在线观看免费视频日本深夜| 亚洲精品成人久久久久久| 亚洲熟妇熟女久久| 午夜亚洲福利在线播放| 在线免费观看不下载黄p国产 | 97超级碰碰碰精品色视频在线观看| 国产又黄又爽又无遮挡在线| 少妇熟女aⅴ在线视频| 亚洲 欧美 日韩 在线 免费| eeuss影院久久| 午夜福利免费观看在线| av女优亚洲男人天堂| 最近最新中文字幕大全电影3| 亚洲av五月六月丁香网| 成人三级黄色视频| 色综合站精品国产| 欧美黑人巨大hd| 99久国产av精品| 久久精品亚洲精品国产色婷小说| 神马国产精品三级电影在线观看| 99热只有精品国产| 免费看十八禁软件| 亚洲国产精品久久男人天堂| 波多野结衣巨乳人妻| 国产免费一级a男人的天堂| 日韩欧美三级三区| 中国美女看黄片| 亚洲avbb在线观看| 天天一区二区日本电影三级| 国产三级在线视频| 久久久精品欧美日韩精品| 黄色片一级片一级黄色片| 高潮久久久久久久久久久不卡| 国产精品亚洲一级av第二区| 国产主播在线观看一区二区| 欧美激情在线99| 免费看日本二区| 成人永久免费在线观看视频| 精品人妻1区二区| 精品久久久久久久毛片微露脸| 操出白浆在线播放| 99精品在免费线老司机午夜| av在线蜜桃| 美女黄网站色视频| 在线观看午夜福利视频| 亚洲精品乱码久久久v下载方式 | 午夜福利在线观看吧| 欧洲精品卡2卡3卡4卡5卡区| www日本在线高清视频| 国产三级在线视频| 免费看美女性在线毛片视频| 日本撒尿小便嘘嘘汇集6| 午夜免费男女啪啪视频观看 | 国产精品亚洲av一区麻豆| 午夜免费男女啪啪视频观看 | 亚洲avbb在线观看| 手机成人av网站| 色噜噜av男人的天堂激情| 两个人视频免费观看高清| 精品午夜福利视频在线观看一区| 搞女人的毛片| 国产成人系列免费观看| 成人午夜高清在线视频| 人妻久久中文字幕网| 国产熟女xx| 女人高潮潮喷娇喘18禁视频| av在线天堂中文字幕| 色视频www国产| 欧美日韩国产亚洲二区| 亚洲真实伦在线观看| 黄色视频,在线免费观看| 成年免费大片在线观看| 国产伦一二天堂av在线观看| 国产午夜精品久久久久久一区二区三区 | 精品国产美女av久久久久小说| 色精品久久人妻99蜜桃| www.熟女人妻精品国产| 欧美日韩综合久久久久久 | 在线a可以看的网站| 国产一区二区在线观看日韩 | 欧美高清成人免费视频www| 又粗又爽又猛毛片免费看| 欧美性猛交黑人性爽| 99热6这里只有精品| 亚洲午夜理论影院| 亚洲美女黄片视频| 偷拍熟女少妇极品色| 国产精品亚洲av一区麻豆| 成人高潮视频无遮挡免费网站| 亚洲国产精品sss在线观看| 天天添夜夜摸| 国产亚洲欧美98| 搞女人的毛片| 最新中文字幕久久久久| 国产精品 国内视频| 高清在线国产一区| 亚洲中文字幕日韩| 天堂√8在线中文| 在线观看日韩欧美| 日韩欧美在线乱码| 女人高潮潮喷娇喘18禁视频| 精品久久久久久久毛片微露脸| 国产69精品久久久久777片| 久久精品国产亚洲av涩爱 | 色综合亚洲欧美另类图片| 午夜福利在线观看吧| 中文字幕人妻熟人妻熟丝袜美 | 日本撒尿小便嘘嘘汇集6| 内射极品少妇av片p| 男女下面进入的视频免费午夜| 一边摸一边抽搐一进一小说| 久久国产精品人妻蜜桃| 日本黄大片高清| 久久久久久国产a免费观看| 一级黄片播放器| 无限看片的www在线观看| 一级a爱片免费观看的视频| 啪啪无遮挡十八禁网站| 天堂影院成人在线观看| 日本黄色视频三级网站网址| 国产一区二区在线av高清观看| 很黄的视频免费| 搞女人的毛片| 草草在线视频免费看| 一本久久中文字幕| 国产av麻豆久久久久久久| 亚洲精品一区av在线观看| 亚洲国产精品成人综合色| 一本精品99久久精品77| 老熟妇仑乱视频hdxx| 日本免费a在线| 成人18禁在线播放| 亚洲成人久久性| 18禁在线播放成人免费| 99热这里只有精品一区| 一级黄色大片毛片| 欧美黑人欧美精品刺激| 在线看三级毛片| 国语自产精品视频在线第100页| 亚洲国产欧洲综合997久久,| 色视频www国产| 色老头精品视频在线观看| 午夜免费观看网址| 最近在线观看免费完整版| 亚洲中文日韩欧美视频| 91久久精品电影网| 国产一区二区在线观看日韩 | 婷婷亚洲欧美| 色吧在线观看| 一进一出抽搐动态| 在线十欧美十亚洲十日本专区| 高清日韩中文字幕在线| 国产成人aa在线观看| 欧美乱色亚洲激情| 欧美午夜高清在线| 亚洲国产色片| 无限看片的www在线观看| 99久久精品热视频| 日韩亚洲欧美综合| 18+在线观看网站| 亚洲美女黄片视频| 欧美成人一区二区免费高清观看| 日日夜夜操网爽| 日韩成人在线观看一区二区三区| 母亲3免费完整高清在线观看| 午夜精品一区二区三区免费看| 成人国产一区最新在线观看| 亚洲欧美一区二区三区黑人| 国产av一区在线观看免费| 国内少妇人妻偷人精品xxx网站| 性色av乱码一区二区三区2| 国产在视频线在精品| 亚洲精品亚洲一区二区| 精品免费久久久久久久清纯| 久久人妻av系列| 亚洲真实伦在线观看| 18禁裸乳无遮挡免费网站照片| 午夜免费男女啪啪视频观看 | 久久伊人香网站| 久久香蕉精品热| av在线天堂中文字幕| 丰满人妻熟妇乱又伦精品不卡| 伊人久久大香线蕉亚洲五| 国产精品 国内视频| 波多野结衣巨乳人妻| 欧美在线黄色| 国产欧美日韩一区二区精品| 国产亚洲精品一区二区www| 色综合亚洲欧美另类图片| 亚洲人成伊人成综合网2020| 黄色成人免费大全| 久久久久久国产a免费观看| 中出人妻视频一区二区| 成人特级黄色片久久久久久久| av在线天堂中文字幕| 国产av一区在线观看免费| 嫩草影视91久久| 久久亚洲真实| 国内精品久久久久精免费| 五月伊人婷婷丁香| 欧美日韩一级在线毛片| 婷婷亚洲欧美| 露出奶头的视频| 中文在线观看免费www的网站| 久久人人精品亚洲av| 一本久久中文字幕| 久久精品国产亚洲av涩爱 | 成人性生交大片免费视频hd| 日韩国内少妇激情av| 欧美黄色片欧美黄色片| 可以在线观看的亚洲视频| 婷婷亚洲欧美| 久久久久久大精品| 国产综合懂色| 97人妻精品一区二区三区麻豆| 真人一进一出gif抽搐免费| 亚洲欧美日韩东京热| 国产一区二区激情短视频| 亚洲不卡免费看| 日本五十路高清| 久久性视频一级片| 亚洲国产精品合色在线| 床上黄色一级片| 成人亚洲精品av一区二区| 又粗又爽又猛毛片免费看| 性色avwww在线观看| 成人永久免费在线观看视频| 哪里可以看免费的av片| 亚洲久久久久久中文字幕| 97碰自拍视频| 欧美乱妇无乱码| 又黄又爽又免费观看的视频| 久久久久久久午夜电影| 色尼玛亚洲综合影院| av中文乱码字幕在线| 欧美成狂野欧美在线观看| 在线十欧美十亚洲十日本专区| 国产精品日韩av在线免费观看| 亚洲第一电影网av| 日韩欧美精品v在线| 国产欧美日韩一区二区精品| 三级男女做爰猛烈吃奶摸视频| 国内毛片毛片毛片毛片毛片| 高潮久久久久久久久久久不卡| 久久这里只有精品中国| xxxwww97欧美| av在线蜜桃| 一区二区三区免费毛片| 国产精品乱码一区二三区的特点| 国产一区二区激情短视频| 欧美大码av| 国产一区在线观看成人免费| 热99re8久久精品国产| 在线a可以看的网站| 狠狠狠狠99中文字幕| 中文字幕av在线有码专区| 岛国视频午夜一区免费看| 精品久久久久久成人av| 亚洲自拍偷在线| 18禁黄网站禁片午夜丰满| 亚洲国产高清在线一区二区三| 女人十人毛片免费观看3o分钟| 亚洲成av人片在线播放无| 五月玫瑰六月丁香| 国内久久婷婷六月综合欲色啪| av天堂在线播放| 麻豆成人午夜福利视频| 一个人看的www免费观看视频| 国产精品女同一区二区软件 | 亚洲在线自拍视频| 国产免费av片在线观看野外av| 久久久国产成人精品二区| 午夜福利高清视频| 高清毛片免费观看视频网站| 国产精品一区二区三区四区免费观看 | 又粗又爽又猛毛片免费看| 97超视频在线观看视频| 日韩欧美在线二视频| 97人妻精品一区二区三区麻豆| 99热6这里只有精品| 每晚都被弄得嗷嗷叫到高潮| 俺也久久电影网| 国语自产精品视频在线第100页| 一a级毛片在线观看| 给我免费播放毛片高清在线观看| 老司机在亚洲福利影院| 精品免费久久久久久久清纯| 搡老熟女国产l中国老女人| 国产一区二区在线av高清观看| 国产精品一区二区免费欧美| 色综合亚洲欧美另类图片| 校园春色视频在线观看| 伊人久久大香线蕉亚洲五| 日韩高清综合在线| 波野结衣二区三区在线 | av黄色大香蕉| 国产伦精品一区二区三区视频9 | 黄色成人免费大全| 亚洲国产精品成人综合色| 日本熟妇午夜| 免费观看人在逋| 久久精品国产综合久久久| 国产精品三级大全| 一级毛片女人18水好多| 国产精品免费一区二区三区在线| 亚洲久久久久久中文字幕| 国产一区在线观看成人免费| 中出人妻视频一区二区| 美女大奶头视频| 99久久精品国产亚洲精品| 日本免费一区二区三区高清不卡| 精品一区二区三区视频在线 | e午夜精品久久久久久久| 亚洲av成人不卡在线观看播放网| 十八禁人妻一区二区| 一个人看的www免费观看视频| 女人十人毛片免费观看3o分钟| 日韩国内少妇激情av| 观看美女的网站| 香蕉久久夜色| 超碰av人人做人人爽久久 | www.999成人在线观看| 午夜久久久久精精品| 18禁黄网站禁片免费观看直播| 国产97色在线日韩免费| 99国产极品粉嫩在线观看| 中文字幕精品亚洲无线码一区| 久久久久久九九精品二区国产| 精品欧美国产一区二区三| 亚洲中文字幕一区二区三区有码在线看| 亚洲五月婷婷丁香| 亚洲国产高清在线一区二区三| 午夜免费激情av| 国产视频一区二区在线看| 在线观看午夜福利视频| 欧美日韩综合久久久久久 | 日本精品一区二区三区蜜桃| 亚洲av成人精品一区久久| 99热只有精品国产| 亚洲av成人精品一区久久| 日韩欧美三级三区| 亚洲在线自拍视频| 亚洲av成人av| 日本一本二区三区精品| 国产精品免费一区二区三区在线| 真实男女啪啪啪动态图| 美女黄网站色视频| 亚洲在线自拍视频| 亚洲熟妇熟女久久| a级一级毛片免费在线观看| 国产成人福利小说| 啦啦啦观看免费观看视频高清| 五月玫瑰六月丁香| 脱女人内裤的视频| 麻豆国产97在线/欧美| 国产激情欧美一区二区| 久久久久久国产a免费观看| 小蜜桃在线观看免费完整版高清| 90打野战视频偷拍视频| 91九色精品人成在线观看| 亚洲性夜色夜夜综合| 亚洲午夜理论影院| 观看美女的网站| 搞女人的毛片| 欧美在线一区亚洲| 一级黄色大片毛片| 久久九九热精品免费| 免费在线观看影片大全网站| 成人性生交大片免费视频hd| 国产三级中文精品| 色视频www国产| 成人三级黄色视频| 亚洲久久久久久中文字幕| 国产免费av片在线观看野外av| 中文字幕人妻熟人妻熟丝袜美 | 欧美乱码精品一区二区三区| 国产黄色小视频在线观看| 国产毛片a区久久久久| 夜夜爽天天搞| 精品不卡国产一区二区三区| 免费观看的影片在线观看| 精品一区二区三区人妻视频| 日韩有码中文字幕| 国产精品影院久久| 首页视频小说图片口味搜索| 精华霜和精华液先用哪个| 亚洲不卡免费看| 日韩欧美一区二区三区在线观看| 精品人妻一区二区三区麻豆 | 国产成人啪精品午夜网站| 国产欧美日韩一区二区三| 国产又黄又爽又无遮挡在线| 亚洲熟妇中文字幕五十中出| 美女被艹到高潮喷水动态| 午夜福利免费观看在线| 熟妇人妻久久中文字幕3abv| 欧美在线一区亚洲| 全区人妻精品视频| 一级作爱视频免费观看| 国产三级中文精品| 欧美色欧美亚洲另类二区| 国产av在哪里看| 蜜桃亚洲精品一区二区三区| 国产一区二区在线观看日韩 | 国产野战对白在线观看| 国产v大片淫在线免费观看| 国产高清激情床上av| 欧美大码av| 亚洲内射少妇av| 国产探花在线观看一区二区| 国产亚洲欧美在线一区二区| 人人妻人人看人人澡| 香蕉丝袜av| www日本在线高清视频| 欧美成人一区二区免费高清观看| 我的老师免费观看完整版| 亚洲人与动物交配视频| av天堂在线播放| 特大巨黑吊av在线直播| 国产精品久久久久久精品电影| 男女那种视频在线观看| 日本精品一区二区三区蜜桃| 真人做人爱边吃奶动态| 亚洲av第一区精品v没综合| av天堂中文字幕网| 欧美色视频一区免费| 中出人妻视频一区二区| 97人妻精品一区二区三区麻豆| 欧美极品一区二区三区四区| 丰满人妻一区二区三区视频av | 国产免费av片在线观看野外av| 国产成人av激情在线播放| 亚洲国产欧美人成| 在线十欧美十亚洲十日本专区| 最后的刺客免费高清国语| 国产精品综合久久久久久久免费| 久久精品国产99精品国产亚洲性色| 日韩 欧美 亚洲 中文字幕| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 变态另类丝袜制服| 久久精品国产亚洲av涩爱 | 最近最新中文字幕大全电影3| 九色成人免费人妻av| 国产伦在线观看视频一区| 免费搜索国产男女视频| 老汉色∧v一级毛片| 国产日本99.免费观看| 最新中文字幕久久久久| 欧美一级a爱片免费观看看| 最近视频中文字幕2019在线8| 少妇高潮的动态图| 久9热在线精品视频| 欧美日本亚洲视频在线播放| 精品久久久久久久久久久久久| 变态另类丝袜制服| 天天一区二区日本电影三级| 亚洲五月婷婷丁香| 免费看a级黄色片| 亚洲国产精品久久男人天堂| 麻豆国产97在线/欧美| 极品教师在线免费播放| 久久国产精品影院| 女警被强在线播放| avwww免费| 一级毛片女人18水好多| 一夜夜www| 国产老妇女一区| 51国产日韩欧美| 久久精品国产自在天天线| 最近最新中文字幕大全电影3| 欧美日韩中文字幕国产精品一区二区三区| 成人高潮视频无遮挡免费网站| 久久性视频一级片| 色在线成人网| 男女那种视频在线观看| 欧美一区二区亚洲| 女同久久另类99精品国产91| bbb黄色大片| 国产不卡一卡二| 欧美一级a爱片免费观看看| 免费大片18禁| 九九在线视频观看精品| 在线免费观看不下载黄p国产 | 国产男靠女视频免费网站| x7x7x7水蜜桃| www国产在线视频色| 亚洲精品在线美女| 日日干狠狠操夜夜爽| 国产熟女xx| 90打野战视频偷拍视频| 亚洲国产欧美网| 午夜福利免费观看在线| 亚洲国产欧美网| 国产一区二区三区在线臀色熟女| 午夜日韩欧美国产| 国产精品自产拍在线观看55亚洲| 欧美高清成人免费视频www| 国产精品精品国产色婷婷| 91在线观看av| 叶爱在线成人免费视频播放| 好看av亚洲va欧美ⅴa在| 观看美女的网站| 国产精品一区二区三区四区免费观看 | 日本 av在线| 亚洲第一电影网av| 国产私拍福利视频在线观看| 99热这里只有是精品50| 亚洲熟妇熟女久久| 波多野结衣巨乳人妻| 国产精品精品国产色婷婷| 两性午夜刺激爽爽歪歪视频在线观看| 国产高清有码在线观看视频| 国产精品99久久99久久久不卡| 国模一区二区三区四区视频| 999久久久精品免费观看国产| 99精品久久久久人妻精品| 又粗又爽又猛毛片免费看| av片东京热男人的天堂| www日本黄色视频网| 日本三级黄在线观看| 人人妻人人看人人澡| 亚洲精品粉嫩美女一区| 国产又黄又爽又无遮挡在线| 91麻豆精品激情在线观看国产| 中文字幕人妻熟人妻熟丝袜美 | 夜夜夜夜夜久久久久| 在线视频色国产色| 一卡2卡三卡四卡精品乱码亚洲| 欧美激情久久久久久爽电影| 我的老师免费观看完整版| 性色avwww在线观看| 亚洲最大成人中文| 黄色视频,在线免费观看| 搞女人的毛片| 欧美zozozo另类| 欧美最新免费一区二区三区 | 在线免费观看的www视频| 免费观看精品视频网站| 欧美性猛交黑人性爽| 天天一区二区日本电影三级| 国产精品久久视频播放| 国产一级毛片七仙女欲春2| 久99久视频精品免费| 国产野战对白在线观看| 精品久久久久久成人av| 久久久久久久久大av| 啦啦啦观看免费观看视频高清| 午夜福利免费观看在线| 内射极品少妇av片p| 中国美女看黄片| 日韩欧美 国产精品| 成年免费大片在线观看| 国产不卡一卡二| xxx96com| 90打野战视频偷拍视频| 无限看片的www在线观看| 国产欧美日韩一区二区精品| 国产高清激情床上av| 中文资源天堂在线| 欧美色视频一区免费| 草草在线视频免费看| 国产熟女xx| 99在线视频只有这里精品首页| 国产高清激情床上av| 91九色精品人成在线观看| 精品欧美国产一区二区三| 久久亚洲精品不卡| 在线观看舔阴道视频| 亚洲欧美日韩高清专用| 国产伦精品一区二区三区视频9 | 国产亚洲欧美在线一区二区| 九九在线视频观看精品| 婷婷丁香在线五月| 国产精品爽爽va在线观看网站| 精品国产亚洲在线| 亚洲男人的天堂狠狠| 精品久久久久久久毛片微露脸| 国产99白浆流出| 国产精品 欧美亚洲| 欧美日韩乱码在线| 日韩人妻高清精品专区| 久久久久精品国产欧美久久久| 搡老妇女老女人老熟妇| 午夜福利高清视频| 观看免费一级毛片| 欧美在线一区亚洲| 日本与韩国留学比较| 免费人成在线观看视频色|