Damage Function of a Quasi-Brittle Material, Damage Rate, Acceleration and Jerk during Uniaxial Compression: Model and Application to Analysis of Trabecular Bone Tissue Destruction

: A diversity of quasi-brittle materials can be observed in various engineering structures and natural objects (rocks, frozen soil, concrete, ceramics, bones, etc.). In order to predict the condition and safety of these objects, a large number of studies aimed at analyzing the strength of quasi-brittle materials has been conducted and presented in publications. However, at the modeling level, the problem of estimating the rate and acceleration of destruction of a quasi-brittle material under loading remains relevant. The purpose of the study was to substantiate the function of damage to a quasi-brittle material under uniaxial compression, determine the rate, acceleration and jerk of the damage process, and also to apply the results obtained to predicting the destruction of trabecular bone tissue. In accordance with the purpose of the study, the basic concepts of fracture mechanics and standard methods of mathematical modeling were used. The proposed model is based on the application of the previously obtained differentiable damage function without parameters. The results of the study are presented in the form of plots and analytical relations for computing the rate, acceleration and jerk of the damage process. Examples are given. The predicted peak of the combined effect of rate, acceleration and jerk of the damage process are found to be of practical interest as an additional criterion for destruction. The simulation results agree with the experimental data known from the available literature.


Research Problem
There is a large class of natural and artificial materials that, under mechanical influences, exhibit to a greater or lesser extent the property of brittleness. Such quasi-brittle materials include rocks, concrete, ceramics, bones, etc. A variety of quasi-brittle materials can be observed in various engineering structures and natural objects. With a decrease in temperature and some other influences, initially non-brittle materials can become quasibrittle [1]. To avoid catastrophes, it is necessary to continue studying the mechanical behavior of quasi-brittle materials under various influences. Depending on the physical and mechanical properties of the material and the effects on it, the appropriate experimental research methods and mathematical models were used. The main focus in this study was given to the substantiation of the damage function of a quasi-brittle material, modeling the kinetics of destruction during uniaxial compression on the example of trabecular bone tissue. Trabecular bone is a multifunctional tissue that maintains a dynamic balance of structure and functions through repeated cycles of remodeling and adaptation to changes in physical activity and physiological needs [2][3][4]. The structure and functions of the trabecular bone tissue have been studied for hundreds of years [5], but a number of problems remain relevant [6]. The relevance of the topic and the practical significance of the presented work are determined by the fact that the solution of this problem is necessary to improve the understanding of the structure and function of the trabecular bone tissue, which will contribute to a more accurate prediction of the risk of bone fractures [7].
Within the framework of the abovementioned problem, the main issue in the presented research was to determine the speed, acceleration and jerk of the process of damage to a quasi-brittle material on the example of trabecular bone tissue under uniaxial compression. An analytical solution to this problem is proposed, which is obtained by standard methods of mathematical modeling of mechanical systems and is based on the use of the damage function justified in [8]. In addition to [8], a more detailed analysis of the damage function was carried out, which allowed us to obtain new analytical relations for determining the rate and other characteristics of the damage process of a quasi-brittle material. A method of applying the obtained results to predict damage to the trabecular bone tissue during uniaxial compression is proposed. Verification of simulation results was performed using experimental data known from the available literature.

Damage Models of Quasi-Brittle Materials
The approach to modeling the accumulation of damage which is used in this work is focused on the application of the methodology of fracture mechanics. The main stages of the evolution of fracture mechanics, starting with the classical work of Griffith (1920), are reflected in [9]. The development of damage models involves the use of certain information about the material that makes up the simulated object; in this case, it is the trabecular bone tissue. The initial building materials of bones are organic and inorganic substances. A functioning bone contains water, mineral and organic substances [10]. Approximately 50% of the dry mass of a bone is accounted for by inorganic material in the form of hydroxyapatite, which, by nature of its mechanical properties, belongs to brittle materials, but is not identical to hydroxyapatite which is contained in rocks [11]. The effect of hydroxyapatite as a component of bones is manifested in the fact that the load displacement diagrams in the tests of trabecular tissue samples for compression are similar (but not identical) to the diagrams obtained in the tests of quasi-brittle rocks; the difference is due to the influence of organic components in the bone material [11]. Accordingly, the trabecular bone tissue is considered in the available literature as a quasi-brittle material [12][13][14], for modeling the state of which some rock models can be adapted [15]. Figure 1 shows typical stress-strain plots for the models of an ideal brittle material and a quasi-brittle material. The diagram in Figure 1a is typical for ultra-high-strength concrete [16]. For a sufficiently small segment, the dependence is almost linear (Figure 1b), which is used in the mathematical description of physical models of destruction of quasi-brittle materials, as shown below in Section 1.3.
sented research was to determine the speed, acceleration and jerk of the proce to a quasi-brittle material on the example of trabecular bone tissue under u pression. An analytical solution to this problem is proposed, which is obtain ard methods of mathematical modeling of mechanical systems and is based the damage function justified in [8]. In addition to [8], a more detailed an damage function was carried out, which allowed us to obtain new analytica determining the rate and other characteristics of the damage process of a qua terial. A method of applying the obtained results to predict damage to the tra tissue during uniaxial compression is proposed. Verification of simulatio performed using experimental data known from the available literature.

