# Modelling a Coupled Thermoelectromechanical Behaviour of Contact Elements via Fractal Surfaces

Authors: G. Mazzucco; F. Moro; M. Guarnieri
Journal: Modelling and Simulation in Engineering
 (2016)

---

## Abstract

A three-dimensional coupled thermoelectromechanical model for electrical connectors is here proposed to evaluate local stress and temperature distributions around the contact area of electric connectors under different applied loads. A micromechanical numerical model has been developed by merging together the contact theory approach, which makes use of the so-called roughness parameters obtained from experimental measurements on real contact surfaces, with the topology description of the rough surface via the theory of fractal geometry. Particularly, the variation of asperities has been evaluated via the Weierstrass-Mandelbrot function. In this way the micromechanical model allowed for an upgraded contact algorithm in terms of effective contact area and thermal and electrical contact conductivities. Such an algorithm is subsequently implemented to construct a global model for performing transient thermoelectromechanical analyses without the need of simulating roughness asperities of contact surfaces, so reducing the computational cost. A comparison between numerical and analytical results shows that the adopted procedure is suitable to simulate the transient thermoelectromechanical response of electric connectors.


---

## Body

## 1. Introduction

Several engineering applications involve connections where the electrical contact relies on a relatively weak pressure, able to ensure a partial adhesion between two elements. The real contact surfaces are not flat but include many asperities [1]; contacts occur in a number n
c of small surfaces nameda-spots, so that the effective contact area A
c is much smaller than the apparent, macroscopic, contact one A
a (Figure 1). The effective contact area can also be reduced due, for example, to natural oxidation and the presence of other superficial contaminant films, that can be present in the contact zone. A direct consequence of the fact that the electrical contact is established only through thesea-spots is that contact resistivity increases, producing the constriction resistance. The total contact resistance can be generally estimated by means of a statistical approach linking the microscopic contact physics to the meso-scale modelling [2–4]; with this approach the total constriction resistance is defined as a function of parameters such as surface roughness, material hardness, normal pressure p
c, and temperature T.Figure 1
Schematic view ofa-spot area.Numerical simulations of electrical contacts can be conducted in agreement with the contact theory [5] which, in general, is able to evaluate only the apparent contact area A
c, possibly overestimating the electrical and the thermal connection; for example, in [6, 7] the electromechanical contact has been evaluated by considering a micro-macro approach where A
c has been statistically defined.To improve the knowledge of the effective interconnection between two bodies in contact, a micromechanical numerical model is here developed to simulate the irregular contact surface. Micromechanical models require a homogenized material response which takes into account microstructural heterogeneity [8]. Analytical models of composite structures at nanoscale can be found in [9–14].More specifically, in the present work the micromechanical model is characterized by means of three-dimensional contact surfaces obtained via the fractal theory [15, 16] which makes use of roughness parameters resulting from experimental measurements on real contact surfaces to simulate the surface asperities (Figure 2).Figure 2
(a) Effective contact area; (b) typical fractal surface.
(a)
(b)Asperities heights, peaks distribution, and so forth are dependent on the contact surface preparation and the load history; in fact residual strains under elastoplastic cycles can modify the contact surfaces. The fractal theory has been followed to represent the initial asperities for a virgin material and to evaluate the effective electrical contact area; the material is modelled to behave elastoplastically and able to undergo large strains.Previous three-dimensional models have been developed [17] to evaluate the global effects in terms of stress and temperature distributions around the contact area; the analyses were supported by experimental tests to retrieve the thermoelectromechanical characteristics of interest.Here an analytical-numerical comparison is additionally carried out to validate this procedure, considering a classic case of a cylinder and a hemisphere in contact under different compression loads.Such problems may undergo instability [18–20], which would require a local investigation of the stress field in large strains [21–24]; nevertheless this topic is not addressed in the present micro-macro approach, under the assumption that the compressive load does not generate buckling.


## 2. The Analytical Model

