Modelling and Simulation in Engineering
Volume 2011, Article ID 469262, 17 pages
http://dx.doi.org/10.1155/2011/469262
Research Article

## Slicing Cuts on Food Materials Using Robotic-Controlled Razor Blade

1Department of Mechanical and Industrial Engineering, University of Minnesota, Duluth, MN 55812, USA
2Food Processing Technology Division, ATAS Lab, Georgia Tech Research Institute, Atlanta, GA 30332, USA

Received 27 May 2011; Revised 10 August 2011; Accepted 17 August 2011

Copyright © 2011 Debao Zhou and Gary McMurray. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Cutting operations using blades can arise in a number of industries, for example, food processing industry, in which cheese, fruit and vegetable, even meat, are involved. Certain questions will rise during these works, such as “why pressing-and-slicing cuts use less force than pressing-only cuts” and “how is the influence of the blade cutting-edge on force”. To answer these questions, this research developed a mathematical expression of the cutting stress tensor. Based on the analysis of the stress tensor on the contact surface, the influence of the blade edge-shape and slicing angle on the resultant cutting force were formulated and discussed. These formulations were further verified using experimental results by robotic cutting of potatoes. Through studying the change of the cutting force, the optimal slicing angle can be obtained in terms of maximum feeding distance and minimum cutting force. Based on the blade sharpness properties and the specific materials, the required cutting force can be predicted. These formulation and experimental results explained the basic theory of blade cutting fracture and further provided the support to optimize the cutting mechanism design and to develop the force control algorithms for the automation of blade cutting operations.

#### 1. Introduction

Food cutting, such as potato or cheese cuts, is different from metal cuts because of the material deformability and shape-variability. Cutting mechanics formulation is constantly a hot research topic since it can provide useful information for cutting operations. In the literature, mainly two methods have been documented, that is, energy formulation method and stress tensor distribution method.

Metal cutting creates plastically deformed offcuts which permanently store energy, while the food offcuts permanently store almost no energy. Thus, many researchers have formulated the food cutting problems using energy method, most notably, using the fracture toughness concept [1, 2]. Using energy concepts, Atkins et al. [3, 4] have been able to explain why the cutting fracture requires smaller force when pressing and slicing compared to pressing only. Kamyab et al. [5] formulated the stress and force distribution in cheese cutting.

Stress distribution formulation can provide an alternative explanation. Two methods can be applied: potential function method and superposition method. In the potential function method [6, 7], the strain and stress tensor are expressed in terms of the space derivatives of certain airy functions in the form of biharmonic equations. Using this method, the closed form expression of the stress distribution generated by a point force acting in either normal or tangential direction to the boundary of the semi-infinite body was derived in [8], the solution to a load applied normal to an infinite half-space was given by [9, 10], and the stress distribution generated by a tangential force applied to a surface was referred as Cerruti problems, and its solution was obtained using reciprocal theorem as summarized in [1113]. Based on the derivation of the integral of the point force airy functions, Love [14] provided the integral for a rectangular area with constant normal pressure. Using the same method, the solution for first-order polynomial load applied to a rectangular surface has been completed by [15], to a triangular region has also been given by [16, 17], and to a circular area by [18]. The second method is to formulate the internal stress distribution, for example, using Hertz contact mechanics [19]. Yoshihara and Matsumoto [20] studied the shearing properties of wood using the stress analysis method. Lucas used fracture mechanics to explain the function of teeth in food cutting based on the detailed investigation of dental structure [21].

Blade sharpness is another factor that affects the cutting forces. Contact between the cutting object and the blade is an area, which can be shown from the microstructure of a knife [22]. Blade sharpness also directly influences the cutting moments and the grip forces applied by an operator such as the research performed by McGorry et al. [23]. Szabo et al. developed a procedure to establish knife-steeling schedules based on increased force due to knife dullness from repetitive use to minimize operator exertions and physical stress associated with work-related musculoskeletal disorders [24].

