Abstract
A new slip-line field model and its associated hodograph of rounded-edge cutting tool were developed for orthogonal micro-cutting operation using matrix technique. The new model considers the existence of dead metal zone in front of the rounded-edge cutting tool. The ploughing forces, chip up-curl radii, chip thicknesses, primary shear zone thicknesses and lengths of bottom side of the dead metal zone are obtained by solving the model depending on the experimental resultant force data. The effects of cutting edge radius, uncut chip thickness, cutting speed and rake angle on these outputs are specified.
Introduction
The effects of cutting edge geometry should be known to understand the precision of micro- and nano-machining operations. In the previous works of Merchant, 1 Lee and Shaffer 2 and Oxley 3 , the cutting tool edge was assumed to be perfectly sharp to simplify the modelling of metal cutting. Besides, Maity and Das4–7 carried out theoretical slip-line field solutions of perfectly sharp cutting tool for orthogonal metal cutting, and researchers considered that the sticking and slipping contact occurred at the chip–tool interface. Fang et al. 8 developed a universal slip-line model of two-dimensional (2D) machining for restricted perfectly sharp cutting tools, and the developed model provided the effects of the chip up-curling and the chip back-flow. Fang9,10 presented a slip-line model for a perfectly sharp cutting tool, and this model considered the tool–chip contact on the tool secondary rake face. The forces, chip thickness and tool–chip contact length can be determined using the proposed slip-line model. Das et al. 11 proposed a slip-line field model for perfectly sharp cutting tools which accepts the adhesion friction at the chip–tool interface by using Kudo’s 12 basic slip-line field. Fang and Dewhurst 13 developed a slip-line model for perfectly sharp cutting tools with built-up edge formation. The model predicted the length and height of the built-up edge, cutting and thrust forces, chip up-curl radius, chip thickness and tool–chip contact length. Das and Dundur 14 developed slip-line field solutions for metal cutting with perfectly sharp tools. Ye et al. 15 presented a slip-line model for perfectly sharp tool of pressure-sensitive material and examined the effect of the material pressure sensitivity on the cutting operation.
However, Albrecht 16 reported that a perfectly sharp cutting tool edge cannot be produced, and the nature radius of the cutting edge was determined to be 0.007 mm. Therefore, various slip-line analyses were carried out using the rounded-edge cutting tools. Abdelmoneim and Scrutton 17 and Abdelmoneim 18 developed a slip-line model for rounded-edge cutting tools, and this model considered the effects of dead metal zone in mathematical formulation. Waldorf et al. 19 presented a slip-line model of the ploughing process for rounded-edge cutting tool in orthogonal machining. The model was based on the assumption of a dead metal zone in front of the finite radius of cutting edge. Fang20,21 proposed a slip-line model to analyse orthogonal metal cutting with a rounded-edge cutting tool. According to the researcher, the tool edge roundness was defined by the tool edge radius, position of the stagnation point and tool–chip frictional shear stress. Fang and Fang 22 investigated the finish cutting with a rounded-edge tool by both theoretical and experimental studies. The two models, including an analytical model and a finite element (FE) model, considered the effects of strain, strain rate, cutting tool edge radius and temperature on the material flow stress. Wang and Jawahir 23 developed a slip-line model for rounded-edge cutting tool–restricted contact grooved tools. Karpat and Ozel 24 presented the cutting mechanics with curvilinear edge cutting tool in the presence of dead metal zone to determine the effects of edge geometry and cutting conditions on the friction and the cutting tool temperature. An analytical slip-line field model was performed to investigate the cutting mechanics and friction at the tool–chip and tool–workpiece interfaces. Jin and Altintas 25 developed a slip-line field model for the micro-cutting process with a rounded-edge cutting tool. Uysal et al. 26 presented a slip-line model of the rounded-edge tool for metal cutting. The model considered the existence of a dead metal zone and consisted of seven sub-regions. Ozturk and Altan 27 presented a slip-line model for a rounded-edge tool. The model was divided into eight sub-regions, and the slip-line angles were calculated depending on the experimental force data. Ozturk et al. 28 metallographically analysed the dead metal zone during metal cutting with the rounded-edge cutting tool. Ozturk and Altan 29 experimentally investigated the effect of the separation point on metal flows and the material behaviour in front of the tool face with a rounded-edge cutting tool.
In the Fang’s20,21 and Wang and Jawahir’s 23 studies, the rounded-edge cutting tool was divided into two straight chords to simplify the models, and the existence of a dead metal zone was neglected. In previous studies, some researchers assumed the cutting tool as perfectly sharp, and some researchers neglected the dead metal zone. It is known that the cutting tools have a little cutting edge radius, and in some cases, they are specially produced with a large cutting edge radius. This cutting edge radius causes the dead metal zone.
In this study, a new slip-line field model and its associated hodograph were developed for a rounded-edge cutting tool by considering the dead metal zone formation. The mathematical formulation of the new model was determined by employing the Dewhurst and Collins’ 30 matrix technique.
The new slip-line field model
The new slip-line model for orthogonal cutting with a rounded-edge cutting tool was developed by assuming a dead metal zone formation. The new slip-line model and its hodograph can be seen in Figures 1 and 2, respectively, and this model has a dead metal zone (GSN) in front of rounded-edge cutting tool. The new slip-line model consists of eight slip-line angles (θ 1, θ 2, η 1, η 2, δ, δ 2, δ 3 and ψ) and two angles of vertices (α 1 and α 2). The angle α 1 is between the free surface of workpiece and the straight A3A2 slip-line, and the angle α2 is between the free surface of chip and the curved IA slip-line. Besides, the new model includes the angle δ 1 between the straight AA3 boundary and the cutting speed. If the slip-line angle ψ = 0, the region SCIJD disappears, and if the slip-line angle δ = 0, the region GKL disappears. In Figure 1, point S is called as separation point, and material flows upwards or downwards from point S.

