# A 4-Mid-Node Plane Model of Base Force Element Method on Complementary Energy Principle

Authors: Yinghua Liu; Yijiang Peng; Lijuan Zhang; Qing Guo
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

Using the base forces as fundamental variables to describe the stress state and the displacement gradients that are the conjugate variables of the base forces to describe the deformation state for the two-dimensional elasticity problems, a 4-mid-node plane model of base force element method (BFEM) based on complementary energy principle is proposed. In this paper, the complementary energy of an element of the BFEM is constructed by using the base forces. The equilibrium conditions are released by the Lagrange multiplier method, and a modified complementary energy principle described by the base forces is obtained. The formulation of the 4-mid-node plane element of the BFEM is derived by assuming that the stress is uniformly distributed on each edge of the plane elements. A procedure of the BFEM on complementary energy principle is developed using MATLAB language. The numerical results of examples show that this model of the BFEM has high precision and is free from mesh sensitivity. This model shows good performances.


---

## Body

## 1. Introduction

The finite element method based on the assumed displacement field has become a method of choice for the solution of a wide variety of problems in structural mechanics. The advantages of the finite element method include the efficient and accurate modeling of domains with complex geometric configurations and varying material properties and the capability of accurately analyzing both geometrically and materially nonlinear problems. Shortcomings of the displacement model of FEM have been observed in the analyses of certain classes of problems, such as the large deformation, the treatment of nearly incompressible materials, the bending of thin plates, and the moving boundary problems. The developments of finite elements and corresponding formulations have been a subject of extensive research, which have been summarized in textbooks by Zienkiewicz [1], Bathe [2, 3], and Cook [4].The attempts to extend and generalize the complementary energy methods that have been proposed over the last almost 50 years for small deformation solid mechanics problems have led to the development of several complementary energy principles and corresponding finite element models for nonlinear elastic solid/structural mechanics problems [5–7].Recently, based on the complementary energy principle, a hybrid stress-function element method is proposed by Fu et al. [8] and Cen et al. [9–12] using the mesh distortion immune technique for developing plane 8-node elements. It starts from the principle of minimum complementary energy and employs the fundamental analytical solutions of the Airy stress function as the trial functions (analytical trial function method). So as long as the element edges keep straight, the resulting 8-node element HSF-Q8 can also produce the exact solutions for quadric displacement fields, and even the element shape degenerates into triangle and concave quadrangle.A new concept was introduced by Gao [13], who used the concept of the “base forces” to replace various stress tensors for the description of the stress state at a point. These base forces can be directly obtained from the strain energy. For large deformation problems, when the base forces were adopted, the derivation of basic formulae was simplified by Gao [14, 15]. Based on the concept of the base forces, precise expressions for stiffness and compliance matrices for the FEM were obtained by Gao [13]. The applications of the stiffness matrix to the plane problems of elasticity using the four-edge plane element and the polygonal element were researched by Peng et al. [16]. Using the concept of base forces as state variables, a three-dimensional formulation of base force element method (BFEM) on complementary energy principle was proposed by Peng and Liu [17] for geometrically nonlinear problems. And the new finite element method based on the concept of base forces was called the base force element method (BFEM) by Peng and Liu [17]. A three-dimensional model of base force element method (BFEM) on complementary energy principle was proposed by Liu and Peng [18] for elasticity problems. The application of 2D base force element method (BFEM) to geometrically nonlinear analysis was proposed by Peng et al. [19].The objective of the present research is to present a new two-dimensional formulation of the base forces element method (BFEM) for elasticity problems. A 4-mid-node plane element model of the BFEM for elasticity problems will be derived. In the present formulations of the BFEM, the “base forces” are treated as unknown variables, and the basic equations are constructed by means of the complementary energy principle. The element equilibrium conditions are fulfilled using the Lagrange multiplier method. Explicit expressions for the control equations are provided, and a procedure of the present method is developed. A number of elasticity problems are solved using the present formulation, and the results are compared with corresponding analytical solutions.


