An Improved Method of Minimizing Tool Vibration during Boring Holes in Large-Size Structures

The paper presents a thoroughly modified method of solving the problem of vibration suppression when boring large-diameter holes in large-size workpieces. A new approach of adjusting the rotational speed of a boring tool is proposed which concerns the selection of the spindle speed in accordance with the results of the simulation of the cutting process. This streamlined method focuses on phenomenological aspects and involves the identification of a Finite Element Model (FEM) of a rotating boring tool only and validating it with a real object, while dispensing with discrete modelling of a completely rigid workpiece. In addition, vibrations in the boring process in all directions were observed, which implies a geometric nonlinearity of the process model. During the simulation, the values of the Root Mean Square (RMS) of the time plots and the dominant values of the “peaks” in the displacement amplitude spectra were obtained. The effectiveness of the method was demonstrated using a selected mechatronic design technique called Experiment-Aided Virtual Prototyping (E-AVP). It was successfully verified by measuring the roughness of the indicated zone of the workpiece surface. The economic profitability of implementing the method in the production practice of enterprises dealing with mechanical processing is also demonstrated.


Introduction
The relative vibrations between the tool and the workpiece are recognized as the main cause of various problems detected during large-size structures machining [1]. Under certain circumstances occurring in the boring process, they can lead to a loss of system stability and the appearance of self-excited chatter vibrations, due to process and structure dynamic interactions resulting in modulated chip thickness [2]. These can cause a reduction in the overall machine tool performance or the quality of the workpiece surface, increased tool wear, especially when boring large diameter holes [3]. In some extreme cases it can even lead to the destruction of the tool or workpiece. To increase the milling efficiency while maintaining the surface quality and minimizing the vibration level, parameters of the machining process should be adjusted [4].
There are many scientific studies devoted to the problems of vibration reduction in the boring process [5]. For example, chatter is a limiting factor during boring of deep holes with long slender boring bars and may be effectively damped by the magnetorheological damper at different machining conditions for boring of Inconel 718 and Al 7075 workpieces, but the method was presented only for a case of small boring tool (length 300 mm, diameter 25 mm) [6]. In [7] Miguélez et al. presented the behavior of boring bars, which were modelled as cantilever Euler-Bernoulli beams, with a passive dynamic vibration absorber for chatter suppression. Only the first mode of vibration of the bars was considered. Bansal and Law presented (in [8]) a receptance coupling-based method to optimally tune and place a tuned mass damper on a similar slender boring bar. Chatter-free depth of cut for boring was increased a few times, which is a significant improvement. The considerations concerned a rotating workpiece, the spindle speed of which has not been identified. In [9] Yadav et al. described a receptance coupling approach in which the substructural receptances of the boring bar modelled as a cantilevered Euler-Bernoulli beam are combined with the substructural receptances of a damper modelled as a rigid mass integrated anywhere within the bar. The improvement of the damping capability of boring tools by using impact dampers was investigated in [10]. Chatter vibration elimination during deep hole boring using slender bars was considered in [11,12]. Composite material may have much higher stiffness and damping in comparison to metal, and tests of such properties for chatter suppression were investigated in [13][14][15][16]. The design and testing of a new dynamic system for boring, based on a tuned holder that increases the dynamic stiffness of slender boring bars by matching the holder natural frequency to the clampedfree boring bar natural frequency, reduced susceptibility to chatter [17]. However, the considerations made there were not verified in the real machining process. In [18] design of the low-cost damped tool holder using different types of high-density materials such as copper and brass to suppress the chatter of boring operations is investigated. The design of an anisotropic boring tool as a chatter suppression method was presented in [19,20]. The obtained theoretical and experimental results showed that using piezoelectric shunt damping it is possible to significantly increase the stability margin in boring operations during tests limited to small boring bars and rotating parts with undefined spindle speed [2]. An adaptive sliding mode control approach is presented to suppress the chatter in the boring process in the presence of uncertainties of model and dynamics [21]. The industrialized version of the boring bar with embedded sensors gives valuable insight into the cutting process during which chatter and excessive insert wear can be detected [22]. The vibration measurement gives an indication of the quality of the machined surface and thus shows the potential of the technology. The results of machining in both Maraging 250 and steel at a cutting speed of 120 m/min, a hole diameter of about 90 mm and a slender bar, whose length is unidentified, are presented. In order to meet the requirements of boring holes about 22 mm diameter with rotating bars, overhang lengths from 160 mm to 420 mm, at high rotational speed (2500 rpm), a boring bar was designed and manufactured with high stiffness pitch-based carbon fiber epoxy composite whose parameters were experimentally determined [23]. Its dynamic characteristics were measured by boring aluminum engine blocks, and the dynamic stiffness of such a bar was increased by approximately 30%. Chatter did not occur when the ratio of length and diameter was about 30% greater than for the tungsten bar. The influence of different inner cores on the dynamic behavior and technical capabilities of non-rotating boring bars (250 mm × 25 mm) during machining of a rotating workpiece with an inner diameter of 60 mm was successfully evaluated [24]. The limit of stability depends not only on the mechanical properties of the boring bar, but also on its fixation and on the machine tool. Performance limiting machining parameters can be raised, thereby improving the productivity of machining operations. In [25] the vibration stability of the boring process is discussed, in which the workpiece rotates, a medium-size bar (0.3 m × 0.020 m) is fixed, and a passive dynamic vibration absorber is used. The boring bar is modelled as an Euler-Bernoulli cantilever beam and the absorber is considered to be attached by a spring and a damper at a certain cross-section of the beam. In order to determine the optimal values of the absorber parameters, the criterion was the maximization of the minimum values of the stability lobe diagram and the classic Nelder-Mead method was used for unconstrained optimization problems. However, no machining test was performed. The stability limit was strongly improved by increasing it by about 15% and the stability lobes were only shown analytically. Active damping of the boring bars with an in-house designed magnetic actuator offers good prospects for practical use in the case of large boring tools [26]. A boring bar with an internal friction damper has also been proposed to reduce chatter vibration of the boring bars [27]. This simple friction damper consisted of several pins axially mounted inside the bar, which caused resistance and energy dissipation during bending vibrations. Unfortunately, users of large-size boring tools have practically no possibility to modify the elastic-damping properties of these tools on their own.
A separate investigation concerns the mechanics of single-point boring operations used in the industry, due to its valuable nature in modelling the chip thickness distribution along the cutting edge [3]. Cutting forces are correlated to the chip area using mechanistic cutting force coefficients. The latter makes it possible to predict them satisfactorily in three Cartesian directions. The model is an essential foundation for examining not only chatter, but also forced vibrations in boring processes. A mathematical model of the cutting force system was developed when boring holes with large diameters and boring head with multiple inserts is spinning [28]. The model includes tool geometry, chip load, cutting edge contact length and process parameters, kinematics and mechanics of the boring process, as well as axial and radial runouts. When runout occurs in the inserts the cutting forces normal to the axis of the hole become periodic at the tooth passing frequency. The model has been experimentally verified with success, so that the surface finish and dimensional quality of the holes are maintained by avoiding excessive forced vibrations. In [29] an analytical model for the stability of turning and boring processes is proposed, which includes the dynamics of the system in a multi-dimensional form with emphasis on the true process geometry and the insert nose radius in a precise manner.
A thorough and exhaustive review of publications in the field of production engineering proves that despite many scientific reports on the dynamics of the boring process, an insufficient number is noted, especially in the case of large-size workpieces. Moreover, scientific research is far away from the requirements in this field [30]. The reason is that academia usually does not have access to the requirements of this niche industry, and the topic, as a cross-sectorial one, is not of particular interest to academia due to the limited possibilities for in-depth research. A high level of investment is also required to have large machine tools in the laboratories. Therefore, most advances in machine tools for large parts have been achieved thanks to solutions proposed by and implemented in close cooperation with the industry.
Large-scale machining is increasingly becoming the domain of small and mediumsized enterprises that do not have their own research facilities. Hence, their interests do not concern the implementation of costly integrated computerized system solutions, but the ad hoc optimization of individual technological activities at the level of process implementation. After all, the production of large-size products is usually estimated at a small number of details and due to the expected profitability-a separate approach is required for each type of product.
The above considerations indicate the only direction of effective solving of the problems of boring large-size objects, which is the performance of basic and applied research in strict cooperation with an industrial partner. The subject of the article is a method of searching for conditions for minimizing the vibration level of a rotating tool during the boring process in large-size products, based on the results of previous computer simulations of the identified boring tool model. The article in question is a creative, progressive and thorough modification of the considerations presented earlier in [31]. However, unlike the previous one, a new computational model of the boring process was proposed in the phenomenological aspect, in which: discrete modeling of an inherently rigid workpiece was abandoned; -the focus was put on the discrete modeling of a boring tool with the use of a set of flexible bars, and on the assessment of model compliance with a real object; -vibrations in the boring process were observed in all directions, and the cutting process itself shows the features of geometric nonlinearity due to the dynamic change in the thickness and width of the layer.

