Modulefpreviousupnextcontentsindex[BIG][Normal][small]
Suiv.: 3 Manipulation de structures et divers Sup.: 2 Manipulation de solutions Préc.: 2.9 Extrema d'une S.D. B Index Table des matières


Normes et résidus

Normes relatives à une S.D. B

Comparer la solution exacte avec la solution calculée consiste à calculer, en chaque noeud:

est la solution exacte et est la solution calculée et a imprimer, pour chaque degré de liberté:

en indiquant de plus le numéro et les coordonnées du noeud où le maximum est atteint.

Le préprocesseur NORMXX permet, dans le cas où on connaît analytiquement la solution d'un problème, de comparer la solution calculée avec la solution exacte. Il appelle le module NORME:

      SUBROUTINE NORME (M,XM,DM,NFMAIL,NIMAIL,NFCOOR,NICOOR,NFB,NIB,
     +   NFBS,NIBS,INDICB,NSM,FONINT,SOLEX,DSOLEX)
C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C    BUT  : IMPRIMER LA SOLUTION EXACTE, LA SOLUTION CALCULEE,
C    ---    LES DIFFERENCES ABSOLUES ET RELATIVES ENTRE ELLES,
C           LES ERREURS L1,L2 ET MAX
C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

Selon la valeur de FONINT, on doit écrire les fonctions SOLEX ou DSOLEX qui ont la forme suivante:

      FUNCTION SOLEX(I,X,Y,Z)

      DOUBLE PRECISION FUNCTION DSOLEX(I,X,Y,Z)

où I est le numéro du degré de liberté du noeud de coordonnées X, Y et Z.

Normes relatives à une S.D. TAE

Le préprocesseur NORMXX permet, dans le cas où on connaît analytiquement les contraintes d'un problème en élasticité, de comparer les contraintes calculées avec cette solution exacte avec les mêmes calculs de normes que pour le cas d'une S.D. B, cas vu ci-dessus. Il appelle le module NORTAE:

      SUBROUTINE NORTAE(M,XM,DM,NFTAE,NITAE,NFTAES,NITAES,INDICB,
     +                  FONINT,SOLEX,DSOLEX,NSM,NC1)
C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C    BUT  : IMPRIMER LA SOLUTION EXACTE, LA SOLUTION CALCULEE,
C    ---    LES DIFFERENCES ABSOLUES ET RELATIVES ENTRE ELLES,
C           LES ERREURS L1,L2 ET MAX  POUR UNE SOLUTION ISSUES DE TAES
C ..................................................................

Selon la valeur de FONINT, on doit écrire les fonctions SOLEX ou DSOLEX qui ont la forme suivante:

      FUNCTION SOLEX(I,X,Y,Z)

      DOUBLE PRECISION FUNCTION DSOLEX(I,X,Y,Z)

où I est le numéro de la composante des contraintes traitée et X,Y et Z sont les coordonnées du point de calcul de cette contrainte.

Norme de l'énergie relative à une S.D. B

Le préprocesseur NORMXX permet, dans le cas où on connaît analytiquement la solution d'un problème, de calculer la norme de l'énergie entre cette solution et la solution obtenue lors de la résolution du problème.

Il appelle le module ENERGI:

      SUBROUTINE ENERGI (M,XM,DM,NFMAIL,NIMAIL,NFCOOR,NICOOR,NFB,NIB,
     +   NFTAE,NITAE,NSM,FONINT,SOLEX,DSOLEX)
C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C    BUT  : CALCULER L'ERREUR EN NORME DE L'ENERGIE
C    ---    |A ( U-UH ; U-UH )|
C ..................................................................

Selon la valeur de FONINT, on doit écrire les fonctions SOLEX ou DSOLEX qui ont la forme suivante:

      FUNCTION SOLEX(I,X,Y,Z)

      DOUBLE PRECISION FUNCTION DSOLEX(I,X,Y,Z)

où I est le numéro du degré de liberté du noeud de coordonnées X, Y et Z.

Erreur de résolution sur une S.D. B

Le module ERREUR permet de calculer l'erreur de résolution d'un système linéaire à partir de la solution de ce système et des matrices et seconds membres élémentaires.

Pour les degrés de liberté dont les valeurs ont été bloquées, le résultat de ce calcul indique la réaction.

      SUBROUTINE ERREUR(M,NFTAE,NITAE,NFB,NIB,NFNDL1,NINDL1,NFBER,
     +                  NIBER,NTBER,NOT1,NOT2)
C  +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C  BUT : CALCUL DU VECTEUR ERREUR.
C  ---
C  +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

Si NTBER est non nul, on doit donner le descriptif associé (cf. ci-dessus).


Modulefpreviousupnextcontentsindex[BIG][Normal][small]
Suiv.: 3 Manipulation de structures et divers Sup.: 2 Manipulation de solutions Préc.: 2.9 Extrema d'une S.D. B Index Table des matières