# Estimating Young’s Modulus of Materials by a New Three-Point Bending Method

Authors: Xiaohu Zeng; Shifeng Wen; Mingxi Li; Gongnan Xie
Journal: Advances in Materials Science and Engineering
 (2014)

---

## Abstract

A new test method based on the three-point bending test is put forward to measure Young’s modulus of materials. The simplified mechanical model is established to make theoretical derivation. This method has not only the advantages of simple specimen preparation and convenient loading device, but also higher precision than the traditional three-point bending method. The method is adopted to obtain Young’s modulus of the aluminum alloy 2024. The feasibility of the method has been demonstrated by comparisons with the corresponding results obtained from the finite element method and experiment method. And the influence of contact friction on the test accuracy is analyzed.


---

## Body

## 1. Introduction

Ultrahigh temperature materials have been used in many fields to protect structural components from the intrinsic high temperatures, for example, on cooled gas turbines blades or vanes and the vehicle nose-tip [1]; therefore how to obtain material properties accurately under the ultrahigh temperature environment is the urgent problem to be solved in current research field [2–6].There are some methods to obtain the mechanical properties of ultrahigh temperature materials especially the film/substrate system (the thermal barrier coating material) at present, such as tensile test, indentation test, and tympanic membrane test. However, all these methods are limited by the temperature and sample processing. The tensile test has a certain requirements of number and size on the sample, and the requirements cannot be met in some cases. First, it invariably runs into problems such as fixture oxidation, strength reduction, and processing difficulties due to high temperature in the fixture design process. Second, the thermal barrier coating (TBC) systems cannot be processed into standard samples. Third, the method is not suitable for ultrahigh temperature environment [7, 8]. Therefore, the indentation test has been applied to some extent. The indentation test has some advantages. First, it is local and nondestructive. Second, the test data can be measured directly in the stress concentration location. Zorzi and Perottoni studied Young’s modulus by indentation test in 2013 [9]. However, this method still has some problems. For example, its precision in ultrahigh temperature environment is hard to control because of the influence of thermal drift produced in the ultrahigh temperature testing process [10]. And the mechanical properties of materials obtained under high pressure stress do not have a good agreement with the results under tensile stresses. The bulge test is designed to determine the mechanical properties of thin film/substrate materials, but there is no report on application of the tympanic membrane test in ultrahigh temperature environment.In the paper, a new test method is put forward to make up for lack of the above tests, namely, simple bending test method. The method has not only the advantages of simple specimen preparation and convenient loading device of the traditional three-point bending test, but also higher precision because the location of maximum stress and maximum deflection is not in the same place. For the traditional three-point and four-point bending test, the bending deflection and failure locations of the sample both occur in the middle position, as well as the contact position of indenter and measuring deflection data position, so the sample’s bending deformation will greatly influence the accuracy of deflection [11]. The emphasis of the paper is to clarify the basic theory of adopting the new method to determine the Young’s modulus of single phase materials. The Young’s modulus is a basic mechanical parameter to characterize elastic deformation properties of solid materials and also an important basis for selecting mechanical components in engineering design. The test principle is described in the following part. And the feasibility of the method has been demonstrated by comparisons with the corresponding results obtained from the finite element method and experiment method. It can provide a theoretical foundation for applying the method to determine the elastic modulus of ultrahigh temperature materials.


## 2. Description of the Test Principle