On the macroscopic scale two contact elements pressed together apparently touch each other on areaA
a that, due to surface roughness, is much larger than area A
c of effective mechanical contact. In fact, on such a scale contact occurs between the asperities of the two surfaces in a number n
c of small surfaces per macro-unit area, nameda-spots, which both increase with the applied contact force P [25]. Current density lines are forced to pass througha-spots, causing a local increase in resistance R
c. For a single circulara-spot of radius a
c and sufficiently thick elements, R
c is analytically known. A more complex question is determining the total constriction resistance corresponding to A
a and given by manya-spots whose number and radius depend on the compressive load P, the surface roughness, and the material properties.A statistical approach that allows linking the macroscopic model needed for numerical simulation with the physics of microscopic contacts is extensively presented in [2–4] and it has been taken as reference for the present work. The statistical characterization of the surface roughness is considered as representative quantities: the mean plane between peaks and valleys of the asperities, the mean absolute asperity slope m, and the root mean square (RMS) of the surface roughness σ
r, whereas the contact between the two surfaces is expressed by the mean plane distance Y. As long as relatively weak contact forces are considered, m and σ
r can be assumed to be constant, so that only Y varies with the applied load P.When the load increases, the asperities in contact are crushed and the mean distanceY between two bodies decreases (Figure 3). We define d
m
a
x as the maximum mean distance between the two bodies (Y
=
d
m
a
x when the two bodies are in contact but the external load is equal to 0). If the external load is applied, the mean distance d is evaluated by (1)
Y
=
d
m
a
x
-
P
K
s
p
o
t
,where K
s
p
o
t is the spot stiffness in the contact surface. The spot stiffness is evaluated in [4, 17] via the relationship(2)
K
s
p
o
t
=
c
1
2
2
π
1
+
c
2
4
σ
r
m
-
c
2
,where c
1, c
2 are two material parameters.Figure 3
Connection between two contact surfaces.Starting from these conditions, the following expression is derived for thea-spot mean radius, their number, and the total contact surface(3)
a
c
=
8
π
σ
r
m
e
x
p
Y
2
2
σ
r
2
e
r
f
c
Y
2
σ
r
,
n
c
=
1
16
m
σ
r
2
A
a
e
r
f
c
Y
/
2
σ
r
e
x
p
Y
2
/
σ
r
2
,
A
c
=
n
c
π
a
c
2
=
1
2
e
r
f
c
Y
2
σ
r
A
a
.In order to relate Y to P, in agreement with most of the available literature, a plastic deformation of the asperities has been assumed, becausea-spots are very small so that weak forces can produce very high local pressures over them, highly exceeding the yield limit [3, 4]. However it must also be noted that hypothesis of elastic deformations is even reported in literature, considering that asperities may behave plastically at first but when thea-spots enlarge enough, loads are elastically supported. As far as the linear assumption is accepted, the plastic condition allows expressing the contact force as P
=
A
c
H
c, with H
c the contact micro-hardness, whereas A
a is related to the macroscopic apparent pressure p
c
=
P
/
A
a, so that (4)
p
c
=
H
c
2
e
r
f
c
d
2
σ
r
.Considering the Vickers hardness test, H
c can be found starting by a Vickers micro-hardness measure H
v
(5)
H
v
=
c
1
d
v
c
2
,where d
v is the mean indentation diagonal; that is, in a similar way(6)
H
c
=
c
1
d
c
c
2
,
d
c
=
2
π
a
c
.Considering the plurality ofa-spots and the proximity among them, each current tube can expand far from its contact point to a radius b (Figure 1). The total contact resistance R
c can be defined as the inverse of the contact conductance R
c
=
1
/
h
c that can be evaluated as (7)
h
c
=
1.25
k
c
m
σ
r
A
c
A
a
0.95in which k
c is the electrical conductivity.Moreover, as the contact surfaces generally present some level of oxidation and often also stray deposits of insulating substances, a film resistance has been added to the constriction resistance from (7). As it strongly depends on the level of cleanness of the surfaces and on the used metal, values taken from literature have been considered.


## 3. The Thermoelectromechanical Model

Under an electrical flux, bodies’ configurations are subjected to modifications caused by electrical energy dissipation. Correspondingly, a coupled thermoelectromechanical model has been accounted for.The electric field is governed by Maxwell’s equation(8)
∫
Γ
J
·
n
d
Γ
=
∫
Ω
r
c
d
Ω
,where Γ is the surface area, Ω is the volume, J is the electrical current density vector (per unit area), and r
c is the internal volumetric current source (per unit volume).J can be defined via Ohm’s law by considering the current density as a function of the electrical conductivity matrix σ
E
=
σ
E
(
T
), with T being temperature, and the gradient of the electrical potential E
E
(9)
J
=
σ
E
·
Ε
E
=
-
σ
E
·
∂
φ
∂
x
,where x is the position vector in the current configuration and φ is the electrical potential. The thermal energy generated by the electrical current can be defined by Joule’s law that evaluates the dissipated electrical energy P
E as(10)
P
E
=
E
E
·
J
.The energy released in the form of internal heat is(11)
r
=
η
P
E
,where r is the heat energy generated during dissipation and η is the energy conversion factor.Thermal variations within a body involve a new mechanical configuration. The stiffness matrix for a coupled thermoelectromechanical model results as(12)
k
m
m
k
m
T
0
k
T
m
k
T
T
k
T
φ
0
k
φ
T
k
φ
φ
,where k
m
m, k
T
T, k
φ
φ are the mechanical, thermal, and electrical parts, respectively, whereas k
m
T, k
T
m, k
φ
T, k
T
φ are the coupling terms.At present the piezoelectric behaviour is not considered and the constitutive relation can be defined as(13)
T
=
C
X
,
T
:
E
e
X
,
T
=
C
X
,
T
:
E
-
E
T
,where the stress tensor T is dependent on the logarithmic elastic strain E
e, obtained by subtracting the thermal strain tensor E
T to total logarithmic strain E, and the constitutive tensor C. Considering an elastoplastic material, the constitutive tensor C is dependent on both the position vector in the reference configuration X and temperature T.


## 4. The Contact Algorithm

