Intrinsic mechanical properties of food in relation to texture parameters

The texture profile analysis test is an imitative test to determine texture properties of food, which quantify the consumer’s perception of eating food. The instrumental texture parameters obtained from this test depend on the specimen size and the nonstandardized test conditions. To overcome this problem, texture properties are here related to intrinsic mechanical properties, which are independent of the test conditions. Two types of materials are used to investigate the effect of viscoelasticity, plasticity and damage on the texture parameters for varying test conditions. Analytical relations between mechanical properties, test conditions, and the instrumental hardness, springiness, cohesiveness, and adhesiveness are determined. The hardness is obtained from the stiffness of the material, which is potentially rate-dependent, and the yield stress of a material in case of plasticity. The springiness is determined by the recoverable or irrecoverable strain in the material, which results from the mechanical properties in combination with the test conditions. Cohesiveness and springiness are found to be strongly related, unless structural damage is present in the material. Adhesiveness is only an indirect measure of the adhesion between the material and compression plate and depends on the test conditions and stiffness of the material. Finite element simulations reveal a decrease of hardness in case of a nonflat top surface, indicating the importance of geometrical effects.

processing (Pascua et al. 2013;Chen 2015). The purpose of the test procedure is to deliver an objective measurement related to perceived sensory attributes. These sensory attributes are evaluated by a trained test panel. The hardness is evaluated as the required force to compress a product by a certain amount, the springiness as the recovery of the deformation, and the cohesiveness as a measure of the deformation before rupture. The deformation is applied by the test panel participant with the teeth or fingers (Di Monaco et al. 2008;Pascua et al. 2013). Besides these three main parameters, many more texture properties are defined, such as adhesiveness, chewiness, and resilience. Many researchers have investigated relations between the sensory evaluation of the test panel and the instrumental texture parameters, obtained by a TPA test. Meullenet et al. (1998) tested a wide variety of food products and found high correlations for the hardness and springiness, but not for cohesiveness and chewiness. However, Wee et al. (2018) did a similar study and found significant correlations for springiness, cohesiveness, chewiness, and resilience but not for the hardness. Breuil and Meullenet (2001) found reasonable predictions for the texture in cheese using an extended instrumental test procedure by considering cone penetration and needle puncture tests next to the uniaxial compression test. Paula and Conti-Silva (2014) used different probes with the texture analyzer as well to study extruded snacks.
Comparing these disparate findings, suggests that the nonstandardized test conditions of the TPA test make the interpretation of results challenging, if not troublesome. The obtained texture parameters differ when using different test parameters (i.e., applied timedisplacement profile) with varying types of probes on different specimen sizes and shapes (Rosenthal 2010;Alvarez et al. 2002;Peleg 2006Peleg , 2019Hyldig and Nielsen 2001;Patel et al. 1992). Clearly, this makes it problematic to compare results. Moreover, it is questionable how meaningful these instrumental texture parameters are in terms of sensory attributes. Consequentially, there is a need for well-defined properties that are ideally independent of the test conditions (Peleg 2019).
Mechanical properties are intrinsic material characteristics that fulfill these requirements. In the TPA test the food material is mechanically loaded under the chosen test conditions. Although instrumental texture parameters are believed to be relevant to the sensory attributes, relating these attributes to mechanical properties that influence the texture parameter may provide novel insights. Di Monaco et al. (2008) related sensory attributes to the results of cyclic compression tests and stress-relaxation experiments, whereas Vincent and Saunders (2002), Kim et al. (2012), and Vliet and Primo-Martín (2011) linked the sensory evaluation of brittle food to fracture mechanics. These studies directly relate mechanical properties to sensory evaluation, without considering the TPA test in relation to mechanical properties.
The TPA test is a test method with a lot of disadvantages and the usefulness of the parameters that are normally derived from the test is limited at best. However, the results of a double compression test are clearly influenced by the mechanical properties of a material. The goal of this study is to analyze the mechanical response during the TPA test and to identify relations between the instrumental texture parameters and intrinsic mechanical properties. This knowledge is useful for multiple purposes. Firstly, it provides insight into the effect of varying test conditions. Secondly, it shows which intrinsic material properties affect specific texture parameters, which also indicates suitable material characterization methods. Moreover, such relations are useful in designing customized food as mechanical properties can be tailored by the microstructure, which is altered by the ingredient composition or processing conditions (Attenburrow et al. 1989;Sozer et al. 2011;Gao et al. 2018), for example, by 3D printing (Le Tohic et al. 2018). To achieve these relations, two different types of materials, namely a linear viscoelastic and an elasto-viscoplastic material, are represented by a suitable constitutive model. The effect of structural damage is included. The modeling approach is, however, limited to food products that remain macroscopically intact during the TPA test. Using these models, the intrinsic material response during the TPA test is described in detail. Each material has its own characteristic mechanical properties, for which analytical relations with the resulting texture parameters are established. The used probe is a flat compression plate, because other types of probes inherently result in heterogeneous deformation that depends on the geometry of the probe. The considered texture parameters are hardness, cohesiveness, and adhesiveness, for which the original definitions (Friedman et al. 1963) are used, and the springiness (Meullenet et al. 1997). The presented methods are easily extended to additional texture parameters such as resilience and chewiness or to modifications of the original definitions (Peleg 1976;Bourne 1978;). In addition, finite element simulations are performed to investigate the effect of friction and surface roughness on the resulting texture properties.
2 Texture profile analysis test 2.1 Test procedure Figure 1 shows the result of a typical texture profile analysis test, in which the applied compressive force is measured as a function of time. First, the compression plate is moving down with constant velocity v to a maximum applied displacement u t , inducing a compressive load on the sample. The time this takes is defined as the test time t t according to (1) Next, the compression plate is moving up with the same velocity, reaching its original position at time 2t t . The compressive force in the sample decreases to zero if there is no adhesion, or to a tensile force if there is adhesion. Subsequently, there is a waiting time t w during which the plate remains at its starting position, after which the loading-unloading step is repeated. During the unloading step, the sample and plate may lose contact after time t lc , counted from the moment at which the unloading started. The compressive strain that is still present in the sample may be recoverable or irrecoverable. The plate starts to move down again at time 2t t + t w . Contact is made again (i.e., the force increases) after time t mc relative to the start of the second cycle. It therefore takes a time interval t rec between losing contact with the sample and touching it again in the second cycle.
From the TPA test, multiple texture parameters can be extracted. In this research, focus is put on the parameters hardness H , adhesiveness A 3 , springiness S, and cohesiveness C. The hardness is the maximum force in the first cycle. The adhesiveness A 3 is equal to the area of the force-time response, for a negative measured force, in the first cycle. The springiness and cohesiveness are defined from Fig. 1 as where s 1 is equal to the test time t t , s 2 is determined by the time between making contact in the second cycle and the test time, that is, s 2 = s 1 − t mc , and A 1 and A 2 are the areas under the loading-unloading curve (for a positive measured force) in cycles 1 and 2, respectively. The test conditions of the TPA test are determined by the maximum applied displacement u t , the velocity v, and the waiting time between the two cycles t w .

