<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing with OASIS Tables v3.0 20080202//EN" "journalpub-oasis3.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:oasis="http://docs.oasis-open.org/ns/oasis-exchange/table" xml:lang="en" dtd-version="3.0">
  <front>
    <journal-meta><journal-id journal-id-type="publisher">MS</journal-id><journal-title-group>
    <journal-title>Mechanical Sciences</journal-title>
    <abbrev-journal-title abbrev-type="publisher">MS</abbrev-journal-title><abbrev-journal-title abbrev-type="nlm-ta">Mech. Sci.</abbrev-journal-title>
  </journal-title-group><issn pub-type="epub">2191-916X</issn><publisher>
    <publisher-name>Copernicus Publications</publisher-name>
    <publisher-loc>Göttingen, Germany</publisher-loc>
  </publisher></journal-meta>
    <article-meta>
      <article-id pub-id-type="doi">10.5194/ms-10-331-2019</article-id><title-group><article-title>Designing optimal automatic cycles of round grinding based on the synthesis of digital twin technologies <?xmltex \hack{\break}?>and dynamic programming method</article-title><alt-title>Designing optimal automatic cycles of round grinding</alt-title>
      </title-group><?xmltex \runningtitle{Designing optimal automatic cycles of round grinding}?><?xmltex \runningauthor{P.~P. Pereverzev et al.}?>
      <contrib-group>
        <contrib contrib-type="author" corresp="no" rid="aff1">
          <name><surname>Pereverzev</surname><given-names>Pavel P.</given-names></name>
          
        <ext-link>https://orcid.org/0000-0002-1416-8918</ext-link></contrib>
        <contrib contrib-type="author" corresp="no" rid="aff1">
          <name><surname>Akintseva</surname><given-names>Aleksandra V.</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="no" rid="aff1 aff2">
          <name><surname>Alsigar</surname><given-names>Masar K.</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="yes" rid="aff1">
          <name><surname>Ardashev</surname><given-names>Dmitrii V.</given-names></name>
          <email>ardashevdv@susu.ru</email>
        </contrib>
        <aff id="aff1"><label>1</label><institution>Faculty of Mechanical Engineering, South Ural State University,
Chelyabinsk, 454080, Russian Federation</institution>
        </aff>
        <aff id="aff2"><label>2</label><institution>Thi-Qar university, College of engineering, Nasiriyah, Republic of
Iraq</institution>
        </aff>
      </contrib-group>
      <author-notes><corresp id="corr1">Dmitrii V. Ardashev (ardashevdv@susu.ru)</corresp></author-notes><pub-date><day>1</day><month>July</month><year>2019</year></pub-date>
      
      <volume>10</volume>
      <issue>1</issue>
      <fpage>331</fpage><lpage>341</lpage>
      <history>
        <date date-type="received"><day>5</day><month>March</month><year>2019</year></date>
           <date date-type="rev-recd"><day>11</day><month>June</month><year>2019</year></date>
           <date date-type="accepted"><day>14</day><month>June</month><year>2019</year></date>
      </history>
      <permissions>
        <copyright-statement>Copyright: © 2019 Pavel P. Pereverzev et al.</copyright-statement>
        <copyright-year>2019</copyright-year>
      <license license-type="open-access"><license-p>This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this licence, visit <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">https://creativecommons.org/licenses/by/4.0/</ext-link></license-p></license></permissions><self-uri xlink:href="https://ms.copernicus.org/articles/10/331/2019/ms-10-331-2019.html">This article is available from https://ms.copernicus.org/articles/10/331/2019/ms-10-331-2019.html</self-uri><self-uri xlink:href="https://ms.copernicus.org/articles/10/331/2019/ms-10-331-2019.pdf">The full text article is available as a PDF file from https://ms.copernicus.org/articles/10/331/2019/ms-10-331-2019.pdf</self-uri>
      <abstract><title>Abstract</title>
    <p id="d1e112">The article presents methodology for designing optimal
feed control cycles resistant to unstable machining conditions for a batch
of parts in round grinding operations performed on CNC machines.</p>
    <p id="d1e115">To improve the quality and reliability of control programs for CNC machines, a digital twin (DT) is proposed to be used. It performs virtual testing of a given grinding cycle for the possibility of defect occurrence at some
combination of variable technological factors. To design the optimal
grinding cycle by calculating the optimum trajectory of the radial feed
change cycle, we propose the dynamic programming method (DPM) used in
solving the classical transport problem, in which a network of roads with
intermediate stations is specified. In the task it is necessary to find the
optimal trajectory of transport. The conditions for the appearance of
defects detected by the DT are introduced into the system of restrictions,
and the feed cycle is again optimized in the DPM optimization system. A new
optimized cycle is again tested by the DT. These iterations are repeated
until the DT fixes the present of a defect for the operation. The
methodology proposed in the article for the synthesis of digital twin and
DPM technologies at the stage of preparation of control programs for CNC
machines makes it possible to guarantee the maximum productivity of the
grinding operation while ensuring the specified quality of the machined
surface under varying processing conditions that vary within the specified
limits.</p>
  </abstract>
    </article-meta>
  </front>
<body>
      

<sec id="Ch1.S1" sec-type="intro">
  <label>1</label><title>Introduction</title>
      <p id="d1e127">To perform grinding operation, all you need is a control program, which
manages the execution of user program to prevent errors and cutting
parameters of the automatic steps cycle in CNC programming radial feed. The
program feed cycle operates in an automatic mode according to the commands
of the active monitoring device and switches the feed-rate depending on the
remaining part of the allowance. Designing grinding cycles for control
programs is an important stage in the development of the technological
process for machining parts in round grinding operations. The quality of
their development depends on the performance of the operation, the full
application of the technological capabilities of the machine, the amount of
defective products.</p>
      <p id="d1e130">However, in automated engineering, there are still no CAM-systems which for
the given grinding conditions would enable us to calculate the optimal
automatic cycle of a step change in feed rate and other operating parameters
that provide, in the aggregate, a grinding cycle.</p>
      <p id="d1e133">The control parameters and the relationship with the cutting parameters with
the productivity of the operation and the quality of the surface to be
treated have been studied by Malkin (1981), Fang et al. (2016), Raoufinia et al. (2019). As follows from these works, as well as the works of Uhlmann et
al. (2016), Pereira et al. (2009), Liu et al. (2017) and Rowe and Ebbrell (2004),
programming the flow parameters of the feed loop is a difficult task,
because it is necessary to set the number of stages of the cycle, the feed
amount at each stage of the cycle, and the distribution of the allowance
along the cycle steps. In addition to step-by-step feed cycles, all other
cutting parameters are programmed in the control program which includes the
speed of the grinding wheel and the<?pagebreak page332?> workpiece, the tool characteristics, its
overrun and other parameters depending on the type of grinding.</p>
      <p id="d1e136">When designing an operation for a CNC machine, the greatest difficulties
occur at the stage when the technologist assigns cutting parameters and
projects a feed switching cycle, depending on the remaining part of the
allowance. In the process of designing the grinding operation, the
technologist should strive to design optimal cycle parameters that will
ensure the maximum productivity of machining the batch of parts, taking into
account the following technological factors and requirements noted in the
works of Petrakov and Chamata (2015), Pereverzev and Akintseva (2017), Maris et al. (1975), Lin et al. (2009) and Nathan et al. (2001):
<list list-type="bullet"><list-item>
      <p id="d1e141">requirements for the accuracy of the treated workpiece surface, specified in its drawing;</p></list-item><list-item>
      <p id="d1e145">presence of limitations on the parameters of the machine;</p></list-item><list-item>
      <p id="d1e149">requirement to ensure the stability of the operation in the presence of
variable technological factors that arise when processing a batch of parts:
variation in the allowance and the dimensions of the workpiece surface,
depending on the initial precision of the billet, the error of basing and
setting, etc.;</p></list-item><list-item>
      <p id="d1e153">presence of a restriction on the strength of the cutting tool (according
to the critical force causing the scattering of the grinding wheel),
depending on the characteristics of the grinding wheel.</p></list-item></list></p>
      <p id="d1e157">The analysis of the above-mentioned requirements for designing the mode part
of the control program for a CNC machine demonstrates that the task of
designing the maximum productive cycle is reduced to the task of designing
the optimal grinding cycle in terms of speed, taking into account all the
restrictions imposed. Such a problem cannot be solved without using models
of the grinding process and mathematical methods for optimizing the cycle.</p>
      <p id="d1e160">A literature review surveys books, scientific articles and other sources of
relevant information, is refers to many scientists have been engaged in
studying the theory of designing optimal grinding cycles for CNC machines of
operational performance management, when designing an operation for a CNC
machine, the greatest difficulties occur at the stage when the technologist
assigns cutting modes and projects a feed switching cycle, depending on the
remaining part of the allowance. In the process of designing the grinding
operation, the technologist should strive to de-sign optimal cycle
parameters that will ensure the maximum productivity of machining the batch
of parts, taking into account the following technological factors and
requirements noted in the works of Cahill et al. (2015), Lur'e (1979),
Pereverzev and Pimenov (2015), Amitay et al. (1989), Phan et al. (2011),
Pereverzev and Akintseva (2015), Akintseva and Pereverzev (2017), Alagumurthi et al. (2007), Lee et al. (2006), Krajnik et al. (2008), Bertsekas (1995), Lee (2008), Nishimura et
al. (1989), and Hashimoto and Lahoti (2004). However, none of the papers considers the improvement of the
design quality of the optimum grinding cycles in terms of their resistance
to variable processing conditions existing in real circular grinding
operations performed on CNC machines.</p>
      <p id="d1e163">Thus, enormous theoretical and practical basis in the design of optimal
cycles of round grinding is not sufficient to solve the task of designing
optimal grinding cycles which would ensure maximum operation efficiency and
stability of meeting the drawing requirements for accuracy and other quality
parameters for any most unfavorable combinations of variable grinding
conditions.</p>
      <p id="d1e166">The aim of the present work is to develop methodology for solving the
problem of designing optimal rounds of circular grinding cycles for CNC
machines that are optimal in terms of speed and stability.</p>
      <p id="d1e169">Let us consider the methodology for designing optimal grinding cycles
resistant to variable processing conditions, for example, round grinding
with radial feed.</p>
</sec>
<sec id="Ch1.S2">
  <label>2</label><title>Model of the allowance removal process in the operation of round outer grinding with radial feed</title>
      <p id="d1e180">In order to optimize the grinding cycle, we need a model of the allowance