Damage Models of Quasi-Brittle Materials
The approach to modeling the accumulation of damage which is used i focused on the application of the methodology of fracture mechanics. The m the evolution of fracture mechanics, starting with the classical work of Griff reflected in [9]. The development of damage models involves the use of mation about the material that makes up the simulated object; in this case, it ular bone tissue. The initial building materials of bones are organic and in stances. A functioning bone contains water, mineral and organic substances imately 50% of the dry mass of a bone is accounted for by inorganic materi of hydroxyapatite, which, by nature of its mechanical properties, belongs to rials, but is not identical to hydroxyapatite which is contained in rocks [11] hydroxyapatite as a component of bones is manifested in the fact that the l ment diagrams in the tests of trabecular tissue samples for compression ar not identical) to the diagrams obtained in the tests of quasi-brittle rocks; the due to the influence of organic components in the bone material [11]. Acc trabecular bone tissue is considered in the available literature as a quasi-br [12][13][14], for modeling the state of which some rock models can be adapted [ Figure 1 shows typical stress-strain plots for the models of an ideal br and a quasi-brittle material. The diagram in Figure 1a is typical for ultraconcrete [16]. For a sufficiently small segment, the dependence is almost l 1b), which is used in the mathematical description of physical models of d quasi-brittle materials, as shown below in Section 1.3.  In fracture mechanics, the state of a quasi-brittle sample is characterized by the proportion of the damaged (and, consequently, nonfunctional, i.e., disconnected) part of the cross-section. The intact part of the cross-section is called the effective area. The proportion of the damaged part of the cross-section is determined by the following ratio: where A 0 is the cross-sectional area in the initial (undamaged) state and the variable value is A-the area of the intact part of the cross-section, i.e., the effective area (0 < A ≤ A 0 ) [17].
The effective area A depends on the sample strain ε, so D is considered to be a function D = D(ε). If at the start of the test ε = 0, then A = A 0 and D = 0. If the axial strain ε increases, the function values D change from 0 to 1. Thus, the value D = 0 corresponds to the undamaged state and the value D = 1 corresponds to complete damage. The main issue in the damage modeling accumulation is the choice, adaptation or justification of the abovementioned damage function D. A number of approaches to the definition of this function are known. It is important to note that in the known research, the damage function was usually chosen at the beginning of the simulation, based on the general understanding of the regularities of the damage accumulation process, in other words, the damage function was postulated.
For example, in [18], using the test results for sandstone, marble and granite, it was shown that the dependence of the value D on the axial strain can be modeled by solving the logistic equation dD/dε = rD(1 − D). As a result, the damage function was obtained in the following form [18]: where ε is the strain, a and r are fitting parameters. An analog of Equation (2) was obtained in [19]. In the Lemaitre model (1987), damage is associated with degradation of Young's modulus [20]: where E 0 and E = (1 − D)E 0 are Young's modulus for the initial and damaged state, respectively. Equation (3) is also used in [14,21]. In the study by Chen et al. (2021), it was assumed that in the initial state, the modeled system consists of N intact elements. If the load increases, then the elements are gradually damaged and turned off from the function. Then, the damage under loading is defined as the ratio of the number of failures of elements (n) to the total number of elements in the initial state (N) [22]: Based on the assumption that the strength of the elements (F) obeys the Weibull distribution law, in [22], using Equation (4), a model of damage development with two parameters F 0 and m was obtained: Equations (1)-(5) model the same process, but differ in form. Differences in approaches to damage analysis can be considered to be a sign of continuation of the search for new models in this area. Some formally different models can be equivalent under certain conditions, as shown in the following subsection.

