Discrete element thermomechanical modelling of rock cutting with valuation of tool wear

The paper presents a thermomechanical discrete element model of rock cutting process. The thermomechanical formulation of the discrete element method considers mechanical and thermal phenomena and their reciprocal influence. The thermal model developed for transient heat conduction problems takes into account conductive heat transfer at the contact between particles and convection on the free surface. The thermal and mechanical problems are coupled by consideration of: (1) heat generated due to friction which is calculated in the mechanical problem and passed to the thermal solution, (2) influence of thermal expansion on mechanical interaction between particles. Estimation of temperature dependent wear has been included into the contact model. The coupled problem is solved using the staggered scheme.The thermomechanical algorithm has been implemented in a discrete element program and applied to simulation of rock cutting with single pick of a dredge cutter head. Numerical results confirm good performance of the developed algorithm.


Introduction
A variety of rock-cutting technologies is used in civil as well as in mining engineering. Figure 1 shows a cutter head of a dredging machine performing rock cutting in underwater conditions. Rock cutting consists in fracturing and disinte-J. Rojek  Performance and efficiency of cutting processes are evaluated using different criteria such as cutting forces and their variation, cutting rate, specific energy and tool wear [2,3,9,17]. The wear leads to significant changes of the shape of rock cutting tools (Fig. 3) and it is an important factor affecting the overall efficiency of excavation works, firstly, because of its adverse effect on cutting forces, excavation rate and energy consumption, and secondly, due to the costs and delays incurred by replacements of the worn-out tools [34].
Prediction of wear is an important element in the design and optimization cutting tools.Tool wear depends on both mechanical and thermal effects in the tool-rock interaction. High temperatures, developed on the tool surface (in some cases temperature on the tool surface can reach as much as 550-800 • C, cf. [20,35]) due to heat generated by friction between tool and rock, reduce the hardness of the surface and increase adhesive and abrasive wear.
Various theoretical models have been developed for rock cutting. Simple analytical models have been proposed by Evans [9] and Nishimatsu [23]. The complex stress state and crack propagation patterns involved in the process of rock cutting, however, makes it difficult for theoretical analysis. The rock failure process during cutting can be traced in more detail by appropriate numerical methods.
A number of numerical studies of rock cutting utilized the finite element method, cf. [16,18,19,30,39]. However, the finite element method based on the continuum mechanics theory of material modelling has serious problems in representing properly discontinuities of the mate-  rial occurring during rock cutting [16]. Special formulations are necessary to introduce the possibility of discontinuum analysis of rock fracture [30,39]. The discrete element method can take into account most kinds of discontinuities and material failure characterized with multiple fracture, which makes it a suitable tool to study rock cutting [14,15,26,28,31,32].
In order to consider the temperature dependent wear of rock cutting tools it is necessary to include thermal effects into the model. Thermal effects, including frictional heating and convective cooling of drag tools during rock cutting have been modelled numerically using 2D thermal finite difference models by Ortega and Glowka [24,25]. Temperature distribution in the tool has been obtained for assumed frictional heating and convective cooling rates. The results have been compared with experimental measurements showing quite a good performance of the simple numerical model. A two-dimensional transient heat transfer model was developed using the finite element method for the study of temperature rise during continuous drag-cutting by Loui and Karanam [20]. The effect of frictional force and cutting speed on the temperature developed at the pickrock interface was compared with the experimental data.
The heat conduction problem can also be analysed using the discrete element method. The discrete element model originally introduced to model mechanical problems [7] can be extended on thermal and thermomechanical problems. The discrete thermal element method (DTEM) proposed by Feng et al. [11] considered heat conduction within the particle using analytical integral solution for the temperature distribution over a circular domain subjected to the Neumann boundary condition. This concept introduced significant complexity into a simple formulation of the discrete element method. A simplified version of the DTEM, termed the pipe-network model, in which each particle is replaced by a simple thermal pipe-network connecting the particle centre with each contact zone associated with the particle has been developed by the same authors [12]. The lattice model of thermal pipes (contact lines) linking heat reservoirs at particle centres is commonly adopted as a thermal model in the discrete element model, cf. [13,22,40]. A discrete element thermal analysis can be performed independently [12,13,40,41] or coupled with a mechanical problem solution [8,22,33,36] as well as with other physical fields [42,43].
The present work extends The discrete element model of rock cutting presented in [26,28] is extended here on thermomechanical analysis. The present contains original elements in theoretical formulation, such as evaluation of temperature dependent wear, and in modelling approach to rock cutting process which has not been analysed before with the discrete element method as the coupled thermomechanical problem.

