Finite-Element-BasedMultiple Normal Loading-Unloading of an Elastic-Plastic Spherical Stick Contact

e repeated normal elastic plastic contact problem of a deformable sphere against a rigid �at under full stick contact condition is investigated with a commercial �nite element soware ANS�S. Emphasis is placed on the effect of strain hardening and hardening model with themaximum interference of load ranging from elastic to fully plastic, which has not yet been reported. Different values of tangent modulus coupled with isotropic and kinematic hardening models are considered to study their in�uence on contact parameters. Up to ten normal loading-unloading cycles are applied with a maximum interference of 200 times the interference required to initiate yielding. Results for the variation of mean contact pressure, contact load, residual interference, and contact area with the increasing number of loading unloading cycles at high hardening parameter as well as for low tangent modulus with two different hardening models are presented. Results are compared with available �nite element simulations and in situ results reported in the literature. It is found that small variation of tangent modulus results in same shakedown behavior and similar interfacial parameters in repeated loading-unloading with both the hardening rules. However at high tangent modulus, the strain hardening and hardening rules have strong in�uence on contact parameters.


Introduction
Multiple repeated normal loading unloading is common in engineering applications.Several researchers have identi�ed the use of repeated normal loading unloading in various �elds of engineering applications such as contact resistance in MEMS micro switches [1,2], head-disk interaction in magnetic storage systems [3], ultrasonic interfacial stiffness measurement [4], stamping mechanism, and in rolling element bearings, gears, cams, and so forth.Earlier the problems of multiple loading-unloading were solved assuming a speci�c pressure distribution.e contact region related to the analysis of Kapoor et al. [5], Williams et al. [6] did not exceed greatly the elastic limit; hence they assumed Hertzian [7] contact.Merwin and Johnson [8], Kulkarni et al. [9], Bhargava et al. [10,11] used Hertzian or modi�ed Herzian pressure distribution even for elastic plastic contacts.In most engineering applications, the contact deformation occurred in elastic plastic as well as in plastic region.e arbitrary selection of pressure distribution where the plasticity effects are dominant eventually would provide inaccurate solution.Beyond the Hertzian assumption of nonadhesive frictionless contact within elastic limit, commercial �nite element soware is best suited to calculate accurately the contact parameters like contact load, contact area, and pressure in the elastic plastic and plastic region [12].e contact analysis of the rough surfaces involves the study of single asperity contact.Kogut and Etsion [13] �rst provided an accurate result of elastic plastic loading of a hemisphere against a rigid �at using commercial �nite element soware ANS�S.ey used �nite element analysis under frictionless condition for a wide range of material elastic properties and sphere size to present generalized empirical relations for contact area, contact load, and mean contact pressure as a function of dimensionless interference.eir investigations in elastic, elastic plastic, and plastic regions also include the effect of plastic properties of the material.Kogut and Etsion inferred that the variation of tangent modulus, plastic property of the material, up to 10 percent of modulus of elasticity yielded in maximum 20 percent difference of contact parameters.Etsion et al. [14] extended their investigation for unloading under perfect slip contact condition with tangent modulus of 2 percent of modulus of elasticity.ey used bilinear isotropic hardening.Jackson et al. [15] studied the effect of elastic properties of material on residual deformation and concentrated on the residual stresses of the sphere under perfect slip contact condition.Kadin et al. [16] offered a multiple loading-unloading model with isotropic hardening using 0.02 tangent modulus under perfect slip contact condition.ey found that the incipient interference required for the secondary plastic �ow in the peripheral zone increases with the increase in tangent modulus.Kral et al. [17] performed repeated normal indentation of a half space by a rigid sphere under perfect slip contact condition using isotropic hardening.It was shown that the effect of plastic properties like strain hardening is more pronounced rather than the elastic properties of the material on the contact parameters during loading unloading in the elastic plastic and plastic region.ey used strain-hardening exponent () up to 0.5.Sahoo et al. [18] studied the effect of strain hardening for elastic plastic loading of a deformable sphere against a rigid �at for frictionless and non-adhesive contact by varying the hardening parameter () up to 0.5.Variation of hardening parameter up to 0.5 enabled them to study the effect of tangent modulus as high as 33 percent of modulus of elasticity.ey inferred that the increase in strain hardening increases the resistance to deformation of a material and the material becomes capable of carrying higher amount of load in a smaller contact area.However Chatterjee and Sahoo [19] observed that the contact loads a�er �rst and second loading were same under perfect slip contact condition for non-hardening material as well as material with hardening parameter of 0.5.
e aforementioned reviews are related to frictionless contact condition.Friction is useful in many engineering applications such as brakes, clutches, bolts, nuts, driving wheels on automobiles, and so forth whereas unproductive friction takes place in engines, gears, cams, bearings, and seals [20].Brizmer et al. [21] found that the pressure distribution and the stress �eld near the contact surfaces at low contact loads are different for perfect slip and full stick contact condition.Chatterjee and Sahoo [22] studied the effect of elastic properties of material for an elastic plastic loading of a deformable sphere against a rigid �at under full stick contact condition.Zait et al. [23] performed the unloading of an elastic plastic spherical contact under full stick contact condition.ey found a substantial variation in load area curve during unloading under full stick contact condition compared to that under perfect slip contact condition.Zait et al. [24] expanded the previous single complete unloading under full stick contact condition to include multiple normal loading unloading.ey used bilinear kinematic hardening with tangent modulus as 2 percent of modulus of elasticity.ey concluded identical shakedown behavior for both isotropic and kinematic hardening.Chatterjee and Sahoo [25] studied the effect of strain hardening on the elastic plastic contact loading of a deformable sphere against a rigid �at under full stick contact condition with two different hardening rules.It was observed that the material with isotropic hardening has higher amount of load carrying capacity than that of kinematic hardening particularly for higher strain hardening.Johnson [26] mentioned that the surface tractions and nonlinear material properties are the factors to in�uence the reversibility of cyclic normal loading in spherical contact problems.us under a speci�c contact condition (slip or stick), the plastic properties of the material and the hardening rule have a particular importance in multiple normal loadingunloading.ough under perfect slip (frictionless) contact condition the effect of strain hardening on loading unloading of a deformable sphere against a rigid �at is available in literature, the effect of plastic properties and hardening model for multiple normal loading-unloading for a deformable spherical contact under full stick contact condition is still missing.
e present investigation therefore aims at studying the effect of strain hardening and hardening models on contact parameters during multiple normal loading-unloading of a deformable sphere against a rigid �at under full stick contact condition using commercial �nite element so�ware ANS�S.