The new slip-line model for a rounded-edge cutting tool.

Hodograph of the new slip-line model.
The new slip-line model follows two traditional slip-line theory assumptions. The first is that the material deformation is under plane-strain conditions. For this reason, the new model is suitable for orthogonal metal cutting with continuous chip formation. The other assumption is that the workpiece material is rigid-plastic. In addition, it is accepted that the static forces are stable in the slip-line field.
Chip formation consists of the primary, secondary and tertiary shear zones; three transition regions; a material pre-flow region and a dead metal zone, as seen in Figure 3. The metal cutting models usually neglect the existence of these transition regions and the dead metal zone, but the shapes, sizes and locations of these transition regions and the dead metal zone depend on cutting conditions, tool geometries and workpiece materials.

Shear zones and regions of the new slip-line model.
All 14 slip-line sub-regions of the new slip-line model of rounded-edge cutting tool for orthogonal machining are shown in Figure 4. Sub-region 1 is a non-uniform field with a straight boundary and formed by tool–chip friction on the tool rake face. In sub-region 1, the material flows from point G to point H and the tool–chip contact length is defined. Sub-region 6 is called as dead metal zone, and it occurs in front of the rounded-edge cutting tool. In sub-region 6, the material does not flow. The non-uniform sub-regions 2, 3 and 4 connect sub-regions 1 and 6. The non-uniform sub-regions 5 and 7 indicate the transition between the primary and secondary shear zones. Sub-region 8 is the tertiary shear zone. The non-uniform sub-region 9 defines the transition between the primary and secondary shear zones. The non-uniform sub-regions 10 and 11 represent the interaction between the primary and secondary shear zones. The non-uniform sub-region 12 provides the chip up-curling effect. The non-uniform sub-region 13 specifies the transition between the primary shear zone and the material pre-flow region, and the uniform sub-region 14 is formed from material pre-flow effect (Figure 4).