Finite element model
Although the aim is obtaining analytical relations between texture parameters and intrinsic mechanical properties, additional simulations of the TPA test are performed using the finite element package MSC.Marc to take the nonlinear material behavior and heterogeneous deformation into account. The considered geometry is a cylinder with a length and diameter of 1 cm. The vertical displacement of the bottom of the cylinder is suppressed, and the compression plate is displaced according to the test parameters u t , v, and t w . The deformation is homogeneous in an ideal test situation, meaning that the sample is flat, the sample and compression plate are perfectly aligned, and friction is negligible. In the homogeneous case, modeling the geometry by a single 2D axisymmetric linear element is sufficient to simulate the TPA experiment, using a time step size of 5.0 · 10 −4 s.
In a real experiment, friction and surface roughness are potential sources that may cause the deformation to be heterogeneous. To investigate these effects, the geometry is modeled using 400 axisymmetric linear elements. The effect of surface roughness is investigated by applying a gradual increase of the sample length from 0.95 cm at the boundary to 1.00 cm at the center of the cylinder. The corresponding mesh of a sample with this arbitrary irregular surface is depicted in Fig. 2.
A smoothed Coulomb friction model is implemented by relating the tangential friction force F t to the normal force F n , the coefficient of friction μ, and the relative sliding velocity v r : A coefficient of friction μ = 0.3, corresponding to friction between starch and steel, is chosen as a reference (Emami and Tabil 2008).

