Fractal Loading Model of the Joint Interface Considering Strain Hardening of Materials

Based on fractal geometry theory, the deformation state of the four stages of the asperity elastic, first elastoplastic, second elastoplastic, and fully plastic deformation and comprehensively considering the hardness of the asperity changes with the amount of deformation in elastoplastic deformation stage due to strain hardening are considered, thereby establishing a single-loading model of the joint interface. By introducing the pushing coefficient and the asperity frequency exponent, each critical frequency exponents of asperity is obtained, and the relationship between the normal contact load and the contact area of the first elastoplastic deformation phase and the second elastoplastic deformation phase of the single asperity in the case of taking into account the change in hardness is inferred, eventually deducing the relationship between the total contact load of the joint interface and the contact area. .e analysis results show that in the elastoplastic deformation stage, when the deformation is constant, the asperity load considering the hardness change is smaller than the unconsidered load, and the difference increases with the increase of the deformation amount. .e establishment of the model provides a theoretical basis for further research on the elastoplastic contact of joint interfaces.


Introduction
e joint formed by assembly of mechanical parts is called the joint interface, which plays an important role in the transmission of motion, load, and energy in the normal operation of mechanical system.e joint interface presents a series of different curvature radius and the height asperities at the microscopic scale.
e contact between the joint interface is discontinuous and only occurs at higher asperity, lead to the real contact area accounts for only a small part of the nominal contact area, resulting in the situation of large load on a small contact area.As a result, researches on the properties of interfacial contact and stress analysis are very complex [1].erefore, to study the deformation behavior and the accurate modeling of joint interface is an important issue for in-depth understanding of the mechanism of friction, wear, lubrication, heat conduction, etc.
Statistical and fractal contact models for solving contact problems on joint interface have been used widely in this field.Statistical contact model was originally put forward by Greenwood and Williamson (GW model) and improved by many subsequent researchers [2][3][4].Zhao et al. deduced a new elastoplastic contact model of joint interface, which describing a long transition period from elastic deformation to fully plastic deformation of joint interface.It is shown that the elastoplastic contact of the asperity plays an important role in the microscopic contact behavior of the joint interface [5].Kogut and Etsion established the contact model between a single asperity and a rigid flat by means of finite element analysis and obtained the relationship between contact area and contact load of a single asperity during loading and unloading [6,7].Kadin et al. applied the conclusion of Etsion to the whole joint interface and got a statistical model of single loading and unloading of joint interface.According to his conclusion, plastic deformation and residual stress may occur in the process of loading and unloading.e actual contact area of the asperity during unloading is larger than the actual contact area of the loading process [8,9].However, the value of statistical parameters depends largely on the filter or resolution of the roughness measuring instrument, so it is not unique for a joint interface.
e fractal model was first proposed by Majumdar and Bhushan (MB model) in 1990.e model holds that the deformation of microconvex body changes from plastic deformation to elastic deformation with the increase of load, contrary to the traditional contact study [10][11][12].Many scholars put forward many kinds of fractal models based on MB fractal model and obtained more accurate contact mechanical properties of joint interface.Wang et al. modified the area distribution density function of asperity in MB model and obtained the modified model of MB elastic and plastic contact [13,14].Morag and Etsion established the elastoplastic contact fractal model of a single asperity and explained the contradiction between the deformation sequence of asperity from plastic deformation to elastic deformation in MB model and the classical Hertz contact theory [15].Tian et al. modified the model further, taking into account the change of material hardness with the change of surface depth in elastoplastic stage and established a new single-loading model of joint interface.However, the model only takes into account the transition from elastic to elastoplastic and elastoplastic to fully plastic deformation stage of the asperity.e description of the elastoplastic deformation stage is seldom involved in the model [16].Yuan et al. proposed an improved model of the fractal elastoplastic contact model of rough surface based on the MB model, so as to deduce a model of the total contact load and the total actual contact area [17].However, the model does not take into account the strain hardening phenomenon of the joint interface material, that is, the hardness of the material is no longer a constant value, but will change with the increase of the amount of deformation.Hardness is an important index to characterize the mechanical properties of materials such as elasticity, plasticity, strength, and toughness.e change of hardness value is directly related to the accuracy of calculation.According to the strain hardening criterion, the average hardness increases with the increase of deformation.e degree of plastic deformation increases, the degree of work hardening and dislocation strengthening increases, and the hardness of the material increases.Based on the above research results and fractal theory, a new hardness change function is expected to be constructed in this paper, considering that the hardness of the material changes with the deformation amount of the asperity in the elastoplastic deformation stage.In this paper, the critical conditions of elastic, elastoplastic, and plastic deformation of asperity are studied, and the four deformation ranges are revised and a fractal theoretical model describing the single loading of the joint interface is proposed.It is expected that the microscopic and macroscopic contact state of the surface of the interface can be more scientifically and reasonably described in order to provide some theoretical basis for the research of contact, friction, wear, and lubrication on the surface of mechanical parts.