Theoretical Background
Figure 1 shows the contact of a deformable sphere and a rigid �at.e dashed line presents the original contour of the sphere, having a radius of , and the rigid �at before the deformation.e solid line shows the loading phase with the interference (), corresponding to the contact radius (), and the contact load ().
Multiple loading unloading consists repeated normal loading and complete unloading.During loading phase, the interference () is increased up to the speci�ed maximum dimensionless interference,  * max .During unloading, the interference () is gradually reduced.At the end of the unloading, under zero contact load and contact area, the sphere has a residual interference ( res ).erefore the original un-deformed spherical geometry is not fully recovered.e original pro�le, the deformed shape a�er loading phase and the sphere pro�le at the end of unloading in single loading unloading are represented in Figure 2. Stick contact condition is applied at the interface of the deformed sphere and the rigid �at.e stick contact condition prevents the contact point of the deformed sphere with the rigid �at from the relative displacement in the radial direction but it allows the axial detachment of the sphere surface from the rigid �at during unloading.
In the present study, multiple normal loading-unloading is simulated considering both the isotropic and kinematic hardening models.Figure 3 describes the difference between isotropic hardening and kinematic hardening by describing the development of the yield surface with progressive yielding.In isotropic or work hardening, the yield surface is uniformly spread out from the center while in kinematic hardening the yield surface translates in stress space with constant size.e variation of different contact parameters is considered with the increase in loading or unloading cycles in the present study.e contact loads, contact area, mean contact pressure in repeated normal loading are normalized by their corresponding values aer �rst loading.Similarly residual interference aer repeated unloading is normalized by its value aer �rst unloading.e contact parameters and the interferences in some cases are normalized with their critical values in a manner as described in Chatterjee and Sahoo [25].e following empirical relations are used to evaluate the critical interference (  ), which initiates the yield inception at �rst loading and the corresponding critical load (  ) and critical contact area (  ) under full stick condition: Here,   is a parameter that depends on Poisson's ratio as   = 1.24 + 1.256.e parameters , , and  are the virgin yield stress, the Young modulus, and Poisson's ratio of the sphere material, respectively, and  is the radius of the sphere.e sphere size used for this analysis is  = 1 m.e material properties used here are Young's Modulus () = 70 GPa, Poisson's Ratio () = 0.3, and Yield stress () = 100 MPa.For full stick contact condition, in�nite friction condition is adopted.