In this paper, mainly applied to the hard and crispy materials such as cheeses and potatoes, we will focus on using shear stresses at fracture to describe the cutting mechanism and blade sharpness. Firstly, the distribution of the external cutting force on the cutting-contact area is described in Section 2.1. From Sections 2.2 to 2.3, cutting force influence factors, such as the shape of the blade cutting edge, slicing angle, and so forth, are discussed. In Section 2.4, the relative sharpness factor concept built in [25] is generalized to be applicable to the pressing and slicing cuts. In Section 3, computation results are provided to illustrate the stress distribution during cutting. The computation results provide several purposes: (1) to show the stress distribution and to verify the analytical model; (2) to find out where the maximum stress is and to use them to evaluate the fracture location; (3) to understand the influence of the blade shape (cutting force shape) and slicing angle to the fracture force; (4) to figure out the fracture mode of blade cuts. Experimental results were provided in Section 4 to verify the formulation. Justification about the application of this algorithm on nonlinear materials are provided in Section 5. Conclusions were drawn in Section 6.

#### 2. Modeling of the Cutting Interaction

##### 2.1. Simplified Model of Pressing and Slicing Cuts on Hard-Crispy Materials

The illustration of potato cutting is shown in Figure 1. The modeling and experiments will focus on using the blade shown in Figure 1(a) to realize the robot-controlled cuts shown in Figure 1(b). The microstructure photo of the edge of a blade is shown in Figure 2(a) with the measured width of the blade cutting edge. For that particular razor blade, the edge of this blade is 850 nm with flat surface. The range of the edge width of a brand new blade is around 500 nm to 1250 nm with the shape shown in Figure 2(b).

Figure 2: Microstructure of the edge of the razor blades used in the modeling and experiments.

The microstructure of potatoes showed that the homogeneity of the particles is 500 times bigger than the blade edge width. Since the target problem is about cutting fracture, the focused problem is the stress distribution at the contact area at the moment just before the initiation of the cutting fracture.

At that moment, enough pressure has been built up beneath the blade, and the blade has a firm and solid contact with the target. Thus, the cutting interaction, such as the potato cutting shown in Figure 1, is simplified as a belt-area force acting on the surface of a semi-infinite body as illustrated in Figure 3(b), where the interaction between the blade and the target (Figure 3(a)) is simplified as a belt-area force, and the target is simplified as a semi-infinite body. Given the contact length between the blade and cutting material is and the width of the blade edge is , the contact area between the blade and the material is belt-like, that is, (, ) in direction and [, ] in direction, where point in frame is at the center of the belt area, plane is on the contact surface, -axis is the blade length direction, and -axis point into the material as shown in Figure 3(b). Normally, the width of a brand new blade () ranges from 500 nm to 1250 nm, and the contact length is normally bigger than one centimeter, roughly , that is, the contact length in direction is much bigger than that in direction. Thus, can be considered as infinite compared to .

Figure 3: Simplified cutting model.

For the symbolic expression in this paper, capital is used to express a cutting force and small is used to express a pressure variable, while subscript and superscript of a variable is used to express the variable is either point-force or point-pressure related. Subscript and superscript stands a variable in normal or direction and subscript and superscript stands a variable in tangential or direction. For example, is an ideal point force acting on the contact surface of a semi-infinite body with normal component and tangential component , while and are the cutting pressure components in the direction and direction, respectively. Subscript and superscript is used to express the variable is line-pressure related.

Force has two components: normal force and tangential force as shown in Figure 4. Then there are where is defined as the slicing angle. If we assume the pressure does not change along a line parallel to the -axis, the line-pressures on a line parallel to -axis are in normal direction and in tangential direction, there are The area pressures and can be used to calculate line pressures and and total force and as follows:

Figure 4: Cutting force acting on a semi-infinite body with slicing angle .

Since the dimension of the blade edge is around 1000 nm (1 um), there is no force sensor in such small size to measure the stress distribution. According to the edge shape of a razor blade, the cutting force distributions in both tangential and normal direction to the cutting surface have to be assumed. In this paper, they are assumed to be trapezoid as shown in Figure 5 and their mathematical expressions are summarized in Table 1.

Table 1: Pressure distribution along a line in -axis direction.
Figure 5: Force distribution profile along -axis at , that is on the contact surface.

There are Then there are

