Télécharger sigma.eso

Retour à la liste

Numérotation des lignes :

  1. C SIGMA SOURCE PV 19/01/25 21:15:10 10085
  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(2)
  42. C
  43. DATA MOEPS2/'LINE','NOER'/
  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. iret1=0
  58. iret2=0
  59. ireps2=0
  60. inoer=0
  61. CALL LIRMOT(MOEPS2,2,IRET1,0)
  62. if (iret1.ne.0) call lirmot(moeps2,2,iret2,0)
  63. if (iret1.eq.1.or.iret2.eq.1) ireps2=1
  64. if (iret1.eq.2.or.iret2.eq.2) inoer=1
  65. C
  66. C ON LIT UN CHAMP/POINT
  67. C
  68. CALL LIROBJ('CHPOINT ',ICHP1,1,IRT1)
  69. IF(IERR.NE.0) RETURN
  70. C
  71. C LECTURE DU MODELE
  72. C
  73. CALL LIROBJ('MMODEL',IPMODL,1,IRT1)
  74. IF(IERR.NE.0) RETURN
  75. CALL ACTOBJ('MMODEL ',IPMODL,1)
  76. C
  77. C LECTURE DU 1ER MCHAML
  78. C
  79. CALL LIROBJ('MCHAML',IPCHA1,1,IRT1)
  80. IF(IERR.NE.0) RETURN
  81. * reduction du mchaml
  82. CALL ACTOBJ('MCHAML ',IPCHA1,1)
  83. call reduaf(ipcha1,ipmodl,ipcha10,0,ier,ker)
  84. if (ier.ne.1) call erreur(ker)
  85. ipcha1=ipcha10
  86. C
  87. C LECTURE DU 2EME MCHAML
  88. C
  89. CALL LIROBJ('MCHAML',IPCHA2,0,IRT1)
  90. if (irt1.ne.0) then
  91. * reduction du mchaml
  92. CALL ACTOBJ('MCHAML ',IPCHA2,1)
  93. call reduaf(ipcha2,ipmodl,ipcha20,0,ier,ker)
  94. if (ier.ne.1) call erreur(ker)
  95. ipcha2=ipcha20
  96. endif
  97. if (ierr.ne.0) return
  98. C
  99. CALL RNGCHA(IPCHA1,IPCHA2,'CARACTERISTIQUES',
  100. 1 'MATRICE DE HOOKE',IPCHE1,IPCHE2)
  101. IF(IERR.NE.0) RETURN
  102. IF (IPCHE2.EQ.0) THEN
  103. IMAT=1
  104. ELSE
  105. IMAT=2
  106. ENDIF
  107. C
  108. C CALCUL DES CONTRAINTES
  109. C
  110. CALL SIGMAP(IREPS2,IPMODL,ICHP1,IPCHE1,IPCHE2,IMAT,IPSTRS,IRET,
  111. > inoer)
  112. C
  113. C ECRITURE DU RESULTAT
  114. C
  115. IF(IRET.EQ.1) THEN
  116. CALL ECROBJ('MCHAML',IPSTRS)
  117. ENDIF
  118.  
  119. RETURN
  120. END
  121.  
  122.  
  123.  
  124.  
  125.  

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