<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Computer Modeling of Forming Surface by Sequential Wire Electrodischarge-Electrochemical Machining</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>khotnyi[</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>ksii Tim</string-name>
          <email>alec.timchenko@gmail.com</email>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>syl Osyp</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Cherkasy State Technological University</institution>
          ,
          <addr-line>Shevchenka blvd., 460, Cherkasy, 18006</addr-line>
          ,
          <country country="UA">Ukraine</country>
        </aff>
      </contrib-group>
      <abstract>
        <p>The mathematical principles of computer modeling of the surface's formation process using the combined sequential technology of electrodischarge cutting and electrochemical machining by a wire electrode are developed. The data of computer modeling of complex surface formation processes using this technology, their analysis and the results of experimental verification are presented in this article. It is proved that the use of computer modeling methods allows to improve the surface forming accuracy and the quality of the parts` surface, which are manufactured by this modern promising technology.</p>
      </abstract>
      <kwd-group>
        <kwd>surface shape accuracy</kwd>
        <kwd>wire electrical discharge machining</kwd>
        <kwd>electrochemical machining</kwd>
        <kwd>EDM/ECM combined process</kwd>
        <kwd>wire-electrode deflection</kwd>
        <kwd>ruled surface formation</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>Wire electrical discharge machining (WEDM) has become one of the most popular
processes for producing precise geometries in hard materials, such as those used in
the tooling industry. Since it is recognized as a precision process, optimization of
different aspects related to dimensional accuracy is a classic research topic [1].</p>
      <p>For WEDM the problem is the reproduction of a given ruled surface by a non-rigid
wire electrode, which is bent due to the forces caused by spark discharges. A
particularly acute problem of surface forming error occurs when thick work-pieces are cut by
a wire-cut EDM, the distance between the upper and lower guide is so wide that the
wire is easily bent.</p>
      <p>The researchers have proposed many effective strategies, including a hybrid
computer-integrated system for the improvement of the corner cutting accuracy [2-4]. At
the end of the processing of one corner`s side plane using different methods, the wire
electrode is straightened and continues the processing of the second side plane.
Unfortunately, such a two-stage strategy is not applicable when continuously cutting in
circular interpolation mode along the circle arc. The bent wire electrode forms a
barrel-shaped surface. The problem with barrel-shaped inaccuracy remains relevant, both
for traditional WEDM and for high-speed WEDM, for workpieces with a thickness
greater than 100 mm [5].</p>
      <p>One of the latest cutting technologies, presented by Japanese researchers [6], is
“High precision cutting for thick work-pieces”. The essence of the technology lies in
the development of an optimum control of the discharge energy, paying attention to
the wire bending. It can significantly reduce the wire bending and effectively improve
the straightness of thick work-pieces. The straightness of 300 mm SKD11 work-piece
is achieved to 7μm accuracy by cutting 3 times on the ROBOCUT machine.
However, the reduction of the wire deflection in this case is achieved by decreasing the
intensity of the discharges, and hence the loss of productivity and increasing of
expensive processing time. Additionally, the entire heat affected zone of WEDM in the
near-surface layers of a workpiece, which is usually detrimental, is not completely
eliminated.</p>
      <p>Significantly more promising is the use of the combined technology of sequential
electrodischarge and electrochemical machining (ECM) with a wire electrode [7].
After roughing electrodischarge cutting, follow 1–2 passes with electrolyte feed to the
processing area. The transverse forces do not act on the wire electrode, it remains
straight and as a result of the anode dissolution, the surface of the part is aligned, the
heat-affected zone is completely eliminated.</p>
      <p>An actual task for the implementation of the proposed combined technology is the
creation of mathematical methods for calculation support and modeling of the
formation of ruled surfaces by a wire electrode. Computer simulation dramatically reduces
the amount of experimental work and the accuracy degree can be acceptable for many
applications.</p>
      <p>The aim of the work is to develop and carry out testing of computer modeling
methods to improve the accuracy of surface forming and surface quality of parts
manufactured using the latest technology of combined electrodischarge and
electrochemical machining with a wire electrode.
2</p>
    </sec>
    <sec id="sec-2">
      <title>Experimental equipment, materials</title>
      <p>Experimental studies of forming processes by electroerosive wire cutting and
subsequent electrochemical machining with a wire electrode were carried out on an
electrodischarge cutting machine SELD-02 (Ukraine) with computer numerical control
(CNC) systems. The rough cutting mode of the generator GKI 300-200 А was used.
Pulses` characteristics: the nominal pulse frequency 22 kHz, pulse duration 3 µs,
current amplitude at load 0.1 Ohm 170 A. The longitudinal tension of the wire varied
within 6 - 8 N, the pressure of the working fluid was 5ꞏ105 Pa. A 0.20 mm Cobra Cut
B (AGIE, Switzerland) wire electrode of a solid brass СuZn37 was used. Pulse
electrochemical treatment was carried out by the same electrode in a 1 M NaCl aqueous
solution in a potentiostatic regime. Workpieces 40, 68, 100 mm height made of tool
stamping steel X12Ф1 (analogue DIN X155CrVMo12-1).</p>
      <p>Shape deviations of the fabricated surface were studied using the Coordinate
Measuring Machine 3D CMM AXIOM CNC ABERLINK. The limit of permissible
error is ± (1.8 + 0.4L/100) μm, where L, mm is the measured length. The software
allows you to set the value of the deviations from the cylindrical shape of the obtained
cloud of points. Surface microgeometry was studied using a TIME 3221 profilometer
and DataView TIME3R Series software.
3</p>
    </sec>
    <sec id="sec-3">
      <title>Methods</title>
      <p>3.1 Mathematical principles of computer modeling of WEDM surface forming
For two-coordinate electroerosive cutting, the wire electrode is flexed due to the
distributed transverse load caused by spark discharges along the entire height of the
workpiece. The computational scheme for determining the shape of the wire electrode
(Fig.1) is a special case of a mathematical model of the wire electrode shape in
fourcoordinate electroerosive machining given in [8]. The experimental-computational
method for determining the transverse load q(z) is presented there, depending on the
cutting regimes and the material being processed from the class of tool steels and
solid alloys.</p>
      <p>For our case, deflection magnitude of a wire electrode in an arbitrary section
(0&lt;z&lt;l) is described by the formula:
f ( z)  q  z (l  z1 )2  z</p>
      <p>2H0  l l
where we state (z – z1) = 0 at z &lt; z1 and (z – z2) = 0 at z &lt; z2.</p>
      <p>To determine the shape of the surface, which is formed by a deflected wire
electrode while moving along the arc of a circle in space, consider the scheme of
motion of individual points of the electrode (Fig. 2). Let the wire fixation nodes (point
М1(x1,y1)) moving along the arc of a given radius r1 (program trajectory). Then point
M2(x2,y2), which is at a height z of a wire electrode with a corresponding deflection

f(z), describes another arc of the radius r2. Vector M1M 2 directed at a tangent to a
 
circle at the point М1 in the opposite directionofthe velocity vector v , then M1M 2 is
perpendicular to the vector r1 . Vector length M1M 2 is equal to f(z). Then
r2 ( z) </p>
      <p>r12  f 2 ( z) .</p>
      <p>
        Or, in the form that is convenient for simulation, we obtain expressions for the
coordinates of the point М2 depending on the position in height z (z1&lt;z&lt;z2):
y x
x2 ( z)  f ( z)  r11  x1, y2 ( z)  y1  f ( z)  r11 . (
        <xref ref-type="bibr" rid="ref3">3</xref>
        )
      </p>
      <p>
        Thus, according to the established form of the wire electrode and taking into
account the correction for half the cut width, which for WEDM consists of the wire
electrode radius and the value of the interelectrode gap (IEG), the shape of the linear
surface obtained as a result of the electrode movement along the arc section of the
given trajectory is built.
(
        <xref ref-type="bibr" rid="ref1">1</xref>
        )
(
        <xref ref-type="bibr" rid="ref2">2</xref>
        )
3.2 Mathematical principles of computer modeling of wire ECM surface forming
The rate of anode dissolution, and hence the change of shape, depends on the density
of the technological current on the surface of the part. For an electrochemical cell
with the following configuration of electrodes, a cylindrical cathode – a cylindrical
anode – the law of distribution must be established of the electric field strength in the
IEG and the current density on the anode surface. The calculation scheme is presented
in Fig.3. Using the method of conformal mappings of the complex variable function
theory [9], we obtain the expressions for the intensity of the electric field:
 x
Ex  A x2  ( y  a)2

      </p>
      <p>x   y  a
x2  ( y  a)2  , Ey  A x2  ( y  a)2
</p>
      <p>y  a 
x2  ( y  a)2  ,
E </p>
      <p>
        Ex2  E y2 ,
(
        <xref ref-type="bibr" rid="ref4">4</xref>
        )
where a  b2  2bR , b 
2(δ  r0  R)
      </p>
      <p>δ(δ  2r0 ) , A  φl1npφ1 0 , p1  aa  bb  aa  δδ  bb ,
here R – radius of the cylindrical surface of the anode, r0 – radius of the cathode, δ
– value of the interelectrode gap, φ0, φ1 – potentials on the cathode and the anode,
respectively.
where s – arc coordinate, which is calculated from point A’ (Fig.3).</p>
      <p>The maximum value of the electric field strength:</p>
      <p>Emax  E</p>
      <p>s  0
jmax 
 2(φ1  φ0 ) .</p>
      <p>a ln p1
2κ(φ1  φ0 ) ,</p>
      <p>a ln p1
The maximum value of current density at the anode`s surface at the point A’: 
where κ – specific electrolyte conductivity.</p>
      <p>It is important to determine the thickness of the dissolved surface layer to
simulate surface formation. To do this, consider the rotational motion of the wire
cathode around the anode with the angular velocity ω (Fig. 3). The linear velocity of
the cathode`s movement in a given trajectory is determined by the expression:
v  ω  (R  δ  r0 ) . (8)
Then the intensity of the electric field on the anode`s surface will take the form:
2a
E  A
(s  R  ω  t)2  a2
.</p>
      <p>By Faraday’s law, the depth of dissolution at the point A’ during time dt:
d   Kv j
s  0</p>
      <p>dt ,
where tk =π/ω. Or
where Kv – coefficient of electrochemical dissolution of anode material, η – current
ratio efficiency, j – current density.</p>
      <p>After integrating (10) in the corresponding boundaries, we find the thickness of the
dissolved layer in one processing run:
m tk dt
 d   4Kv κAa , (11)
0 0 (R  ω  t)2  a2</p>
      <p>1  R
m  4Kv κA arctg( ) .</p>
      <p>Rω a</p>
      <p>For potentiostatic pulse regime of ECM at φ1 = 0, φ0 = -U, U – voltage, we finally
get:
(12)
m  4Kv κ</p>
      <p>U ti arctg(
R  ω  ln( p1) tn
 R
a
) ,
where ti – pulse on-time, tn – period of pulses.</p>
      <p>
        After the electroerosive cutting on the arc of radius R with the deformed wire
electrode, the resulting surface will have a convexity in the middle of the workpiece along
the z axis (Fig. 4). In the case of the subsequent ECM using nondeformed wire
electrode, minimal and maximal IEGs will be δ1 and δ2, respectively. From the formulas`
analysis (
        <xref ref-type="bibr" rid="ref4">4</xref>
        ), (7), (13) it follows that the intensity of the electric field, the current
density, and, as a result, the intensity of the anode dissolution increases with the decrease
of IEG. Therefore, in the middle of the workpiece material removal will be bigger and
the process will tend to leveling the surface in height. That is, when wire moves in a
circle with radius R + δ + r0 there will be removed thicknesses Δm1, Δm2, considering
Δm1&gt; Δm2, as δ1 &lt;δ2.
(6)
(7)
(9)
(10)
(13)
      </p>
      <p>A convex surface will become flat if equality is satisfied:</p>
      <p>δ1  m1  δ2  m2 .</p>
      <p>
        We find the thickness of the dissolved layer in accordance with (13), then formula
(14) will take form:
(14)
(15)
Based on the mathematical dependences obtained, a series of computational
experiments on computer simulation of surface forming processes using the WEDM and
ECM with a wire electrode with the experimental verification of individual results
was performed using the Mathcad package. Computational experiments have shown
that the wire electrode deflection is significantly increased with increasing distance
between fixation points – electrode guides. The parameter l is included in (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) at the
power 2. This explains the emergence of the unacceptable surface formation error
problem in the processing of thick parts.
      </p>
      <p>Computer modeling of shaping of ruled surfaces (Fig. 5) using WEDM by
deformed wire electrode during cutting in circular interpolation regime revealed a
significant dependence of the error from the trajectory arc radius. The determinant factor
is the relationship between the thickness of the workpiece and the magnitude of the
arc radius – the generating line of the cylindrical surface.
tion error of the cylinders with a radius less than 0.5 mm becomes significant (curve
2). Curve 3 shows that surface formation error, due to the wire electrode deflection,
will exceed 20 microns when cutting a cylindrical part with a radius less than 1 mm in
diameter. With an increase in the thickness of the workpiece, the formation error will
increase dramatically. The theoretically calculated curve 4 for the thickness of the
workpiece of 150 mm, shows, that for radii of 5 mm and less, error becomes
unacceptably large.</p>
      <p>The following simulation is performed for the aforementioned workpiece materials
and WEDM cutting regime. Experimental verification of computer simulation with
the coordinate-measuring machine control of the convexity of the obtained cylindrical
surfaces, confirmed the satisfactory accuracy of the model. The maximum measured
difference between the results of the computational and physical experiment does not
exceed 8% (Fig.6).</p>
      <p>For other WEDM technological regimes, it is advisable to use the method of
determining the power factors influencing the wire electrode described in [8] and to
conduct a similar simulation according to the proposed method.</p>
      <p>
        The task of the next stage of computer simulation of the cylindrical surface
aligning process using the wire ECM technology is based on a reasonable choice of the
amplitude-time parameters of the pulsed current and the velocity of the electrodes`
movement along the trajectory. Theoretical calculations, simulations and particular
experimental studies have shown that already at the beginning, anodic dissolution
predominantly aligns the difference between the convex part in the middle of the
height of the cylindrical surface and the less convex regions from above and below
(Δm1 and Δm2 respectively, Fig. 4). Then, with the time of dissolution, the difference
between Δm1 and Δm2 decreases, and more intense removal of the thickness of the
entire surface layer becomes more predominant than the surface leveling in terms of
the height of the workpiece. Indeed, from formulas (
        <xref ref-type="bibr" rid="ref4">4</xref>
        ), (13) it can be seen that the
dissolution rate depends on the intensity of the electric field, which is higher on the
humps. Over time, the humps become smaller and the unevenness of the distribution
of tension is leveled. In practice, this means that it is not necessary to continue the
process until the full alignment of the surface, but it is advisable to stop when
achieving the required precision of reproduction of geometric parameters (Fig.5). This once
again proves the importance of computer simulation in the design process of the
combined technology. Its use allows to predict the magnitude of the surface formation
error obtained after WEDM and, accordingly, the amplitude-time parameters of the
pulsed current and the velocity of the electrode movement along the trajectory to
reduce deviation level to the acceptable result using subsequent ECM with a wire
electrode.
      </p>
      <p>Increasing the intensity of the predominant dissolution of humps is possible by
reducing the gap between the wire electrode and the surface of the part. The electric
field strength concentration increase becomes more noticeable on the humps
compared with the caverns. However, setting the IEG of less than 0.1 mm is not allowed
due to disruption of a uniform electrolyte flow throughout the height of the part. The
uneven and insufficient renewal of the electrolyte in the intercolectrode space leads to
a decrease in the current efficiency (η in formulas (10-13), (15)), which significantly
influences surface shaping and, accordingly, the accuracy of its modeling.</p>
      <p>In addition to increasing accuracy, due to electrochemical machining with a wire
electrode, surface quality is also improved. The heat affected zone, obtained due to
spark discharge, is completely removed. Surface roughness is reduced.
Experimentally obtained decrease in roughness from Ra 2.7 – 3.2 μm after WEDM to Ra 0.65 –
1.2 μm after the ECM stage. It should be noted that the possibilities of improving
surface properties by finishing electrochemical machining, in particular, reducing the
surface roughness, are not limited to the results obtained in this paper. By choosing
the amplitude-time characteristics of pulses, the electrolyte composition, modern
technologies of electrochemical surface polishing achieve even more significant
results [10,11].
5</p>
      <p>Conclusion
1. There are created mathematical principles of computer simulation of the process of
surface forming using combined WEDM and ECM with a wire electrode.</p>
      <p>2. The implemented complex of computational experiments allowed broadening
the idea of complex surface forming processes using the technology of a combined
sequential WEDM and ECM with a wire electrode. Experimental verification of
computer simulation results confirmed the adequacy of the proposed mathematical
models.</p>
      <p>3. Usage of computer simulation technologies allowed to improve the shaping
accuracy and surface quality of the parts manufactured with the combined technology
WEDM and ECM with a wire electrode.</p>
      <p>4. The obtained results are an element of the scientific basis for the further
development of the modern promising technology.
The International Journal of Advanced Manufacturing Technology, vol.
74(58), pp.917-929 (2014) doi:10.1007/s00170-014-6002-5
5. Deng, C., Liu, Z., Zhang, M., Ji, Y., Pan, H.: Minimizing drum-shaped
inaccuracy in high-speed wire electrical discharge machining after multiple cuts. The
International Journal of Advanced Manufacturing Technology, pp. 1-11
(2019) doi:10.1007/s00170-018-3068-5
6. Takayama, Y., Makino, Y., Niu, Y., Uchida, H.: The Latest Technology of
Wire-cut EDM. Procedia CIRP, vol. 42, pp. 623-626 (2016)
doi:10.1016/j.procir.2016.02.259
7. Osipenko, V.I., Plakhotny, A.P., Denisenko, A.Yu.: Improved methodology
for calculating the processes of surface anodic dissolution of spark eroded
recast layer at electrochemical machining with wire electrode. Pratsi Odeskoho
politekhnichnoho universytetu, vol. 1(43), pp.55-60 (2014)
8. Osypenko, V., Plakhotnyi, O.: Determination of force factors influencing
surface forming in wire electrical discharge machining. Elektronnaya obrabotka
materialov, vol. 54(6), pp.9-18 (2018) (in Russian) doi:
10.5281/zenodo.1968973. The English version of the journal: Surface
Engineering and Applied Electrochemistry.
9. Zhytnykov, V.P., Zaitsev, A.N.: Pulse precision electrochemical machining.</p>
      <p>
        Mashynostroenye, Moscow (2008) (in Russian).
10. Taylor, E.J., Inman, M.: Electrochemical Surface Finishing. The
Electrochemical Society Interface, vol. 23(
        <xref ref-type="bibr" rid="ref3">3</xref>
        ), pp. 57-61 (2014)
doi:10.1149/2.F05143if
11. Yang, G., Wang, B., Tawfiq, K., Wei, H., Zhou, S., Chen, G.: Electropolishing
of surfaces: theory and applications. Surface Engineering, vol. 33(
        <xref ref-type="bibr" rid="ref2">2</xref>
        ), pp.
149166 (2016) doi:10.1080/02670844.2016.1198452
      </p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Sanchez</surname>
            ,
            <given-names>J.A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Plaza</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          , Lopez de Lacalle,
          <string-name>
            <given-names>L.</given-names>
            ,
            <surname>Lamikiz</surname>
          </string-name>
          ,
          <string-name>
            <surname>A.</surname>
          </string-name>
          :
          <article-title>Computer simulation of wire-EDM taper-cutting</article-title>
          .
          <source>International Journal of Computer Integrated Manufacturing</source>
          , vol.
          <volume>19</volume>
          (
          <issue>7</issue>
          ), pp.
          <fpage>727</fpage>
          -
          <lpage>735</lpage>
          (
          <year>2006</year>
          ) doi:10.1080/09511920600628855
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <surname>Sanchez</surname>
            ,
            <given-names>J.A.</given-names>
          </string-name>
          , Lopez de Lacalle,
          <string-name>
            <given-names>L.N.</given-names>
            ,
            <surname>Lamikiz</surname>
          </string-name>
          ,
          <string-name>
            <surname>A.</surname>
          </string-name>
          :
          <article-title>A computer-aided system for the optimization of the accuracy of the wire electro-discharge machining process</article-title>
          .
          <source>International Journal of Computer Integrated Manufacturing</source>
          , vol.
          <volume>17</volume>
          (
          <issue>5</issue>
          ), pp.
          <fpage>413</fpage>
          -
          <lpage>420</lpage>
          (
          <year>2004</year>
          ) doi:10.1080/09511920310001626590
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <surname>Han</surname>
            ,
            <given-names>F.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Zhang</surname>
          </string-name>
          , J.,
          <string-name>
            <surname>Soichiro</surname>
            <given-names>I.</given-names>
          </string-name>
          :
          <article-title>Corner error simulation of rough cutting in wire EDM</article-title>
          .
          <article-title>Precision engineering</article-title>
          , vol.
          <volume>31</volume>
          (
          <issue>4</issue>
          ), pp.
          <fpage>331</fpage>
          -
          <lpage>336</lpage>
          (
          <year>2007</year>
          ) doi:10.1016/j.precisioneng.
          <year>2007</year>
          .
          <volume>01</volume>
          .005
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4.
          <string-name>
            <surname>Chen</surname>
            ,
            <given-names>Z.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Huang</surname>
            ,
            <given-names>Y.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Zhang</surname>
            ,
            <given-names>Z.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Li</surname>
            ,
            <given-names>H.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ming</surname>
            ,
            <given-names>W.</given-names>
          </string-name>
          , Zhang, G.:
          <article-title>An analysis and optimization of the geometrical inaccuracy in WEDM rough corner cutting</article-title>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>