Slip-line sub-regions of the new slip-line model.
The new slip-line field model considers the cutting edge radius and the dead metal zone, and this model is adapted for cutting tools used in machining operations because perfectly sharp cutting tools cannot be produced and they have a little cutting edge radius. In some cases, cutting tools are specially produced with a large cutting edge radius. This cutting edge radius causes the dead metal zone. In most of the classic models, the chip formation occurs on a single shear plane. In the presented model, the chip deformation occurred in a shear zone differently from the classic models. This approach is much closer to practical metal cutting operations than a shear plane model. The new model considers the chip up-curling effect caused by the convex IA slip-line and the material pre-flow effect (slip-line sub-region AA 2 A 3, Figure 3) in the shear zone in contrast to the previous models. The new model clearly presents a strong interaction between the primary and secondary shear zones. As an example, a change in the slip-line sub-regions 1–5 and 7 in the secondary shear zone causes a change in the slip-line sub-regions 10 and 11 in the primary shear zone (Figure 4).
Mathematical formulation of the new slip-line field model
The matrix technique developed by Dewhurst and Collins 30 has been employed for numerically solving the mathematical formulation of the new slip-line model. This technique was expressed in Appendix 3. Also, Fang’s20,21 work has been used to establish and solve the new slip-line field model.
Determination of the tool–chip friction
The tool–chip friction occurs on the tool rake face and is determined according to equation (1) 20
where ζ is the factor for tool–chip friction on the tool rake face, τ is the corresponding tool–chip frictional shear stress and k is the material shear flow stress.
Determination of the slip-lines in the secondary shear zone
The slip-lines HL and KJ seen in Figure 1 are two base slip-lines, and their radii are indicated as column vectors σ
1 and σ
2, respectively. In the following equations,
Relationships among the slip-lines in the proposed new slip-line model
In the HLG slip-line region, GL is calculated as
In the LKG slip-line region, KG and KL are given by
In the KJDG slip-line region, GD and JD are expressed as
In the GDS slip-line region, GS and SD are represented by
In the DCS slip-line region, CD and CS are calculated as
In the JICD slip-line region, DJ, DC, IJ and IC are determined by
The convex IA slip-line is calculated by
where ρ is the magnitude of total velocity jump across the slip-lines AICSN, ω is the angular velocity of a curled chip rotation and c is a column vector representing a unit circle as given in equation (16)
Relationships among the slip-lines in the hodograph of the proposed new slip-line model
Equations (17)–(22) are determined by Green 31
In the i′jd′c′ slip-line region, jd′ and c′d′ are represented by
In the c′d′s′s slip-line region, c′s, d′s′ and ss′ are determined by
In the d′gs′ slip-line region, s′d′, gs′ and gd′ are calculated as
In the gd′jk slip-line, d′j, d′g, kg and kj are given by
In the lg′gk slip-line region, g′l and g′g are expressed as
In the hg′l slip-line region, hl is determined by
Equations (2)–(37) have been employed to determine two base slip-lines HL and KJ. After slip-lines HL and KJ are determined, all slip-lines in the secondary shear zone JDSGHLKJ and the transition region ICSDJI and the IA slip-line can be determined by matrix technique.
Determinations of the slip-lines in the primary and tertiary shear zones
The slip-lines in the primary and tertiary shear zones can be determined from relationships between the slip-lines in the new model (Figure 1) and its hodograph (Figure 2), when the slip-lines in the secondary shear zone are determined. The A3A2 and B2N slip-lines are assumed as straight, and the curved slip-lines A 2 A 1, I 1 A 1, I 1 C 1, C 1 B 1 and B 1 B 2 can be calculated according to equations (38)–(42)
In the above equations, ρ′ is the magnitude of total velocity jump across the slip-lines A 3 A 2 A 1 I 1 C 1 B 1 B 2 N, ΔS is the thickness of the primary shear zone and the SN slip-line is the length of the bottom side of the dead metal zone and calculated according to equation (47) when it is assumed as straight. Besides, it is assumed that ΔS is equal to SB2. Equation (44) used to determine the length between the machined surface and the end point of the dead metal zone (separation point S shown in Figure 5) has been taken from Ozturk’s 32 work

The dead metal zone and the separation angle (θS ) in the new slip-line model.
where r is the cutting edge radius, a is the length between the machined surface and the end point of the dead metal zone (separation point S), t 1 is the uncut chip thickness, V is the cutting speed and θS is the separation angle shown in Figure 5.
Determination of the resultant and ploughing forces
The forces and the bending moments are resolved at point I and transmitted through the slip-lines IA, IJ, KJ, KL and HL as shown in Figure 6. Equations (48)–(50) must be satisfied to provide the force and bending moment equilibriums