In order to simulate the thermoelectric interconnection when two bodies are touching, the contact theory has been applied [5].Considering two bodies,Ω
1 and Ω
2 (Figure 4), for example, representative of connection components, two surfaces can be identified, Γ
1 (with Γ
1
∈
Ω
1, namedslave) and Γ
2 (with Γ
2
∈
Ω
2 namedmaster), where contact is possible. Theclosed contactcondition is achieved and the two bodies are in contact if the contact surface Γ
c
=
Γ
1
∩
Γ
2
≠
0.Figure 4
Contact scheme.Contact is defined when the fundamental conditions are set:(i)
Nonpenetration conditions:(14)
u
m
-
u
s
·
n
+
g
≥
0
on
Γ
c
g
=
X
-
m
-
X
s
·
n
on
Γ
c
,
whereu
i (with i equal to m, master, and s, slave) are the displacement vectors, X
i are the position vectors in the reference configuration, g is the gap function (the distance between two points in contact), and n is the normal vector (Figure 4).
(ii)
Action-reaction conditions:(15)
t
m
+
t
s
=
0
on
Γ
c
,
wheret
i are the stress vectors.
(iii)
Kuhn-Tucker conditions:(16)
g
≥
0
,
p
c
≤
0
,
p
c
·
g
=
0
,
on
Γ
c
,
wherep
c is the normal pressure: p
c
=
t
·
n.Generally, friction effects between two bodies in contact are dependent on the roughness of the contact surface; in our case, if one considers the adopted material asperity model, friction is a direct consequence of the model itself and it comes from the transversal asperity connections, so that a normal condition only has been considered (Figure 5).Figure 5
Normal contact definition in the fractal surface.The normal contact is characterized by a normal pressure, defined based on the penalty method; that is,(17)
p
c
=
k
n
·
g
n
α
,where g
n is the gap function along the normal direction and k
n is the penalty coefficient.In the developed three-dimensional numerical model here presented, characterized by tetrahedral elements, master and slave surfaces have been defined on the contact faces of the elements (Figure6). The closed contact condition has been considered in thecontact pair, defined through a slave node and a master point. On these surfaces only the gap function g is evaluated.Figure 6
Master and slave surfaces in finite elements.The master pointx
m in a contact pair is chosen as the point with the minimum distance from the slave node x
s (minimal distance rule); generally it does not coincide with a masternode but with a generic point belonging to surface Γ
2.The minimum distance between slave nodes and master points can be defined by considering that the master point is obtained by the orthogonal projection of the slave node onto the master surface (see Figure7), where the normal direction n is defined by(18)
n
=
a
1
×
a
2
a
1
×
a
2and a
i are the tangent vectors to the master surface in x
m (Figure 4) [26].Figure 7
Definition of the normal direction in a point and location of a master point.A “master element” is characterized by the elements faces, and the tangent vector at a generic point x
m can be consequently defined taking into account the shape functions derivatives in the master element. A master point can generally be represented by referring to two different reference systems: a global reference system (g.r.s.), x
m
=
x
m
(
x
1
,
x
2
,
x
3
), and a local curvilinear reference system (l.r.s.) belonging to the master surface, x
m
=
x
m
ξ
=
x
m
(
ξ
1
,
ξ
2
). The shape functions of a master element related to point x
m can be defined as(19)
x
1
m
=
∑
i
=
1
n
N
i
ξ
1
,
ξ
2
x
1
i
,
x
2
m
=
∑
i
=
1
n
N
i
ξ
1
,
ξ
2
x
2
i
,
x
3
m
=
∑
i
=
1
n
N
i
ξ
1
,
ξ
2
x
3
iin which n is the number of element nodes, N
i is the shape function for node i, and x
i
=
(
x
1
i
,
x
2
i
,
x
3
i
) are the coordinates of the master element nodes referred to as the g.r.s.


## 5. Thermoelectromechanical Conditions

The contact surfaces shown before must additionally transfer thermal and electrical potential fluxes.The heat flux per unit areaq defined between a contact pair has been assumed as(20)
q
=
h
T
m
,
T
s
,
p
c
T
m
-
T
s
,where the thermal conductivity h is dependent on the temperature T
i (with i equal to m, master, and s, slave) and the normal contact pressure p
c.The electric flux densityJ has been defined as(21)
J
=
h
c
T
a
v
g
,
p
c
φ
m
-
φ
s
=
h
c
Δ
φ with Δ
φ being the electric potential between master and slave surfaces and h
c the gap of electrical conductance or electrical conductivity in the contact zone, which is dependent on the average temperature in the contact zone T
a
v
g
=
T
m
-
T
s and on the normal pressure p
c.By assuming that the effective contact areaA
c can be represented as a circular area, the contact resistance R
c can be obtained following Yovanovich’s resistance relation [3, 4] and considering that, in general, the contact size is much larger than the mean free path of electrons; correspondingly the contact resistance can be obtained as(22)
R
c
=
ρ
m
+
ρ
s
2
D
,where ρ
i is the resistivity of the contact surface i and D is the contact diameter. Experimental tests have been conducted on samples in aluminium alloy with assumed resistivity equal to 3.4
·
10
-
8 Ωm [27]. A thermoelectrical conductance due to the interstitial gas layers [16] has been here neglected.


## 6. Roughness Tests