Thermomechanical model of rock cutting: basic assumptions
A rock cutting process is considered here as a coupled thermomechanical problem, in which mechanical phenomena are accompanied by correlated thermal effects. Thermomechanical model of the tool-rock system is shown schematically in Fig. 4. The main physical phenomenon analyzed in the mechanical problem is the interaction of the tool with a rock leading to a failure of the rock characterized with discontinuous material behaviour. Interaction forces between tool and rock are determined. Simultaneously heat generated by friction between the tool and rock is calculated and passed to the thermal model. The thermal model takes into account heat absorption and conduction through the tool and rock as well as convective heat transfer on the free surface of the rock and tool and heat transfer in the contact surface between the rock and tool. Changing temperature distribution is determined in the solution of the thermal problem Change of temperature cause material expansion and may produce thermal stresses. Increase of temperature influences material properties and tool wear.
In the model of rock cutting assumed in this work, the tool is considered rigid, assuming that its stiffness is sufficient to produce rock failure and its deformation is irrelevant for the purposes of modelling of rock failure. Thus no stress analysis is performed for the tool.

Thermomechanical discrete element method formulation
The numerical model of the tool-rock system allowing us to simulate a process of rock cutting as coupled thermomechanical problem has been developed within the framework of the discrete element method (DEM) employing spherical (in 3D) or cylindrical (in 2D) particles. Discrete element formulation using spherical or cylindrical particles was first proposed by Cundall and Strack [7]. The present model has been developed within the author's own implementation of DEM [26,29].

Equations of motion
The translational and rotational motion of rigid spherical or cylindrical elements (particles) is governed by the standard equations of rigid body dynamics. For the ith element we have where u i is the element centroid displacement in a fixed (inertial) coordinate frame X, ω i -the angular velocity, m i -the element mass, J i -the moment of inertia, F i -the resultant force, and T i -the resultant moment about the central axes.
The form of the rotational Eq. (2) is valid for spheres and cylinders (in 2D) and is simplified with respect to a general form for an arbitrary rigid body with the rotational inertial properties represented by a second order tensor. Vectors F i and T i are the sums of all forces and moments applied to the ith element due to external load, F ext i and T ext i , respectively, the contact interactions with neighbouring spheres, F cont i j , j = 1, . . . , n c i , where n c i are the number of elements being in contact with the ith discrete element, and the forces and moments resulting from external damping, F damp i and T damp i , respectively: where s c ij is the vector connecting the centre of mass of the ith element with the contact point with the jth element.

