Télécharger elasti.eso

Retour à la liste

Numérotation des lignes :

  1. C ELASTI SOURCE CB215821 16/12/05 21:39:21 9237
  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. -INC CCOPTIO
  24. *
  25. IPCHE1=0
  26. IPCHE2=0
  27. IPCHE3=0
  28. IPSTRS=0
  29. IRETOU=0
  30. IRET =0
  31. *
  32. * LECTURE DU MODELE
  33. *
  34. CALL LIROBJ('MMODEL',IPMODL,1,IRETOU)
  35. IF(IERR.NE.0) RETURN
  36. *
  37. * ON LIT LE 1-ER CHAMP/ELEMENT
  38. *
  39. CALL LIROBJ('MCHAML',IPIN,1,IRETOU)
  40. IF(IERR.NE.0) RETURN
  41. CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  42. IF(IR .NE. 1) CALL ERREUR(KER)
  43. IF(IERR .NE. 0) RETURN
  44. *
  45. * LECTURE DU MCHAML DE MATERIAU
  46. *
  47. CALL LIROBJ('MCHAML',IPIN,1,IRETOU)
  48. IF(IERR.NE.0) RETURN
  49. CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
  50. IF(IR .NE. 1) CALL ERREUR(KER)
  51. IF(IERR .NE. 0) RETURN
  52. C
  53. C LECTURE DU MCHAML DE VARIABLES INTERNES
  54. C
  55. CALL LIROBJ('MCHAML',IPIN,0,IRETOU)
  56. IF (IERR.NE.0) RETURN
  57. IPCHE3 = 0
  58. IF (IRETOU .EQ. 1) THEN
  59. CALL REDUAF(IPIN,IPMODL,IPCHE3,0,IR,KER)
  60. IF(IR .NE. 1) CALL ERREUR(KER)
  61. IF(IERR .NE. 0) RETURN
  62. ENDIF
  63. *
  64. * CALCUL DES CONTRAINTES OU DES DEFORMATIONS
  65. *
  66. CALL ELAS1 (IPMODL,IPCHE1,IPCHE2,IPCHE3, IPSTRS,IRET)
  67. IF (IERR .NE. 0) RETURN
  68. *
  69. IF (IRET.EQ.1) CALL ECROBJ('MCHAML',IPSTRS)
  70.  
  71. RETURN
  72. END
  73.  
  74.  
  75.  

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