Télécharger la plaquette Cast3M


fr en

Télécharger energi.eso

Retour à la liste

Numérotation des lignes : 

energi

  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