MATHEMATICAL SIMULATION OF KNIFE PROFILE RESISTANCE FORCE DURING FISH CUTTING
Abstract and keywords
Abstract (English):
The article focuses on studying the fish cutting process and modeling forces of harmful resistance. The fish muscular tissue rheological properties are described by a Maxwell-Thomson model. The conditions of constrained compression of the material across the width and the absence of constrained compression in the direction of movement of the knife are accepted. On the basis of the energy approach, the profile resistance force of the double-edged knife has been interpreted as deformational force of the friction at the macroscopic scale level, provided that the surface of the faces is smooth. The mathematical models for dimensional and dimensionless profile resistance forces of the knife without side edges have been developed. The dependence of the dimensional force on the sharpening angles, knife thickness, rheological properties and cutting speed has been established. The dependence of the dimensionless force on the dimensionless cutting speed and measure of the muscle tissue elasticity has been shown. The profile resistance forces of flat-back knife and double-edged knives have been analyzed. With sharpening angle of back edges = 5°; 10°; 20°; 50°, force maximums are 0.317; 0.306; 0.288; 0.274, respectively. When the values of instantaneous modulus of elasticity 1.5·105; 2.0·105; 2.5·105; 3.0·105 N/m2, the maximums of the specified force are 0.310; 0.411; 0.513; 0.614 N, respectively. With the values of elasticity = 4; 7; 11; 15, dimensionless force maximums of flat-back knife are 1.959; 3.166; 4.774; 6.381 and without side edges - 1.193; 1.864; 2.764; 3.663, respectively.

Keywords:
fish, cutting, force, resistance, profile, knife, edge, rheology, viscoelasticity
Text
Publication text (PDF): Read Download