Uniaxial material response
The test conditions are converted into relevant quantities for material characterization. The engineering strain ε and strain rateε are calculated for a deformation u, velocity v, and initial specimen length L 0 using The true strain definition, denoted byε, is accounting for the deformation at each instant of the deforming sample length and is applied for large deformations: The maximum applied deformation u t is converted into the maximum applied engineering strain ε t and true strainε t by substituting u = u t in Equations (4) and (5), respectively. Because the velocity is kept constant in the TPA test, the true strain rate is not constant during the measurement. The engineering strain rateε t is therefore used as the relevant test parameter. Note that ε t andε t are negative in the TPA test, because the material is loaded in compression. The strain rateε t is negative during loading and positive during unloading. Thereforeε t is defined positive.
During the test, the compressive force F is measured. The engineering stress σ and the true stressσ are calculated using the initial cross-sectional area A 0 and the corresponding deformed area A: The difference between the engineering and true stress is important because changes in the cross-sectional area during deformation are generally hard to measure. This prevents the direct measurement of the true stress. The adhesion between the material and the compression plate is characterized by an engineering stress σ a or true stressσ a , which implies that they are separated when the tensile stress in the material reaches this adhesion stress.
To investigate the relation between mechanical properties and texture properties, two types of materials are considered that exhibit a different characteristic mechanical response,  Figure 3 sketches the characteristic force response of these materials for a displacement that is first increasing and subsequently kept constant in time. The key characteristic of a linear viscoelastic material is its time-dependent response. The stiffness depends on the applied deformation rate, and the force increase is, in general, nonlinear, whereas the force is decreasing if the displacement is kept constant. The stiffness of an elasto-viscoplastic material does not depend on the deformation rate. An elastic regime is followed by a plastic regime in which a small increasing force is measured due to an increasing cross-sectional area. The plastic strain in the material is irrecoverable, whereas the strain in the linear viscoelastic material is always recoverable. By introducing damage the force decreases with increasing plastic deformation.
First, the constitutive models of the two materials are presented, after which the constitutive equations are solved for the different steps in the TPA test. From this, analytical relations are obtained to link the texture properties to the mechanical properties.

Tensor notations
The constitutive models are formulated in a tensorial form in the three-dimensional space. Using a Cartesian vector basis {e 1 , e 2 , e 3 } and summation over repeated indices, the secondorder tensor A and the transpose A T are denoted as The second-order tensor C results from the inner product of tensors A and B according to and the scalar c from the double inner product of two tensors is The hydrostatic part A h and deviatoric part A d of a tensor are given by where I = e i e i is the unit tensor, and the trace of a tensor is The material time derivative is denoted asȦ.

1D constitutive model
The standard solid model is a time-dependent model to describe linear viscoelastic materials. The mechanical analog of the model consists of a Maxwell element, with Young's modulus E 1 and viscosity η, and a parallel spring with Young's modulus E ∞ ; see Fig. 4. The characteristic time of the model is the relaxation time λ given by The formulation of this standard solid model is restricted to small strains. Therefore analysis is limited to a maximum compressive strain of 0.1, at which the engineering stress and strain are approximately equal to the true stress and strain. The total stress is calculated from the stress in the parallel spring σ ∞ and the stress in the Maxwell element σ M using whereε v is the strain rate in the viscous part of the Maxwell element, obtained by subtracting the strain rateε e in the elastic part as follows: By substituting Equation (14) into Equation (13) the differential equation that relates the stress and strain is given by By introducing the total stiffness E and the fraction of the material response that is viscoelastic φ, the stress-strain relation given by Equation (15) is rewritten as The material parameters are easily determined from a stress-relaxation experiment; see Fig. 5. The stress is measured as a function of time for an applied strain ε, which is kept constant during the experiment. The time t step , in which this strain is applied, is typically smaller than the relaxation time to accurately measure the mechanical properties. The material reacts initially with a total stiffness E. After loading, the stress is exponentially decreasing toward a plateau value that corresponds to the stress in the parallel spring. The relaxation time λ is obtained by linearizing the curve around t = t step and determining its intersection point with the plateau stress E ∞ ε.

3D constitutive model
The 3D integral formulation of Equation (15), which is a built-in material model in the finite element package MSC.Marc, is withε ε ε h andε ε ε d the hydrostatic and deviatoric parts of the strain-rate tensor. The timedependence of the shear modulus G and bulk modulus K is expressed as of which the shear and bulk moduli are related to the corresponding Young's moduli and Poisson's ratio ν by The influence of the Poisson's ratio on the TPA test results is small. A standard value of the initial Poisson ratio's ν = 0.15 is therefore chosen in the analysis. The 3D formulation is used in finite element simulations to include the effect of friction and surface roughness. Without these effects, a uniaxial stress state is recovered. This implies that exact analytical relations between the TPA parameters, test conditions, and material properties are then derived by solving Equation (15).  The loading response of the linear viscoelastic material for compression cycle i with a strain rateε = −ε t is given by

