Nonlinear finite elements/Homework 3/Solutions
Problem: The Patch Test
Taylor, R.L., Simo, J.C., Zienkiewicz, O.C., and Chan, A.C.H, 1986, The patch test - a condition for assessing FEM convergence, International Journal for Numerical Methods in Engineering, 22}, pp. 39-62.
What is the patch test and why is it used?
- What is the patch test?
The standard patch test involves setting up a mesh of elements as shown in Figure 1(a) and checking the predicted solution against a known solution. This mesh is called a patch.
The elements should be arbitrary quadrilaterals rather than rectangles. No body forces should be applied. Material properties should be uniform and linear elastic.
The displacements of nodes , , , and are prescribed. The prescribed value depends on the order of the polynomial that is to be represented exactly by the elements. A similar process is followed for three-dimensional elements.
Suppose we wish to have both constant and linear displacements be represented exactly by the elements. Then the prescribed displacement field is
where () and () are constants set by the user. These constants should be nonzero.
Therefore, the prescribed nodal displacements are
where are the coordinates of node .
The finite element solution for these prescribed displacements should be given by equation (2) throughout the patch.
Hence the displacements at nodes , , , and should be given by equation (3)}.
The element strains should be constant and be given by
The element stresses should also be constant.
These results should be satisfied with a high degree of precision.
The paper talks about an additional test that can be used to assess stability (see Figure 1(b)).
- Why is it used?
For a finite element analysis, we seek elements (and therefore shape functions) that behave such that as the mesh is refined the approximate solution converges to the exact solution.
The patch test is a necessary condition for convergence. If a finite element solution is to converge with refinement, the finite elements must pass the patch test.
In practice, the patch test checks whether an element passes the completeness test.
Suppose the approximate solution over an element is given by
Suppose that we would like the element to represent exactly linear and constant functions. Then, in two dimensions, the shape functions () are said to be complete if
where , and are arbitrary constants.
This means that a complete element can represent exactly all rigid motions (constant displacements) and all constant strain states.
The patch test is also used to check the correct implementation of a finite element program.
What is the consistency requirement?What is the stability condition?
The term consistency condition comes from the finite difference literature where it refers to the formal consistency of the difference equations with the correct differential equation.
In the paper, the consistency requirement implies that as the size of the elements tends to zero, the approximate equation will represent the governing differential equation exactly.
By stability, the authors mean that the solution of the finite element system of equations is unique. For linear systems, this means that the stiffness matrix has to be nonsingular.
What do the authors mean by convergence?
Let the finite element approximation be
where is the vector of shape functions, and is the vector of nodal displacements.
By convergence, the authors mean that the approximate solution as , where is the exact solution, and is the element size.
List the 2-D structural elements available in ANSYS and describe
what shape functions and degrees of freedom each uses.
Listed below are some example of the 2-D structural elements used by ANSYS. The descriptions and the shape functions of each of these elements are from the ANSYS Element Manual.
- PLANE42: 2-D 4-Node Structural Solid
- The PLANE42 element can be used either as a plane element (plane stress or plane strain) or as an axisymmetric element. The element is defined by four nodes (see Figure 4) with two translational degrees of freedom at each node. The shape functions for this element are
- where and .
- However, if incompatible elements are used (KEYOPT(2)=0), then the approximate solution is of the form
- where nodes and are incompatible nodes. The forces at the incompatible nodes are zero. Hence the element stiffness matrix can be statically condensed and and do not have to correspond to any physical nodes.
- PLANE82: 2-D 8-Node Structural Solid
- The PLANE82 element is defined by eight nodes (see Figure 5) with two translational degrees of freedom at each node. The element may be used as a plane elementor as an axisymmetric element. The shape functions of this element are
- where and .
- PLANE182: 2-D 4-Node Structural Solid
- PLANE182 is a four-noded quadrilateral element and has the same displacement shape functions as PLANE42. In addition, the element allows for pressure degrees of freedom if KEYOPT(6) = 1 is chosen (mixed u-p formulation). The hydrostatic pressure has an interpolation function one order lower than the displacement interpolation. If KEYOPT(1) = 0 is chosen, a method with full integration is used. One pressure degree of freedom isused per element with this option. If KEYOPT(1) = 1 is chosen, a method with selective reduced integration (and hourglass control) is used in the computations. If KEYOPT(1) = 2 is chosen, an enhanced strain method isused and the pressure is interpolated with linear shape functions. There are three pressure degrees of freedom in the element.
- PLANE183 - 2-D 8-Node Structural Solid
- PLANE183 is a two-dimensional eight-node element similar to PLANE82. This element can also include a mixed u-p formulationwhen KEYOPT(6) is set to 1.
Check whether each of these elements passes the patch test. Show details of your calculations.
Patch test for PLANE182 elements
In this test we assume that the displacements on the patch have the form
Then, assuming plane strain, the strains in the patch are
The stresses are given by
If and , the stresses are
If and , the stresses are
Let us apply the patch test to the patch shown in Figure 6.
The exact solution for the displacements at the nodes is
We apply the exact displacements at nodes 1, 2, 3, and 4 and compute the remaining displacements and the stresses and strains using ANSYS. The ANSYS script we have used is shown below:
/prep7 ! ! Element type ! et,1,182 ! ! Keyopt(1) - 0 = full integration, 1 = uniform reduced integration ! 2 = enhanced strain ! !keyopt,1,1,0 !keyopt,1,1,1 keyopt,1,1,2 ! ! Keyopt(3) - 0 = plane stress, 2 = plane strain ! keyopt,1,3,2 ! ! Keyopt(6) - 0 = disp (u) formulation, 1 = mixed (u-p) formulation ! !keyopt,1,6,0 keyopt,1,6,1 ! ! Material properties ! mp,ex,1,1000.0 mp,prxy,1,0.499999 ! ! Create nodes ! n,1,0,0 n,2,2,0 n,3,2,3 n,4,0,2 n,5,0.4,0.4 n,6,1.4,0.6 n,7,1.5,2 n,8,0.3,1.6 ! ! Create elements ! e, 1, 2, 6, 5 e, 2, 3, 7, 6 e, 7, 3, 4, 8 e, 1, 5, 8, 4 e, 5, 6, 7, 8 ! ! Apply displacements to the boundary nodes ! d,1,ux,0.001 d,1,uy,-0.001 d,2,ux,0.005 d,2,uy,-0.005 d,3,ux,0.014 d,3,uy,-0.014 d,4,ux,0.007 d,4,uy,-0.007 fini ! ! Solve (Ku = f) ! /solu solve fini ! ! Write out results ! /post1 /output,plane182Patch,out prnsol, u, comp presol, s, comp presol, epel, comp /output,, fini
The results from ANSYS are compared with the exact values in the tables below.
These results clearly show that the element is passing the patch test. Similar results are obtained by changing the keyopt values in the input file.
Patch test for PLANE42 The patch test (Test B) is applied to the model shown in Figure 5a in the paper. The same solution given in the handout is considered for this analysis.
Listed below are the inputs for the 4-node element. To analyze the incompressible materials, is set to 0.4999 and the element type is changed to PLANE182. The keyopt command should be changed to reflect the options available of PLANE182.
/prep7 !change this line to et,1,182 for PLANE182 et,1,42 !change this line to keyopt,1,2,0 to include the extra displacement keyopt,1,2,1 mp,ex,1,1e3 mp,prxy,1,0.3 k,1,0,0 k,2,2,0 k,3,2,3 k,4,0,2 k,5,0.4,0.4 k,6,1.4,0.6 k,7,1.5,2 k,8,0.3,1.6 !create lines and set their number of division to 1 for each line l,1,2,1 l,2,3,1 l,3,4,1 l,4,1,1 l,5,6,1 l,6,7,1 l,7,8,1 l,8,5,1 l,1,5,1 l,2,6,1 l,3,7,1 l,4,8,1 a,1,2,6,5 a,3,4,8,7 a,1,5,8,4 a,5,6,7,8 a,2,3,7,6 !AMAP is used instead of AMESH in this case to ensure the quad elements AMAP,1,1,2,6,5 AMAP,2,3,4,8,7 AMAP,3,1,5,8,4 AMAP,4,5,6,7,8 AMAP,5,6,2,3,7 !apply displacement conditions to the surrounding nodes dk,1,all,0 dk,2,ux,0.004 dk,2,uy,0 dk,3,ux,0.004 dk,3,uy,-0.0018 dk,4,ux,0 dk,4,uy,-0.0012 fini
Patch test for PLANE82
For an 8-node element, the similar model is created. However, the displacements must also be applied to the midnodes. Only the inputs for PLANE183 are shown here. To run the analysis using PLANE82, the ET line should be changed to et,1,82, and the keyopt line should be changed to reflect the options available for PLANE82. The input for the model using PLANE183 is shown below
/prep7 !change this line to et,1,82 for PLANE82 et,1,183 !sets the element for using the mixed formulation (u-p). keyopt,1,6,1 mp,ex,1,1e3 mp,prxy,1,0.4999 k,1,0,0 k,2,2,0 k,3,2,3 k,4,0,2 k,5,0.4,0.4 k,6,1.4,0.6 k,7,1.5,2 k,8,0.3,1.6 !create lines and set their number of division to 1 for each line l,1,2,1 l,2,3,1 l,3,4,1 l,4,1,1 l,5,6,1 l,6,7,1 l,7,8,1 l,8,5,1 l,1,5,1 l,2,6,1 l,3,7,1 l,4,8,1 a,1,2,6,5 a,3,4,8,7 a,1,5,8,4 a,5,6,7,8 amesh,all !this area is meshed separately to force a quad mesh a,2,3,7,6 AMAP,5,6,2,3,7 !apply displacement conditions to the surrounding nodes d,1,all,0 d,2,ux,0.004 d,2,uy,0 d,9,ux,0.004 d,9,uy,-0.0018 d,10,ux,0 d,10,uy,-0.0012 d,3,ux,0.002 d,3,uy,0 d,20,ux,0.004 d,20,uy,-0.0006*1.5 d,11,ux,0.002 d,11,uy,-0.0006*2.5 d,18,ux,0 d,18,uy,-0.0006 fini
Results The displacements obtained from ANSYS are compared to (10). In order pass the patch test, the displacement at each node must match (10), and stresses have to be constant everywhere.
Similarly for PLANE182 and PLANE183, the results shown in the table below use to simulate the incompressible material. The verification of the results is also similar to the previous cases.
Hence, we conclude that these four elements used by ANSYS, PLANE42, PLANE82, PLANE182, and PLANE183 passed the patch test.