Simulation Model
The process of boring holes in a perfectly rigid large-size workpiece mounted on a machine tool table ( Figure 1) using a rotating flexible boring tool ( Figure 2) can be considered in the convention of a non-stationary Finite Element Model (FEM) [4,32,33]. When boring holes of large diameters, the design of the tool is more complicated and cannot be modelled with a single bar. Thus, in the appropriate machining process model (Figure 3), it is possible to distinguish structural system, i.e., non-stationary discrete model of rotating indexable tool with a desired spindle speed n and linearly displacing with a desired feed speed v f , as a set of 4 flexible finite elements, the Euler-Bernoulli Bars (E-BBs) [32], Spring-Damping Elements (SDEs) [34], and cutting process, i.e., Coupling Element (CE) no. l [32,35]. Hence, the behavior of the structural system obtained in this way, unlike the hybrid system considered earlier [31], is described by the vector of its generalized coordinates q in local Cartesian coordinate systems x e1 , x e2 , x e3 of every E-BB no. e. The relatively small (much below 1000 rpm) values of the boring tool rotational speed allow for disregarding the influence of the gyroscopic effects related to spinning on its dynamic properties [32,36].   Moreover, the scheme of the boring process shows: rake angle γ 0 and clearance angle α 0 , as elements of cutting edge geometry in the orthogonal plane of the edge, -cutting edge angle κ r , -force F yl1 , acting in the direction of the cutting speed v c , -the time-varying thickness of the cutting layer h l , and the force acting in its direction-F yl2 , -the time-varying width of the cutting layer b l , and the force acting in its direction-F yl3 .
Within the range of the cutting speed of 100-240 m/min and the changed depth of cutting in the range of 0.625-1.750 mm, the cutting forces in boring have a linear dependency with the chip area [3]. The latter allows for assumptions about the proportional model of the cutting process [37,38]. Moreover, the cutting forces are separated into tangential and friction directions, but the friction force is further projected into the radial and feed directions [3]. Hence, taking into account the changes in the thickness h l (t) and the width b l (t) of the cutting layer over time, the components of the cutting forces acting in the direction of the Cartesian coordinates y l1 , y l2 , y l3 of CE no. l, were obtained in the forms: where and: b D -desired cutting layer width (b D = a e /sin κ r ); ∆b l (t)-dynamic change in cutting layer width for CE no. l; h D -desired cutting layer thickness (h D = f z sin κ r ); ∆h l (.)-dynamic change in cutting layer thickness for CE no. l; k dl -average dynamic specific cutting pressure for CE no. l; µ l2 , µ l3 -cutting force ratios for CE no. l, respectively as the quotient of the forces F yl2 and F yl1 , and the forces F yl3 and F yl1 ; τ l -time-delay between the same position of CE no. l and of CE no. l-1; a e -radial depth of cutting; f z -feed per tooth (f z = v f /(nz)); z-number of boring tool teeth. Dependencies (1)- (3) show that due to the presence of nonlinear constraints in the form of inequalities, the components of cutting forces have the character of strong geometric nonlinearities.
Description of cutting forces for CE no. l in six-dimensional space takes the form: where: and: q zl (t)-relative displacement of edge tip and workpiece along direction y l1 at instant of time t; q zl (t − τ l )-relative displacement of edge tip and workpiece along direction y l1 at instant of time t-τ l . Components of generalized displacements ∆w l (t) in the Cartesian coordinate system y l1 , y l2 , y l3 of CE no. l are related to the vector of generalized displacements q of the structural system using the time-dependent constraints equation [32]: where: T l (t)-transformation matrix of displacements vector q from the x e1 , x e2 , x e3 coordinates of E-BBs, e = 1, . . . , 4, to the coordinate system y l1 , y l2 , y l3 of CE no. l [32,35,39]. The equation of the dynamics of the non-stationary model of the cutting process, after taking into account the expression (11), has the final form: where :M, L, K-matrices of inertia, damping and stiffness of the set of E-BBs and accompanying SDEs, i l -number of "active" coupling elements, i.e., cutting edges currently being in contact with the workpiece. The way of determining these matrices is shown in the Appendix A. The dependence of the matrices D Pl and D Ol on the vector of generalized displacements q results in another argument proving the geometric nonlinearity of the model of dynamics of the boring process, described by Equation (12). It should be noted that the latter equation describes the behavior of the structural system in generalized coordinates and shows significant advantages over the description of the behavior of the system considered in hybrid coordinates [31].

