12.10      Error Recovery

 

RECOVER_ERROR

 

            RECOVER_ERROR   step_number = ns, etc....

 

Determine the magnitude of the error at any given time step.

 

Note    Variable Name        Type        Default      Description

 

(1)        Element_group(s)   string          [none]     Element group(s) for which the error is to be
                                                                                   estimated.  Name(s) must be enclosed in
                                                                                   quotation marks.

 

            Step_number         integer            [0]        Step number at which error recovery is to be
                                                                                   initiated.

 

            Frequency              integer            [0]        Error recovery frequency

 

 

EXAMPLE

            Recover_error /

                        element_group(s) = "group_1", "group_2" /

                        step_number = 1 /

                        frequency = 1

 

Notes/

(1)        List all element group names or numbers for which the error is to be estimated.

 

(2)        The finite element provides an approximation to the exact solution of a mathematically well-posed problem.  The difference between the exact and approximate solutions decreases as the size of the subdivision h gets smaller, or as the order of the interpolation polynomial functions used increases.  In order to determine exactly the magnitude of the error at a given stage of subdivision the exact solution must be available, a situation not common in practice.

 

            Error estimates can be obtained by projecting the field quantities onto the nodes to obtain a continuous field and by comparing the projected quantities with the elemental values, e.g., the error in stresses in continuum type problems defined as

 

                                                                                                                      (1)

 

where  = "exact solution" and  = "finite element" stress, can be estimated with good accuracy as:

 

                                                                                                                     (2)

 


 

where  = "projected" finite element stress.  The projection is performed by using the same interpolation functions as used for representing the displacement field , viz.,

 

                                                                                                                            (3)

 

where   = vector of nodal displacements, and N = nodal interpolation function.  Then

 

                                                                                                                        (4)

 

The approximating equation is achieved by a weighted residual requirement for equality between and  and , viz.,

 

                                                                                                      (5)

 

giving rise to the following matrix problem:

 

                                                                                                       (6)

 

where

 

                                                                                                               (7)

 

The above computation is made particularly simple if a lumped or diagonal approximation is made to A by using nodal quadrature.  With both  and  fields available, the error in stress can be estimated with good accuracy as:

 

                                                                                                                     (8)

 

The error computed in this manner is an a posterior  error estimate, since it can only be  computed after the solution of the problem has been obtained. The following error norms are computed:

 

            = L2 norm of stress error                      (9)

 

            = RMS error in stress                                                   (10)

 

            These error norms are evaluated over the whole domain and over each finite element subdomain.  We note that:

 

                                                                                                              (11)

 

where i refers to the individual finite element subdomains.  A posterior error estimates can be used to refine the finite element approximation locally to achieve results of a given desired accuracy economically and optimally.  For that purpose, the error is projected onto the nodes and may be displayed graphically.

 


 

 

Notes . .