## 2. Basic Equations

Consider a two-dimensional domain of solid medium, letxα (α=1,2) denote the Lagrangian coordinate system, where P and Q are the position vectors of a material point before and after deformation, respectively. Two triads for original and current configurations can be defined as
(1)Pα=∂P∂xα,Qα=∂Q∂xα.
Let u denote the displacement of a point; then
(2)u=Q-P.The gradient of displacementuα can be written as
(3)uα=∂u∂xα=Qα-Pα.In order to describe the stress state at a pointQ, a parallelogram with the edges dx1Q1 and dx2Q2 is shown in Figure 1. Define
(4)Tα=dTαdxα+1,dxα⟶0,
where 3=1 for indexes. Quantities Tα (α=1,2) are called the base forces at point Q in the two-dimensional coordinate system xα.Figure 1
Base forces on plane problem.TheTα and uα can directly be expressed by strain energy. Thus, uα is just the conjugate variable of Tα. It can be seen that the mechanics problem can be completely established by means of Tα and uα.According to the definitions of Cauchy stress tensors, the relation between the base forces and Cauchy stress tensors can be given as(5)σ=AQ-1Tα⊗Qα,
where ⊗ is the dyadic symbol, the summation rule is implied, and AQ is the current base area
(6)AQ=|Q1×Q2|.For the isotropic material, the complementary energy density can be expressed as follows:(7)WC=(2E)-1[(1+ν)J2T-νJ1T2]
in which E is Young’s modulus, ν is Poisson’s ratio, and J1T and J2T are the invariants of Tα, so
(8)J1T=AP-1Tα·Pα,J2T=(AP2)-1(Tα·Tβ)pαβ,
where AP is the original base area, pαβ is a metric tensor in the original configurations, and there are the following relations:
(9)AP=|P1×P2|,pαβ=Pα·Pβ.


## 3. Compliance Matrix

Now, consider a 4-mid-node plane element as shown in Figure2. Let I, J, K, and L denote its edges and TI, TJ, TK, and TL the force vectors acting on each of the edges.Figure 2
Four-mid-node plane element.Substituting (8) into (7), the complementary energy of an element can be reduced as follows:
(10)WCe=1+ν2EA[(TI·TJ)pIJ-ν1+ν(TI·PI)2],
where A is the original area of the element, TI is the force vectors acting on the center of the boundary edge I of the element, pIJ is the dot product of position vectors PI and PJ of points I and J, and there are
(11)TI=TI1e1+TI2e2,TJ=TJ1e1+TJ2e2(12)pIJ=PI·PJ,PI=PI1e1+PI2e2,PJ=PJ1e1+PJ2e2
in which e1 and e2 are the unit vectors.From (10), we can obtain the generalized displacement corresponding to TI(13)δI=∂WCe∂TI=CIJ·TJ,
where CIJ is a second-order tensor that is called the compliance matrix
(14)CIJ=1+νEA(pIJU-ν1+νPI⊗PJ)(I,J=1,2,3,4),
where U is the unit tensor, and it can be expressed as
(15)U=e1⊗e1+e2⊗e2.Substituting (12) and (15) into (14), the compliance matrix of an element can be reduced as follows:
(16)CIJ=1+νEA[(11+νPI1PJ1+PI2PJ2)e1⊗e1-ν1+νPI1PJ2e1⊗e2-ν1+νPI2PJ1e2⊗e1+(PI1PJ1+11+νPI2PJ2)e2⊗e2],
or
(17)CIJ=CI1J1e1⊗e1+CI1J2e1⊗e2+CI2J1e2⊗e1+CI2J2e2⊗e2,
or
(18)CIJ=1+νEA×[11+νPI1PJ1+PI2PJ2-ν1+νPI1PJ2-ν1+νPI2PJ1PI1PJ1+11+νPI2PJ2](I,J=1,2,3,4).For a strain problem, it is necessary to replaceE with E/(1-ν2) and ν with ν/(1-ν) in (16)–(18).


