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

    Investigation of anode shaping process during co-rotating electrochemical machining of convex structure on inner surface

    2023-09-05 09:45:06ShuofangZHOUDengyongWANGWenjianCAODiZHU
    CHINESE JOURNAL OF AERONAUTICS 2023年7期

    Shuofang ZHOU, Dengyong WANG, Wenjian CAO, Di ZHU,*

    a College of Mechanical and Electrical Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China

    b Jiangsu Key Laboratory of Precision and Micro-Manufacturing Technology, Nanjing 210016, China

    KEYWORDS Convex structure forming process;Co-rotating motion;Current density;Electrochemical machining;Revolving part

    Abstract A novel co-rotating electrochemical machining method is proposed for fabricating convex structures on the inner surface of a revolving part.The electrodes motion and material removal method of co-rotating electrochemical machining are different from traditional electrochemical machining.An equivalent kinematic model is established to analyze the novel electrodes motion,since the anode and cathode rotate in the same direction while the cathode simultaneously feeds along the line of centres.According to the kinematic equations of the electrodes and Faraday’s law, a material removal model is established to simulate the evolution of the anode profile in co-rotating electrochemical machining.The simulation results indicate that the machining accuracy of the convex structure is strongly affected by the angular velocity ratio and the radius of the cathode tool.An increase of the angular velocity ratio can improve the machining accuracy of a convex structure.A small difference in the radius of the cathode tool will cause changes in the shape of the sidewalls of the convex structure.The width of the cathode window affects only the width of the convex structure and the inclination α of the sidewall.A relation between the width of the cathode window and the width of the convex structure was obtained.The formation process for a convex structure under electrochemical dissolution was revealed.Based on the simulation results,the optimal angular velocity ratio and cathode radius were selected for an experimental verification,and 12 convex structures were simultaneously fabricated on the inner surface of a thin-walled revolving part.The experimental results are in good agreement with the simulation results, which verifies the correctness of the theoretical analysis.Therefore, inner surface co-rotating electrochemical machining is an effective method for fabricating convex structures on the inner surface of a revolving part.

    1.Introduction

    Aviation and aerospace applications use many thin-walled revolving parts made of difficult-to-cut materials, such as engine casings and jet engine rings, whose surfaces have complex convex structures.1,2The cutting force can easily cause the deformation of revolving parts machined by conventional mechanical methods, and repairing the deformation is extremely difficult.3,4Electrochemical machining (ECM) is a noncontact form of machining, and it uses electrochemical anodic dissolution to remove material from the workpiece.5,6Because of the advantages such as regardless of material hardness, no relevant tool wear, no machining deformation and good surface quality, ECM has been widely used for the machining of difficult-to-cut materials.7–9.

    There has been progress in the electrochemical fabrication of surface structures on a revolving part.Liu et al.10used a rotating dentate cathode to machine a micro-inner annular groove in a small diameter metallic tube, and the machining mechanisms were analyzed.Zhang et al.11optimized the cathode structure and successfully fabricated small inner-walled ring grooves in an injector nozzle with ECM.Tang et al.12conducted research on ECM of a special-shaped inner spiral tube,and various flow models were compared and analyzed.Patel et al.13fabricated the microstructure on the curved surface by through-mask electrochemical micromachining(TMECMM), and the shaping mechanisms were analyzed numerically.Mahdavinejad and Hatami14applied ECM to polish the inner surface of a gun barrel chambers, achieving good performance.Sheng et al.15used blocky cathode tools to machine thin-walled revolving parts.However, multiple block cathode tools need to be prepared according to the shape of the convex structures on the surface of the revolving part.Zhu et al.16proposed counter-rotating electrochemical machining (CRECM) for machining the outer surface of a thin-walled revolving part.CRECM can simultaneously fabricate all the convex structures on the anode surface, which significantly improves the machining accuracy and efficiency.Wang et al.17fabricated intensive cylindrical pillar array by counter-rotating electrochemical machining, and the cathode tool was manufactured using additive manufacturing technology.

    For the internal structure of a thin-walled revolving part,conventional machining tools cannot effectively reach the machining area and the machining performance is limited.Because of the large amount of material that needs to be removed from the inside of the workpiece, conventional machining requires long processing cycles and has a high machining cost.3In CRECM, the anode workpiece and cathode tool have similar radii and rotate synchronously in opposite directions, which enables the method to simultaneously fabricate all the convex structures on the outer surface of a revolving part.However, CRECM cannot machine the internal structure of a revolving part.At present, there is no effective method for fabricating convex structures on the inner surface of large and thin-walled revolving parts,such as engine casings.For these kinds of special parts, inner surface co-rotating electrochemical machining (ICRECM) is proposed.ICRECM can simultaneously fabricate all the convex structures on the inner surface of a revolving part, which can significantly improve the machining efficiency.The electrodes motion and material removal method of co-rotating electrochemical machining are different from traditional electrochemical machining and counter-rotating electrochemical machining (CRECM).

    The ECM process is quite complex,in which the machining accuracy is affected by many factors.18,19In order to obtain a desired machining accuracy,the evolution of anode profile has been studied in ECM.20Jain and Gehlot21established a mathematical model to simulate the evolution of the anode profile in through-mask electrochemical micromachining (TMECMM).Zhitnikov et al.22applied the theory of functions of a complex variable for investigation of the electrochemical machining process.The stepwise function of current efficiency was used to determine the movement velocity of the anode boundary.In order to predict the anode profile in wire electrochemical micromachining (Wire-ECMM), Sharma et al.23established a two-dimensional simulation model using finite element method(FEM).Predictions from the developed model agreed well with the experimental results with a minimum error of 3%.Chen et al.24presented a multi-physics model and the quasi steady state shortcut (QSSSC) approach to simulate the changes of workpiece shape in pulse electrochemical machining (PECM).To improve the surface quality and the machining accuracy, low frequency vibrations was added to the cathode tool.Hewidy et al.25presented an analytical approach to assess the mechanism of metal removal for ECM assisted by low frequency vibrations.

    This paper analyzes the formation of convex structures on the inner surface of a revolving part.A mathematical model based on the kinematic equations of the electrodes and anodic material dissolution is established.The trajectories of contour points on the cathode window are plotted according to the mathematical equations.The influence of the angular velocity ratio, cathode tool radius and cathode window width on the machining accuracy of a convex structure are analyzed in detail.It is found that an increase of the angular velocity ratio can improve the machining accuracy of a convex structure.A small difference in the radius of the cathode tool will cause changes in the shape of the sidewalls of the convex structure.The width of the cathode window affects only the width of the convex structure and the inclination α of the sidewall.The formation process for a convex structure under electrochemical dissolution was revealed.The excess dissolution of material results in differences between the simulated profile under electrolysis and the trajectory envelope.Based on the simulation results,the optimal angular velocity ratio and cathode radius were selected for an experimental verification, and 12 convex structures were simultaneously fabricated on the inner surface of a thin-walled revolving part, indicating that ICRECM has good machining performance.The experimental results are in good agreement with the simulation results,which verify the theoretical analysis.

    Fig.1 Principle of ICRECM.

    2.Mathematical model

    2.1.Inner surface co-rotating electrochemical machining

    Fig.1 illustrates ICRECM.The cathode tool is a cylindrical electrode with one or more windows.The interior of each window is insulated.The anode workpiece rotates at a constant angular velocity (ω).The cathode tool simultaneously rotates inside the anode workpiece in the same direction but at n times the angular velocity of the workpiece.The cathode tool is fed toward the anode workpiece at a constant velocity (v).The high-speed electrolyte flushes through the narrow electrode gap from the side.The anode material is dissolved rapidly when a certain voltage is applied between the anode and the cathode.Thus, convex structures can be fabricated on the inner surface of the anode workpiece.ICRECM is quite different from conventional ECM because of the unique and complex movement of the electrodes.

    2.2.Kinematic model

    To accurately simulate the motion of ICRECM,an equivalent kinematic model is established.In ICRECM, the anode and cathode rotate in the same direction while the tool is simultaneously fed linearly toward the workpiece, which makes the overall motion extremely complex.We consider a coordinate system where the anode workpiece is assumed to be stationary.Fig.2 shows the equivalent two-dimensional decomposition of the machining process in this workpiece coordinate system.The center of the anode is the origin O1(0,0), and the center of the cathode is O2(X0, Y0).The motion of the cathode tool can be decomposed into three components:26,27.

    (a) The cathode tool rotates counterclockwise around the origin O1.

    (b) The cathode tool rotates clockwise around its own center O2at n times the angular velocity of the workpiece.

    (c) The cathode tool is linearly fed toward the anode workpiece.

    Fig.2 Equivalent kinematic model for ICRECM.

    In the present study,we determine the equations of motion for each step and derive the equivalent kinematic equations.The coordinates of any point M (X, Y) on the tool cathode at any time can then be obtained from the equivalent kinematic equations.

    First,the cathode tool rotates counterclockwise around the center O1of the anode workpiece with angular velocity ω for time t.The coordinates of any point on the cathode surface M1(X1, Y1) are determined as follows:

    2.3.Material removal model

    The electric field distribution on the irregular revolving surface is difficult to directly obtain an analytical solution,so the finite element method is adopted.The domains and boundaries for material removal model by ICRECM are shown in Fig.3.These include the electrolyte domain Ω, anode boundary Γ1,and cathode boundaries Γ2,Γ3,Γ4,and Γ5.The anode boundary Γ1is also the inner surface of the electrolyte tank.In order to improve the continuity of the physical fields in specific regions, local mesh refinement strategies are adopted, including boundary-layer refinement and adaptive refinement.The motion of the cathode boundaries follows Eqs.(7) and (8).In addition, the coordinates of the anode boundary Γ1as it changes during machining can be obtained from the electrochemical dissolution rate of anode surface points.Thus, the overall geometric shape of the convex structure can be predicted.

    In this simulation,the material removal model of ICRECM is simplified with the following assumptions:28–31(1)The cathode tool and the anode workpiece are equipotential surfaces with different potential values.(2)The current density distribution on each electrode surface obeys Ohm’s law.(3) The electrolyte flow within the inter-electrode gap is ideal, and the electrical conductivity k is approximately constant.

    The potential in the electrolyte domain Ω satisfies the Laplace equation:32.

    The boundaries of the insulated windows Γ3, Γ4, Γ5are as follows:

    Fig.3 Boundaries and computational mesh for ICRECM.

    where u is the normal coordinates at the boundary surface.

    The normal current density on the electrode surface is proportional to the local potential gradient:33.

    where E→xand E→yare the components of the electric field intensity along x and y directions respectively.

    During the anode shaping process, the whole processing time is assumed to be T.The whole process can be divided into N cycles:

    where(xj,yj)and(xj+1,yj+1)are the coordinates of any point on the anode boundary.

    According to Eq.(14), the actual volume electrochemical equivalent ηW has an important role in the simulation.In this study, the anode workpiece in ICRECM is made of 304 stainless steel(304SS),and the actual volume electrochemical equivalent can be obtained by weighting loss measurements.34The function curve between the actual volume electrochemical equivalent ηW and the current density (i) for 304SS in 20 wt%NaNO3at 30°C is shown in Fig.4.The function can be fitted based on data points as follows:35.

    Fig.4 Function curve between actual volume electrochemical equivalent and current density (i) for 304SS in 20wt% NaNO3 at 30 °C.35

    Fig.5 Schematic diagram of material removal.

    Fig.5 shows the material removal process in the simulation.The anode workpiece is assumed to be stationary.It is connected to the positive pole of the power supply.The cathode tool has a complex motion.It is connected to the negative pole of the power supply.The anode boundary Γ1moves outward at velocity v→eunder the action of the electric field, so that the inner radius of the anode gradually increases.Over time, the amount of material dissolved in each subsequent time interval increases.

    3.Simulation of the formation of a convex structure

    3.1.Determination of angular velocity ratio

    3.1.1.Trajectories of contour points on the cathode window

    In ICRECM, the formation of a convex structure mainly depends on the trajectories of significant contour points on each cathode window.The coordinates of the contour points can be calculated by Eqs.(7)and(8)at different times.The trajectories of the contour points can then be obtained by connecting these coordinates.In this study, we consider different angular velocity ratios n between the anode workpiece and the cathode tool.The other simulation parameters are listed in Table 1.

    Fig.6 shows the trajectories of two contour points, A and B, for different angular velocity ratios n = 3, 6, or 9.Theradius of the cathode tool has to be correspondingly reduced with an increase of the angular velocity ratio.The number of convex structures fabricated depends on the angular velocity ratio.For n = 3, points A and B touch the inner surface of the anode workpiece three times, so three convex structures are produced, as shown in Fig.6(a).With an increase of the angular velocity ratio,the trajectories of points A and B make contact more frequently with the inner surface of the circular anode,and so the window on the cathode surface will produce more structures.It can be seen from Fig.6(c)that nine convex structures are manufactured when n=9,and the convex structures are evenly distributed.

    Table 1 Simulation parameters.

    To describe and analyze the simulation results, a series of structural parameters are defined, as shown in Fig.7.Based on the characteristics of a convex structure, some accuracy indexes are defined to evaluate the convex structure,including:the fillet RDat the top,the fillet RGat the root,the inclination angle α of the sidewall,the height H,the width D and the angle β.

    The width of the convex structure is defined as the chord length of the circle with radius R3, where R3is calculated as follows:

    Fig.6 Trajectories of contour points A and B for different angular velocity ratios n.

    Fig.7 Simulation structure parameters in ICRECM.

    where R2is the inner radius of the anode workpiece after processing, and R4is the initial inner radius of the anode workpiece.

    The width D and the angle β of the convex structure satisfy the following formula:

    where RCis the cathode tool radius.

    3.1.2.Anode profiles under electrochemical dissolution

    Fig.8 Simulated profiles of convex structures for different angular velocity ratios.

    The machining accuracy of a convex structure strongly depends on the angular velocity ratio.To explore the formation of convex structures for different angular velocity ratios,electrochemical dissolution is simulated for a range of angular velocity ratios (n = 2–10).The width of the cathode window DCis fixed at 18 mm, and the other simulation parameters are listed in Table 1.

    Fig.9 Accuracy indexes of a convex structure for different angular velocity ratios.

    Fig.10 Current density distribution in ICRECM.

    Fig.11 Current density distribution on anode surface for various angular velocity ratios.

    Fig.12 Trajectories of cathode tool for different radii when angular velocity ratio n = 6.

    Fig.8 shows the simulated profiles of convex structures when n=3,6,or 9.There are clearly differences.The accuracy indexes are plotted in Fig.9 for the different angular velocity ratios.Fig.9(a) shows the trends for D, H, RD, and RG.The width D of the convex structure increases with an increase of the angular velocity ratio.The height H decreases linearly with an increase of angular velocity ratio over a certain range.With an increase of the angular velocity ratio, the top fillet RDand root fillet RGfirst decrease and then remain unchanged.As shown in Fig.9(b), the inclination angle α of the sidewall decreases rapidly with an increase of the angular velocity ratio.

    To explore why the angular velocity ratio affects the machining accuracy of a convex structure, the current density distribution in the electrolyte domain Ω is calculated with simulation software(COMSOL 6.0).As shown in Fig.10,there is a high current density region only in the narrow gap near the processing area.In addition, the current density on the anode surface in the processing area is calculated for various angular velocity ratios, as shown in Fig.11(a).The peak current density is the same for the different angular velocity ratios.Fig.11(b)shows the angle of the region where the current density is greater than 40 A/cm2for different angular velocity ratios.With an increase of the angular velocity ratio,the angle of the high current density region decreases sharply at first and then decreases slowly.The smaller the angle of the high current density region, the better the localization and the higher the machining accuracy of the convex structure.However, if the angular velocity ratio is too high, the cathode radius will be small, which may overly limit the width of the cathode window.

    3.2.Determination of cathode tool radius

    Here, the angular velocity ratio n = 6 is used to study the influence of the cathode tool radius on the formation of a convex structure.The other simulation parameters are listed in Table 1.Fig.12 shows close-ups of the trajectories of cathode tools with different radii.When RC= 16.5 mm, both sides of the convex structure have a conical profile.For RC=17 mm,the sidewalls of the convex structure are V-shaped.When RC-≥17.5 mm, the sidewalls are straight.The inclination angle α and the fillet RGof the convex structure both increase as the radius of the cathode tool increases.Regardless of the cathode radius,the trajectories of contour points A and B always have a non-zero inclination angle, which suggests that the inclination angle α of a convex structure cannot be zero.A small difference in the radius of the cathode tool changes the trajectories of contour points A and B, and so the contour of convex structure also changes.

    The profile of the convex structure is simulated under electrochemical dissolution when the angular velocity ratio n = 6 and the cathode radius RC= 17.5 mm.Fig.13(a) shows the evolution of the convex structure on the inner surface of a revolving part at different machining times.The cathode tool is fed continuously toward the anode workpiece,and the material of the anode workpiece is rapidly removed under the high current density.The convex structure becomes narrower over time, and the fillet RDand the fillet RGgradually decrease.The inclination angle α of the sidewall becomes significantly smaller.The corrosion caused by stray current has a great impact on the machining accuracy of the convex structure.The error between the simulated profile under electrochemical dissolution and the trajectory envelope is large at the top of the convex structure, as shown in Fig.13(b).In Section 3.4, the causes of the error are discussed in detail.

    3.3.Determination of cathode window width

    The machining accuracy of the convex structure is also affected by the width D of the cathode window.The angular velocity ratio n = 6 and the cathode radius RC= 17.5 mm are used in the simulation.The other simulation parameters are listed in Table 1.Fig.14 shows the simulated profiles of convex structures with different values for the angle of the cathode window.As the angle θ increases, the angle β of the convex structure also increases.The angle θ of the cathode window and the central angle γ for the trajectories of the contour points satisfy the following:

    To quantitatively analyze the profiles of the convex structures for different values of the angle of the cathode window,the accuracy indexes are plotted in Fig.15.It can be seen from Fig.15(a)that the fitting curves of the simulation results for H,RD, and RGare horizontal, which indicates that they do not change with an increase of the angle of the cathode window.Fig.15(b) shows that as the angle of the cathode window increases in a certain range, the inclination angle α and the angle β increase linearly.The angle θ of the cathode window and the angle β of the convex structure satisfy the following:

    Fig.13 Simulated profiles of a convex structure.

    Fig.14 Simulated profiles of convex structures with different values for angle of cathode window.

    Therefore, the width of the cathode window can be calculated according to the desired width of the convex structure.

    3.4.Formation of a convex structure under electrochemical dissolution

    In ICRECM, the error between the simulated profile under electrochemical dissolution and the trajectory envelope is large at the top of the convex structure, as shown in Fig.13(b).To explore the causes of this error,the current density on the surface of the convex structure is researched during the 500th rotation.With an angular velocity of 6 r/min, it takes 8 s for the cathode tool to move from position Ⅰto position Ⅱ, as shown in Fig.16.In this paper, 304SS is used as the material of the anode workpiece.The minimum current density to break the passive film of 304SS is 3.96 A/cm2.34,35When the current density is lower than 3.96 A/cm2, the material will not dissolve.

    Fig.15 Accuracy indexes of a convex structure for different values of the angle of the cathode window.

    Fig.16 Cathode movement during the 500th rotation.

    Fig.17 Current density distribution on surface of convex structure at different machining times.

    Fig.17 shows the current density distribution on the surface of the convex structure at different times.The motion can be divided into two stages according to the position of the cathode tool relative to the convex structure.In the first stage,from 0 to 4 s, the cathode rotates close to the convex structure, as shown in Fig.17(a)–(d).The stray current sweeps through points F, G, and P in turn.The current density gradually increases during this time.In the second stage, from 4 to 8 s,the cathode rotates away from the convex structure, as shown in Fig.17(d)–(g).The right profile of the convex structure is attacked again by the stray current in this time.Similarly,the left profile of the convex structure is also attacked twice by the stray current.In summary, the convex structure is attacked twice by the stray current during each rotation.

    The current densities of characteristic points at different times are shown in Fig.18.The figure gives the start and end times of material dissolution for each point.The current density curves for symmetrically positioned characteristic points are symmetrical,which means that the convex structure fabricated is also symmetrical.

    There are two peaks in the current density curve of point S in time 0–8 s, as shown in Fig.18(a).However, the peak current density is lower than 3.96 A/cm2, so the material near point S does not dissolve.In Fig.18(b) and (c) there are also two peaks in the current density curves, and the peak current density exceeds 3.96 A/cm2.The material near points F, J,G, and K is dissolved twice in time 0–8 s.The material near points F and J dissolves to form the fillet RDat the top of the convex structure.As shown in Fig.18(d), the current density curves of points P and Q have only one peak,and the peak current densities are as high as 72 A/cm2.The period when the current density exceeds 3.96 A/cm2lasts for 1.73 s.In Fig.17(c) and(e), the high current density regions are at the roots of the convex structure and include points P and Q.A large amount of material in these regions is removed so that the height of the convex structure increases.

    In summary,the sidewall of the convex structure is attacked by stray current many times during each rotation, resulting in the excess dissolution of the material.In the fabrication of a convex structure, the duration of when the stray current attacks the sidewall increases from the top to the root, which leads to the error in the profile shown in Fig.13(b).Therefore,the error between the simulated profile under electrochemical dissolution and the trajectory envelope is inevitable.

    4.Experimental validation

    4.1.Experimental procedure

    As shown in Fig.19, a self-designed experimental set-up was used to study the formation of convex structures on the inner surface of a revolving part.The electrolyte distributor was used to redistribute the electrolyte, so that the pressure and flow in the jet nozzle were stable.The high-speed electrolyte was sprayed from the jet nozzle and flushed through the narrow electrode gap from the side.The anode workpiece and cathode tool were, respectively, mounted on two rotating shafts, which rotated in the same direction.The cathode tool was simultaneously fed along the line of centres.

    Fig.19 Experimental set-up for ICRECM.

    Workpieces of stainless steel 304 with an inner surface radius of 90 mm were used in the experiment.Two cylindrical cathode tools were used to machine the anode workpieces with ICRECM.One had a radius of 17.5 and the other 31 mm.The angular velocity ratios were 6 and 3,respectively.For the cathode tool with a radius of 17.5 mm, as shown in Fig.20(a)and(b), two rectangular windows of widths 12 and 18 mm were fabricated in it.The windows were insulated by ceramic.The other tool had a single insulated window of width 18 mm, as shown in Fig.20(c).The applied voltage was 25 V, and the electrolyte was 20wt% NaNO3, which was kept at a constant temperature of 30 °C during the experiment.The feed depth of the cathode tool was 8 mm at a constant speed of 0.012 mm/min.

    Fig.20 Photographs of cathode tools.

    4.2.Experimental results

    As shown in Fig.21, multiple convex structures were successfully fabricated on the inner surface of thin-walled revolving parts using cathode tools of different radii.Fig.21(a) shows that the cathode with a single insulated window makes three convex structures on the inner surface of the anode workpiece when the angular velocity ratio was 3.In contrast, Fig.21(b)shows that with the other tool with two windows, when the angular velocity ratio was 6,12 convex structures with two different widths were fabricated on the anode workpiece.

    The coordinates of different points on the profiles of the convex structures were measured.The accuracy indexes of the convex structure were obtained by fitting these points with Matrix Laboratory (MATLAB 2020).

    Fig.21 Pictures of workpieces machined using ICRECM.

    Fig.22 Photographs showing the profiles of fabricated convex structures.

    Table 2 Accuracy indexes for the fabricated convex structures.

    Fig.23 Comparison of simulated and experimental convex profiles (top) and errors (bottom).

    Fig.24 Errors between simulated and experimental results for accuracy indexes.

    Fig.22 shows the profiles of the convex structures fabricated with different experimental parameters.The values of the accuracy indexes for the fabricated convex structures are presented in Table 2.The fillet radii RDand RGand the height H of the convex structures I and II are approximately equal,but the inclination angle α is larger for DC=18 mm.In addition, the widths of the convex structures are 7.95 and 13.82 mm, respectively, which are consistent with Eq.(27).Comparing structures Ⅱand III, the height H, the fillet radii RDand RG, and the inclination angle α are larger for n = 3.In conclusion,the trends for the accuracy indexes of the fabricated convex structures are consistent with the simulation results.

    Fig.23 shows that the simulated and experimental convex profiles are in good agreement and that the errors are within the permissible range.The profiles of the convex structures are symmetrical, so the errors between the simulated and experimental results are also approximately symmetric.Fig.24 shows the errors between the simulated and experimental results for accuracy indexes over the three types of convex structure.The errors for the height H and the width D are small, whereas the errors are larger for the fillet radii.However,the errors between the simulated and experimental results are within the allowable error range for the process,so that the machining performance is acceptable for industrial applications.

    5.Conclusions

    In this paper,mathematical equations are derived for the complex relative motion of the electrodes and the material removal model is established based on the kinematic equations for the electrodes and Faraday’s law.Multiple convex structures are fabricated on the inner surface of a revolving part using ICRECM.The main conclusions can be summarized as follows:

    (1) The material removal process is simulated numerically,and the evolution of anode profile is investigated.Over time, the amount of material dissolved in each time interval increases.

    (2) The simulation results indicate that the machining accuracy of the convex structure strongly depends on the angular velocity ratio, the radius of the cathode tool and the width of the cathode window.As the angular velocity ratio n increases, the inclination angle α decreases significantly.When n ≥10, the inclination angle α is less than 9°.When RC≥17.5 mm, the sidewalls of the convex structure are straight.As the angle θ of the cathode window increases, the angle β of the convex structure also increases.

    (3) The sidewall of the convex structure is attacked by stray current many times during each rotation,resulting in the excess dissolution of the material.The excess dissolution of material during electrolysis causes errors between the simulated profile and the trajectory envelope of the cathode tool.

    (4) Multiple convex structures were simultaneously fabricated on the inner surface of a thin-walled revolving part,indicating that ICRECM has good machining performance.The error between simulated and experimental results for accuracy indexes is less than 7%, which verifies that the proposed simulation model is valid for predicting convex shapes fabricated using ICRECM.

    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 work is supported by the National Natural Science Foundation of China (No.52175414), National Natural Science Foundation of China for Creative Research Groups(No.51921003), the China Postdoctoral Science Foundation(No.2019M661833).

    少妇人妻 视频| 国产欧美亚洲国产| svipshipincom国产片| 成人国产av品久久久| 久久天躁狠狠躁夜夜2o2o | 久久久精品94久久精品| 在线观看免费视频网站a站| 国产精品欧美亚洲77777| videos熟女内射| 一本—道久久a久久精品蜜桃钙片| 一个人免费看片子| 极品人妻少妇av视频| 999精品在线视频| 老鸭窝网址在线观看| 日韩 亚洲 欧美在线| 丝袜人妻中文字幕| 超碰97精品在线观看| 婷婷色综合大香蕉| 好男人电影高清在线观看| 午夜老司机福利片| 亚洲精品一区蜜桃| 亚洲国产av新网站| 久久精品国产a三级三级三级| 国产一区二区在线观看av| 亚洲av成人不卡在线观看播放网 | 亚洲精品美女久久av网站| 国产免费视频播放在线视频| 91精品三级在线观看| 久久国产亚洲av麻豆专区| 亚洲中文av在线| 国产成人啪精品午夜网站| 亚洲人成网站在线观看播放| 久久综合国产亚洲精品| 啦啦啦 在线观看视频| 91精品伊人久久大香线蕉| 国产高清不卡午夜福利| 欧美亚洲 丝袜 人妻 在线| 天堂俺去俺来也www色官网| 别揉我奶头~嗯~啊~动态视频 | 一级毛片我不卡| 你懂的网址亚洲精品在线观看| 中文字幕精品免费在线观看视频| 男女下面插进去视频免费观看| 久热爱精品视频在线9| 免费高清在线观看视频在线观看| 热re99久久精品国产66热6| 观看av在线不卡| 涩涩av久久男人的天堂| 午夜久久久在线观看| 欧美老熟妇乱子伦牲交| 精品少妇久久久久久888优播| videosex国产| 2021少妇久久久久久久久久久| 99国产精品一区二区蜜桃av | www.精华液| 中文字幕人妻丝袜一区二区| 亚洲成人手机| 大型av网站在线播放| 亚洲五月色婷婷综合| 如日韩欧美国产精品一区二区三区| 成人国产一区最新在线观看 | 超碰97精品在线观看| 国产成人91sexporn| 国产淫语在线视频| 国产成人啪精品午夜网站| 亚洲精品美女久久av网站| 国产在视频线精品| 日韩大码丰满熟妇| 成人国产一区最新在线观看 | 悠悠久久av| 少妇的丰满在线观看| 亚洲视频免费观看视频| 真人做人爱边吃奶动态| 国产成人精品久久二区二区91| 亚洲精品在线美女| 成人国产av品久久久| 国产日韩欧美视频二区| 无限看片的www在线观看| 伦理电影免费视频| 色婷婷av一区二区三区视频| 日韩av在线免费看完整版不卡| 欧美97在线视频| 精品卡一卡二卡四卡免费| 又大又爽又粗| 新久久久久国产一级毛片| 欧美大码av| 精品免费久久久久久久清纯 | 丰满迷人的少妇在线观看| 夫妻性生交免费视频一级片| 国产精品麻豆人妻色哟哟久久| 亚洲精品美女久久av网站| 一区二区三区精品91| 男女高潮啪啪啪动态图| 91九色精品人成在线观看| 久久精品国产综合久久久| 亚洲精品中文字幕在线视频| 欧美97在线视频| 99国产综合亚洲精品| 国产色视频综合| 男女边吃奶边做爰视频| 91九色精品人成在线观看| 亚洲国产成人一精品久久久| 国产又色又爽无遮挡免| 亚洲欧美精品自产自拍| 一级毛片 在线播放| 亚洲,一卡二卡三卡| 午夜激情av网站| 国产又爽黄色视频| 亚洲av电影在线进入| 国产成人精品无人区| 亚洲综合色网址| 天堂8中文在线网| 不卡av一区二区三区| 美女视频免费永久观看网站| 国产主播在线观看一区二区 | netflix在线观看网站| 97精品久久久久久久久久精品| 自线自在国产av| 99re6热这里在线精品视频| cao死你这个sao货| 亚洲国产av新网站| 久久综合国产亚洲精品| 亚洲欧美激情在线| 亚洲,欧美精品.| 欧美日韩成人在线一区二区| 国产精品一区二区在线不卡| 亚洲午夜精品一区,二区,三区| 日韩制服丝袜自拍偷拍| h视频一区二区三区| 亚洲欧美一区二区三区黑人| 菩萨蛮人人尽说江南好唐韦庄| av天堂久久9| 久久久精品区二区三区| 免费在线观看完整版高清| 日韩人妻精品一区2区三区| 一边亲一边摸免费视频| 天天躁日日躁夜夜躁夜夜| 男女午夜视频在线观看| 国产一区二区激情短视频 | 美女大奶头黄色视频| 18禁黄网站禁片午夜丰满| 免费一级毛片在线播放高清视频 | 这个男人来自地球电影免费观看| 国产欧美日韩综合在线一区二区| 久久久久久久久免费视频了| 国产亚洲一区二区精品| 欧美激情 高清一区二区三区| 国产女主播在线喷水免费视频网站| 国产成人免费观看mmmm| 我要看黄色一级片免费的| 无遮挡黄片免费观看| 亚洲av日韩精品久久久久久密 | 色婷婷久久久亚洲欧美| 男女无遮挡免费网站观看| 99国产综合亚洲精品| 国产成人91sexporn| 天天躁日日躁夜夜躁夜夜| 涩涩av久久男人的天堂| 国产成人精品久久久久久| 性少妇av在线| av电影中文网址| 午夜影院在线不卡| 99re6热这里在线精品视频| 在线观看免费午夜福利视频| 高潮久久久久久久久久久不卡| 妹子高潮喷水视频| a级片在线免费高清观看视频| 男女免费视频国产| e午夜精品久久久久久久| 一本久久精品| 国产在线观看jvid| 免费高清在线观看日韩| 久久免费观看电影| 少妇的丰满在线观看| 久久精品国产综合久久久| 晚上一个人看的免费电影| 蜜桃在线观看..| 男女边摸边吃奶| 国产亚洲欧美精品永久| 一本久久精品| 久久久国产欧美日韩av| 国产亚洲精品久久久久5区| 成人三级做爰电影| av不卡在线播放| 久久精品aⅴ一区二区三区四区| 一级片免费观看大全| 国产成人啪精品午夜网站| 亚洲欧美日韩另类电影网站| 午夜福利,免费看| 久久99精品国语久久久| 青春草亚洲视频在线观看| 国产精品一区二区精品视频观看| 精品一区二区三卡| 亚洲欧美一区二区三区黑人| 婷婷丁香在线五月| 少妇 在线观看| 国产av一区二区精品久久| 男的添女的下面高潮视频| 日本猛色少妇xxxxx猛交久久| 老司机亚洲免费影院| 天天影视国产精品| 欧美人与善性xxx| 亚洲精品日韩在线中文字幕| 18禁裸乳无遮挡动漫免费视频| 丝袜美足系列| 亚洲欧美激情在线| 国产主播在线观看一区二区 | 午夜91福利影院| av线在线观看网站| 国产av一区二区精品久久| 狠狠婷婷综合久久久久久88av| 两人在一起打扑克的视频| av国产久精品久网站免费入址| 天天躁夜夜躁狠狠躁躁| www.av在线官网国产| 日本色播在线视频| 久久国产精品人妻蜜桃| 日韩中文字幕欧美一区二区 | av又黄又爽大尺度在线免费看| 日韩av不卡免费在线播放| 制服诱惑二区| 黄色a级毛片大全视频| 精品国产超薄肉色丝袜足j| 亚洲视频免费观看视频| 久久久亚洲精品成人影院| av又黄又爽大尺度在线免费看| 好男人视频免费观看在线| 国产成人系列免费观看| 亚洲七黄色美女视频| 丝袜喷水一区| 国产淫语在线视频| 久久人妻福利社区极品人妻图片 | 一级毛片我不卡| 新久久久久国产一级毛片| 亚洲av电影在线进入| 女警被强在线播放| 久久久久久久精品精品| 丝袜脚勾引网站| 男女午夜视频在线观看| 午夜福利在线免费观看网站| 少妇人妻 视频| 9191精品国产免费久久| 天天躁夜夜躁狠狠躁躁| 久久精品人人爽人人爽视色| 18在线观看网站| 高潮久久久久久久久久久不卡| 超碰97精品在线观看| 性高湖久久久久久久久免费观看| 两人在一起打扑克的视频| 免费日韩欧美在线观看| 精品一品国产午夜福利视频| 亚洲激情五月婷婷啪啪| 国产高清videossex| 国产视频一区二区在线看| 欧美变态另类bdsm刘玥| 黄频高清免费视频| 国产成人免费观看mmmm| 日日摸夜夜添夜夜爱| 一区二区三区激情视频| 国产精品国产av在线观看| 中国美女看黄片| 又粗又硬又长又爽又黄的视频| 国产老妇伦熟女老妇高清| 午夜影院在线不卡| 男女下面插进去视频免费观看| 美女扒开内裤让男人捅视频| 国产精品99久久99久久久不卡| 一区福利在线观看| 国产黄频视频在线观看| 亚洲精品久久午夜乱码| 自线自在国产av| 少妇裸体淫交视频免费看高清 | 国产一区二区三区av在线| 中国国产av一级| 青青草视频在线视频观看| 欧美成人精品欧美一级黄| 欧美黄色片欧美黄色片| 只有这里有精品99| 久久久久精品国产欧美久久久 | 啦啦啦视频在线资源免费观看| 91字幕亚洲| 丝袜美腿诱惑在线| 男女高潮啪啪啪动态图| 国产无遮挡羞羞视频在线观看| 在线精品无人区一区二区三| 别揉我奶头~嗯~啊~动态视频 | 免费av中文字幕在线| bbb黄色大片| 久久久久久久精品精品| 免费少妇av软件| 日本色播在线视频| 成年人午夜在线观看视频| 欧美激情极品国产一区二区三区| 99久久精品国产亚洲精品| 一级黄片播放器| 天天添夜夜摸| av网站在线播放免费| 色婷婷久久久亚洲欧美| 男女午夜视频在线观看| 大话2 男鬼变身卡| 一级片'在线观看视频| 中文字幕精品免费在线观看视频| 精品亚洲成国产av| 性色av乱码一区二区三区2| 爱豆传媒免费全集在线观看| 亚洲国产最新在线播放| 在线亚洲精品国产二区图片欧美| 91精品三级在线观看| 大香蕉久久成人网| 美女福利国产在线| 91成人精品电影| 国产高清不卡午夜福利| 老司机影院毛片| 欧美另类一区| 久久精品国产a三级三级三级| 日韩人妻精品一区2区三区| 久久精品久久精品一区二区三区| 亚洲第一av免费看| 最新的欧美精品一区二区| 久久 成人 亚洲| 国产精品久久久av美女十八| 大陆偷拍与自拍| 国产成人欧美在线观看 | 日本欧美国产在线视频| 国产淫语在线视频| 成年人黄色毛片网站| 国产1区2区3区精品| 精品亚洲成a人片在线观看| 90打野战视频偷拍视频| 久久久久视频综合| 久热爱精品视频在线9| 亚洲成人国产一区在线观看 | 亚洲伊人色综图| 亚洲欧美一区二区三区黑人| 国产精品偷伦视频观看了| 日本午夜av视频| 久久精品久久精品一区二区三区| 欧美av亚洲av综合av国产av| 男女之事视频高清在线观看 | 精品少妇久久久久久888优播| 久久精品国产a三级三级三级| 国产有黄有色有爽视频| 大码成人一级视频| 精品福利观看| kizo精华| 亚洲熟女精品中文字幕| 亚洲av片天天在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 日韩av不卡免费在线播放| 日日爽夜夜爽网站| 黑人巨大精品欧美一区二区蜜桃| 久久精品aⅴ一区二区三区四区| 一区福利在线观看| 欧美日韩一级在线毛片| 一级黄色大片毛片| 色网站视频免费| 亚洲情色 制服丝袜| 国产成人系列免费观看| 日韩伦理黄色片| 性色av一级| 丝袜美腿诱惑在线| 狂野欧美激情性xxxx| 亚洲欧美日韩高清在线视频 | 午夜福利视频在线观看免费| 人妻一区二区av| 菩萨蛮人人尽说江南好唐韦庄| 欧美人与性动交α欧美软件| 亚洲av成人不卡在线观看播放网 | 午夜免费男女啪啪视频观看| 成人手机av| 日韩免费高清中文字幕av| av欧美777| 亚洲精品日韩在线中文字幕| 国产精品一区二区免费欧美 | 午夜激情av网站| 亚洲少妇的诱惑av| 中国美女看黄片| 午夜福利视频精品| avwww免费| 亚洲精品日韩在线中文字幕| 校园人妻丝袜中文字幕| 国产精品一二三区在线看| 久久精品国产亚洲av涩爱| 可以免费在线观看a视频的电影网站| 999精品在线视频| 日韩一本色道免费dvd| 欧美 亚洲 国产 日韩一| 捣出白浆h1v1| 久久人人爽av亚洲精品天堂| 成年动漫av网址| 亚洲欧洲国产日韩| 性色av一级| 人成视频在线观看免费观看| 男女边吃奶边做爰视频| 各种免费的搞黄视频| 一边亲一边摸免费视频| 日韩中文字幕欧美一区二区 | 亚洲av成人精品一二三区| 亚洲专区国产一区二区| 亚洲人成77777在线视频| 9色porny在线观看| 女人高潮潮喷娇喘18禁视频| 日韩一本色道免费dvd| 国产成人av激情在线播放| 亚洲av成人不卡在线观看播放网 | 日韩一本色道免费dvd| 亚洲,一卡二卡三卡| 精品国产一区二区三区四区第35| 美女主播在线视频| 国产精品久久久久成人av| 国产在线视频一区二区| 成人三级做爰电影| 欧美成人精品欧美一级黄| 男的添女的下面高潮视频| 91精品伊人久久大香线蕉| 久久亚洲精品不卡| 久久精品久久久久久噜噜老黄| 亚洲五月色婷婷综合| 大陆偷拍与自拍| 国产av国产精品国产| 亚洲中文av在线| av国产久精品久网站免费入址| 99国产综合亚洲精品| 免费日韩欧美在线观看| 丝袜人妻中文字幕| 午夜视频精品福利| 日韩制服骚丝袜av| 天天影视国产精品| 欧美大码av| 国产一区二区三区综合在线观看| netflix在线观看网站| 人体艺术视频欧美日本| www日本在线高清视频| 交换朋友夫妻互换小说| 99热全是精品| 国产极品粉嫩免费观看在线| videos熟女内射| 久久久国产一区二区| 80岁老熟妇乱子伦牲交| 看免费av毛片| 丝袜美腿诱惑在线| 一级a爱视频在线免费观看| 国产色视频综合| 91字幕亚洲| 成人国产av品久久久| 亚洲国产中文字幕在线视频| 日韩 欧美 亚洲 中文字幕| 久久精品亚洲av国产电影网| 亚洲五月婷婷丁香| 欧美精品亚洲一区二区| 国产免费福利视频在线观看| 9色porny在线观看| 美女午夜性视频免费| 色综合欧美亚洲国产小说| 国产在线一区二区三区精| 精品熟女少妇八av免费久了| 成年人午夜在线观看视频| 成人亚洲欧美一区二区av| 国产成人欧美在线观看 | 搡老岳熟女国产| 国产精品亚洲av一区麻豆| 极品少妇高潮喷水抽搐| 一区二区av电影网| 晚上一个人看的免费电影| 精品福利永久在线观看| 在线 av 中文字幕| 免费高清在线观看日韩| 脱女人内裤的视频| 91麻豆精品激情在线观看国产 | 夫妻性生交免费视频一级片| 欧美日韩视频高清一区二区三区二| 午夜免费男女啪啪视频观看| 高清av免费在线| 丝袜喷水一区| 99热网站在线观看| 精品亚洲成a人片在线观看| 又大又黄又爽视频免费| 国产成人欧美| 久9热在线精品视频| 女人被躁到高潮嗷嗷叫费观| 亚洲人成电影免费在线| 欧美亚洲日本最大视频资源| 欧美日韩成人在线一区二区| 国产亚洲av高清不卡| 亚洲黑人精品在线| 人成视频在线观看免费观看| av不卡在线播放| 免费观看av网站的网址| 中文字幕高清在线视频| 国产在线免费精品| 亚洲色图综合在线观看| 亚洲,欧美,日韩| 中文字幕最新亚洲高清| 国产午夜精品一二区理论片| 国产男女内射视频| 天天躁狠狠躁夜夜躁狠狠躁| 首页视频小说图片口味搜索 | 亚洲国产av新网站| 中文字幕亚洲精品专区| 国产一区二区三区av在线| 国产精品久久久久久精品古装| 色婷婷av一区二区三区视频| 午夜免费男女啪啪视频观看| 99久久综合免费| 亚洲国产精品一区二区三区在线| 黄网站色视频无遮挡免费观看| 精品一区二区三区四区五区乱码 | 亚洲免费av在线视频| 国产午夜精品一二区理论片| 精品少妇内射三级| 欧美精品人与动牲交sv欧美| 国产欧美日韩一区二区三区在线| 一级片'在线观看视频| www.999成人在线观看| 婷婷色综合www| 免费看十八禁软件| 久久久久久久大尺度免费视频| 欧美黄色淫秽网站| 国产一卡二卡三卡精品| 丰满少妇做爰视频| 超色免费av| 精品人妻一区二区三区麻豆| av天堂久久9| 精品久久久久久久毛片微露脸 | 一区二区三区激情视频| 亚洲免费av在线视频| 宅男免费午夜| 精品少妇内射三级| 亚洲精品自拍成人| www.自偷自拍.com| 欧美中文综合在线视频| 日韩熟女老妇一区二区性免费视频| 色婷婷久久久亚洲欧美| 欧美人与性动交α欧美精品济南到| 亚洲国产毛片av蜜桃av| 久久久久久久久免费视频了| 成年美女黄网站色视频大全免费| 婷婷成人精品国产| 两人在一起打扑克的视频| 国产91精品成人一区二区三区 | 久久久久久久精品精品| 男女免费视频国产| 色婷婷久久久亚洲欧美| 老鸭窝网址在线观看| 亚洲欧美激情在线| 成年人黄色毛片网站| 亚洲专区中文字幕在线| 国产成人91sexporn| 欧美av亚洲av综合av国产av| 亚洲精品乱久久久久久| 国产又爽黄色视频| 自拍欧美九色日韩亚洲蝌蚪91| 午夜福利乱码中文字幕| 男男h啪啪无遮挡| 久久久亚洲精品成人影院| 悠悠久久av| 两人在一起打扑克的视频| 青春草视频在线免费观看| 国产精品欧美亚洲77777| 18禁国产床啪视频网站| 日日夜夜操网爽| 你懂的网址亚洲精品在线观看| a级毛片黄视频| 久久天躁狠狠躁夜夜2o2o | 人成视频在线观看免费观看| 精品少妇一区二区三区视频日本电影| 国产黄色免费在线视频| 婷婷色av中文字幕| 日本一区二区免费在线视频| 亚洲国产精品国产精品| 日韩中文字幕视频在线看片| 人妻一区二区av| 欧美国产精品va在线观看不卡| 久9热在线精品视频| h视频一区二区三区| 丝袜美足系列| 国产精品偷伦视频观看了| 国产又爽黄色视频| 在线观看免费高清a一片| 久久精品成人免费网站| 女人高潮潮喷娇喘18禁视频| 天堂8中文在线网| 男女下面插进去视频免费观看| 自线自在国产av| 国产高清国产精品国产三级| 叶爱在线成人免费视频播放| 国产一级毛片在线| 夫妻午夜视频| 久久女婷五月综合色啪小说| 久久精品久久久久久久性| 又粗又硬又长又爽又黄的视频| videosex国产| 韩国高清视频一区二区三区| av网站在线播放免费| 人人妻人人澡人人爽人人夜夜| 精品少妇一区二区三区视频日本电影| 亚洲成国产人片在线观看| www.av在线官网国产| 久久人人爽人人片av| 在线av久久热| 国产精品九九99| 一级黄片播放器| 大码成人一级视频| 看十八女毛片水多多多| 一级黄片播放器| 国产男人的电影天堂91| 人人妻人人澡人人爽人人夜夜| 丰满饥渴人妻一区二区三| 午夜激情久久久久久久| 国产男女超爽视频在线观看| 妹子高潮喷水视频| √禁漫天堂资源中文www| 天天添夜夜摸| 每晚都被弄得嗷嗷叫到高潮| 夫妻午夜视频|