Selecting the Best Spindle Speed
The proposed method of supervising the boring process by Experiment-Aided Virtual Prototyping (E-AVP) determines the best rotational speed of the tool based on the simulation of the boring process, carried out for a computational model of the tool tuned to the results of modal tests, along with the adopted model of the cutting process ( Figure 4).  [4,32] was not applied here. The reason is that the positions of the E-BBs nodes are, in general, different from the locations of the accelerometers used in the experiment. After taking into account the time-varying positions of the edges (CE no. l) of the boring tool, a non-stationary computational model of the machining process is created. It is impossible to estimate the parameters of the computational model of the boring process based on the assessment of compliance of the Root Mean Square (RMS) of the simulated plots with the values obtained during the measurements in the real process. However, due to the simplicity of the mathematical description of the adopted model of proportional cutting dynamics, which is its advantage, for the correct estimation of the parameters of the cutting process, it is enough to actually choose only 3 values of abstract meaning, i.e., average dynamic specific cutting pressure k dl and the cutting force coefficients µ 2 and µ 3 . Determining the values of these parameters is only a means to an end. The exact determination of the values sought is time-consuming, and at the same time does not significantly translate into improved results. Hence, these cutting process parameters are estimated using the so-called "mechatronic" procedure [40]. It allows, in contrast to the method presented in [4], to resign from the simulation cycle for the case of cutting with the standard spindle speed in order to adjust the values of these parameters. It is sufficient to perform a single simulation, the results of which form the basis for the selection of the best speed later. Subsequently, the permissible range of tested spindle speeds is selected and for chosen speeds from this range, the cutting process simulations are performed. It should be noted that instead of the previously considered abstract hybrid model [31], this time a non-stationary structural model of the machining process with deep physical meaning was used for the simulation. The spindle speed is selected as the best if it assures the lowest vibration level or the lowest dominant amplitudes in the spectra.
The key elements of the proposed method are the identification of model parameters of the boring tool itself (before machining), estimation of the necessary parameters of the cutting process, computer simulations made for a finite set of spindle speeds, and finally the machining process with parameters derived from the simulation results. It is an innovative model of decision support in the design of boring processes in large-size objects and is addressed to a wide group of recipients interested in the results of research on the quality and efficiency of production.