The simplified geometry model of new test device is shown in Figure1. Taking the point O as the origin of coordinates, the coordinate system is established.The simplified geometry model of test device to measure Young’s modulus of materials. (a) Schematic diagram of geometry; (b) schematic diagram of the cross section of the sample; (c) force analysis diagram of the cross section of the sample.
(a)
(b)
(c)The support reactions are obtained by equilibrium conditions∑
M
O
=
0 and ∑
M
A
=
0; that is,
(1)
F
O
y
=
-
F
a
l
F
A
y
=
F
(
l
+
a
)
l
,
where F
O
y and F
A
y are the support reactions of O and A, respectively, F is the applying load, l is the distance between two supports, and a is the distance between the loading point and the support A. Then the shear forces and bending moment equations of the sample can be obtained, as follows:
(2)
F
s
(
x
1
)
=
-
F
a
l
M
(
x
1
)
=
-
F
a
l
x
1
(
0
≤
x
1
≤
l
)
F
s
(
x
2
)
=
F
M
(
x
2
)
=
F
(
x
2
-
(
l
+
a
)
)
(
l
≤
x
2
≤
l
+
a
)
,
where F
s
(
x
1
) and F
s
(
x
2
) are the shear forces of O
A and A
B, respectively, and M
(
x
1
) and M
(
x
2
) are bending moments of O
A and A
B, respectively. From the theory of Mechanics of Materials [12], the equation of bending deformation is obtained:
(3)
1
ρ
(
x
)
=
M
(
x
)
E
I
,
where ρ is the curvature radius of deflection curve, E is the Young’s modulus of the sample, and I is the inertia moment of z-axis. Because the cross section height is far smaller than the span, the influence on deformation of the shear stress is very small and can be ignored. Equation (3) can be rewritten as
(4)
1
ρ
(
x
)
=
M
(
x
)
E
I
.From the theory of Higher Mathematics, the following equation can be obtained:(5)
1
ρ
=
d
2
y
/
d
x
2
[
1
+
(
d
y
/
d
x
)
2
]
3
.The deflection curve is nearly flat in the elastic deformation stage, so the value ofd
y
/
d
x is very small. The following equation is obtained:
(6)
1
ρ
≈
d
2
y
d
x
2
,
(7)
∴
d
2
y
d
x
2
=
M
E
I
.Then the following approximate differential equations of the deflection curve can be obtained from (2) and (7):
(8)
E
I
d
2
y
1
d
x
1
2
=
-
F
a
l
x
1
,
(9)
E
I
d
y
1
d
x
1
=
-
F
a
2
l
x
1
2
+
C
1
,
(10)
E
I
y
1
=
-
F
a
6
l
x
1
3
+
C
1
x
1
+
D
1
,
(11)
E
I
d
2
y
2
d
x
2
2
=
F
(
x
2
-
(
l
+
a
)
)
,
(12)
E
I
d
y
2
d
x
2
=
F
2
(
x
2
-
(
l
+
a
)
)
2
+
C
2
,
(13)
E
I
y
2
=
F
6
(
x
2
-
(
l
+
a
)
)
3
+
C
2
x
2
+
D
2
,
where y
1 and y
2 are the deflections of O
A and A
B, respectively and C
1, C
2, D
1, and D
2 are unknown coefficients. Consider
(14)
I
=
b
(
2
h
)
3
12
=
2
b
h
3
3
.The boundary conditions and continuity and smoothness conditions are described as follows:(15)
y
1
(
0
)
=
0
(
x
1
=
0
)
d
y
1
d
x
1
=
d
y
2
d
x
2
y
1
=
y
2
=
0
(
x
1
=
x
2
=
l
)
.
Substituting (15) into (9), (10), (12), and (13), the following equations are obtained:
(16)
C
1
=
F
6
l
a
D
1
=
0
(17)
C
2
=
-
F
6
(
2
l
a
+
3
a
2
)
D
2
=
F
6
(
a
3
+
2
a
l
2
+
3
l
a
2
)
(18)
E
I
y
2
=
F
6
(
x
2
-
(
l
+
a
)
)
3
-
F
6
(
2
l
a
+
3
a
2
)
x
2
+
F
6
(
a
3
+
3
l
a
2
+
2
a
l
2
)
V
V
V
V
V
V
V
V
V
V
I
I
(
l
≤
x
2
≤
l
+
a
)
(19)
∴
E
I
y
2
=
-
F
3
(
a
3
+
a
2
l
)
(
x
2
=
l
+
a
)
.As can be seen from the above equation, the bending deflection ofB is proportional to the applying load. From (19), the following equation is obtained:
(20)
E
=
-
F
3
I
y
2
(
a
3
+
a
2
l
)
.Therefore, as long as the ratio of applying load and bending deflection in the elastic deformation stage is obtained, Young’s modulus of the sample can be calculated.