Mechanical contact model
The overall mechanical behaviour of the system is determined by the contact laws assumed for the particle interaction. The contact law can be seen as the material model on the microscopic level. The formulation of the contact model employs the decomposition of the contact force between two elements 1 F c into the normal and tangential components, F n and F t , respectively: where n is the unit normal vector at the contact point.
Modelling of rock or other cohesive materials requires contact models with cohesion allowing tensile interaction force between particle [4]. In the present formulation, rock materials are modelled using the elastic perfectly brittle model of contact interaction, in which initial bonding between neighbouring particles is assumed. The bonds are established between particles i and j satisfying the condition: where r 0 i and r 0 j are the particle radii at the initial temperature, g 0 is the overlap or gap at the contact point, and g 0 max is a length parameter, being a tolerance in contact verification. The cohesive bonds can be broken under excessive load which allows us to simulate initiation and propagation of material fracture. After decohesion, the contact is treated assuming the standard contact model with Coulomb friction. The rheological schemes of the contact interaction of bonded and unbonded particles are shown in Fig. 5.
The normal and tangential particle interaction of bonded particles (Fig. 5a) is represented by the springs connected in parallel with the dashpots providing a mechanism to dissipate contact oscillations. Thus, the normal and tangential contact forces are decomposed into the elastic, F ne and F te , 1 In the next part of this section indices denoting the elements will be omitted. and damping parts, F nd and F td , respectively: The elastic normal force is calculated from the linear constitutive relationship: where k n is the interface stiffness in the normal direction, and u n is the change of the distance between the particles evaluated with respect to the distance when the cohesive bond has been established Considering the distance g 0 in the formula for u n allows us to eliminate initial stresses in the specimen due to overlaps or gaps between bonded particles produced by a generation procedure. Thermal expansion is taken into account in calculation of the particle radii: where α i and α j are the linear thermal expansion coefficients, and θ i , θ j , θ 0 i , θ 0 j are the particle current and initial temperatures.
The tangential elastic force is given by the linear equation k t -interface stiffness in the tangential direction, u trelative displacement at the contact point in tangential direction. The relative tangential displacement u t is evaluated incrementally: where u old t is the vector of the relative tangential displacement from the previous time step rotated to the present contact plane and Δu t is the incremental relative tangential displacement with v t being the relative tangential velocity at the contact point determined as where v c r is the relative velocity at the contact point and v rn is the normal relative velocity at the contact point The contact damping forces in the normal and tangential directions are given by respectively, where c n and c t are the damping coefficients. The damping coefficients, c n and c t , can be taken as appropriate fractions, ξ n and ξ t , of the critical damping in the normal and tangential direction, c cr n and c cr t , respectively: For the system of two rigid bodies with masses m i and m j , connected with a spring of stiffness k, the critical damping c cr is given by, cf. [4]: By taking k = k n or k t in Eq. (22) the critical damping c cr n and c cr t is obtained. The normal tensile and tangential elastic contact forces are limited by the tensile and shear interface strengths, R n and R t , respectively: Cohesive bonds are broken instantaneously when the interface strength is exceeded either by the tangential contact force or by the tensile contact force.
In the absence of cohesion the normal contact force can be compressive only (F n ≤ 0) and tangential contact force can be nonzero due to friction. The frictional contact can be treated as a problem analogous to that of elastoplasticity. This analogy allows us to calculate the friction force employing the standard radial return algorithm. First a trial state is calculated and then the slip condition is checked If φ trial ≤ 0, we have stick contact and the friction force is assigned the trial value otherwise (slip contact) a return mapping is performed giving

Background damping
The contact damping previously described is a function of the relative velocity of the contacting bodies. It is sometimes necessary to apply damping, referred here as background damping, for non-contacting particles to dissipate their energy. The background damping is represented by the terms F  (1) and (2). Two types of the background damping have been implemented: -non-viscous damping where ε vt , ε vr , ε nvt and ε nvr are damping constants. The viscous damping is proportional to the magnitude of the velocity, while the non-viscous damping is proportional to the magnitude of the resultant force and resultant moment (Fig. 6).

Thermal formulation of the discrete element method
Thermal formulation of the discrete element method introduced in this work is based on the assumption that the temperature difference inside particles is negligible and the temperature can be considered uniform within particles. Following this assumption heat conduction inside particles is neglected whereas heat transfer to and from particles through their boundary is considered. Such a simplification of the heat conduction problem is typical for the lumped capacitance model, also called the lumped system analysis [4]. This assumption Fig. 6 Scheme of the heat conduction problem for a discrete element is justified for the discrete element model employing relatively small particles and it is consistent with the formulation of the mechanical problem. Similarly as the mechanical problem is governed by the contact interaction, the heat conduction problem is governed by the conductive heat transfer at the particle contacts. The scheme of the heat transfer for a single particle is shown in Fig. 6. The thermal model is expressed mathematically by the heat balance equation, which can be written for a single particle in the following form: where: m i -particle mass, c-specific heat, θ i -particle temperature, Q i -heat sources or heat fluxes per single particle. Q i includes externally supplied heat source Q ext i , heat generated through friction dissipation and absorbed by the particle Q gen i j , heat conducted through the contact interface Q cont i j , and convective and radiative heat transfer between particles and environment on the free surface, Q conv i and Q rad where n c is the number of particles being in contact with the ith particle.