Modal Identification of Boring Tool
During the experimental investigation, the Sandvik CoroBore ® 825 XL boring tool, manufactured by Sandvik Coromant, Sandviken, Sweden, an international supplier of the cutting tools was used [41]. Modal identification of the boring tool mounted in a damped adaptor at its minimal extension (see, Figure 2) was performed in order to validate the parameters of the tool FEM, in accordance with the procedure presented in the previous section. Positions and directions of the accelerometers used during modal tests are shown in Figure 5. Impact modal tests were performed with a modal hammer and the object responses were measured with 15 accelerometers (marked as A17-A31 in Figure 5). There were 4 sets of 40 impacts made in the direction of accelerometers A18, A21, A22 and A29. For each set, the Frequency Response Function (FRF) was calculated using the H3 estimator. Examples of the FRF are illustrated in Figure 6. From the point of view of the analysis of vibrations occurring during machining, the most important are the modes of vibration in the low-frequency range, especially those characterized by significant values of displacement. For the tested boring tool, the frequencies up to 2000 Hz were analyzed. One should also pay attention to the coherence function. Seemingly, it looks low on the collective chart. However, if you look at individual channels, it turns out that the low values of coherence relate primarily to those channels that were weakly excited at a given point. For example, when excitation was acting in the direction of A18, low coherence concerns accelerometers A22 and, especially, A21, which measure accelerations in planes perpendicular to excitation direction. Figure 6 shows FRFs measured only for A18, 21 and 22 to make plots easier to analyze. During tests, FRFs for all of the 15 accelerometers were calculated and later used for modal identification. Identification of the parameters of the modal model of the boring tool was performed using the ERA (Eigenvalue Realisation Algorithm) [4,32,42] and the p-LSCFD (polyreference-Least Squares Complex Frequency Domain) [32,43,44] methods. Measurements were made by uniaxial accelerometers that were mounted to measure vibration for three orthogonal directions at 5 different points ( Figure 5). An important element facilitating the assessment of the correctness of the estimated figures is the symmetrical configuration of the structure. The displacements on one side of the boring tool should correspond to the displacements on the other side.
The values of the natural frequencies and dimensionless damping coefficients (the so-called modal damping), determined using both methods in the frequency range from 0 to 2000 Hz, are summarized in Table 1. The convergence of their values proves the correctness of the identified modal parameters.
The determined normal modes of natural vibrations in the frequency range up to 1000 Hz, using both identification methods, are of very good quality. This is evidenced by the comparison of the results of both methods with the use of the MAC criterion (Figure 7).
Graphs of the normal modes of the boring tool in the frequency range up to 1000 Hz are presented using the wireframe model in Figure 8. The presented modes are normalized in such a way that for each mode the maximum deflection from the equilibrium position is equal to one. The tool holder is located in the upper part of the drawing. The first two figures illustrate bending modes with respect to clamping the boring tool in the X-Z and Y-Z planes. The third form is the bending of the boring tool itself in the X-Z plane. The fourth figure is a "winglet" mode. The fifth and sixth forms are very similar and are bending ones in the Y-Z plane. The first torsional mode has not been identified, which may be due to its effective damping by the propulsion system.