Example of the Equivalence of Two Approaches to Damage Determination
Taking into account the need for model validation, we paid attention to the consistency of the results obtained using different approaches. In this case, consistency confirms the adequacy of the models or their equivalence under certain conditions. For example, one of the abovementioned approaches to substantiating the damage function [22] is based on the initial assumption that the strength of the elements obeys the Weibull distribution law.
Using a different approach [8], it can be shown that for uniaxial compression, an analog of equation (5) can be obtained based on the assumption that the change in the effective area (d A) of the sample with the initial height H 0 and the initial cross-sectional area A 0 is proportional to the change in the elastic potential energy if the strain increases from ε to (ε + dε). In this case, the height of the sample varies from εH 0 to (ε + dε)H 0 . The height change is equal to H 0 dε. Accordingly, the change in volume is dV = AH 0 dε. The elastic potential energy per unit volume is w = ε 2 E/2. The change in the elastic potential energy is dW = wdV = (ε 2 E/2) AH 0 dε. Let us assume that the change in the effective area A by an infinitesimal amount d A is proportional to the change in the elastic potential energy by dW : The effective modulus of elasticity E is the same for all particles of the material and does not change during loading (which is explained below in Section 2). The value proportionality coefficient β is a constant within the framework of solving a specific problem and is determined by Equation (7): The value ε test extr corresponds to the peak (extreme value) of the apparent stress σ test extr = F test extr /A 0 and is determined experimentally as a result of standard tests for uniaxial compression. The values ε test extr and σ test extr correspond to the extremum points in Figure 1b. Using (7), we transformed Equation (6) to (8): where the minus sign means that the effective area A decreases with increasing strain ε. Let us divide both sides of Equation (8) by A 0 and denote Θ = A/A 0 . Then, Equation (8) can be rewritten in the dimensionless form (9): We integrated Equation (9): We found the integration constant C 1 as follows: if ε = 0, then A = A 0 and Θ = 1. Then, C 1 = 0. Thus, the function Θ = Θ(ε) is determined by Equation (10): Function Θ = A/A 0 , where 0 ≤ A ≤ A 0 , determines the fraction of undamaged particles and, therefore, is a function of the residual resource of the cross-sectional area. Effective cross-sectional area is A = A 0 Θ. Taking into account Equation (10): The share of damaged particles D = 1 − Θ is a damage measure and defined as the function D = D(ε): Symmetry 2021, 13, 1759

of 17
Formal comparison of damage models (5) and (12) leads to the condition of their equivalence (13): The plots of the remaining resource (10) and damage (12) functions are shown in Figure 2.
The plots of the remaining resource (10) and damage (12) functions are shown in Figure 2. Remark 1: Note that Equation (12) was justified without postulating the law of damage evolution at the development start. Nevertheless, using the logic of fracture mechanics, a model of damage accumulation of quasi-brittle materials under uniaxial compression was substantiated, in which the final Equation (12) corresponds to the Weibull distribution law, which is confirmed by the equivalence of Equations (5) [22] and (12) [8] under condition (13).
Remark 2: To assess the damage of a quasi-brittle material during uniaxial compression according to Equation (12), the results of standard tests and ( Figure 1b) are sufficient.

A Model in Which the Elastic Modulus Does Not Change under Loading
As it was noted in the comments to Equations (6) and (7), in this study, it was assumed that the effective modulus of elasticity, like Young's modulus, is the same for the material of all particles and does not change during loading. However, damage and disconnection of the weakest particles lead to a decrease in sample stiffness during loading step by step. As a basis, the natural assumption that strength of the particles is not the same and that with increasing load the weakest of the particles remaining undamaged are damaged and turned off from functioning was used. The proportion of undamaged particles ( ) was determined using Equation (10). The material was considered to be a mechanical system of particles in equilibrium under the action of external and internal forces. Schematically, the damage process is shown in Figure 3, where the strain and apparent stress, respectively, are = / and = / . The understanding of the damage process in accordance to Figure 3 was used in the development of model (6)- (12).

Remark 1.
Note that Equation (12) was justified without postulating the law of damage evolution at the development start. Nevertheless, using the logic of fracture mechanics, a model of damage accumulation of quasi-brittle materials under uniaxial compression was substantiated, in which the final Equation (12) corresponds to the Weibull distribution law, which is confirmed by the equivalence of Equations (5) [22] and (12) [8] under condition (14).

Remark 2.
To assess the damage of a quasi-brittle material during uniaxial compression according to Equation (12), the results of standard tests ε test extr and σ test extr ( Figure 1b) are sufficient.

