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    JK148537  25/12/12    21:15:04     12418          
  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.       MODORI = IPMODL
  41.       CALL PIMODL(MODORI,IPMODL,MAILDG,1)
  42.       if (ierr.ne.0) return
  43.       IF (IPMODL.EQ.0) then
  44.         call erreur(21)
  45.         return
  46.       ENDIF
  47. *
  48. *     ON LIT LE 1-ER CHAMP/ELEMENT
  49. *
  50.       CALL LIROBJ('MCHAML  ',IPIN,1,IRETOU)
  51.       IF(IERR.NE.0) RETURN
  52.       CALL ACTOBJ('MCHAML  ',IPIN,1)
  53.       CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  54.       IF(IR   .NE. 1) CALL ERREUR(KER)
  55.       IF(IERR .NE. 0) RETURN
  56. *
  57. *     LECTURE DU MCHAML DE MATERIAU
  58. *
  59.       CALL LIROBJ('MCHAML  ',IPIN,1,IRETOU)
  60.       IF(IERR.NE.0) RETURN
  61.       CALL ACTOBJ('MCHAML  ',IPIN,1)
  62.       CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
  63.       IF(IR   .NE. 1) CALL ERREUR(KER)
  64.       IF(IERR .NE. 0) RETURN
  65. C
  66. C       LECTURE DU MCHAML DE VARIABLES INTERNES
  67. C
  68.       CALL LIROBJ('MCHAML  ',IPIN,0,IRETOU)
  69.       IF (IERR.NE.0) RETURN
  70.       IPCHE3 = 0
  71.       IF (IRETOU .EQ. 1) THEN
  72.         CALL ACTOBJ('MCHAML  ',IPIN,1)
  73.         CALL REDUAF(IPIN,IPMODL,IPCHE3,0,IR,KER)
  74.         IF(IR   .NE. 1) CALL ERREUR(KER)
  75.         IF(IERR .NE. 0) RETURN
  76.       ENDIF
  77. *
  78. *     CALCUL DES CONTRAINTES OU DES DEFORMATIONS
  79. *
  80.       SEGACT,MCOORD
  81.       CALL ELAS1 (IPMODL,IPCHE1,IPCHE2,IPCHE3, IPSTRS,IRET)
  82.       SEGDES,MCOORD
  83.       IF (IERR .NE. 0) RETURN
  84. *
  85.       IF (IRET.EQ.1) THEN
  86.         CALL ACTOBJ('MCHAML  ',IPSTRS,1)
  87.         CALL ECROBJ('MCHAML  ',IPSTRS)
  88.       ENDIF
  89.  
  90.       END
  91.  
  92.  
  93.  

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