Télécharger invari.eso

Retour à la liste

Numérotation des lignes :

  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