A Model in Which the Elastic Modulus Does Not Change under Loading
As it was noted in the comments to Equations (6) and (7), in this study, it was assumed that the effective modulus of elasticity, like Young's modulus, is the same for the material of all particles and does not change during loading. However, damage and disconnection of the weakest particles lead to a decrease in sample stiffness during loading step by step. As a basis, the natural assumption that strength of the particles is not the same and that with increasing load the weakest of the particles remaining undamaged are damaged and turned off from functioning was used. The proportion of undamaged particles (Θ) was determined using Equation (10). The material was considered to be a mechanical system of particles in equilibrium under the action of external and internal forces. Schematically, the damage process is shown in Figure 3, where the strain and apparent stress, respectively, are ε = ∆L/L 0 and σ = F/A 0 . The understanding of the damage process in accordance to Figure 3 was used in the development of model (6)- (12). mmetry 2021, 13, x FOR PEER REVIEW Figure 3. Cyclic growth and drop of the load as a result of sequential micropart step, the weakest particle among the functioning particles is damaged and d creases by the value of the bearing capacity of that particle. Consequently, the the system of interacting particles corresponds to a lower load value. Then, the strength of the next weakest particle, etc. The number of particles gradually de tive cross-sectional area decreases, resulting in a decrease in the sample's stif capacity. The figure shows a simplified scheme without taking into account t the growth and drop off the load during particle damage take place at the m the material and a corresponding decrease in the cross-sectional area from crease in the sample's stiffness from to . The elastic modulus of th change. Figure 3 shows that at the pre-peak stage, the increase in the load e In the post-peak state, the increase in the load is less than its decrease plained by a decrease in the effective area (11), as a result of which the loses the ability to withstand a continuously increasing load.
As for the realism of the scheme in Figure 3, the following shou scheme corresponds to a real process, then there should be a saw-tooth perimental plot. Since the growth and drop of the load during the dest is carried out at the microlevel, this pattern can be fixed if the sensitivit ment is high enough. The saw-tooth pattern can be seen in the experime ment plot in [23], which confirms the adequacy of the scheme presente

External and Internal Damage Processes
In the framework of the manifestation of cause-and-effect relation cal tests, we deal with two processes, external and internal, in relation external process is characterized by the rate of change of the load autom by the test equipment in accordance with the specified program. The ch . Cyclic growth and drop of the load as a result of sequential microparticle damage. At each step, the weakest particle among the functioning particles is damaged and disabled. The load decreases by the value of the bearing capacity of that particle. Consequently, the equilibrium state of the system of interacting particles corresponds to a lower load value. Then, the load increases to the strength of the next weakest particle, etc. The number of particles gradually decreases, so the effective cross-sectional area decreases, resulting in a decrease in the sample's stiffness and its bearing capacity. The figure shows a simplified scheme without taking into account the scale. In practice, the growth and drop off the load during particle damage take place at the microlevel. Damage to the material and a corresponding decrease in the cross-sectional area from A 0 to A leads to a decrease in the sample's stiffness from EA 0 to E A. The elastic modulus of the material does not change. Figure 3 shows that at the pre-peak stage, the increase in the load exceeds its decrease. In the post-peak state, the increase in the load is less than its decrease. This feature is explained by a decrease in the effective area (11), as a result of which the sample gradually loses the ability to withstand a continuously increasing load.
As for the realism of the scheme in Figure 3, the following should be noted. If the scheme corresponds to a real process, then there should be a saw-tooth pattern in the experimental plot. Since the growth and drop of the load during the destruction of particles is carried out at the microlevel, this pattern can be fixed if the sensitivity of the test equipment is high enough. The saw-tooth pattern can be seen in the experimental load displacement plot in [23], which confirms the adequacy of the scheme presented in Figure 3.