## 4. Governing Equations

The total complementary energy of the elastic system can be written as(19)ΠC=∑n(WCe-u-I·TI).According to the complementary energy principle,ΠCne takes stationary value under equilibrium conditions and stress boundary conditions for the real displacement and stress state.The equilibrium conditions can be released by the Lagrange multiplier method. For the stress boundary condition and the corresponding condition of the face forces between elements, computational techniques can be applied to achieve a rotation. Using the Lagrange multiplier method, a new function for an element can be introduced as follows:(20)ΠCe*(T,λ,λ3)=ΠCe(T)+λ(∑I=14TI)+λ3(TI×PI)
in which arbitrary vectors λ, λ3 are the Lagrange multipliers, and λ can be expressed as
(21)λ=λ1e1+λ2e2.For the elastic system, there is(22)ΠC*=∑n[ΠCe*(T,λ,λ3)].
And by means of the modified complementary energy principle
(23)δΠC*=∑n[δΠCe*(T,λ,λ3)]=0.Further, (23) can be expressed as
(24)∂ΠC*(T,λ,λ3)∂T=0,∂ΠC*(T,λ,λ3)∂λ=0,∂ΠC*(T,λ,λ3)∂λ3=0.
Equations (24) are the compatibility equations and displacement boundary conditions for the elastic system. These are the governing equations of the BFEM.An explicit expression of stress for the 4-mid-node element, according to (5), can be written as
(25)σ=A-1TI⊗PI
or
(26)σ=A-1∑I=14[TI1PI1e1⊗e1+TI1PI2e1⊗e2=A-1+TI2PI1e2⊗e1+TI2PI2e2⊗e2]
or
(27)σ=1A[∑I=14TI1PI1∑I=14TI1PI2∑I=14TI2PI1∑I=14TI2PI2].An explicit expression for displacement of node, based on the governing equation of element, can be developed. For a 4-mid-node plane element, the governing equation can be written as(28)δI=∂ΠCe*(T,λ,λ3)∂TI.Substituting (20) into (28) and considering (19), (10), and (14), the explicit expression for displacement of nodes can be written as
(29)δI=CIJ·TJ+λ+λ3ε·PI,
where ε is the alternating tensor and ε can be written as
(30)ε=εαβeα⊗eβ
or
(31)ε=e1⊗e2-e2⊗e1.Substituting (11), (17), (21), and (31) into (29), the explicit expression for displacement of nodes can be written as
(32)δI=(CI1J1TJ1+CI1J2TJ2+λ1+λ3PI2)e1+(CI2J1TJ1+CI2J2TJ2+λ2-λ3PI1)e2
or
(33)δI={CI1J1TJ1+CI1J2TJ2+λ1+λ3PI2CI2J1TJ1+CI2J2TJ2+λ2-λ3PI1}.


## 5. Numerical Examples

