Finite Element Analysis. Barna Szabó

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

Читать онлайн книгу Finite Element Analysis - Barna Szabó страница 37

Автор:
Жанр:
Серия:
Издательство:
Finite Element Analysis - Barna Szabó

Скачать книгу

(QoIs) such as, for example, the maximum and minimum values of u or u prime on upper I equals left-parenthesis 0 comma script l right-parenthesis. 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

      (1.95)pi left-parenthesis u Subscript upper F upper E Baseline right-parenthesis minus pi left-parenthesis u Subscript upper E upper X Baseline right-parenthesis almost-equals StartFraction script upper C squared Over upper N Superscript 2 beta Baseline EndFraction

      where script upper C equals Overscript def Endscripts upper C double-vertical-bar u Subscript upper E upper X Baseline double-vertical-bar Subscript upper E left-parenthesis upper I right-parenthesis. There are three unknowns: pi left-parenthesis u Subscript upper E upper X Baseline right-parenthesis, script upper C and β. Assume that we have a sequence of solutions corresponding to the hierarchic sequence of finite element spaces upper S Subscript i minus 2 Baseline subset-of upper S Subscript i minus 1 Baseline subset-of upper S Subscript i. Let us denote the corresponding computed potential energy values by pi Subscript i minus 2, pi Subscript i minus 1, πi and the degrees of freedom by upper N Subscript i minus 2, upper N Subscript i minus 1, Ni. We will denote the estimate for pi left-parenthesis u Subscript upper E upper X Baseline right-parenthesis by π∞. With this notation we have:

      and, repeating with i minus 1 substituted for i, it is possible to eliminate 2 beta to obtain:

      where

upper Q equals log StartFraction upper N Subscript i minus 1 Baseline Over upper N Subscript i Baseline EndFraction left-parenthesis log StartFraction upper N Subscript i minus 2 Baseline Over upper N Subscript i minus 1 Baseline EndFraction right-parenthesis Superscript negative 1 dot

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

      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

Скачать книгу