External and Internal Damage Processes
In the framework of the manifestation of cause-and-effect relationships, in mechanical tests, we deal with two processes, external and internal, in relation to the sample. The external process is characterized by the rate of change of the load automatically generated by the test equipment in accordance with the specified program. The characteristics of this process were automatically measured and displayed as a load displacement plot (Figure 4). Thus, we obtained the primary data for continuing the study.
The function F = F(∆L) (Figure 4) can be represented analytically as an equation [8]:  The experimental data ( Figure 4) can be transformed into the characteristics of the internal process, which include the damage of the sample material D (12) and the dependence of the apparent (σ) and effective ( σ) stresses on strain ε [8]:  The experimental data ( Figure 4) can be transformed into the c internal process, which include the damage of the sample material pendence of the apparent ( ) and effective ( ) stresses on strain [8] = ( ( ) ) = / ≈ 1.396 Equation (17) can be rewritten in the form of Hooke's law for un = Taking into account (17) and (18), the effective modulus of elas using Equation (19): In experiments, only part of the damage is visible on the surface although more complete data on the structure and properties of the Equation (17) can be rewritten in the form of Hooke's law for uniaxial loading: Taking into account (17) and (18), the effective modulus of elasticity is determined using Equation (19): In experiments, only part of the damage is visible on the surface of the sample, and although more complete data on the structure and properties of the material can be obtained, for example, using tomography, experimental data alone are not sufficient to better understand the process of damage accumulation [24][25][26]. Additional data can be obtained using mathematical modeling [27]. Nevertheless, despite the rapid development of mathematical models predicting the accumulation of damage during loading of quasibrittle materials, a number of questions remain open [7]. The most pressing issues are related to the determination of the rate and other kinetic characteristics of damage to quasi-brittle materials, primarily trabecular bone tissue. The solution of these issues is necessary to improve the assessment of the risk of bone fractures [14].

Damage Rate, Acceleration and Jerk
To determine the damage rate, we used Equation (12), which models damage as a function of strain ε. The strain ε is the result of an external influence, characterizes an external process and depends on the time t: ε = ε(t). The rate of the external impact can be defined as the dependence of the dimensionless strain ε on time (v load,ε ) or as the dependence of displacement, expressed in meters, on time (v load ). The rate dimensions in the first and second cases are, respectively, 1/s and m/s. If the height of the sample is equal to L 0 (Figure 4), then the ratio between v load,ε and v load is determined by Equation (20): v load,ε = v load /L 0 (20) Let the loading rate of the sample and the rate of damage accumulation be determined, respectively, as v load,ε = dε/dt and v damage = dD/dt. The function D is determined using Equation (12). Then, the following equation is true: For further comparison of the results of our modeling with the experimental data known from the available literature, we assumed that the rate v ε is a constant value since the experimental data known from the available literature are obtained, as a rule, at a constant strain rate. Then, the acceleration of damage accumulation (a damage ) and the jerk (j damage ) were determined using the following equations: A jerk j damage (1/s 3 ) shows how quickly the acceleration changes a damage (1/s 2 ). A rapid change in acceleration means a rapid strain that can lead to the destruction of the material. Therefore, the predicted breakthrough can be considered to be one of the harbingers of damage and destruction. The impact of the jerk can be partially reduced by using shock-absorbing materials or devices [28].
Using Equation (12), we transformed Equations (21)-(23) to the form (24)-(26), respectively: Equation (24) predicts that the damage rate is proportional to the loading rate. The influence of the loading rate on the crack opening rate using an experimental method was investigated in [29], where it was shown that if the load rate increases from 17.1 mm/s to 1750 mm/s (i.e., by about 100 times), then the crack opening rate increases at the same ratio (from 0.02 m/s to 2 m/s). A similar dependence is predicted by Equation (24).
The plots of Equations (24)- (26) for the examples in Figure 2 are shown in Figure 5.

Example
Let us consider the application of Equations (10)- (26) to the modeling of a sam trabecular bone tissue. As the initial data, we used the results of experiments pu in [26]. In the experiments in [26], the samples had the shape of a cylinder with a of 8 mm and a diameter of 4 mm. The samples were loaded until destruction a = 0.005 mm/s that, according to (20), corresponds to ,ε = 0.000625 1 0.0625% strain per second). The experimental curve from [26] is shown in Figure  line

Example
Let us consider the application of Equations (10)- (26) to the modeling of a sample of trabecular bone tissue. As the initial data, we used the results of experiments published in [26]. In the experiments in [26], the samples had the shape of a cylinder with a height of 8 mm and a diameter of 4 mm. The samples were loaded until destruction at a rate v load = 0.005 mm/s that, according to (20), corresponds to v load,ε = 0.000625 1/s (or 0.0625% strain per second). The experimental curve from [26] is shown in Figure 6                  A comparison of the plots in Figures 6-9 leads to the following conclusions.
2. Figures 6-9 show the data corresponding to the download rate v load = 0.005 mm/s = 5·10 −6 m/s. This is a quasi-static load, in which the absolute value of acceleration (a damage ) at the failure point is a small value that can be ignored (a damage = 257·10 −6 s −1 ). Nevertheless, if the rate of external influence increases, for example, to 5 m/s (18 km/h), that is, by 10 6 times, then, theoretically, the acceleration (25) will increase by 10 12 times (then, a damage = 257·10 6 s −1 ), and this will inevitably lead to the bone fracture. The failure point is located near the negative peak of function a damage = a damage (ε) (Figure 8).