Analytical solution
with the initial state given by time t 0,i , strain ε 0,i , and stress σ 0,i = 0. In the first compression cycle (i.e., i = 1) it holds that t 0,i = 0, ε 0,i = 0 and the loading step is ended at time t 1,1 = t t . The second loading step starts at t 0,2 = 2t t + t w + t mc and is ended at time t 1,2 = 3t t + t w . The time it takes to retrieve contact t mc is obtained by calculating the material response throughout the whole first cycle and during the waiting time, which also provides the initial condition ε 0,2 .
The uniaxial unloading response of the material withε =ε t is given by in which the initial conditions for the first and second cycle are t 1,1 = t t and σ 1,1 = σ (t t ) and t 1,2 = 3t t + t w and σ 1,2 = σ (3t t + t w ), respectively. Equation (22) is only valid when the material is in contact with the compression plate, which is not the case during the entire Fig. 8 Schematic representation of the different regimes during unloading. The loading step is finished at time t t (a). In the adhesion-free case the material is in contact with the plate until σ = 0 is reached upon unloading (b). In case of adhesion, contact is lost when the adhesion tensile stress σ a is reached (c), followed by an instantaneous decrease of the sample height (d), after which the strain recovers toward zero (e) unloading step. Figure 8 shows the different regimes. At t t the plate starts moving upward: the sample is initially in contact with the plate, that is,ε =ε t . Upon unloading, the material reaches a stress equal to zero after a time interval t 0 lc,1 . The sample is still in contact with the plate until the tensile stress in the material reaches the adhesion stress σ a at t t + t lc,1 .
Contact with the compression plate is lost when exceeding σ a in tension, and at this moment the stress drops instantaneously to zero. The time at which contact is lost depends on the material properties, the adhesion stress, the applied strain, and applied strain rate; t lc,i is the time interval it takes until the stress in Equation (22) reaches σ a , which results in in which (x) denotes the principle branch of the Lambert-W function of the argument x (Corless et al. 1996). Note that the moment of losing contact in the adhesion-free case equals t 0 lc,i , which is found by substituting σ a = 0 into Equation (23), and t a is the additional time it takes to reach the adhesion stress σ a starting from σ (t t + t 0 lc,1 ) = 0: Figure 7 reveals the strain and stress at the moment of losing contact. The strain in the material at an infinitesimally small time before the moment of losing contact is because the strain in the material is prescribed by the compression plate. In the case of adhesion, this is followed by an instantaneous change in stress and strain, exactly at the moment of losing contact. The stress becomes equal to zero because there is no external force on the material. The change in compressive strain is characterized by the linear elastic response of the material due to the infinitely fast stress rate. The strain in the material at an infinitesimally small time after the moment of losing contact is From the moment contact is lost the material is unconstrained, and the strain will recover toward zero. This material relaxation between cycles 1 and 2 results from solving Equation (17) with boundary conditions σ = 0 andσ = 0: The springiness and cohesiveness depend on the deformed sample height at the time contact is re-established. The time between the moments of losing contact and making contact is the time the strain recovers toward zero and is defined as t rec in Fig. 1. To determine the time t rec , the total recovery time is split into two parts, the time between losing contact in the first cycle and the time the plate starts to move to perform the second cycle results in ε(2t t + t w ) using Equation (27). From this moment the plate starts to move while the material is still relaxing according to Equation (27). To determine the time t mc , the equation is solved, obtaining The second cycle is completed by repeating the steps taken to derive the first cycle, starting with Equation (21). The initial conditions are given by ε 0,2 = −t mcεt and t 0,2 = 2t t + t w + t mc . The analytical derivations allow to establish an expression of the TPA parameters in terms of the mechanical properties and test parameters.