To evaluate the electrical resistance between two contact bodies under different compression load levels, experimental tests have been conducted with an aluminium cylinder and a hemisphere (Figure8(a)). In this simplified case the apparent contact surface A
a remains circular and can be also analytically estimated. The effective contact area, that is, the zone where the thermoelectric connections are defined, is dependent on contact surface asperities; the surface roughness in the contact area has consequently been characterized.Figure 8
Sample of electrical contact (a); roughness surface test zone (b).
(a)
(b)For each contact surface of cylinders and hemispheres, 20 roughness profiles have been evaluated, in 20 different positions (the measurements have been taken in the neighbourhood of the contact zone between the cylinder and the hemisphere).Profile lengths are equal to 4.0 mm considering a cut-off length of 0.8 mm. The sampling measures have been taken every 2μm.Typical examples of roughness profiles evaluated in different positions of the sample are shown in Figure9.Figure 9
Typical roughness profiles in two different positions.
(a)
(b)The RMS surface roughnessσ
r in the contact zone has been defined as(23)
σ
r
=
1
2
σ
r
,
1
2
+
σ
r
,
2
2
,where σ
r
,
i is the RMS value in the contact surface i. The RMS for the two contact surfaces has been obtained based on the roughness tests, reaching an average value of σ
r
=
0.105 μm. The mean absolute asperity slope m
- has been evaluated equal to 13%.


## 7. The Fractal Surface

The roughness parameters obtained via the roughness tests are dependent on the resolution of available measuring instruments because of the multiscale character of the roughness and its nonstationary features. The fractal theory can be adopted [27] to evaluate the invariant roughness parameters for all scale levels.At this purpose the numerical representations of the effective contact area have been conducted in agreement with the fractal approach. A micromechanical model has been developed to evaluate where the thermoelectrical connection occurs, taking into account the material asperity in the contact surface.So, starting from the theory of the fractal geometry, the topology description of the rough surface has been possible. The asperity heightz
(
x
,
y
) in a plane (
x
,
y
) has been calculated with the modified Weierstrass-Mandelbrot (WM) function [16] (24)
z
x
,
y
=
L
G
L
D
-
2
ln
⁡
γ
M
1
/
2
∑
m
=
1
M
∑
n
=
0
N
γ
n
D
-
3
cos
⁡
φ
m
,
n
-
cos
⁡
2
π
γ
n
x
2
+
y
2
0.5
L
cos
⁡
tan
-
1
⁡
y
x
-
π
m
M
+
φ
m
,
n
,where the mechanical characteristics are L, sample length, G, fractal roughness, D, fractal dimension; considering a two-dimensional surface, the limits of parameter D are 2 < D < 3. D represents the extent of space occupied by the rough surface; γ is the scaling parameter (equal to 1.5); M is the number of superimposed ridges used in constructing the surface profile; n is a frequency index; and φ
m
,
n is a random phase with interval [
0
–
2
π
].The only unknown variables in (24) are G and D, which can be experimentally estimated.


### 7.1. Fractal Surface Generation

To define the fractal surface two different codes have been developed: the first written in C++ language is able to define the WM coefficients after reading different 1D roughness profiles experimentally obtained, with roughness surface tests at different positions in the same area of the analysed contact surface. This code evaluates the surface conditions with different parameters such as maximum and minimum height peak present in the profiles, average height of asperities, profile slope, and asperity density distribution.These surface characteristics are necessary to obtain an equivalent fractal surface developed via a second code, written in Fortran 90, where the three-dimensional WM surface (Figure10) is defined.Figure 10
Example of a generated fractal surface (a); solid geometry (b); mesh of the numerical model (c).
(a)
(b)
(c)After adopting such a method for generating a fractal surface, a solid geometry has been then defined and afterwards a numerical model has been created (Figures10(b) and 10(c)). The irregular surface has been meshed via a triangular discretization.As shown in (24), the fractal surface can be defined provided that a “random phase” parameter, representative of an independent random variable, is known. Hence, a specific Fortran function has been called to generate randomly a phase parameter as shown in Box 1.Box 1:Random phase function in Fortran 90.
Real
∗
8 function  r
a
n
d
o
m
P
h
a
s
e
(
)
implicit none
real
∗
8      :: phi
real
∗
8, parameter :: pi = 3.141592653589793d0
callR
A
N
D
O
M
_
S
E
E
D
(
)
callR
A
N
D
O
M
_
N
U
M
B
E
R
p
h
i
randomPhase =a
b
s
2
.
d
0
∗
p
i
∗
p
h
i
return
end function randomPhase


## 7.1. Fractal Surface Generation

To define the fractal surface two different codes have been developed: the first written in C++ language is able to define the WM coefficients after reading different 1D roughness profiles experimentally obtained, with roughness surface tests at different positions in the same area of the analysed contact surface. This code evaluates the surface conditions with different parameters such as maximum and minimum height peak present in the profiles, average height of asperities, profile slope, and asperity density distribution.These surface characteristics are necessary to obtain an equivalent fractal surface developed via a second code, written in Fortran 90, where the three-dimensional WM surface (Figure10) is defined.Figure 10
Example of a generated fractal surface (a); solid geometry (b); mesh of the numerical model (c).
(a)
(b)
(c)After adopting such a method for generating a fractal surface, a solid geometry has been then defined and afterwards a numerical model has been created (Figures10(b) and 10(c)). The irregular surface has been meshed via a triangular discretization.As shown in (24), the fractal surface can be defined provided that a “random phase” parameter, representative of an independent random variable, is known. Hence, a specific Fortran function has been called to generate randomly a phase parameter as shown in Box 1.Box 1:Random phase function in Fortran 90.
Real
∗
8 function  r
a
n
d
o
m
P
h
a
s
e
(
)
implicit none
real
∗
8      :: phi
real
∗
8, parameter :: pi = 3.141592653589793d0
callR
A
N
D
O
M
_
S
E
E
D
(
)
callR
A
N
D
O
M
_
N
U
M
B
E
R
p
h
i
randomPhase =a
b
s
2
.
d
0
∗
p
i
∗
p
h
i
return
end function randomPhase