The Calculation Model of a Boring Tool
Based on the characteristic geometrical dimensions of the CoroBore ® 825 XL boring tool, a FEM was created containing four 2-node finite elements of the Euler-Bernoulli Bar (E-BB) type (6 degrees of freedom in each node), connected by Spring-Damping Elements (SDE) [31,32] (see, Figure 3).
Using the proprietary programs developed in Fortran, the frequencies and modes of natural vibrations of the created computational model were calculated. The calculation results for the first 6 normal modes are graphically illustrated in Figure 9. These results were compared with the boring tool modal parameters obtained from the identification (Table 1, Figure 8). Despite slight differences in the values of frequencies 3 and 4, the obtained results were found to be consistent. This also applies to the calculated modes of natural vibrations. With the exception of some differences, mainly noticeable in forms 3rd and 5th, in other cases, the obtained figures generally reflect the image consistent with the experimentally determined. The above allows us to take into account the obtained computational model of the boring tool in computer simulations of the boring process.

Computer Simulations
The process of boring a hole (∅ 733.44 mm) in rigid large-size cast iron (EN-GJS-400-18-LT) workpiece, from the standard production program of the industrial partner PHS HYDROTOR Inc., Tuchola, Poland, was simulated using proprietary software developed in Fortran. Taking into account the previously postulated "mechatronic" procedure [40], on the basis of the considerations presented in the literature [3], the values of the parameters of the calculation model of the boring process were estimated as: k dl = 20 daN/mm 2 , µ l2 = 0.4, µ l3 = 0.2. Selected simulation results for 5 different pairs of technological parameter values are summarized in Table 2. Contrary to the considerations presented in [31], this time the vibrations were observed in three characteristic directions of the components of the cutting force [3], i.e., radial, feed speed and cutting speed. Thanks to creative modifications of the simulation software, the average simulation time was reduced to about 20% of the main machining time in the case of the simulation run on the Intel Core i7 6700 processor, compared to the calculations in [31]. The machining was first considered according to the standard technology of the cooperating company (i.e., n = 105 rpm, v f = 9.6 mm/min), and then the obtained results were compared with the proposed method. Then, the hole was machined for standard parameters, and subsequently-in a series of repetitive simulations-the predicted values of the best spindle speeds were determined, for which the boring process was planned to be carried out.
Thus, in Figure 10a, we can observe time plots and spectra of vibrations of the tool edge tip in the radial direction (i.e., normal to the bored hole surface), for the best spindle speed, and in Figure 10b for an adverse spindle speed. In Figure 11a, we can see time plots and spectra of vibrations of the tool edge tip in the feed speed direction, for the best spindle speed, and in Figure 11b for an adverse spindle speed. In Figure 12a, we can observe time plots and spectra of vibrations of the tool edge tip in the cutting speed direction, for the best spindle speed, and in Figure 12b for an adverse spindle speed. The starting times of the observation of the spectra were 300 s, with 2 17 samples for each plot. It is easy to see that at an adverse spindle speed, the vibrations that occur are much greater. In addition, relatively low excitation frequencies from 1.75 to 2.08 Hz following the entry of the cutting edge into the material and resulting from the adopted range of rotational speeds of the boring tool, cause-in all cases-vibrations to be stimulated with a frequency close to the first mode of natural vibrations, with the predominant vibration amplitudes in the feed speed direction.    (Table 2) in the radial direction and in the feed speed direction, we obtain the best configuration of parameters for n = 120 rpm, but the adverse configuration of parameters for n = 125 rpm. Although in the cutting speed direction, the best configuration of parameters is obtained for n = 115 rpm, the direction of these vibrations, unlike the first two, is not decisive from the point of view of machining quality (geometric accuracy, surface roughness). Hence, it can be assumed that the best-simulated spindle speed is n = 120 rpm.

