In order to measure the residual stress components in an elastic thin plate, the hole-drilling strain-gage method has been used. This method enables to determine the relation between the magnitudes and directions of the principal stresses and the strain relaxation about the hole. In the existing analytical models based on stress field, the formulations associated with the hole-drilling method are based on the assumption of an infinite plate, this may cause some errors for a finite plate and it’s difficult to validate these solutions by FE methods. Furthermore, in the composite, the displacement field is continuous but the stress field is not necessarily continuous, the displacement field based method has to be used. In the present paper an analyt-ical model based on a displacement field described by a function with coefficients to determine for a finite round thin plate is presented. The coefficients used in the displacement field are independent on the three residual stress components, and they are determined by minimization of the internal strain energy during the hole-drilling process. Once the coefficients in the dis-placement field are determined, three strains measured in three radial directions are utilized to determine the three residual stress components. The proposed analytical model can be also adapted to infinite plate by assuming that the diameter of the round plate tends to infinite.
Residual stresses represent stresses that remain in a solid material when the original cause of the stresses has been removed. They are generated in the solid material after manufacturing processes such as casting, welding, machining, heat treatment, etc.
In order to measure the residual stresses in a thin plate, the hole-drilling method has been used experimentally [
In the present paper, a displacement field based analytical model for a finite round thin plate is proposed. The proposed displacement field is expressed by three residual stress components and some coefficients to determine. These coefficients are independent of the three residual stress components. The minimization of the internal strain energy enables to determine these coefficients. Once these coefficients are determined, three strains in three radial directions are used to determine the residual stress components.
The minimization of internal strain energy replaces the force equilibrium resolution and a flexible displacement field function can adapt very well to the boundary conditions. This displacement based model has been validated by FE method in ABAQUS and it is also adapted for infinite plate mathematically.
Hole-drilling method is to determine the residual stress in a thin plate by measuring the strains in three different radial directions occurred during the hole-drilling presented in
due to two facts: the plate is thin and the dimension of the plate is much larger than the drilled hole. All the stresses in the out-plane direction are negligible.
The plate is supposed large enough, thus a round plate is taken as an example to be treated even for other different forms of the plate. Indeed, all the specimen for the strain measurement can be fixed by a large circle perimeter and the outer bound becomes a round surface.
The residual stress components in the plate are given by:
〈 σ 0 〉 = 〈 σ x 0 σ y 0 σ x y 0 〉 (1)
In order to determine the above three residual stress components, the relation between the three measured strains in three radial directions and the three residual stress components will be investigated. In the literature [
In the experimental measurement, the strains are measured in three radial directions, thus the residual stress vector is expressed in the cylindrical coordinates as follows by using coordinates transformation:
{ σ 0 } P = { σ r 0 σ θ 0 σ r θ 0 } = { σ c + σ n σ c − σ n σ τ } (2)
where σ c = σ x 0 + σ y 0 2 , σ n = σ Δ cos 2 θ + σ x y 0 sin 2 θ , σ τ = − σ Δ sin 2 θ + σ x y 0 cos 2 θ , σ Δ = σ x 0 − σ y 0 2 .
By the analysis and observation of the results obtained by using FE method in ABAQUS®, the displacement field is proposed by a polynomial of ξ − ( 2 i − 1 ) with coefficients to determine Ai, Bi and Ci and with three residual stress components, as well as the angle θ .
{ u r = ∑ i = 0 N C C i a E σ c ξ − j + ∑ i = 0 N A A i a E σ n ξ − j u θ = ∑ i = 0 N B B i a E σ τ ξ − j with ξ = r a , j = 2 i − 1 (3)
where a is the drilled hole’s radius, E is the Yong’s Module; NA, NB and NC define the degrees of the polynomial in the displacement; Ci, Ai and Bi are coefficients to determine and they are supposed to be independent of the residual stress components, Additionally σ c and σ n in ur are independent one another as ur holds for any choice of the residual stresses.
The three cylindrical components of the strain field { ε } P are given by the follows:
ε r = ∂ u r ∂ r ; ε θ = 1 r ∂ u θ ∂ θ + u r r ; γ r θ = 1 r ∂ u r ∂ θ + ∂ u θ ∂ r − u θ r (4a)
By using Equations (3) and Equation (4a), the strain in the radial direction ε r can be expressed as follows:
ε r = − ∑ i = 0 N C j C i σ c E ξ − 2 i − ∑ i = 0 N A j A i σ n E ξ − 2 i , with j = 2 i − 1 (4b)
Once the coefficients used in the displacements field are determined, Equation (4b) can be used to determine the three residual stress components by using three measured strain in three radial directions.
The hole-drilling process is supposed to be a plane stress problem, thus the final stress state after drilling can be calculated by the following:
{ σ 1 } P = { σ 0 } P + [ C ] { ε } P (5)
with [ C ] = E 1 − ν 2 [ 1 ν 0 ν 1 0 0 0 ( 1 − ν ) / 2 ]
By using Equation (4) and Equation (5), the components of the final stress can be expressed as follows:
{ σ r 1 = σ c + 1 1 − ν 2 ∑ i = 0 N C ( − j + ν ) C i ξ − 2 i σ c + σ n + 1 1 − ν 2 ( ∑ i = 0 N A ( − j + ν ) A i ξ − 2 i − 2 ν ∑ i = 0 N B B i ξ − 2 i ) σ n σ θ 1 = σ c + 1 1 − ν 2 ∑ i = 0 N C ( − ν j + 1 ) C i ξ − 2 i σ c − σ n + 1 1 − ν 2 ( ∑ i = 0 N A ( − ν j + 1 ) A i ξ − 2 i − 2 ∑ i = 0 N B B i ξ − 2 i ) σ n σ r θ 1 = σ τ + 1 ( 1 + ν ) ∑ i = 0 N A A i ξ − 2 i σ τ − 1 2 ( 1 + ν ) ∑ i = 0 N B ( 1 + j ) B i ξ − 2 i σ τ
with j = 2 i − 1 (6)
In the thin elastic plate, two boundary conditions are used:
1) For the assumption that the outer bound is very large with respect to the drilled hole, at the outer bound of the circle plate ( ξ = λ ), the displacements u r and u θ are null. Furthermore, σ c and σ σ are independent, with Equation (3) the following relation is held:
∑ i = 0 N C C i λ − 2 i + 1 = 0 ; ∑ i = 0 N A A i λ − 2 i + 1 = 0 ; ∑ i = 0 N B B i λ − 2 i + 1 = 0 (7)
For a finite plate, λ is given by the plate’s dimension. For infinite plate, λ → ∞ and Equation (7) gives:
C 0 → 0 ; A 0 → 0 ; B 0 → 0 (8)
It is worth mentioning that the values the above three coefficients tend to zero does not mean that they can be directed eliminated from the displacement field function, because their derivatives about other coefficients are not zero which are used in the internal strain energy.
2) In the inner free boundary at the hole, the radial stress and the shearing stress in the plane are null.
σ r 1 = σ r θ 1 = 0 ( ξ = r / a = 1 ) (9)
which gives the following with the independence between σ c and σ σ :
{ ∑ i = 0 N C ( − j + ν ) C i = − ( 1 − ν 2 ) ∑ i = 0 N A ( − j + ν ) A i − 2 ν ∑ i = 0 N B B i = − ( 1 − ν 2 ) ∑ i = 0 N A A i − 1 2 ∑ i = 0 N B ( 1 + j ) B i = − ( 1 + ν ) (10)
These two boundary conditions provide six equations which allows the six coefficients C 0 , C 1 , A 0 , A 1 , B 0 and B 1 to be dependent with other coefficients. And by using Equation (10) their derives by C i , A i and B i can be calculated. These derives will be always kept in the calculations even for the infinite case λ → ∞ , C 0 → 0 ; A 0 → 0 ; B 0 → 0 .
1) Choice of Ci
The coefficients Ci concerns σ c that is the average residual normal stress. Due to its smooth change of the average residual normal stress, NC = 1 is proposed. The parts concerning Ci in Equation (7) and Equation (10) given in the two boundary conditions can be expressed as follows:
{ C 0 λ + C 1 λ − 1 = 0 ( 1 + ν ) C 0 + ( − 1 + ν ) C 1 = − ( 1 − ν 2 ) (11)
The above equation gives the coefficients C0 and C1. And for the infinite plate, λ → ∞ which gives the same expression as in the Kirsch’s equation [
2) Choice of Ai and Bi
N A = 2 , N B = 2 are the minimum for the calculation. That doesn’t avoid other choices of these two parameters. A i and B i are set as unknown variables ( i > 1 ). By using the parts concerning A i and B i in Equations ((7) and (10)) given by the two boundary conditions, the coefficient derives about A i and B i ,
that is, ∂ A k ∂ A i , ∂ B k ∂ A i , ∂ A k ∂ B i , ∂ B k ∂ B i ( k = 0 , 1 ) can be calculated:
{ ∑ i = 0 N A A i λ − 2 i + 1 = 0 ; ∑ i = 0 N B B i λ − 2 i + 1 = 0 ∑ i = 0 N A ( − j + ν ) A i − 2 ν ∑ i = 0 N B B i = − ( 1 − ν 2 ) ; ∑ i = 0 N A A i − 1 2 ∑ i = 0 N B ( 1 + j ) B i = − ( 1 + ν ) (12)
The coefficients Ai and Bi to determine in the displacement field can be obtained the minimizing the internal strain energy of the circle plate during the hole-drilling process. The internal strain energy for a determined circle plate can be given by follows:
ϕ int = 1 2 ∭ V ( 〈 σ 1 〉 P + 〈 σ 0 〉 P ) { ε 1 } d V = t 2 ∫ θ = 0 θ = 2 π ∫ r = a r = λ a ( 〈 σ 1 〉 P + 〈 σ 0 〉 P ) [ C ] − 1 ( { σ 1 } P − { σ 0 } P ) r d r d θ (13)
where t is the thickness of the plate and the minimization of the internal strain energy is translated by:
∂ ϕ int ∂ A i = 0; ∂ ϕ int ∂ B i = 0 ( i > 1 ) (14)
Equation (12) and Equation (14) form a linear system of the coefficients Ai and Bi. A linear resolution is used to solve the equations.
λ → ∞ in Equation (13) and Equation (14) gives the solution for the infinite plate. As there are some infinite-value terms in the linear system of equations, the limits of the expression possessing infinite-value terms have been calculated and simplified in the deductions.
The results obtained by FE simulations in ABAQUS using the brick elements are used to validate the analytical model. The investigation of the displacement field used in the analytical model and that in the FE simulations in ABAQUS is presented hereafter.
The numerical simulation in ABAQUS of the hole-drilling process has been fulfilled by two steps:
1) Step 1: initial settings
The residual stress field < σ x 0 , σ y 0 , τ x y 0 > in all the plate is predefined and will remain unchanged during the simulation; the displacements in the vertical surface of the plate’s outer bound and the drilled hole are set zero.
2) Step 2: hole drilling process settings
Disable the fixed displacements setting on the vertical surface of the drilled hole. In this phase, the plate will be deformed with the constraints of the outer bound. The equilibrium of the plate represents the end of the hole drilling process.
The above settings have shown the differences of displacement based model and stress based model. In the displacement based model, the displacements are used as unknown variables and it’s easy to be fulfilled by FE method.
Two cases of residual stress have been taken under investigation: 1) σ x 0 = 20 Mpa; 2) τ x y 0 = 20 Mpa in order to simulate the tensile and shearing effects of the residual stresses individually. All the other cases can be considered as combined cases of the tensile and shearing effects. The geometrical parameters of plate used in
Numerical model results vs. analytical model
1) Unidimensional horizontal tension test σ x 0 = 20 Mpa;
The strain values in three radial directions have been investigated in the experimental hole drilling method. As the displacement field determines the strain and stress field, only the displacement field is necessary to be validated. For the circular structure, as the θ = π / 2 is a special case of θ = 0 , only two directions are necessary to be taken for the validation as shown in
The cylindrical displacements u r , u θ for different value of Poisson’s ratio ν in θ = 0 direction have been shown in
The cylindrical displacements u r , u θ for different value of Poisson’s ratio ν in θ = π / 4 direction have been shown in
All the non-zero displacements in
2) Shearing test τ x y 0 = 20 Mpa
The displacement distributions in two directions are shown in
Based the above observation, the best place to put the strain gages are near the hole.
3) Unidimensional horizontal tension test σ x 0 = 20 Mpa with different diameter
In order to show the present analytical solution is suitable for thin plate to satisfy the plane stress hypothesis, different diameters of the plate are studied. The unidimensional tension test σ x 0 = 20 Mpa is used.
A displacement field based analytical model to determine the residual stress components in a finite elastic thin plate with hole-drilling method has been developed.
The proposed analytical model uses the finite boundary conditions that are realist and it has been validated by FE simulations. The obtained displacement fields in different directions for pure tensile and pure shearing are in very good agreement with those obtained by shell FE simulations (Abaqus). Theoretically this model keeps the capacity to calculate the residual stress components for an infinite plate.
This proposed model enables the determination of residual stress components in a thin finite elastic plate to be fulfilled with an efficient analytical solution without using FE simulations that are time-consuming.
The future work will be carried out by using incremental hole-drilling method in a composite plate including different displacement field in z direction where the stress based model cannot be used.
The authors declare no conflicts of interest regarding the publication of this paper.
Li, Y.M., Abbes, B., Abbes, F., Guo, Y.Q. and Gong, X.L. (2018) A Displacement Based Analytical Model to Determine Residual Stress Components in a Finite Elastic Thin Plate with Hole-Drilling Method. Journal of Applied Mathematics and Physics, 6, 1571-1580. https://doi.org/10.4236/jamp.2018.68133