To verify the validity and accuracy of the 4-mid-node plane model of base force element method proposed in this paper, a program on the BFEM is made, and a number of the numerical examples on the BFEM for the typical plane elasticity problems are presented in this section. The results obtained from the 2D BFEM are compared with corresponding theoretical solutions. The responses obtained using the traditional displacement finite element method are also given for some problems in order to assess the advantages of the 2D BFEM. All calculations of the traditional displacement finite element method are performed using the four-node isoparametric element (Q4 model).Example 1 (pure shear of rectangular plate).
Consider the problem of a rectangular plate subjected to uniformly distributing shearing forces as shown in Figure3. Parameters about the calculated specimen are taken as a=2.0 m, b=1.0 m, E=105 MPa, ν=0.3, and τ0=1.0 MPa for a plane stress problem. The calculation is done using two different element meshes with the center nodes of each of the edges as shown in Figure 4 successively.Figure 3
Rectangular plate subjected to shearing forces.Two kinds of meshes for the pure shear problem.
(a)
Uniform 4-side elements
(b)
Arbitrary 4-side elementsThe numerical solutions are that the stress component of various elements isτxy=1 for two kinds of different meshes. The results are consistent with the theoretical results [1].Example 2 (thick cylinder with fixed outer surface under internal pressure).
Consider a thick cylinder with inner radiusa and outer radius b under inner pressure pi as shown in Figure 5.Figure 5
Thick cylinder with fixed outer edge under internal pressure.The calculated parameters are taken to bea=1.0m, b=5.0m, E=106MPa, ν=0.3, and pi=1.0MPa for a plane strain problem. The calculation is done using three different element meshes with the center nodes of each of the edges as shown in Figure 6 successively.Three kinds of meshes for cylinder.
(a)
20 × 10 elements
(b)
20 × 7 elements
(c)
20 × 5 elementsThe values of stress componentsσr, σθ, and ur for various radiuses r are listed in Tables 1, 2, and 3, respectively. Comparisons of the results from the theoretical solution [1] and traditional 4-node quadrilateral isoparametric element (Q4 model) are also given in Tables 1, 2, and 3, respectively. The numerical results of the present model are consistent with those of the theoretical solutions and have given very good performance compared with Q4 model.Table 1
The values ofσr.
Meshes
Location
σ
r/(MPa)
Q4
BFEM
Theoretical
2
0
×
10
r
=
2.8000
−0.2066
−0.2072
−0.2069
20
×
7
r
=
3.0000
−0.1914
−0.1926
−0.1919
20
×
5
r
=
3.0000
−0.1912
−0.1935
−0.1919Table 2
The values ofσθ.
Meshes
Location
σ
θ/(MPa)
Q4
BFEM
Theoretical
2
0
×
10
r
=
2.8000
0.0274
0.0253
0.0250
20
×
7
r
=
3.0000
0.0113
0.0107
0.0101
20
×
5
r
=
3.0000
0.0124
0.0116
0.0101Table 3
The values ofur.
Meshes
Location
u
r/(×10-6 m)
Q4
BFEM
Theoretical
2
0
×
10
r
=
1.0000
1.1274
1.1340
1.1346
20
×
7
r
=
1.0000
1.1212
1.1340
1.1346
20
×
5
r
=
1.0000
1.1113
1.1340
1.1346Example 3 (a cantilever under uniformly distributed load).
It consists of a straight cantilever beam under uniformly distributed load, as shown in Figure7. The beam has length L=10m, cross-sectional dimensions b=1m and h=1m, elastic modulus E=109 MPa, and Poisson’s ratio ν=0, respectively. While the applied total load level is p=10 MPa. The calculated specimen was divided into the 4-mid-node base force elements with the center nodes of each of the edges, as shown in Figure 8.Figure 7
Cantilever beam under uniformly distributed load.Meshes of cantilever beam under uniformly distributed load.
(a)
7 × 7 elements
(b)
14 × 7 elements
(c)
20 × 7 elementsThe vertical displacement solutions of BFEM at the tip of the beam with different meshes are listed in Table4. The vertical displacement of BFEM at the tip of the beam is compared with that provided by traditional 4-node quadrilateral isoparametric element (Q4 model) and the analytical solution.Table 4
Vertical displacement of the cantilever beam under uniformly distributed load.
Meshes
v (×10-4)
BFEM
Q4 model
Analytical
2
0
×
7
1.4441
1.3447
1.5000
14
×
7
1.4041
1.2062
1.5000
7
×
7
1.2802
0.7528
1.5000Table4 shows that the results of the present formulation agree well with those of the analytical solution, and the 4-mid-node element of BFEM has given very good performance compared with Q4 model for the large aspect ratio of elements.Example 4 (a cantilever with large aspect ratio under a concentrated end load).
A clamped cantilever with a concentrated end load is carried out as shown in Figure9. The beam has length L=5 m, cross-sectional dimensions b=1 m and h=0.1 m, elastic modulus E=109 N/m2, and Poisson’s ratio ν=0. While the applied total load level is P=1 N. The calculated specimen was divided into the 4-mid-node base force elements with the center nodes of each of the edges, as shown in Figure 10.Figure 9
Cantilever beam under a concentrated end load.Meshes of cantilever beam under a concentrated end load.
(a)
60 × 9 elements
(b)
30 × 9 elements
(c)
15 × 9 elementsThe vertical displacement of BFEM at the tip of the beam is compared with that provided by traditional 4-node quadrilateral isoparametric element (Q4 model) and the analytical solution in Table5. The results show that the 4-mid-node element of BFEM has given very good performance compared with Q4 model for the large aspect ratio of elements.Table 5
Vertical displacement of the cantilever beam with a concentrated end load.
Meshes
v (×10-4)
BFEM
Q4 model
Analytical
6
0
×
9
−5.0633
−3.7122
−5.0000
30
×
9
−5.0623
−2.0935
−5.0000
15
×
9
−5.0581
−0.7629
−5.0000Example 5 (a cantilever with mesh distortion under a concentrated end load).
A clamped cantilever with a concentrated end load is carried out as shown in Figure9. The beam has length L=10 m, cross-sectional dimensions b=1 m and h=1 m, elastic modulus E=106 Pa, and Poisson’s ratio ν=0, while the applied total load level is P=1 N. The calculated specimen was divided into the 4-mid-node base force elements with the center nodes of each of the edges, as shown in Figure 11.Distortion meshes of cantilever beam.
(a)
Mesh distortion parametersd=0
(b)
Mesh distortion parametersd=0.5
(c)
Mesh distortion parametersd=1.0
(d)
Mesh distortion parametersd=1.5The normalized vertical displacement of BFEM at the tip of the beam is compared with that provided by traditional 4-node quadrilateral isoparametric element (Q4 model) and the analytical solution in Table6. The results show that the 4-mid-node element of BFEM has given very good performance compared with Q4 model for the distortion mesh problems.Table 6
Normalized vertical displacement with the effect of mesh sensitivity.
Mesh distortion parametersd
0.0
0.5
1.0
1.5
BFEM
1.0058
1.0014
0.9868
0.9654
Q4
0.0744
0.0463
0.0380
0.0361


