Abstract
This study analyzed the potential of modelling meat mastication by using pork and poultry meat as food with different physical properties under different grilling temperatures. For the purpose of modelling oral processing, temporal dominance of sensations and finite element methods were employed. A panel with ten subjects was trained and used for oral processing analysis and temporal dominance of sensations revealing in-mouth sensations and mastication characteristics. In parallel, the second aim was to evaluate the mechanical properties of the samples and explore the potential of simulating the first bite using the finite element method. Based on the textural parameters, a 3D model of grilled meat was created and a first-bite simulation was performed. A higher level of differences was observed comparing the number of chews for pork meat compared to poultry meat. The chewing rate showed a statistical difference with values in the range of 1.31 chews/s to 1.46 chews/s for pork meat and between 1.36 chews/s and 1.42 chews/s for poultry meat. Firmness was the predominant sensory attribute recognized by panelists at the beginning of mastication, which confirmed our approach used for first-bite modelling. Simulation results show the growth of internal stress following the jaw’s path. Presented models demonstrate that the highest values are around teeth pressure and lead to a conclusion that upon biting, the meat structure will suffer irreversible damage dividing the grilled meat into two pieces, as it happens during the first bite. The main conclusion of this study is that by combining results from oral processing and testing of mechanical properties of the grilled products, it is possible to simulate the first bite.
1. Introduction
Analysis of different structural changes, including in-mouth sensations that occur during mastication, is referred to as oral processing [1]. The complexity of understanding food oral processing consists of three types of evaluations, depending on the time the measurements occur, namely: (i) ex ante, (ii) ongoing, and (iii) ex post. Ex ante analyses mainly target different physical and mechanical properties measured prior to mastication that can help in understanding deformations and fracturing of food that occur in the mouth from the first bite until swallowing and may support explaining eating and sensory perception [2]. The majority of ongoing studies that analyze oral processing during mastication confirm that basic characteristics such as the number of chews, bite size, consumption time, and eating rate are correlated with food characteristics [3–5]. In parallel, the use of temporal dominance of sensation (TDS) helps in understating in-mouth sensation changes [6]. This sensory tool is capable of screening sensations associated with textural and flavor attributes that occur during mastication from the first bite to swallowing [7]. Food structural changes, particle size distribution, and saliva incorporation also help in understanding the structural behavior of food [7–9]. Finally, ex post studies focus on bolus characteristics, swallowing, and effects of food consumption such as energy and nutrient intake [10] and satiety [11]. The majority of these methods are used in terms of interpreting certain oral processing phenomena, rather than in modelling certain mastication activities.
The finite element method (FEM) is an engineering tool that found its application in food science focusing on mechanical properties and predicting the behavior of materials, such as in the work of Celik et al. [12] for the purpose of simulating dropping of apples, Yousefi et al. [13] for drop testing of pears, Haiyan et al. [14] for dropping eggs, or Guessasma et al. [15] for simulating the behavior of starch/protein composite. All these studies used FEM to determine the deformation behavior of food [13]. When it comes to mastication, FEM was used in analyzing stress and loading on teeth [16], occlusal static and dynamic loading during biting or chewing [17], or analyzing simulating mastication associated with different materials used as substitutes for molars [18]. However, all these studies used FEM from a dental rather than a food science perspective. In meat science, the use of FEM supported by mechanical properties is not new but was mostly associated with mass/heat transfer in terms of cooling and/or heating [19–24]. However, there were no studies that simulated the first bite of grilled meat using FEM.
Regarding culinary methods associated with meat, the two basic methods are cooking and grilling [9] with sous-vide [25] and ohmic heating [26] promoted as new promising techniques. Meat that is used as food is a postmortem skeletal muscle tissue of animals [27] as it undergoes a variety of physiological and biochemical changes after slaughtering [28]. The composition of meat depends on numerous characteristics such as species/breeds, age, and also the muscle position in the carcass [29], where the complex structure of meat comprises fibers, intramuscular connective tissue, and intramuscular fat [30]. The complexity of biological materials influences mechanical properties and, as such, cannot be easily explained like in the case of less complex materials such as metals and alloys [31]. When it comes to solid food, there is a certain barrier in connecting the field of material sciences with biological sciences [32], mainly due to the matrix of food tissues and certain limitations in using mechanical tests per se ignoring a wider (food science) perspective. In that sense, meat is considered one of the most complex types of food with interacting components operating at the microscale [33]. Therefore, modelling meat as a material and meat consumption as a process is considered a scientific challenge.
This paper had two main objectives. The first was to characterize sensory and oral processing characteristics of pork and poultry meat as food with different physical properties under different grilling temperatures grilled at three temperatures and identify parameters of interest for simulating mastication from the first bite. The second aim was to evaluate the mechanical properties of the samples and explore the potential of simulating the first bite using the finite element method.
2. Materials and Methods
2.1. Pork and Poultry Meat
Pork (m. longissimus dorsi) and poultry (m. pectoralis major) meat was purchased locally from a butcher store in Belgrade, Serbia. After purchasing, meat was cut into uniform pieces of 200 gr (pork) and 100 gr (poultry), vacuumed, and stored in the refrigerator at 4°C. Before grilling, meat was placed at room temperature for half an hour.
2.2. Cooking Loss
Prepared uniform pieces of meat were grilled using a Tefal grill (OptiGrill+) on three defined grilling regimes (T1 is the shortest grilling time; T2 is the medium grilling time; and T3 is the longest grilling time) in accordance with the instrument’s manual for pork and poultry meat (time duration between four and seven minutes). Grilling temperatures on the surfaces show a range between 263.47°C and 272.1°C for pork meat and 238.42°C and 249.24°C for poultry meat, with medium grilling time having the lowest temperature. The achieved thermal temperature in the center of the product was from 67.17°C to 90.63°C (pork) and 78.77°C to 99.3°C (poultry), confirming its safe consumption.
The indicator consists of a sensor measuring the grilling meat thickness, based on which the instrument indicates the doneness level by changing the LED bulb light color. To control the process, temperatures in the samples’ center during grilling were measured with a digital thermometer (Trotec GmbH - Model BT20, Germany), while the temperature of the grilling surface was measured using an infrared thermometer (TES-1327KUSB) (Table 1). Grilled samples were settled to rest at room temperature and then sealed in vacuum bags and stored in the refrigerator at 4°C.
The weight of the samples before and after grilling was measured using an analytic balance (OHAUS Adventurer-Model AR2140, USA), and the cooking loss was calculated by measuring the weight differences [25].
2.3. Density
Upon grilling, density was measured by determining the sample mass and volume [34] and measuring the external dimensions (length, width, and thickness). Grilled samples were cut into cubical samples (20 × 20 × 20 mm) using a thin-bladed sharp knife. Assuming a constant cubic shape, based on 30 replicates, for both types of meat and three grilling temperatures, density was calculated as shown in Table 2. Dimensions of the samples were measured with a digital vernier caliper while mass was measured using an analytical balance (OHAUS Adventurer-Model AR2140, USA).
2.4. Sensory and Oral Processing Evaluation
The panel consisted of 10 panelists experienced in both sensory and oral processing evaluation participating in the study. All panelists were in good general health condition, with no reported dental problems and with a normal range for BMI of 18-25 kg/m2, as suggested by Forde et al. [35]. All panelists are well-trained, have already participated in sensory and oral processing studies associated with animal-origin food, and have passed basic tests for taste, odor, and texture in line with ISO 8586 [36]. Two training sessions (duration 2 hours each) were performed for sensory evaluation, while for oral processing methods used, one 2-hour session was organized as suggested by Djekic et al. [37].
2.5. Temporal Dominance of Sensations
For performing TDS evaluation, six sensory characteristics have been predefined, as follows: firm, meat flavor, juicy, dry, fibrous, and soft, modified from References [9, 38, 39]. The meanings of each sensation were explained to the panelists [40] during the first training session for sensory evaluation. Definitions that were used were as follows: firm/soft-strong/weak force required to compress food between molars; juicy/dry-high/low quantity of juice released from food during mastication; fibrous—the presence of fibers in meat samples; meat flavor—the strength of all flavors typical for pork/poultry meat. Reference food used during the training session was pork and poultry meat prepared in three culinary methods—grilling, cooking in boiling water, and sous-vide. Also, panelists were instructed to select a dominant sensation only when they perceived one, and they could choose the same attribute more than once.
The grilled meat was cut into cubical samples (20 × 20 × 20 mm) and provided to panelists for TDS. This test was completed using an internally developed application and in line with the method explained in Djekic et al. [9]. When mastication started, the panelists selected “start,” and dominant sensations were recorded. At the moment of swallowing, panelists were instructed to select “stop.” Testing of all grilled samples (two types of meat, three grilling temperatures) was performed in three replications. Samples, coded with three-digit codes, were served to the panelists randomly, one at a time on two consecutive days—the first day for pork meat and the second day for poultry meat. To clean the palates and remove aftertastes between two mastications, panelists used white bread and tap water.
To calculate TDS curves, the following rules have been applied: for each point on the time scale, the value of the dominance rate (overall share of a specific attribute) was determined by dividing the number of mentions of a sensation divided by the overall number of sampling (from all replications and panelists) [41]. The “rule of the thumb” in interpreting the curves is that “the higher the proportion, the higher the agreement among panelists about the dominating sensation over the period.” All plotted results were standardized from the first bite to swallowing [9, 42].
2.6. Oral Processing
The grilled meat was cut into cubical samples (20 × 20 × 20 mm) and presented to panelists for mastication. The mass of every sample was measured with a technical balance of 0.01 g accuracy. A digital video camera was positioned in front of the panelist, so the upper part of the person chewing was visible, as proposed by Forde et al. [35]. All members of the panel chewed all grilled samples (two types of meat, three grilling temperatures) in three replications. Upon completion, video clips have been analyzed using a stopwatch, resulting in the number of chews and total oral exposure time [43, 44].
Based on these data, the following oral processing parameters have been calculated: the number of chews, consumption time per sample, chewing cycle duration (s/chew), chewing rate (chews/s), eating rate (g/s), and average bite size (g) [4, 35].
2.7. Mechanical Properties
Mechanical properties of two types of meat samples have been determined under atmospheric temperature/pressure conditions (20°C, 1 atm) using Brookfield CT3 Texture Analyser for compression and TA.XT plus Texture Analyser, Stable Micro System for Warner–Bratzler shear tests. Compression test parameters were as follows: test type—compression, test speed—1 mm/s, trigger load—10 g, target mode—30%, load cell—10 kg, probe selection—cylindrical probe (plate of 50.8 mm diameter), and sample shape—rectangular. This test was performed on all three planes of the cubic sample (Figure 1).

