Télécharger la plaquette Cast3M


fr en

Télécharger invari.eso

Retour à la liste

Numérotation des lignes : 

invari

  1. C INVARI    SOURCE    CB215821  19/07/30    21:16:56     10273          
  2. *
  3.       SUBROUTINE INVARI
  4. *
  5. *---------------------------------------------------------------------
  6. *
  7. *         CALCUL DES 3 INVARIANTS D'UN TENSEUR D'ORDRE 2
  8. *
  9. *   SYNTAXE :
  10. *   ---------
  11. *
  12. *     CH1 CH2 CH3  = INVA  MOD CHEL (CAR) ( MOT );
  13. *
  14. *           MOD   MODELE DE CALCUL, TYPE MMODEL
  15. *           CHEL  OBJET DE TYPE MCHAML DE SOUS TYPE CONTRAINTES
  16. *                 OU DEFORMATIONS
  17. *           CAR   OBJET DE TYPE MCHAML DE SOUS TYPE CARACTERISTIQUES
  18. *           MOT   MOT QUI INDIQUE POUR LES COQUES OU ON CALCULE LES
  19. *                 CONTRAINTES
  20. *           CH1   OBJET DE TYPE MCHAML DE SOUS TYPE STRESSES
  21. *           CH2   OBJET DE TYPE MCHAML DE SOUS TYPE STRESSES
  22. *           CH3   OBJET DE TYPE MCHAML DE SOUS TYPE STRESSES
  23. *
  24. *     PASSAGE AUX NOUVEAU CHAMELEM PAR JM CAMPENON LE 12/90
  25. *
  26. *--------------------------------------------------------------------
  27. *
  28.       IMPLICIT INTEGER(I-N)
  29.       IMPLICIT REAL*8(A-H,O-Z)
  30.  
  31.  
  32. -INC PPARAM
  33. -INC CCOPTIO
  34.  
  35.       CHARACTER*4 MOCLE(3)
  36.       DATA MOCLE/'SUPE','MOYE','INFE'/
  37. *
  38.       IPCHE1=0
  39.       IPCHE2=0
  40. *
  41.       CALL LIROBJ('MMODEL  ',IPMODL,1,IRET)
  42.       IF(IERR.NE.0) RETURN
  43.       CALL ACTOBJ('MMODEL  ',IPMODL,1)
  44. *
  45.       CALL LIROBJ('MCHAML  ',IPIN,1,IRET)
  46.       IF(IERR.NE.0) RETURN
  47.       CALL ACTOBJ('MCHAML  ',IPIN,1)
  48.       CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  49.       IF(IR   .NE. 1) CALL ERREUR(KER)
  50.       IF(IERR .NE. 0) RETURN
  51. *
  52.       CALL LIROBJ('MCHAML',IPIN,0,IRET)
  53.       IF (IERR.NE.0) RETURN
  54.       IPCHE2=0
  55.       IF (IRET .EQ. 1) THEN
  56.         CALL ACTOBJ('MCHAML',IPIN,1)
  57.         CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
  58.         IF(IR   .NE. 1) CALL ERREUR(KER)
  59.         IF(IERR .NE. 0) RETURN
  60.       ENDIF
  61. *
  62. *     LECTURE D'UN MOT SI NECESSAIRE
  63. *
  64.       IMIL=0
  65.       CALL LIRMOT(MOCLE,3,IMOT,0)
  66.       IF (IMOT.EQ.1) IMIL=1
  67.       IF (IMOT.EQ.3) IMIL=-1
  68. *
  69.       CALL INVACA(IPMODL,IPCHE1,IPCHS1,IPCHS2,IPCHS3,IPCHE2,IMIL,IRET)
  70.       IF(IERR .NE. 0 .OR. IRET.EQ.0) RETURN
  71. *
  72.       CALL ACTOBJ('MCHAML',IPCHS3,1)
  73.       CALL ACTOBJ('MCHAML',IPCHS2,1)
  74.       CALL ACTOBJ('MCHAML',IPCHS1,1)
  75.       CALL ECROBJ('MCHAML',IPCHS3)
  76.       CALL ECROBJ('MCHAML',IPCHS2)
  77.       CALL ECROBJ('MCHAML',IPCHS1)
  78.  
  79.       END
  80.  
  81.  
  82.  

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