Télécharger energi.eso

Retour à la liste

Numérotation des lignes :

  1. C ENERGI SOURCE CB215821 16/12/05 21:39:21 9237
  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. -INC CCOPTIO
  25. IPMODL = 0
  26. IPCHE1 = 0
  27. IPCHE2 = 0
  28. IPCHR = 0
  29.  
  30. *
  31. * LECTURE DU MODELE
  32. *
  33. CALL LIROBJ('MMODEL',IPMODL,1,IRETOU)
  34. IF (IERR.NE.0) RETURN
  35. *
  36. * LECTURE DU PREMIER MCHAML
  37. *
  38. CALL LIROBJ('MCHAML',IPIN,1,IRETOU)
  39. IF (IERR.NE.0) RETURN
  40. CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  41. IF(IR .NE. 1) CALL ERREUR(KER)
  42. IF(IERR .NE. 0) RETURN
  43. *
  44. * LECTURE DU DEUXIEME MCHAML
  45. *
  46. CALL LIROBJ('MCHAML',IPIN,1,IRETOU)
  47. IF (IERR.NE.0) RETURN
  48. CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
  49. IF(IR .NE. 1) CALL ERREUR(KER)
  50. IF(IERR .NE. 0) RETURN
  51. *
  52. CALL ENERCA(IPMODL,IPCHE1,IPCHE2,IPCHR)
  53. IF(IERR .NE. 0) RETURN
  54. *
  55. CALL ECROBJ('MCHAML',IPCHR)
  56. *
  57. RETURN
  58. END
  59.  
  60.  
  61.  

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