Finite Element Model
A commercial �nite element package ANSYS 11.0 is used in the present study in order to get the accurate solution of a complicated subject like multiple loading-unloading under full stick contact condition.e sphere is represented by a quarter of a circle, due to its axisymmetry.A line models the rigid �at.e elements as shown in Figure 4 consists the mesh of maximum 18653 number of plane183 �nite elements.Plane 183 is six-node triangular axisymmetric element.e mesh density at the bottom of the sphere is coarsest one and is made gradually �ner towards the sphere summit.e �nest mesh density near the contact region simultaneously allows the sphere's curvature to be captured and accurately simulated during deformation with a reduction in computation time.Window 2 of Figure 4 presents the enlarged view of the �nest mesh density at sphere summit.e detail description and boundary condition of FE model can be found in [25].

Results and Discussion
Engineering stress-strain curves are used within elastic limit.e dimension of the specimen changes substantially in the region of plastic deformation.e increment of strain in conjunction with true stress can be termed as strain hardening.
Strain hardening causes an increase in strength and hardness of the metal.Strain hardening is expressed in terms of tangent modulus (  ), which is the slope of the stress-strain curve.
Below the proportional limit the tangent modulus is the same as the Young's modulus ().Above the proportional limit the tangent modulus varies with the strain.e tangent modulus is useful in describing the behavior of materials that have been stressed beyond the elastic region.In elastic perfectly plastic cases, the tangent modulus becomes zero.Very few material exhibit elastic perfectly plastic behaviors, generally all the materials follow the multilinear behavior with some tangent modulus.is multi-linear behavior can be assumed as bilinear behavior for analysis purpose in elastic-plastic cases.In this analysis a bilinear material property, as shown in Figure 5, is used for the deformable sphere.Shankar and Mayuram [27] mentioned that the tangent modulus for the most practical materials is less than 0.05, whereas Kadin et al. [16] found the tangent modulus for most practical materials below 0.02.However both the authors used tangent modulus up to 0.1 for analytical purpose.On the other hand, Ovcharenko et al. [28] used stainless steel specimen with tangent modulus of 0.26 (Figure 6(b)) in their in situ investigation).It is also available in literature that structural steel, aluminum alloys have signi�cant amount of strain hardening.tangent modulus,   = 0.025 using kinematic hardening.is simulation was performed to validate the present results with Zait et al. [24].Zait et al. furnished the results (Figure 4) with maximum dimensionless interference of 60 using kinematic hardening.ey have shown that with small tangent modulus the materials result in elastic shakedown even under the in�uence of kinematic hardening.e agreement between the present simulation and the results of Zait et al. [24] is appreciably good.Figures 6(b) and 6(c) are the details of Figure 6(a).Figure 6(b) shows the decrease of maximum contact load with the increase in loading cycles for multiple loading unloading for the same dimensionless interference of loading.Figure 6(c) presents the increase of residual interferences with the increase in unloading cycles during ten loading unloading cycles.Zait et al. [24] also observed the decrease of maximum contact load and increase of dimensionless interference during multiple loading unloading.