Thermomechanical contact
Except for mechanical effects thermal phenomena, heat transfer and frictional heat generation, are included in the interaction model of contacting particles (Fig. 7). Particleto-particle conductive heat transfer rate Q cont i j is given as follows: where θ i and θ j are the contacting particle temperatures, A cont is the contact area, and h cont is the thermal contact conductance coefficient defining the conductivity through the contact interface between two contacting surfaces. In general case of a contact of different bodies or particles, the thermal contact conductance, is dependent on the surface roughness, material properties, interface temperature and interface pressure [5,21]. It is sometimes convenient to express heat transfer at the contact in terms of one parameter H cont = h cont A cont , then Eq. (34) can be rewritten as [40]: Heat generation through frictional dissipation is calculated using the following formula where F t is the friction force, v ir rt is the irreversible part of the relative tangential velocity, and 0 ≤ χ ≤ 1 is the part of the friction work converted to heat. Heat generated at the contact point is absorbed by the contacting particles Q gen = Q gen i j + Q gen ji (37) If the particles are of the same materials, it is assumed that the heat is absorbed equally by the particles: Q gen i j = Q gen ji = 0.5Q gen (38) In a general case, when the particle thermal properties can be different, heat absorbed by each of the contacting particle is assumed to be proportional to the thermal effusivity of the particle materials with the effusivity e being defined as e = λρc (40) where λ is the thermal conductivity, c is the specific heat capacity, and ρ is the density. The heat partition coefficient α is given by the following equation

Convective heat transfer
Convective heat transfer between particles and the environment can be determined from Newton's law of cooling as [4]: where θ a is the ambient temperature, h conv is the convective heat transfer coefficient, andĀ i is the surface area assigned to the boundary particle as with ξ being the coefficient given by where summation is performed over all the particles on the convective boundary, and L bound is the length of the boundary with the convection condition.

Thermomechanical coupling in the discrete element method
Thermo-mechanical problem defined by the system of coupled equations formed by the equations of motion (1) and (2) and the heat balance Eq. (32) with appropriate initial conditions. Equations (45) and (46) are coupled with Eq. (47) by consideration of the following effects: -frictional heat generation, -thermal expansion of the particles and its effect on particle interaction.
The system of coupled equations is solved using the staggered solution scheme, in which the mechanical and thermal problems are analysed separately.
1. Solution of the mechanical problem.
Equations (45) and (46) are integrated in time using the explicit central difference scheme. Employing the equations for the known configuration at the time t n the solution for the time t n+1 is obtained in the following way: Thermal expansion of the particles and resulting thermal components of the interaction forces are considered in the solution of the mechanical problem. In the considered problem, thermally induced stresses in rock are of minor importance, but in other problems they may be the main factor leading to rock damage [36]. Heat generated by friction is evaluated in Eqs. (45) and (46) and passed to Eq. (47).

Solution of the thermal problem
Heat conduction Eq. (47) is integrated in time using the explicit forward Euler scheme The thermal problem is solved on the modified particle configuration determined in the solution of the mechanical problem. Particle temperatures evaluated in the solution of the thermal problem are passed to the solution of the mechanical problem.
The explicit time integration scheme of the coupled thermomechanical problem is conditionally stable. The time integration step is limited by the critical step Δt cr : where Δt mech cr is the critical time step for the solution of the mechanical problem depending on the highest eigenfrequency of the discrete system max and Δt term cr is the critical time step for the solution of the thermal problem. The critical time step for the solution of the thermal problem can be estimated by the critical value for the one-dimensional heat conduction problem [4].
where l min is the minimum particle centre distance and a is the thermal diffusivity a = λ ρc . (58)

Wear model
Wear of rock cutting tools is dominated by abrasive and adhesive mechanisms [35]. Abrasive wear is one of the most important mechanisms in cutting of hard rocks especially in the presence of quartzite. Scraping of the rock surface leads to high temperatures, which softens the steel body and hardmetal tip, resulting in increasing wear of adhesive character. Adhesive wear is typical for heat generating rock of low abrasivity. Different mechanisms often act in parallel. Wear of rock cutting tools depends on many factors, including rock/soil properties, tool characteristics and cutting process parameters [35]. Based on the experimental data different empirical wear prediction methods for rock cutting tools have been developed, cf. [6,34]. In the approach adopted in this paper, the wear evaluation is included in the simulation of rock cutting process. The main motivation has been to increase possibilities to predict wear of rock cutting tools under different process conditions taking into account both mechanical and thermal effects.

Archard model of wear
Among the quantitative models of wear the classical formula of Archard [1] is still widely used. Its validity for rock cutting has been confirmed recently by experimental studies [37]. The Archard model assumes that the wear rateẇ is proportional to the contact pressure p n and to the slip velocity v Ṫ where H denotes the hardness of worn surface and k is a dimensionless wear parameter. The Archard model was derived originally for adhesive wear, the same form of equation, however, can be obtained for abrasive wear, cf. [27]. Values of adhesive and abrasive wear constants k for different combinations of materials can be determined in laboratory tests.