The transmitting of forces and bending moment through slip-lines in the new model.
where the variables f 1, f 2 and f 3 are symbolized three non-linear functions. Powell’s 33 algorithm is performed to solve equations (48)–(50). Equation (51) must be satisfied to provide the force and bending moment equilibriums and to achieve the force-free chip
where w is the width of cut.
If the forces transmitted through the slip-lines A
3
A
2, A
2
A
1, I
1
A
1, I
1
C
1, C
1
B
1, B
1
B
2 and B
2
N are denoted by the vectors
The cutting force and thrust force can be determined by decomposing the resultant force in the parallel and normal directions to the cutting speed. The ploughing force (P) can be determined according to equation (53)
Constraints of stresses and determination of the chip up-curl radius and the chip thickness
According to Hill’s34,35 overstressing theory, the extension of a stress field into a rigid field is possible for a limited range of solutions of the model for each value of the tool rake angle. For an acceptable solution, the two angles of vertices (α 1 and α 2) seen in Figure 1 should not be overstressed, and this condition is mathematically expressed in equations (54) and (55)
where PA is the hydrostatic pressure at point A. The two angles of vertices (α 1 and α 2) can be determined according to equations (56)–(57)
where γ is the tool rake angle and Rui is the intermediate variable of chip up-curl radius and calculated according to equation (58)
Besides, Ru is calculated according to equation (59), and t2 is the chip thickness and determined according to equation (60)
where Vh is the chip speed.
Constraints of slip-line angles
All the eight slip-line angles (θ 1, θ 2, δ, δ 2, δ 3, ψ, η 1 and η 2) must be positive in their numerical value as given in equation (61)
In addition, the slip-line angles have the following relationships
Inputs and outputs of the new slip-line field model
Inputs and outputs of the new slip-line field model are given in Figure 7. The input parameters were used in the equations to determine the slip-lines and to calculate the resultant and ploughing forces. The tool rake angle was used in equations (56), (62) and (65), and the cutting edge radius was used in equations (43)–(45). Machining parameters such as the cutting speed, uncut chip thickness and width of cut were used in equations (44), (52), (53), (57) and (58). Experimental results were performed in equations (44)–(46), (52) and (56), and the stress boundary conditions were used in equations (1), (51)–(55) and (64).

Inputs and outputs of the new slip-line field model.
Experimental study
In this study, a shaping-based quick stop device (QSD) was used in experiments (Figure 8). The QSD was developed to collect the chip-root samples. It is necessary to effectively freeze the cutting action without substantially disturbing the state of the chip. There are two basic motions in this QSD. One of them is cutting and the other is sensitive depth-of-cut motions. Brass (CuZn30) was chosen as workpiece material. Chemical composition of the samples was Cu% = 69.73, Zn% = 30.216, Pb% = 0.006, Fe% = 0.008 and Ni% = 0.04, and the yield strength, the tensile strength, the elongation and the elastic modulus were 320 MPa, 440 MPa, 11% and 114 GPa, respectively. The dimensions of the workpieces were set as 32 × 30 × 1.5 mm3 by laser cutting and the depth of heat-affected zone was measured as 60 μm by Bulut Machine HVS-1000 micro-hardness device. Therefore, the uncut chip thicknesses were chosen more than 60 μm. Uncoated tungsten carbide (WC) cutting tool inserts were chosen as TPGN160308.

Shaping-based quick stop device.
All 54 experiments were performed with various cutting edge radii, uncut chip thicknesses, cutting speeds and rake angles. Cutting parameters are given in Table 1. Cutting edge radii were measured by the optical measuring system.
Cutting parameters.
Brass samples were machined by the shaping-based QSD. The workpiece was not cut along the length and the chip was not broken to investigate the dead metal zone. The sanding and polishing operations of the samples were performed using MetaServ 2000 Sander/Polisher device. Then, the samples were etched by etchant including 100 mL H2O, 96 mL ethanol and 59 g FeCl3, as given in Voort’s 36 study. Investigation of these etched brass samples was carried out by SOIF XJP-2 optical microscope. Dead metal zone microphotographs were observed for the rounded-edge cutting tools (Figure 9).

Microphotographs of the dead metal zone for rounded-edge cutting tools of V = 0.5 m/min, t 1 = 150 μm and γ = 2°: (a) r = 50 μm and (b) r = 100 μm.
Results and discussions
Mathematical formulations of the new slip-line field model were solved by MATLAB® software. The values of PA /k and τ/k were assumed to be 1.1 and 0.9, respectively, from Fang’s 21 and Fang and Jawahir’s 37 studies. The separation angle (θS ) was chosen as 60° from Fang’s 21 work and Ozturk’s 32 experimental study. Besides, the length between the machined surface and the end point of the dead metal zone (a) was calculated using equation (44) based on Ozturk’s 32 work. The approach to solve the new slip-line field model was taken to be 1%. All variables were solved by equalizing the resultant force obtained from the computer-controlled QSD to the resultant force obtained from the new model. The slip-line angles were obtained by solving the model and the maximum and minimum values are given in Table 2.
Maximum and minimum values of the slip-line angles.
The effects of cutting edge radius and uncut chip thickness on the ploughing force are given in Figures 10 and 11, respectively. Besides, the change in the ratio of the resultant force to the ploughing force according to the cutting edge radius and uncut chip thickness can be seen in Figures 12–15.