## 8. Mechanical Characteristics

An aluminium alloy has been tested. The mechanical characteristics in terms of elastic modulusE and yield stress have been obtained at room temperature via tensile tests. A thermal variation of such parameters has been accounted for in agreement with Eurocode 9 prescriptions [28]. When temperature increases, the elastic modulus and yield stress decrease (Figure 11(a)).Figure 11
Mechanical characteristics versus temperature (a); stress-strain relation at different temperatures (b).
(a)
(b)For the sake of simplicity, a bilinear temperature-dependent elastoplastic relation has been assumed for the numerical analyses (Figure11(b)).The hardness parameters have been experimentally defined by considering the Vickers measure; namely, a valueH
v
=
315 MPa at room temperature was obtained for the tested alloy.


## 9. Numerical Analyses

Two different models have been carried out in the following: a micromechanical model, to evaluate the effective contact area between two contact surfaces (and consequently upgrading the contact algorithm), and a global model where the thermoelectromechanical characteristics have been modified in agreement with the micromechanical results.


### 9.1. The Micromechanical Model

To evaluate the effective contact zone, where the real electrothermomechanical connections are established, fractal interfaces associated with the contact algorithm have been considered. In the real contact surfaces the irregular asperities generate stress concentrations that can be estimated through the fractal surfaces themselves.The WM equation used to generate the roughness surface considers a random phase parameter that allows the creation of different fractal surfaces with the same mechanical characteristics (derived by experimental tests). The contact surfaces used in the numerical models have been built by taking an average value of 10 different surfaces obtained via the WM equation to simulate the cylinder contact surface and 10 for the hemispheric one. The maximum in-plane extension of the fractal surface has been assumed 1.0 × 1.0 mm2 and the average asperity height is equal to about 0.2 μm.The independent variablesG and D have been calculated by taking into account that maximum, minimum, and average asperity heights must be the same, when comparing fractal surfaces and experimental measurements.The micro-model of the cylinder in contact with the hemisphere consists of 123610 nodes and 513833 tetrahedral elements; circular macro asperities have been created to simulate the effective contact surface, by means of surface preparation techniques.During compression of the two electrical joint surfaces, the different asperity heights come into contact at different times and the effective contact area depends on the external load as shown in Figure12.Figure 12
Effective contact surface at 79% of the applied load (a); at 89% of the applied load (b); at 100% of the applied load (c).
(a)
(b)
(c)Figures12(b) and 12(c) show that contact occurs first in the circular macro asperities and, subsequently, contact involves the internal part of the surface. As shown in Figure 13(a), the joint stiffness increases when surface asperities enter into contact but stiffness decreases at higher temperatures. The relation in terms of contact stiffness, temperature, and gap function has been used to evaluate the mechanical penalty coefficient in the global model.Figure 13
Contact stiffening at different temperatures (a); Von Mises stress during contact closure.
(a)
(b)When two asperities are touching, the transferred mechanical stress rapidly increases above the yield stress (Figure13(b)).As already stated, when compression occurs, at the contact surface the effective contact areaA
c changes with temperature, but as reported also in [27], the relation between the applied load and the effective contact area results to be linear, with decreasing slope at increasing temperatures (see Figure 14). At thermal melting point (about 550°C) the curve becomes almost horizontal.Figure 14
Normal pressure versusa-dimensional contact area.If the effective contact area is equal to zero the two surfaces are not in contact and the contact resistanceR
c is infinite. In agreement with (22), R
c decreases as A
c increases (Figure 15(a)).Figure 15
Contact resistanceR
c versusa-dimensional contact area (a); R
c versus normal pressure (b).
(a)
(b)As evidenced in Figure14, by varying normal pressure and temperature, the effective contact area A
c changes, which allows obtaining R
c law in terms of pressure and temperature variations (Figure 15(b)).Hence, the contact algorithm adopted in the macro scale model has been modified in light of the thermoelectromechanical relations coming from the micromechanical model just illustrated. In the global one, roughness contact surfaces have not been represented to reduce the computational cost, and the contact area itself gives the apparent contact areaA
a related to the pair, once the effective contact area A
c is evaluated following the relationship shown in Figure 16. Consequently, the contact resistance R
c has been defined for each contact pair.Figure 16
Relation between temperature, contact area, and normal pressure.


### 9.2. Global Model and Comparison of the Results