Texture parameters
The stress in the material after the loading step in the first compression cycle σ (t t ) is given by Equation (21), from which the hardness results, making use of the fact that the test time t t = −ε t /ε t : To calculate the springiness, the times s 1 and s 2 indicated in Fig. 1 are required. The resulting springiness equals in which the time to make contact t mc is given by Equation (29). The areas A 1 and A 2 in Fig. 7 are equivalent to A 1 and A 2 in Fig. 1 and are calculated by integrating the stress as a function of time, whereby the stress is negative. The area A i of cycle i is split into a loading part A l,i and an unloading part A u,i , which are found by integrating Equations (21) and (22): Then the cohesiveness C is calculated as The adhesiveness A 3 is found by integrating the unloading curve, given by Equation (22), in the interval [t t + t 0 lc , t t + t lc ] and multiplying it with the initial cross-sectional area A 0 :

Constitutive model
The constitutive model that is used to describe an elasto-viscoplastic material is a specific case of the model presented by Jonkers et al. (2020) and is implemented in MSC.Marc using the user subroutine HYPELA2. To separate elastic and plastic deformation, the total deformation gradient tensor is multiplicatively decomposed (Lee 1969) in an elastic part F e and a plastic part F p : The Cauchy stress tensor that results from this deformation is calculated as with the shear modulus G and the bulk modulus K. The volume change ratio J and the isochoric elastic finger tensorB e are defined as The plastic rate of deformation tensorD p is calculated from the equivalent plastic shear strain rateγ eq p , the equivalent shear stress τ eq , and the deviatoric part of the elastic second Piola-Kirchhoff-like stress tensorŜ The deviatoric part of the elastic second Piola-Kirchhoff-like stress tensor is defined aŝ and the equivalent shear stress is given by The equivalent plastic strain rate is obtained from the equivalent shear stress aṡ with the reference shear rateγ 0 , the characteristic stress τ c , and the exponent m. The value m = 50 is chosen to ensure that the yield stress is approximately independent of the applied strain rate, andγ 0 = 10 −3 s −1 is chosen arbitrarily as the reference shear rate for τ c . The resulting uniaxial intrinsic material response is shown in Fig. 9. The material parameters that influence this response are the elastic moduli G and K and the characteristic stress τ c , which defines the maximum stress, that is, the yield stressσ y at the corresponding yield strainε y . Note that the material parametersσ y andε y are equivalent measures (i.e., defined positive), whereas in compression it holds thatσ = −σ y andε = −ε y at the yield point. Figure 3 reveals a small force increase when reaching the yield stress, whereas the true stress remains constant in the plastic regime. This is due to the increase of the cross-sectional area during plastic deformation. By introducing damage in the material the stress does not remain equal but is decreasing after reaching the yield stress. This effect is also known as strain softening. The softening regime is characterized by a damage parameter ξ , which decreases the initial, undamaged shear modulus G 0 and characteristic stress τ 0 of the material: In the undamaged case the parameters retain their initial values, that is, τ c = τ 0 and G = G 0 . The damage parameter, which is initially zero, evolves with the equivalent plastic strain rate according toξ The damage parameter d determines the rate at which the damage develops, and ξ ∞ defines the limit value of ξ . The influence of these parameters on the intrinsic material response is depicted in Fig. 9. The addition of damage in the material may have mesh-dependent results as a consequence if no precautions are taken. Here this effect is regularized by the viscoplastic response of the material, in particular, by the model parameter m in Equation (42).

Analytical approximation
The material is nonlinear, and a numerical method is required to obtain the exact response of the material to the prescribed deformation in the TPA test. However, by considering a uniaxial stress state and by making some approximations on the effective material response, analytical relations between the texture parameters and intrinsic mechanical properties can be recovered, which closely describe the actual material response. The specific form of Equation (42) implies that the plastic deformation is continuously increasing. Therefore the identified yield strainε y is not purely elastic, but has a small plastic contribution. By neglecting this plastic contribution and, in addition, the nonlinearity of the elastic regime, the yield strainε y can be approximated bȳ and the irrecoverable plastic strainε p byε p =ε t +ε y .
Moreover, all analytical relations are expressed in true stress and strain definitions because the purpose is to relate the TPA parameters to the intrinsic material parameters. If damage is included, then the plastic deformation induced by the applied strain affects the mechanical properties. The equivalent plastic shear strain γ eq p,t , resulting from ε t , is calculated as where the factor √ 3 converts the uniaxial strain to an equivalent shear strain. The damage parameter ξ is obtained by integrating Equation (44), which results in The damaged Young's modulus is calculated by combining Equations (20) and (43), resulting in in which ξ t is the induced damage, found by substituting γ eq p = γ eq p,t into Equation (48). Since here the time-dependence of the material behavior is neglected, which is a valid assumption for a high power-law exponent m in Equation (42), the only relevant test parameter is the applied strainε t . The standard values of t w = 1 s andε t = 2 s −1 , which correspond to a velocity of 20 mm/s (Vliet and Primo-Martín 2011) for a specimen length L 0 = 10 mm, are chosen.  Figure 10a shows the result of a typical TPA experiment, obtained by a finite element simulation of an elasto-viscoplastic material without damage, and Fig. 10b shows the corresponding stress-strain response. The applied strain exceeds the yield point, resulting in plastic deformation. The hardness is found from the deformed cross-sectional area and the stress in the material as