Effect of cutting edge radius on ploughing force, γ = 2°.

Effect of uncut chip thickness on ploughing force, γ = 6°.

Effect of cutting edge radius on the ratio of resultant force to ploughing force, V = 0.5 m/min.

Effect of cutting edge radius on the ratio of resultant force to ploughing force, t1 = 100 μm.

Effect of uncut chip thickness on the ratio of resultant force to ploughing force, V = 0.25 m/min.

Effect of uncut chip thickness on the ratio of resultant force to ploughing force, r = 50 μm.
The ploughing force increased when the cutting edge radius increased as shown in Figure 10. It was specified that the effect of ploughing force on the resultant force increased as the cutting edge radius increased (Figures 12 and 13) because higher cutting edge radius made the ploughing process of the cutting tool difficult. The change in the ratio of the resultant force to the ploughing force was affected by the cutting speed very much as the cutting edge radius decreased. It was seen that the effect of the ploughing force on the resultant force decreased slightly with an increase in the cutting speed (seen in Figure 15) and the rake angle (seen in Figures 12–15). For given cutting conditions, the ratio of resultant force to ploughing force decreased when the cutting speed increased due to the fact that the increase in the cutting speed reduced the dead metal zone formation and this caused the decrease in the ploughing force.
The ploughing force increased with a cutting speed of 0.25 and 0.5 m/min as given in Figure 11. It decreased with a cutting speed of 0.75 m/min when the uncut chip thickness increased. In addition, it was specified that the effect of ploughing force on the resultant force increased with an increase in uncut chip thickness (Figures 14 and 15).
Effects of the cutting edge radius and the uncut chip thickness on the primary shear zone (ΔS) and the length of the bottom side of the dead metal (SN) are given in Figures 16–19.

Effect of cutting edge radius on the primary shear zone thickness, V = 0.5 m/min.

Effect of uncut chip thickness on the primary shear zone thickness, r = 150 μm.

Effect of cutting edge radius on the length of the bottom side of the dead metal, V = 0.75 m/min.

Effect of uncut chip thickness on the length of the bottom side of the dead metal, r = 50 μm.
The primary shear zone thickness (ΔS) increased with the increase in the cutting edge radius, decreased slightly with the increase in the uncut chip thickness, and it did not change with the cutting speed and rake angle very much as shown in Figures 16 and 17.
The length of the bottom side of the dead metal zone (SN) increased as the cutting edge radius increased and decreased a little when the uncut chip thickness increased (Figures 18 and 19). Besides, it did not change with cutting speed very much. The size of the dead metal zone grew as the length of the bottom side of the dead metal zone increased. Therefore, it can be pointed that the reason of increasing the ploughing force with the increase in the cutting edge radius was the growing of dead metal zone size.
The effects of the cutting edge radius and the uncut chip thickness on the chip up-curl radius (Ru ) and the chip thickness (t2 ) are given in Figures 20–23.

Effect of the cutting edge radius on the chip up-curl radius and chip thickness, V = 0.5 m/min, γ = 2°.

Effect of the cutting edge radius on the chip up-curl radius and chip thickness, V = 0.75 m/min, γ = 2°.

Effect of uncut chip thickness on the chip up-curl radius and chip thickness, V = 0.75 m/min, γ = 2°.