## 3. The Test Device and Materials

The new test system of determining Young’s modulus proposed in the paper is shown in Figure2. The test system is refitted by INSTRON8871 electrohydraulic servo fatigue testing machine. The system is equipped with a temperature control system and a cooling system, and it can create a high temperature environment above 800°C. The displacement meter comes with INSTRON8871 electrohydraulic servo fatigue testing machine. Its accuracy is ±0.1% which can meet the requirements of the test. The fatigue loading system of ISTRON8871 fatigue machine is adopted to apply loads.Figure 2
Schematic diagram of the test system.The test device adopted in the paper is shown in Figure3. One end of the above transitional connection rod 7 is connected directly to the electrohydraulic servo fatigue testing machine. And the other end is connected with indenter 9 by bolt connection. The two nuts 8 and 15 are used to fasten indenter 9 and the lower transitional connection rod 16. The fixing device 12 connected to the lower transitional connection rod 16 by bolt connection is used to fix sample 10. And there are removable grooves 14 below the fixing device. It can change the contact position of indenter and sample to ensure the adequacy of test results.Figure 3
Schematic diagram of the test device.The material of sample is aluminum alloy 2024. Its Young’s modulus (E
0
=
71
GPa) is known by consulting document. And its composition and content are shown in Table 1 [13]. The sample is in a rectangle shape with length of 50 mm, width of 20 mm, and height of 2 mm. The material of indenter is cast nickel-base superalloy K403. The deformation of the indenter relative to the bending deflection of the sample is too small, so its deformation is not taken into consideration during the test. As the influence of contact friction between the indenter and sample on bending deflection is very small, it is also ignored.Table 1
The composition and content of aluminum alloy 2024 (wt.%).
Si
Fe
Cu
Mn
Mg
Cr
Zn
Ti
Others
Al
Single
Total
0.50
0.50
3.8~4.9
0.30~0.9
1.2~1.8
0.10
0.25
0.15
0.05
0.15
allowance


## 4. The Finite Element Model

According to the characteristics of the new test device, the finite element model can be simplified as shown in Figure4. The finite element software ABAQUS [14] is adopted to carry on the numerical simulation in the paper. In order to accurately simulate the stress and deformation gradient in the contact region, the grids of contact region are refined. As the deformation of the indenter relative to the bending deflection of the sample can be ignored, the indenter and supports are simplified as a rigid body without deformation in the finite element analysis in order to simplify the calculation. And the influence of friction is not taken into consideration in the process of numerical simulation. The research object in this paper is Young’s modulus of materials, so the ambient temperature is not taken into consideration temporarily. The boundary condition applied to the model is that the two supports are fixed, the indenter is limited to move along the Y direction only, and the initial state of the sample is in horizontality. The load is applied on the indenter.The finite element model for numerical calculations performed in the present study. (a) Model of indentation fatigue; (b) the reference point is highlighted.
(a)
(b)


## 5. Results and Discussions

