Télécharger la plaquette Cast3M


fr en

Télécharger elasti.eso

Retour à la liste

Numérotation des lignes : 

elasti

  1. C ELASTI    SOURCE    CB215821  20/07/31    14:32:23     10671          
  2.  
  3.       SUBROUTINE ELASTI
  4.  
  5. *_______________________________________________________________________
  6. *
  7. *                     OPERATEUR ELASTICITE
  8. *
  9. *      SIG = ELAS MODL EPS MATER ( VA1 )
  10. *                OU                         EN ELASTICITE LINEAIRE
  11. *      EPS = ELAS MODL SIG MATER ( VA1 )
  12. *
  13. *            MODL   MODELE DE CALCUL,TYPE MMODEL
  14. *            EPS    CHAMP PAR ELEMENT DE DEFORMATIONS,TYPE MCHAML
  15. *            SIG    CHAMP PAR ELEMENT DE CONTRAINTES,TYPE MCHAML
  16. *            MATER  CHAMP PAR ELEMENT DE CARACTERISTIQUE,TYPE MCHAML
  17. *            VA1    CHAMP PAR ELEMENT DE VARIABLES INTERNES,TYPE MCHAML
  18. *_______________________________________________________________________
  19. *
  20.       IMPLICIT INTEGER(I-N)
  21.       IMPLICIT REAL*8(A-H,O-Z)
  22. *
  23.  
  24. -INC PPARAM
  25. -INC CCOPTIO
  26. -INC SMCOORD
  27. *
  28.       IPCHE1=0
  29.       IPCHE2=0
  30.       IPCHE3=0
  31.       IPSTRS=0
  32.       IRETOU=0
  33.       IRET  =0
  34. *
  35. *     LECTURE DU MODELE
  36. *
  37.       CALL LIROBJ('MMODEL  ',IPMODL,1,IRETOU)
  38.       IF(IERR.NE.0) RETURN
  39.       CALL ACTOBJ('MMODEL  ',IPMODL,1)
  40. *
  41. *     ON LIT LE 1-ER CHAMP/ELEMENT
  42. *
  43.       CALL LIROBJ('MCHAML  ',IPIN,1,IRETOU)
  44.       IF(IERR.NE.0) RETURN
  45.       CALL ACTOBJ('MCHAML  ',IPIN,1)
  46.       CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  47.       IF(IR   .NE. 1) CALL ERREUR(KER)
  48.       IF(IERR .NE. 0) RETURN
  49. *
  50. *     LECTURE DU MCHAML DE MATERIAU
  51. *
  52.       CALL LIROBJ('MCHAML  ',IPIN,1,IRETOU)
  53.       IF(IERR.NE.0) RETURN
  54.       CALL ACTOBJ('MCHAML  ',IPIN,1)
  55.       CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
  56.       IF(IR   .NE. 1) CALL ERREUR(KER)
  57.       IF(IERR .NE. 0) RETURN
  58. C
  59. C       LECTURE DU MCHAML DE VARIABLES INTERNES
  60. C
  61.       CALL LIROBJ('MCHAML  ',IPIN,0,IRETOU)
  62.       IF (IERR.NE.0) RETURN
  63.       IPCHE3 = 0
  64.       IF (IRETOU .EQ. 1) THEN
  65.         CALL ACTOBJ('MCHAML  ',IPIN,1)
  66.         CALL REDUAF(IPIN,IPMODL,IPCHE3,0,IR,KER)
  67.         IF(IR   .NE. 1) CALL ERREUR(KER)
  68.         IF(IERR .NE. 0) RETURN
  69.       ENDIF
  70. *
  71. *     CALCUL DES CONTRAINTES OU DES DEFORMATIONS
  72. *
  73.       SEGACT,MCOORD
  74.       CALL ELAS1 (IPMODL,IPCHE1,IPCHE2,IPCHE3, IPSTRS,IRET)
  75.       SEGDES,MCOORD
  76.       IF (IERR .NE. 0) RETURN
  77. *
  78.       IF (IRET.EQ.1) THEN
  79.         CALL ACTOBJ('MCHAML  ',IPSTRS,1)
  80.         CALL ECROBJ('MCHAML  ',IPSTRS)
  81.       ENDIF
  82.  
  83.       END
  84.  
  85.  

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