Interfacial Parameters with Low Tangent Modulus.
In the �rst part of present analysis, materials are chosen as elastic perfectly plastic material and the materials with tangent modulus of 0.025, 0.05, and 0.09 for studying the effect of strain hardening with both isotropic and kinematic hardening rule.Figure 7(a) presents the normalized contact load,  0 during ten loading-unloading cycles with bilinear isotropic hardening. 0 is the contact load aer the completion of �rst loading.e maximum loading was performed up to a dimensionless interference of  * max = 100.It is found that the contact load decreases with each loading-unloading cycle.Chatterjee and Sahoo [19] found identical contact load in second loading (Figures 9(a) and 9(b)) with bilinear isotropic hardening under perfect slip contact condition.So this decrease in contact load is due to the contact conditions, which was also inferred by Zait et al. [24].e decrease of contact load with increase in loading cycles depends on strain hardening.e decrease in contact load is less with increase in tangent modulus.e decrease of contact loads aer ten loading-unloading cycles are 9.5, 7.73, 7, and 5.1 percent for materials with tangent modulus of 0.0, 0.025, 0.05, and 0.09, respectively.interference of  * max = 100.e decrement trend of contact load is same for both the hardening rule.e decrease of contact loads aer ten loading-unloading cycles are 9.5, 7.55, 6.5 and 2.8 percent for materials with tangent modulus of 0.0, 0.025, 0.05 and 0.09 respectively.It is seen from the results of both the hardening rules that the decrease of contact load is less with kinematic hardening than that of isotropic hardening.Elastic perfectly plastic materials are insensitive to the hardening rule; however the contact load using isotropic and kinematic hardening with the increase in loading cycle has signi�cant variation for higher tangent modulus.
Figure 8(a) is the plot of normalized residual displacement,  * res / * max , during ten loading unloading cycles with bilinear isotropic hardening.e residual displacement or interferences are found aer complete unloading.e contact load and area becomes zero aer complete unloading.e residual interferences are normalized by the maximum dimensionless interference of loading.e maximum dimensionless loading interference is  * max = 100.It can be seen from the plot that the normalized residual interference during ten loading-unloading cycles is dependent on tangent modulus.e normalized residual displacement decreases with the increase of tangent modulus.e normalized residual displacement in case of elastic perfectly plastic material is 3.2, 6.4, and 11.52 percent higher aer ten loading-unloading when compared with the materials having tangent modulus of 0.025, 0.05, and 0.09, respectively.Kadin et al. [16] also observed that higher hardening results in lower residual strain.Figure 8(b) shows the normalized residual displacement,  * res / * max , as a function of loading cycles with bilinear kinematic hardening.It is revealed from the �gure that aer ten loading unloading, the normalized residual displacement of elastic perfectly plastic material is 3.61, 8.26, and 20.57percent higher than that with the materials having tangent modulus of 0.025, 0.05, and 0.09, respectively.Comparing the results with isotropic and kinematic hardening, it is observed that using kinematic hardening yields less residual displacement than that with isotropic hardening.is is more pronounced for the materials with high tangent modulus.Zolotarevskiy et al. [29] studied the elastic plastic spherical contact under cyclic tangential loading in presliding with 2% hardening and found the tangential displacement at the completion of the �rst cycle is less in the case of kinematic hardening compared to that in the isotropic hardening.e agreement between Zolotarevskiy et al. and present results is excellent.
Figure 9 shows the increment of residual interference of  * res with respect to  * res0 (residual interference aer �rst unloading) during ten unloading.We have studied the effect of hardening type also.It can be seen from the �gure that the increase of residual interference aer tenth unloading with respect to the residual interference aer �rst unloading ranges from 4 to 4.5 percent irrespective of hardening type and tangent modulus except the simulation with tangent modulus of 0.09 using kinematic hardening.e simulation of material with tangent modulus of 0.09 and kinematic hardening resulted in 1.8 percent increase of residual interference aer ten unloading with respect to residual interference aer �rst unloading.e study of contact load and residual interferences during ten loading unloading with tangent modulus up to 0.09 and two different hardening type clearly indicates that the interfacial parameters depend predominantly on the extent of strain hardening and hardening rule with the increase in tangent modulus.