Other simplifications in this paper are as follows. (1) Relaxation and creep are ignored since we consider the instantaneous cuts. (2) The cutting force intensity is zero at the edge of the contact area. (3) The offcuts move away from the knife and no friction force acts on the side of the blade; the contact area between the blade edge and the material does not change once full contact is established.

##### 2.2. Modeling of the Stress Distribution in the Cutting Materials

The modeling starts from a stress field generated by a point force on a semi-infinite solid body. Superposition method is used to obtain the stress field generated by a belt-shape area force on a semi-infinite solid body. The procedure is Table 2. In the last row of Table 2, the functions , , , , , and can be explicitly expressed using (A.8) to (A.10) and (A.12) to (A.14) using the variables . The details are shown in the appendix.

Table 2: Area external force.

Although the modeling method uses a standard superposition which can be found in any material mechanics textbook, it is first time in the literature that the stress distribution due to a belt-shaped area force was obtained in a closed and manageable form.

Substituting and in (5) into (*) in Table 2, the stress tensor at point generated by cutting force is obtained as

The principle stresses , , can then be obtained through solving where the solutions to are the principle stresses , , and ; the solutions to (3 by 1 vector) express the direction of the principle stresses. The maximum shear stress () can then be obtained as where and are one of the numbers 1, 2, or 3.

##### 2.3. Fracture and Initialization during Blade Cutting

Tresca’s failure criterion [26] was adopted to identify fracture initialization during blade cutting. This criterion can be expressed as where is the ultimate shear strength, is the ratio of the fracture initialization force over continuous cutting force, the function returns the absolute value of the expression , and the function returns the biggest value among , , and .

During the cutting process, there are two sets of fracture: fracture initialization and postfracture. Prior to material fracture, all of the stresses in the material are due to deformation. At the moment of fracture initialization, (9) can be applied with . After the fracture initiation, postfracture criterion will apply. During postfracture, the initial failure criterion () is no longer valid. The stress concentration factor needs be determined experimentally using where and are the ultimate shear stress at fracture initialization (material property) and the ultimate stress during postfracture (determined experimentally), respectively, and is the cutting force during continuous cutting (post-fracture), is the cutting force at the moment of fracture initiation. From (6)–(9), there is where function expresses a function of variable . can be explicitly expressed using (A.8) to (A.10) and (A.12) to (A.14). At the moment of fracture initiation, the required force is obtained by substituting (11) into (9), where is the cutting force during cutting fracture and () represents the and coordinates of the fracture location, and function can be explicitly expressed using (9) and (11). is the normal cutting force and is the tangential cutting force.