Texture parameters
Changes in cross-sectional area in the elastic regime result from the Poisson's ratio ν. The deformed area A is, neglecting nonlinearities in the elastic regime, given by In the plastic regime, the change in cross-sectional area is based on the incompressibility of plastic deformation, that is, the volume remains constant. The actual hardness H is calculated using the deformed cross-sectional area. If the area changes are small, then the hardness can be approximated from the initial cross-sectional area: If the deformation is purely elastic, then the springiness and cohesiveness are equal to one. In the case of plastic deformation, at times s 1 and s 2 in Fig. 1, the corresponding strain values at these time instances can be extracted from Fig. 10b. The strain corresponding to s 1 is the applied strainε t , and that corresponding to s 2 is the recovered elastic strainε y , resulting in the springiness Similarly, the areas A 1 and A 2 of Fig. 1 have equivalent energies per unit volume, which are found by extracting W 1 and W 2 in Fig. 10b. The required energy per unit volume for the loading of the first compression cycle is given by W 1 = 1 2ε yσy − (ε t +ε y )σ y . The recovered energy during unloading, which is also equal to the loading and unloading energy in the second compression cycle, is given by W 2 = 1 2ε yσy . From this the cohesiveness is calculated as Clearly, the springiness and cohesiveness are both measures of plastic deformation in the material and provide exactly the same information for the considered material. Figure 11a shows the stress-strain and force-time responses of a damaging elastoviscoplastic material for different values ofε t . Figure 11a reveals the interplay between the damage, causing the force to decrease and the increase in cross-sectional area, which result in a force increase. Once ξ = ξ ∞ is reached, only the drive for an increasing force prevails. Depending onε t , this increasing force may result in a second peak that is larger than the first peak. In that case the force of the first peak is known as the fracturability F r , and the force of the second peak corresponds to the hardness H .
The approximations of the hardness (52) and springiness (53) still hold. The calculation of cohesiveness requires the energy per unit volume in the softening regime W d , which is approximated by Analogously to the calculation of cohesiveness for the undamaged case, the cohesiveness of the damaging material is approximated by Due to the nearly time-independent material behavior, adhesion is not affecting any TPA parameter except for the adhesiveness. The time t a to reach the adhesion stressσ a is calcu- from which the adhesiveness is approximated by

Homogeneous deformation
The analytical derivations enable investigating the relation between texture parameters, mechanical properties, and test conditions. Analytical relations for the linear viscoelastic material are exact and show the effect of time-dependent properties. We compare the analytical results of the elasto-viscoplastic material with the numerical solution to investigate the accuracy of the approximations made. An overview of the two materials, including the analyzed mechanical properties and test parameters, is presented in Table 1.