Experimental Research of the Boring Process
Experimental research to evidence the accuracy of computer predictions, concerned boring a hole in a rigid large-size cast iron (EN-GJS-400-18-LT) workpiece, which was mounted on a milling and boring WHN 13-15 CNC machine, produced by TOS VARNS-DORF A.S. in Varnsdorf, Czech Republic (see, Figure 1). Boring was performed using Sandvik CoroBore ® 825 XL boring tool (see, Figures 2 and 5). In Table 3, the roughness parameters R a and R z of the hole surface at the lowest place ( Figure 1) are compared, obtained in accordance with the standard technology of the cooperating industry and on the basis of 3 modifications of technological parameters, as a consequence of previously conducted simulations (Table 2). Compared to the standard technology (i.e., n = 105 rpm, v f = 9.6 mm/min), the vibration level decreased noticeably in the case of the best simulation ( Table 2). The RMS values were reduced accordingly: in the radial direction-by 9%, in the feed speed direction-by 9%, and in the cutting speed direction-by 4%. The dominant amplitudes in the spectra remained almost the same comparing to standard machining technology. In addition, the evident reduction in vibration of the boring tool is accompanied by an increase in the spindle speed from 105 to 120 rpm (i.e., 14%), which results in a reduction of the main machining time by 0.89 min (11%) and a simultaneous increase in machining efficiency. At the same time, the measured values of the R a and R z parameters did not deteriorate (Table 3), which confirms that the requirements for maintaining the required surface quality were assured.
Although a further increase in the spindle speed to 125 rpm results in a reduction of the roughness parameters, doing so is in fact not recommended. This is because it adversely affects the durability of the cutting inserts, significantly shortening their service life. Thereby, it unreasonably increases processing costs.

