Télécharger energi.eso

Retour à la liste

Numérotation des lignes :

  1. C ENERGI SOURCE CB215821 19/08/01 21:15:45 10279
  2. SUBROUTINE ENERGI
  3. *_______________________________________________________________________
  4. *
  5. * OPERATEUR DENSITE D'ENERGIE
  6. *
  7. * CALCULE LE CHAMELEM PRODUIT SIGMA * EPS
  8. *
  9. * CHEL3 = ENER MODL CHEL1 CHEL2 ;
  10. *
  11. * MODL MODELE DE CALCUL ,TYPE MMODEL
  12. * CHEL1 CHAMP PAR ELEMENT DE SOUS TYPE CONTRAINTE OU
  13. * DEFORMATION ( TYPE MCHAML )
  14. * CHEL2 CHAMP PAR ELEMENT DE SOUS TYPE CONTRAINTE OU
  15. * DEFORMATION ( TYPE MCHAML )
  16. * CHEL3 CHAMP PAR ELEMENT DE SOUS TYPE SCALAIRE ( TYPE
  17. * MCHAML )
  18. *
  19. * PASSAGE AUX NOUVEAUX MCHAML PAR JM CAMPENON LE 01/91
  20. *_______________________________________________________________________
  21. IMPLICIT INTEGER(I-N)
  22. IMPLICIT REAL*8 (A-H,O-Z)
  23.  
  24.  
  25. -INC PPARAM
  26. -INC CCOPTIO
  27. IPMODL = 0
  28. IPCHE1 = 0
  29. IPCHE2 = 0
  30. IPCHR = 0
  31.  
  32. *
  33. * LECTURE DU MODELE
  34. *
  35. CALL LIROBJ('MMODEL ',IPMODL,1,IRETOU)
  36. CALL ACTOBJ('MMODEL ',IPMODL,1)
  37. IF (IERR.NE.0) RETURN
  38. *
  39. * LECTURE DU PREMIER MCHAML
  40. *
  41. CALL LIROBJ('MCHAML ',IPIN,1,IRETOU)
  42. CALL ACTOBJ('MCHAML ',IPIN,1)
  43. IF (IERR.NE.0) RETURN
  44. CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  45. IF(IR .NE. 1) CALL ERREUR(KER)
  46. IF(IERR .NE. 0) RETURN
  47. *
  48. * LECTURE DU DEUXIEME MCHAML
  49. *
  50. CALL LIROBJ('MCHAML ',IPIN,1,IRETOU)
  51. CALL ACTOBJ('MCHAML ',IPIN,1)
  52. IF (IERR.NE.0) RETURN
  53. CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
  54. IF(IR .NE. 1) CALL ERREUR(KER)
  55. IF(IERR .NE. 0) RETURN
  56. *
  57. CALL ENERCA(IPMODL,IPCHE1,IPCHE2,IPCHR)
  58. IF(IERR .NE. 0) RETURN
  59. *
  60. CALL ACTOBJ('MCHAML',IPCHR,1)
  61. CALL ECROBJ('MCHAML',IPCHR)
  62.  
  63. END
  64.  
  65.  
  66.  

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