Linear viscoelastic material
The hardness given by Equation (30) is determined by the product of the relaxation time and the strain rate λε t and the viscoelastic fraction φ of the material as presented in Fig. 12.
If λε t is high or if the fraction φ is low, then the material responds linear elastically and reaches the hardness value H e given by For very small values of λε t , the hardness takes the value (1 − φ)H e corresponding to the elastic response of the parallel spring with Young's modulus E ∞ . A combination of material parameters and test conditions that gives a normalized hardness H/H e close to one also results in springiness and cohesiveness values close to one because it is dominated by the elastic component. Figure 13a shows the effect of the material parameters on the springiness (Equation (31)) and cohesiveness (Equation (34)) for the adhesion-free case and constant test conditions ε t = −0.1,ε t = −2.0 s −1 , and t w = 1.0 s. Two characteristic regimes are visible. A small relaxation time means that the viscous component dominates the elastic response. Contact is lost in a short time period (i.e., t lc is small), but the material relaxes fast after losing contact. A large relaxation time implies that within the considered timescale, the elasticity is more important. The compressive strain decreases with the test rate for a longer time (i.e., t lc is closer to t t ). In comparison to a smaller relaxation time, the unconstrained material relaxation given by Equation (27) is slower, but this is compensated by the fast constrained unloading.
The interplay between the two regimes results in a minimum for the springiness and cohesiveness (as a function of relaxation time). These minimum values decrease for an increasing fraction φ, meaning that the viscoelastic response becomes more important, as depicted in Fig. 13a. The value of the total stiffness E does not affect the results, only the fraction φ matters. For a large relaxation time, φ has only a small influence on the TPA parameters, which corresponds to the elastic regime, also clearly visible in the hardness; see Fig. 12.
To show the effect of the test parameters in Fig. 13b, a material with relaxation time λ = 0.05 s and fraction φ = 0.995 is taken. The springiness and cohesiveness are functions of the test time t t = −ε t /ε t : a small test time corresponds to an elastic response and thus the texture parameter close to one. Larger values of t t imply a viscous response of the material. By increasing the test time the time to relax, while the plate is moving, becomes larger, which results in higher springiness and cohesiveness values. The same effect is obtained by increasing the time between the two compression cycles t w , as presented in Fig. 13b. Figure 14a shows the increasing adhesiveness A 3 given by (35) for an increasing adhesion stress σ a or a decreasing Young's modulus E. A determining factor in the adhesiveness is the time t a it takes to reach the adhesion stress. This time depends on the combination of σ a and E, as derived in Equation (23). Moreover, the adhesiveness depends on the re- laxation time λ, the fraction φ, and the strain rateε t . Consequently, for this material it is questionable whether the adhesiveness is a meaningful parameter, which is also discussed by Fiszman and Damasio (2000). Figure 14b shows the effect of adhesion on the springiness and cohesiveness. The effect is negligible when the adhesion stress is much smaller than the compressive stress at the applied strain ε t . This is illustrated by decreasing the deformation rate, which results in a lower maximum compressive stress. Subsequently, the springiness and cohesiveness are more affected by adhesion.

Elasto-viscoplastic material
The analytical solutions for the elasto-viscoplastic material without damage are compared to the nonlinear numerical results in Fig. 15 for different combinations of Young's modulus E, yield stressσ y , and applied strainε t . The validity of neglecting changes in cross-sectional area to calculate the hardness (52) depends onε t . The numerical solutions for a large compressive strain (ε t < −0.5) deviate significantly from the analytical solution, as depicted in Fig. 15a. A large compressive strain results in a large increase of the cross-sectional area, which is not taken into account in the analytical approximation. Next to this, neglecting the nonlinearity of the elastic regime causes a deviation between the analytical and numerical results if the yield strain ε y is large.
The analytical approximations for the springiness (53) and cohesiveness (54) are in good agreement with the numerical results. These texture parameters are only sensitive to the nonlinearities in the constitutive model and insensitive to variations in cross-sectional area. In the undamaged case, the springiness and cohesiveness provide the same information on the material response, in particular, a measure of the plastic, irrecoverable deformation in the material.
This similarity does not hold if damage, resulting in the characteristic strain softening regime, is added. Figure 16a compares the analytical and numerical results of the springiness for different values of ξ ∞ . The numerical results confirm that the damage parameters do not affect the springiness value. In contrast to the springiness, the cohesiveness depends on the damage parameters. Because the effect of changing d turns out to be negligible, Fig. 16b only shows the results of the cohesiveness for different values of the damage limit value ξ ∞ . All cohesiveness values are below the prediction of the undamaged material as a result of the softening, which becomes more pronounced for larger values of ξ ∞ . Whereas the springiness only provides information on the deformation in the sample, the cohesiveness is a measure of structural damage in the material due to this deformation.
The adhesiveness A 3 of an elasto-viscoplastic material is depicted in Fig. 17. The analytical relation (58) is in adequate agreement with the actual, numerically obtained adhesiveness. If there is no damage, then the adhesiveness is directly obtained from the adhesion stressσ a , strain rateε t , cross-sectional area A 0 , and Young's modulus E and does not depend on the applied strainε t or yield stressσ y . This is not the case if damage is added to the material. The decrease of the stiffness induced by damage results in a higher adhesiveness. In Fig. 17, this is shown by increasing the limit value ξ ∞ , which reveals the same trends for a material with a lower yield strain ε y =σ y /E or for a higher applied compressive strainε t .