Warner–Bratzler shear test had the following parameters: test speed—1 mm/sec, target mode—distance (21 mm), sample shape–rectangular, selected probe—HDP/WBV (Warner-Bratzler V slot blade), and load cell—50 kg. This test was performed normally on the axis of the muscle fibers [9]. In such a way, the tool cuts the samples, so the shearing is perpendicular to the longitudinal positioning of the fibers [45].
All meat samples that were cut off along the axis of the muscle fibers had cubic dimensions (20 × 20 × 20 mm), and tests were performed in 15 replicates. Dimensions of the samples were manually measured with a digital vernier caliper. Preparation of the samples was performed using thin-bladed sharp knives to minimize the damage to the fibers and take into account the direction of the fibers [46]. As an output of this test, firmness (N) was obtained.
True stress and strain calculations have been obtained according to the research of Vallespir et al. [47] and Vicente [48] (equations (1) and (2)). Rupture stress (σR, MPa) and strain (εR) have been extracted from the first peak of the stress-strain curve using TexturePro CT V1.9 Build 35 software. Young’s modulus (Ed) was calculated for the common linear part of the stress-strain curve (equation (3)).
Legend: F(t) is the force at time t; Ho is the initial sample height; ΔH is the height difference; H(t) is the height at time t; and Ao is the sample area.
From material science, it is known that depending on the mechanical properties under load, isotropic materials are the ones directionally independent as opposed to orthotropic (interchangeable across the three orthogonal axes) or anisotropic having different properties in all directions [32]. As meat is a complex material, for this study, it was considered an orthotropic material with three planes/axes of symmetry (Figure 1). Therefore, the calculation of all the abovementioned parameters has been performed for all three planes.
To calculate Poisson’s ratio, it has been assumed that during compression tests, expansion of meat in the direction perpendicular to the specific loading direction is equal in the other two planes with constant volume before and after loading, equations (4) and (5):where Vb is the volume of the sample before compression; Va is the volume of the sample after compression; x, y, and z are the lengths of the cube in three planes; assumed compression axial to x-axis.
Poisson’s ratio as the ratio between the transversal (lateral) strain and the longitudinal strain, in a tensile or compressive test, is presented as follows:where νij is corresponding to an expansion in direction when compression is applied in direction.
2.8. Simulation of the First Bite
In this study, a 3D solid model of grilled meat was created by using solid modelling software. The mesh has been developed in the SolidWorks Simulation FEM code. The tetrahedral solid element type was chosen with 34,146 elements and 40,228 nodes used for mesh construction of the 3D model. A similar mesh was used in the work of Wang and Sun [21] that modeled roasted meat with four-node tetrahedral elements.
For the purpose of performing this simulation, the following assumptions apply: (i) 3D solid meat is presented as a cube (20 × 20 × 20 mm); (ii) first bite occurs perpendicular to the longitudinal positioning of the fibers; (iii) first bite occurs parallel to the x-axis; (iv) first-bite force was assumed as the value obtained from WB test for all types of grilled meat (Table 3) divided by two (upper and lower jaw); and (v) first bite was calculated as line pressure (20 mm) assuming teeth (incisors) width of 1 mm. To simulate the first bite, values for density, Young’s modulus, and Poisson’s ratio for the three orthogonal axes will be used.
2.9. Statistical Processing
One-way ANOVA and Tukey’s HSD post hoc tests were used to distinguish statistical differences between grilling time as the heating factor, while meat type was not considered as factor for statistical analysis. The level of statistical significance was set at 0.05.
3. Results and Discussion
3.1. Oral Processing Characteristics
Tables 2 and 3 give an overview of oral processing parameters for two types of meat. For pork meat, (Table 2) the number of chews was between 27 and 30, with an average mastication time between 19 and 22 seconds. The chewing cycle was around 0.75 s/chew, with the chewing rate up to 1.46 chews/s. The lowest average bite size was for medium grilling time. Oral processing characteristics of poultry meat (Table 3) fall in a narrower range: number of chews between 22 and 24; chewing cycle between 0.73 s/chew and 0.76 s/chew, and chewing rate up to 1.42 chews/s. More pronounced differences between the number of chews were observed for pork meat (Table 2) grilled at the highest temperature compared to poultry meat (Table 3), but with no statistical difference (). This pattern is expected as higher temperatures trigger crust formation and heat-related toughening associated with harder meat [25].
The chewing rate was the only parameter with an observed statistical difference (only for pork meat) with values in the range of 1.31 chews/s-1.46 chews/s (pork meat, Table 2) and 1.36 chews/s-1.42 chews/s (poultry meat, Table 3). These values are slightly below the mean chewing rate of 1.53 chews/s outlined by Farooq and Sazonov [49]. It needs to be mentioned that meat as a very complex material influences oral processing parameters to fall out of typical, expected patterns observed in solid food such as cheese or chocolate [4]. Mastication pattern is influenced by the mechanical properties of meat, a mixture of saliva and liquid released from the meat, and culinary treatment applied to meat [9, 50, 51].
3.2. Temporal Dominance of Sensations
When the aim of performing sensory analysis is to examine textural and flavor changes during mastication, TDS is recognized as one of the most promising tools and has been used in meat science [7, 9]. It is a dynamic sensory method with the aim of capturing “dominant sensation(s)” while consuming food [52]. Grilled meat is unique in its flavor and taste, and grilling is one of the most popular culinary methods used at home and/or restaurants [53]. A similar approach in employing TDS with pork loins was performed by Watanabe et al. [54].
Figure 2 presents the results of TDS applied on evaluating pork and meat products treated under three grilling temperatures, showing the proportion of citations for each sensory attribute. The firmness of the samples dominated in pork meat products during the beginning of mastication, similar to the results of Watanabe et al. [54]. In parallel, juiciness (for T1 and T2) and dry attribute (for T3) dominate until half of the consumption time. These patterns may have been expected as culinary methods using temperatures above 140°C enable the development of odor and meat flavors [55], preventing juice losses and textural damages [56]. In the second half of the mastication process, the softness of the sample (and bolus) was predominant. For all three grilling times, in the last third of the mastication period, no sensation had over 30% of dominance rate.