Ensuring resource saving when cutting fish provides for a thorough analysis of the resistance forces. Reduction of harmful resistance force involves the reduction of its components - the profile resistance forces and the friction force. Mathematical modeling of the food cutting process is a relevant scientific direction and attracts the attention of researchers. An experimental relationship between the resistance force of the fibrous material destruction and the motion speed of the knife has been established [1]. The effect of the blade sharpening angle on the parameters of this process was analyzed [2]. The regularities of high-speed processing of viscoelastic materials in a wide range of speeds: from 0.001 m /s to 10 m /s, have been studied [3, 4]. Finite element modeling of the viscoelastic product destruction has been performed [5]. An approach to determine the optimal geometry of a knife by minimizing the normal component of the resulting resistance (maximizing “sliding” cutting) has been developed [6, 7]. However, currently there is no analytical description of the harmful resistance forces acting on the working tool when cutting fish. At the same time, to optimize the geometry of the knife according by the criterion of minimum resistance, mathematical modeling of the forces acting on its faces is required. Material The choice of rheological models of fish muscle tissue has been substantiated. Differential equations of models with their solutions for three different loading conditions have been considered. The results of experimental tests of muscular tissue of scab, mackerel, Atlantic sardinella for direct creep, relaxation and reverse creep are presented. It has been established that the muscle tissue of the fish before destruction shows a limited flow under load, relaxes under constant load to an equilibrium state, and is fully restored at full unloading. It is shown that the results of the experimental tests approximately correspond to the three-element Maxwell - Thomson rheological model. Methods The purpose of mathematical modeling is to determine the profile resistance force that occurs when cutting fish with a two-edged knife without side edges. According to modern concepts of tribology, the profile resistance force is interpreted as the deformation friction force at the macroscopic scale level [8]; Popov, 2017 [9]. At this level, it is necessary to solve the problem of determining the sliding friction force of the stamp (knife as a completely solid smooth body) on the deformable viscoelastic material (cut muscle tissue of the fish). In this case, the roughness of the knife surface is not taken into account, since the macrogeometry of the profile of the cutting working tool is the determining factor. The task of finding the deformation force of friction when sliding a knife on a viscoelastic material at the macroscopic level (the profile resistance force) is advisable to solve on the basis of an energy approach [10, 11]. This approach is more general and is based on the determination of energy losses due to the viscoelastic deformation of the material. The energy approach is valid for a stamp of arbitrary shape, which allows it to be applied to knives with different profiles. Mathematical simulation According to the energy approach, when the stamp slides, the energy dissipation in a viscoelastic material occurs due to viscous friction in the Kelvin-Voigt element damper of the Maxwell-Thomson model. The deformation friction force of the stamp when sliding on a viscoelastic base is determined as follows: , (1) where p(y) - the normal contact pressure on the stamp from the material; g(y) - the function that determines the stamp geometry (the shape of the knife profile); g׳(y) - derivative of the specified function; a, b - coordinates along the axis (0, y) of the extreme points of the arbitrary shape stamp. The function g(y) determines the difference between the coordinates of the point of the stamp profile with the coordinate and the profile point with zero coordinate, so that g(0) = 0. Consider the movement in the fish muscle tissue of an elementary knife, shown in Fig. 1. The center of the moving Cartesian coordinate system Oxy is at a point. Fig. 1. The scheme of movement of double-edged knife without side edges during fish cutting In a quasistatic setting, we believe that at any time the knife is in a state of equilibrium. When cutting, the blade performs a rectilinear uniform motion in the material. The material is in conditions of constrained compression along the axis Ox, and the cutting is carried out without compressing the material in the direction of knife movement along the axis Oy (the height of the material sample is much larger than its width). In Fig. 1 are marked: - the contact pressure of the horizontal elementary fiber dy at the point A, normal to the direction of knife movement; - the counter material resistance to the movement of the point A, due to the deformation of the specified fiber along the axis Ox; - the contact pressure of the material on the front inclined face at a point A; - the normal contact pressure of the horizontal elementary fiber dy at the point I of the back inclined edge of the knife; - associated contact pressure at the point I of the back inclined edge of the knife; - the contact pressure of the material on the back inclined edge at a point I. We introduce the following notation: E0 - instantaneous modulus of elasticity of fish muscle tissue (Young's modulus); E1 - high-elasticity modulus of fish muscle tissue; η - the coefficient of dynamic viscosity of the fish muscle tissue; e01 = E0/E1 - a measure of fish muscle tissue elasticity; lk - blade edge length; δ - half the thickness of the knife; α - half the angle of sharpening the front edge of the knife; β - half the angle of sharpening the back edge of the knife; +; hm - height of the inclined front edge; l - half the width of the material in the equilibrium state; ν - knife speed; (quasistatic modulus of elasticity); ; ; ; - dimensional profile resistance force; - dimensionless profile resistance force. The right half of a two-edging (double-edged) knife (Fig. 1) has two extreme contact points: a point O(0, 0) and a point G , as well as one corner point B. Thus, in expression (1) a = 0, . In this regard, the profile resistance force of a two-edging knife according to (1) is determined by the following expression: . (2) Let us write the boundary conditions of the knife contact with the material in the following form p(0) = 0; . Let us represent the integral in expression (2) as the sum of the integrals. Then, taking into account the boundary conditions, the expression (2) has the following form: (3) The profile shape of the two-edging knife without side edges determines the type of functions g(y) and : Let us substitute into expression (3) the expressions for normal contact pressures p2(y) and p6(y) at the front and back inclined faces and the coordinates of the extreme contact point G, obtained by successively solving the first order differential equation of the Maxwell - Thomson model in regions OB and BG [12, 13]. We use the formula of Newton - Leibniz and get the expression for the profile resistance force of a two-edging knife without side edges: . Let us introduce the dimensionless knife speed: . Then taking into account ; ; ; (dimensionless coordinate of a point G); the expression for the dimensionless profile resistance force of a two-edging knife has the following form: Results Fig. 2-6 show the simulation results of the profile resistance force of a two-edging knife without side edges. The force F1 corresponds to a knife with a straight butt, - a knife with side edges, - a knife without side edges. 1 2 3 4 4 3 2 1 а b Fig. 2. Dependence of the profile resistance force on the cutting speed at different rheological properties of the material (α = 5°; β = 15°; δ = 3 mm; l = 50 mm; L = hm = 12 mm): a - at low speeds; b - at high speeds; 1 - E0 = 1.5 ∙ 105 N/m2, E1 = 0.1 ∙ 105 N/m2; 2 - E0 = 2 ∙ 105 N/m2, E1 = 0.3 ∙ 105 N/m2; 3 - E0 = 2.5 ∙ 105 N/m2, E1 = 0.6 ∙ 105 N/m2; 4 - E0 = 3 ∙ 105 N/m2, E1 =0.8 ∙ 105 N/m2 4 3 2 1 4 3 2 1 а b Fig. 3. Dependence of the profile resistance force on the cutting speed at different sharpness half-angle β of the back inclined edge: a - at low speeds; b - at high speeds; E0 = 3 ∙ 105 N/m2; e01 = 3.75; δ = 3 mm: 1 - β = 5°; 2 - β = 10°; 3 - β = 20°; 4 - β = 50° 1 2 3 4 1 2 3 4 a b Fig. 4. Dependencies of the dimensionless profile resistance force of double-edged knife with side edges (a) and knife without side edges (b) on the dimensionless cutting speed at different values of elasticity measure (): 1 - e01 = 4; 2 - e01 = 7; 3 - e01 = 11; 4 - e01 = 15 3 2 1 3 2 1 a b Fig. 5. Dependencies of the dimension (a) and dimensionless (b) profile resistance force on the dimensionless cutting speed: a - E0 = 3 ∙ 105 N/m2, E1 = 0.8 ∙ 105 N/m2; η = 1.5 ∙ 107 N∙s/m2; α = 5°; δ = 3 mm; 1 - F1; 2 - ; 3 - ; b - e01 = 3.75; tαβ = 0.8; 1 - ; 2 - ; 3 - 3 2 1 3 2 1 a b Fig. 6. Dimension (a) and dimensionless (b) profile resistance force: a - dependence of dimension forces on coefficient of dynamic viscosity of the material (ν = 1 mm/s): 1 - F1; 2 - ; 3 - ; b - dependence of dimensionless forces on elasticity measure (; ; ): 1 - ; 2 - ; 3 - Discussion of results To reduce energy costs for cutting it is necessary to strive to reduce the resistance forces. In this regard, let us analyze the influence of various parameters on the dimensional and dimensionless profile resistance forces of a two-edged (double-edged) knife without side edges (Fig. 1). In this case, the process of reverse creep of the material proceeds without intermediate stress relaxation with constant deformation of the muscle tissue. Fig. 2 shows that with a monotonous increase in the elastic properties of the material, the maximum force increases significantly and is achieved at very low cutting speeds. With a further increase in cutting speed with increased stiffness of muscle tissue, the profile resistance force of the knife decreases. Fig. 3 shows that in the absence of the knife side faces, the influence on the indicated force of the sharpening angle of the back inclined edges is essential. With an increase in the specified angle of sharpening at low cutting speeds, the force maximum decreases, and at high speeds, the force also decreases. At angles of sharpening 5°; 10°; 20°; 50° the maximums of the considered force are 0.317; 0.306; 0.288; 0.274 N, respectively. Fig. 4 shows that the absence of stress relaxation with a constant deformation of the material by the side faces leads to a significant decrease in the dimensionless profile resistance force, ceteris paribus. With the values of the elasticity measure 4; 7; 11; 15 the maximums of the dimensionless force profile resistance force of the knife with the side faces are 1.959; 3.166; 4.774; 6.381, without side faces - 1.193; 1.864; 2.764; 3.663, respectively. In general, Fig. 2-4 show that as the cutting speed increases, the force tends asymptotically to zero. Fig. 5, 6 illustrate that the greatest profile resistance force occurs in the case of cutting fish with a straight butt knife, and the smallest force corresponds to a two-edging knife without side edges, which, with accepted assumptions, is the most effective tool. The dependences of dimensional and dimensionless forces on the elastic modulus, elasticity measures and coefficient of dynamic viscosity of a material are monotonic. When the values of the instantaneous modulus of elasticity 3·105 N/m2; lagging modulus of elasticity 0.8·105 N/m2; coefficient of dynamic viscosity 1.5·107 N∙s/m2; the half angle of sharpening the front inclined edges and the half thickness of the knife is 3 mm; then the limit of the profile resistance force of the knife with a straight butt is 1.345 N; the maxima of the resistance forces of the knife with and without side edges are 1.025 N and 0.646 N, respectively. When the measure of material elasticity is 3.75 and the ratio of tangent angles of sharpening is 0.8, the limit of the dimensionless profile resistance force of the knife with a straight butt is 2.367, the maxima of the dimensionless profile resistance force of knives with and without side edges are 1.804 and 1.136, respectively. With increasing material viscosity, the profile resistance force of the knife with a straight butt increases and asymptotically tends to the limit value, while the profile resistance forces of the two-edged knives decrease and asymptotically tend to zero. At the same time from Fig. 4-6, we see that the minimum dimensional force occures when cutting fish with a knife without side edges. With an increase in the measure of material elasticity, the smallest dimensionless force also corresponds to the specified two-edged knife. Conclusion 1. The profile resistance force of the knife essentially depends on the rheological properties of the fish - the elastic modulus and the dynamic viscosity coefficient. With an increase in the rigidity of the raw material, this force is significantly reduced at medium and high cutting speeds. 2. The dependence of the profile resistance force of the knife on the cutting speed is non-monotonic with a pronounced maximum. With increasing speed, this force asymptotically tends to a zero value, which distinguishes it from the resistance force of a single-edged knife, which, under the same conditions, monotonously increases and asymptotically tends to its limit value. 3. The knife geometry affects the profile resistance force of the knife. Half-angles of sharpening the front and back inclined edges affect the dimensional force, which determines the expediency of setting and solving the problem of optimizing the knife profile by the criterion of the minimum force of harmful resistance. The half thickness of the knife affects the instant-elastic and high-elastic components of the resistance force. 4. The absence of side edges leads to the absence of viscoelastic energy losses during stress relaxation, which leads to a significant reduction in the profile resistance force. This allows us to conclude that for knives with straight edges, ceteris paribus, the lowest energy costs for cutting are provided when cutting with a thin two-edging knife without side edges. This is confirmed by current trends in the development of working tools of filleting machines.
References

1. 1. Dowgiallo A. Cutting force of fibrous materials. Journal of Food Engineering, 2005, no. 66, pp. 57-61.

2. 2. Boisly M., Schuldt S., Kaestner M. G., Schneider Y., Rohm H. Experimental characterisation and numerical modelling of cutting processes in viscoelastic solids. Journal of Food Engineering, 2016, no. 191, pp. 1-9.

3. 3. Schuldt S., Arnold G., Kowalewski J., Schneider Y., Rohm H. Analysis of the sharpness of blades for food cutting. Journal of Food Engineering, 2016, no. 188, pp. 13-20.

4. 4. Schuldt S., Schneider Y., Rohm H. High-speed cutting of foods: Cutting behavior and initial cutting forces. Journal of Food Engineering, 2018, no. 230, pp. 55-62.

5. 5. Pagani M., Perego U. Explicit dynamics simulation of blade cutting of thin elastoplastic shells using “directional” cohesive elements in solid-shell finite element models. Computer methods in applied Mechanics and Engineering, 2015, no. 285, pp. 515-541.

6. 6. Atkins T. Optimum blade configurations for the cutting of soft solids. Engineering Fracture Mechanics, 2006, no. 73, pp. 2523-2531.

7. 7. Atkins T. Prediction of sticking and sliding lengths on the rake faces of tools using cutting forces. International Journal of Mechanical Sciences, 2015, no. 91, pp. 33-45.

8. 8. Goryacheva I. G. Mekhanika frikcionnogo vzaimodejstviya [Mechanics of friction interaction]. Moscow, Nauka Publ., 2001. 478 p.

9. 9. Popov V. L. Contact Mechanics and Friction: Physical Principles and Applications. Berlin, Springer Verlag GmbH, 2017. 391 p.

10. 10. Soldatenkov I. A. K raschetu deformacionnoj sostavlyayushchej sily treniya dlya standartnogo vyazkouprugogo osnovaniya [On calculating deformation component of friction force for a standard viscoelastic base]. Trenie i iznos, 2008, no. 1, vol. 29, pp. 12-21.

11. 11. Soldatenkov I. A. Raschet treniya indentora s fraktal'noj sherohovatost'yu o vyazkouprugoe osnovanie [Calculation of friction indenter with fractal roughness on a viscoelastic base]. Trenie i iznos, 2015, no. 3, vol. 36, pp. 257-262.

12. 12. Ageev O. V., Fatyhov Yu. A., Samojlova N. V. Vybor i identifikaciya reologicheskoj modeli strukturno-mekhanicheskih svojstv myshechnoj tkani ryby [Selection and identification of rheological model of structural-mechanical properties of fish muscle tissue]. Izvestiya Kaliningradskogo gosudarstvennogo tekhnicheskogo universiteta, 2018, no. 49, pp. 75-91.

13. 13. Ageev O. V., Naumov V. A., Fatyhov Yu. A., Samojlov N. V. Analiz sootvetstviya reologicheskih modelej strukturno-mekhanicheskim svojstvam ryby [Analysis of rheological models compliance with structural and mechanical properties of fish]. Nauchnyj zhurnal Sankt-Peterburgskogo nacional'nogo issledovatel'skogo universiteta informacionnyh tekhnologij, mekhaniki i optiki. Seriya: Processy i apparaty pishchevyh proizvodstv, 2018, no. 2 (36), pp. 34-43. DOI: 10.17586/2310-1164-2018-11-2-34-43.


Login or Create
* Forgot password?