Fractal Model of a Single Asperity
Majumdar et al. show that the contours of joint interface topography in practical engineering have fractal characteristics; mathematical characteristics are continuity, nondifferentiability, and self-affinity [10,11].e joint interface profile can be described by the Weierstrass-Mandelbrot (W-M) function, which is expressed as where x is the horizontal coordinate of the profile function of the joint interface, and the corresponding function value is the height of the profile; D is the fractal dimension of the surface profile (for a physically continuous surface, 1 < D < 2); G is the length scale parameter of the surface, which reflects the amplitude of Z(x), and is the measurement constant; n min is the lowest frequency exponent corresponding to the profile; and c n determines the spectrum of surface roughness, which is the frequency density control parameter, c > 1. e actual surface profile has an unstable randomness [18], and its lowest frequency is related to the sample length, which is given by c n min � 1/L.In order to satisfy the requirements for high spectral density and for phase randomization, c � 1.5.

Elastic Deformation of a Single
Asperity.On the microscopic scale, the contact between the two joint interfaces is essentially a contact between the asperity and the asperity, which can be simplified as a contact between a series of equivalent asperities on the joint interface and a rigid flat surface.Assuming that the joint interface is isotropic, there is no interaction between the asperity and the asperity during the contact process, and no large deformation will occur.e equation before deformation of the asperity with frequency exponent n is obtained as follows: Figure 1 shows an asperity in equivalent joint interface contacts with a rigid flat surface.e height of the asperity is h n , the interference of the asperity is ω n during the loading process, and the size of the substrate of the asperity is l n .According to equation (2), the curvature radius of an asperity with frequency exponent n at any point x is obtained as follows: When x � 0, the curvature radius of asperity is minimum: 2 Advances in Materials Science and Engineering e height before deformation of the asperity is In the loading process of asperity, the deformation will increase with the increase of normal contact load.Accordingly, asperity will change from elastic deformation to elastoplastic deformation and then to fully plastic deformation.
e elastic critical interference of asperity at initial yield is [16] where K � 0.454 + 041υ, υ is the Poisson ratio of the softer material, H is the hardness of the softer material, H � 2.8Y, and E ′ is the equivalent elastic modulus, ). E 1 and E 2 are, respectively, elastic modulus of two objects in contact with each other.υ 1 ′ and υ 2 ′ are, respectively, the corresponding Poisson ratios.When ω n < ω nec , a n < a nec , the asperity is in a state of elastic deformation.According to Hertz theory, the contact area of the asperity is Substituting equation (6) in equation ( 7), the critical contact area of the elastic asperity is According to the Hertz contact theory, the normal load on a single asperity is Substituting equations ( 4) and (7) in equation ( 9), we can obtain According to equations ( 6), (8), and (9), we can get the critical contact load of the elastic asperity:

Elastic-Plastic Deformation of a Single Asperity.
Literature [6] through the finite element analysis of a single asperity, it is concluded that when the asperity actual deformation is greater than the elastic critical interference (ω n > ω nec ), the yield phenomenon begins to appear, and the elastoplastic deformation of the asperity occurs.According to the results of [6], the elastoplastic deformation of asperity can be divided into two different stages according to the ratio ω n /ω nec , namely, the first elastoplastic deformation stage when 1 < ω n /ω nec ≤ 6 and the second elastoplastic deformation stage when 6 < ω n /ω nec ≤ 110.Define ω nepc � 6ω nec as the first elastoplastic critical interference, where the actual contact area is a nepc .e actual deformation of asperity ω npc � 110ω nec is defined as the second elastoplastic critical interference, and the actual contact area is a npc .e relationship between contact area-deformation and contact load-deformation in the elastoplastic deformation stage of asperity is [6] a n a nec � 0.93 From the above equations, we can get Advances in Materials Science and Engineering where f nec is contact load for ω � ω nec and f nep1 and f nep2 are contact loads in the first elastoplastic stage and the second elastoplastic stage, respectively.Both f nep1 and f nep2 obtained above are related to the hardness (H) of the material.However, according to the plastic strengthening principle, the hardness is not a constant when the material yields, but a function related to the deformation, that is, it changes with the deformation.erefore, it is not accurate to describe elastoplastic deformation by the above formula.
In order to express the characteristics of elastoplastic deformation more accurately, the concept of limit mean geometric hardness is introduced.
According to equations ( 13) and ( 14), H G (a) is fitted into the following segmented relations.
e first elastoplastic deformation stage is e second elastoplastic deformation stage is where c 1 , c 2 , c 3 , and c 4 are the coefficients to be solved.
(1) Equation ( 15) should satisfy two limiting conditions: where p ea (a) is the average contact pressure of the asperity in elastic stage, which is given by p ea (a) � (f ne /a); p epa1 (a) is the average contact pressure of the asperity in the first elastoplastic deformation stage and is given by p epa1 (a) � (f nep1 /a).Substituting equations ( 11) and ( 15) in equation ( 17), we can obtain Substituting equations ( 13) and (15) in equation ( 18), we can obtain Derived from equation ( 21), Considering the change of hardness, the normal contact load of a single asperity in the first elastoplastic stage is Substituting equations ( 15), (20), and (22) in equation ( 23), new equations are yielded: (2) Equation ( 16) should satisfy two limiting conditions: where p epa2 (a) � (f nep2 /a) is the average contact pressure of the asperity in the second elastoplastic stage.Substituting equations ( 13) and ( 16) in equation ( 25), we can obtain Substituting equations ( 14) and ( 16) in equation ( 26 (29) Substituting equation (29) in equation ( 27), we can obtain Considering the change of hardness, the normal contact load of a single asperity in the second elastoplastic stage is (31)

Full Plastic Deformation of a Single Asperity.
As the deformation continues to increase, when ω n > 110ω nec , the contact area a n > a npc and the asperity enters the stage of full plastic deformation.At this stage, the hardness of the material is no longer affected by the deformation and can be regarded as a constant.When the hardness of the material is 4 Advances in Materials Science and Engineering given, according to literature [7], the contact load and contact area of the asperity at this stage can be expressed as (32) In conclusion, with the increase of load and deformation, the contact area of the same asperity increases gradually, i.e., a nec < a nepc < a npc .With the increase of the load and contact area, the asperity underwent elastic deformation, first elastoplastic deformation, second elastoplastic deformation, and full plastic deformation successively.Under constant load and deformation, the actual contact area of the asperity is related to the radius of curvature at the vertex of the asperity.

Asperity's Frequency
Exponent n.When using W-M function to describe the surface profile of an asperity, the profile function is related to the asperity's frequency exponent.In other words, the radius of curvature at the vertex of the asperity and the height of the asperity vary with the frequency exponent when the load is constant.According to the equations ( 5)-( 7), it was found that the value of h n , R n , and ω nec correlated with the frequency exponent.When the frequency exponent is constant, the deformation of the asperity is not greater than the height of the asperity under the action of the load.In order to obtain the critical value of the frequency exponent, we take h n � ω nec , i.e., e elastic critical frequency exponent can be obtained as follows: where int{} is the integer part of the value in the parenthesis.Similarly, the first elastoplastic critical frequency exponent can be obtained: e second elastoplastic critical frequency exponent can be obtained: From the above, when the asperity frequency exponent is n min < n ≤ n ec , elastic deformation only takes place in these asperities under contact load.When n ec < n ≤ n epc , elastic deformation or the first elastoplastic deformation can take place in these asperities.When n epc < n ≤ n pc , elastic deformation, the first elastoplastic deformation, or the second elastoplastic deformation can take place in these asperities, and full plastic deformation never occur.When n pc < n ≤ n max , elastic deformation, elastoplastic deformation, or full plastic deformation can take place in these asperities.

Actual Contact Area and Normal Contact Load of Joint Interface
According to reference [10], when the asperity frequency exponent is n, the area distribution density function of the asperity on the joint interface is defined as where a nl represents the largest contact area when the asperity's frequency exponent is n.
In order to simplify equation (36), we define the area distribution function of the asperity of any frequency exponent as n n (a) � Mn(a).According to reference [17], the actual contact area of joint interface is

When the Frequency Exponent Belongs to
When the frequency exponent belongs to n min < n ≤ n ec , even if these asperities are completely deformed, only elastic deformation will occur, and a nl < a nec .In this case, the actual contact area of the joint interface is defined as A r1 : In this case, the contact load of the joint interface is as follows: Substituting equation (11) in equation (39), we can obtain

When the Frequency Exponent Belongs to
When the frequency exponent belongs to n ec < n ≤ n epc , for the case a nec < a nl ≤ a nepc , elastic deformation or the first elastoplastic deformation may take place in these asperities.At this point, the actual contact area of the joint interface consists of two parts, the elastic deformation stage and the first elastoplastic deformation stage: Advances in Materials Science and Engineering For the determined frequency exponent, the maximum actual contact area of the asperity appears at the maximum deformation amount ω n , where the maximum value of the elastic deformation phase ω n appears at ω nec , whereupon formula (42) is simplified to e contact load is given by Substituting equations ( 24) and (36) in equation ( 46), we can obtain

When the Frequency Exponent Belongs to
When the frequency exponent belongs to n epc < n ≤ n pc , for the case a nepc < a nl ≤ a npc , elastic deformation, the first elastoplastic deformation, or the second elastoplastic deformation may take place in these asperities.At this point, the actual contact area of the joint interface consists of three parts: the elastic deformation stage, the first elastoplastic deformation stage, and the second elastoplastic deformation stage: Mn(a)ada In this case, the contact load of the joint interface is as follows: When the second elastoplastic deformation occurs, the normal contact load of the joint interface is as follows: 6 Advances in Materials Science and Engineering Substituting equations ( 31) and (38) in equation ( 46), we can obtain (51)

When the Frequency Exponent Belongs to n pc < n.
When the frequency exponent belongs to n pc < n, elastic deformation, elastoplastic deformation, or full plastic deformation may take place in these asperities.e actual contact area of the joint interface can be evaluated as (52) In this case, the contact load of the joint interface is as follows: For all frequency exponents, the total actual contact area of the joint interface is e total contact load of the joint interface is e total real contact area and the total contact load in a nondimensional form can be written as follows: where A a is the nominal contact area and is given by Advances in Materials Science and Engineering

Results Analysis
In order to further analyze the above calculation results, the parameters of equivalent joint interface are taken as shown in Table 1 [16].
Figure 2 shows the relation between all critical contact areas and frequency exponents of single asperity when D � 1.5.It can be seen from the figure that as for one definite asperity, when frequency exponent n is certain, elastic critical contact area is minimum, followed by the first elastoplastic critical contact area, and the second elastoplastic critical contact area is maximum.With gradual increase of contact load, the contact area increases.e single asperity is firstly subject to elastic deformation, followed by the first elastoplastic deformation, the second elastoplastic deformation, and fully plastic deformation successively, which is consistent with typical contact mechanics theory.As for different asperities, with increase of frequency exponent, all critical contact areas decrease correspondingly, which shows that elastic critical contact area, the first elastoplastic critical contact area, and the second elastoplastic critical contact area are all related to frequency exponent n.
Figure 3 shows the relation curve between fractal dimension D and critical frequency exponent n of asperity.When fractal dimension is definite, elastic critical frequency exponent n ec , the first elastoplastic critical frequency exponent n epc , and the second elastoplastic critical frequency exponent n pc increase gradually.As shown in Figure 3, when D < 1.06, n ec , n epc , and n pc are all negative.As for asperities with minimum value of frequency exponent being 0, elastic deformation, elastoplastic deformation, and fully plastic deformation will all occur.When D � 1.13, n ec and n epc are negative and n pc is positive.At this time, as for asperities with minimum value of frequency exponent being 0, elastic deformation, elastoplastic deformation will occur, except fully plastic deformation.
For D � 1.5, G � 2.5 × 10 −9 m, H � 5.5 × 10 9 N/m 2 , we can obtain the elastic critical frequency exponent n ec � 32, the first elastoplastic critical frequency exponent n epc � 36, and the second elastoplastic critical frequency exponent n pc � 43.ese asperities whose frequency exponents range from 20 to 32 are only under elastic deformation.Elastic deformation and the first elastoplastic deformation can occur in these asperities whose frequency exponents range from 33 to 36.Elastic deformation, the first elastoplastic deformation, and the second elastoplastic deformation can occur in these asperities whose frequency exponents range from 37 to 43.When frequency exponents range from 43 to 50, all deformations types can occur in these asperities.
Figure 4 shows relation comparison diagram between contact load and contact area of single asperity with and without hardness change at the first elastoplastic stage.e comparison diagram is simulation result when n � 33.It can be seen from the figure that with gradual increase of contact area, with contact area of single asperity over 3.2 × 10 −13 m 2 , contact load of the same asperity with hardness change will be less than that without hardness change.In addition, as the amount of deformation increases, the difference between them tends to increase.Advances in Materials Science and Engineering definite, contact load of the same asperity with hardness change will be less than that without hardness change.In addition, with increase of deformation amount, the difference between them tends to increase, which is consistent with the change trend at the first elastoplastic stage.
Figure 6 shows the relation between limit mean geometric hardness and contact area (logarithm) of single asperity at the first elastoplastic deformation stage.Figure 6(a) shows relation curve that D � 1.1, 1.3, 1.5, 1.7 when n � 34. Figure 6(b) shows relation curve that n � 32, 33, 34, 35 when D � 1.5.It can be seen from Figure 6 that the limit mean geometric hardness of single asperity is related to contact area, fractal dimension, and frequency exponent in the first elastoplastic deformation stage.e limit mean geometric hardness increases with increase of contact area.When n is definite, the relation between limit mean geometric hardness and contact area of asperity is related to fractal dimension D.
e larger the D is, the more obvious the relation curve between them changes; when D is definite, the relation between limit mean geometric hardness and contact area of asperity is related to frequency exponent n. e smaller n is, the more obvious the relation curve between them changes.
Figure 7 shows the relation between limit mean geometric hardness and contact area (logarithm) of single asperity at the second elastoplastic deformation stage.During loading, the contact area increases with the increase of deformation of a single asperity.e ratio of deformation to the natural height of the asperity is defined as the pushing coefficient, namely, the pushing coefficient k � ω n /h n , 0 ≤ k ≤ 0.9.When fractal dimension is 1.5, we will research the relation between contact load and contact area of single asperity with frequency exponent n being 30, 35, and 40, respectively, during loading.
When n � 30, the asperity will only be subject to elastic deformation.During loading, even the pushing coefficient k is maximum, no plastic deformation will occur.e relation between contact area and contact load is f ∼ a 1.5 approximately, as shown in Figure 8(a).
As is shown in Figure 8(b), when n � 35, elastic deformation and the first elastoplastic deformation may take place in the asperity during loading.When the pushing coefficient k is less than 0.247, the asperity will under elastic deformation.At this time, the relation between contact area and contact load is f ∼ a 1.5 approximately; when the pushing coefficient is over 0.247, the first elastoplastic deformation occurs.At this time, the relation between contact area and contact load is f ∼ a 1.1093 approximately.As is shown in Figure 8(c), when n � 40, elastic deformation, the first elastoplastic deformation, and the second elastoplastic   Advances in Materials Science and Engineering deformation may take place in the asperity during loading.
When the pushing coefficient is greater than 0.1954, the asperity begins to enter the second elastoplastic deformation, the relation between contact area and contact load is f ∼ a 1.0977 approximately.When n � 45 and the pushing coefficient is greater than 0.472, the asperity begins to enter fully plastic deformation, the relation between contact area and contact load is f ∼ a approximately.
Figure 9 shows that when the minimum frequency exponent is 20 and the maximum value is 32, the actual contact area of the joint interface increases with the increase of the total contact load, and the relation between them is F * r ∼ A * 1.5 r approximately.During the whole deformation process, the joint interface appears to be of elastic property.
Figure 10 shows the relation between dimensionless total real contact area and the dimensionless total contact loading   Advances in Materials Science and Engineering of the joint interface when minimum value and maximum value of frequency exponent of asperities are, respectively, 20 and 36.When the pushing coefficient k is 0.1648, the asperity begins to be subject to the first elastoplastic deformation, and the relation between dimensionless total real contact area and the dimensionless total contact loading is approximately. e joint interface appears to be of elastic property.As the load continues increasing, when the pushing coefficient k is 0.5564; at this point, F * r > 6.9023 × 10 −3 , the relation between above load and area is approximately r .e joint interface presents elastoplastic properties.At this time, the first elastoplastic deformation takes place in these asperities whose frequency exponents range 33 ∼ 36 .
Figure 11 shows the relation between dimensionless total real contact area and the dimensionless total contact loading of the joint interface when minimum value and maximum value of frequency exponent of asperities are, respectively, 44 and 50.When the pushing coefficient k is 0.0621, the asperity begins to be subject to the second elastoplastic deformation, at this point, 0.0016 < F * r < 0.023, and the relation between dimensionless total real contact area and the dimensionless total contact loading is F * r ∼A * 1.0977 r approximately.When the pushing coefficient k is 0.7076, the asperity begins to be subject to the fully plastic deformation; at this point, F * r > 0.023, and the relation between dimensionless total real contact area and the dimensionless total contact loading is approximately F * r ∼ A * r .

Conclusions
(1) e contact mechanical properties of a single asperity on a joint interface are related to the frequency exponent of the asperity, while the frequency exponent of an asperity is related to the fractal dimension and the profile scale parameter.In this paper, the critical frequency exponents of each deformation stage of a single asperity are obtained, and Advances in Materials Science and Engineering the deformation characteristics of the asperity under different frequency exponents are obtained.(2) e normal contact load of a single asperity in the elastoplastic stage is related to the hardness of the material.When the material yields, the hardness H is not a constant, but a function related to the amount of deformation.In this paper, the limit mean geometric hardness is introduced to express the normal contact load of a single asperity in the first and second elastoplastic deformation stages considering the change of hardness.(3) e relationship between the contact load and contact area of a single asperity in the first and second elastoplastic deformation stages considering the change of material hardness and not considering the change of hardness is compared, respectively.When the contact area of the asperity is the same, the contact load with the change of hardness is smaller than that without the change of hardness, and the difference between them increases with the increase of deformation.(4) e limit mean geometric hardness of a single asperity is related to the contact area, fractal dimension, and frequency exponent during the elastoplastic deformation stage, and the limit mean geometric hardness increases with the increase of the contact area.When the frequency exponent is constant, the relationship between the limit mean geometric hardness and the contact area of the asperity is related to the fractal dimension.e larger the fractal dimension is, the more obvious the relationship curve between them is.When the fractal dimension is constant, the relationship between the limit mean geometric hardness of the asperity and the contact area is related to the frequency exponent.e smaller the frequency exponent is, the more obvious the change of the curve is.
(5) e relationship between dimensionless total real contact area and the dimensionless total contact loading of joint interface at each deformation stage is obtained by introducing the pushing coefficient, and the critical pushing coefficient exists.When the pushing coefficient exceeds this value, some asperities begin to deform in the next stage, and the inflection point appears on the relationship curve.(6) e current model still has some limitations on the applicable scope of materials, and further research and improvement are needed to make it more applicable in the future.Nomenclature ω n : Interference of the asperity ω nec : Elastic critical interference of the asperity ω nepc : First elastoplastic critical interference of the asperity ω npc : Second elastoplastic critical interference of the asperity a n : Contact area of the asperity a nec : Elastic critical contact area of the asperity a nepc : First elastoplastic critical contact area of the asperity a npc : Second elastoplastic critical contact area of the asperity f ne : Normal load in the elastic deformation of a single asperity f nec : Normal contact load for ω � ω nec f nep1 : Normal contact load of a single asperity in the first elastoplastic stage f nep2 : Normal contact load of a single asperity in the second elastoplastic stage f np : Normal contact load of a single asperity in the full plastic deformation stage f nep1 ′ : Normal contact load of a single asperity in the first elastoplastic stage considering the change of hardness f nep2 ′ : Normal contact load of a single asperity in the second elastoplastic stage considering the change of hardness H G1 (a n ): Limit mean geometric hardness in the first elastoplastic deformation stage H G2 (a n ): Limit mean geometric hardness in the second elastoplastic deformation stage n ec : Elastic critical frequency exponent n epc : First elastoplastic critical frequency exponent n pc : Second elastoplastic critical frequency exponent A r : Actual contact area of the joint interface A r1 : Actual contact area of the joint interface for n min < n ≤ n ec A r2 : Actual contact area of the joint interface for n ec < n ≤ n epc A r3 : Actual contact area of the joint interface for n epc < n ≤ n pc A r4 : Actual contact area of the joint interface for n pc < n F r1 : Actual contact load of the joint interface for n min < n ≤ n ec F r2 : Actual contact load of the joint interface for n ec < n ≤ n epc F r3 : Actual contact load of the joint interface for n epc < n ≤ n pc F r4 : Actual contact load of the joint interface for n pc < n p ea (a): Average contact pressure of the asperity in elastic stage p epa1 (a): Average contact pressure of the asperity in the first elastoplastic deformation stage p epa2 (a): Average contact pressure of the asperity in the second elastoplastic stage.

Figure 5
Figure 5 shows relation comparison diagram between contact load and contact area of single asperity with and without hardness change at the second elastoplastic stage.e comparison diagram is simulation result when n � 37.It can be seen from the figure that when the deformation is

Figure 2 :
Figure 2: e relationship between critical contact area and frequency exponent of a single asperity.

Figure 3 :
Figure 3: e relationship between fractal dimension D and critical frequency exponent n of a single asperity.

2 )
Contact areas of a single asperity (a/m Change of hardness is considered Change of hardness is not considered

Figure 4 : 2 )
Figure 4: e relationship between contact load and contact area of single asperity in the first elastoplastic deformation stage.

Figure 5 :
Figure 5: e relationship between contact load and contact area of single asperity in the second elastoplastic deformation stage.

2 ) 2 ) 2 )Figure 8 :
Figure 8: e relationship between contact load and contact area of a single asperity during loading.

Figure 10 : 3 nFigure 9 :
Figure 10:   e relationship between dimensionless total real contact area and the dimensionless total contact loading for 20 ≤ n ≤ 36.

Table 1 :
e parameters of equivalent joint interface.