Interfacial Parameters with High Tangent Modulus.
It is necessary to investigate the effect of strain hardening and hardening rule on interfacial parameters with higher tangent modulus to understand the response of the materials with signi�cant amount of strain hardening such as stainless steel, structural steel and aluminum alloys in repeated loading unloading.In our forth-coming simulations, the tangent modulus (  ) is varied according to a hardening parameter ().e hardening parameter is de�ned as  =   /(    ).
We have considered four different values of , covering wide range of tangent modulus to depict the effect of strain hardening in single asperity multiple loading unloading contact analysis with other material properties being constant.e values of  used in this analysis are within range 0 ≤  ≤ 0. as most of the practical materials falls in this range [30].e value of  equals to zero indicates elastic perfectly plastic material behavior, which is an idealized material behavior.e hardening parameters used for this analysis and their corresponding   values are shown in Table 1.[21].e mean contact pressures increase with the increase of tangent modulus for both the isotropic and kinematic hardening rule.For the tangent modulus ranging from 0.09 to 0.33, the increase of mean contact pressures with isotropic hardening are in the range of 46.56 to 118 percent than that elastic perfectly plastic material aer ten loading cycles whereas the increase of mean contact pressures with kinematic hardening ranges from 42 to 113 percent for the same range of tangent modulus aer ten loading cycles.is indicates that the effect of tangent modulus or strain hardening on the mean contact pressure is signi�cantly more pronounced than the hardening rule.Figure 10(b) shows the dimensionless mean contact pressure during ten loading cycles for the maximum dimensionless interference of  * max = 100 It is clear from the �gure that the mean contact pressure intensi�ed with the increase in dimensionless interference of loading.Mean contact pressure using isotropic hardening increased from 15.98 to 34.5 percent aer ten loading cycles for the four indicated tangent modulus in the �gure with the increase of dimensionless interference of loading from 50 to 100.It is also observed from the �gure that a marginal deviation of mean contact pressure is resulted with kinematic hardening.As already mentioned earlier that the elastic perfectly plastic material is insensitive to the hardening rule, the evolution of mean contact pressure depends predominantly on the strain hardening characteristics and the intensity of loading rather than the hardening rule.

Mean
Figure 10(c) represents the evolution of dimensionless mean contact pressure at the end of each seven loading cycles for the maximum dimensionless interference of  * max = 200.e evolution of mean contact pressure clearly indicates the same qualitative trend of Figures 10(a) and 10(b).It can be seen from the �gure that the evolution of mean contact pressure with isotropic hardening is higher than that with kinematic hardening.e dimensionless mean contact pressures with isotropic hardening are 6.4,8.9, and 8.8 percent higher aer �rst loading than that with kinematic hardening for the materials having tangent modulus of 0.09, 0.23, and 0.33, respectively.
Kral et al. [17] simulated repeated indentation of a half space by a rigid sphere with isotropic hardening.ey used three strain hardening exponent as 0, 0.3, and 0.5 and ere is no scope of comparison as we have simulated with a deformable sphere against a rigid �at but the �ualitative natures of present results (Figures 10(a)-10(c)) are in good agreement with the results of Kral et al. e deviation of mean contact pressure with the increase of loading cycles for varying tangent modulus in case of isotropic and kinematic hardening is shown in Figure 11.e maximum dimensionless interference of loading is  * max = 100,  0 represents the mean contact pressure aer �rst loading.It is found from the �gure that the variation of mean contact pressure aer every cyclic loading depends on the tangent modulus or the strain hardening characteristics and hardening type.It can also be seen from the �gure that the mean contact pressure decreases with the increasing loading cycle for the materials with isotropic hardening irrespective of the extent of tangent modulus though the decrease is not so prominent with kinematic hardening.Kral et al. [17] also mentioned that the maximum contact pressure is lower in second load half cycle in case of an indentation of a half space by a rigid sphere with isotropic hardening.e decrease of mean contact pressures with isotropic hardening aer second loading compared to �rst loading are �.3, �.8, 3.8, and 0.2 percent for the materials with tangent modulus of 0, 0.09, 0.23, and 0.33, respectively.e decrement of mean contact pressure with isotropic hardening is consistent with increasing loading cycles, aer tenth loading the decreases are 15.74, 12.5, 6.5, and 0.3 percent compared to the mean contact pressure aer �rst loading for the materials with tangent modulus of same order as mentioned for second loading.
It is therefore imperative that the maximum decrement rate of mean contact pressure is occurred aer second loading cycle and the decreases in subsequent cycles are insigni�cant.Kral et al. also observed the same pressure pro�le distribution aer second load half cycle.However the variations of mean contact pressure with the increasing loading cycles are not pronounced in case of kinematic hardening.

