Télécharger la plaquette Cast3M


fr en

Télécharger errare.eso

Retour à la liste

Numérotation des lignes : 

errare

  1. C ERRARE    SOURCE    PV        22/04/25    21:15:06     11344          
  2.       SUBROUTINE ERRARE(IPMODL,IPCHE1,IPCHE2,XERR,IPCHRR)
  3. *_______________________________________________________________________
  4. *
  5. *      CALCUL D'ERREUR SUR LES CONTRAINTES
  6. *
  7. *      METHODE ZIENKIEWICZ - VERSION BARZIC ET RICHARD
  8. *
  9. *_______________________________________________________________________
  10. *
  11. *    ENTREES :
  12. *    ---------
  13. *        IPMODL   POINTEUR SUR UN MMODEL
  14. *        IPCHE1   POINTEUR SUR UN MCHAML DE CONTRAINTES
  15. *        IPCHE2   POINTEUR SUR UN MCHAML DE MATERIAU
  16. *
  17. *
  18. *    SORTIES :
  19. *    ---------
  20. *        XERR     ERREUR GENERALE SCALAIRE
  21. *        IPCHRR   MCHAML D'ERREUR PAR ELEMENT
  22. *                 = 0 SI ERREUR
  23. *
  24. *    PASSAGE AUX NOUVEAUX CHAMELEM PAR JM CAMPENON LE 02/91
  25. *
  26. *------------------------------------------------------------------
  27. *
  28.       IMPLICIT INTEGER(I-N)
  29.       IMPLICIT REAL*8 (X)
  30. *
  31. -INC PPARAM
  32. -INC CCOPTIO
  33. -INC SMCOORD
  34.       real*8 flo
  35. *
  36.       XERR = -1.E+9
  37.       IPCHRR=0
  38.       SEGACT,MCOORD
  39.  
  40.       CALL SIGMOY(IPMODL,IPCHE1,IPOSGM)
  41.       IF (IPOSGM.EQ.0) GOTO 310
  42. *
  43.       CALL ADCHEL(IPCHE1,IPOSGM,IPODSI,-1)
  44.       IF (IPODSI.EQ.0) GOTO 300
  45. *
  46. *   CHAMP DE VARIABLES INTERNES : IPCH3 = 0 ici (a lire dans ERRE ?)
  47. *
  48.       IPCHE3 = 0
  49. *
  50.       CALL ELAS1(IPMODL,IPCHE1,IPCHE2,IPCHE3,IPOEPS,IRET)
  51.       IF (IRET.NE.1) GOTO 290
  52. *
  53.       CALL ELAS1(IPMODL,IPODSI,IPCHE2,IPCHE3,IPODEP,IRET)
  54.       IF (IRET.NE.1) GOTO 280
  55.  
  56. *
  57.       CALL ENERCA(IPMODL,IPCHE1,IPOEPS,IPOENE)
  58.       IF (IPOENE.EQ.0) GOTO 270
  59. *
  60.       CALL ENERCA(IPMODL,IPODSI,IPODEP,IPODEN)
  61.       IF (IPODEN.EQ.0) GOTO 260
  62. *
  63.       CALL INTGCA(IPMODL,IPOENE,0,1,IRET,XRET,IPORET)
  64.       IF (IRET.EQ.0) GOTO 250
  65. *
  66.       CALL INTGCA(IPMODL,IPODEN,0,1,IRET,XDRE,IPODRE)
  67.       IF (IRET.EQ.0) GOTO 240
  68. *
  69.       XERR=(XDRE/(XRET+XDRE))**(0.5D0)
  70.  
  71.       CALL ADCHEL(IPORET,IPODRE,IPOSOM,1)
  72.  
  73. *     Pour la PUISSANCE : IOPERA=1
  74.       IOPERA= 1
  75.       IARGU = 1
  76.       I1    =-1
  77.       CALL OPCHE1(IPOSOM,IOPERA,IARGU,I1,FLO,IPOFAC,IRET)
  78.       IF(IRET.EQ.0) GOTO 220
  79. *
  80.       IPLMOT = -1
  81.       IPLREE = -1
  82.       CALL MUCHSC(IPMODL,IPOFAC,IPODRE,IPOCAR,
  83.      &                   IPLMOT,IPLMOT,IPLMOT,IPLREE,1)
  84.       IF (IPOCAR.EQ.0) GOTO 210
  85. *
  86.       I1    =1
  87.       CALL OPCHE1(IPOCAR,IOPERA,IARGU,I1,FLO,IPCHRR,IRET)
  88.       IF (IRET.NE.0) GOTO 200
  89. *
  90.       IPCHRR=0
  91.  200  CALL DTCHAM(IPOCAR)
  92.  210  CALL DTCHAM(IPOFAC)
  93.  220  CALL DTCHAM(IPOSOM)
  94.  230  CALL DTCHAM(IPODRE)
  95.  240  CALL DTCHAM(IPORET)
  96.  250  CALL DTCHAM(IPODEN)
  97.  260  CALL DTCHAM(IPOENE)
  98.  270  CALL DTCHAM(IPODEP)
  99.  280  CALL DTCHAM(IPOEPS)
  100.  290  CALL DTCHAM(IPODSI)
  101.  300  CALL DTCHAM(IPOSGM)
  102.  310  CONTINUE
  103.       SEGDES,MCOORD
  104.  
  105. *      RETURN
  106.       END
  107.  
  108.  
  109.  
  110.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales