Télécharger energi.eso

Retour à la liste

Numérotation des lignes :

  1. C ENERGI SOURCE CB215821 18/09/21 21:15:41 9930
  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. CALL ACTOBJ('MMODEL ',IPMODL,1)
  36. *
  37. * LECTURE DU PREMIER MCHAML
  38. *
  39. CALL LIROBJ('MCHAML',IPIN,1,IRETOU)
  40. IF (IERR.NE.0) RETURN
  41. CALL ACTOBJ('MCHAML ',IPIN,1)
  42. CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  43. IF(IR .NE. 1) CALL ERREUR(KER)
  44. IF(IERR .NE. 0) RETURN
  45. *
  46. * LECTURE DU DEUXIEME MCHAML
  47. *
  48. CALL LIROBJ('MCHAML',IPIN,1,IRETOU)
  49. IF (IERR.NE.0) RETURN
  50. CALL ACTOBJ('MCHAML ',IPIN,1)
  51. CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
  52. IF(IR .NE. 1) CALL ERREUR(KER)
  53. IF(IERR .NE. 0) RETURN
  54. *
  55. CALL ENERCA(IPMODL,IPCHE1,IPCHE2,IPCHR)
  56. IF(IERR .NE. 0) RETURN
  57. *
  58. CALL ECROBJ('MCHAML',IPCHR)
  59. *
  60. RETURN
  61. END
  62.  
  63.  
  64.  
  65.  

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