(a)

(b)

(c)

(d)

(e)

(f)
As for poultry meat, the flavor was very dominant at the beginning of the chewing process for T1, with firmness being firstly observed by panelists for T2 and T3. Juiciness was recognized as the most dominant in the middle of the mastication period. During the last period of mastication, similar to pork meat, softness was prevailing. Knowing that poultry meat is slightly softer than pork meat, it is not surprising that flavor and juiciness are observed at the beginning of the consumption process, as chewing affects the perception of juiciness [57] and the release of aroma [52].
For both types of meat, at T3, the dry sensation was also observed by panelists due to the fact that the final internal temperature of pork and poultry at the longest grilling time was the highest, corresponding to previous research that internal temperature above 80°C decreases tenderness [58, 59]. As the first dominant attribute, “firmness” is common for all samples; this finding is valuable for understanding how to simulate the first bite.
3.3. Mechanical Properties
From Figure 1, it is obvious that the grilling surface was set on the plane perpendicular to fibers (plane xOy). A higher level of cooking loss was observed for poultry meat (up to 40%) than that of pork meat (almost 30%) (Table 4). Firmness obtained by the Warner–Bratzler shear force test shows that the values increase with the increase in grilling temperature. This test measures the maximum force needed to cut off/shear a meat sample [45], and this value was used for modelling the first bite.
Values of Young’s modulus increase with the increase in grilling temperature. On the other side, due to the complexity of meat, values of this modulus differ on all three planes with no clear pattern. The Poisson’s ratio was between 0.2460 and 0.4913 for pork and between 0.2296 and 0.4909 for poultry meat. Young’s modulus in the study of Jahanbakhshian et al. [60] for chicken nugget crust and the crumb was 0.26 and 0.37, respectively. Kim and Hung [61] reported that Poisson’s ratio of meat products was in the range of 0.2 to 0.4, while in the work of Nowak et al. [62], the average value of Poisson’s ratio of various meat products was 0.49. Most soft tissues are considered roughly incompressible materials with a Poisson ratio up to 0.49 [63]. This means that if the tissue is compressed axially (in one direction), then it must expand laterally in the other two directions [64]. This assumption was used in our study to calculate Poisson’s ratio values in all three directions.
3.4. First Bite Simulation
For simulating the first bite, calculated values for density, Young’s modulus, and Poisson’s ratio for the three orthogonal axes (Tables 2–4) have been used to set the parameters of meat as an orthotropic material. Six models have been developed (pork meat T1, pork meat T2, pork meat T3, poultry meat T1, poultry meat T2, and poultry meat T3). Pressures applied at the first bite were in the range of 1.46-1.98 N/mm2 (pork) and 0.85-1.07 N/mm2 (poultry). Values were calculated using firmness values for the Warner–Bratzler test divided by 40 (length of incisors that adhere to the surface 20 mm, incisors width 1 mm, 2 jaws—upper and lower and assuming both jaws have equal pressure, so we have pressure area of 20 mm2 for the upper jaw and the same value for the lower jaw, in total 40 mm2).
According to the FEM simulation results, maximum values for von Mises stress were obtained for the highest grilling temperatures (1.98 N/mm2 for pork meat T3 and 1.07 N/mm2 for poultry meat T3). The von Mises stress is often used in biomechanical calculations as the stress value enables predicting the food damage behavior [65].
A comparison of results shows that higher grilling time corresponds to higher values of van Mises stress for both types of meat. Values of maximum von Mises stress rise from 1.181 N/mm2 to 1.450 N/mm2 for pork and from 0.686 N/mm2 to 0.909 N/mm2 for poultry meat. These computed values by using FEM enable prediction of how deformation will occur during the first bite [12], knowing that the average duration of the chewing cycle (and the first bite) is between 0.73 and 0.78 s/chew for both types of meat (Tables 2 and 3). The 3D models (Figure 3) show the development of internal stress in the first-bite direction leading to crack propagation and breakage of the sample by the pressure of the upper and lower jaw as expected during the first bite. The highest values for all six samples are in the area of teeth pressure, indicating that upon impact, the structure will most probably experience irreversible damage leading to diving the sample into two pieces. These results show that the use of FEM in combination with experimental data associated with food materials may enable fair prediction of crack initiation and propagation [66]. This is the first step in modelling the food mastication cycle from first bite breakdown, the interactions between comminuted food and saliva during bolus formation, to the release and transport of taste and aromas [67]. Future challenges in modelling oral digestion should focus on (i) widening the perspective of understanding mechanic of food to how they flow/deform during mastication and (ii) understanding physicochemical reactions supported by the role of enzymes in the dissolution of taste compounds into saliva and release of aromatics in the mouth [68].