## 6. Conclusions

In this paper, a 4-mid-node plane model of base force element method (BFEM) on complementary energy principle is proposed for elasticity problems. Based on the complementary energy principle, the equilibrium conditions are released by the Lagrange multiplier method, and a modified complementary energy principle described by the base forces is obtained. The 4-mid-node plane element model is derived by assuming that the stress is uniformly distributed on each edge of a plane element. A BFEM procedure is developed using MATLAB language. The chief features of the method are that the model does not introduce an interpolating function and can be used in any coordinate system and is not necessary to introduce Gauss’ integral for calculating the stiffness coefficient at a point.The present model is applied for the solution of some elasticity problems. The comparisons are made with the corresponding analytical solutions, and a good agreement of results is observed. A series of numerical experiments are performed in order to assess the relative performances of the present model; based on these studies, it is found that the 4-mid-node plane model of BFEM with complementary energy principle provides reliable predictions for elasticity problems. A good performance of the 4-mid-node plane element of BFEM is observed for some problems.The 4-mid-node plane element model of BFEM is also compared with the standard displacement method. It has been demonstrated by the results of numerical calculation that the new BFEM formulations proposed in this paper are effective and accurate. The 4-mid-node plane element model can be used efficiently for the elasticity analysis with increasing element aspect ratios and distortion meshes. The results presented here revealed a better overall performance of the BFEM in calculations. These results provide further confirmation that the BFEM can be used successfully and efficiently in the structural analysis, thus serving as an alternative method, and certain applications, being the method of choice for the analysis.



---

*Source: 706759-2013-09-24.xml*