The global model has been carried out as illustrated in Guarnieri et al. [29], where a hemisphere and a cylinder are put in contact at different compressive load levels and different electrical potentials. These types of shapes ensure a closed solution for the analytical approach, because the resulting contact surface is always circular. The contact characteristics have been compared considering the analytical formulation and the previously described fractal characterization.The schematic view of contact between the two surfaces is shown in Figure8, where the maximum diameter for the samples d is assumed to be equal to 90 mm and the material is conceived to behave thermoelastoplastically.A current densityJ
=
1.0 A/mm2 is applied to the hemisphere after the load application, at an initial temperature of 20°C.The three-dimensional numerical model consists of 730000 tetrahedral elements (Figure17(a)) and linear shape functions. Due to symmetry, only 1/4 of the geometry has been considered.Figure 17
Three-dimensional model of the electric connector (a); sphere indentation (b).
(a)
(b)After the application of the compressive load, the sphere touches the cylinder in the contact zone and the contact pressure increases following the elastoplastic relationship. A permanent indentation occurs if the elastic field is exceeded as shown in Figure17(b) [30].The geometric configuration of the sample allowed for obtaining a circular apparent contact surface at different load levels (see Figure18), which is proved in the following to correspond to that analytically evaluated [1].Figure 18
Apparent contact surface at different load levels (at room temperatureT
= 20°C): 2% (a); 11.6% (b); 18.4% (c); 35% (d); 56% (e); 100% (f).
(a)
(b)
(c)
(d)
(e)
(f)In the contact algorithm the electrical conductanceh
c
i at the contact pair i has been obtained via the apparent contact area A
a
,
i and the effective contact area A
c
,
i, based on Yovanovich’s formulation [2, 4](25)
h
c
i
=
1.25
k
c
m
σ
A
c
,
i
A
a
,
i
0.95
,where k
c is the average electrical conductance of the two bodies in contact(26)
1
k
c
=
2
1
k
c
,
m
+
1
k
c
,
s
.
A
c
i has been calculated, via the micromechanical models, as a function of contact pressure and temperature. In the analytical and numerical models k
c has been assumed to be equal to 2.9
×
10
4 S/mm. Similarly, the thermal conductance has been evaluated and the thermal conductivity h has been assumed to be equal to 0.237 W/(K·mm). This assumption allowed evaluating the variation of the electrical and the thermal conductance at different contact positions, under different contact pressures.Analytical and numerical results have been compared by considering the contact surface at 20°C.A comparison in terms of real contact areaA
c and h
v is reported in Table 1 and in Figure 19.Table 1
Comparison between analytical and numerical results.
% load
A
a (mm2)
A
c
,
a
n
a
l
y
t
i
c
a
l (mm2)
A
c
,
f
r
a
c
t
a
l (mm2)
P (MPa)
h
c
,
a
n
a
l
y
t
i
c
a
l (S/mm2)
h
c
,
f
r
a
c
t
a
l (S/mm2)
0.00E + 00
0.00E + 00
0.00E + 00
0.00E + 00
0.00E + 00
0.00E + 00
0.00E + 00
1.00E − 02
1.26E − 01
1.02E − 02
1.27E − 03
5.06E + 01
3.98E + 03
5.48E + 02
1.00E − 01
2.89E − 01
2.63E − 02
1.24E − 02
2.20E + 02
4.44E + 03
2.17E + 03
2.50E − 01
5.15E − 01
4.97E − 02
4.02E − 02
3.09E + 02
4.69E + 03
3.83E + 03
4.00E − 01
8.49E − 01
8.14E − 02
7.65E − 02
3.00E + 02
4.69E + 03
4.39E + 03
5.50E − 01
1.32E + 00
1.24E − 01
1.17E − 01
2.64E + 02
4.69E + 03
4.39E + 03
7.50E − 01
1.99E + 00
1.81E − 01
1.84E − 01
2.40E + 02
4.69E + 03
4.50E + 03
1.00E + 00
2.90E + 00
2.64E − 01
2.74E − 01
2.20E + 02
4.69E + 03
4.60E + 03Figure 19
Apparent and real contact area at different load levels (20°C) (a); comparison between analytical and numerical electrical conductivities (b).
(a)
(b)As shown, the numerical and the analytic results forA
c are pretty similar at different load levels, which confirms that fractal surfaces are good to evaluate the real contact surface in these types of problems.A slight discrepancy is encountered in the first load steps: such difference can be explained by the fact that the analytical method is based on a statistical approach and the correspondent results are averaged, whereas the fractal approach takes into account the effective contact area found through the micromechanical model during the overall compressive load histories; that is, the analytical method is not able to locally catch the real contact area at low load values. Indeed, the deviation between the results gradually decreases when the load increases, due to the fact that the contact spots defined by the fractal model are comparable to those obtained via the statistical approach. In fact, the resulting contact areas must finally coincide, getting closer to the real one.According to the micromechanical model, the real contact areaA
c is consequently a function dependent on load and temperature, A
c
=
A
c
P
,
T. This allows evaluating the interaction between two bodies during an electrical connection, by taking into account the different thermal deformations at different positions of the contact zone during transient load scenarios.The global model is also able to evaluate the displacement variations and the stress fields (Figures20(a) and 20(b)) when an electrical current is applied at different compressive loads.Figure 20
Vertical displacements versus position at different loads levels (a); Von Mises stress (b).
(a)
(b)The contact discontinuity generates a sudden deviation in terms of electrical potential and temperature (Figure21); when the contact area increases, at higher loads, the peak of temperature and electrical potential decreases, being locally higher than the thermal and the electrical contact conductivities.Figure 21
Comparison between electrical potential versus position at different loads levels (a); comparison between temperature versus position at different loads levels (b).
(a)
(b)