(a)

(b)

(c)

(d)

(e)

(f)
TDS curves give an overall perception of major sensory attributes of grilled pork meat during mastication. It is obvious that differences in the dynamic sensory perception of the samples exist in terms of whether textural or taste/odor attributes prevail, depending on the type of meat and grilling parameters. However, as firmness was the predominant sensory attribute recognized by a panelist at the beginning of mastication, simulation of the first bite using mechanical characteristics may pave the way for simulating the entire grilled meat mastication process. Overall, chewing cycles mainly serve for food breakdown initiating it from the first bite [52]. In parallel, when mastication starts, it triggers saliva incorporation in the mouth acting as a solvent for releasing different taste-active compounds [69]. Awareness of in-mouth changes obtained during TDS may be useful not only for interpreting TDS graphs but also for designing new types of food associated with recognized sensory profiles [70]. However, it must be mentioned that as in the case of all types of organic materials, the complexity of the material, variability of different impacts, and sample orientation influence results [12].
4. Implications for Sensory Science
This paper provides new insights into explaining changes in mechanical properties, mastication, and dynamic sensory perception of grilled meat ending with the modelling of the first bite. In addition to the increased need for understanding food oral processing and digestion [71, 72], authors believe that the application of such modelling can help in revealing not only food breakdown [73, 74] but also flavor release [1], important in gastronomy and sensory science. Our results send a message that a variety of mastication aspects should be considered for understanding different gastronomic factors associated with grilled meat. By considering the information provided in this paper, chefs may find ideas for adapting grilling methods for target groups of meat consumers. Finally, results may be also of interest to various gastronomy stakeholders in the future development of tailor-made food for different groups of consumers such as the elderly, dental wearers, or children.
5. Conclusions
This study highlights the potential of using finite element modelling in oral processing studies and is one of the first of its kind that tried to connect mechanical properties in modelling the first bite. This approach has been used as a first step in simulating meat behavior during mastication, as TDS confirmed that the firmness of grilled meat is a sensation attributed to the first bite while oral processing showed that regardless of the type of grilled meat and grilling temperature used, duration of the first bite is between 0.73 s and 0.78 s. This model may be a starting point in defining more accurate methods for simulating the mastication of complex solid materials such as meat.
The limitation of the study may be the fact that the simulation of the first bite was based only perpendicular to the fibers considering that the direction of fibers during the first bite may not be relevant to average consumers. Also, realistic jaw motions as well as the geometry and rigidity of the teeth have been simplified.
Future research on 3D modelling of the mastication process of pork meat could include mechanical properties of meat bolus at different stages of mastication, including saliva incorporation that affects changes observed during TDS.
Data Availability
Data are available upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The results within this research originated from a Proof of Concept project #5229 “Design of artificial masticator for modelling food oral processing” financed by the Innovation Fund from the budget of the Government of the Republic of Serbia, the Ministry of Education, Science and Technical Development through the “Competitiveness and Jobs Project.”