Télécharger errare.eso

Retour à la liste

Numérotation des lignes :

  1. C ERRARE SOURCE FANDEUR 12/04/11 21:15:18 7350
  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 CCOPTIO
  32. *
  33. IPCHRR=0
  34. CALL SIGMOY(IPMODL,IPCHE1,IPOSGM)
  35. IF (IPOSGM.EQ.0) GOTO 310
  36. *
  37. CALL ADCHEL(IPCHE1,IPOSGM,IPODSI,-1)
  38. IF (IPODSI.EQ.0) GOTO 300
  39. *
  40. * CHAMP DE VARIABLES INTERNES : IPCH3 = 0 ici (a lire dans ERRE ?)
  41. *
  42. IPCHE3 = 0
  43. *
  44. CALL ELAS1(IPMODL,IPCHE1,IPCHE2,IPCHE3,IPOEPS,IRET)
  45. IF (IRET.NE.1) GOTO 290
  46. *
  47. CALL ELAS1(IPMODL,IPODSI,IPCHE2,IPCHE3,IPODEP,IRET)
  48. IF (IRET.NE.1) GOTO 280
  49. *
  50. CALL ENERCA(IPMODL,IPCHE1,IPOEPS,IPOENE)
  51. IF (IPOENE.EQ.0) GOTO 270
  52. *
  53. CALL ENERCA(IPMODL,IPODSI,IPODEP,IPODEN)
  54. IF (IPODEN.EQ.0) GOTO 260
  55. *
  56. CALL INTGCA(IPMODL,IPOENE,0,1,IRET,XRET,IPORET)
  57. IF (IRET.EQ.0) GOTO 250
  58. *
  59. CALL INTGCA(IPMODL,IPODEN,0,1,IRET,XDRE,IPODRE)
  60. IF (IRET.EQ.0) GOTO 240
  61. *
  62. XERR=(XDRE/(XRET+XDRE))**(.5)
  63. CALL ADCHEL(IPORET,IPODRE,IPOSOM,1)
  64. CALL PUCHEL(IPOSOM,-1,IPOFAC,IRET)
  65. IF(IRET.EQ.0) GOTO 220
  66. *
  67. IPLMOT = -1
  68. IPLREE = -1
  69. CALL MUCHSC(IPMODL,IPOFAC,IPODRE,IPOCAR,
  70. & IPLMOT,IPLMOT,IPLMOT,IPLREE,1)
  71. IF (IPOCAR.EQ.0) GOTO 210
  72. *
  73. CALL PUCHEL(IPOCAR,1,IPCHRR,IRET)
  74. IF (IRET.NE.0) GOTO 200
  75. *
  76. IPCHRR=0
  77. 200 CALL DTCHAM(IPOCAR)
  78. 210 CALL DTCHAM(IPOFAC)
  79. 220 CALL DTCHAM(IPOSOM)
  80. 230 CALL DTCHAM(IPODRE)
  81. 240 CALL DTCHAM(IPORET)
  82. 250 CALL DTCHAM(IPODEN)
  83. 260 CALL DTCHAM(IPOENE)
  84. 270 CALL DTCHAM(IPODEP)
  85. 280 CALL DTCHAM(IPOEPS)
  86. 290 CALL DTCHAM(IPODSI)
  87. 300 CALL DTCHAM(IPOSGM)
  88. 310 RETURN
  89. END
  90.  
  91.  
  92.  

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