Wear model accounting for thermal effects
Influence of temperature on wear can be taken into account by adaptation of the law of Archard [1] given by Eq. (59). Thermal effects can be captured approximately by taking hardness H changing with temperature θ Taking this into account in Eq. (59) we havė

Numerical implementation of wear evaluation algorithm
The tool wear w, i.e. the amount (depth) of removed material is obtained by integrating the wear rateẇ in time The Archard law given by Eq. (61) has been implemented in the thermomechanical discrete element model presented earlier. Distribution of wear on the tool surface is estimated and the tool shape is modified according to the calculated wear. Wear is a relatively slow process. In the numerical algorithm developed wear is accelerated using scaled wear constants.
The developed numerical algorithm of wear evaluation takes advantage of the tool discretization with discrete elements. This will allow us to easily modify the shape of the tool. The tool shape is changed eliminating particles if the accumulated wear exceeds the particle diameter. The tool shape is modified continuously during the analysis.

Problem definition
2D simulation of rock cutting with one pick of a dredge cutterhead has been performed in order to verify performance of the discrete element model developed. Cuttersuction dredgers are equipped with cutter head (Fig. 1), which can be used for disintegration of hard rocks. Rock cutting during dredging is performed under water and accurate modelling of rock cutting should take into account influence of the hydrostatic and pore water pressures [35]. Here, these effects will be neglected assuming working conditions for relatively shallow water. Another effect, water cooling of cutting tools, has been considered in the solution of the thermal problem. Figure 8 shows a work scheme of a dredge cutter head. The picks follow cycloidal trajectories resulting from the swing motion and rotation of the cutter head. The teeth cut in a rock during a part of the revolution and cool off in water the rest of the revolution when the teeth rotate in water.
Simulation of a complete cutter head would require excessively large computational cost. Numerical analysis has been performed for one pick in an undercutting operation using   Fig. 9. Three cycles consisting of rock cutting with tool heating and subsequent non-cutting rotation with tool cooling have been analysed. The swing velocity of the cutter head was assumed 0.2 m s −1 , and angular velocity 1.6204 s −1 , which with the distance of the tooth from the axis of rotation 0.7 m gives circumferential velocity 1.134 m s −1 .
Mechanical properties of the rock are assumed the following: Young modulus E = 18 GPa, Poisson ratio ν = 0.18, unconfined compression strength σ c = 60 MPa, and tensile strength σ t = 7 MPa. Thermal properties have been taken as follows: for the rock heat capacity c = 1,970 J/(kg K), heat conductivity λ = 5 W/(m K), and linear thermal expansion α = 10 × 10 −6 K −1 , and for the steel heat capacity c = 450 J/(kg K), and heat conductivity λ = 60 W/(m K). The densities of the rock and steel have been taken,ρ = 7,830 kg m −3 andρ = 2,500 kg m −3 , respectively. Normal values of the heat transfer coefficient for water cooling are between 5,000 and 25,000 W/(m 2 K) [24,35], here a low value of 5,000 W/(m 2 K) has been assumed.