3.
The jerk (j damage ) shows how fast the acceleration changes. In the considered example, at the failure point, j damage = 19.108·10 −6 s −3 . This is a negligible value. However, if v load = 5 m/s, then the jerk (26) will grow 10 18 times and will be j damage = 19.108·10 12 s −3 . The failure point is located near the negative peak of function j damage = j damage (ε) (Figure 9). 4.
It is important to note that the fracture point is close to the peak values of the damage rate, as well as of the acceleration and jerk (Figures 8 and 9). This indicates that the destruction occurred with the greatest combined effect of the damage rate, acceleration and jerk. Thus, a quantitative assessment of the overall impact of the damage rate, as well as of the acceleration and jerk, can be taken into account in predicting the fracture risk.

The Damage Function and Its Derivatives
Even without reference to the example discussed above (Figures 6 and 7), engineering intuition suggests that the rate, acceleration and jerk of the damage process are somehow related to the strength of the quasi-brittle material. How can we obtain quantitative assessment of the influence of velocity and other kinetic characteristics of the destruction process of a quasi-brittle material on its strength under uniaxial compression? To solve this problem, we turned again to the physical sense of the derivatives of the damage function D (12). All components of this function have a clear physical meaning and can be determined as the result of standard uniaxial compression tests. The function D shows the change in the proportion of damaged material depending on the strain. It is important to note that Equation (12) is differentiable. The first derivative (v damage ) shows how quickly the proportion of the damaged material (D) increases depending on time t (21). The second derivative (a damage ) shows how fast the rate (v damage ) changes depending on time (22). The third derivative (jerk) shows how fast the second derivative (acceleration) changes with time, etc. Higher derivatives are not often used in engineering practice, so the terms for denoting their physical content still required consensus. Formally, all the derivatives of the damage function D (12) can be represented in the following form: The practical significance of Equations (27)- (30) lies in the possibility of their use in models for predicting destruction of a quasi-brittle material under uniaxial compression, which is mentioned in Section 5 and discussed in more detail in the next subsection.

Application of the Damage Function and Its Derivatives to the Prediction of Destruction
Let us consider the possible application of Equations (27)- (30) for predicting the destruction of a quasi-brittle material under uniaxial compression. The analysis of the plots of damage rate, acceleration and jerk in the abovementioned example (Figures 6 and 7) leads to confidence in the reality of such a possibility. From the physical point of view, it is logical to assume that the probability of material destruction is proportional to the combined influence of rate, acceleration and jerk. To take into account the joint influence of rate, acceleration and jerk, we transformed these functions to a dimensionless, normalized to one, form. To do this, we divided both parts of each of Equations (24)- (26), respectively, by max(abs(v damage )), max(abs(a damage )) and max(abs(j damage )). Thus, we obtained functions that we denoted as − v damage , − a damage and − j damage . Accordingly, the values of these functions were equal to the dimensionless estimates of the influence of rate, acceleration and jerk on material destruction. The plots of the functions are shown in Figure 10.
The practical significance of Equations (27)-(30) lies in the possibility of the models for predicting destruction of a quasi-brittle material under uniaxial comp which is mentioned in Section 5 and discussed in more detail in the next subsecti