The load spectrum adopted in numerical simulation is a monotonically increasing spectrum, as shown in Figure5. Figure 6 shows the relationship between the bending deflection and the applying load of the reference point of the sample. From the figure we can see that the slope of the load spectrum has no effect on the numerical simulation results; that is to say, the load’s increase speed with time has no effect on the bending deflection of the sample. And the deformation of the sample has two obvious stages. (
1
) The first stage is stage AB: the stage where the bending deflection of the reference point monotonically increases with the applying load, that is, the elastic deformation stage. (
2
) The second stage is stage BC: the stage where the bending deflection of the reference point rapidly increases with the applying load, that is, the plastic deformation stage. From the second section we know that as long as the ratio of applying load and bending deflection in the elastic deformation stage is obtained, Young’s modulus of the sample can be calculated. Therefore, the focus of this paper is the elastic deformation stage, that is, stage AB as shown in the figure. Young’s modulus of the sample in finite element simulation is calculated. Consider E
f
=
70.24
GPa. This result has a good agreement with the Young’s modulus obtained by consulting document.Figure 5
The load spectrum adopted in numerical simulation.The results of finite element simulation. (a) The relationship between the bending deflection and the applying load of the reference point. The load spectrum isF
=
t; (b) the relationship between the bending deflection and the applying load of the reference point. The load spectrum is F
=
2
t; (c) the relationship between the bending deflection and the applying load of the reference point. The load spectrum is F
=
3
t.
(a)
(b)
(c)


## 6. Experiment Results and Discussions

The load spectrum adopted in the test is a monotonically increasing spectrum(
F
=
t
) in the same as the numerical simulation. It is applied to the sample by the loading system of ISTRON8871 fatigue machine. The test data are set to record every 0.02 s. All the tests are carried out at room temperature. The above load spectrum is applied on the sample to test three times without difference. The relationship between the bending deflection and the applying load of the sample is shown in Figure 7.The relationship between the bending deflection (d) and the applying load (F) of the sample. (a) The d-F curve of sample 1; (b) the d-F curve of sample 2; (c) the d-F curve of sample 3.
(a)
(b)
(c)From Figure7, we can see that the deformation of the sample has also two obvious stages: (
1
) stage AB, the stage where the bending deflection of the reference point monotonically increases with the applying load, that is, the elastic deformation stage; (
2
) stage BC, the stage where the bending deflection of the reference point rapidly increases with the applying load, that is, the plastic deformation stage. The focus of this paper is stage AB. The test data processing is shown in Table 2. The method adopted in calculation process is linear regression method. The test results show that Young’s modulus of materials obtained by the new test method proposed in the paper corresponds to the Young’s modulus obtained by consulting document well. And the feasibility of the method has been demonstrated.Table 2
Test data processing.
Sample number
Young’s modulus
1
E
1
=
72.6 GPa
2
E
2
=
71.8 GPa
3
E
3
=
73.5 GPa
E
t
=
(
1
/
3
)
(
E
1
+
E
2
+
E
3
)
=
72.6 GPa


## 7. Error Analysis

The test error is calculated ase
=
(
E
t
-
E
0
)
/
E
0
=
2.3
%. The factors that produce the error may be that the influences of deformation of the indenter and contact friction are not taken into consideration. And the influence of contact friction is studied as follows.Figure8 shows the relationship between the bending deflection and the applying load of the reference point of the sample under different friction forces by finite element simulation. The Young’s modulus is calculated, as shown in Table 3. From the figure, we can see that the contact friction has influence on determination of Young’s modulus and the bigger contact friction coefficient is, the greater Young’s modulus will be. The bigger contact friction coefficient will cause bigger shear stress to prevent the sample continually deforming, so the bending deflection will be smaller and Young’s modulus will be greater.Table 3
The calculation results.
Contact friction coefficient
Young’s modulus
0
70.24 GPa
0.25
75.1 GPa
0.5
79.3 GPaThe relationship between the bending deflection (d) and the applying load (F) under different friction forces. (a) The contact friction coefficient is 0. (b) The contact friction coefficient is 0.25. (c) The contact friction coefficient is 0.5.
(a)
(b)
(c)


## 8. Concluding Remarks

In the paper, a new test method and test device based on a new three-point bending test are put forward to measure Young’s modulus of materials, and the main results are summarized as follows.(1)
The new method has advantages of simple principle, convenient operation, and wide field of application. The feasibility of the test method and test device has been demonstrated through finite element simulation and test on aluminum alloy 2024.
(2)
The contact friction has influence on determination of Young’s modulus of materials and the bigger contact friction coefficient is, the greater Young’s modulus will be.



---

*Source: 189423-2014-07-17.xml*