Dimensionless Load with Loading
Cycles.e variation of normalized contact load up to ten loading cycles for both isotropic and kinematic hardening rule is plotted in Figure 12(a). 0 is the contact load aer �rst loading and the maximum dimensionless loading for the simulation of repeated loading unloading is  * max = 100.e contact load decreases with the increase in loading cycle and the rate of decrease is the function of hardening rule and tangent modulus.e decrease of contact load with the increase of loading cycle is more pronounced in isotropic hardening than the kinematic hardening.e material with less tangent modulus produced more variation in contact load with the increase in each load cycle.
Figure 12(b) shows the variation of normalized contact load up to seven cyclic loading for the maximum loading of dimensionless interference  * max = 200.Figures 12(a) and 12(b) demonstrate that the variation of contact load with each loading cycle signi�cantly depends on the extent of maximum loading interference.e decrease of contact load with maximum dimensionless loading interference of 200 aer seventh cyclic loading is 2% more even aer the tenth cyclic loading with maximum dimensionless loading interference of 100 for the elastic perfectly plastic material.However the materials with high tangent modulus (  = 0.33 and kinematic hardening produced no variation of contact load with the increase of loading cycle.4.2.3.Residual Displacements. Figure 13(a) shows normalized residual interference,  * res / * max , aer each unloading during ten loading unloading cycles for maximum dimensionless interference of loading,  * max = 100 ( * max = 459 for elastic perfectly plastic material).It reveals from the �gure that the residual displacement depends predominantly on the tangent modulus and hardening rule.e residual displacement for the materials with the largest tangent modulus (  = 0.33 ) using isotropic hardening is about 61% that of the elastic perfectly plastic material aer �rst unloading.Kral et al. [17] found half residual displacement for material with strain hardening exponent 0.5 compared with the elastic perfectly plastic material aer �rst unloading under perfect slip contact condition.ey simulated with maximum dimensionless load,  * max = 300 and inferred that above  * max = 50, the contact load increases significantly with the increase of strain hardening for the same indentation depth.Chatterjee and Sahoo [25] reported the same qualitative results while studying the effect of strain hardening during the loading of a deformable sphere against a rigid �at under full stick contact condition.e residual displacement with isotropic hardening is signi�cantly higher than that of with kinematic hardening for higher tangent modulus.�er �rst unloading the residual displacements with isotropic hardening are 9, 52, and 86 percent higher than the residual displacement with kinematic hardening for the materials with tangent modulus of 0.09, 0.23, and 0.33, respectively.
Figure 13(b) shows the normalized residual displacement aer each unloading for seven loading unloading cycles.e maximum dimensionless interference of loading for this simulation is  * max = 200.e normalized residual displacements aer �rst unloading increased up to 12% for the increase of dimensionless interference of loading from 100 to 200.Hardening rule and the extent of loading has a signi�cant effect on residual displacement.
Figure 14(a) presents the increment of residual interference with the increase of unloading cycles compared with the residual interference aer �rst unloading. * res0 is the residual interference aer �rst unloading.It reveals from the �gure that the increases in residual interferences with isotropic hardening are signi�cantly more than that with kinematic hardening.e �gure also demonstrates that the increment of residual interference with the increase in unloading cycles depends on tangent modulus.Higher tangent modulus yields slightly less increment of residual interference with the increase in unloading cycles. of residual interference with increasing loading unloading cycles.is may cause early shakedown at higher maximum dimensionless interference of loading.and 116) are used in this simulation and the results are compared with the available literature.Ovcharenko et al. [28] found in their experiments an increase of around 5.3% and 8% of  0 with elastic perfectly plastic copper sphere ( = 400) using maximum dimensionless loading ( * ) of 48 and 116 at the end of tenth loading.Zait et al. [24] estimated 8-9% growth of  0 aer sixth loading with the deformable sphere ( = 1000) against a rigid �at using maximum dimensionless loading ( * ) of 40.We have used deformable sphere of E/Y ratio 700, our simulation resulted 12% and 23.5% increase of  0 for the maximum loading ( * ) of 48 and 116, respectively, aer tenth loading.e agreement between the present results and the �ndings of �ait et al. [24] is appreciably good.e variation of our results with that of Ovcharenko et al. [28] is due to the fact they have used elastic �at contrary to our assumption of rigid �at.�oreover the variation can be attributed to some extent to the discrete meshing of contact elements with �nite resolution and the magnitude of the contact element stiffness used in the �nite element soware during evaluating contact radius.It is also evident from the �gure that the growth of normalized contact area increases with the increase of maximum loading.Figure 15(b) shows the growth of normalized contact area aer each loading cycle with four different tangent modules using isotropic and kinematic hardening.It is evident from the plot that the contact area increases slightly aer each loading cycles indicating marginal effect of strain hardening in isotropic hardening.e materials with kinematic hardening exhibit very less increase of contact area during ten loading cycles with maximum dimensionless interference of loading  * = 100.Small discrepancy was observed for the results of higher tangent modulus (  = 0.23 and 0.33) in kinematic hardening, probably due to the discretization of contact elements as discussed earlier.