Effect of uncut chip thickness on the chip up-curl radius and chip thickness, V = 0.75 m/min, γ = 6°.
The chip-breaking process and the chip formation were affected by the chip up-curl radius (Ru ). It was observed that the chip up-curl radius increased when the cutting edge radius increased. Thus, longer chips were obtained with an increase in the chip up-curl radius (Figures 20 and 21). Additionally, it was seen that the chip up-curl radius decreased with an increase in the uncut chip thickness, and the increase in uncut chip thickness made the chip curling difficult (Figures 22 and 23). Besides, the chip up-curl radius increased as the rake angle increased and decreased as the cutting speed increased. It was specified that the chip thickness (t2 ) decreased slightly when the cutting edge radius rose and increased when the uncut chip thickness increased. In addition, the chip thickness increased with the increase in the rake angle and decreased with the increase in the cutting speed (Figures 20–23).
Conclusion
In this study, a new slip-line field model for orthogonal machining with a rounded-edge cutting tool was developed. The proposed slip-line model simultaneously predicted the resultant force, ploughing force, primary shear zone, length of the bottom side of the dead metal zone, chip up-curl radius and chip thickness. The mathematical formulations of the new model and hodograph were developed based on the basic matrix algorithm of Dewhurst and Collins. 30 The results obtained from solving the new slip-line model and hodograph are as follows:
The mathematical formulation of the new model was calculated by equalizing the resultant force obtained from the slip-line field model to the experimental resultant force. The approach to solve the new slip-line field model was taken to be1%. Therefore, the new model results were consistent with the experimental results.
It was specified that the ploughing force efficiency on the resultant force increased with the increase in the cutting edge radius.
The ratio of the resultant force to the ploughing force was much more affected by the cutting speed with the reduction in the cutting edge radius.
It was seen that the ploughing force efficiency decreased as the rake angle and the cutting speed increased. Besides, it was specified that the ploughing force was most affected by the cutting edge radius in machining with a rounded-edge cutting tool.
The ploughing force usually increased with the increase in the uncut chip thickness. In addition, it was seen that the effect of ploughing force on the resultant force increased as the uncut chip thickness increased.
It was specified that the primary shear zone thickness (ΔS) increased with the increase in the cutting edge radius, decreased a little with the increase in the uncut chip thickness, and it did not change with the cutting speed and rake angle very much.
It was defined that the length of the bottom side of the dead metal zone (SN) increased as the cutting edge radius increased and decreased slightly when the uncut chip thickness increased. It was specified that it did not change with the cutting speed very much.
It was seen that the chip up-curl radius (Ru ) increased with an increase in the cutting edge radius and decreased with an increase in the uncut chip thickness. Furthermore, the chip up-curl radius increased as the rake angle increased and decreased as the cutting speed increased.
It was specified that the chip thickness (t2 ) decreased a little when the cutting edge radius increased and increased when the uncut chip thickness increased. In addition, the chip thickness increased with the increase in the rake angle and decreased with the increase in the cutting speed.
Footnotes
Appendix 1
Appendix 2
The computational flow chart of the new slip-line model is given in Figure 24.
Appendix 3
In this study, P, P*, Q, Q*, R and G are the basic matrix operators defined by Dewhurst and Collins 30 which were used to solve the developed slip-line model. According to the researchers’ study, a slip-line is to be thought of as being represented by the column vector of coefficients in the series expansion of its radius of curvature. Thus, the curves and flow directions forming the slip-lines can be defined by mathematical expressions.
The column vector coefficients of the curves are expressed as equation (66)
These latter components are related to those of σ 1 and σ 2 and can be written in matrix form as equation (67)
where
and
where
The significance of the P* and Q* matrix operators can be seen by considering the special case in which one of the base slip-lines (σ 2) degenerates into a point (σ 2 = 0 in equation (67)). P* and Q* are hence the matrix operators which generate the singular field on the convex side of a given slip-line. From equation (67), therefore, it is evident that the curves σ 3 and σ 4 can be regarded as the sum of certain curves in the singular fields constructed on the two base curves σ 1 and σ 2.
We now define the unstarred P and Q matrices by equation (71)
where Rθ and Rψ are reversion matrix operators. For example, Rψ reverses the intrinsic direction of a given slip-line with angular span ψ
When the intrinsic directions of σ 3 and σ 4 are opposite to those in equation (67), the mathematical expressions are represented as equation (73)
When a slip-line makes an angle λ with a rough straight boundary, the G matrix operator (equation (74)) is used
where I is the unit matrix, and J is expressed as equation (75)
Declaration of conflicting interests
The authors declare that there is no conflict of interest.
Funding
This research has been supported by Yildiz Technical University Scientific Research Projects Coordination Department. Project Number: 29-06-01-DOP01.
