Abstract
We study a hypoplastic model for soil and granular materials stemming from geomechanical engineering which further incorporates effects of degradation of the granular hardness, therefore allowing for the description of environmental weathering. The governing system is described by a nonlinear system of transcendental-differential equations for stress and strain rate, which is investigated with respect to its long-time dynamic. Under deviatoric stress control, two different solutions of the underlying, implicit differential equations are constructed analytically. The spherical components of stress and strain rate converge asymptotically to an attractor and lead to the sparsification of material states. Whereas under cyclic loading-unloading carried out in a numerical simulation, finite ratcheting of the deviatoric strain rate is observed in the form of a square spiral.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
Introduction
We study stress–strain rate constitutive relations describing granular materials like cohesionless soils or broken rocks within the hypoplastic theory proposed by Kolymbas [21] and developed further by [34, 36, 37]. Unlike hyper- and hypoelastic material laws, the hypoplastic response differs for loading and unloading, thus corresponding to inelastic materials [19]. On the other hand, the strain is not decomposed into elastic and plastic parts unlike classical elastoplastic concepts [3, 20]. The hypoplastic model is incrementally nonlinear and of the rate type. See other incrementally nonlinear constitutive equations in [4, 13], the rate-independence concept and hysteresis in [31, 32, 35]. Under cyclic loading, one can naturally obtain limit circles with barodesy [22, 23]. For mathematical modeling of granular and multiphase media, we cite [15,16,17, 28,29,30], and for well-posedness and stability analysis of the relevant differential equations, we refer to [11, 12, 14, 33].
Our study relies on a simplified version of the hypoplastic constitutive relation that was introduced independently by Bauer [5] and Gudehus [18]. In previous works [8, 9, 24, 26, 27], we omitted the pressure and density dependence of soil materials and considered the strain–stress law as a nonlinear differential equation for the stress under a given proportional strain rate, the so-called strain control. The existence of an analytical solution made it possible to describe explicitly stress paths obtained for the specific cases of monotonic compression, extension, and isochoric deformations. Recently, the case of unknown strain rate that should be derived from a given proportional stress, called stress control, was recovered within implicit differential equations in [10, 25] together with the pressure and density dependence of model parameters in the hypoplastic law. In such studies, the density of a granular material, which represents the distribution of pores and solid grains, is expressed with the help of the corresponding void ratio.
Our current mathematical modeling is based on the degradation of granular hardness suggested by Bauer [6, 7], which is suitable for the description of mechanical weathering of soils and granular materials under dry and wet states. This can be caused by long-time phenomena in various environmental events. We investigate the dynamic of the hypoplastic model subjected to the exponential degeneration for special cases of strain and stress control. Under the strain control, the spherical stress is constant. Whereas in the stress control case, the spherical components of stress and strain rate are described by implicit differential equations, which yield to two different solutions, and consequently, two degradation scenarios. It is shown that both solutions converge asymptotically to an attractor implying sparsification of the material state. When cyclic loading-unloading was applied, then finite ratcheting occurred, that is, a shift of the hysteresis loops was observed in a numerical simulation carried out by the standard 4th order Runge–Kutta method. The corresponding hysteresis curve of the deviatoric strain rate factor has, in this case, the form of a square spiral.
Modeling
For the reader’s convenience, the notations used throughout the paper are collected alphabetically in Table 1.
In this study, we consider a deformable granular body with constant volume and whose granular hardness, \(h_\textrm{s}\), can be subjected to degradation according to the following phenomenological equation found in [7]:
where the parameter c has the dimension of time and \(h_\textrm{s}^\infty\) denotes the final value of the degraded solid hardness. Therefore, the dependence of the granular hardness on the time t is expressed for prescribed \(0\le h_\textrm{s}^\infty <h_\textrm{s}^0\) by an exponential function
where \(h_\textrm{s}^0 =h_\textrm{s}(t_0)\) for some \(t_0\ge 0\), and \(h_\textrm{s}(t)\rightarrow h_\textrm{s}^\infty\) as \(t\rightarrow \infty\). In particular, if \(h_\textrm{s}^\infty =0\), then Eq. 1 follows a constant proportion
It is worth mentioning that the degradation of the solid hardness as described in Eq. 1 is reasonable for modeling the cases when the material is subjected, for example, to variable weather conditions. Indeed, in these cases, even in the absence of applied external forces, we may still observe a degradation of the material caused by progressive or instantaneous weathering; see [7].
To describe the constitutive stress–strain relation for hypoplastic granular materials, we consider an extension of the model from [10] that takes into account the degradation of the solid hardness. More precisely, considering the Cauchy stress \(\varvec{\sigma }\) to be a three-by-three matrix and assuming that its orthogonal decomposition after normalization is of the form
the constitutive response between the stress and the three-by-three matrix of the strain rate \(\dot{\varvec{\varepsilon }}\) can be expressed as follows
using the dyadic product of tensors, where \(a>0\) is the yield strength, weight \(\kappa \in (0,1)\), \(\textbf{I}\) is the three-by-three identity matrix, and \(\Vert \dot{\varvec{\varepsilon }}\Vert =\sqrt{\dot{\varvec{\varepsilon }}:\dot{\varvec{\varepsilon }}}\) stands for the Frobenius norm. Concerning the physical interpretation of the right-hand side of the constitutive relation Eq. 4, the two terms which are linear with respect to the strain rate \(\dot{\varvec{\varepsilon }}\) account for the hypoelastic behavior, whereas the term with the norm \(\Vert \dot{\varvec{\varepsilon }}\Vert\) describes the inelastic behavior and allows for different signs of the strain rate, therefore distinguishing the loading and unloading cycles. Note that the inelastic term is weighted by the density factor \(f_\textrm{d}>0\), while the stiffness factor, \(f_\textrm{s}<0\), influences both the hypoelastic and the inelastic processes of the system. Finally, a non-constant granular hardness \(h_\textrm{s}\) entering Eq. 4 allows degradation phenomena.
For simplicity, the governing relations are given in time \(t\ge t_0\ge 0\) and we omit the spatial dependence for all physical variables. Note that, unlike the model in [7], the constitutive relation Eq. 4 is rate dependent. This is due to the perturbation originated by the term accounting for the degradation which, as seen in Eq. 2, shows a strong dependence on the time. There is a number of studies devoted to the investigation of rate-dependent processes, and in the past few years, this theory became more popular in the engineering community for applications in modeling hysteresis in multi-functional materials, see, e.g., [1, 2].
For cohesionless granular materials only non-positive principal stresses and not all equal to zero
are physically relevant; this yields, in particular, that the pressure \(p =-(\sigma _1 +\sigma _2 +\sigma _3)/3\) is positive. Like in [7], the factors \(f_\textrm{d}\) and \(f_\textrm{s}\) depend on the void ratio e. Indeed, in our investigation, while the volume of the grains in the body is assumed to be a constant \(V_\textrm{s}\), the void space has variable volume \(V_\textrm{v}\), hence the void ratio e is considered as a state quantity
The bounds \(e_\textrm{d}>0\) and \(e_\textrm{i}>0\) depend proportionally on the mechanical pressure according to the following identities found in [6] (which are related to Eq. 15)
for fixed values \(0< e_\textrm{min} <e_\textrm{max}\), by means of the auxiliary variable
Notably, the variable y appears in the compression law proposed by Bauer with typically \(n\in (0,1)\) and is related to the point of inflection observed in experiments in the compression curve, see [7].
Now, observe that \(\textrm{tr} \dot{\varvec{\varepsilon }}\) is the relative volume increment and, as mentioned before, the solid volume \(V_\textrm{s}\) is constant during the evolution. Denoting by V the total volume, the physical consistency of the system relies on a volume balance equation which reads as follows:
which in terms of the void ratio e introduced in Eq. 6 can be interpreted as
where \(\delta =\textrm{tr}\varvec{\varepsilon }\) is the volume strain (dilatation). Further, the solution of Eq. 9 for prescribed \(\delta _0\) and \(e_0\) is
Note that Eq. 9 can be also understood as a mass balance equation. Indeed, recalling that the material’s density is expressed in terms of the mass m of the granular body by
we get the continuity equation
In order to obtain a more approachable formulation of the problem, let us consider the decomposition of the strain rate into deviatoric and spherical parts
By differentiating Eq. 3, we get
Recalling that \(p =-\textrm{tr}\varvec{\sigma }/3\), from Eq. 4, it follows the following scalar differential equation for the trace
Using the identities above, we rewrite Eq. 4 in terms of the deviatoric part
while by differentiating Eq. 8, we obtain
Besides accounting for the effects of degradation, in our study of the constitutive equation for stress–strain relation, we assume that the void ratio changes in time according to the equation
The identity above has been proposed in [7] for modeling under isotropic compression and it allows us to simplify the system by restoring constant density and stiffness factors in Eq. 13. Similar technique has been used in [25] for studying the long-term dynamic of such type of system. More precisely, for prescribed \(y_0\ge 0\) and \(e_0\in (e_\textrm{min}, e_\textrm{max}) \exp (- y_0)\), the positive solution to Eq. 14 is expressed in terms of the auxiliary variable y as
In turn, the density and stiffness factors are functions of the current void ratio
for a model parameter \(b >0\), and the critical void ratio \(e_c\) is described, as counterpart of Eq. 7, by the identity
with constant \(e_\textrm{crt}\in (e_\textrm{min}, e_\textrm{max})\). As a consequence of Eqs. 15 and 16 the dependence of the void ratio can be suppressed and, thanks to Eq. 7, both the density and the stiffness correspond to the following constant factors
In summary, the coupled system we want to investigate comprehends the algebraic Eq. 8, and the differential Eqs. 12 and 13, for unknowns \(\varvec{\varepsilon }^*, p, \hat{\varvec{\sigma }}^*\). Our investigation of system Eqs. 12 and 13 focuses on its behavior when subjected to periodically oscillating loading-unloading cycles, therefore its analysis needs to be carried out in the corresponding monotonicity intervals of the cycle. Since the system is positive 1-homogeneous, by choosing a loading parameter s(t) to be a monotone differentiable function in some interval \(t\in (t_0,t_1)\), assuming that the prime denotes the derivative with respect to s, from Eq. 12, we get by the chain rule
where \(\dot{s}\) stands for the derivative of s with respect to t; and analogously for Eq. 13. Thus, dividing by this term, we finally obtain the corresponding system of differential Eqs. 12 and 13 in terms of s. More precisely:
where \(s_\pm =1\) if \(\dot{s}>0\), and \(s_\pm =-1\) if \(\dot{s}<0\). With this in mind, we identify a loading cycle for s increasing, i.e., \(\dot{s}>0\), while an unloading cycle refers to s decreasing, \(\dot{s}<0\).
The governing system comprehends two equations: one tensor equation for the deviator, Eq. 19, and one scalar equation for the trace, Eq. 20, which is implicitly given in terms of the auxiliary variable y introduced in Eq. 8, where the granular hardness \(h_\textrm{s}\) represents an inhomogeneous coefficient given by Eq. 2. Note that both equations are implicit through the norm \(\Vert \varvec{\varepsilon }^\prime \Vert\). Besides, the system is under-determined for unknown \(p, (\varvec{\varepsilon }^*)^\prime , \hat{\varvec{\sigma }}^*\) and it can be interpreted in terms of the stress if the strain is given, and vice versa. More precisely, the strain-controlled case consists in choosing a fixed tensor \(\textbf{D}\) normalized by \(\textrm{tr} \textbf{D} =1\), and assuming that the strain tensor has the form \(\varvec{\varepsilon } =\delta (s(t)) \textbf{D}\) as presented in the “Strain control” section. Analogously, in the stress-controlled case, we assume that a fixed tensor \(\textbf{T}\) is normalized by \(\textrm{tr} \textbf{T} =1\), and consider the stress given by \(\varvec{\sigma } =-3 p(s(t)) \textbf{T}\) in the “Stress control” section. The study of the long-time behavior for periodical cycling loading-unloading under stress control is contained in the “Long-time behavior of the stress-controlled case” section.
As mentioned before, in our investigation, we will also take into account equations Eqs. 9 and 14 which describe the changes of the void ratio. Combining these two identities, together with Eq. 15, and considering the time transformation, we get
We can see that the monotonicity intervals of y and \(\delta\) coincide; therefore, it is natural to choose the loading parameter s so that \(\dot{s}>0\) for \(\dot{\delta }>0\), while \(\dot{s}<0\) for \(\dot{\delta }<0\). The opposite relation holds for the variable y, that is, \(\dot{s}>0\) for \(\dot{y}<0\), and \(\dot{s}<0\) for \(\dot{y}>0\). Since y is associated with the pressure \(p=-\textrm{tr} \varvec{\sigma }\) via identity Eq. 8, the rate \(\dot{y}<0\) corresponds to compressive loading, and \(\dot{y}>0\) to unloading related to the extension phase.
Throughout the paper, we consider the initial conditions at \(t=t_0\) corresponding to the loading parameter value \(s(t_0)\) prescribed by given \(e_0, p_0, \delta _0, (\varvec{\varepsilon }^*)^\prime _0, \hat{\varvec{\sigma }}^*_0\) and \(y_0 =(3 p_0/h_\textrm{s}^0)^n\).
Analysis
The solvability of the initial value problem in the general case is open. We study the special cases under strain and stress control conditions.
Strain control
For a prescribed three-by-three symmetric second-order tensor \(\textbf{D}\) such that \(\textrm{tr} \textbf{D} =1\), let us assume that
hence its deviatoric part reads \((\varvec{\varepsilon }^*)^\prime =\delta '\textbf{D}^*\), and
Recalling that the loading parameter s follows the monotonicity intervals of \(\delta\), from the identities above and Eq. 19, we get the following differential equation with respect to \(\hat{\varvec{\sigma }}^*\) coupled with \(\delta ^\prime\):
Similarly, from Eq. 20, we obtain
this together with Eq. 21 implies
Equation 24 is trivially satisfied if \(y^\prime =0\), that is, \(y=y_0\) is constant and
In this case, by Eq. 21, we have \(\delta ^\prime =0\) and consequently the differential Eq. 23 for the deviatoric stress corresponds to
hence its solution can be obtained explicitly and has the form
Note that, for \(\delta '=y'=0\), Eq. 22 indicates that there is no change of strain; therefore, in accordance with the identity above, the stress forces are due only to the degradation of the granular material.
The other possible solution of Eq. 24 is expressed by the transcendental equation for y coupled with \(\hat{\varvec{\sigma }}^*\), namely,
Although the solvability of Eq. 26 follows from the general theory of ordinary differential equations, there is no explicit formula for its solution. The investigation of such a problem can thus rely on numerical techniques and it is a question for a future work. The theoretical and numerical analysis of a related hypoplastic problem like Eq. 4 for the strain-controlled case has been carried out in [8, 9, 24, 26, 27], however, disregarding the degradation effects.
Stress control
For a prescribed three-by-three symmetric second-order tensor \(\textbf{T}\) with non-positive eigenvalues, and \(\textrm{tr} \textbf{T} =1\), let
such that Eq. 5 holds and \(\textrm{tr} \varvec{\sigma }<0\). Inserting Eq. 27 into Eqs. 19 and 20, dividing by \(f_\textrm{s}\) and \(a^2\), using \((\hat{\varvec{\sigma }}^*)^\prime =\textbf{0}\), we get the coupled system for \((\varvec{\varepsilon }^*)^\prime\) and \(\delta ^\prime\):
Note that, the first equation is implicit as the norm \(\Vert \varvec{\varepsilon }^\prime \Vert =\sqrt{\Vert (\varvec{\varepsilon }^*)^\prime \Vert ^2 +(\delta ^\prime )^2/3}\). On the other hand, recalling the definition of y, we can see that Eq. 29 gives an equation for the pressure p, while the assumption of Eq. 27 indicates the relation between the stress components and the pressure, \(\varvec{\sigma } =-3 p \textbf{T}\).
We will show that, under certain conditions, we can find an explicit solution for the system. The first step towards such an analytical solution is to derive an expression for the norm \(\Vert \varvec{\varepsilon }^\prime \Vert\).
Theorem 1
Let Eq. 27 hold and assume that
Then
with the positive discriminant
Proof
Recalling the decomposition of the strain rate Eq. 11 (with the corresponding time transformation), from Eq. 28 it follows the identity
This leads to a quadratic equation for \(\Vert \varvec{\varepsilon }^\prime \Vert\), more precisely:
whose discriminant is given by Eq. 32 and is positive thanks to Eq. 30. To see that Eq. 31 is indeed the positive solution of the equation above it is enough to observe that Eq. 30 implies \(\frac{f_\textrm{d}}{a} \Vert \textbf{T}^*\Vert ^2 \le \Vert \textbf{T}^*\Vert\), while \(\sqrt{\textrm{Disc}}\ge \Vert \textbf{T}^*\Vert |\delta ^\prime _\textrm{s}|\). \(\square\)
Theorem 1 allows us to write the deviatoric part of the strain rate, Eq. 28, in terms of the function \(\delta\) as follows
Using this and Eq. 21, Eq. 29 can also be reduced to a quadratic differential equation in \(y^\prime\) whose analysis is the content of the following theorem.
Theorem 2
(Analytical solution) If Eqs. 27 and 30 hold, then the auxiliary variable y in Eq. 8 satisfies the quadratic differential equation
with the inhomogeneous coefficients depending on y as follows
where for brevity
and \(c_0, B, C, D\) are constant factors given by
Assume that
then Eq. 35 has two solutions, \(y_+\) and \(y_-\), which satisfy respectively the following two differential equations:
Proof
Using Eqs. 31 and 34, from Eq. 29, we get the equation
Recalling that by Eq. 21\(\delta ^\prime =-y^\prime /(1+\frac{\exp (y-y_0)}{e_0})\), we rewrite the identity above using the notation in Eq. 37 to obtain
where \(\textrm{Disc}\), given by Eq. 32, corresponds to
Squaring Eq. 41, after gathering alike terms, we get the quadratic Eq. 35 for \(y^\prime\) with the corresponding coefficients \(c_2(y), c_1(y), c_0\) from Eqs. 36-37. Note that
after some simplification yields to
Therefore, assumption Eq. 38 ensures that the discriminant of Eq. 35 is positive. Further, Eq. 38 together with Eq. 30 implies \(c_0\ge 0\), hence the quadratic Eq. 35 has two solutions, Eqs. 39 and 40, and the proof is complete. \(\square\)
Finally, by solving the differential Eqs. 39 or 40 with respect to y and recalling the relation \(\delta ^\prime =-y^\prime /(1+\frac{\exp (y-y_0)}{e_0})\), we obtain an explicit formula for the strain rate via Eq. 34.
Long-time behavior of the stress-controlled case
In the stress-controlled case analyzed in Theorem 2, we can make the following statements about the long-time behavior of the model:
Theorem 3
(Attractor) Under the stress control Eq. 27 and assumption Eq. 30, as \(t\rightarrow \infty\) the system Eqs. 9, 14, 28, and 29 tends exponentially to a stationary state given by constant
Proof
Indeed, as \(t\rightarrow \infty\) in Eq. 2, we have \(h_\textrm{s}(t)\rightarrow h_\textrm{s}^\infty\) and \(\dot{h}_\textrm{s} =0\) in the limit. Then \(h^\prime _\textrm{s} =0\) and, by Theorem 1, the norm \(\Vert \varvec{\varepsilon }^\prime \Vert\) is given by Eq. 31 with \(\delta ^\prime _\textrm{s} =\delta ^\prime\). Moreover, for \(h^\prime _\textrm{s} =0\), the discriminant in Eq. 42 corresponds to \(\sqrt{\textrm{Disc}} =\sqrt{D} |\delta ^\prime |\) and we can factorize Eq. 41 as follows
This implies that, in the limit, y is constant. We claim that this convergence is exponential, and as a consequence of Eqs. 10 and 15, we conclude Eq. 44. Indeed, note that equations Eqs. 39 and 40 are of the form \(y^\prime =(h^\prime _\textrm{s}/ h_\textrm{s}) F(y)\) for some bounded function F. Observing that \(h^\prime _\textrm{s}/ h_\textrm{s}\) behaves as the exponential function \(\frac{\exp (-t)}{\gamma +\exp (-t)}\), with \(\gamma >0\) constant, we need to distinguish two cases: \(F(y_0)<0\) and \(F(y_0)>0\), for a prescribed \(y(0)=y_0\). The case \(F(y_0)=0\) is trivial. If \(F(y_0)<0\) the exponential decay is clear. On the other hand, if \(F(y_0)>0\), then the convergence is driven by the exponential term which multiplies F. In summary, in either case, one can show that y converges to a constant at an exponential rate which proves the result. \(\square\)
For an illustrative example, we consider the normalized stress matrix in Eq. 27 of the shear form:
such that the principal eigenvalues of \({\textbf {T}}\) are \(\{0.8, 0.2, 0\}\), and \(\Vert {\textbf {T}}^*\Vert \approx 0.588\). The material parameters are taken from [6, 7] as follows: \(h_\textrm{s}^\infty =78.5\) MPa, \(h_\textrm{s}^0=120\) MPa and \(c=4\) hours for the degradation in Eq. 1; the exponents \(n =0.82\) in Eq. 8, and \(\alpha =0.18\), \(\beta =1.05\), \(b=1\) in Eq. 16; the void ratio bounds \(e_\textrm{min}=0.1\) and \(e_\textrm{max}=0.3\) in Eq. 7, \(e_\textrm{crt}=0.24\) in Eq. 17; \(a=0.6\) and \(\kappa =0.6\) in the hypoplastic law of Eq. 4. The initial values at \(t_0=0\) are prescribed as \(e_0=0.2\), \(p_0=10\) MPa, and \(\delta _0=0\) such that the factors \(f_\textrm{d}\approx 1.561\) and \(f_\textrm{s}\approx -0.844\) in Eq. 16. These data yield the coefficient \(c_0\approx 0.136\) in Eq. 36 and satisfy the solvability condition of Eq. 38 of Theorem 2.
We discretize the system of equations based on equidistant time points. To compute solutions \(y=y_+\) to the first-order differential Eq. 39 and \(y=y_-\) to Eq. 40, the standard fourth-order Runge–Kutta method is applied with the uniform step of 1 (hour).
A monotone behavior as \(t\rightarrow \infty\) in the hypoplastic model is described by the unified factor \(s_\pm =\pm 1\) implying the loading parameter \(s(t)=s_\pm t\) versus time t. This means that \(s_\pm =1\) when \(\delta\) is increasing, and \(s_\pm =-1\) when \(\delta\) decreases. In this case, a numerical solution to the system Eqs. 30–40 is plotted in Fig. 1.
The plots from (a) to (f) show respectively: the granular hardness \(h_\textrm{s}\), auxiliary variable y, pressure p, void ratio e, volume strain \(\delta\), and the deviatoric strain rate factor \(k_\textrm{d}\) for \(t\in [0,50]\). As seen in the figures, the spherical components \(y, p, e, \delta\) coincide for the loading/unloading, and they demonstrate two different scenarios of degradation corresponding either \(y_+\) or \(y_-\). In contrast, the deviatoric component \(k_\textrm{d}\) is the same for \(y=y_+\) and \(y=y_-\), whereas two curves in the plot (f) distinguish either loading as \(s_\pm =1\) or unloading as \(s_\pm =-1\) in Eq. 34.
Cyclic loading-unloading
In this section, we aim to investigate the model behavior in the stress-controlled case when the loading-unloading cycles are given by a continuous function s(t) such that
with an increasing sequence of turning points \(t_{k+1}>t_k\) for \(k =0,1,2,\ldots\). For instance, starting with \(s(t_0)=0\), we set \(t_{k+1}=2k-1\) and the piecewise-linear periodic function
such that \(\dot{s}=1\) for loading, and \(\dot{s}=-1\) for unloading, which is illustrated in the plot (a) of Fig. 2 during 10 cycles.
After solving the hypoplastic system Eqs. 30–40 numerically as described before, in Fig. 2b, we depict the corresponding void ratio e versus pressure p, while the stress–strain curve presenting the trace of the stress tensor \(\textrm{tr}\varvec{\sigma } =-3p\) versus volume strain \(\delta\) is plotted in Fig. 2c. Recalling that assumption of Eq. 38 yields two solutions, if \(y_+\not = y_-\), we get distinct \(\delta (y_+)\not =\delta (y_-)\) in Eq. 21. As a consequence Fig. 2b and c describe two possible scenarios of the granular body sparsification corresponding to the two roots \(y_+\) and \(y_-\) of the quadratic differential Eq. 35. Indeed, in Fig. 2b, we can see that the curves identifying the void ratio in either case, \(e(y_+)\) and \(e(y_-)\), approach one another as the pressure increases. Similarly, Fig. 2c shows that at higher pressure (therefore, smaller values of \(-3p\)) the volume strain has basically the same curvature for both \(\delta (y_+)\) and \(\delta (y_-)\).
To study the long-time behavior of the deviatoric strain rate factor \(k_\textrm{d}\), we consider the process of loading-unloading with an increasing number of cycles. The deviatoric strain rate factor \(k_\textrm{d}\) during 2, 3, and 10 loading-unloading cycles is shown within the plots (a), (b), and (c) in Fig. 3, respectively. The dynamic presents a square-spiral curve, in which hysteresis loops are not closed and obey ratcheting. This limit cyclic behavior can be explained by the fact that \(h_s^\infty >0\), which is in accordance with the expected limiting behavior of the granular degradation of Eq. 2 which is rate dependent. Despite that, when the number of cycles increases, the loops converge to an attractor set in accordance with Theorem 3.
In order to highlight the role of degradation, we vary the creep time c in the exponential function Eq. 2 for the granular hardness \(h_\textrm{s}\). Compared to \(c =4\) hours from Fig. 3, increasing the creep time to \(c =20\) hours decreases the deviation of the solution and suppresses its hysteresis ratcheting as presented in Fig. 4.
From our numerical tests, we found that an increase in the value of c delays the degradation, leading to a kind of limiting cycle which can be observed, while the convergence to a stationary constant state becomes visible just with a higher number of cycles.
Data availability
This manuscript has no data from the repository.
References
M. Al Janaideh, P. Krejčí and G.A. Monteiro, Inverse rate-dependent Prandtl–Ishlinskii operators and applications, Appl. Math. (Prague) 68 (2023), 713–726.
M. Al Janaideh and M. Rakotondrabe, Precision motion control of a piezoelectric cantilever positioning system with rate-dependent hysteresis nonlinearities, Nonlinear Dyn. 104 (2021), 3385–3405.
B.D. Annin, V.A. Kovtunenko and V.M. Sadovskii, Variational and hemivariational inequalities in mechanics of elastoplastic, granular media, and quasibrittle cracks, in: “Analysis, Modelling, Optimization, and Numerical Techniques”, G.O. Tost and O. Vasilieva, editors, pp. 49–56, Springer Proc. Math. Stat. 121, 2015.
P.J. Armstrong and C.O. Frederick, A mathematical representation of the multiaxial Bauschinger effect, C.E.G.B. Report RD/B/N 731, 1966.
E. Bauer, Calibration of a comprehensive hypoplastic model for granular materials, Soils Found. 36 (1996), 13–26.
E. Bauer, Analysis of shear band bifurcation with a hypoplastic model for a pressure and density sensitive granular material, Mech Mater. 31 (1999), 597–609.
E. Bauer, Long-term behavior of coarse-grained rockfill material and their constitutive modeling, in: “Dam Engineering - Recent Advances in Design and Analysis”, Z. Fu and E. Bauer, editors, IntechOpen, 2021.
E. Bauer, V.A. Kovtunenko, P. Krejčí, N. Krenn, L. Siváková and A.V. Zubkova, Modified model for proportional loading and unloading of hypoplastic materials, in: “Extended Abstracts Spring 2018. Singularly Perturbed Systems, Multiscale Phenomena and Hysteresis: Theory and Applications”, A. Korobeinikov, M. Caubergh, T. Lázaro and J. Sardanyés, editors, Trends in Mathematics 11, pp. 201–210, Birkhäuser, Ham, 2019.
E. Bauer, V.A. Kovtunenko, P. Krejčí, N. Krenn, L. Siváková and A.V. Zubkova, On proportional deformation paths in hypoplasticity, Acta Mechanica 231 (2020), 1603–1619.
E. Bauer, V.A. Kovtunenko, P. Krejčí, G.A. Monteiro and J. Runcziková, Stress-controlled ratchetting in hypoplasticity, Acta Mechanica, 234 (2023), 4047–4093.
M. Brokate, Newton and Bouligand derivatives of the scalar play and stop operator, Math. Model. Nat. Phenom. 15 (2020), 51.
M. Brokate and P. Krejčí, Wellposedness of kinematic hardening models in elastoplasticity, RAIRO Modél. Math. Anal. Numér. 32 (1998), 177–209.
D. Caillerie and R. Chambon, Existence and uniqueness for rate problems of geomechanics, Revue Française de Génie Civil 8 (2004), 537–561.
M. Eleuteri, C. Gavioli and J. Kopfová, Fatigue and phase transition in an oscillating elastoplastic beam, Math. Model. Nat. Phenom. 15 (2020), 41.
K. Fellner and V.A. Kovtunenko, A singularly perturbed nonlinear Poisson–Boltzmann equation: uniform and super-asymptotic expansions, Math. Meth. Appl. Sci. 38 (2015), 3575–3586.
K. Fellner and V.A. Kovtunenko, A discontinuous Poisson–Boltzmann equation with interfacial transfer: homogenisation and residual error estimate, Appl Anal. 95 (2016), 2661–2682.
J.R. González Granada and V.A. Kovtunenko, Entropy method for generalized Poisson–Nernst–Planck equations, Anal. Math. Phys. 8 (2018), 603–619.
G. Gudehus, A comprehensive constitutive equation for granular materials, Soils Found. 36 (1996), 1–12.
J. Hron, J. Málek and K.R. Rajagopal, Simple flows of fluids with pressure dependent viscosities, Proc. R. Soc. A 457 (2001), 1603–1622.
A.M. Khludnev and V.A. Kovtunenko, Analysis of Cracks in Solids, WIT-Press, Southampton, Boston, 2000.
D. Kolymbas, Introduction to Hypoplasticity, A.A. Balkema, Rotterdam, 2000.
D. Kolymbas, A Primer on Theoretical Soil Mechanics, Cambridge Univ. Press, 2022.
D. Kolymbas and G. Medicus, Genealogy of hypoplasticity and barodesy, Int. J. Numer. Anal. Methods Geomech. 40 (2016), 2532–2550.
V.A. Kovtunenko, E. Bauer, J. Eliaš, P. Krejčí, G.A. Monteiro and L. Straková (Siváková), Cyclic behavior of simple models in hypoplasticity and plasticity with nonlinear kinematic hardening, J. Sib. Fed. Univ. - Math. Phys. 14 (2021), 756–767.
V.A. Kovtunenko, J. Eliaš, P. Krejčí, G.A. Monteiro and J. Runcziková, Stress-controlled hysteresis and long-time dynamics of implicit differential equations arising in hypoplasticity, Arch. Math. (Brno) 59 (2023), 275–286.
V.A. Kovtunenko, P. Krejčí, E. Bauer, L. Siváková and A.V. Zubkova, On Lyapunov stability in hypoplasticity, in: “Proc. Equadiff 2017 Conference”, K. Mikula, D. Ševčovič and J. Urbán, editors, pp. 107–116, Slovak University of Technology, Bratislava, 2017.
V.A. Kovtunenko, P. Krejčí, N. Krenn, E. Bauer, L. Siváková and A.V. Zubkova, On feasibility of rate-independent stress paths under proportional deformations within hypoplastic constitutive model for granular materials, Mathematical Models in Engineering 5 (2019), 119–126.
V.A. Kovtunenko and A.V. Zubkova, Mathematical modeling of a discontinuous solution of the generalized Poisson–Nernst–Planck problem in a two-phase medium, Kinet. Relat. Mod. 11 (2018), 119–135.
V.A. Kovtunenko and A.V. Zubkova, Homogenization of the generalized Poisson–Nernst–Planck problem in a two-phase medium: correctors and estimates, Appl Anal. 100 (2021a), 253–274.
V.A. Kovtunenko and A.V. Zubkova, Existence and two-scale convergence of the generalised Poisson–Nernst–Planck problem with non-linear interface conditions, Eur. J. Appl. Math. 32 (2021b), 683–710.
M.A. Krasnosel’skiǐ and A.V. Pokrovskiǐ, Systems with Hysteresis, Springer, Berlin, Heidelberg, 1989.
P. Krejčí, Hysteresis, Convexity and Dissipation in Hyperbolic Equations, Gakkotosho, Tokyo, 1996.
H. Lamba, P. Krejčí and D. Rachinskii, The global stability of a class of history-dependent macroeconomic models, Math. Model. Nat. Phenom. 15 (2020), 49.
D. Mašín, Modelling of Soil Behaviour with Hypoplasticity: Another Approach to Soil Constitutive Modelling, Springer Nature, Switzerland, 2019.
A. Mielke and T. Roubíček, Rate-Independent Systems: Theory and Application, Springer, New York, 2015.
A. Niemunis and I. Herle, Hypoplastic model for cohesionless soils with elastic strain range, Mech. Cohes.-Frict. Mat. 2 (1997), 279–299.
W. Wu, E. Bauer and D. Kolymbas, Hypoplastic constitutive model with critical state for granular materials, Mech. Mater. 23 (1996), 45–69 .
Acknowledgements
The authors thank Erich Bauer for the discussion and the anonymous referees for their suggestions to improve the manuscript.
Funding
Open access funding provided by University of Graz. The work is supported by the OeAD Scientific & Technological Cooperation (WTZ CZ18/2020 and MULT 06/2023) financed by the Austrian Federal Ministry of Science, Research and Economy (BMWFW), by the Czech Ministry of Education, Youth and Sports (MŠMT) in the frame of the Danuble multilateral project 8X23001. JR also acknowledges the support by the Grant Agency of the Czech Technical University in Prague, Project No. SGS23/092/OHK1/2T/11. VAK acknowledges the financial support by the University of Graz. The Institute of Mathematics of the Czech Academy of Sciences is supported by RVO:67985840.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare no competing interest.
Rights and permissions
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/.
About this article
Cite this article
Kovtunenko, V.A., Krejčí, P., Monteiro, G.A. et al. HYSTERESIS OF IMPLICIT EQUATIONS IN HYPOPLASTICITY FOR SOIL MATERIALS WITH GRANULAR HARDNESS DEGRADATION. J Math Sci 280, 453–467 (2024). https://doi.org/10.1007/s10958-024-07089-x
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10958-024-07089-x