Finite Element Analysis of Pin-on-Disc Tribology Test: Sumit Khot, Utpal Borah
Finite Element Analysis of Pin-on-Disc Tribology Test: Sumit Khot, Utpal Borah
Scientific Officer-F, Metal Forming & Tribology Programme, Materials Technology Division, IGCAR, Kalpakkam-603102, Tamilnadu,
India
Abstract: The main objective of this work was to evaluate the state of stress and strain in Pin-on-Disc (POD) tribology test setup with of
316LN austenitic stainless steel as the tribo elements under self-mated conditions using finite element method (FEM). Type 316LN
stainless steel is a major core and structural material in the Prototype Fast Breeder Reactor (PFBR. In PFBR there are many in-core and
out-of-core component involving contact pairs of 316LN stainless steel and undergoing sliding wear during operation as well as
maintenance. Estimation of wear during operation of the machine would lead to developing appropriate wear mitigation approaches.
However, measurement of in-situ wear in machine components is very difficult, if not impossible. Finite element method (FEM) based
numerical modeling of machine operation with appropriate wear models would enable estimation of wear a-priori. As accuracy of
calculated wear values strongly depends on the state of stress and strain in the components, accurate modelling of the state of stress and
strain is essential.
Keywords: FEA, Tribology, ABAQUS, Pin-On-Disc, WEAR, state of stress, state of strain, Elastoplastic Analysis.
δε. σdΩ - Ω δd. f dΩ - ni=1 δdi .pi = 0 (6) displacement vector, and η the penalty parameter vector.
Ω
Where, Ω is the Finite Deformation Analysis volume, δε is the
Contact identification is carried according to three rules. First,
Virtual rate of deformation tensor, ζ is the Cauchy stress
a position vector is identified between the current location of a
tensor, δd is the Virtual displacement vector, f is the Body
node, Φ, and a base point in the contact sphere, Φc
force vector, and p_i is an External force vector acting on
v = Φ – Φc (12)
points n. For the Continuum element (Pin on Disc), finite
If |v| > R, where R is the radius of the sphere, no contact
element formulation is given by,
occurs. The analyses were carried out using multiple values of
m the penalty parameter in order to verify the independency of
. Fe )j + m
T
j=1(δue
T n T
j=1 (δue . Me . ue )j - i=1 δdi . pi = 0 (7) results on the selected values.
Where, m is the number of elements, ue is the Element
displacement vector, Fe the Element force vector, Me is the Usually the friction between the contacting surfaces must be
Element mass matrix. taken into the account and it must be represented by Coulomb
friction model because it influences structural behavior and
Newton’s method with updating of the tangential stiffness is the problem may be quasi-static and time independent. So the
used during each iteration. The linear equation system is equivalent frictional stress is defined by,
solved using a preconditioned conjugate gradient solver and ηe=( τ2𝑖 )1/2 , i=1,2,…,
as convergence criteria, norms for residual forces and
increment displacement fields. In an elastic-plastic problem, where the frictional shear components are proportional to the
the concentrated loads and point supports are associated with coefficient of friction and the contact pressure, i.e.
the high local stresses. To solve this condition at these Ɣ
τ𝑖 = µp. 𝑒𝑖 (13)
locations, it takes much iteration which is wasted. The system Ɣ
is preconditioned using a Crout element-by-element scheme. Where, Ɣ𝑖 is the slip rate in direction I,Ɣ𝑒 is the slip velocity,
The finite deformation analysis is formulated using a polar µ is the coefficient of friction (given a constant results
decomposition approach and the Jaumann objective stress independent value), p is the contact pressure. [15]
rate. The three essential ingredients of elastic and plastic
analysis are yield criterion, a flow rule and hardening rule. The commercial finite element (FE) software ABAQUS can
The yield criterion relates the state of stress to the onset of handle several material and structural non-linearties, such as
yielding. The flow rule relates the state of stress {ζ} with the plasticity, visco-elasticity, friction, etc. The coupled-field
corresponding increment of plastic strain {dέp} when an analyses, for instance thermal–structural, can be performed as
increment of plastic flow occurs. The hardening rule describes well.
how the yield criterion modified by straining beyond initial
yield. 8. Methodology
Elastoplasticity of the deformable material is specified Methodology includes the modeling and FEA of pin on disc
according to a bilinear Von-Mises model using isotropic model by ABAQUS v.6.10, which includes the
hardening. Since the deformable material has a very low pre-processing, analysis and post-processing stages. It is
hardening tendency, a bilinear model with a constant possible to make an exact copy of the configuration of the pin
tangential modulus is a suitable fit for the material parameters. on disc test, but this will require a very long CPU-time
The plastic modulus for isotropic hardening isgiven by, (Central Processing Unit time). This is not necessary to obtain
E.E
H’ = T (8) good results. A simplified model will be made. The model that
E−E T
is made contains as less elements as possible in order to make
Where E is the elastic modulus, the simulation run very fast. Difficulties with the settings like
ET isthe tangent modulus. materials, boundary conditions and loads can be solved easily.
The evolution equation for the yield surface by isotropic The main goal is to get the simulation running. After this step,
hardening is given by, the model will be refined in order to get more accurate results.
2 1
k = k0 + 3
H’.λ + 3
H’εpe (9) A proper friction formulation will be chosen and the number
Where, of elements will be increased to be able to check the real
λ is proportionality coefficient pressure distribution between pin and disc, due to the different
p 2
materials of the pin. In the property module the different
εe = Strain in elastic-plastic material= 𝜀 𝑝 . 𝜀 𝑝 (10) material properties and body sections are assigned to specified
3
The contact problem is kind of geometrically non-linear regions of the pin and the disc. The structure with spherical
problem that arises when different structures, or different ended pin with radius of R=5mm was thus represented in
surfaces of a single structure comes into contact or slide on 3-dimensional (3-D) sphere on plane contact model. Both the
one another with friction. Contact forces either gained or lost pin and disc are made up of 316(LN) Stainless Steel. Young’s
must be determined in order to calculate the structural modulus and Poisson’s ratio of material are 197GPa and 0.29
behavior. Contact is enforced using a penalty function respectively. Density of material is 7850 kg/m3. To find
method, i.e. the potential energy of the finite element model is plasticity, using Ludwingson Equation,
added a so-called penalty term as ζ = K1. έn1 + e(K2 + n2.έ)(14)
1 1 Where,
Π = 2 U T . K.U - U T . F + 2 P T η. P (11)
K1 = 1232.8, n1=0.312, K2 = 5.1, n2 = -70.6
Where, U is the nodal displacement vector, K the stiffness Initial yield stress = 320.27MPa at έ = 0
matrix, F the nodal force vector, P the penetration
Volume 4 Issue 4, April 2015
www.ijsr.net
Paper ID: SUB153450 1478
Licensed Under Creative Commons Attribution CC BY
International Journal of Science and Research (IJSR)
ISSN (Online): 2319-7064
Index Copernicus Value (2013): 6.14 | Impact Factor (2013): 4.438
Table 1: Plastic Properties (From Ludwingson Equation)
Yield Stress Plastic Strain From the above given table, the correct element size of
320.27 MPa 0 contact between pin and disc is 100 micron and it will be used
488.94 MPa 0.05 for further analysis of pin on disc test.
601.6 MPa 0.1
682.08 MPa 0.15
746.129 MPa 0.2
9.2Elastoplastic Analysis of Pin on Disc
799.92 MPa 0.25
846.7 MPa 0.3 In Elastoplastic analysis, both elastic and plastic properties
888.46 MPa 0.35 [12] must be put in the FEA model and check the stress-strains
properties at the load of 10N. A very fine mesh is employed in
9. Results and Discussion the contact zone of 100 micron element size of pin and disc
where the deformation gradients are the highest. This fine
mesh region would have to extend the full length of the unit
In the beginning elastic FEA was carried out to evaluate the
cell to obtain an accurate solution. The normal load of 10N is
maximum contact pressure between the pin and the disc. The
applied through the axis of the rigid pin head of radius 5mm
results from elastic FEA were compared with analytically
and the different coefficients of friction values such as 0, 0.1,
calculated maximum contact pressure from Hertz contact
and 0.3 are used in contact between pin and disc.
theory and the optimum FEA mesh size was determined.
Subsequently Elastoplastic FEA of pin on disc was done with
constant load on the pin and various friction conditions. The
state of stress and state of strain results obtained are discussed
in the following. Normalized Elasticity Modulus is 1.07x1011
N/m2.Contact area between pin and disc is 7.05 x 10-5 N/m2.
And maximum contact pressure between pin and disc is 9.6 x
108 N/m2.
The mesh sizes of pin on disc configuration were analyzed Figure 4: Elastoplastic Analysis of Pin on Disc Visualization
with the FEA approach. The plastic deformations and the
influence of friction on the contact pressure distribution were 9.3 State of Stress and State of Strain in Pin on Disc
considered to be negligible in this case.
Von Mises Stress(S,Mises):-These are the equivalent
The concentrated load was applied on pin was F=10N and the tensile stresses between Pin and Disc material when both
sizes of the elements used in mesh validation were 100 μm, materials are in sliding contact and material start to yield.
120 μm, 140 μm, 200 μm. The Von Mises stress satisfies the property that two stress
states with equal distortion energy have equal Von Mises
stress.
Table 2: Comparison between Analytical and FEA Results of Logarithmic Strain(LE):- The logarithmic strain is true
different mesh sizes strain to provide the correct measure of the final strain when
deformation takes place in a series of increments, taking
into account the influence of the strain path. It is necessary
to verify that the disc material experiences the plain strain
conditions in the central region of the wear track. For this to
be the case, the logarithmic strain components LE12, LE13
and LE22 must be equal or close to zero in the central
region. The FE results from 3D non-linear contact
calculations for a pin load of 10N and using coefficient of
friction of 0, 0.1 and 0.3.
10. Conclusion
In this work the pin-on-disc (POD) tribology test setup was
modeled using elastoplastic finite element method (FEM) for
316LN austenitic stainless steel with 10N applied load and
100 rpm speed using the commercial FEM code ABAQUS
Figure 5: State of Stress version 6.10. The FEM model was analytically validated and
optimum mesh size for the contact elements were determined
by comparing the maximum contact pressure from elastic
analysis results with analytical results from Hertz contact
theory applied to ball on flat surface. Using this optimum
mesh size elastoplastic FEM analysis was carried out to
determine the state of stress and strain in the pin and the disc
for different friction conditions between the pin and the disc.
For this type of contact problem (ball on flat surface) the FEM
mesh size for accurate maximum contact pressure is strongly
dependent on the applied load and material properties. The
optimum mesh size for 316LN austenitic stainless steel with
10N applied load is 100 micron for which the error in
maximum contact pressure is 2.7%.The von-Mises stress
reduces for elastoplastic model as compared to elastic model.
With increasing friction, both equivalent plastic strain and
von-Mises stress increase. A state of plane strain condition is
established along the mean diameter of the contact ring as
depicted by negligible ε33 below the contact point.
References
[1] Yunfeng Zhang and Zhe Yuan, ―Wear Characteristics
Analysis of Pin-on-Disc pairs based on Finite Element
Model‖, Applied Mechanics and Materials, Volumes
148-149, 2012, pp 806-809.W.-K. Chen, Linear
Figure 6: State of Strain Networks and Systems (Book style). Belmont, CA:
Wadsworth, 1993, pp. 123–135.
Table 3: State of Stress and State of Strain in Pin-on-Disc [2] Stanislaw Kucharski and ZenonMroz, ―Identification of
wear process parameters in reciprocating ball-on-disc
tests‖, Tribology International, Volume 44, 2011,
pp154–164.
[3] ArkadiuszStachowiak and WieslawZwierzycki,
―Tribocorrosion modeling of stainless steel in a sliding
pair of pin-on-plate type‖, Tribology International,
Volume 44, September 2011, pp 1216-1224.
[4] Kris De Moerlooze and Farid Al-Bender, ―An
experimental study of ball-on-flat wear on a newly
developed rotational tribometer‖, Wear, Volume 271,
2011, pp 1005-1016.
[5] I.I. Argatov, ―Asymptotic modeling of reciprocating
sliding wear with application to local interwire contact‖,
Wear, Volume 271, 2011, pp1147-1155.
[6] A. Torres Perez and S. Austen, ―A model of friction for a
pin-on-disc configuration with imposed pin rotation‖,
Mechanism and Machine Theory, Volume 46, 2011,
pp1755- 1772.