Télécharger elocrak1.eso

Retour à la liste

Numérotation des lignes :

  1. C ELOCRAK1 SOURCE BP208322 17/03/01 21:17:23 9325
  2. SUBROUTINE ELOCRAK1(WRK52,WRK53,WRK54,NVARI,Iecou)
  3. C
  4. C====&===1=========2=========3=========4=========5=========6=========7==
  5. C Commentaires : Subroutine permettant de mettre en oeuvre le
  6. C modele ELOCRACK (CEA) pour representer
  7. C le comportement cyclique du béton sous
  8. C sollicitations sismiques
  9. C
  10. C Traits : -
  11. C -
  12. C -
  13. C -
  14. C -
  15. C
  16. C Auteur : B. Richard (Dr - Ing.) - CEA/DEN/DANS/DM2S/SEMT/EMSI
  17. C====&===1=========2=========3=========4=========5=========6=========7==
  18. C
  19. C----DECLARATION GENERALES----------------------------------------------
  20. C
  21. IMPLICIT INTEGER(I-N)
  22. IMPLICIT REAL*8(A-H,O-Z)
  23. C
  24. C----LES INCLUDES-------------------------------------------------------
  25. C
  26. -INC CCOPTIO
  27. -INC DECHE
  28. C
  29. C----DECLARATION PARAMETRES---------------------------------------------
  30. C
  31. SEGMENT IECOU
  32. INTEGER icow1,icow2,icow3,icow4,icow5,icow6,icow7,
  33. 1 icow8,icow9,icow10,icow11,icow12,icow13,icow14,icow15,icow16,
  34. 2 icow17,icow18,icow19,icow20,icow21,icow22,icow23,icow24,
  35. 3 icow25,icow26,icow27,icow28,icow29,icow30,ICARA,
  36. 4 icow32,icow33,NSTRS1,MFR1,icow36,icow37,icow38,
  37. 5 icow39,icow40,icow41,icow42,icow43,icow44
  38.  
  39. INTEGER icow45,icow46,icow47,icow48,icow49,icow50,
  40. 1 icow51,icow52,icow53,icow54,icow55,icow56
  41. 2 icow57,icow58
  42. ENDSEGMENT
  43.  
  44. SEGMENT WRKK2
  45. REAL*8 EPSILI(NSTRSV)
  46. END SEGMENT
  47.  
  48. SEGMENT WRK6
  49. REAL*8 SIG0S(NSTRS),DEPSTS(NSTRS)
  50. END SEGMENT
  51.  
  52. INTEGER NVARI,KCAS,IRTD,ISTRS
  53.  
  54. REAL*8 PREC,EPAI,FAC,AUX,AUX1,AUX2,YOUN,XNU,UN
  55.  
  56. PARAMETER (UN=1.D0)
  57.  
  58. KERRE=0
  59. C
  60. C----MISE EN DONNEES----------------------------------------------------
  61. C
  62. YOUN = XMAT(1)
  63. XNU = XMAT(2)
  64. C
  65. C----CAS DE LA FORMULATION COQUES MINCES--------------------------------
  66. C
  67. IF (MFR1 .EQ. 9) THEN
  68.  
  69. EPAI=XCARB(1)
  70. SEGINI WRK6
  71.  
  72. DO 101 ISTRS=1,NSTRS1
  73. SIG0S(ISTRS)=SIG0(ISTRS)
  74. DEPSTS(ISTRS)=DEPST(ISTRS)
  75. 101 CONTINUE
  76.  
  77. NSTRSV=4
  78. IFOUR2=-2
  79. JDIM =2
  80.  
  81. DO 102 ISTRS=1,2
  82. SIG0(ISTRS)=SIG0S(ISTRS)/EPAI
  83. 102 CONTINUE
  84.  
  85. DEPST(3)=0.D0
  86. DEPST(4)=DEPSTS(3)
  87.  
  88. SIG0(3)=0.D0
  89. SIG0(4)=SIG0S(3)/EPAI
  90. C
  91. C----CAS DE LA FORMULATION MASSIVE--------------------------------------
  92. C
  93. ELSE IF (MFR1 .EQ. 1) THEN
  94.  
  95. NSTRSV=NSTRS1
  96. IFOUR2=IFOUR
  97. JDIM = IDIM
  98.  
  99. ELSE
  100.  
  101. PRINT*,'MFR1=',MFR1
  102. KERRE=57
  103. RETURN
  104.  
  105. END IF
  106.  
  107. SEGINI WRKK2
  108.  
  109. CMATE = 'ISOTROPE'
  110. KCAS=1
  111.  
  112. CALL DOHMAS(XMAT,CMATE,IFOUR2,NSTRSV,KCAS,DDHOOK,IRTD)
  113.  
  114. IF ( IRTD .EQ. 1) THEN
  115.  
  116. CALL MATVE1 (DDHOOK,DEPST,NSTRSV,NSTRSV,DSIGT,1)
  117.  
  118. PREC=1.D-08
  119.  
  120. CALL DOHMAS(XMAT,CMATE,IFOUR2,NSTRSV,2,DDHOOK,IRTD)
  121.  
  122. CALL INVALM(DDHOOK,NSTRSV,NSTRSV,IRTD,PREC)
  123.  
  124. IF (IRTD.EQ.0)THEN
  125.  
  126. CALL MATVE1 (DDHOOK,SIG0,NSTRSV,NSTRSV,EPSILI,1)
  127.  
  128. icarbi=icara
  129. C
  130. C----CALCUL DES CONTRAINTES DE CAUCHY DANS LE CAS 3D--------------------
  131. C
  132. CALL ELOCRAK2(WRK52,WRK53,WRK54,WRKK2,NSTRSV,NVARI,
  133. 1 ICARbi,JDIM,IFOUR2)
  134.  
  135. icara=icarbi
  136. C
  137. C----CALCUL DES CONTRAINTES VRAIES DANS LE CAS DES COQUES MINCES--------
  138. C
  139. IF (MFR1 .EQ. 9) THEN
  140.  
  141. FAC=(EPAI**3)/12.D0
  142. AUX =FAC*YOUN/(1.D0-XNU*XNU)
  143. AUX1=FAC*YOUN*.5D0/(1.D0+XNU)
  144. AUX2=EPAI*YOUN*.5D0/(1.d0+XNU)/1.2d0
  145.  
  146. DO 103 ISTRS=1,NSTRS1
  147. SIG0(ISTRS)=SIG0S(ISTRS)
  148. DEPST(ISTRS)=DEPSTS(ISTRS)
  149. 103 CONTINUE
  150.  
  151. DO 104 ISTRS=1,2
  152. SIGF(ISTRS)=SIGF(ISTRS)*EPAI
  153. 104 CONTINUE
  154.  
  155. SIGF(3)=SIGF(4)*EPAI
  156. SIGF(4)=SIG0(4)+AUX*(DEPST(4)+XNU*DEPST(5))
  157. SIGF(5)=SIG0(5)+AUX*(DEPST(5)+XNU*DEPST(4))
  158. SIGF(6)=SIG0(6)+AUX1*DEPST(6)
  159. SIGF(7)=SIG0(7)+AUX2*DEPST(7)
  160. SIGF(8)=SIG0(8)+AUX2*DEPST(8)
  161. SEGSUP WRK6
  162. END IF
  163.  
  164. ELSE
  165. PRINT*,'ERREUR DANS INVALM'
  166. KERRE=56
  167. END IF
  168. ELSE
  169. PRINT*,'ERREUR DANS DOHMAS'
  170. KERRE=56
  171. END IF
  172. SEGSUP WRKK2
  173. RETURN
  174. END
  175.  
  176.  
  177.  
  178.  
  179.  
  180.  
  181.  
  182.  
  183.  
  184.  
  185.  
  186.  
  187.  
  188.  

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