## 9.1. The Micromechanical Model

To evaluate the effective contact zone, where the real electrothermomechanical connections are established, fractal interfaces associated with the contact algorithm have been considered. In the real contact surfaces the irregular asperities generate stress concentrations that can be estimated through the fractal surfaces themselves.The WM equation used to generate the roughness surface considers a random phase parameter that allows the creation of different fractal surfaces with the same mechanical characteristics (derived by experimental tests). The contact surfaces used in the numerical models have been built by taking an average value of 10 different surfaces obtained via the WM equation to simulate the cylinder contact surface and 10 for the hemispheric one. The maximum in-plane extension of the fractal surface has been assumed 1.0 × 1.0 mm2 and the average asperity height is equal to about 0.2 μm.The independent variablesG and D have been calculated by taking into account that maximum, minimum, and average asperity heights must be the same, when comparing fractal surfaces and experimental measurements.The micro-model of the cylinder in contact with the hemisphere consists of 123610 nodes and 513833 tetrahedral elements; circular macro asperities have been created to simulate the effective contact surface, by means of surface preparation techniques.During compression of the two electrical joint surfaces, the different asperity heights come into contact at different times and the effective contact area depends on the external load as shown in Figure12.Figure 12
Effective contact surface at 79% of the applied load (a); at 89% of the applied load (b); at 100% of the applied load (c).
(a)
(b)
(c)Figures12(b) and 12(c) show that contact occurs first in the circular macro asperities and, subsequently, contact involves the internal part of the surface. As shown in Figure 13(a), the joint stiffness increases when surface asperities enter into contact but stiffness decreases at higher temperatures. The relation in terms of contact stiffness, temperature, and gap function has been used to evaluate the mechanical penalty coefficient in the global model.Figure 13
Contact stiffening at different temperatures (a); Von Mises stress during contact closure.
(a)
(b)When two asperities are touching, the transferred mechanical stress rapidly increases above the yield stress (Figure13(b)).As already stated, when compression occurs, at the contact surface the effective contact areaA
c changes with temperature, but as reported also in [27], the relation between the applied load and the effective contact area results to be linear, with decreasing slope at increasing temperatures (see Figure 14). At thermal melting point (about 550°C) the curve becomes almost horizontal.Figure 14
Normal pressure versusa-dimensional contact area.If the effective contact area is equal to zero the two surfaces are not in contact and the contact resistanceR
c is infinite. In agreement with (22), R
c decreases as A
c increases (Figure 15(a)).Figure 15
Contact resistanceR
c versusa-dimensional contact area (a); R
c versus normal pressure (b).
(a)
(b)As evidenced in Figure14, by varying normal pressure and temperature, the effective contact area A
c changes, which allows obtaining R
c law in terms of pressure and temperature variations (Figure 15(b)).Hence, the contact algorithm adopted in the macro scale model has been modified in light of the thermoelectromechanical relations coming from the micromechanical model just illustrated. In the global one, roughness contact surfaces have not been represented to reduce the computational cost, and the contact area itself gives the apparent contact areaA
a related to the pair, once the effective contact area A
c is evaluated following the relationship shown in Figure 16. Consequently, the contact resistance R
c has been defined for each contact pair.Figure 16
Relation between temperature, contact area, and normal pressure.


## 9.2. Global Model and Comparison of the Results

