Télécharger sigma.eso

Retour à la liste

Numérotation des lignes :

  1. C SIGMA SOURCE FANDEUR 16/01/07 21:16:04 8756
  2. SUBROUTINE SIGMA
  3. *_______________________________________________________________________
  4. *
  5. * OPERATEUR CONTRAINTES LINEAIRES
  6. *
  7. * Ancienne syntaxe :
  8. * __________________
  9. *
  10. * SIG1=SIGMA ( II ) MAT1 ( CAR1 ) CHP1 ;
  11. *
  12. * II mot cle indiquant que l'on veut EPS2
  13. * MAT1 champ par {l{ment de mat{riau ou de HOOKE
  14. * CAR1 champ par {l{ment de caract{ristiques (facultatif)
  15. * CHP1 CHPOINT de d{placement
  16. *
  17. * JUILLET 84 EBERSOLT
  18. * FEVRIER 87 BROCHARD :
  19. * - LE CALCUL EFFECTIF EST FAIT DANS SIGMAP
  20. *
  21. *
  22. * Nouvelle syntaxe :
  23. * __________________
  24. *
  25. * SIG1=SIGMA ( II ) MODL1 | CAR1 CHP1 ;
  26. * | HOO1 (CAR1)
  27. *
  28. * II mot cle indiquant que l'on veut EPS2
  29. * MODL1 mod}le de calcul, type MMODEL
  30. * CAR1 champ par {l{ment de caract{ristiques, type MCHAML
  31. * HOO1 champ par {l{ment de matrice de HOOKE, type MCHAML
  32. * CHP1 CHPOINT de d{placement
  33. *
  34. * Passage aux nouveaux CHAMELEMs par I.Monnier le 6 6 1990
  35. *
  36. *_______________________________________________________________________
  37. IMPLICIT INTEGER(I-N)
  38. IMPLICIT REAL*8(A-H,O-Z)
  39. -INC CCOPTIO
  40. C
  41. CHARACTER*4 MOEPS2(1)
  42. C
  43. DATA MOEPS2/'LINE'/
  44. C
  45. C IPMODL MODELE MMODEL
  46. C ICHP1 CHAMPOINT DEPLACEMENTS
  47. C IPCHA1 CHAMELEM MATERIAU OU DE HOOKE
  48. C IPCHA2 CHAMELEM CARACTERISTIQUE
  49. C IPSTRS CHAMELEM CONTRAINTES
  50. IPCHE1=0
  51. IPCHE2=0
  52. IPCHA1=0
  53. IPCHA2=0
  54. C
  55. C VEUT ON CALCULER LES GRANDES DEFORMATIONS ?
  56. C
  57. CALL LIRMOT(MOEPS2,1,IREPS2,0)
  58. C
  59. C ON LIT UN CHAMP/POINT
  60. C
  61. CALL LIROBJ('CHPOINT ',ICHP1,1,IRT1)
  62. IF(IERR.NE.0) RETURN
  63. C
  64. C LECTURE DU MODELE
  65. C
  66. CALL LIROBJ('MMODEL',IPMODL,1,IRT1)
  67. IF(IERR.NE.0) RETURN
  68. C
  69. C LECTURE DU 1ER MCHAML
  70. C
  71. CALL LIROBJ('MCHAML',IPCHA1,1,IRT1)
  72. IF(IERR.NE.0) RETURN
  73. * reduction du mchaml
  74. call reduaf(ipcha1,ipmodl,ipcha10,0,ier,ker)
  75. if (ier.ne.1) call erreur(ker)
  76. ipcha1=ipcha10
  77. C
  78. C LECTURE DU 2EME MCHAML
  79. C
  80. CALL LIROBJ('MCHAML',IPCHA2,0,IRT1)
  81. if (irt1.ne.0) then
  82. * reduction du mchaml
  83. call reduaf(ipcha2,ipmodl,ipcha20,0,ier,ker)
  84. if (ier.ne.1) call erreur(ker)
  85. ipcha2=ipcha20
  86. endif
  87. if (ierr.ne.0) return
  88. C
  89. CALL RNGCHA(IPCHA1,IPCHA2,'CARACTERISTIQUES',
  90. 1 'MATRICE DE HOOKE',IPCHE1,IPCHE2)
  91. IF(IERR.NE.0) RETURN
  92. IF (IPCHE2.EQ.0) THEN
  93. IMAT=1
  94. ELSE
  95. IMAT=2
  96. ENDIF
  97. C
  98. C CALCUL DES CONTRAINTES
  99. C
  100. CALL SIGMAP(IREPS2,IPMODL,ICHP1,IPCHE1,IPCHE2,IMAT,IPSTRS,IRET)
  101. C
  102. C ECRITURE DU RESULTAT
  103. C
  104. IF(IRET.EQ.1) THEN
  105. CALL ECROBJ('MCHAML',IPSTRS)
  106. ENDIF
  107.  
  108. RETURN
  109. END
  110.  
  111.  
  112.  

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