Heterogeneous deformation
The effect of friction and surface roughness are analyzed for the linear viscoelastic material as the effect on the elasto-viscoplastic material is similar. The effect of friction turns out to be negligible. The hardness increases 1% due to friction, whereas the springiness and cohesiveness remain unaffected.
We investigate the effect of surface roughness for a varying relaxation time λ. The other material parameters (E = 20 MPa and φ = 0.995) and test parameters (ε t = −0.1,ε t = 2.0 s −1 , and t w = 0.2 s) are constant. The effect of roughness on the TPA tests with λ = 0.03 s is presented in Fig. 18a. The stiffness in the first loading cycle is smaller because the initial contact area is small. This effect is less pronounced in the second loading cycle as the contact surface is flattened in the first cycle. The relative deviation of a parameter x is given by and depicted for the hardness, springiness, and cohesiveness in Fig. 18b. The hardness is most affected, especially if λ is large, resulting in a more elastic response. The relative deviation of cohesiveness exceeds 15% for small relaxation times. However, the absolute deviation is below 0.06 in the whole interval and becomes significant only in the region where the cohesiveness is small, similar to the minimum cohesiveness value in Fig. 13a. The relative deviation of the springiness is below 3% and therefore negligible.

Conclusion
In this research, constitutive models are used to predict the effective mechanical response of the materials. However, an appropriate constitutive model is not necessary to establish relations between texture parameters and intrinsic mechanical properties. Careful experimental characterization of the mechanical behavior of a food material provides all the information that is essential to predict its texture properties. A double compression test, the TPA test, can be used to determine the mechanical characteristics of a particular food product, but the analysis should involve the mechanical properties instead of the currently used instrumental texture parameters. Analytical relations between instrumental texture parameters, test conditions, and mechanical properties have been derived, and analyzed by numerical modeling. Two types of representative materials were chosen and characterized by a suitable constitutive model that is used to simulate the texture profile analysis test.
For the time-dependent linear viscoelastic material, exact relations are obtained by solving the differential equation for all steps of the TPA test. The characteristic relaxation time of the material in combination with the applied deformation rate is the governing factor in the interplay between the elastic and viscous components of the material. Combining these two characteristic time measures with the time-dependent stiffness of the material provides the instrumental texture parameters of a viscoelastic material. Although the values of springiness and cohesiveness are not equal for a specific set of material and test parameters, the same dependencies were identified. Consequently, the springiness and cohesiveness qualitatively provide the same structural information on the material when plasticity is not taken into account.
The large-deformation response of an elasto-viscoplastic material is studied to include the effect of plastic deformation. By neglecting some nonlinearities in the constitutive model, simple analytical relations are found between texture and material properties. Clearly, the springiness and cohesiveness are a direct result of plasticity in the material and provide exactly the same qualitative information. The plastic deformation, resulting from the applied strain, is easily found from the stiffness and yield stress of the material. By adding damage to the material, that is, a reduction of mechanical properties, it is found that the springiness is independent of the damage parameters for the considered damage model. The springiness only provides information on the current strain in the material, which may be recoverable (viscoelastic) or irrecoverable (plastic). However, the cohesiveness decreases when the damage in the material is more severe and contains therefore relevant structural information.
Two nonideal contact situations that cause heterogeneous deformation are investigated. The effect of friction on the results of the TPA test is negligible. Surface roughness has a major influence on the hardness but a minor effect on the other texture parameters as it mainly involves the first compression cycle.
Using the maximum compressive stress instead of the maximum force as the definition of hardness, which is currently affected by the specimen size, is suggested. An accurate calculation of the compressive stress requires the deformed cross-sectional area, which however can be approximated by the initial cross-sectional area if volumetric changes in the material are large, which may be the case, for instance, for porous materials. The springiness and cohesiveness provide important structural information, which is something that is definitely relevant for sensory attributes. However, the relevance of the instrumental adhesiveness is questionable for the used materials, because it does not only depend on the adhesion between the material and the probe. Clearly, the instrumental texture parameters are a result of the mechanical characteristics of the food product. A better option is therefore to link the intrinsic mechanical properties directly to sensory attributes, as the dependence on test conditions is inherently taken into account.
Funding The authors acknowledge TNO (Netherlands Organisation for Applied Scientific Research) for the financial support.

Data availability
The data generated or analyzed in this study are available upon reasonable request.
Code availability Numerical simulations have been conducted using MSC.Marc 2013.

Conflicts of interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.