Mechanothermodynamic Entropy and Analysis of Damage State of Complex Systems

Abstract: Mechanics from its side and thermodynamics from its side consider evolution of complex systems, including the Universe. Created classical thermodynamic theory of evolution has one important drawback since it predicts an inevitable heat death of the Universe which is unlikely to take place according to the modern perceptions. The attempts to create a generalized theory of evolution in mechanics were unsuccessful since mechanical equations do not discriminate between future and past. It is natural that the union of mechanics and thermodynamics was difficult to realize since they are based on different methodology. We make an attempt to propose a generalized theory of evolution which is based on the concept of tribo-fatigue entropy. Essence of the proposed approach is that tribo-fatigue entropy is determined by the processes of damageability conditioned by thermodynamic and mechanical effects causing to the change of states of any systems. Law of entropy increase is formulated analytically in the general form. Mechanothermodynamical function is constructed for specific case of fatigue damage of materials due to variation of temperature from 3 K to 0.8 of melting temperature basing on the analysis of 136 experimental results.


Introduction
As is known, any scientific discipline serves and aims at understanding and describing these or those regularities and features of certain phenomena, situations, events caused by the existence of some real or thinkable objects that manifest themselves their specific properties [1,2]. Based on the considerations that the study of a new object, as a rule, generates a new scientific discipline as applied to Mechanics, the hierarchy of objects can be built. Successive complexity of objects studied by mechanics is shown in works [1][2][3][4]. The next new object is multiphase-mechanothermodynamics system. Methods of mechanics alone are insufficient for its study as well as methods of thermodynamics alone.
These methods and models addressing coupled problems of both stress-strain states and energy states of the complex systems working under mechanical and thermodynamic loading are discussed in well-known monographs [5][6][7][8][9][10][11]. The concepts of entropy and damage are important in order to construct models of the mechanothermodynamical system.
Basic ideas about mechanical behavior of materials under the fracture process are discussed in [9]. The work [12] addresses the main aspects of Damage Mechanics as the branch of Fracture Mechanics and some of its applications. Fundamentals of Physical Mesomechanics of heterogeneous media that develops at the boundary of Continuum Mechanics, Physics of Plasticity and Strength of Materials and studies the stressed and damaged material at linked micro-, meso-and macro levels are given in [13].
Constitutive relations for strain-induced damage in terms of thermodynamic consistency and also applications of Continuum Damage Mechanics to failures of mechanical and civil engineering v i , i " 1, 2, 3 (1) where σ ij are stresses; ρ is density; f i are volumetric forces; v i are velocities.
The law of conservation of mechanical energy for the continuum of the volume V, with regard to the repeated index summation rule is obtained by multiplying scalar Equation (1) by the velocity vector v i : The right hand-side of Equation (2) is the change in the kinetic energy K of the continuum of volume V: The equation for mechnical energy of continuum [6] is derived basing on the known transformations with regard to Gauss-Ostrogradsky's theorem: or: where ε ij is the strain rate; Π is the continuum surface; the l are the direction cosines at the continuum surface; δU/dt is the power of internal forces; δA/dt is the power of internal surfaces and volumetric forces. The symbol δ in expression (4) is used to underline that the increment in the general case cannot be an accurate differential.
In the thermomechanical statement the rate of change in internal energy U [6] is usually given by the integral: where u " lim ∆mÑ0 up∆mq ∆m is the specific internal energy (internal energy density) of the elementary volume of the mass ∆m.
The rate of heat transfer to the continuum is expressed as follows: where c i is the characteristic of the heat flux per unit area of continuum surface per unit time due to heat conduction; z is the constant of heat radiation per unit mass per unit time.
The law of change in the energy of thermomechanical continuum then assumes the form: Entropy 2016, 18, 268 4 of 34 In (7), the transformation of surface integrals into volume integrals allows the local form of the energy equation to be obtained: If the scalar product of Equation (1) and the velocity vector v i is subtracted from Equation (8), then the following form of the local energy equation is obtained: where dq is the heat flux per unit mass. According to expression (9), the rate of change in the internal energy is equal to the sum of the stress power and the heat flux to the continuum.
As applied to the thermodynamic system, two characteristic functions of its state are defined: absolute temperature T and entropy S that can be interpreted as the characteristic of ordering (or chaotic state) of the thermodynamic system. Usually it is assumed that the entropy possesses the property of additivity, i.e.,: S " In continuum mechanics [6,7], the specific entropy per unit mass is considered: The specific entropy increment ds can be due to the interaction with the environment (the increment ds (e) ) or inside the system itself (the increment ds (i) ) [6,7]: ds " ds peq`d s piq .
The increment ds (i) is equal to zero in reversible processes and is greater than zero in irreversible processes.
If the heat flux per unit mass is expressed through dq, then for reversible processes the increment will be: Tds " dq.
According to the second law of thermodynamics, the rate of change in the total entropy S of the continuum of the volume V cannot be smaller than the sum of the heat flux through the volume boundary and the entropy produced inside the volume by external sources (Clausius-Duhem's inequality) [6,7]: d dt where e is the power of local external entropy sources per unit mass. The equality in Formula (14) is valid for reversible processes and the inequality-for irreversible processes.
Transforming the surface integral into the volume integral in expression (14) can yield the relation for the rate of the internal entropy production per unit mass: Entropy 2016, 18, 268 5 of 34 In continuum mechanics, it is assumed that the stress tensor can be decomposed into two parts: the conservative part σ pCq ij for reversible processes (elastic deformation, liquid pressure) and the dissipative part σ pDq ij for irreversible processes (plastic deformation, liquid viscous stresses): The expression for the rate of change in energy (9) can then be presented in the following form: If it is assumed that relation (13) is valid for irreversible processes, then the rate of the total entropy production is: Expression (18) for the rate of a total change of local entropy in the elementary volume of the continuum can be very convenient in practice.
In view of assumption (10) on the entropy additivity the sum in (18) can be supplemented by other terms considering the internal entropy production in the liquid (gas) volume due to different mechanisms. Similarly, for the continuum volume dV, for example, internal chemical processes can be considered [5]: If dV is considered not as a finite volume but as an elementary volume of continuum, then changes in its specific energy and entropy based on Equation (17), Equations (19) and (20) can be written in the following differential form: ds " where n k is the number of mols per unit mass. For the continuum of the volume V, expressions (21) and (22) on the basis of relations (5) and (11) will assume the form: dS " The introduction of the chemical entropy component (the last terms in (21)-(23)) allowed one not only to obtain a more complete picture of the continuum state, but also to describe self-organization processes that result in initiating stable structures when the heat flux to the continuum is increased.
The above-presented known models of energy and entropy states of continuum (17)- (24), being rather general, do not nevertheless permit one to satisfactorily describe some processes occurring in such a continuum as a deformable solid. Convenient representation of the additivity of energy and entropy components (11) used, for example, for modeling elastic deformation is not suitable for the description of non-linear processes. The available models do not also allow for the entropy growth due to damageability of solids as a specific characteristic of changes in the structural organization. According to the tribo-fatigue concepts [1][2][3][4][67][68][69], the damageability is interpreted as any irreversible change in structure, continuity, shape, etc. of a deformable solid that leads to its limiting state (separation into pieces under mechanical fatigue, occurrence of the pitting of critical density or sizes under contact fatigue, limiting wear under sliding fatigue). Although, for example, the elasticity limit is implicitly taken into account during plasticity modeling, the damageability, for example, during mechanical or contact fatigue occurs in the conditions of linear elastic deformation and requires a particular approach for its description with regard to limiting fatigue characteristics of material. The above drawbacks are overcome in the approach below.

Main Principles
According to [2,4,67], the mechanothermodynamical (MTD) system in the general case represents the thermodynamic continuum with solids distributed (scattered) within it, interacting with each other and with the continuum. Consider its fragment of limited size ΩpX, Y, Zq shown in Figure 1.
The continuum has a temperature θ and a chemical composition Ch. Here are two interacting solid elements (A and B) that can move relatively to each other at the contact zone S(x, y, z). Arbitrary mechanical loads applied to one of them (for example, to element A) in x, y, z coordinate system are transformed into the internal transverse forces Q x , Q y , Q z , longitudinal forces N x , N y , N z and also into the bending moments M x , M y , M z . Element B is pressed to element A by the loads that are transformed into the distributed normal pressure p(x, y) and the tangential tractions q(x, y). The origin of the coordinates is placed at the point of original contact O of the two elements (prior to deformation). It is easy to see that the elements A and B together form the Tribo-Fatigue system [4] which could be reduced to the friction pair [2] in the absence of internal forces (N i " 0, Q i " 0, M i " 0, i = x, y, z). Thus, the Tribo-Fatigue system is the friction pair in which at least one of the elements perceives non-contact loads and, consequently, undergoes volumetric deformation. This representation of the MTD system has an advantage that the analysis of the states of a solid and the components of a system can adopt the appropriate solutions known in mechanics of deformable solid, in contact mechanics, in mechanics of tribo-fatigue systems (tribo-fatigue) and in tribology.
Our main task is to describe the energy state of the MTD system under the action of mechanical and thermodynamic loads with regard to the environmental influence.
The energy state of any system is very interesting in itself. However, as applied to the MTD system it is very important to study its damage and, as a result, to study the conditions of reaching the limiting state (fatigue fracture, wear etc.). Of special interest is the analysis of translimiting or supercritical conditions [2].  The origin of the coordinates is placed at the point of original contact О of the two elements (prior to deformation). It is easy to see that the elements А and В together form the Tribo-Fatigue system [4] which could be reduced to the friction pair [2] in the absence of internal forces ( x, y, z). Thus, the Tribo-Fatigue system is the friction pair in which at least one of the elements perceives non-contact loads and, consequently, undergoes volumetric deformation. This representation of the MTD system has an advantage that the analysis of the states of a solid and the components of a system can The main ideas, which are the fundamentals of the given theory, can be formulated with regard to [2][3][4] as follows: I. Due to the fact that the elements of the MTD system are subject to the loads of different nature-mechanical, thermal and electrochemical, the traditional analysis of their damage and limiting state under the action only of mechanical stresses or strains [70][71][72][73][74][75][76][77][78] etc. can be the basis for research. However it appears insufficient and, as a result, is ineffective. This means that there is a need to analyze MTD system states using more general energy concepts.
II. Considering that the damage of MTD system solids is determined by mechanical, thermodynamic [79][80][81][82][83] etc. and electrochemical loads, it is needed to introduce the generalized representation of its complex damage that is caused by these loads acting at a time. Call such damage any irreversible changes in shape, size, volume, mass, composition, structure, continuity and, as a result, physical and mechanical properties of the system elements. This means the corresponding changes in the functions of the system as the integrity.
III. The onset and development of complex damage is mainly determined by means of four particular phenomena: mechanical fatigue, friction, wear, thermodynamic and electrochemical processes. These phenomena are called particular in the sense that each of them can be realized as independent and separate. This leads to the corresponding energy state and damage in terms of particular (separate) criteria.
IV. In the general case, all these particular phenomena and processes in the MTD system appear simultaneously and within one area. The states of such a system are then caused not by one of any mentioned phenomena but by their joint (collective) development and, consequently, by their interaction.
V. If the physical state of the MTD system is described by means of all input energy u ř , then its  IX. The energy u 0 is considered to be a fundamental constant for a given material. It shouldn't depend on testing conditions, input energy types, damage mechanisms.
X. Dangerous volume is the 3D area V ij Ă V 0 of the deformable solid (V 0 is its working volume) with the critical state of material at all its points, of which it consists.
XI. In the general case, the limiting (critical) state of the MDT system is reached not due to a simple growth of effective energy components and, hence, due to the accumulation of irreversible damages that are caused by separate actions (loads of different nature), but as a result of their dialectical interaction, whose direction is characterized by the development of spontaneous phenomena of hardening-softening of materials in the given operating or testing conditions. In such a way, taking into consideration function (26), the hypothesis of the limiting (critical) state of the MTD system can be represented in the following general form: Φpu where m k , k = 1, 2, . . . , are some characteristic properties (hardening-softening) of contacting materials, Λ kzlzn ¡ 1 are the functions (parameters) of dialectic interactions of effective energies (irreversible damages) that are caused by loads of different nature. This means that at Λ k > 1, the damage increase is realized, at Λ l < 1-its decrease, and at Λ n = 1-its stable development. XII. Taking into consideration Item III, from the physical viewpoint, hypothesis (27) should be multi-criterion, i.e., it should describe not only the states of the system as the integrity but its separate elements in terms of different criteria of performance loss (wear, fatigue damage, pitting, corrosion damage, thermal damage, etc.). In particular cases, it is possible to reach the corresponding limiting (critical) states in terms of one or two, three or several criteria at a time.
XIII. Reaching the limiting state: u means the full loss of the integrity of the MTD system, i.e., of all its functions. At the same time damageability of its elements: reach the critical value: ψ e f f u´ψσpchq , ψ τpchq , ψ Tpchq , Λ kzlzn , m k¯" 1.
XIV. If t = t 0 is the time of origination of the system and T ' is the time of reaching the limiting state, then the failure time of its functions corresponds to the relative lifetime (longevity) t{T ' " 1. But the system lifetime T˚as the material object is longer than its lifetime as the functional integrity (T˚ąą T ' ) since at the time moment t ą T ' the process of degradation-disintegration is realized by forming a great number of remains, pieces, fragments, etc. This process develops under the influence of not only possible mechanical loads but mainly of the environment-up to the system death as the material object at the time moment t = T˚. The system death means its complete disintegration into an infinitely large number of ultimately small particles (for example, atoms). The translimiting existence of the system as a gradually disintegrating material object can then be described by the following conditions: where d ψ is the average size of disintegrating particles and the natural relation ψ e f f u pd ψ q should exist between ψ Σ and d ψ . Then the condition for the system death is: XV. The particles of the "old system" disintegration are not destructed but are spent for the formation and growth of a number of "new systems". This is the essence of the MTD system evolution hysteresis.
Finally it should be noted that not irreversible but effective energy causing damages production is considered in the present work. Determination of parameters denoting effective (damaging) part of energies of different nature (due to stresses, friction, temperature) is explained in Section 4 "Energy Theory of Damage and Limiting States".
Our idea is that failure depends on fundamental physical property of the material-critical energy u 0 which is independent of the conditions of deformation, shapes and sizes of objects. Limiting energy u 0 is considered as the initial activation energy of the disintegration process (critical energy) which approximately corresponds both to the sublimation heat for metals and crystals with ionic bonds and to the activation energy of thermal destruction for polymers. It is independent of conditions and ways of reaching the limiting state. Longevity of a gear or a shaft in the given conditions of operation would be different but determined by fundamental parameter u 0 .
On the other hand effective energy is of course dependent on all these conditions and external loads and corresponding stress-strain state. Limiting state occurs when effective energy reaches u 0 . In cases of different sizes, shapes, number of cycles and stress concentrators effective energy is different. Therefore it may need different conditions and time for effective energy to reach fundamental material constant u 0 .
In our case entropy is defined as the ratio between the effective energy and critical energy u 0 (see Section 5 "Mechanothermodynamical states"). Therefore the proposed formulation of entropy simultaneously contains acting and critical energy. The limiting state occurs when this ratio is equal to 1. Latter studies could show effectiveness of this in certain cases.

Energy Theory of Damage and Limiting States
Let us specify function (25). To determine the effective energy, consider the work of internal forces in the elementary volume dV of tribo-fatigue systems (in Figure 1). In the general case, the differential of the work of the internal forces and the temperature dT ř can be written with regard to the rule of disclosing the biscalar product of the stress and strain tensors σ and ε: du " σ ij¨¨d ε ij`k dT Σ "¨σ xx σ xy σ xz σ yx σ yy σ yz σ zx σ zy σ zz‹ ‚¨¨¨d ε xx dε xy dε xz dε yx dε yy dε yz dε zx dγ zy dε zz‹ ‚`kdTΣ " " σ xx dε xx`σyy dε yy`σzz dε zz`σxy dε xy`σxz dε xz`σyz dε yz`k dT Σ ; (34) here k is the Boltzmann constant.
We proceed from the idea that in the general case, according to [2,4], the main role in forming wear-fatigue damage is played by the normal and shear stresses that cause the processes of shear (due to friction) and tear (due to tension-compression).
In this case, it is reasonable to divide the tensor σ into two parts: σ τ is the tensor of friction-shear stresses, or, briefly, the shear tensor and σ n is the tensor of normal stresses (tension-compression), or, briefly, the tear tensor. So in (28), the tear part σ n and shear part σ τ of the tensor σ will be set as: According to items III and IV, the tensors σ and ε should be represented as follows: Here the stress and strain tensors with the superscript V are caused by the action of volume loads (the general cases of 3D bending, torsion, tension-compression) and those with the superscript W-by the contact interaction of the system elements. Expression (35) with regard to (36) can be given as follows: In the case of the linear relationship between the stresses and strains, expression (36) will assume the form: and (37) will be as follows: From (40) it is seen that the tear part σ n of the tensor σ is the sum of the tear parts of the tensors at the volume strain σ pVq n and the surface load (friction) σ pWq n , whereas the shear part σ τ is the sum of the shear parts σ pVq τ and σ pWq τ . This means the vital difference of the generalized approach to the construction of the criterion for the limiting state of the MTD system. From total energy (40), its effective part is separated according to Items V and VIII with regard to [2,3]. To do this, introduce the coefficients of appropriate dimensions A n (V), A τ (V) and A T (V) that determine the fraction of the absorbed energy: or : du where Λ M\T (V) and Λ τ\n (V) are the functions of interaction between energies of different nature. The subscript τ\n means that the function Λ describes the interaction between the shear (τ) and tear (σ) components of the effective energy, and the subscript M\T means that the function Λ describes the interaction between the mechanical (M) and thermal (T) parts of the effective energy. That fact that the coefficients A can be, generally speaking, different for different points of the volume V, enables one to take into account the inhomogeneity of environment. Taking into consideration (42), criteria (27) can be specified with no regard to the environmental influence: In the case of the linear relationship between the stresses and strains, expressions (41) and (42) will be as follows: or : With regard to expression (36), criterion (43) can be represented as follows: When the time effects should be taken into consideration, criterion (46) will assume the form: Thus, expression (45) is the concretization of function (25) and formula (46) is the concretization of criterion (27) for that case when the environmental influence is not taken into account.
Criterion (27) in the form of (46) and (47) says: when the sum of interacting effective energy components caused by the action of force, frictional, and thermal (thermodynamic) loads reach the critical (limiting) quantity u 0 , the limiting (or critical) state of the MTD system (of both separate elements of the system and the system as the integrity) is realized. Physically, this state is determined by many and different damages.
The fundamental character of the parameter u 0 has been mentioned above. According to [82][83][84][85][86][87][88][89][90][91][92][93][94][95], the parameter u 0 will be interpreted as the initial activation energy of the disintegration process. It is shown that the quantity u 0 approximately corresponds both to the sublimation heat for metals and crystals with ionic bonds and to the activation energy of thermal destruction for polymers: On the other hand, the quantity u 0 is determined as the activation energy for mechanical fracture: In such a way, the energy U 0 can be considered to be the material constant: Taking into consideration the physical-mechanical and thermodynamic representations of the processes of damage and fracture [84,85,87], write down (48) in the following form: where s k is the reduction coefficient, σ th is the theoretical strength, E is the elasticity modulus, C a is the atom heat capacity, α V is the thermal expansion of the volume, k is the Boltzmann constant, T S is the melting point, θ D is the Debye temperature, h is the Planck constant. According to (49), it can be taken approximately [84]: where ε˚« 0.6 is the limiting strain of the interatomic bond. Calculations according to (50) are not difficult. Methods of experimental determination of the quantity u 0 have also been developed [85]. From equality (49) it follows that u 0 is the activation energy of a given material, which is by the order of magnitude equal to 1-10 eV per one particle or molecule (~10 2 -10 3 kJ/mol), i.e., the value that is close to the energy of interatomic bond rupture in the solid [88]. Its level doesn't depend on how the rupture is -mechanically, thermally or by their simultaneous action. In [85], it is possible to find the tables containing the u 0 values for different materials.
From (49) it is possible to find the thermomechanical constant of the material [2]: The constant θ σ characterizes the strength loss per 1 K. Criterion (46) is written in the absolute values of physical parameters-the values of the effective and critical energy components. This criterion can be easily made dimensionless by diving it by the quantity u 0 . Then it can be represented in terms of irreversible (effective) damage: or in detailed form: According to (54), it is also possible to determine particular energy measures of damage: due to the effective energies of different nature that are determined by the force (the subscript n), frictional (the subscript τ), and thermodynamic (the subscript T) loads, respectively. Now criterion (52) can be written in dimensionless form: "´ψ According to (58), the limiting state of the MTD system is reached when the sum of interacting damages (0 < ψ < 1) at mechanical and thermodynamic loads is equal to 1. Criterion (46) in form (58) is convenient because all damage measures are dimensionless and are within the same range 0 ď ψ ď 1.
Since numerous and infinite actions, as well as the interaction effects of physical damages of many types (dislocations, vacancies, non-elastic deformations, etc.) cannot be described and predicted exactly, when analyzing the MTD system, one introduces the concept of the interaction of dangerous volumes [2] that contain a real complex of damages (defects generated by the action of the corresponding fields of stresses (strains)). According to the statistical model of the deformable solid with the dangerous volume [88], such a volume should depend on the geometric parameters of the solid responsible for its working volume V 0 , on the distribution function parameters p(σ´1) and p(σ) of the durability limit σ´1 and the effective stresses σ considering both the effective stress probabilities P and γ 0 and gradients G σ : Here ϑ V describes how the limiting durability is influenced by the body shape and the scheme of body loading during fatigue tests.
Thus, the dangerous volume can serve as the equivalent of the complex of damages, as its value is proportional, in particular to the level of effective stresses and, hence, to the number (concentration) of defects (damages).
The boundary between the volumes of dangerous and safe, as it follows from expression (59), is generally blurred and probabilistic in nature. As the damage probability P of the solid increases, the dangerous volume V Pγ is growing. At a given value of P the volume can vary depending on the confidence probability γ 0 . It means that at P = const: if γ min ď γ 0 ď γ max . Here γ min , γ max form the permissible range. If it is accepted that γ 0 = const, then the dangerous volume will have a single value associated with the damage probability P. Scattered damage within the dangerous volume is characteristic not only for the so-called smooth bodies but also for the elements with the structural stress concentrators [88].
Thus, if in the uniaxial stress state, the stress distribution σ (x, y, z) in x, y, z coordinates is known, then the dangerous volume should be calculated by the formula: where σ´1 min being the lower boundary of the range of the durability limit σ´1 statistical distribution is such that if σ´1 ă σ´1 min , then P = 0. From expression (61) it follows that the generalized condition for fatigue fracture is of the form: with some probability P under the confidence probability γ 0 . If: then the fatigue damage cannot occur physically (because in this case, σ ă σ´1 min ); hence, (63) is the generalized condition of non-fracture.
The methods for calculation of dangerous volumes V ij for friction pairs and tribo-fatigue systems are developed similar to (59): ( 64) and outlined in [4,[90][91][92]. Here σ pV, Wq lim is the limiting stress based on the assigned criterion of damage and fracture.
Further, the following dimensionless characteristics of damage can be introduced: integral energy damage within the dangerous volume: and the average energy damage (at each point of the dangerous volume) : The accumulation of energy damage in time within the dangerous volume is described by the formulas: Ψ e f f u pV, tq " Having used (63)-(68), the MTD system damage can be described and analyzed using the most general representations-the energy concepts with regard to the influence of numerous and different factors taken into account by (59), including the scale effect, i.e., the changes in the size and shape (mass) of system elements.
According to [2,95], the function Λ k\l\n for damage interactions in the MTD system is determined by the parameters ρ of the effective energy ratio (Λ k\l\n ): The quantities Λ calculated by (69) describe the influence of the load parameter ratio on the character and direction of interaction of irreversible damages [2][3][4]. If Λ > 1, then the system is self-softening because at the balance of hardening-softening phenomena, softening processes are dominant. If Λ < 1, the system is self-hardening, because at the balance of hardening-softening phenomena, hardening processes are dominant. At Λ = 1, the system is stable-the spontaneous hardening-softening phenomena are at balance in it. The general analysis of damage interactions in the MTD systems will be given on account of its fundamental importance in a separate paper.
After criterion (27) has been basically formalized, the action of electrochemical loads (damages) should be taken into consideration according to Item VII. It should be said at once that it is difficult to do in the strict mechanothermodynamical statement: electrochemical reactions are very diverse and complex, when the environment interacts with a deformable solid body, and are insufficiently studied. That's why, the approach proposed in [2,3] is adopted: the simplification is introduced, according to which the damage of solids in the environment is determined by corrosion-electrochemical processes. In addition, the hypothesis is put forward, following to which the effective energy of corrosion-electrochemical damage is proportional to the square of the corrosion speed, i.e., : If according to Item VII, 0 ď D ch ď 1 is the parameter of corrosion-electrochemical damage of the body, then based on [2,4,92], criterion (26) with regard to its shape will be as follows: Here v ch is the corrosion speed in this environment, v chpTq , v chpσq , v chpτq is the corrosion speed in the same environment under thermal, force, and friction actions, respectively; b e are the coefficients responsible for corrosive erosion processes; m V(‚) are the parameters responsible for the electrochemical activity of materials at force (the subscript σ), friction (the subscript τ), and thermodynamic (the subscript T) loads, wherein m V(‚) = 2/A ch and the parameter A ch ¡ 1.
In [95], other methods for assessment of the parameter D ch can be found. As seen, Equation (72) is the specification of criterion (27). According to this criterion, the limiting state of the MTD system is reached when the sum of dialectically interacting irreversible damages at force, friction, and thermodynamic loads (including electrochemical damage when acted upon by stress, friction, temperature) becomes equal to unity.
Further, consider the particular case when in (46) it is assumed that Firstly, the stress state is caused by volume deformation, for which all components of the stress tensor, except for one component σ (one-dimensional tension-compression, pure bending), can be neglected. Secondly, the stress state is caused by surface friction, for which all components of the stress tensor, except for one component τ w , can be neglected. Then (40) assumes the following form: or in accordance with (72): (77) is the simplest form of the energy criterion of the limiting state that is nevertheless of great practical importance [2].
If there is no electrochemical influence of the environment (D ch = 0), then: Equation (78) is the simplest form of the energy criterion of the limiting state, which, nevertheless, is of great practical importance [2,92,95]. It serves particularly for the development of methods of assessing the parameters a T , a σ , a τ . In fact, at Λ MzT " Λ τzn " 1, the boundary conditions are the following: T Σ " 0 , τ w " 0 : a n σ 2 d " u 0 , a n " u 0 {σ 2 d ; where σ d , τ d are the force and friction limiting stresses as T Ñ 0. These are called the limits of (mechanical) destruction, T d is the destruction temperature (when σ = 0, τ w = 0) or the thermal destruction limit. The effective ("dangerous") part of total energy of strain can also be determined from the following physical considerations. It shall be assumed that the strain energy flux U generated in the material sample during its cyclic strain (ε " ε max sinωt) in the homogeneous (linear) stress state is to a certain extent similar to the light flux. In fact, it is continuously excited when the loading cycle is repeated with the speed ω " 1{λ. This enables one to consider it as a wave (with the length λ). Some part of the energy u generated in such a way can be absorbed by material atoms and structural formations, followed by damage of material. Denote the absorbed part of the energy by u eff . Then the generated energy u is equal to: u " u e f f`u cons , where u cons is the non-absorbed part (it is called the conservative part) of the generated energy u.
If the analogy of light and energy strain is justified, then the strain absorption law may be similar to Bouguer's light absorption law. Consequently, the equation relating the energy u cons passed through the material strain volume V and the generated energy u is of the form: or, in accordance with Lambert, in differential form: Here as in Bourguer-Lambert's equation, the coefficient χ ε independent of u is the energy absorption parameter.
Taking into account (81) into (80), the strain energy absorption law is obtained: and hence if u = 0 or V = 0 then u eff = 0. If V Ñ 8 it appears that according to (81) u cons = u, i.e., all input energy is dissipated within such a volume. ‚ changes in (internal) Tribo-Fatigue entropy (wear-fatigue damage mechanics [2]).
It should be noted that approach (83) can also be extended to the case of friction, since any indenter drives a strain wave upstream in the thin surface layer of the solid to which indenter is pressed to. Energy absorption parameter in this case will be χ γ where the subscript γ denotes the shear strain. Similarly, heat absorption in the deformable solid body can also be considered. Finally, the problem of strain energy absorption in the non-uniform (including complex) stress state can be easily solved by putting the dangerous volume V " V Pγ into (81)- (83).
It should be noted that although criterion (78) is special, it is fundamental and general in nature. Its general nature is caused by the fact that in this case, all four particular phenomena responsible for the MTD system state (in the statement simplified in terms of the stress-strain state) are taken into account (in accordance with Item III). Its fundamental nature is that here, as in complete solution (46), Λ n\τ takes into account the interaction of effective mechanical energy components caused by friction τ w and normal σ stresses, whereas Λ M\T takes into account the interaction of the thermal and mechanical components of the effective energy. The effective energy thermal component is determined by the variations of the total temperature T ř = T 2´T1 in the bodies contact zone caused by all heat sources, including the heat released during mechanical (spatial and surface) strain, structural changes, etc.

Mechanothermodynamical States
Within the framework of mechanothermodynamics a special approach is being developed to assess the entropy in terms of a generalized energy state. Following this approach and formula (77), out of the total energy (specific) due to some particular loads (force, temperature, etc.), its effective part directly spent for the damage production is defined by the experimentally found coefficients A l in formulas (41), (42), (77)  : u e f f l " A l u l , where the u l are the specific internal energies at tear (u n ), shear (u τ ), thermal action (u T ). The total specific energy of an elementary volume and a rate of its change are then given as: Moreover, the Λ-functions are used to take into consideration a complex (nonadditive) character of interactions between effective energies of different nature expressed by formula (42). This allows the total effective energy of the system to be assessed: .
where the Λ α are the possible combinations of interaction of effective energies (irreversible damages).
The specific feature of Λ-functions is such that: and hence: u Thus, using coefficients A l and Λ-functions it is possible to assess energy interaction due to different-nature loads. Such interaction can cause both a sharp growth and a substantial decrease of effective energy, resulting in damages and limiting states, as compared to the one calculated by the ordinary additivity model of type (17): The total effective energy of volume V and its accumulation in time with regard to formula (87) are of the form: and : The principal moment of the mechanothermodynamical model is the account of the limiting state (limits of plasticity, strength, fatigue, etc.) according to Item XIII: where u 0 is the limiting density of the internal energy interpreted as the initial activation energy of the disintegration process.
A relationship between the current state (mechanical, thermomechanical, energy) of an elementary volume of a solid (medium) and its limiting state enables one to construct the parameter of local energy damageability: dimensionless: Local energy damageability (94) or (95) is most general out of the damageability parameters constructed in terms of different mechanical (thermomechanical) states ϕ : ; eqv is the equivalent mechanical state; the ij are the components of the tensor φ; the i are the main components of the tensor φ; S and D ij and the sphere and deviator parts of the tensor φ; n and τ are the normal and tangential components of the tensor φ; int is the intensity φ; u is the specific potential strain energy (internal energy density); the indices at u mean: n u and τ u are the specific potential strain energy at tension-compression and shear; e f f u is the effective specific potential strain energy. Integral damageability measures can be built on the basis of local measures (96) with the use of the model of a deformable solid with a dangerous volume (64)-(68) [89].
By the dangerous volume is understood the spatial region of a loaded solid, at each point of which the value of local damageability is smaller than the limiting one : or: Dangeous volumes are calculated by the following general formula: The integral condition of damageability of a solid or a system can be written in the form: where V 0 is the working volume of the solid. To analyze at a time dangerous volumes local damageability distributed within them, the function of damageability of unit volume is introduced: The function of damageability of the entire volume V will then be as: The simplest functions of damageability accumulation in time for unit volume and the enrire volume will be of the following form, respectively: The indices of volume-mean damageability: and its accumulation in time can be used: The analysis of formulas (94), (100), (102) comes to the conclusion that conceptually, they are related to the concept of entropy as a difference (or relations) between two states (configurations) of a system, the degree of its organization (chaotic state). As applied to damageability, such states are current and limiting. Now using local energy damageability (94), construct specific (per unit mass) tribo-fatigue entropy (up to a constant): or: where A ψ is the dimensional parameter (J¨mol´1¨K´1).
On the basis of the expressions for entropy (18), as well as of formulas (85), (86) the local entropy and the rate of its change within an elementary volume will be: and: From formulas (108) and (109), it is seen that unlike the thermomechanical model, the state indices of the mechanothermodynamical system u and s are not equivalent. This is due to the fact that the calculation of the specific tribo-fatigue entropy s TF by formula (106) is supplemented by the limiting state in the form of the limiting density of the internal energy u 0 .
The tribo-fatigue entropy S TF is calculated not within the entire volume V, but only within its damageable part, i.e., within the energy effective dangerous volume V e f f u : On the basis of formulas (11), (106) and (110), the tribo-fatigue entropy of volume V will be: where: The principal feature of tribo-fatigue total S TF and specific s TF entropies should be emphasized. They allow the difference between two states to be assesses not only quantitatively (as thermomechanical entropy), but also qualitatively, as the value of the limiting density of the internal energy u 0 is explicitly introduced into the calculation of the specific entropy s TF . Thus, s TF and S TF allow one to answer the question how much the current state of a solid or a system is dangerous in comparison with limiting states.
The total entropy and the rate of its change for a solid of a system with regard to (111) and (113) assume the form: S " and: Based on formulas (106) (117) In this respect, bearing in mind the limiting states of a solid or a system, models (115)-(117) permit one to answer the question whether the current state is the point of a qualitative jump in the system, i.e., whether the current state is close to the limiting one. A similar (dialectical as a matter of fact) qualitative transition differs, for example, from the bifurcation point having the uncertainty in a further development of events and the possibility to predict the system behavior after a transition on the basis of the analysis of s TF and S TF . Particular limiting states (limit of strength, mechanical or contact fatigue, etc.) enable one to predict the situation after transition the given point: principal changes in the system properties and behavior or the formation of a new system based on the previous one.
As an example, there can be non-linear deformation or generation of microcracks in the solid (or the system) that cause the changes in its strength and fatigue properties, and, hence, to its response to loads. In turn, formed macrocracks lead to local continuum violation-formation of new free surfaces (possibly, of new solids-destruction products), i.e., a new system.
It should be noted that models (115)-(117) were built using a traditional concept of entropy additivity (10) although with regard to substational improvements. These models also contain reversible processes described by the entropy components s l not yielding primary damages and, hence, the limiting states-the points of a qualitative change of the system.
It is therefore more advisable for a qualitative and quantitative analysis of evolution of systems (whose states are traditionally defined as bifurcation branches) that the entropy state should be determined using the mechanothermodynamical model of the solid using only tribo-fatigue entropy. In this case, formulas (111)-(113) for entropy and its accumulation will be of the form: and: To identify the points of qualitative change in the limiting states of solids (systems), the indices of relative integral entropy and its accumulation built on the basis of the concept of integral condition of damageability of a solid (99) can be used: The indices S TF , S TF ptq, ω S , ω S ptq can grow infinitely, allowing not only the limiting states of type (93), but also different transmitting states to be described; in essence, they "provide" a quantitative description of the law of increase of entropy. Now based on formulas (24), (115), (117) and (119) let us construct generalized expresssions for entropy, a rate of its change as well as its accumulation in the mechanothermodynamical system consisting of a liquid (gas) medium of volume V and a solid of volume V ψ : Entropy state indices can be built similarly for a system composed of many media. It should be noted that the interaction (contact) of two media in formulas (122)-(125), which can be complex in nature, is taken into account only implicitly in terms of medium state parameters (stresses, strains, temperature). It is obvious that this is only the first step to a comprehensive (generalized) solution of the problem stated.
The simplified writing of expression (123) for the entropy increment of the mechanothermodynamical system composed of finite volumes dV and dV ψ given in [68] can be re-written in the following form: Expression (125) can also be represented in terms of specific quantities as: or on the basis of (123) as: In formulas (111)-(113) for calculation of the tribo-fatigue entropy S TF and its accumulation S TF ptq, the specific entropy s TF is assumed to be integrated in terms only of the damageable part of the solid-the dangerous volume. However the influence of undamagable regions can also be allowed for by integrating S TF over the entire volume: where: From expression (131), it is seen that ψ e f f u < 0 is observed outside the dangerous volume (at u e f f Σ ă u 0 ). This means that the specific tribo-fatigue entropy s TF also appears to be negative (or less than unity for its alternative definition) outside the dangerous volume where the limiting state is not reached. Negative values of ψ e f f u and s TF can then be interpreted as the absence of damageability or in other words as the retention of structure and/or properties of the solid.
As follows from the above-stated, the assumption on the entropy additivity is wrong in the general case for a system composed of both a solid and a liquid (gas) where chemical reactions can occur. By analogy with Λ-functions, interaction functions of different-nature energy (179) it is necessary to introduce interaction functions of different-nature entropy by adding them to expression (125) in effort to determine total effective entropy: or : where the subscripts Q and Ch denote the thermodynamic and chemical entropy components.
Formulas (132)-(133) are supplemented by the generalized interaction functions Λ pSq TzTF , Λ pSq QzCh , Λ pSq TzTFzCh in mechanothermodynamical systems. This means that the hypothesis on the additivity of thermodynamic and tribo-fatigue entropy is not accepted. The appropriate Λ-functions of interactions should then be specified and introduced into Equations (132)-(133).

Laws of Mechanothermodynamics
From the above-stated three main conclusions follow:

1.
Damages are the fundamental physical property (and the functional duty) of any system and all of its elements.

2.
Damage of each object (any existing one) inevitably grows up to its breakdown-decomposition (disintegration) into a set of particles of arbitrarily small size, i.e., it is the unidirectional process of time.

3.
Not only the unity and struggle of opposites but also the directivity of various and complex physical processes of hardening-softening (depending on the level of loads and time) are typical of the evolution of the system by damage. It means that the Λ-function of damage interactions (of all kinds) can take three classes of values: (a) Λ < 1 when the hardening process is dominant; (b) Λ > 1 when the softening process is dominant; (c) Λ = 1 when a stable hardening to softening process ratio is found.
Thus, the First Law of Mechanothermodynamics (see Figure 2) states that for evolution of any system the unidirectional process of its damage and disintegration, finally, into an infinitely large number of small components (fragments, atoms, elementary particles, etc.) is inevitable. In fact, it is equivalent to the recognition of the thesis on the unending of evolution if it is taken into account that disintegration products of any system become a construction material for new systems. Thus, the evolution hysteresis is formed. Generalizing, it may be said that our Universe is indestructible as it evolves by damage. This corresponds to the philosophical concept that matter and motion are eternal, and damage is the fundamental property (and a duty) of all systems, including living and intelligent ones [1,67,96,97].
The Second Law of Mechanothermodynamics states Figure 3): Λ-functions of interactions must take three classes of values (Λ ¡ 1) to describe not only the unity and struggle, but also the directivity of physical hardening-softening processes in the system, i.e., the system evolution by damage [1,67].
It can be seen that the state of a system can be equivalently described in terms of either energy or entropy. The main drawback of such descriptions is the known unreality of energy and, hence, of entropy: physical energy carriers are not detected and, apparently, do not exist. As Feynman [98] said figuratively, they cannot be touched. Damages are quite a different matter: they are physically real, can be touched, actually define any of the conceivable states of material bodies and systems. Kinetic process of their accumulation, as well as the time stream is inevitable and unidirectional. If Mechanothermodynamics considers the damage of a system as its fundamental physical property (and duty), one can hope that based on it, the consistent general theory of evolution of any systems, including living and intelligent, can be created. For instance, the idea of life as of a special method of damage accumulation (biological, mechanical, intellectual, etc.) is being developed in [68,97,99]. The Second Law of Mechanothermodynamics states (see Figure 3): Λ-functions of interactions must take three classes of values (Λ 1) to describe not only the unity and struggle, but also the directivity of physical hardening-softening processes in the system, i.e., the system evolution by damage [1,67].
It can be seen that the state of a system can be equivalently described in terms of either energy or entropy. The main drawback of such descriptions is the known unreality of energy and, hence, of entropy: physical energy carriers are not detected and, apparently, do not exist. As Feynman [98] said figuratively, they cannot be touched. Damages are quite a different matter: they are physically real, can be touched, actually define any of the conceivable states of material bodies and systems. Kinetic process of their accumulation, as well as the time stream is inevitable and unidirectional. If Mechanothermodynamics considers the damage of a system as its fundamental physical property (and duty), one can hope that based on it, the consistent general theory of evolution of any systems, including living and intelligent, can be created. For instance, the idea of life as of a special method of damage accumulation (biological, mechanical, intellectual, etc.) is being developed in [68,97,99].
Thus, the attempt was made above to formulate the basic provisions of a new (or, better to say, integrated) physical discipline-Mechanothermodynamics with the use of the energy principles. This discipline combines two branches of Physics in effort not to argue or not to compete with each other, but to take a fresh look at the MTD system evolution (Figure 4).  The second principle states the driving force and the main source of the emergence and development of processes of internal damageability of any system -these are dialectic Λ-interactions of effective energy components in the absorbing medium. Λ-interaction functions must take three classes of values 1) in order to reflect not only (Λ unity and struggle direction physical hardening-, but also of softening processes in the system. Since hardening is always finite, whereas the softening intensity can be infinitely high, the interaction of such processes inevitably leads a system to or the critical damageability-limited state

Conclusion
Effective energy absorbed in the system when it reaches a critical (limiting) state is identically equal to released (and dissipated) energy during its degradation up to decomposition Effective energy conservation law , Thus, the attempt was made above to formulate the basic provisions of a new (or, better to say, integrated) physical discipline-Mechanothermodynamics with the use of the energy principles. This discipline combines two branches of Physics in effort not to argue or not to compete with each other, but to take a fresh look at the MTD system evolution (Figure 4).       Thus, Tribo-Fatigue has become a bridge for transition from Mechanics and Thermodynamics to Mechanothermodynamics. The fact that the both ways are leading to one purpose and, finally, yield the same (unified) result, means that the above-mentioned two methodologies of analysis are valid, correct and do not contradict each other. From Figure 1 it is clear that the MTD system should be followed by an object that is somewhat more complex (however simpler) than the real system, for example, the MTD system with some "elements of intelligence". The first works in this area of research are already available [100][101][102].

Analysis of Experimental Data
Experimental verification of generalized criterion (72) of the limiting state of a MTD system is extremely difficult because of the absence of relevant experimental data. Their acquisition is though very relevant but at the same time is very difficult and expensive. Therefore, in this paper, the analysis of the particular solution of criterion (77) in form of (78) is given.
From (78) it is easy to obtain a number of formulas important for application. So, the conditions of purely thermal (or thermodynamic) damage (when σ = 0 and τ w = 0) or purely mechanical damage [when T ř Ñ 0)] will be as follows: Λ nzτ´a n σ 2`a τ τ 2 w¯" u 0 .
For isothermal mechanical fatigue (when τ w = 0) we have: and for isothermal frictional fatigue (when σ = 0) we obtain: The general analysis of the above-described partial criteria allows three main conclusions to be reached:

1.
The growth of loading parameters (σ, τ w , T ř , D) results in the corresponding acceleration of reaching the limiting state (u 0 ).

2.
The limiting state of the system can also be reached by increasing only one (any) of the loading parameters (when maintaining the same values of other parameters). 3.
If Λ > 1, the damageability of the system accordingly enhances (i.e., the processes of its softening are dominant), and if Λ < 1 it slows down (i.e., the processes of its hardening appear to be preferable) in comparison with the damage due to the joint action of loading parameters alone (with no regard to the dialectic interaction of irreversible damages).
The last conclusion is also the result of a fundamentally new approach to the construction of the criterion of the limiting state of MTD systems [102]. According to this approach, not the mutual influence of the factors, but the interaction (Λ ¡ 1) of phenomena determines the damage processes in the MTD system [1,[67][68][69]102]. In this regard, the paper synthesized the results of more than 600 diverse experimental data. This permitted the generalized MTD function of critical by damage states to be revealed.
Refer to one of the special cases of criterion (78)-isothermal mechanical fatigue. From (136) it follows that: According to (138), the dependence of limiting stresses on the parameter of thermomechanic resistance C T in the double logarithmic coordinates is to be a straight line with the angular coefficient (1/2). The general regularity is as follows: the higher the value of the parameter C T , the greater is the quantity σ´1 T . Figure 6 shows a satisfactory evidence of this dependence for numerous different-grade steels tested for fatigue in different conditions [93,95,[103][104][105]. It is seen that the C T value varied by more than two orders, i.e., by a factor of 100 or more, and the values of the endurance limit σ´1 T -by more than two orders, i.e., by a factor of 10 or more, thus the testing temperature varied in the range from the helium temperature to 0.8 T s (T S is the temperature of melting). As shown in Figure 6 is the quantity Figure 6 shows a satisfactory evidence of this dependence for numerous different-grade steels tested for fatigue in different conditions [93,95,[103][104][105]. It is seen that the CT value varied by more than two orders, i.e. by a factor of 100 or more, and the values of the endurance limit T 1 − σ -by more than two orders, i.e., by a factor of 10 or more, thus the testing temperature varied in the range from the helium temperature to 0.8 Ts ( S T is the temperature of melting). As shown in Figure 6, Equation (138) adequately describes the results of 136 experiments. In Figure 7, the similar analysis is made on the basis of the test results for tension at different temperatures (σиT is the stress limit). In this case, it is taken that σ-1 = σиT in Equation (138). It is obvious: the correlation coefficient is very high-not less than r = 0.722 (very occasionally), but in the most cases it exceeds r = 0.9; the analysis includes more than 300 test results. Works [103,104] and may others contain other examples of successful experimental approbation of criterion (138). This allows us to hope that even more general criteria (for example, Equations (77) and (78)) will appear to be practically acceptable. In our opinion, further studies should confirm our hope.
Note that the attempts to construct the explicit temperature dependence of limiting stresses in uniform, semilogarithmic and logarithmic coordinates for different materials and different testing conditions are quite ineffective (Figure 8).  In Figure 7, the similar analysis is made on the basis of the test results for tension at different temperatures (σ uT is the stress limit). In this case, it is taken that σ -1 = σ uT in Equation (138). It is obvious: the correlation coefficient is very high-not less than r = 0.722 (very occasionally), but in the most cases it exceeds r = 0.9; the analysis includes more than 300 test results. Works [103,104] and may others contain other examples of successful experimental approbation of criterion (138). This allows us to hope that even more general criteria (for example, Equations (77) and (78)) will appear to be practically acceptable. In our opinion, further studies should confirm our hope. limit T 1 − σ -by more than two orders, i.e., by a factor of 10 or more, thus the testing temperature varied in the range from the helium temperature to 0.8 Ts ( S T is the temperature of melting). As shown in Figure 6, Equation (138) adequately describes the results of 136 experiments. In Figure 7, the similar analysis is made on the basis of the test results for tension at different temperatures (σиT is the stress limit). In this case, it is taken that σ-1 = σиT in Equation (138). It is obvious: the correlation coefficient is very high-not less than r = 0.722 (very occasionally), but in the most cases it exceeds r = 0.9; the analysis includes more than 300 test results. Works [103,104] and may others contain other examples of successful experimental approbation of criterion (138). This allows us to hope that even more general criteria (for example, Equations (77) and (78)) will appear to be practically acceptable. In our opinion, further studies should confirm our hope.
Note that the attempts to construct the explicit temperature dependence of limiting stresses in uniform, semilogarithmic and logarithmic coordinates for different materials and different testing conditions are quite ineffective (Figure 8).  Note that the attempts to construct the explicit temperature dependence of limiting stresses in uniform, semilogarithmic and logarithmic coordinates for different materials and different testing conditions are quite ineffective (Figure 8).

Conclusions
It is shown, that the creation of Mechanothermodynamics-a generalized physical discipline is possible by constructing two bridges. The first is tribo-fatigue entropy which paved the way from Thermodynamics to Mechanics. The second is a fundamental understanding in Tribo-Fatigue of irreversible damage of everything that exists which paved the way from Mechanics to Thermodynamics (Figures 4 and 5). The first and the second principles of Mechanothermodynamics are presented and discussed.
Main provisions (I-XV) which constitute the foundation of the general theory of evolution of mechanothermodynamical systems are formulated.
The following theories and models are developed: • energy theory of limiting states, • energy theory of damage, • foundations of the theory of electrochemical damage.
Methods and procedures for calculating effective (dangerous) energy which is spent on generation, motion and accumulation of irreversible damages are developed (see formulas (79)-(83) and corresponding text).
Fundamentals of the theory of Λ-interaction between damages caused by loads of different nature (mechanical, thermodynamic, etc.; see formulas (69), (70), (136) and corresponding text) are given. This theory allows considering the influence of the spontaneous processes of hardening-softening on the state of damage of mechanothermodynamical systems.
In practical terms, a single MTD function of critical by damage (limiting) states of metal and polymer materials operating in different conditions (formula (138)) is obtained in this work. The analysis of 136 experimental results ( Figure 6) showed that this function is fundamental: it is valid for low-, average-and high-strength states of pure metals, alloys, and polymers over a wide range of temperatures of medium (from helium temperature to 0.8 TS, where TS is the material melting temperature) and mechanical stresses (up to the strength limit for single static loading) while the fatigue life was of the order of 10 6 -10 8 cycles.
Finally, it should be noted that the research in Mechanothermodynamics is just beginning. The deepening and broadening of the scope of works in this new and promising area of knowledge are expected in the near future. In the authors' opinion, the utility of works is so great that it is difficult to foresee their results.

Conclusions
It is shown, that the creation of Mechanothermodynamics-a generalized physical discipline is possible by constructing two bridges. The first is tribo-fatigue entropy which paved the way from Thermodynamics to Mechanics. The second is a fundamental understanding in Tribo-Fatigue of irreversible damage of everything that exists which paved the way from Mechanics to Thermodynamics (Figures 4 and 5). The first and the second principles of Mechanothermodynamics are presented and discussed.
Main provisions (I-XV) which constitute the foundation of the general theory of evolution of mechanothermodynamical systems are formulated.
The following theories and models are developed: Methods and procedures for calculating effective (dangerous) energy which is spent on generation, motion and accumulation of irreversible damages are developed (see formulas (79)-(83) and corresponding text).
Fundamentals of the theory of Λ-interaction between damages caused by loads of different nature (mechanical, thermodynamic, etc.; see formulas (69), (70), (136) and corresponding text) are given. This theory allows considering the influence of the spontaneous processes of hardening-softening on the state of damage of mechanothermodynamical systems.
In practical terms, a single MTD function of critical by damage (limiting) states of metal and polymer materials operating in different conditions (formula (138)) is obtained in this work. The analysis of 136 experimental results ( Figure 6) showed that this function is fundamental: it is valid for low-, average-and high-strength states of pure metals, alloys, and polymers over a wide range of temperatures of medium (from helium temperature to 0.8 T S , where T S is the material melting temperature) and mechanical stresses (up to the strength limit for single static loading) while the fatigue life was of the order of 10 6 -10 8 cycles.
Finally, it should be noted that the research in Mechanothermodynamics is just beginning. The deepening and broadening of the scope of works in this new and promising area of knowledge are expected in the near future. In the authors' opinion, the utility of works is so great that it is difficult to foresee their results.