Application of the Damage Function and Its Derivatives to the Prediction of Destruc
Let us consider the possible application of Equations (27)- (30) for predicting struction of a quasi-brittle material under uniaxial compression. The analysis of t of damage rate, acceleration and jerk in the abovementioned example (Figures leads to confidence in the reality of such a possibility. From the physical point of is logical to assume that the probability of material destruction is proportional to t bined influence of rate, acceleration and jerk. To take into account the joint influ rate, acceleration and jerk, we transformed these functions to a dimensionless, nor to one, form. To do this, we divided both parts of each of Equations (24)-(26), resp by ( abs( )), ( abs( ))and ( abs( )). Thus, we o functions that we denoted as ̄ , ̄ and ̄ . Accordingly, the v these functions were equal to the dimensionless estimates of the influence of rate, ation and jerk on material destruction. The plots of the functions are shown in Fig   Figure 10. Dimensionless normalized functions ̄ , ̄ and ̄ .
The total effect of rate, acceleration and jerk on material damage can be det as the average of three estimates: = ( )=sum(abs( ̄ ) + abs( ̄ ) + abs( ̄ ))/3 It is logical to assume that, theoretically, the maximum of Equation (31) corr to fracture. The graph of Equation (31) is shown in Figure 10 with a double gre Two extrema on this line are interesting: at = 3.3% and = 4.5%. The absciss of the extrema ( = 4.5%) almost coincides with the experimental value of the ( = 4.6%), which corresponds to the fracture point [26].
In Figure 10, it can be seen that failure of the sample corresponds to a chang sign of acceleration, that is, the compressive strain decrease. Consequently, the f force in the contact zone of the particles also decreases, which can lead to shea particles, for example, in concrete or rocks. As for the trabecular bone tissue, hy cally, the change in the sign of acceleration and jerk (Figure 8) may be associated appearance of local stretching of the trabecula, which does not contradict the expe of Turunen et al. [30], who found that, under loading, the structure of cancello tissue demonstrates an alternation of positive and negative local volumetric defo Thus, rate, acceleration and jerk (27)  The total effect of rate, acceleration and jerk on material damage can be determined as the average of three estimates: It is logical to assume that, theoretically, the maximum of Equation (31) corresponds to fracture. The graph of Equation (31) is shown in Figure 10 with a double green line. Two extrema on this line are interesting: at ε = 3.3% and ε = 4.5%. The abscissa of one of the extrema (ε = 4.5%) almost coincides with the experimental value of the abscissa (ε = 4.6%), which corresponds to the fracture point [26].
In Figure 10, it can be seen that failure of the sample corresponds to a change in the sign of acceleration, that is, the compressive strain decrease. Consequently, the frictional force in the contact zone of the particles also decreases, which can lead to shear of the particles, for example, in concrete or rocks. As for the trabecular bone tissue, hypothetically, the change in the sign of acceleration and jerk (Figure 8) may be associated with the appearance of local stretching of the trabecula, which does not contradict the experiments of Turunen et al. [30], who found that, under loading, the structure of cancellous bone tissue demonstrates an alternation of positive and negative local volumetric deformities. Thus, rate, acceleration and jerk (27)-(29) make it possible to obtain an adequate prognosis of destruction of the trabecular bone tissue. Maximum (31), as one of the criteria for failure under uniaxial loading, can be taken into account in further studies.
It should be noted that only the first, second and third derivatives are taken into account in Equation (31). More derivatives can be used to account for more factors. For this, using Equation (30), we rewrote Equation (31) in the general form (32): Equations (27)- (30) can be determined analytically or numerically. Equation (31) shows the simplest case of using rate, acceleration and jerk to predict fracture, that is, it is assumed that dimensionless estimates of the rate effect, acceleration and jerk equally contribute to material fracture. It would be logical and more realistic to introduce weighting factors. However, the question of justifying the weighting coefficients remains open. Therefore, the author does not offer Equation (31) as a replacement for the generally accepted criteria and estimates of damage, but only shows another possible application of the damage function and its derivatives.

Methodological Aspects of Modeling Trabecular Bone Tissue
As noted above, trabeculae are periodically renewed, as a result of which the age of various trabeculae and, consequently, their mechanical characteristics do not coincide. For example, according to [31], the strength of the cancellous bone tissue can range from 0.15 to 13.7 MPa. Due to the difference in the age of trabeculae and other factors, the real structure of trabeculae is heterogeneous, and each of its fragments is unique, as discussed in [32]. The heterogeneity of the trabecular bone tissue significantly complicates the modeling, which is necessary, first of all, to assess the risk of bone fractures.
One of the methods of modeling structures as mechanical analogs of trabecular bone tissue is the finite element method (FEM) [32,33]. This is a computer-based method, the practical application and capabilities of which positively correlate with the development of computer technologies. New FEM models, which take into account the inhomogeneity of trabeculae, anisotropy and nonlinear behavior of samples of trabecular bone tissue, allow us to study the functioning of the trabecular bone in more detail [34]. However, the limiting factor in nonlinear FEM models is the relatively high cost of their computer implementation [32,35]. Therefore, in some cases, the use of simplified FEM models continues, in which the inhomogeneity and asymmetry of the elastic characteristics of real bone tissue is ignored since it is assumed that the finite element material is homogeneous, isotropic and linearly elastic [36]. In simplified models, the elastic properties of the material are characterized by two constants (Young's modulus and Poisson's ratio). In more complex models, models of anisotropic materials are used, the elastic properties of which are characterized by a large number of constants-from three in the simplest case to 21 in the case of the most general type of anisotropy [37]. Nevertheless, the use of various models and a comparative analysis of the simulation results make it possible to overcome difficulties, the reasons for which are that the experiments contain errors and the models use simplifications, and therefore reasonable questions always arise, which, at the same time, motivate the search for new solutions [32,38,39].
An analysis of the available literature shows that a new and more complex model does not always guarantee a better result compared to a simple but physically adequate model. In this regard, in the context of this article, a recent study [34] is of scientific and practical interest, in which, in particular, it is showed that "heterogeneous anisotropic models accurately predicted the apparent elastic modulus but were no better than a simple homogeneous isotropic model" [34]. In other words, the values of the apparent elastic modulus calculated using two different models are almost the same. Such a coincidence confirms the adequacy of both the complex model and the simple model. Of course, the anisotropic model allows us to get a more detailed idea of the features of the sample behavior during loading. However, in some cases, it is advisable to use a simpler model.
In this case, the adequacy of a simple isotropic model can be explained with the following considerations. For example, if both models are adequate, then the experimental load displacement plots for the trabecular tissue sample should coincide with the theoretical plots for both the anisotropic model and the isotropic model.
The experimental load displacement plot contains the primary data on the behavior of the sample during loading. If the experimental load displacement plot coincides with the plot for both the simplified and anisotropic models, it means that the isotropic model directly or indirectly takes into account the same factors as the anisotropic model. Thus, in the study cited above [34], we found confirmation that it is not always necessary to strive for a formal complication of the model. In other words, the methodological guideline in this case may be the principle concerning the choice of theory, known as Ockham's razor: do not multiply entities beyond the need [40,41]. Essentially, the development of a model of damage to a quasi-brittle material for determining static and kinematic characteristics (10)-(26) during uniaxial compression is based on this principle. Namely, Section 1.3 was fulfilled, and characteristics (10)- (26) were obtained without postulating the damage function and without using the law of softening of a quasi-brittle material on the descending branch of the load displacement plot under uniaxial compression.
Let us pay attention to the manifestations of symmetry and asymmetry of the object of research and the process of its functioning. Uniaxial compression can be considered to be a special case of an axisymmetric problem of mechanics. The trabecular bone tissue fills the proximal part of the human femur. At the same time, the right and left femurs are symmetrical, but the external symmetry hides the asymmetry and uniqueness of the structure of the trabecular tissue [32].
The model proposed above can be used to analyze not only trabecular bone tissue, but also some other quasi-brittle materials. For example, the previous version of the proposed model (similar to the Furamura model [42,43]) was used for the analysis of low-and highstrength concrete (30 and 70 MPa) [44]. However, the details of the answer require separate consideration in another article.

Remark 3.
It should be noted that this article considers the quasi-static uniaxial compression of not too rigid materials, whereas in practice, there is three-dimensional dynamic loading of various materials. Accordingly, the application scope of the obtained results is limited.

Conclusions
Based on the results of the study presented above, the following conclusions can be made: 1.
The model of damage accumulation of a quasi-brittle material under uniaxial compression is substantiated, in which the final Equation (12) corresponds to the Weibull distribution law, which is confirmed by the equivalence of Equations (5) [22] and (12) [8] under condition (13).

2.
A model of the behavior of a quasi-brittle material under loading is described, in which the elastic modulus of the material does not change under loading at the prepeak and post-peak loading stages (despite the decrease in the stiffness of the sample with its gradual damage) (Section 2). 3.
The analytical relations for calculating the rate, acceleration and jerk of the damage process are justified (Section 3.2). 4.
The possibility of predicting material destruction using the damage function and its derivatives is investigated. It is established that the predicted peak of the combined effect of rate, acceleration and jerk of the damage process is of practical interest as an additional criterion of destruction (Section 5).

5.
Application limitations to the obtained results are noted (Section 6, Remark 3).
Funding: This research received no external funding.
Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Acknowledgments: The author highly appreciates the consultations of Professor Rudolf Meltzer in Petrozavodsk State University concerning the biomechanics of trabecular bone tissue. Special thanks should be expressed to each of the three anonymous reviewers for their valuable comments on the manuscript. The staff of the MDPI editorial office include English editing service deserve full gratitude.

Conflicts of Interest:
The author declares no conflict of interest.