Télécharger defsig.eso

Retour à la liste

Numérotation des lignes :

  1. C DEFSIG SOURCE BP208322 17/03/01 21:16:59 9325
  2. SUBROUTINE DEFSIG(MFR,NDEF,
  3. . INPLAS,IND,WRK1,WRK5,WTRAV,
  4. . IVASTF,IVARIF,IVADEP,COB,XMOB,
  5. . IB,IGAU,CMATE,MATE,MELE,KERRE)
  6. *
  7. **************************************************
  8. * ENTREES
  9. **************************************************
  10. *
  11. * MFR : formulation de l'élément
  12. * NSTRS : nombre de composantes des contraintes
  13. * NVARI : nombre des variables internes
  14. * NDEF : nombre des déformations
  15. * INPLAS : numéro du matériau inélastique
  16. * IND :
  17. * SIGF : contraintes à la fin du pas (WRK1)
  18. * SIG0 : contraintes au début du pas (WRK1)
  19. * DSIGT : incrément de contrainte au cours du pas (WTRAV)
  20. * EPINF : déformations inélastiques à la fin du pas (WRK5)
  21. * EPIN0 : déformations inélastiques au début du pas (WRK5)
  22. * EPST0 : déformations totales au début du pas (WRK5)
  23. * VARF : variables internes à la fin du pas (WRK1)
  24. * IVASTF : pointeur sur un segment mptval de contraintes
  25. * IVARIF : pointeur sur un segment mptval de variables internes
  26. * IVADEP : pointeur sur un segment mptval de deformations inelastiques
  27. * COB : porosité
  28. * IB : numéro de l'élément
  29. * IGAU : numéro du point de Gauss
  30. * DEFP : incrément de déformations inélastiques (WRK1)
  31. * CMATE : nom du matériau
  32. * MATE : numéro du matériau élastique
  33. * MELE : numéro élément fini
  34. *
  35. **************************************************
  36. * SORTIES
  37. **************************************************
  38. *
  39. * KERRE : indice d'erreur
  40. * SEGMENT COMPRENANT :
  41. * - les contraintes finales
  42. * - les variables internes finales
  43. * - les incréments de déformations inélastiques
  44. *
  45. **************************************************
  46. *
  47. IMPLICIT INTEGER(I-N)
  48. IMPLICIT REAL*8(A-H,O-Z)
  49. *
  50. -INC CCOPTIO
  51. -INC SMCHAML
  52. *
  53. SEGMENT MPTVAL
  54. INTEGER IPOS(NS) ,NSOF(NS)
  55. INTEGER IVAL(NCOSOU)
  56. CHARACTER*16 TYVAL(NCOSOU)
  57. ENDSEGMENT
  58. *
  59. SEGMENT WRK1
  60. REAL*8 DDHOOK(LHOOK,LHOOK),SIG0(NSTRS),DEPST(NSTRS)
  61. REAL*8 SIGF(NSTRS),VAR0(NVARI),VARF(NVARI)
  62. REAL*8 DEFP(NSTRS),XCAR(ICARA)
  63. ENDSEGMENT
  64. *
  65. SEGMENT WRK5
  66. REAL*8 EPIN0(NSTRS),EPINF(NSTRS),EPST0(NSTRS)
  67. ENDSEGMENT
  68. *
  69. SEGMENT WTRAV
  70. REAL*8 DDAUX(LHOOK,LHOOK),VALMAT(NUMAT)
  71. REAL*8 VALCAR(NUCAR),DSIGT(NSTRS)
  72. REAL*8 TXR(IDIM,IDIM),DDHOMU(LHOOK,LHOOK)
  73. REAL*8 XLOC(3,3),XGLOB(3,3)
  74. REAL*8 D1HOOK(LHOOK,LHOOK),ROTHOO(LHOOK,LHOOK)
  75. ENDSEGMENT
  76. *
  77. CHARACTER*8 CMATE
  78. *
  79. * rearrangement pour milieu poreux
  80. *
  81. NSTRS=SIGF(/1)
  82. NVARI=VARF(/1)
  83. KERRE=0
  84. IF(MFR.EQ.33) THEN
  85. *
  86. *ZZZZZ CAS DES JOINTS NON LINEAIRES : A VOIR !
  87. *
  88. * calcul des contraintes totales
  89. *
  90. SIGF(NSTRS) = SIG0(NSTRS)+DEPST(NSTRS)/XMOB
  91. DO 1994 IC=1,3
  92. IF(IFOUR.EQ.-2.AND.IC.EQ.3) GO TO 1994
  93. SIGF(IC) = SIGF(IC)
  94. & -COB*(EPST0(NSTRS)+DEPST(NSTRS))
  95. *
  96. * ce qui suit a ete pompe plus bas
  97. *
  98. IF(IND.EQ.1.OR.INPLAS.EQ.31.OR.INPLAS.EQ.30
  99. & .OR.INPLAS.EQ.37.OR.INPLAS.EQ.66
  100. & .OR.INPLAS.EQ.141)THEN
  101. SIGF(NSTRS) = SIGF(NSTRS)
  102. & +COB*(DEPST(IC)- (EPINF(IC)-EPIN0(IC)))
  103. ELSE
  104. SIGF(NSTRS) = SIGF(NSTRS)
  105. & +COB*(DEPST(IC)- DEFP(IC))
  106. ENDIF
  107. 1994 CONTINUE
  108. ENDIF
  109. *
  110. MPTVAL=IVASTF
  111. DO 1116 IC=1,NSTRS
  112. MELVAL=IVAL(IC)
  113. VELCHE(IGAU,IB)=SIGF(IC)
  114. 1116 continue
  115. c
  116. c et les variables internes finales
  117. c
  118. MPTVAL=IVARIF
  119. IF(MFR.EQ.7.AND.CMATE.EQ.'SECTION')THEN
  120. DO IC=1,NVARI
  121. MELVAL=IVAL(IC)
  122. IELCHE(IGAU,IB)=NINT(VARF(IC))
  123. END DO
  124. *
  125. ELSE IF(INPLAS.EQ.74) THEN
  126. *
  127. * Modele MAXWELL - composante EPSE
  128. *
  129. MELVAL=IVAL(1)
  130. VELCHE(IGAU,IB)=VARF(1)
  131. *
  132. * les autres composantes
  133. *
  134. DO IC=2,NVARI
  135. MELVAL=IVAL(IC)
  136. IELCHE(IGAU,IB)=NINT(VARF(IC))
  137. END DO
  138. *
  139. ELSE
  140. DO 1117 IC=1,NVARI
  141. MELVAL=IVAL(IC)
  142. VELCHE(IGAU,IB)=VARF(IC)
  143. 1117 continue
  144. ENDIF
  145. c
  146. c et les increments de deformations plastiques
  147. c
  148. MPTVAL=IVADEP
  149. DO 1118 IC=1,NDEF
  150. MELVAL=IVAL(IC)
  151. IF(IND.EQ.1.OR.INPLAS.EQ.31.OR.INPLAS.EQ.30
  152. & .OR.INPLAS.EQ.37.OR.INPLAS.EQ.66.OR.INPLAS.EQ.118
  153. & .OR. INPLAS.EQ.141)THEN
  154. VELCHE(IGAU,IB)=EPINF(IC)-EPIN0(IC)
  155. ELSE
  156. VELCHE(IGAU,IB)=DEFP(IC)
  157. ENDIF
  158. 1118 continue
  159. c
  160. 1000 RETURN
  161. END
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  
  172.  
  173.  

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