Télécharger elasti.eso

Retour à la liste

Numérotation des lignes :

  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