Conclusion
e effect of strain hardening and hardening models on the contact parameters during the multiple normal loading unloading under full stick contact were analyzed with �nite element soware ANSYS.e sphere was deformed against a rigid �at with the maximum interference up to 200 times the interference necessary for initial yield.Maximum ten repeated loading unloading cycles were performed with different plastic properties using isotropic and kinematic hardening models.e contact load decreases and the residual interference increases with the increase in loading unloading cycles.e variations of contact parameters with the increase in loading unloading cycles are predominantly due to the contact conditions.Under full stick contact condition, the effect of hardening model on the variation of contact parameters is more pronounced at large tangent modulus.us small variation of tangent modulus resulted same shakedown behavior and similar interfacial parameters in repeated loading unloading with both the hardening rule as reported in literature.However at high tangent modulus, the strain hardening and hardening rules have strong in�uences on contact parameters like mean contact pressure, contact load, residual interference and contact area.
It was found that the evolution of mean contact pressure depends predominantly on the strain hardening characteristics and the intensity of loading.e mean contact pressure increases with the increase in tangent modulus and intensity of loading.However the variation of mean contact pressure owing to the change of hardening model is not signi�cant even at large load.e mean contact pressure decreases with the increase in loading cycles in isotropic hardening.is decrease is prominent aer second loading although due to the increased interfacial conformity, the decrease is not so pronounced in subsequent cycles.Kinematic hardening produced insigni�cant variation of mean contact pressure with increasing loading cycles.e contact load decreases with the increase in loading cycles.e rate of decrease increases with tangent modulus and the intensity of maximum loading.e decrease of contact load is large enough in isotropic hardening for the materials with low tangent modulus compared to kinematic hardening for the materials with high tangent modulus.
e residual interference varies with both tangent modulus and hardening model.e residual interference decreases with the increase in tangent modulus and it is signi�cantly higher in isotropic hardening compared to that with kinematic hardening.e increase in residual interference with the increasing number of unloading cycles was severely affected by hardening model at high tangent modulus.e contact area increases slightly aer each loading cycles in isotropic hardening whereas the materials with kinematic hardening exhibit very less increase of contact area.e growth of contact area was found to be independent of strain hardening but increases with the increase in intensity of loading.e effect of hardening model on contact parameters at high tangent modulus clearly indicates different shakedown behavior for isotropic and kinematic hardening for repeated normal loading unloading under full stick contact condition.e analysis of hysteretic loops at high tangent modulus is an interesting �eld for future study.