Discrete element model
The rock sample of the shape shown in Fig. 9 is represented by an assembly of randomly distributed 92,000 discs of radii 1-1.5 mm. The compaction of the assembly is characterized by the average porosity n = 0.12. The rigid tool is modelled with 5,400 discs of equal radius of 1 mm. The porosity of the particle discretization amounts to 0.13.
The key issue in the discrete element method simulation is determination of micromechanical parameters yielding required macroscopic properties. There are different theoretical and numerical methods to find micro-macro relationships. A convenient method based on dimensional analysis have been proposed in [14] and successfully employed in [10,28,38]. Using the dimensionless micro-macro relationships given in [28] the following discrete element parameters have been evaluated: k n = 20 GPa, k t = 4 GPa, R n = R t = 0.1 MN m −1 . The Coulomb friction coefficient characterizing the frictional particle interaction after debonding has been assumed μ = 0.8.
The densitiesρ given above for the rock and steel are bulk densities. The particle material ρ density should be evaluated taking into account the porosity: The main micromechanical parameter employed in the thermal formulation is the thermal contact conductance H cont . In case of rock modelling, the discrete element method is used as a discretization method-discrete elements do not represent separate particles which come into contact. The thermal contact conductance for the rock model can be estimated assuming that the heat flux transferred in the contact is equivalent to the heat flux in a bar connecting the particle centres shown in Fig. 10: 2r 2r x i x j Fig. 10 Equivalence of two particles to a bar in a heat conduction problem where A bar is the bar cross-section area and l its length. Taking A = 2rb and l = 2r (r being the arithmetic mean of the contacting particle radii and b-the thickness) we have the thermal contact conductance given by In our case, the thickness b for the 2D model is assumed 1 m. It must be remarked that the micro-macro relationship (65) is approximate only and it does not take into account characteristics of the particle assembly (density of connections or coordination number). Similarly to the analysis performed in [28] for mechanical parameters it would be possible to investigate the dimensionless parameter H cont /(λd) with d being a length parameter (d = b for 2D, d =r for 3D). Convective heat transfer between the tool and water is evaluated using Eq. (42) assuming the water temperature θ w = 283 K. For the rock-tool mechanical interaction the following set of parameters has been assumed: k n = k t = 50 GPa, μ = 0.4. The coefficient χ which determines the part of the friction work converted to heat in Eq. (36) is assumed 0.85. The conductive heat transfer at the contact between the tool and rock is defined analogously to Eq. (42) as: where θ tool and θ rock are the tool and rock surface temperatures and the areaĀ i is defined similarly as in Eq. (42). A constant value of the thermal contact conductance h t−r = 1,500 W/(m 2 K) has been assumed. The set of the data is completed with the parameters of the wear model, the temperature dependent hardness H and dimensionless Archard coefficient k. The rock cutting tool is made of steel whose hardness measured in pressure units (Brinell hardness) as a function of temperature is given in Table 1. The dimensionless Archard coefficient k has been estimated in laboratory tests (Fig. 11) in the range 0.002-0.003. The average value k = 0.0025 has been scale 10,000 times in order to obtain significant shape changes in short time.

Numerical results
Results of the analysis are shown in Figs. 12, 13, 14 Fig. 13. An increase of tool temperature due to heat generated by friction in the contact between the rock and tool can be observed. We can see high temperature at the tool surface in the contact area and a fast drop of temperature with the distance from the surface. Steep thermal gradients are also predicted by other theoretical and numerical studies [20,24,35]. As expected, during the cooling phase the maximum temperature and temperature gradients decrease due to convective heat transfer at the tool surface and heat conduction through the tool body, which can be seen in Fig. 14. Variation of the maximum tool temperature during repeated heating and cooling in water is shown in Fig. 15. The temperature rises during cutting and drops during cooling in water in each cycle, the temperature peaks at the end of subsequent cycles are higher and mean temperature increases. If we continued our analysis, after certain number of cycles, a stationary state in which mean temperature and temperature peaks would stabilize. In this state the amount of energy absorbed by the tool during heating would be equal to that released to the cooling water. Figure 16 shows accumulated wear on the tooth surface at different instants of the first tooth pass of cutting. Maximum wear amounts are localized in the areas with the most intensive friction and heating. Since the wear coefficients have been upscaled the change of the tool shape due to wear is visible even after the first cut (Fig. 17). The tool wear is quantified by the curve showing the tooth mass loss in Fig. 18.

Conclusions
Extension of the discrete element formulation on thermal and thermomechanical problems has allowed us to analyse a complex practical problem of rock cutting taking into account both mechanical and thermal effects and their reciprocal influence. Thermomechanical model of rock cutting enables more realistic treatment of tool wear which depends both on mechanical as well as thermal effects. The thermal and thermomechanical discrete element method formulation contains original elements. The thermal formulation is based on the assumption of the uniform temperature in the discrete element. Heat conduction is modelled by the heat transfer at the contact between particles, which is consistent with the formulation of the mechanical problem in the discrete element method in which macroscopic material behaviour is governed by the contact interaction.
Numerical results confirm a correct performance of the developed thermomechanical formulation of the discrete element method. The results of rock cutting simulation show qualitative agreement with theoretical expectations and practical observations.  Fig. 17 Thermomechanical simulation of rock cutting with a pick of a dredge cutter head-change of the tool shape due to wear. a Initial shape, b shape after the first cut In order to get quantitatively correct results the discrete element modelling should be based on experimentally measured physical parameters and should be extended on threedimensional geometry.