removal process first of all during the whole processing cycle, with a
stepwise changing feed rate on the commands of the active control device.
The model of the removal process of allowance is based on two models: the
model of the cutting force and the model of the relationship between the
cutting force and elastic deformations of the technological system, with
program and actual feeds. Therefore, in the preparation of control programs
for CNC machine tools, an analytical model of cutting force is needed to
calculate the cutting conditions and grinding cycles in a wide range of
various process parameters and processing conditions. Therefore, empirical
models that have a number of significant drawbacks cannot be taken as a
model of the cutting force, namely: they have a narrowly focused range of
application, the coefficients in the empirical formulas are obtained
experimentally for the specific conditions of the operation, they do not
establish an analytical relationship between the main technological
parameters of the process.</p>
      <?pagebreak page333?><p id="d1e183">A theoretical wide-range model of cutting force which takes into account
most of the main technological factors was proposed by Pereverzev and Akintseva (2017), Pereverzev et al. (2015) and Pereverzev and Pimenov (2016). It served as a model of the radial component of the
cutting force for circular grinding. The model is obtained by summing the
cutting forces of individual grains in the cutting area of the grinding
wheel with the workpiece and by applying the principle of equality of the
working cutting forces and resistance forces of the processed metal to
plastic deformation, formulated by Korchak (1974):
          <disp-formula id="Ch1.E1" content-type="numbered"><label>1</label><mml:math id="M1" display="block"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo><mml:mn mathvariant="normal">9</mml:mn><mml:mi mathvariant="italic">σ</mml:mi><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mi mathvariant="italic">π</mml:mi><mml:mtext>d</mml:mtext><mml:msub><mml:mi>n</mml:mi><mml:mtext>det </mml:mtext></mml:msub><mml:mi>B</mml:mi><mml:msub><mml:mi>S</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow><mml:mi>V</mml:mi></mml:mfrac></mml:mstyle><mml:mo>+</mml:mo><mml:mi mathvariant="italic">η</mml:mi><mml:mi>B</mml:mi><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mi mathvariant="italic">σ</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:mfrac></mml:mstyle><mml:msqrt><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mtext>d</mml:mtext><mml:mi>D</mml:mi><mml:msub><mml:mi>S</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>n</mml:mi><mml:mtext>det </mml:mtext></mml:msub><mml:mo>(</mml:mo><mml:mi>D</mml:mi><mml:mo>+</mml:mo><mml:mi>d</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mfrac></mml:mstyle></mml:msqrt><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        where <inline-formula><mml:math id="M2" display="inline"><mml:mi>B</mml:mi></mml:math></inline-formula> – width of the working surface of the wheel, mm; <inline-formula><mml:math id="M3" display="inline"><mml:mi>d</mml:mi></mml:math></inline-formula> – diameter of the
machined surface, mm; <inline-formula><mml:math id="M4" display="inline"><mml:mi>D</mml:mi></mml:math></inline-formula> – diameter of grinding wheel, mm; <inline-formula><mml:math id="M5" display="inline"><mml:mrow><mml:msub><mml:mi>n</mml:mi><mml:mtext>det</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula> –
rotational speed, RPM; <inline-formula><mml:math id="M6" display="inline"><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mtext>f</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula> – actual minute feed, mm min<inline-formula><mml:math id="M7" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>; <inline-formula><mml:math id="M8" display="inline"><mml:mi>V</mml:mi></mml:math></inline-formula> –
circumferential speed of rotation of the wheel; <inline-formula><mml:math id="M9" display="inline"><mml:mi>i</mml:mi></mml:math></inline-formula> – revolution number of the
workpiece; <inline-formula><mml:math id="M10" display="inline"><mml:mi>k</mml:mi></mml:math></inline-formula> – radius number in the section of the treated surface; <inline-formula><mml:math id="M11" display="inline"><mml:mi>z</mml:mi></mml:math></inline-formula> –
radial feed stage number; <inline-formula><mml:math id="M12" display="inline"><mml:mi mathvariant="italic">η</mml:mi></mml:math></inline-formula> – degree of blunting of the grinding
wheel; <inline-formula><mml:math id="M13" display="inline"><mml:mi mathvariant="italic">σ</mml:mi></mml:math></inline-formula> – stress intensity in the moving volume of the deformable
metal, which characterizes the resistance of the metal to the plastic flow,
kg mm<inline-formula><mml:math id="M14" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>.</p>
      <p id="d1e422">Degree of bluntness of the grinding wheel (<inline-formula><mml:math id="M15" display="inline"><mml:mi mathvariant="italic">η</mml:mi></mml:math></inline-formula>), equals to the ratio of
the total area of the blunting areas of all the grains of the grinding wheel
on its entire surface to the geometric area of the entire working surface of
the grinding wheel (that is, the value of <inline-formula><mml:math id="M16" display="inline"><mml:mi mathvariant="italic">η</mml:mi></mml:math></inline-formula> determines the relative
support surface of the grinding wheel along the areas of blunting the
grinding wheel grains).</p>
      <p id="d1e439">The presented model of cutting force for round external grinding (1) is
based on the fundamental laws of the mechanics of metal plastic deformation
in the cutting zone, takes into account the kinematics and especially the
process of round external grinding, connects the cutting forces with
machining modes (radial and axial feedrate, rotational speed of the
workpiece and the wheel, etc.), and also covers most of the main
technological parameters in a wide range of their variation – physical and
mechanical properties of the metal being ground, geometric parameters of the
contact area of the grinding wheel and the workpiece, the characteristics of
the wheel, bluntings of abrasive grains of the grinding wheel, etc.</p>
      <p id="d1e443">To establish the relationship between the cutting force <inline-formula><mml:math id="M17" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>, the
elastic deformations <inline-formula><mml:math id="M18" display="inline"><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> with the accumulated program
<inline-formula><mml:math id="M19" display="inline"><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>p</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> and factual <inline-formula><mml:math id="M20" display="inline"><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> feeds, let us consider the
calculation scheme shown in Fig. 1. The diagram shows the graphs and
dimensional chains that establish the mathematical relationship between the
elastic deformations <inline-formula><mml:math id="M21" display="inline"><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> of the technological system, the
accumulated displacements for <inline-formula><mml:math id="M22" display="inline"><mml:mi>i</mml:mi></mml:math></inline-formula> of the program <inline-formula><mml:math id="M23" display="inline"><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>p</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> and actual <inline-formula><mml:math id="M24" display="inline"><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> feeds with the current actual values feedrate <inline-formula><mml:math id="M25" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> and the current values of the radii <inline-formula><mml:math id="M26" display="inline"><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> of the machined surface at each <inline-formula><mml:math id="M27" display="inline"><mml:mi>i</mml:mi></mml:math></inline-formula>turn of the workpiece, for the case when a workpiece having the initial radial run-out is ground <inline-formula><mml:math id="M28" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>R</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F1" specific-use="star"><?xmltex \currentcnt{1}?><label>Figure 1</label><caption><p id="d1e691">Calculation scheme of the relationship between feeds,
elastic deformations and current radii at round grinding with radial feed: (1) trajectory of the radius of the part and the actual position i while
turning; (2) trajectory of the radius of the part and the actual position
<inline-formula><mml:math id="M29" display="inline"><mml:mrow><mml:mi>i</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula> while turning; (3) interrupted cutting zone; (4) programmed feed path.</p></caption>
        <?xmltex \igopts{width=341.433071pt}?><graphic xlink:href="https://ms.copernicus.org/articles/10/331/2019/ms-10-331-2019-f01.jpg"/>

      </fig>

      <p id="d1e712">In Fig. 1, two abscissa-axes are used to synchronize <inline-formula><mml:math id="M30" display="inline"><mml:mi>k</mml:mi></mml:math></inline-formula> – the radius number, <inline-formula><mml:math id="M31" display="inline"><mml:mi>i</mml:mi></mml:math></inline-formula> – revolution number of the workpiece; and <inline-formula><mml:math id="M32" display="inline"><mml:mi>z</mml:mi></mml:math></inline-formula> – stage number of the feed cycle
four radius in one stage of the cycle. The actual feed is shown on the graph
in the form of two nonlinear trajectories. The lower trajectory reflects the
change in the accumulated actual value <inline-formula><mml:math id="M33" display="inline"><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> during (<inline-formula><mml:math id="M34" display="inline"><mml:mrow><mml:mi>i</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula>) component revolutions, with <inline-formula><mml:math id="M35" display="inline"><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> being equal to the part of the allowance ground from the <inline-formula><mml:math id="M36" display="inline"><mml:mi>k</mml:mi></mml:math></inline-formula>th radius in <inline-formula><mml:math id="M37" display="inline"><mml:mrow><mml:mi>i</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula> revolutions of the part. At the same time, the lower trajectory shows the current value of the <inline-formula><mml:math id="M38" display="inline"><mml:mi>k</mml:mi></mml:math></inline-formula>th radius <inline-formula><mml:math id="M39" display="inline"><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> of the part that came from the previous turn.</p>
      <p id="d1e855">On the <inline-formula><mml:math id="M40" display="inline"><mml:mi>i</mml:mi></mml:math></inline-formula>th turn from the radius <inline-formula><mml:math id="M41" display="inline"><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> the layer of metal
<inline-formula><mml:math id="M42" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> , which is equal to the depth of cutting or the value of the actual feed per revolution of the workpiece, will be ground.
Therefore, the upper trajectory of the graph shows the change in the value
of the part of the allowance ground for <inline-formula><mml:math id="M43" display="inline"><mml:mi>i</mml:mi></mml:math></inline-formula> of the component revolutions or
accumulated in the actual feedrate per revolution. At the same time, the
upper trajectory shows the current change in the radius <inline-formula><mml:math id="M44" display="inline"><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>of the
part after the <inline-formula><mml:math id="M45" display="inline"><mml:mi>i</mml:mi></mml:math></inline-formula>th turn, i.e.

              <disp-formula specific-use="align" content-type="numbered"><mml:math id="M46" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E2"><mml:mtd><mml:mtext>2</mml:mtext></mml:mtd><mml:mtd><mml:mstyle displaystyle="true" class="stylechange"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:msub><mml:mi>Y</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E3"><mml:mtd><mml:mtext>3</mml:mtext></mml:mtd><mml:mtd><mml:mstyle displaystyle="true" class="stylechange"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>p</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>R</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E4"><mml:mtd><mml:mtext>4</mml:mtext></mml:mtd><mml:mtd><mml:mstyle class="stylechange" displaystyle="true"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>p</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mi>k</mml:mi></mml:munderover><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mi>i</mml:mi></mml:munderover><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mi>z</mml:mi></mml:munderover><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>p</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E5"><mml:mtd><mml:mtext>5</mml:mtext></mml:mtd><mml:mtd><mml:mstyle class="stylechange" displaystyle="true"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mi>k</mml:mi></mml:munderover><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mi>i</mml:mi></mml:munderover><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mi>z</mml:mi></mml:munderover><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E6"><mml:mtd><mml:mtext>6</mml:mtext></mml:mtd><mml:mtd><mml:mstyle displaystyle="true" class="stylechange"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:msub><mml:mi>y</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi mathvariant="italic">γ</mml:mi><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:msub><mml:mi>Y</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E7"><mml:mtd><mml:mtext>7</mml:mtext></mml:mtd><mml:mtd><mml:mstyle displaystyle="true" class="stylechange"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:msub><mml:mtext>w</mml:mtext><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E8"><mml:mtd><mml:mtext>8</mml:mtext></mml:mtd><mml:mtd><mml:mstyle displaystyle="true" class="stylechange"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>R</mml:mi><mml:mi>k</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>R</mml:mi><mml:mo>max⁡</mml:mo></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:msub><mml:mtext>w</mml:mtext><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>

          where <inline-formula><mml:math id="M47" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>R</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> – initial radial runout of workpiece, mm; <inline-formula><mml:math id="M48" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> – actual feed, mm/revolution; <inline-formula><mml:math id="M49" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mi>p</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> –
programmed feed, mm/double stroke; <inline-formula><mml:math id="M50" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> – radial component of
cutting force, kg; <inline-formula><mml:math id="M51" display="inline"><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> – current size of the radius in the
section of the machined surface, mm; <inline-formula><mml:math id="M52" display="inline"><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mo>max⁡</mml:mo></mml:msub></mml:mrow></mml:math></inline-formula> – maximum radius of the
workpiece, mm; <inline-formula><mml:math id="M53" display="inline"><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:msub><mml:mtext>w</mml:mtext><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> – initial radius of the workpiece on the <inline-formula><mml:math id="M54" display="inline"><mml:mi>k</mml:mi></mml:math></inline-formula>th radius, mm; <inline-formula><mml:math id="M55" display="inline"><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> – sum of the actual radial feeds accumulated at each billet radius for all the workpiece turns made, mm; <inline-formula><mml:math id="M56" display="inline"><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>p</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> – sum of
program feeds accumulated at each billet radius over all perfect double
strokes, mm; <inline-formula><mml:math id="M57" display="inline"><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> – elastic deformation of the technological system,
mm; <inline-formula><mml:math id="M58" display="inline"><mml:mi mathvariant="italic">γ</mml:mi></mml:math></inline-formula> – compliance of the technological system, mm kg<inline-formula><mml:math id="M59" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>.</p>
      <?pagebreak page334?><p id="d1e1621">Solving the set of the obtained Eqs. (1)–(8) with respect to the actual
feed-rate <inline-formula><mml:math id="M60" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:math></inline-formula> we obtain the following formula which makes
it possible to calculate current values of the actual feed-rate <inline-formula><mml:math id="M61" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> at each radius and the turn of the balance during the whole grinding cycle.
          <disp-formula id="Ch1.E9" content-type="numbered"><label>9</label><mml:math id="M62" display="block"><mml:mtable columnspacing="1em" class="split" rowspacing="0.2ex" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mo mathsize="2.5em">[</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mo>-</mml:mo><mml:mi mathvariant="italic">γ</mml:mi><mml:msub><mml:mi>K</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:mn mathvariant="normal">2</mml:mn><mml:mo>(</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>+</mml:mo><mml:msub><mml:mi>K</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mi mathvariant="italic">γ</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mfrac></mml:mstyle></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mfenced open="" close="]"><mml:mrow><mml:mo>+</mml:mo><mml:msqrt><mml:mrow><mml:msup><mml:mfenced open="(" close=")"><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mi mathvariant="italic">γ</mml:mi><mml:msub><mml:mi>K</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:mn mathvariant="normal">2</mml:mn><mml:mo>(</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>+</mml:mo><mml:msub><mml:mi>K</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mi mathvariant="italic">γ</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mfrac></mml:mstyle></mml:mfenced><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>p</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:msub><mml:mtext>w</mml:mtext><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mo>+</mml:mo><mml:msub><mml:mi>K</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mi mathvariant="italic">γ</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mrow></mml:msqrt></mml:mrow></mml:mfenced><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
        where <inline-formula><mml:math id="M63" display="inline"><mml:mrow><mml:msub><mml:mi>K</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M64" display="inline"><mml:mrow><mml:msub><mml:mi>K</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> – analytical coefficients characterizing the
interrelation of various technological parameters of the grinding process.
          <disp-formula id="Ch1.E10" content-type="numbered"><label>10</label><mml:math id="M65" display="block"><mml:mrow><mml:msub><mml:mi>K</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mi mathvariant="italic">π</mml:mi><mml:mtext>d</mml:mtext><mml:mi>B</mml:mi><mml:mi mathvariant="italic">σ</mml:mi><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo><mml:mn mathvariant="normal">9</mml:mn></mml:mrow><mml:mi>V</mml:mi></mml:mfrac></mml:mstyle><mml:mo>,</mml:mo><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:msub><mml:mi>K</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mi>B</mml:mi><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mi mathvariant="italic">σ</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:mfrac></mml:mstyle><mml:msqrt><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mi>D</mml:mi><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mi>n</mml:mi><mml:mtext>det</mml:mtext></mml:msub><mml:mo>(</mml:mo><mml:mi>D</mml:mi><mml:mo>+</mml:mo><mml:mi>d</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mfrac></mml:mstyle></mml:msqrt></mml:mrow></mml:math></disp-formula>
        Substituting the value of the actual feed per revolution into the given
above equations, calculated according to Eq. (9), we can calculate the
current values of the cutting force, elastic deformation and the radius of
the part on each turn during the whole grinding cycle.</p>
      <p id="d1e1956">Given that the feed per turn is functionally (Eq. 11) related to the minute feed
          <disp-formula id="Ch1.E11" content-type="numbered"><label>11</label><mml:math id="M66" display="block"><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mrow><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>f</mml:mi><mml:mo>,</mml:mo><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:msub><mml:mo>⋅</mml:mo><mml:msub><mml:mi>n</mml:mi><mml:mtext>det </mml:mtext></mml:msub></mml:mrow></mml:math></disp-formula>
        The graphs of the minute feed for each radius and revolution of the part
during the entire grinding cycle, synchronized in time, is shown in Fig. 2.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F2" specific-use="star"><?xmltex \currentcnt{2}?><label>Figure 2</label><caption><p id="d1e2014">The graphs of minute feeds (programmed and actual) in a
two-step grinding cycle. Points of actual feed rate are calculated by four
radiuses on every round of workpiece</p></caption>
        <?xmltex \igopts{width=369.885827pt}?><graphic xlink:href="https://ms.copernicus.org/articles/10/331/2019/ms-10-331-2019-f02.png"/>

      </fig>

      <p id="d1e2023">Then, after substituting Eq. (9) into Eq. (11), it is possible to calculate the
current minute feed values for each radius and revolution of the part during
the entire grinding cycle. The calculation of the grinding cycle time <inline-formula><mml:math id="M67" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mtext>c</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula> is carried out by summing the time of all the revolutions of the
part per grinding cycle using Eq. (12)
          <disp-formula id="Ch1.E12" content-type="numbered"><label>12</label><mml:math id="M68" display="block"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mtext>c</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mi>z</mml:mi></mml:munderover><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mi>i</mml:mi></mml:munderover><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        where <inline-formula><mml:math id="M69" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> – time of one turn of the workpiece.</p>
      <p id="d1e2096">Thus, the combination of Eqs. (1)–(12) is a wide-range analytical model
of the allowance removal process for the operation of circular external
grinding with a radial feed that enables the calculation of the current
values of all types of feeds, forces, elastic deformations and radii of the
surface for any given grinding cycle and processing conditions.</p>
</sec>
<sec id="Ch1.S3">
  <label>3</label><title>Methodology for calculating optimal feed rates for specified grinding conditions</title>
      <p id="d1e2107">Limitations of models are specifically developed by Pereverzev and Akintseva (2015, 2017), Akintseva and Pereverzev (2017), Bertsekas (1995), Lee (2008) to support objective function of
accuracy, roughness and hardness of the treated surface as well as
technological limitations on the technical data of the machine and the
strength of the grinding wheel.</p>
      <p id="d1e2110">To calculate the optimal step-by-step grinding cycle, it is necessary to
assign optimum minute feed rates for each stage of the cycle and assign the
optimum values of the allowance parts taken at each stage of the cycle.
According to the active monitoring device, the CNC system switches the cycle
stages during grinding. The minimum grinding cycle time is assumed as the
target function.</p>
      <?pagebreak page335?><p id="d1e2113">With optimum feed and allowance values at each stage of the cycle, it is
necessary to ensure that the cycle will be maximum in speed from all
possible cycles and provide a given quality of the treated surface. For this
purpose, we developed models of the objective function restrictions for
accuracy, non-brushing, roughness of the machined surface, as well as
technological limitations on the technical data of the machine and the
crumbling of the grinding wheel. These models described by Pereverzev and Akintseva (2015, 2017), and Akintseva and Pereverzev (2017).<?xmltex \hack{\newpage}?></p>
      <p id="d1e2117">This formulation of the problem of calculating the parameters of the
grinding cycle is of an optimization nature and requires the choice of a
mathematical optimization method.</p>
      <p id="d1e2121">The analysis of suitable methods described by Pereverzev and Akintseva (2015, 2017), and Akintseva and Pereverzev (2017)
showed that the most suitable for optimizing cycles is the dynamic
programming method (DPM) developed by Bellman (1960). The optimization of
the DPM cycle is carried out by analogy with the optimization of the
transport task, in which a network of roads with intermediate stations is
specified, it is required to find the optimal transport route between two
points A and B. With regard to optimizing the processing cycle, the road
diagram is similar to the “grinding process state options” for different
feeds when different parts of the allowance are being removed. In our case,
the starting point is the initial state of the process (procurement
parameters, etc.). The final point is the final state of the process (the
parameters of the finished part in accuracy and quality, the main time,
etc.). As a result, the problem arises of finding a cycle (the “optimal
trajectory of removing the allowance”), which has a minimal basic time among
various possible variants of cycles, provided all the technological
restrictions are met.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F3" specific-use="star"><?xmltex \currentcnt{3}?><label>Figure 3</label><caption><p id="d1e2126">Grid for optimizing the grinding cycle (<inline-formula><mml:math id="M70" display="inline"><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mtext>pm</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula> – nominal program minute feed in the DPM optimization system, mm min<inline-formula><mml:math id="M71" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>).</p></caption>
        <?xmltex \igopts{width=369.885827pt}?><graphic xlink:href="https://ms.copernicus.org/articles/10/331/2019/ms-10-331-2019-f03.png"/>

      </fig>

      <p id="d1e2158">When optimizing the grinding cycle with the use of DPM instead of the road
network and stations, a coordinate grid with the axes “Program feed
<inline-formula><mml:math id="M72" display="inline"><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mtext>pm</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula>” – “Allowance <inline-formula><mml:math id="M73" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mtext>n</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula>” (Fig. 3) is used. Since the DPM is a discrete
optimization method, the feed and overflow axes are divided into small
discrete parts. Each intersection of the coordinate grid has an information
cell containing the data necessary for the calculation (including the
optimal time to reach the state and the coordinates from which the optimum
move is made).</p>
      <p id="d1e2183">Each move is checked for validity of accuracy, roughness and hardness of the
treated surface as well as technological limitations on the technical data
of the machine and the strength of the grinding wheel. The feed level number
<inline-formula><mml:math id="M74" display="inline"><mml:mi>m</mml:mi></mml:math></inline-formula> varies in the feed range of the machine from <inline-formula><mml:math id="M75" display="inline"><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mo>≤</mml:mo><mml:mi>m</mml:mi><mml:mo>≤</mml:mo><mml:mi>M</mml:mi></mml:mrow></mml:math></inline-formula>, where <inline-formula><mml:math id="M76" display="inline"><mml:mi>M</mml:mi></mml:math></inline-formula> is the maximum feed level (number). The number of the allowance part <inline-formula><mml:math id="M77" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula> varies from <inline-formula><mml:math id="M78" display="inline"><mml:mrow><mml:mn mathvariant="normal">0</mml:mn><mml:mo>≤</mml:mo><mml:mi>n</mml:mi><mml:mo>≤</mml:mo><mml:mi>N</mml:mi></mml:mrow></mml:math></inline-formula>, where <inline-formula><mml:math id="M79" display="inline"><mml:mi>N</mml:mi></mml:math></inline-formula> is the maximum number of the remaining part of the allowance. Moreover, for <inline-formula><mml:math id="M80" display="inline"><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula> the allowance is maximum, and for <inline-formula><mml:math id="M81" display="inline"><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mi>N</mml:mi></mml:mrow></mml:math></inline-formula> the allowance is zero, i.e. allowance was removed completely. The level
<inline-formula><mml:math id="M82" display="inline"><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mi>N</mml:mi></mml:mrow></mml:math></inline-formula> is the reference point of the removed allowance from the center of the
tolerance field of the radius of the part to the maximum value of the radius
of the workpiece.</p>
      <p id="d1e2283">In contrast to the classical transport problem, in which advance time (i.e.,
time of travel between stations) and permissible strokes is known, then in
the optimization of grinding cycles these key values are not known in
advance. Therefore, the determination of the travel time (i.e., the removal
time of the allowance discrepancy between adjacent grid nodes) is carried
out by modulating the removal of the metal, equal to the magnitude of the
allowance discrepancy between the coordinates of the nodes of the
optimization grid. Simulation is carried out using the removal model of the
allowance (1)–(12). In the process of modeling the removal of the part of
the allowance, the parameters of the state of the grinding process (yield,
cutting depth, cutting force, travel time, surface radius, elastic
deformation) are determined. Each stroke is checked for tolerance by four
main limitations of the operation: accuracy, roughness, roughness and
incidence of the wheel, ensuring the required quality and precision of
processing, as well as the main limitations of the technological system
(machine drive power and feed ranges).</p>
      <p id="d1e2287">Admissibility of each move is checked step by step on the principle whether
after the next move it is possible to obtain a suitable part with minimum
cutting conditions due to the technical capabilities of the machine. If not,
the move is considered to be unacceptable.</p>
      <p id="d1e2290">The method of dynamic programming allows you to organize a direct search of
options and avoid a full search. Each<?pagebreak page336?> station is characterized by two state
coordinates, which will be written in the form [<inline-formula><mml:math id="M83" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula>, <inline-formula><mml:math id="M84" display="inline"><mml:mi>m</mml:mi></mml:math></inline-formula>].</p>
      <p id="d1e2307">The process of optimizing the grinding cycle is carried out in full
compliance with the DPM. The search of the moves (Fig. 3) is carried out by
calculating the time of all admissible moves to the state [<inline-formula><mml:math id="M85" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula>, <inline-formula><mml:math id="M86" display="inline"><mml:mi>m</mml:mi></mml:math></inline-formula>] from the
previous level [<inline-formula><mml:math id="M87" display="inline"><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M88" display="inline"><mml:mi>l</mml:mi></mml:math></inline-formula>], where <inline-formula><mml:math id="M89" display="inline"><mml:mi>l</mml:mi></mml:math></inline-formula> – feedrate discrepancy number at <inline-formula><mml:math id="M90" display="inline"><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula> level. In
accordance with the MDP procedure, from all competing moves, the optimal
move with the minimum optimal time <inline-formula><mml:math id="M91" display="inline"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> was calculated using Eq. (13):
          <disp-formula id="Ch1.E13" content-type="numbered"><label>13</label><mml:math id="M92" display="block"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mi>m</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:munder><mml:mo movablelimits="false">min⁡</mml:mo><mml:mi>l</mml:mi></mml:munder><mml:mfenced close="}" open="{"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo><mml:mi>l</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msubsup><mml:mi mathvariant="italic">τ</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo><mml:mi>l</mml:mi></mml:mrow></mml:msubsup></mml:mrow></mml:mfenced><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        where <inline-formula><mml:math id="M93" display="inline"><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">τ</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo><mml:mi>l</mml:mi></mml:mrow></mml:msubsup></mml:mrow></mml:math></inline-formula> – time of travel to the state with coordinates
[<inline-formula><mml:math id="M94" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula>, <inline-formula><mml:math id="M95" display="inline"><mml:mi>m</mml:mi></mml:math></inline-formula>] from the previous level of allowed states [<inline-formula><mml:math id="M96" display="inline"><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M97" display="inline"><mml:mi>l</mml:mi></mml:math></inline-formula>] at all values l from 1 till <inline-formula><mml:math id="M98" display="inline"><mml:mi>M</mml:mi></mml:math></inline-formula>. on the previous <inline-formula><mml:math id="M99" display="inline"><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula> level, min; <inline-formula><mml:math id="M100" display="inline"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mi>l</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> – optimal time to reach the
station at the previous level, min; <inline-formula><mml:math id="M101" display="inline"><mml:mi>m</mml:mi></mml:math></inline-formula> – number of programmed supply
discrepancies in the DPM optimization system; <inline-formula><mml:math id="M102" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula> – number of the allowance
discrepancy in the DPM optimization system. It should be noted that the
travel time <inline-formula><mml:math id="M103" display="inline"><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">τ</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo><mml:mi>l</mml:mi></mml:mrow></mml:msubsup></mml:mrow></mml:math></inline-formula> is the time necessary to remove the part
of the allowance between the <inline-formula><mml:math id="M104" display="inline"><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M105" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula> levels when the program feed is turned on at level <inline-formula><mml:math id="M106" display="inline"><mml:mi>m</mml:mi></mml:math></inline-formula> from the process state [<inline-formula><mml:math id="M107" display="inline"><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M108" display="inline"><mml:mi>l</mml:mi></mml:math></inline-formula>]. The travel time value <inline-formula><mml:math id="M109" display="inline"><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">τ</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo><mml:mi>l</mml:mi></mml:mrow></mml:msubsup></mml:mrow></mml:math></inline-formula> is determined by calculation, using the allowance removal model (1)–(12).</p>
      <p id="d1e2654">After determining the optimal move, the information cell of state [<inline-formula><mml:math id="M110" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula>, <inline-formula><mml:math id="M111" display="inline"><mml:mi>m</mml:mi></mml:math></inline-formula>]
records the level number from which the optimal move is made. The recording
is made in the variable <inline-formula><mml:math id="M112" display="inline"><mml:mrow><mml:msup><mml:mi>m</mml:mi><mml:mo>*</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula> and is recorded in the information box together
with the value of the optimal time and technological parameters of the
grinding process.</p>
      <p id="d1e2682">A directional search of variants from the calculation of the optimal time
values is carried out at all allowance discrepancies.</p>
      <p id="d1e2685">Thus, the optimization process begins with the execution of the set of
straight strokes on the first allowance discrepancy. Then a transition is
made from the previous (<inline-formula><mml:math id="M113" display="inline"><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula>)th level of allowance, at the <inline-formula><mml:math id="M114" display="inline"><mml:mi>m</mml:mi></mml:math></inline-formula>th program feed rate to the <inline-formula><mml:math id="M115" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula>th level of allowance. Each turn is checked for the
admissibility according to the given restrictions of the objective function.
When the [<inline-formula><mml:math id="M116" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula>, <inline-formula><mml:math id="M117" display="inline"><mml:mi>m</mml:mi></mml:math></inline-formula>]th state of the grinding process is reached by one or more
permissible moves, the choice of the optimal stroke ensuring the minimum
time <inline-formula><mml:math id="M118" display="inline"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> to reach the state [<inline-formula><mml:math id="M119" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula>, <inline-formula><mml:math id="M120" display="inline"><mml:mi>m</mml:mi></mml:math></inline-formula>] should be made. The coordinates
of the previous state of the process in the variable <inline-formula><mml:math id="M121" display="inline"><mml:mrow><mml:msup><mml:mi>m</mml:mi><mml:mo>*</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula> are remembered, the transition is considered optimal. After reaching the last level of allowance <inline-formula><mml:math id="M122" display="inline"><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mi>N</mml:mi></mml:mrow></mml:math></inline-formula> (and hence the final state of the process), at the last level of the grid allowance, there may be several cells located at different numbers <inline-formula><mml:math id="M123" display="inline"><mml:mi>m</mml:mi></mml:math></inline-formula> of program feed rates and having different time <inline-formula><mml:math id="M124" display="inline"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>. This means that it is
possible to obtain a suitable workpiece surface for different values of the
actual feed rate at the end of the cycle and the feed rates corresponding to
them, but for different processing times. Therefore, out of all available
options at the last level of the allowance, we select one having the minimum
time <inline-formula><mml:math id="M125" display="inline"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>N</mml:mi><mml:mo>,</mml:mo><mml:mi>M</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> of reaching the final state, which is the minimum time of the
<inline-formula><mml:math id="M126" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mtext>c</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula> cycle, i.e.
          <disp-formula id="Ch1.E14" content-type="numbered"><label>14</label><mml:math id="M127" display="block"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mtext>c</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:munder><mml:mo movablelimits="false">min⁡</mml:mo><mml:mi>m</mml:mi></mml:munder><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>N</mml:mi><mml:mo>,</mml:mo><mml:mi>m</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        where <inline-formula><mml:math id="M128" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mtext>c</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula> – grinding cycle time, min;</p>
      <p id="d1e2875">To fix the optimum control path for the grinding cycle, it is necessary to
perform a retraction procedure that starts from this single final state
having <inline-formula><mml:math id="M129" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mtext>c</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula>. By analogy with the<?pagebreak page337?> procedure for the back stroke in
the optimization of the route, we determine the previous state of the
process <inline-formula><mml:math id="M130" display="inline"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>, using the stored coordinates of the previous state of
the process <inline-formula><mml:math id="M131" display="inline"><mml:mrow><mml:msup><mml:mi>m</mml:mi><mml:mo>*</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula>, going from the final level of the allowance to the initial,
memorizing the program speeds of the feeds and allowances, after which the
program speeds of the feeds were changed by the commands of the active
control device.</p>
      <p id="d1e2920">Figure 2 shows the procedure for the back stroke. In cells on the
intersections of the coordinate grid, only numerical values of <inline-formula><mml:math id="M132" display="inline"><mml:mrow><mml:msup><mml:mi>m</mml:mi><mml:mo>*</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula> are
conditionally given, which is obtained after optimization of the grinding
cycle. Empty cells indicate that these intersections have no admissible
moves.</p>
      <p id="d1e2934">The reverse move is made on the basis of the use of the reverse mark <inline-formula><mml:math id="M133" display="inline"><mml:mrow><mml:msup><mml:mi>m</mml:mi><mml:mo>*</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula>, which
stores the feedrate number, from which the optimal move is made. Using this
data, the optimal cycle path is restored, which ensures the minimum cycle
time. Therefore, for demonstrating the procedure for the return stroke Fig. 3 presents arrows which in accordance with the value of <inline-formula><mml:math id="M134" display="inline"><mml:mrow><mml:msup><mml:mi>m</mml:mi><mml:mo>*</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula> show the coordinate grid node at the previous level, from which the optimal move is made.</p>
      <p id="d1e2959">Thus, we have developed a method for designing optimal round grinding cycles
and makes it possible to calculate optimal radial feed rates at all stages
of the cycle and optimal distribution of the removed allowance in the
control cycle steps which need the minimum main cycle time of grinding.</p>
</sec>
<sec id="Ch1.S4">
  <label>4</label><title>Synthesis of the digital twin technology and dynamic programming for testing the grinding cycle for resistance to variable technological factors</title>
      <p id="d1e2970">The parameters of the grinding cycle on the CNC machine are constant for any
variable batch processing conditions: blunting the grinding wheel grains,
allowance swing, oscillation of the initial radial runout of the workpiece,
size wear of the grinding wheel and reduction of the cutting speed, variable
hardness of the technological system in different sections of the surface to
be treated, the vibration of the physical and mechanical properties of the
material of the part, the variation in the parameters of the characteristics
of the grinding wheel, the operation error of the active control device. The
number of combinations of variable factors, even for three values of each,
is measured in hundreds. But the grinding cycle remains unchanged with any
aggregate combination of variable factors, which causes a fluctuation in the
accuracy of processing and other quality parameters of the grinding surface.</p>
      <p id="d1e2973">The optimum grinding cycle, designed on the basis of the average
deterministic grinding conditions, has low reliability and stability to the
influence of variable technological factors on the accuracy of processing
and the provision of other quality parameters. Therefore, the
above-mentioned procedure for designing the optimum grinding cycle does not
take into account unforeseen situations associated with unstable grinding
conditions, blunt grinding of wheel grains, the variation in the allowance
or the initial radial runout of the workpiece, etc., which can lead to the
appearance of defects in grinding.</p>
      <p id="d1e2976">To design an optimal grinding cycle that is resistant to the combined effect
of various variable technological factors, a system for diagnosing the
grinding cycle for stability and the optimal cycle design system is
necessary.</p>
      <p id="d1e2979">When developing a diagnostic system for the stability of the grinding cycle,
the concept of a “digital twin” (DT) is applied to the cumulative effect of
the constantly changing variables that arise when processing a batch of
parts. Within the framework of this concept, a system model DT was designed
for: (1) Preventing the defect and the causes of its occurrence in the
operation of round grinding; (2) Increasing the reliability and stability of
grinding cycle to the combined effect of variable factors; (3) Providing full
automation of the design of the control program at the stage of calculation
and programming of cutting parameters; (4) Forecasting fluctuations in
accuracy, roughness, hardness of the machined surface after machining a
batch of parts; (5) Fixing an array of variable factors which caused defects;
(6) Forming the array of restrictions of the objective function with variable
grinding conditions for the cycle optimization system; (7) Automatic transfer
of the array of the objective function restrictions according to the
changing grinding conditions in the cycle optimization system; (8) Automatic
loading of the tested cycle into the diagnostic system; (9) Automatic
completion of grinding cycle diagnostics.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F4" specific-use="star"><?xmltex \currentcnt{4}?><label>Figure 4</label><caption><p id="d1e2985">Scheme of interaction of the system for optimizing the MDP
cycle with the digital twin of the cycle testing system: (1) defect
prevention; (2) improving the reliability of the grinding cycle; (3) forecasting the accuracy and quality of processing; (4) formation of
restrictions for variable conditions; (5) generation of variable grinding
conditions; (6) cycle testing.</p></caption>
        <?xmltex \igopts{width=341.433071pt}?><graphic xlink:href="https://ms.copernicus.org/articles/10/331/2019/ms-10-331-2019-f04.png"/>

      </fig>

      <p id="d1e2994">The developed system “DT-CicleStab” provides synthesis of diagnostic system
and optimization system according to the following algorithm (Fig. 4): (1) The first version of the optimum grinding cycle is designed based on
averaged grinding conditions and checking restrictions under these
conditions; (2) The optimum grinding cycle is transferred to the cycle
diagnostics system; (3) An array of combinations of variable factors is
generated; (4) diagnosis of the optimum grinding cycle is carried out for all variants
given by an array of variable factors combinations. In the process of
diagnostics, the optimal cycle is tested by modeling the grinding process
during the whole processing cycle in hundreds of combination variants of
variable technological factors (formed in point 3 of the algorithm)which are
variable in the given ranges of their variation. After each testing, the
estimated values of machining accuracy, roughness, surface hardness and
other specified quality parameters are evaluated; (5) All types of defects
and the conditions for their occurrence are identified. If there is no
defect, this algorithm is implemented; (6) Forming an array of additional
restrictions of the objective function on the basis of the revealed
conditions for the appearance of defects; (7) Transfer of additional
restrictions of the objective function to the system for designing the
optimum grinding cycle; (8) Optimization of the grinding cycle with an
additional array of the objective<?pagebreak page338?> function restrictions. A transition is
made to point 2 of this algorithm; (9) Forming a report and documentation on
the design of the optimum grinding cycle; (10) Transfer of data for the
optimum grinding cycle to the control system design for the CNC machine
tool. Thus, we have developed a system for designing optimal cycles with
improved reliability and resistance to variable grinding conditions.</p>
</sec>
<sec id="Ch1.S5">
  <label>5</label><title>Experimental studies</title>
      <p id="d1e3006">This article presents experimental studies results to proving the validity
and accuracy of the predicted model of cutting forces (1) and removal
allowance model (10) in the grinding cycle. Other models have been
experimentally proven in the works of Malkin (1981), Rowe and Ebbrell (2004),
Petrakov and Chamata (2015), Pereverzev and Akintseva (2017), Lur'e (1979). Several
aspects of the models used for this purpose has been developed, based on the
prediction of the performance of grinding cycles.</p>
      <p id="d1e3009">The radial component of the cutting force was adopted as the main object for
testing the adequacy of grinding models, because this parameter is the most
sensitive one to the action of all technological factors.</p>
      <p id="d1e3012">Since the grinding process is accompanied by a constant change in the
blunting of the grains of the wheel, and the process is non-deterministic in
connection with the continuous processes of blunting the grains and tearing
them out of the bundle, then the cutting force is also constantly changing.
Therefore, in carrying out the experiments, the theoretical value of the
cutting force was calculated for the worst processing conditions,
corresponding to grinding by a blunt wheel. Thus, the procedure for
performing the experiments is as follows: (1) A circular mortise grinding of
the sample is carried out according to a given grinding cycle; (2) During the
grinding process, the numerical values of the radial component of the
cutting force, the actual feed rate and the cycle time are measured and
recorded; (3) For a given cycle, a theoretical calculation is made of the
variation in the radial component of the cutting force, the actual feed rate
and the cycle time, depending on the program feed rates taken in the
experiment; (4) The estimation of coincidence of theoretical and experimental
values of controlled parameters is made; (5) Experiments are conducted with
various combinations of parameters of the grinding wheel characteristics and
the processed steel grades; (6) The processing cycle is assigned from the
working condition of the grinding wheel in the blunt mode; (7) The allowance
for processing is assigned sufficiently large to show different stages of
blunting of the grinding wheel grains.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F5" specific-use="star"><?xmltex \currentcnt{5}?><label>Figure 5</label><caption><p id="d1e3018">Functional diagram of the computer Center.</p></caption>
        <?xmltex \igopts{width=398.338583pt}?><graphic xlink:href="https://ms.copernicus.org/articles/10/331/2019/ms-10-331-2019-f05.png"/>

      </fig>

      <?xmltex \floatpos{t}?><fig id="Ch1.F6" specific-use="star"><?xmltex \currentcnt{6}?><label>Figure 6</label><caption><p id="d1e3029">Graphs of complex experiment No. 1: (1) experimental
graph of the change in the program feed rate in time; (2) experimental
graph of the change in the actual speed of delivery in time; (3) experimental graph of the variation of the radial cutting force in time; (4) theoretical graph of the change in the radial component of the cutting
force in time.</p></caption>
        <?xmltex \igopts{width=341.433071pt}?><graphic xlink:href="https://ms.copernicus.org/articles/10/331/2019/ms-10-331-2019-f06.png"/>

      </fig>

      <p id="d1e3038">The experiments were carried out on a measuring and computing stand (MCS)
(Fig. 5) created on the basis of a CNC machine of model 3M151F2. The bench
is designed for measuring and recording various parameters of the grinding
process, entering measurement results in the computer, primary processing,
sorting, storage and graphical display of information, automatic adjustment
of number of process parameters, program control of processing modes. To
perform these functions, the MCS is equipped with number of sensors which
together with standard instrumentation are connected to the common computer
bus through the modules of the standard CAMAC interface (Fig. 5).</p>
      <p id="d1e3041">The stand provides two options for measuring the radial cutting force. In
the first variant, the measurement is performed using displacement sensors
(inductive transducers) mounted in the flange of the bearings of the
grinding head spindle. The springs of inductive converters constantly press
the levers against the spindle.</p>
      <p id="d1e3044">Grinding of the samples was carried out with a wheel of <inline-formula><mml:math id="M135" display="inline"><mml:mrow><mml:mn mathvariant="normal">600</mml:mn><mml:mo>×</mml:mo><mml:mn mathvariant="normal">40</mml:mn><mml:mo>×</mml:mo></mml:mrow></mml:math></inline-formula> at the frequency of the wheel rotation and the<?pagebreak page339?> sample,
respectively, 1480 and 180 rpm. 1 % soda solution was used as a
technological liquid. The compliance of the technological system was not
changed and according to the measurement results it was 0.0002 mm N<inline-formula><mml:math id="M136" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>. The
experimental values of the monitored parameters were displayed graphically
on the display screen and on the printing device.</p>
      <p id="d1e3074">Figure 6 shows the results of two experiments in the form of a set of the
following graphs visually monitored on the display screen: (1) Experimental
graph for changing the program feedrate in time (shown as a step graph); (2) Experimental graph of the change in the actual feed rate in time; (3) Experimental graph of the variation of the radial component of the cutting
force in time; (4) Theoretical graph of the change in the value of the radial
component of the cutting force. The theoretical calculation of the cycle
time differs from the experimental one by 0.1 % … 0.5 %.</p>
      <p id="d1e3077">The theoretical models of the metal removal process during the grinding
cycle are adequate to the results of the<?pagebreak page340?> experiments, which allowed us to
use them in calculating the regulatory and reference database for the design
of grinding cycles and the creation of automated cycle design systems.</p>
</sec>
<sec id="Ch1.S6" sec-type="conclusions">
  <label>6</label><title>Conclusions</title>
      <p id="d1e3088"><list list-type="order">
          <list-item>

      <p id="d1e3093">The developed methodology for the design of automatic grinding cycles with
optimal parameters of cutting conditions for CNC machines allows to ensure
the stable quality of the surface to be processed in terms of accuracy,
roughness and hardness, under varying conditions for processing a batch of
parts.</p>
          </list-item>
          <list-item>

      <p id="d1e3099">The methodology of designing optimal cycles is based on the synthesis of:
<list list-type="bullet"><list-item>
      <p id="d1e3104">wide-range analytical model of cutting forces;</p></list-item><list-item>
      <p id="d1e3108">removal allowance model based on the relationship of the elastic
deformations of the technological system with the parameters of cutting
conditions and processing conditions;</p></list-item><list-item>
      <p id="d1e3112">Error of processing model can be calculated for given grinding conditions;</p></list-item><list-item>
      <p id="d1e3116">optimization technique of grinding cycle based on dynamic programming
method;</p></list-item><list-item>
      <p id="d1e3120">stability of Diagnostic and reliability of the grinding cycle to varying
conditions;</p></list-item></list></p>
          </list-item>
          <list-item>

      <p id="d1e3126">results and practical optimization also involves planning to increasing the
level of Industrial Automation programming for CNC machines in digital
production.</p>
          </list-item>
          <list-item>

      <p id="d1e3132">The developed model of forming is analytical, due to the fact that it was
developed based on mathematical interrelation of cutting force, elastic
deformations of a technological system, cutting depth, programmed and actual
feeds and treated surface radiuses sizes.</p>
          </list-item>
          <list-item>

      <p id="d1e3138">The obtained model of allowance removal is a wide-range in terms of
variability of technological factors. The model has an analytical character
due to the fact that it was obtained on the basis of mathematical
interrelation between a radius size of a treated surface and actual feeds
with analytical model of cutting force, elaborated on the basis of
fundamental mechanics conformity of cutting process and the theory of
plastic deformation of metal in a cutting area and which makes a connection
with major technological factors, which are:
<list list-type="bullet"><list-item>
      <p id="d1e3143">physical and mechanical properties of a grinded metal (stress intensity);</p></list-item><list-item>
      <p id="d1e3147">geometrical parameters of the contact area of a wheel and a workpiece
(actual speed of feed, diameters of a workpiece and a wheel, rotational
speed of a workpiece, width of a workpiece treated area);</p></list-item><list-item>
      <p id="d1e3151">elastic properties of the technological system (elastic deformation and
compliance of the technological system);</p></list-item><list-item>
      <p id="d1e3155">characteristic of a grinding wheel and the grade of wheel grains blunting.</p></list-item></list></p>
          </list-item>
          <list-item>

      <p id="d1e3162">Interrelation between the model of allowance removal with cutting force
gives an opportunity to estimate the influence of different technological
factors on productivity, accuracy and quality of processing during variation
of the control parameters of the cylindrical exterior grinding cycle. These
parameters are cycle stage quantity, program speeds of radial speed on each
cycle stage, distribution of allowance on cycle stages.</p>
          </list-item>
          <list-item>

      <p id="d1e3168">The grinding cycle diagnostics system allows you to design an optimal
grinding cycle, ensuring a stable quality of processing a batch of parts.</p>
          </list-item>
        </list></p>
</sec>

      
      </body>
    <back><notes notes-type="dataavailability"><title>Data availability</title>

      <p id="d1e3177">All data included in this study are available upon request by contacting the corresponding author.</p>
  </notes><notes notes-type="authorcontribution"><title>Author contributions</title>

      <p id="d1e3183">PP has made substantial contributions to the conception and design of the work. AV has made the acquisition, analysis and the interpretation of data for the work. DV has drafted the work or revised it critically for important intellectual content. MA has made a contribution to the acquisition of simulation experimental data and data collation.</p>
  </notes><notes notes-type="competinginterests"><title>Competing interests</title>

      <p id="d1e3189">The authors declare that they have no conflict of interest.</p>
  </notes><ack><title>Acknowledgements</title><p id="d1e3195">South Ural State University is grateful for the financial support of the
Ministry of Education and Science of the Russian Federation.</p></ack><notes notes-type="financialsupport"><title>Financial support</title>

      <p id="d1e3200">This research has been supported by the Ministry of Education and Science of the Russian Federation (grant no. 9.5589.2017/8.9).</p>
  </notes><notes notes-type="reviewstatement"><title>Review statement</title>

      <p id="d1e3206">This paper was edited by Xichun Luo and reviewed by two anonymous referees.</p>
  </notes><ref-list>
    <title>References</title>

      <ref id="bib1.bib1"><label>1</label><mixed-citation>Akintseva, A. V. and Pereverzev, P. P.: Prospects for the development
of the theory of designing optimal cycles of machining in a multidimensional
space of control parameters, EDP Sciences, 129, 01018, <ext-link xlink:href="https://doi.org/10.1051/matecconf/201712901018" ext-link-type="DOI">10.1051/matecconf/201712901018</ext-link>, 2017.</mixed-citation></ref>
      <ref id="bib1.bib2"><label>2</label><mixed-citation>Alagumurthi, N., Panairadja, K., and Soundararajan, V.: Optimization of
grinding process through design of experiment (DOE) – a comparative study,
Mater. Manuf. Proc., 21, 19–21, <ext-link xlink:href="https://doi.org/10.1080/AMP-200060605" ext-link-type="DOI">10.1080/AMP-200060605</ext-link>, 2007.</mixed-citation></ref>
      <ref id="bib1.bib3"><label>3</label><mixed-citation>
Amitay, G., Malkin, S., and Koren, Y.: Adaptive control optimization of
grinding, J. Eng. Ind., 103, 103–108, 10.31399/asm.hb.v16.a0002175, 1989.</mixed-citation></ref>
      <ref id="bib1.bib4"><label>4</label><mixed-citation>
Bellman, R.: Dynamic programming, Foreign Literature Publishing House,
Moscow, Russia, 1960.</mixed-citation></ref>
      <ref id="bib1.bib5"><label>5</label><mixed-citation>
Bertsekas, D.: Dynamic programming and optimal control, Prentice Hall, Reading, New
Jersey, USA, 1995.</mixed-citation></ref>
      <ref id="bib1.bib6"><label>6</label><mixed-citation>Cahill, M. J., Bechtold, M. J., Fess, E., Wolfs, F. L., and Bechtold,
R.: Ultrasonic precision optical grinding technology, P. SPIE, 96330, <ext-link xlink:href="https://doi.org/10.1117/12.2195977" ext-link-type="DOI">10.1117/12.2195977</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bib7"><label>7</label><mixed-citation>Fang, F., Cheng, K., Ding, H., Chen, S., and Zhao, L.: Sustainable
design and analysis of CNC machine tools: sustainable design index based
approach and its application perspectives, ASME Proceedings, 3, MSEC2016-8730, <ext-link xlink:href="https://doi.org/10.1115/MSEC2016-8730" ext-link-type="DOI">10.1115/MSEC2016-8730</ext-link>, 2016.</mixed-citation></ref>
      <ref id="bib1.bib8"><label>8</label><mixed-citation>Hashimoto, F. and Lahoti, G. D.: Optimization of set-up conditions for
stability of the centerless grinding process, CIRP Ann.-Manuf. Techn.,
53, 1271–274, <ext-link xlink:href="https://doi.org/10.1016/S0007-8506(07)60696-9" ext-link-type="DOI">10.1016/S0007-8506(07)60696-9</ext-link>,
2004.</mixed-citation></ref>
      <ref id="bib1.bib9"><label>9</label><mixed-citation>
Korchak, S. N.: The Productivity of the grinding of steel parts,
Mashinostroenie, Moscow, Russia, 1974.</mixed-citation></ref>
      <ref id="bib1.bib10"><label>10</label><mixed-citation>Krajnik, P., Drazumeric, R., Meyer, D., Kopac, J., and Zeppenfeld, C.:
Simulation of workpiece forming and centre displacement in plunge
centreless grinding, Int. J. Mach. Tool. Manu., 48, 824–831, <ext-link xlink:href="https://doi.org/10.1016/j.ijmachtools.2007.12.008" ext-link-type="DOI">10.1016/j.ijmachtools.2007.12.008</ext-link>, 2008.</mixed-citation></ref>
      <ref id="bib1.bib11"><label>11</label><mixed-citation>Lee, C. W.: Dynamic optimization of the grinding process in batch
production, J. Manuf. Sci. E.-T. ASME, 131, 485–494, <ext-link xlink:href="https://doi.org/10.1115/msec_icmp2008-72212" ext-link-type="DOI">10.1115/msec_icmp2008-72212</ext-link>, 2008.</mixed-citation></ref>
      <ref id="bib1.bib12"><label>12</label><mixed-citation>Lee, T. S., Ting, T. P., Lin, Y. J., and Htay, T.: A particle swarm
approach for grinding process optimization, J. Adv. Manuf.
Technol., 33, 1128–1135, <ext-link xlink:href="https://doi.org/10.1007/s00170-006-0538-y" ext-link-type="DOI">10.1007/s00170-006-0538-y</ext-link>, 2006.</mixed-citation></ref>
      <ref id="bib1.bib13"><label>13</label><mixed-citation>
Lin, L., Yong, L., Liwei, W., and Jun, X.: PMAC-based tracking control
system for 8-axis automated tape-laying machine, Chinese J. Aeronaut., 22, 558–563, 10.1016/S1000-9361(08)60141-7, 2009.</mixed-citation></ref>
      <ref id="bib1.bib14"><label>14</label><mixed-citation>
Liu, Y. M., Yang, T. Y., He, Z., and Li, J. Y.: Analytical modelling of
grinding process in rail profile correction cnsidering grinding pattern,
Arch Civ. Mech. Eng., 18, 669–678, 10.1016/j.acme.2017.10.009, 2017.</mixed-citation></ref>
      <ref id="bib1.bib15"><label>15</label><mixed-citation>
Lur'e, G. B.: Optimizing the grinding cycle by adaptive control,
Mashinostroitel, 3, 12–14, 1979.</mixed-citation></ref>
      <ref id="bib1.bib16"><label>16</label><mixed-citation>Malkin, S.: Grinding cycle optimization, CIRP Ann., 30, 223–226, <ext-link xlink:href="https://doi.org/10.1016/S0007-8506(07)60930-5" ext-link-type="DOI">10.1016/S0007-8506(07)60930-5</ext-link>, 1981.</mixed-citation></ref>
      <ref id="bib1.bib17"><label>17</label><mixed-citation>Maris, M., Snoeys, R., and Peters, J.: Analysis of plunge grinding
operations, CIRP Ann., 241, 225–230, <ext-link xlink:href="https://doi.org/10.1007/978-1-4613-1965-8_6" ext-link-type="DOI">10.1007/978-1-4613-1965-8_6</ext-link>, 1975.</mixed-citation></ref>
      <ref id="bib1.bib18"><label>18</label><mixed-citation>Nathan, R. D., Vijayaraghavan, L., and Krishnamurthy, R.: Intelligent
estimation of burning limits to aid in cylindrical grinding cycle planning,
Heavy Veh. Syst., 80, 48–59, <ext-link xlink:href="https://doi.org/10.1504/IJHVS.2001.001154" ext-link-type="DOI">10.1504/IJHVS.2001.001154</ext-link>, 2001.
</mixed-citation></ref><?xmltex \hack{\newpage}?>
      <ref id="bib1.bib19"><label>19</label><mixed-citation>Nishimura, T., Inasaki, I., and Yamamoto, N.: Study on optimization of
internal grinding cycle, T. Jpn. Soc. Mech.
Eng., 55, 1808–18013, <ext-link xlink:href="https://doi.org/10.1299/kikaic.55.1808" ext-link-type="DOI">10.1299/kikaic.55.1808</ext-link>, 1989.</mixed-citation></ref>
      <ref id="bib1.bib20"><label>20</label><mixed-citation>Pereira, W. X., Diniz, A. E., and Hassu, A.: Comparing different plunge
cylindrical grinding cycles based on workpiece roughness and process
vibration, J. Braz. Soc. Mech. Sci., 31, 161–163, <ext-link xlink:href="https://doi.org/10.1590/S1678-58782009000200009" ext-link-type="DOI">10.1590/S1678-58782009000200009</ext-link>, 2009.</mixed-citation></ref>
      <ref id="bib1.bib21"><label>21</label><mixed-citation>Pereverzev, P. P. and Akintseva, A. V.: Automatic cycles
multiparametric optimization of internal grinding, Elsevier, Procedia
Engineer., 129, 121–126, <ext-link xlink:href="https://doi.org/10.1016/j.proeng.2015.12.019" ext-link-type="DOI">10.1016/j.proeng.2015.12.019</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bib22"><label>22</label><mixed-citation>Pereverzev, P. P. and Akintseva A. V.: The problem resolution of
designing of cycles in the conditions of modern automated production, Russian Engineering Research,
37, 523–529, <ext-link xlink:href="https://doi.org/10.3103/S1068798X1706017X" ext-link-type="DOI">10.3103/S1068798X1706017X</ext-link>, 2017.</mixed-citation></ref>
      <ref id="bib1.bib23"><label>23</label><mixed-citation>Pereverzev, P. P. and Pimenov D. Y.: Optimization of control programs
for numerically controlled machine tools by dynamic programming, Russian Engineering Research, 35,
135–142, <ext-link xlink:href="https://doi.org/10.3103/S1068798X15020197" ext-link-type="DOI">10.3103/S1068798X15020197</ext-link>,
2015.</mixed-citation></ref>
      <ref id="bib1.bib24"><label>24</label><mixed-citation>Petrakov, Y. and Chamata, S.: Control of grinding the mandrel working
surface of cold-rolling mills, East. Eur. J. Enterpr.
Technol., 2, 55–63, <ext-link xlink:href="https://doi.org/10.15587/1729-4061.2015.39042" ext-link-type="DOI">10.15587/1729-4061.2015.39042</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bib25"><label>25</label><mixed-citation>Pereverzev, P. P. and Pimenov, D. Y.: A grinding force model allowing
for dulling of abrasive wheel cutting grains in plunge cylindrical grinding,
J Frict. Wear., 37, 60–65, <ext-link xlink:href="https://doi.org/10.3103/S106836661601013X" ext-link-type="DOI">10.3103/S106836661601013X</ext-link>, 2016.</mixed-citation></ref>
      <ref id="bib1.bib26"><label>26</label><mixed-citation>Pereverzev, P. P., Popova, A. V., and Pimenov, D. Yu.: Relation between
the cutting force in internal grinding and the elastic deformation of the
technological system, Russian Engineering Research, 35, 215–217, <ext-link xlink:href="https://doi.org/10.3103/S1068798X15030156" ext-link-type="DOI">10.3103/S1068798X15030156</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bib27"><label>27</label><mixed-citation>Phan, A. M., Summers, M. P., and Parmigiani, J. P.: Optimization device
for grinding media performance parameters, Int. Mechan. Eng. Congr. Expos
(IMECE), 3, 915–923, <ext-link xlink:href="https://doi.org/10.1115/IMECE2011-64210" ext-link-type="DOI">10.1115/IMECE2011-64210</ext-link>, 2011.</mixed-citation></ref>
      <ref id="bib1.bib28"><label>28</label><mixed-citation>Raoufinia, M., Petrakov, Y. V., Ataei, A., Parand, R., and
Abou-El-Hossein, K.: Error compensation of complex three-dimensional
surfaces machined on computer-numeric-control grinding machine tools, J.
Appl. Sci., 9, 1356–1361, <ext-link xlink:href="https://doi.org/10.3923/jas.2009.1356.1361" ext-link-type="DOI">10.3923/jas.2009.1356.1361</ext-link>,
2009.</mixed-citation></ref>
      <ref id="bib1.bib29"><label>29</label><mixed-citation>Rowe, W. B. and Ebbrell, S.: Process requirements for cost-effective
precision grinding, CIRP Ann., 53, 255–258,
<ext-link xlink:href="https://doi.org/10.1016/S0007-8506(07)60692-1" ext-link-type="DOI">10.1016/S0007-8506(07)60692-1</ext-link>, 2004.</mixed-citation></ref>
      <ref id="bib1.bib30"><label>30</label><mixed-citation>Uhlmann, E., Koprowski, S., Weingaertner, W. L., and Rolon, D. A.:
Modelling and simulation of grinding processes with mounted points: Part I
of II – Grinding tool surface characterization, Proc. CIRP, 46, 599–602, <ext-link xlink:href="https://doi.org/10.1016/j.procir.2016.03.205" ext-link-type="DOI">10.1016/j.procir.2016.03.205</ext-link>, 2016.</mixed-citation></ref>

  </ref-list></back>
    <!--<article-title-html>Designing optimal automatic cycles of round grinding based on the synthesis of digital twin technologies and dynamic programming method</article-title-html>
<abstract-html><p>The article presents methodology for designing optimal
feed control cycles resistant to unstable machining conditions for a batch
of parts in round grinding operations performed on CNC machines.</p><p>To improve the quality and reliability of control programs for CNC machines, a digital twin (DT) is proposed to be used. It performs virtual testing of a given grinding cycle for the possibility of defect occurrence at some
combination of variable technological factors. To design the optimal
grinding cycle by calculating the optimum trajectory of the radial feed
change cycle, we propose the dynamic programming method (DPM) used in
solving the classical transport problem, in which a network of roads with
intermediate stations is specified. In the task it is necessary to find the
optimal trajectory of transport. The conditions for the appearance of
defects detected by the DT are introduced into the system of restrictions,
and the feed cycle is again optimized in the DPM optimization system. A new
optimized cycle is again tested by the DT. These iterations are repeated
until the DT fixes the present of a defect for the operation. The
methodology proposed in the article for the synthesis of digital twin and
DPM technologies at the stage of preparation of control programs for CNC
machines makes it possible to guarantee the maximum productivity of the
grinding operation while ensuring the specified quality of the machined
surface under varying processing conditions that vary within the specified
limits.</p></abstract-html>
<ref-html id="bib1.bib1"><label>1</label><mixed-citation>
Akintseva, A. V. and Pereverzev, P. P.: Prospects for the development
of the theory of designing optimal cycles of machining in a multidimensional
space of control parameters, EDP Sciences, 129, 01018, <a href="https://doi.org/10.1051/matecconf/201712901018" target="_blank">https://doi.org/10.1051/matecconf/201712901018</a>, 2017.
</mixed-citation></ref-html>
<ref-html id="bib1.bib2"><label>2</label><mixed-citation>
Alagumurthi, N., Panairadja, K., and Soundararajan, V.: Optimization of
grinding process through design of experiment (DOE) – a comparative study,
Mater. Manuf. Proc., 21, 19–21, <a href="https://doi.org/10.1080/AMP-200060605" target="_blank">https://doi.org/10.1080/AMP-200060605</a>, 2007.
</mixed-citation></ref-html>
<ref-html id="bib1.bib3"><label>3</label><mixed-citation>
Amitay, G., Malkin, S., and Koren, Y.: Adaptive control optimization of
grinding, J. Eng. Ind., 103, 103–108, 10.31399/asm.hb.v16.a0002175, 1989.
</mixed-citation></ref-html>
<ref-html id="bib1.bib4"><label>4</label><mixed-citation>
Bellman, R.: Dynamic programming, Foreign Literature Publishing House,
Moscow, Russia, 1960.
</mixed-citation></ref-html>
<ref-html id="bib1.bib5"><label>5</label><mixed-citation>
Bertsekas, D.: Dynamic programming and optimal control, Prentice Hall, Reading, New
Jersey, USA, 1995.
</mixed-citation></ref-html>
<ref-html id="bib1.bib6"><label>6</label><mixed-citation>
Cahill, M. J., Bechtold, M. J., Fess, E., Wolfs, F. L., and Bechtold,
R.: Ultrasonic precision optical grinding technology, P. SPIE, 96330, <a href="https://doi.org/10.1117/12.2195977" target="_blank">https://doi.org/10.1117/12.2195977</a>, 2015.
</mixed-citation></ref-html>
<ref-html id="bib1.bib7"><label>7</label><mixed-citation>
Fang, F., Cheng, K., Ding, H., Chen, S., and Zhao, L.: Sustainable
design and analysis of CNC machine tools: sustainable design index based
approach and its application perspectives, ASME Proceedings, 3, MSEC2016-8730, <a href="https://doi.org/10.1115/MSEC2016-8730" target="_blank">https://doi.org/10.1115/MSEC2016-8730</a>, 2016.
</mixed-citation></ref-html>
<ref-html id="bib1.bib8"><label>8</label><mixed-citation>
Hashimoto, F. and Lahoti, G. D.: Optimization of set-up conditions for
stability of the centerless grinding process, CIRP Ann.-Manuf. Techn.,
53, 1271–274, <a href="https://doi.org/10.1016/S0007-8506(07)60696-9" target="_blank">https://doi.org/10.1016/S0007-8506(07)60696-9</a>,
2004.
</mixed-citation></ref-html>
<ref-html id="bib1.bib9"><label>9</label><mixed-citation>
Korchak, S. N.: The Productivity of the grinding of steel parts,
Mashinostroenie, Moscow, Russia, 1974.
</mixed-citation></ref-html>
<ref-html id="bib1.bib10"><label>10</label><mixed-citation>
Krajnik, P., Drazumeric, R., Meyer, D., Kopac, J., and Zeppenfeld, C.:
Simulation of workpiece forming and centre displacement in plunge
centreless grinding, Int. J. Mach. Tool. Manu., 48, 824–831, <a href="https://doi.org/10.1016/j.ijmachtools.2007.12.008" target="_blank">https://doi.org/10.1016/j.ijmachtools.2007.12.008</a>, 2008.
</mixed-citation></ref-html>
<ref-html id="bib1.bib11"><label>11</label><mixed-citation>
Lee, C. W.: Dynamic optimization of the grinding process in batch
production, J. Manuf. Sci. E.-T. ASME, 131, 485–494, <a href="https://doi.org/10.1115/msec_icmp2008-72212" target="_blank">https://doi.org/10.1115/msec_icmp2008-72212</a>, 2008.
</mixed-citation></ref-html>
<ref-html id="bib1.bib12"><label>12</label><mixed-citation>
Lee, T. S., Ting, T. P., Lin, Y. J., and Htay, T.: A particle swarm
approach for grinding process optimization, J. Adv. Manuf.
Technol., 33, 1128–1135, <a href="https://doi.org/10.1007/s00170-006-0538-y" target="_blank">https://doi.org/10.1007/s00170-006-0538-y</a>, 2006.
</mixed-citation></ref-html>
<ref-html id="bib1.bib13"><label>13</label><mixed-citation>
Lin, L., Yong, L., Liwei, W., and Jun, X.: PMAC-based tracking control
system for 8-axis automated tape-laying machine, Chinese J. Aeronaut., 22, 558–563, 10.1016/S1000-9361(08)60141-7, 2009.
</mixed-citation></ref-html>
<ref-html id="bib1.bib14"><label>14</label><mixed-citation>
Liu, Y. M., Yang, T. Y., He, Z., and Li, J. Y.: Analytical modelling of
grinding process in rail profile correction cnsidering grinding pattern,
Arch Civ. Mech. Eng., 18, 669–678, 10.1016/j.acme.2017.10.009, 2017.
</mixed-citation></ref-html>
<ref-html id="bib1.bib15"><label>15</label><mixed-citation>
Lur'e, G. B.: Optimizing the grinding cycle by adaptive control,
Mashinostroitel, 3, 12–14, 1979.
</mixed-citation></ref-html>
<ref-html id="bib1.bib16"><label>16</label><mixed-citation>
Malkin, S.: Grinding cycle optimization, CIRP Ann., 30, 223–226, <a href="https://doi.org/10.1016/S0007-8506(07)60930-5" target="_blank">https://doi.org/10.1016/S0007-8506(07)60930-5</a>, 1981.
</mixed-citation></ref-html>
<ref-html id="bib1.bib17"><label>17</label><mixed-citation>
Maris, M., Snoeys, R., and Peters, J.: Analysis of plunge grinding
operations, CIRP Ann., 241, 225–230, <a href="https://doi.org/10.1007/978-1-4613-1965-8_6" target="_blank">https://doi.org/10.1007/978-1-4613-1965-8_6</a>, 1975.
</mixed-citation></ref-html>
<ref-html id="bib1.bib18"><label>18</label><mixed-citation>
Nathan, R. D., Vijayaraghavan, L., and Krishnamurthy, R.: Intelligent
estimation of burning limits to aid in cylindrical grinding cycle planning,
Heavy Veh. Syst., 80, 48–59, <a href="https://doi.org/10.1504/IJHVS.2001.001154" target="_blank">https://doi.org/10.1504/IJHVS.2001.001154</a>, 2001.

</mixed-citation></ref-html>
<ref-html id="bib1.bib19"><label>19</label><mixed-citation>
Nishimura, T., Inasaki, I., and Yamamoto, N.: Study on optimization of
internal grinding cycle, T. Jpn. Soc. Mech.
Eng., 55, 1808–18013, <a href="https://doi.org/10.1299/kikaic.55.1808" target="_blank">https://doi.org/10.1299/kikaic.55.1808</a>, 1989.
</mixed-citation></ref-html>
<ref-html id="bib1.bib20"><label>20</label><mixed-citation>
Pereira, W. X., Diniz, A. E., and Hassu, A.: Comparing different plunge
cylindrical grinding cycles based on workpiece roughness and process
vibration, J. Braz. Soc. Mech. Sci., 31, 161–163, <a href="https://doi.org/10.1590/S1678-58782009000200009" target="_blank">https://doi.org/10.1590/S1678-58782009000200009</a>, 2009.
</mixed-citation></ref-html>
<ref-html id="bib1.bib21"><label>21</label><mixed-citation>
Pereverzev, P. P. and Akintseva, A. V.: Automatic cycles
multiparametric optimization of internal grinding, Elsevier, Procedia
Engineer., 129, 121–126, <a href="https://doi.org/10.1016/j.proeng.2015.12.019" target="_blank">https://doi.org/10.1016/j.proeng.2015.12.019</a>, 2015.
</mixed-citation></ref-html>
<ref-html id="bib1.bib22"><label>22</label><mixed-citation>
Pereverzev, P. P. and Akintseva A. V.: The problem resolution of
designing of cycles in the conditions of modern automated production, Russian Engineering Research,
37, 523–529, <a href="https://doi.org/10.3103/S1068798X1706017X" target="_blank">https://doi.org/10.3103/S1068798X1706017X</a>, 2017.
</mixed-citation></ref-html>
<ref-html id="bib1.bib23"><label>23</label><mixed-citation>
Pereverzev, P. P. and Pimenov D. Y.: Optimization of control programs
for numerically controlled machine tools by dynamic programming, Russian Engineering Research, 35,
135–142, <a href="https://doi.org/10.3103/S1068798X15020197" target="_blank">https://doi.org/10.3103/S1068798X15020197</a>,
2015.
</mixed-citation></ref-html>
<ref-html id="bib1.bib24"><label>24</label><mixed-citation>
Petrakov, Y. and Chamata, S.: Control of grinding the mandrel working
surface of cold-rolling mills, East. Eur. J. Enterpr.
Technol., 2, 55–63, <a href="https://doi.org/10.15587/1729-4061.2015.39042" target="_blank">https://doi.org/10.15587/1729-4061.2015.39042</a>, 2015.
</mixed-citation></ref-html>
<ref-html id="bib1.bib25"><label>25</label><mixed-citation>
Pereverzev, P. P. and Pimenov, D. Y.: A grinding force model allowing
for dulling of abrasive wheel cutting grains in plunge cylindrical grinding,
J Frict. Wear., 37, 60–65, <a href="https://doi.org/10.3103/S106836661601013X" target="_blank">https://doi.org/10.3103/S106836661601013X</a>, 2016.
</mixed-citation></ref-html>
<ref-html id="bib1.bib26"><label>26</label><mixed-citation>
Pereverzev, P. P., Popova, A. V., and Pimenov, D. Yu.: Relation between
the cutting force in internal grinding and the elastic deformation of the
technological system, Russian Engineering Research, 35, 215–217, <a href="https://doi.org/10.3103/S1068798X15030156" target="_blank">https://doi.org/10.3103/S1068798X15030156</a>, 2015.
</mixed-citation></ref-html>
<ref-html id="bib1.bib27"><label>27</label><mixed-citation>
Phan, A. M., Summers, M. P., and Parmigiani, J. P.: Optimization device
for grinding media performance parameters, Int. Mechan. Eng. Congr. Expos
(IMECE), 3, 915–923, <a href="https://doi.org/10.1115/IMECE2011-64210" target="_blank">https://doi.org/10.1115/IMECE2011-64210</a>, 2011.
</mixed-citation></ref-html>
<ref-html id="bib1.bib28"><label>28</label><mixed-citation>
Raoufinia, M., Petrakov, Y. V., Ataei, A., Parand, R., and
Abou-El-Hossein, K.: Error compensation of complex three-dimensional
surfaces machined on computer-numeric-control grinding machine tools, J.
Appl. Sci., 9, 1356–1361, <a href="https://doi.org/10.3923/jas.2009.1356.1361" target="_blank">https://doi.org/10.3923/jas.2009.1356.1361</a>,
2009.
</mixed-citation></ref-html>
<ref-html id="bib1.bib29"><label>29</label><mixed-citation>
Rowe, W. B. and Ebbrell, S.: Process requirements for cost-effective
precision grinding, CIRP Ann., 53, 255–258,
<a href="https://doi.org/10.1016/S0007-8506(07)60692-1" target="_blank">https://doi.org/10.1016/S0007-8506(07)60692-1</a>, 2004.
</mixed-citation></ref-html>
<ref-html id="bib1.bib30"><label>30</label><mixed-citation>
Uhlmann, E., Koprowski, S., Weingaertner, W. L., and Rolon, D. A.:
Modelling and simulation of grinding processes with mounted points: Part I
of II – Grinding tool surface characterization, Proc. CIRP, 46, 599–602, <a href="https://doi.org/10.1016/j.procir.2016.03.205" target="_blank">https://doi.org/10.1016/j.procir.2016.03.205</a>, 2016.
</mixed-citation></ref-html>--></article>