F 1 :
A deformable sphere pressed by a rigid �at.

F 2 :
ree di�erent pro�les of the sphere.

Figure 6 (F 3 :
a) represents the load interference hysteretic loop during ten repeated loading unloading.e maximum dimensionless interference for loading is  * = 100, with (a) Isotropic and (b) kinematic hardening models for two-dimensional stress �eld.

F 4 :F 5 :
Finite element mesh of a sphere generated by ANSYS.StressStrain Stress-strain diagram for a material with bilinear properties.

Figure 7 (
b) represents the normalized contact load,   during ten loading unloading in case of bilinear kinematic hardening with maximum loading up to the dimensionless s e q u e n t l o a d i n g u n l o a d i n g c y c l e s Dimensionless normal interference ( / ) Dimensionless contact load ( / ) = 0.025 , kinematic interference ( / ) Dimensionless contact load ( / ) = 0.025 , kinematic (c) F 6: (a) e load-interference hysteretic loop during ten loading-unloading cycles following maximum loading to a dimensionless interference of  * = 100 (b) Increase of residual interferences and decrease of contact load.(c) Increase of residual interferences.

F 7 :F 8 :
Normalized contact load,  0 , as a function of loading cycles with (a) bilinear isotropic hardening, (b) bilinear kinematic hardening.Normalized residual interference,  * res / * max , as a function of loading cycles with (a) bilinear isotropic hardening, (b) bilinear kinematic hardening.

F 10 :
Dimensionless contact pressure, , as a function of loading cycles for maximum loading, (a)  * max = 50, (b)  * max = 100, (c)  * max = 200.indentation load up to 300 times the load necessary for initial yield.Kral et al. inferred that the distribution of contact pressure over the contact area increases with increasing load and strain hardening characteristics (Figures 4(b) and 8).

F 11 :
Normalized contact pressure,  0 , as a function of loading cycles for maximum loading,  * max = 100.

Figure 14 (F 12 :F 13 :
b) shows the normalized residual interference during seven loading unloading cycles for the maximum dimensionless interferences of loading  * max = 200.Figures 14(a) and 14(b) indicates that higher dimensionless interference of loading produces marginally less increment Normalized contact load,  0 , as a function of loading cycles for maximum loading, (a)  * max = 100, (b)  * max = 200.Normalized residual interference,  * res / * max , as a function of loading cycles for maximum loading, (a)  * max = 100, (b)  * max = 200.

Figure 15
(a)   presents the normalized contact area at the end of each loading. 0 is the contact area at the end of �rst loading.e present simulation consists ten loading unloading cycles.Two different maximum dimensionless loadings ( * = 48

F 14 :F 15 :
Normalized residual interference,  * res / * res0 , as a function of loading cycles for maximum loading, (a)  * max = 100, (b)  * max = 200.Normalized contact area,  0 , as a function of the number of loading cycles; (a) comparison of different studies; (b) comparison of hardening models for maximum loading,  * max = 100.

Nomenclature𝑎𝑎:
Contact area radius : Modulus of elasticity of the sphere : Yield Strength of the sphere material : Real contact area : Radius of the sphere : Contact load : Interference : Poisson's ratio of sphere : Mean contact pressure   : Tangent modulus of the sphere  * : Dimensionless contact load,   , in stick contact  * : Dimensionless contact area,   , in stick contact  * : Dimensionless interference, /  , in stick contact.Subscripts : Critical values res: Residual values following unloading max: Maximum values during loading-unloading process.Superscripts * ∶ Dimensionless.
1: Different  and   values used for the study of strain hardening effect.
Contact Pressure Distribution.Figure 10(a) is the plot of dimensionless mean contact pressure (/) as a function of ten loading cycles.e maximum dimensionless interference for loading is  * max = 0 ( * max = 196 for T