The global model has been carried out as illustrated in Guarnieri et al. [29], where a hemisphere and a cylinder are put in contact at different compressive load levels and different electrical potentials. These types of shapes ensure a closed solution for the analytical approach, because the resulting contact surface is always circular. The contact characteristics have been compared considering the analytical formulation and the previously described fractal characterization.The schematic view of contact between the two surfaces is shown in Figure8, where the maximum diameter for the samples d is assumed to be equal to 90 mm and the material is conceived to behave thermoelastoplastically.A current densityJ
=
1.0 A/mm2 is applied to the hemisphere after the load application, at an initial temperature of 20°C.The three-dimensional numerical model consists of 730000 tetrahedral elements (Figure17(a)) and linear shape functions. Due to symmetry, only 1/4 of the geometry has been considered.Figure 17
Three-dimensional model of the electric connector (a); sphere indentation (b).
(a)
(b)After the application of the compressive load, the sphere touches the cylinder in the contact zone and the contact pressure increases following the elastoplastic relationship. A permanent indentation occurs if the elastic field is exceeded as shown in Figure17(b) [30].The geometric configuration of the sample allowed for obtaining a circular apparent contact surface at different load levels (see Figure18), which is proved in the following to correspond to that analytically evaluated [1].Figure 18
Apparent contact surface at different load levels (at room temperatureT
= 20°C): 2% (a); 11.6% (b); 18.4% (c); 35% (d); 56% (e); 100% (f).
(a)
(b)
(c)
(d)
(e)
(f)In the contact algorithm the electrical conductanceh
c
i at the contact pair i has been obtained via the apparent contact area A
a
,
i and the effective contact area A
c
,
i, based on Yovanovich’s formulation [2, 4](25)
h
c
i
=
1.25
k
c
m
σ
A
c
,
i
A
a
,
i
0.95
,where k
c is the average electrical conductance of the two bodies in contact(26)
1
k
c
=
2
1
k
c
,
m
+
1
k
c
,
s
.
A
c
i has been calculated, via the micromechanical models, as a function of contact pressure and temperature. In the analytical and numerical models k
c has been assumed to be equal to 2.9
×
10
4 S/mm. Similarly, the thermal conductance has been evaluated and the thermal conductivity h has been assumed to be equal to 0.237 W/(K·mm). This assumption allowed evaluating the variation of the electrical and the thermal conductance at different contact positions, under different contact pressures.Analytical and numerical results have been compared by considering the contact surface at 20°C.A comparison in terms of real contact areaA
c and h
v is reported in Table 1 and in Figure 19.Table 1
Comparison between analytical and numerical results.
% load
A
a (mm2)
A
c
,
a
n
a
l
y
t
i
c
a
l (mm2)
A
c
,
f
r
a
c
t
a
l (mm2)
P (MPa)
h
c
,
a
n
a
l
y
t
i
c
a
l (S/mm2)
h
c
,
f
r
a
c
t
a
l (S/mm2)
0.00E + 00
0.00E + 00
0.00E + 00
0.00E + 00
0.00E + 00
0.00E + 00
0.00E + 00
1.00E − 02
1.26E − 01
1.02E − 02
1.27E − 03
5.06E + 01
3.98E + 03
5.48E + 02
1.00E − 01
2.89E − 01
2.63E − 02
1.24E − 02
2.20E + 02
4.44E + 03
2.17E + 03
2.50E − 01
5.15E − 01
4.97E − 02
4.02E − 02
3.09E + 02
4.69E + 03
3.83E + 03
4.00E − 01
8.49E − 01
8.14E − 02
7.65E − 02
3.00E + 02
4.69E + 03
4.39E + 03
5.50E − 01
1.32E + 00
1.24E − 01
1.17E − 01
2.64E + 02
4.69E + 03
4.39E + 03
7.50E − 01
1.99E + 00
1.81E − 01
1.84E − 01
2.40E + 02
4.69E + 03
4.50E + 03
1.00E + 00
2.90E + 00
2.64E − 01
2.74E − 01
2.20E + 02
4.69E + 03
4.60E + 03Figure 19
Apparent and real contact area at different load levels (20°C) (a); comparison between analytical and numerical electrical conductivities (b).
(a)
(b)As shown, the numerical and the analytic results forA
c are pretty similar at different load levels, which confirms that fractal surfaces are good to evaluate the real contact surface in these types of problems.A slight discrepancy is encountered in the first load steps: such difference can be explained by the fact that the analytical method is based on a statistical approach and the correspondent results are averaged, whereas the fractal approach takes into account the effective contact area found through the micromechanical model during the overall compressive load histories; that is, the analytical method is not able to locally catch the real contact area at low load values. Indeed, the deviation between the results gradually decreases when the load increases, due to the fact that the contact spots defined by the fractal model are comparable to those obtained via the statistical approach. In fact, the resulting contact areas must finally coincide, getting closer to the real one.According to the micromechanical model, the real contact areaA
c is consequently a function dependent on load and temperature, A
c
=
A
c
P
,
T. This allows evaluating the interaction between two bodies during an electrical connection, by taking into account the different thermal deformations at different positions of the contact zone during transient load scenarios.The global model is also able to evaluate the displacement variations and the stress fields (Figures20(a) and 20(b)) when an electrical current is applied at different compressive loads.Figure 20
Vertical displacements versus position at different loads levels (a); Von Mises stress (b).
(a)
(b)The contact discontinuity generates a sudden deviation in terms of electrical potential and temperature (Figure21); when the contact area increases, at higher loads, the peak of temperature and electrical potential decreases, being locally higher than the thermal and the electrical contact conductivities.Figure 21
Comparison between electrical potential versus position at different loads levels (a); comparison between temperature versus position at different loads levels (b).
(a)
(b)


## 10. Conclusions

A three-dimensional coupled thermoelectromechanical model for electrical connectors has been here proposed to evaluate local stress and temperature distributions around the contact area of the connectors under different applied loads.The electrical resistance arising in the contact area varies with the effective contact areaA
c, where thermal and electrical contacts are established. It can be assumed that A
c depends on the surface roughness and the contact pressure.A micromechanical numerical model has been developed by merging together the contact theory approach [5], which makes use of roughness parameters obtained from experimental measurements on real contact surfaces, with the topology description of the rough surface via the theory of fractal geometry. Particularly, the asperities variation has been evaluated via the Weierstrass-Mandelbrot function [16].In this way the micromechanical model has allowed for obtaining an upgraded contact algorithm in terms of effective contact areaA
c, as well as thermal and electrical contact conductivities, qualifying a small region of the total contact zone.Such an algorithm has been subsequently implemented in a global model for performing transient thermoelectromechanical analyses without the need of simulating roughness asperities of contact surfaces, so reducing the computational costs.A comparison between numerical and analytical results proved that the adopted procedure is suitable to simulate the transient thermoelectromechanical response of electric connectors.



---

*Source: 5219876-2016-08-25.xml*