Finite Element Analysis. Barna Szabó

Чтение книги онлайн.

СКАЧАТЬ EndRoot EndFraction dot"/>

      The error of approximation over the entire domain is:

(1.130)

      By Theorem 1.2, the exact value of the potential energy is

      (1.131)

      and the relative error in energy norm on the entire domain is:

      (1.132)

Remark 1.13 In estimating the local error we used . It can be shown that in the special case of this problem () this relationship holds and therefore using the equal sign in eq. (1.128) is justified. In the general case () however, and eq. (1.128) will be an estimate of the local error in the finite element solution. Therefore the equal sign in eq. (1.128) has to be replaced by the approximately equal () sign and the first equal sign in eq. (1.130) has to be replaced with the less or equal () sign.

Example 1.13 This example illustrates the distribution of the relative error among the elements for a fixed mesh and polynomial degree for selected fractional values of α. Uniform mesh on the domain with and for is used. The exact solution for is shown in Fig. 1.12. The percent relative error in energy norm associated with the kth element, given by eq. (1.129), is shown in Table 1.4 and the relative error for the entire domain is shown in the last column.

      It is seen that for all values of α the maximum error is associated with the first element.

Example 1.14 This example illustrates the distribution of the relative error among the elements for a fixed mesh and polynomial degree for selected integer values of α. Uniform mesh on the domain with and for is used. The percent relative error in energy norm associated with the kth element, given by eq. (1.129), is shown in Table 1.5 and the relative error for the entire domain is shown in the last column.

      The error of approximation for is zero. This follows directly from Theorem 1.4: The exact solution is a polynomial of degree 2. Therefore it lies in the finite element space and hence the finite element solution is the same as the exact solution.

Remark 1.14 In the foregoing discussion it was tacitly assumed that all data computed by numerical integration were accurate and the coefficient matrices of the linear equations were such that small changes in the right‐hand‐side vector produce small changes in the solution vector. This happens when the condition number of the coefficient matrix is reasonably small. In the finite element method the condition number depends on the choice of the shape functions, the mapping functions and the mesh. In one‐dimensional setting the mapping is linear and the shape functions are energy‐orthogonal, therefore round‐off errors are not significant. This is not the case in two and three dimensions, however.

Table 1.4 Example: Element‐by‐element and total relative errors in energy norm (percent) for selected fractional values of α.

	Element number
α	1	2	3	4	СКАЧАТЬ В начало < 40 41 42 43 44 45 46 47 48 49 > В конец