Implementation Profitability Assessment
The purpose of these considerations, which make a significant contribution compared to [31], is to develop the basis for managerial insights that facilitate the accuracy of decisionmaker actions in the industrial practice of small and medium-sized enterprises that are not interested in implementing advanced system solutions, but in meeting the requirements of the profitability criterion at the technological process level of the product.
In order to obtain information on the profitability resulting from the application of the presented method, the standards of technological times were compared [45]. Let us reconsider the process of boring a hole in a rigid large-size cast iron EN-GJS-400-18-LT workpiece on the WHN 13-15 CNC table milling and boring machine in one pass. According to the standard technology, such machining has so far been performed in a unit time (approximately equal to the execution time), which is: where: l p0 -number of inspection cuts for pass, -t pk0 -time of inspection cut of pass, -t g0 -main time of pass.
For the sake of simplicity, in Formula (13), other components of the auxiliary time t p were omitted. Assuming l pk0 = 1, t pk0 = t g0 = 8.02 min., we get t j0 = 16.04 min.
In the solution for the best combination of technological parameters, the expected main time of pass is t g = 7.13 min., and there are no inspection passes (l p = 0). This gives the total unit time t j = 7.13 min, which is 56% shorter. The total time per unit for boring a hole is therefore: for standard technology-16.04 min, -in the case of technology based on the best combination of parameters-7.13 min.
In the case of the proposed technology, there is no need to increase the machining time by the theoretical and experimental modal analysis times, and by the boring simulation series time, because the above concern only the boring tool, regardless of the number of rigid workpieces mounted on the machine.
The profitability of applying the obtained results is therefore seen in the category of: minimizing the vibration level of the boring tool, thereby maintaining the required product quality, and shortening the total standard of time by 8.91 min, i.e., by 56%.

Conclusions
The results obtained during the research confirmed, in the production conditions of the industrial partner, the effectiveness of the proposed method of suppressing the vibrations of a rotating boring tool. Because of the use of the Experiment-Aided Virtual Prototyping (E-AVP), the best spindle speed was selected when boring holes in large-size structures. Due to technical impossibility, the vibrations of the workpiece were not measured on-line during machining. Nevertheless, the results of roughness measurements (i.e., R a , R z ) of the bored surface made with a profilometer after boring confirmed the accuracy of the prediction for the best selection of technological parameters n, v f for the machining process.
The present study makes it possible to predict the best conditions of the process of boring large workpieces only on the basis of a simulation of the computational model of the process in which the parameters of the boring tool were validated with a real object. The approach presented in the article, due to its uncomplicated nature, can be easily applied in the economic practice of many companies dealing with the machining of large-size objects-even those who do not have their own research facilities and their own hardware and software infrastructure.
The applied method does not require interference in the structure of the machine tool and, apart from identifying the computational model of the boring tool, no previous experimental tests are needed to simulate the boring process. Computational models of many boring tools can be prepared off-line for the selected configurations (e.g., overhang, boring diameter), and then the appropriate model can be selected to perform the simulation.
The assessment of the profitability of implementing the proposed method should be seen in the category of minimizing the vibration level of the boring tool, resulting in an improvement in the quality of workmanship, as well as a significant reduction in the production standard of the unit execution time. In addition, the cost reduction of the material to be removed cannot be overestimated due to the lack of inspection cuts that actually exist in standard technology.
k and axis x ej of E-BB no. e, i = 1, . . . , 3, j = 1, . . . , 3, N ek -matrix of shape functions of E-BB no. e, determined for the coordinates of the connection point of SDE no. k and E-BB no. e.