Finite Element Analysis. Barna Szabó

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

СКАЧАТЬ (QoIs) such as, for example, the maximum and minimum values of u or on . Since finite element solutions are approximations to an exact solution, it is not sufficient to report the value of a QoI computed from the finite element solution. It is also necessary to provide an estimate of the relative error in the QoI, or present evidence that the relative error in the QoI is not greater than an acceptable value.

      In this section we will use the a priori estimates described in Section 1.5.2 to obtain a posteriori estimates of error in energy norm. It is possible to obtain very accurate estimates for a large class of problems which includes most problems of practical interest.

      Error estimation based on extrapolation

For most practical problems the estimate (1.92) is sufficiently sharp so that the less than or equal sign (≤) can be replaced by the approximately equal sign (≈) and this a priori estimate can be used in an a posteriori fashion.

The computed values of the potential energy corresponding to a sequence of finite element spaces can be used for estimating the error in energy norm by extrapolation. Sequences of finite element spaces that have this property are called hierarchic sequences. By Theorem 1.5 and eq. (1.92) we have:

      (1.95)

      where . There are three unknowns: , and β. Assume that we have a sequence of solutions corresponding to the hierarchic sequence of finite element spaces . Let us denote the corresponding computed potential energy values by , , πi and the degrees of freedom by , , Ni. We will denote the estimate for by π∞. With this notation we have:

(1.96)

(1.97)

On dividing eq. (1.96) with eq. (1.97) and taking the logarithm we get

(1.98)

      and, repeating with substituted for i, it is possible to eliminate to obtain:

(1.99)

      where

Equation (1.99) can be solved for π∞ to obtain an estimate for the exact value of the potential energy.

      The relative error in energy norm corresponding to the ith finite element solution in the sequence is estimated from

(1.100)

      Usually the percent relative error is reported. This estimator has been tested against the known exact solution of many problems of various smoothness. The results have shown that it works well for a wide range of problems, including most problems of practical interest; however, it cannot be guaranteed to work well for all conceivable problems. For example, this method would fail if the exact solution would happen to be energy‐orthogonal to all basis functions associated with (say) odd values of i.

      Remark 1.12 From equation (1.92) we get

СКАЧАТЬ