In general, “sharpness of a blade” is an approximate measurement of the magnitude of the applied force in different cuts by keeping the cutting material () and cutting manners ( and ) unchanged. It may be possible to explicitly express the relationship between the cutting force and the factors, , , , , and () using (12). However, since the magnitude of and is in several hundred nanometers, it will not be practical to put a microscope in a workshop to measure the width for each blade in order to estimate its sharpness. So a relative sharpness factor is defined as follows: When setting , there is Let be the cutting force of the sharpest knife, be the cutting force of the dullest knife, and be any other cutting force. In the sharpest case, such as a brand new blade from factory, the blade relative sharpness factor is defined as The relative sharpness of certain blade can be expressed as In the dullest case (in which the maximum allowable force has to be applied in order to realize cutting), there is Then for convenience, another parameter, knife relative sharpness level, is defined as where is a user-defined integer which is used to distinguish the sharpness level of a blade and int(#) is a function to round the number # to the nearest integer. The function int(#) forces the knife sharpness level value to an integer. Substituting (15)–(17) into (18), there is It can be easily seen that the knife relative sharpness level can be determined using the cutting forces at different conditions. Since the cutting forces are measureable, a knife can be categorized into -level sharpness starting from level 1 by defining in the sharpest case.

#### 3. Simulation Analysis

In the computation, the change of the stresses with the parameters (contact length), (half width of the blade cutting edge), (load shape), and (different locations in material) are discussed, and the stress distributions in the material under different conditions have been visualized for further parametric study. By changing the cutting force profile, the influence from blade sharpness can be observed and the fracture modes are obtained. In the calculation, the coordinate values or distance data are normalized by (half width of the blade edge), the stresses generated by normal and tangential forces are normalized by the maximum pressure and , respectively, (or set = 1 mm, = 1 pa, and = 1 pa).

##### 3.1. Stress Variation with Depth

The results shown in this section are to illustrate the change of the stress distribution with the depth () in the material at different locations. The results were obtained using (*) in Table 2 by assuming , constant pressure, = 1 pa, and = 1 pa. The results of the stresses , , , , , have been calculated and only the results for stress are shown in Figure 7. In Figure 7, each curve represents the stresses change with (from 0 to or ) on the same coordinates ( from to ).

It is observed from the results of , , , , , that at the contact surface, the maximum stress may not exceed the cutting stress from the blade and the stress will decrease dramatically as increases. At the location where the coordinate (depth) is about twice of the half width of the cutting blade (at ), the force decreases roughly to half of the maximum stress. In general, the magnitude of the width of a sharp blade is about 1 um and the size of the inhomogeneity in potatoes and cheese is bigger than 0.1 mm (100 um). This means that at the depth there are inhomogeneity, the stress already reduced to very small value.

Moreover, this paper is only interested in the stress to generate the fracture. From the calculation, it is found when slicing angle is smaller than 10°, the fracture happens at the depth of or at about = 0.06–0.32 um; when slicing angle is bigger than 10°, the fracture happens at the contact surface or at the depth of = 0 um. Thus, in the microscale sense, the material can be considered as homogenous although in macrosense, the cutting material may be anisotropic or heterogeneous. When we only consider the stress to generate the cutting fracture, the material can be considered as homogenous.

Here we define the zone with at least tenth of the applied maximum stress from the blade as the effective zone. Since the depth of the effective zone in the cutting modeling is much smaller than the size of the fiber, or inhomogeneous, or anisotropic zone of the cutting materials, although the material itself is inhomogeneous or anisotropic, only in the cutting effective zone can the cutting modeling still be approximated as homogenous material.

##### 3.2. Change of the Stress Distribution with Cutting Width

This calculation is to show the variation of the stresses with different . The results are shown in Figure 8 where the abscissa is y coordinate with and the ordinate is the stress magnitude at different . In Figure 8, each curve represents the stress change with different coordinates between 0 and with step size at . The values have been marked on each curve. For clarity, the different stress components are drawn with different line styles.

The subplots in Figure 8 show the stresses under three different external force distribution profiles, or rectangle-shaped profile is shown in Figure 8(a), or a trapezoid shaped profile is in Figure 8(b) and or a triangle-shaped profile is in Figure 8(c). It is observed that all the stress components, except , have the maximum magnitude when . They then decrease as increases. The maximum happens inside the material at . Since is very small (in the magnitude of 10−9 m), the location of maximum is very close to the surface. The maximum values of , and are in the same magnitude of or , respectively.

##### 3.3. Fracture Force via Blade Shape

The knife shape can be defined using the parameters , , and . If and increase, the maximum pressure will decrease when the total external force keeps unchanged. The influence of the edge shape can be roughly expressed using w (please refer to Figure 5). Using the expression in (6), the stress distributions for the different values of the on the blade-material contact plane are obtained and shown in Figure 9. From Figure 9, it is observed that by keeping the external force unchanged, when the force distribution changes from constant intensity to linear intensity, the maximum magnitude of all the stresses increases. This leads to an increase in the magnitude of the maximum shear stress. Thus, by assuming is constant, the external force to realize cutting fracture will decrease. Or in order to realize cutting on the same material with the same constraints, the blade with the linear edge shape will use less force than any other blades. It will be the sharpest one.

##### 3.4. Fracture Force via Slicing Angle

In this section, the knife shape and the external force intensity profile are fixed with . The obtained stress distribution is shown in Figure 8(b). It is observed that there are four possible locations with maximum stress, and they are summarized in Table 3.

Table 3: Locations and values of the possible maximum stresses.

The changes of the corresponding stresses distribution, the principle stresses and the maximum shear stress, with slicing angle, are shown in Figures 6(a), 6(b), and 6(c), respectively. In Figure 6, only the stresses at location (ii) () = () are illustrated. Similar results at locations (i), (iii), and (iv) can be obtained accordingly. The results in Figure 6(a) are obtained using (6) by assuming , that is, = 1 when α = 0° and = 1 when α = 90°. The results in Figures 6(b) and 6(c) are obtained using (7) and (8), respectively. Using Tresca’s fracture criterion [26], if is generated by , according to (9) by setting , the cutting force is obtained using (12) as shown in Figure 6(d). Using the same method, the maximum shear stresses and the corresponding external cutting forces at locations (i), (iii), and (iv) can also be obtained. The results are shown in Figure 10(a).

Figure 6: Stresses and external forces change with slicing angle at location (ii).
Figure 7: Stress tensor changes along for different .
Figure 8: Comparison of the normalized tensor for different cutting force profiles.
Figure 9: Stresses change with on the blade-material contact plane with same cutting force.
Figure 10: Maximum shear stresses and required minimum cutting force.

Since the external force is obtained using the same ultimate shear stress , the smallest one among all the at the four locations from (i) to (iv) will initiate cutting fracture. For clarity, the largest maximum shear stress and its required smallest external force are redrawn in Figure 10(b). It can be observed from Figure 10(b) that when is from 0° to 10°, the shear stress at location (iv) will initialize fracture, and when is from 10° to 90°, the shear stress at location (ii) will initialize fracture.

The total cutting force and its tangential and normal components are shown in Figure 11. From Figure 11, it can be seen that the fracture at slicing angle from 0° to 10° is due to which is just beneath the surface ( or roughly  nm). By considering the stress direction, it can be seen that it is the mode II fracture, edge-sliding fracture [1]. The fracture from 10° to 90° is due to and which is just on the surface (). It is the mode III fracture, out-of-plane tearing [1]. From the external force profile shown in Figure 11, the influence of the slicing angle can be clearly seen. The required force to cut by pressing-only is far larger than the force required by pressing-and-slicing cuts. Note also is that force does not change a lot from 10° to 90° in the mode III fracture and its value is just the total force when the slicing angle is 90°.

Figure 11: Total cutting force and the force components in direction () and direction ().

#### 4. Experimental Verification

Raw Russet Baking Potatoes, purchased from a local Kroger Store in Atlanta Ga, USA, were selected as the testing materials in the experiments. The potatoes were purchased in fresh. They are firm and smooth without dark spots, green areas, mold, or cuts. They are in the size of 125 mm–150 mm long and about 75 mm in diameter. The moisture content is about 4 g water per 5 g fresh potato solid. Moisture content of the sample was determined by drying thin potato slabs (10 mm by 20 mm by 45 mm) at 70°C under vacuum for 24 hours. The experiments were carried out in an air conditioning controlled room at temperature about 22°C and moisture about 70%. Once a fresh potato was peeled and cut to a roughly 20 mm by 40 mm × 100 mm brick, the experiments were performed without delay. The experimental setup is shown in Figure 12. The system consists of an ABB robot IRB 140 [27] for motion generation, speed control and distance measurement, and an ATI Force/Torque sensor ISA F/T-16 Mini40 [28] for force measurement. The robot is commanded to move at speed of 0.5 mm/s for cutting. The force data are saved on a central computer. The forces during the cut of potatoes are shown in Figure 13.

Figure 12: Experiment setup.
Figure 13: Cutting force changes with contact length.
##### 4.1. Evaluating the Influence of Slicing Angle and Cutting Fracture Modes

When changing the slicing angle, the same shapes of the cutting force changing profiles have been obtained. The average data with standard deviation less than 0.2 lb from 5 sets of cuts using different slicing angles on potatoes are summarized in Figure 14. According to the maximum external force when is near 90°, where the tear mode dominates, the ultimate shear stress () can be estimated as  lb. Then, using this and the procedure to obtain and in Figure 11, the theoretical fracture forces can be estimated using (12) by assuming the same , , , , and since only the slicing angle () changed during different trails. Both of the experimental data and the theoretically estimated data are shown in Figure 14. Good match is observed.

Figure 14: Experimental data and theoretical estimation of the cutting forces during potato cuts.
##### 4.2. Evaluating Relative Blade Sharpness Factor

These experiments were used to demonstrate how to implement the sharpness definition in this paper. In the experiments, blades with different sharpness were used to cut the same material to study the effect of blade sharpness. The materials (potatoes) were in rectangular shape with the same width (13 mm) and thickness (30 mm) to make sure the same contact length of each cut. The blades were prepared such that the sharpest blade and the dullest blade defined in (15)–(17) were included.

Table 4: Blade relative sharpness level.
Figure 15: Cutting results using different sharp blades (only blades A, D, E are shown).

By using this method, any blade can be assigned to a certain sharpness level. When the cutting force reaches certain level by which the knife is decided to be dull, the knife needs to be resharpened or disposed.

#### 5. Justification

It can be seen that the computation results from the model agree well with the potato cutting results. Noted is that the results are from the cuts of hard-crispy material. The modeling method should applicable.

Then the question is how about large deformation nonlinear materials. The mathematical model is clearly not applicable to nonlinear materials, such as meat. However, our everyday experience told us that the cuts on meat follow the rule: slicing is easier that pressing. Here, the author just provides a weak explanation of the cuts on meat: Type II fracture happens at the depth () about 0.37 times of half of the blade width or about = 185 nm if = 1000 nm. It is almost at the contact surface when you compare the whole thickness of the cutting material. Type III fracture happens at the contact surface, that is, = 0. In this small distance or on the contact surface, the inhomogeneity, inelasticity, and anisotropy of the cutting material will not affect the stress distribution since the react force/pressure will balance the external force statically based on Newton’s third law of motion. From Figure 8, it has been checked that the stress distribution in these locations is balanced by the external force. Thus, at the moment of cutting initialization, although the model will not be valid at bigger for big deformation nonlinear materials, the model would be still valid at and . Since we only care about the stresses to generate the possible fracture (the results in Figure 7 are just to visualize the stress distribution yielded by the model), the stress when is not within our consideration and does not relate to our conclusion. Thus, the model is applicable to large deformation nonlinear materials only when . We further constrain our application range for hard-crispy materials within , that is, not all .

#### 6. Conclusions

In order to understand the question: why pressing-and-slicing cuts use less force than pressing-only cuts during food cutting, such as potatoes, this paper formulated and studied the change of the stress distribution with various influence factors. The cutting interaction was modeled as a belt-area force acting on the surface of a semi-infinite body. The cutting force was assumed to be in certain profile based on the observed shape of blade cutting edges. The closed form expression of the cutting interaction has been developed using the direct integration method. Compared to the knowledge in the current literature, the improvements of the modeling are twofolds: (1) this work originally expressed the cutting force (force at cutting fracture) and its change with slicing angle using mathematical equations. (2) for the first time, this work mathematically expressed the blade sharpness with blade parameters and used relative sharpness concepts based on mathematical descriptions (not based on experience) to describe blade sharpness.

Regarding the future work, the stress intensity factors and will be analyzed and quantified. The influence of the relative moving speed will also be investigated. On the long term, the research work will try to understand and model the behavior of interesting biomaterials or hybrid materials during robot-controlled cuts. Then robots can be intelligent to adapt themselves for any deformation of material and any varieties of material’s structure.

#### A. Derivation of Stress Distribution dut to Various Forces

##### A.1. Stress Distribution due to Point Force

(1) Normal Point Force Only
When there is only normal point force , for material with Poisson’s ratio , the normal stress and the shear stress at point are given as where The superscript (stands for normal direction and on a point) in (A.1a) to (A.1g) represents the stresses generated by external force .

(2) Tangential Point Force Only
When there is only tangential point force , the stress distribution is known as the Cerruti solution and is given in [8] as where The superscript in (A.2a) to (A.2g) represents the stresses generated by external force .

##### A.2. Stress Distribution due to Line Force

(1) Normal Line Force Only
The closed form solution of the stresses at any point () in the semi-infinite body generated by line force acting along () are shown in (A.3).When and , the stress tensor is

(2) Tangential Line Force Only
The closed-form expression of the stress distribution by a tangential force on line segment [] is given in (A.5). where
As and go to and , respectively, that is, a line force acting on the boundary of a semi-infinite body, the stress distribution changes to

##### A.3. Stress Distribution due to Area Force

(1) Normal Area Force Only
The stresses generated by the normal force intensity between [] are where , and .
The stresses generated by the normal force intensity between () are obtained as
The stresses generated by the normal force intensity between [] are obtained as
The stress tensor generated by normal force can then be expressed as

(2) Tangential Area Force Only
The stresses generated by the tangential force intensity between [] are obtained as
The stresses generated by the tangential force intensity between () are obtained as
The stresses generated by the tangential force intensity between [] are obtained as
The final stress tensor generated by tangential force can then be expressed as

(3) Force with Both Normal and Tangential Components
Using (A.11) and (A.15), the stress distribution generated by both and is obtained as Equation (A.16) can be rewritten as where the functions , , , , , and can be explicitly expressed using (A.8) to (A.10) and (A.12) to (A.14) using the variables .

#### Symbols

 𝑎: Half of the width of the blade edge (see Figure 3) 𝐴(𝑥,𝑦,𝑧): A point in the material with coordinate (𝑥,𝑦,𝑧) in frame 𝑜𝑥𝑦𝑧 (see Figure 3) 𝐾𝑠: Ratio of the fracture initialization force over continuous cutting force 𝑙: Contact length between the blade and the material (see Figure 3) 𝑝𝑛 and 𝑝𝑡: Cutting pressure components in the 𝑧 (normal) direction and 𝑥 (tangential) direction, respectively (see Table 2) 𝑝𝑙𝑛 and 𝑝𝑙𝑡: Normal and tangential components of line distributed cutting force 𝑃: Total cutting force (see Figure 4 and Table 2) 𝑃𝑐: Force during continuous cutting fracture 𝑃𝑢: Required force to initialize cutting fracture 𝑃𝑛 and 𝑃𝑡: Normal and tangential component of total applied external force, respectively 𝑃𝑛𝑢 and 𝑃𝑡𝑢: Normal and tangential component of total cutting force during fracture, respectively 𝑃𝑝𝑛 and 𝑃𝑝𝑡: Normal and tangential component of a point applied external force, respectively 𝑞𝑛 and 𝑞𝑡: Maximum of the area force intensity between [−𝑤,𝑤] in the normal and tangential direction, respectively 𝑟: Distance between point 𝑜 and point 𝐵 (see Figure 3(b)) 𝑅: Distance between point 𝐴 and point 𝐵 (see Figure 3(b)) 𝑤: Width of the top side of the trapezoid profile of force intensity distribution (see Figure 5) 𝛼: Slicing angel, that is, the angle between 𝑃 and 𝑃𝑛 (see Figures 3(b) and 2) 𝜎 and 𝜏: 𝜎 is the tensile stress or just an expression of any of the stress components and 𝜏 is only the shear stress 𝜏𝑖𝑗: Maximum shear stress, 𝑖 or 𝑗=1, 2, or 3 𝜏max: Maximum of the maximum shear stresses 𝜏𝑖𝑗, 𝑖 or 𝑗=1, 2, or 3 𝜏𝑢: Ultimate shear stress 𝜇: Poisson’s ratio 𝜂: Blade relative sharpness factor 𝜂0: Relative sharpness factor of the sharpest blade 𝜂𝑓: Relative sharpness factor of the dullest blade 𝜅: Blade relative sharpness level.

#### References

1. K. Hellan, Introduction to Fracture Mechanics, McGraw-Hill, New York, NY, USA, 1984.
2. M. Mahvash and Hayward V., “Haptic rendering of cutting: a fracture mechanics approach,” Haptics-e, vol. 2, no. 3, pp. 1–12, 2001.
3. A. G. Atkins, X. Xu, and G. Jeronimidis, “Cutting, by 'pressing and slicing,' of thin floppy slices of materials illustrated by experiments on cheddar cheese and salami,” Journal of Materials Science, vol. 39, no. 8, pp. 2761–2766, 2004.
4. A. Atkins and Y. Mai, Elastic and Plastic Fracture: Metals, Polymers, Ceramics, Composites, Biological Materials, Ellis Horwood, Chichester, UK, 1985.
5. I. Kamyab, S. Chakrabarti, and J. G. Williams, “Cutting cheese with wire,” Journal of Materials Science, vol. 33, no. 11, pp. 2763–2770, 1998.
6. S. Way, “Some observations on the theory of contact pressures,” Journal of Applied Mechanics, vol. 7, pp. 147–157, 1940.
7. M. Sadd, Elasticity: Theory, Applications, and Numerics, Elsevier, Burlington, Mass, USA, 2005.
8. Z. Xu, Elasticity Mechanics, Higher-level Education, 5th edition, 1996.
9. A. E. H. Love, “The stress produced in a semi-infinite solid by pressure on part of the boundary,” Philosophical Transactions of the Royal Society A, vol. 667, pp. 377–420, 1929.
10. J. Boussinesq, Application des Potentiels a I’Etude de l’Euqilibre et du Mouvement des Solides Elastiques, Gauthier-Villars, Paris, France, 1885.
11. F. F. Ling, M. Lai, and D. Lucc, Fundamentals of Surface Mechanics: with Applications, John Wiley Sons, Nwe York, NY, USA, 1973.
12. S. Timoshenko, Theory of Elasticity, McGraw-Hill, New York, NY, USA, 3rd edition, 1970.
13. K. L. Johnson, Contact Mechanics, Cambridge University, New York, NY, USA, 1985.
14. A. E. H. Love, A Treatise on the Mathematical Theory of Elasticity, Cambridge University, New York, NY, USA, 4th edition, 1927.
15. J. R. Dydo and H. R. Busby, “Elasticity solutions for constant and linearly varying loads applied to a rectangular surface patch on the elastic half-space,” Journal of Elasticity, vol. 38, no. 2, pp. 153–163, 1995.
16. O. J. Svec and G. M.L. Gladwell, “An explicit boussinesq solution for a polynomial distribution of pressure over a triangular region,” Journal of Elasticity, vol. 1, no. 2, pp. 167–170, 1971.
17. J. Li and E. J. Berger, “A Boussinesq-Cerruti solution set for constant and linear distribution of normal and tangential load over a triangular area,” Journal of Elasticity, vol. 63, no. 2, pp. 137–151, 2001.
18. G. M. Hamilton and L. E. Goodman, “The Stress field created by a circular sliding contact,” Journal of Applied Mechanics, vol. 33, pp. 371–376, 1966.
19. A. C. Fischer-Cripps, Introduction to Contact Mechanics, Springer, New York, NY, USA, 2000.
20. H. Yoshihara and A. Matsumoto, “Measurement of the shearing properties of wood by in-plane shear test using a thin specimen,” Journal of International Academy of Wood Science, vol. 39, no. 2, pp. 141–152, 2005.
21. P. Lucas, Dental Functional Morphology: How Teeth Work, Cambridge University Press, New York, NY, USA, 1st edition, 2007.
22. Tech Edge 2005, http://www.furitechnics.com.au/information/Tech%20Edge/Cutting%20edge%20ang les.htm.
23. R. W. McGorry, P. C. Dowd, and P. G. Dempsey, “The effect of blade finish and blade edge angle on forces used in meat cutting operations,” Applied Ergonomics, vol. 36, no. 1, pp. 71–77, 2005.
24. R. L. Szabo, R. G. Radwin, and C. J. Henderson, “The influence of knife dullness on poultry processing operator exertions and the effectiveness of periodic knife steeling,” American Industrial Hygiene Association Journal, vol. 62, no. 4, pp. 428–433, 2001.
25. D. Zhou and G. McMurray, “Modeling of blade sharpness and compression cut of biomaterials,” Robotica, vol. 28, no. 2, pp. 311–319, 2010.
26. C. Lipson and R. Juvinall, Handbook of Stress and Strength: Design and Material Applications, Macmillan, New York, NY, USA, 1963.
27. ABB Robotics, “Product Manual IRB 140,” ABB Robotics Products AB publication, no. M2000, article 3HAC 7564-1, 2000.
28. ATI Industry Automation, “Installation and operations manual for ISA F/T-16,” ATI Industry Automation Inc., 1998http://www.ati-ia.com/app_content/documents/9610-05-1019%20ISA%20FT-16%20%28ALL%29.pdf.