Télécharger ecbpel.eso

Retour à la liste

Numérotation des lignes :

  1. C ECBPEL SOURCE BP208322 17/03/01 21:17:05 9325
  2. SUBROUTINE ECBPEL(SIG0,NSTRS,DEPST,VAR0,XMAT,NCOMAT,XCAR,ICARA,
  3. 1NVARI,SIGF,VARF,DEFP,MFR1,DDAUX,CMATE,VALMAT,VALCAR,N2EL,N2PTEL,
  4. 2 NBPGAU,IFOU,IB,IGAU,EPAIST,MELE,NPINT,NBGMAT,NELMAT,SECT,
  5. 3 LHOOK,TXR,XLOC,XGLOB,D1HOOK,ROTHOO,DDHOMU,CRIGI,DSIGT,KERRE,DT)
  6. C
  7. C ==================================================================
  8. C CE SOUS-PROGRAMME EST APPELE DANS "ECOUL2" VIA ECOU25.
  9. C IL PREPARE L'INTEGRATION DE LA FORMULE DE RELAXATION
  10. C DE CABLES DE PRECONTRAINTE SELON LA METHODE DU TEMPS
  11. C EQUIVALENT PRECONISEE PAR LE BPEL 91
  12. C
  13. C ENTREES:
  14. C -------
  15. C NSTRS = NBR. DE COMPOSANTES DES CONTR. OU DES DEFORM.
  16. C SIG0(NSTRS) = CONTR. AU DEBUT DU PAS D'INTEGRATION
  17. C DEPST(NSTRS) = INCREMENT DES DEFORM. CALCULE ELASTIQUEMENT
  18. C A PARTIR DE L'INCREMENT DES DEFORM. TOTALES
  19. C NVARI = NBR. DE VARIABLES INTERNES
  20. C VAR0(NVARI) = VARIABLES INTERNES AU DEBUT DU PAS D'INTEGRATION
  21. C CE TABLEAU NE CONTIENT QU'UN INDICE CORRESPONDANT A
  22. C - LA DEFORM. INELASTIQUE CUMULEE
  23. C NCOMAT = NBR. DE CARACTERISTIQUES MECANIQUES DU MATERIAU
  24. C XMAT(NCOMAT) = CARACTERISTIQUES MECANIQUES DU MATERIAU
  25. C MFR = INDICE DE LA FORMULATION MECANIQUE
  26. C ICARA = NBR. DE CARACT. GEOMETRIQUES DES ELEMENTS FINIS
  27. C XCAR(ICARA) = CARACT. GEOMETRIQUES DES ELEMENTS FINIS
  28. * DDAUX = MATRICE DE HOOKE ELASTIQUE
  29. * NSTRS = NBRE DE COMPOSANTES DES DEFORMATIONS
  30. * CMATE = NOM DU MATERIAU
  31. * VALMAT= TABLEAU DE CARACTERISTIQUES DU MATERIAU
  32. * VALCAR= TABLEAU DE CARACTERISTIQUES GEOMETRIQUES
  33. * N2EL = NBRE D ELEMENTS DANS SEGMENT DE HOOKE
  34. * N2PTEL= NBRE DE POINTS DANS SEGMENT DE HOOKE
  35. * IFOU = OPTION DE CALCUL
  36. * IB = NUMERO DE L ELEMENT COURANT
  37. * IGAU = NUMERO DU POINT COURANT
  38. * EPAIST= EPAISSEUR
  39. * NBPGAU= NBRE DE POINTS DE GAUSS
  40. * MELE = NUMERO DE L ELEMENT FINI
  41. * NPINT = NBRE DE POINTS D INTEGRATION
  42. * NBGMAT= NBRE DE POINTS DANS SEGMENT DE CARACTERISTIQUES
  43. * NELMAT= NBRE D ELEMENTS DANS SEGMENT DE CARACTERISTIQUES
  44. * SECT = SECTION
  45. * LHOOK = TAILLE DE LA MATRICE DE HOOKE
  46. * TXR,XLOC,XGLOB,D1HOOK,ROTHOO,DDHOMU,CRIGI = TABLEAUX UTILISES
  47. * POUR LE CALCUL DE LA MATRICE DE HOOKE
  48. C
  49. C SORTIES:
  50. C -------
  51. C SIGF(NSTRS)= CONTR. A LA FIN DU PAS D'INTEGRATION
  52. C VARF(NVARI)= VARIABLES INTERNES A LA FIN DU PAS D'INTEGRATION
  53. C DEFP(NSTRS)= INCREMENT DES DEFORM. PLASTIQUES A LA FIN DU PAS
  54. C D'INTEGRATION
  55. C KERRE = INDICE QUI REGIT LES ERREURS
  56. C = 77 SI LA DEFORM. PLAST. CUMULEE ENDOMMAGEE (2IEME VAR.
  57. C INT.) EST EN DEHORS DE LA COURBE DE TRACTION, DS.
  58. C LE CAS DE L'ECROUISSAGE ET DE L'ENDOMM. ISOTROPES.
  59. C CECI PEUT SE PRODUIRE SUITE A L'APPEL A "CRIDAM"
  60. C = 99 SI LA FORMULATION MECANIQUE N'EST PAS DISPONIBLE
  61. C POUR LE MODELE CONSIDERE OU S'IL Y A INCOMPATIBILITE
  62. C ENTRE MFR ET IFOUR
  63. C ==================================================================
  64. C ICI IL FAUT PROGRAMMER EN FORTRAN PUR
  65. C ===================================================================
  66. C
  67. IMPLICIT INTEGER(I-N)
  68. IMPLICIT REAL*8(A-H,O-Z)
  69. -INC CCOPTIO
  70. DIMENSION SIG0(*),DEPST(*),VAR0(*),XMAT(*),XCAR(*),SIGF(*),
  71. & VARF(*),DEFP(*)
  72. DIMENSION RSIG0(6),RDSIGT(6),RDEFP(6),RSIGF(6),DSIGT(*)
  73. DIMENSION VALMAT(*),VALCAR(*)
  74. DIMENSION TXR(IDIM,*),CRIGI(12)
  75. DIMENSION DDAUX(LHOOK,*),DDHOMU(LHOOK,*)
  76. DIMENSION XLOC(3,3),XGLOB(3,3)
  77. DIMENSION D1HOOK(LHOOK,*),ROTHOO(LHOOK,*)
  78. CHARACTER*8 CMATE
  79. REAL*8 K1,K2,MU
  80. C
  81. C CALCUL DU TEMPS EQUIVALENT
  82. C
  83. EACIER = XMAT(1)
  84. AIRE = XCAR(1)
  85. TENS0 = SIG0(1)
  86. SIGP0 = (TENS0/AIRE) + (EACIER*DEPST(1))
  87. EPIN0 = VAR0(1)
  88. DSIGP0 = EPIN0*EACIER
  89. RO1000 = XMAT(5)
  90. FPRG = XMAT(6)
  91. TCHA = XMAT(7)
  92.  
  93. IF (RO1000 .LE. 2.5D0) THEN
  94. K1 = 6.D-3
  95. K2 = 1.1D0
  96. ELSE
  97. K1 = 8.D-3
  98. K2 = 1.25D0
  99. IF (RO1000 .GE. 8.D0) THEN
  100. K2 = 1.5D0
  101. ENDIF
  102. ENDIF
  103.  
  104. * WRITE(6,*) 'DT = ',DT
  105. * WRITE(6,*) 'DEPST(1) = ',DEPST(1)
  106. * WRITE(6,*) 'SIG0(1) = ',SIG0(1)
  107. * WRITE(6,*) 'VAR0(1) = ',VAR0(1)
  108. * WRITE(6,*) 'SIGP0 = ',SIGP0
  109. * WRITE(6,*) 'DSIGP0 = ',DSIGP0
  110.  
  111. ** WRITE(6,*) 'EACIER = ',EACIER
  112. ** WRITE(6,*) 'AIRE = ',AIRE
  113. ** WRITE(6,*) 'TENS0 = ',TENS0
  114. ** WRITE(6,*) 'SIGP0 = ',SIGP0
  115. ** WRITE(6,*) 'EPIN0 = ',EPIN0
  116. ** WRITE(6,*) 'DSIGP0 = ',DSIGP0
  117. ** WRITE(6,*) 'RO1000 = ',RO1000
  118. ** WRITE(6,*) 'FPRG = ',FPRG
  119. ** WRITE(6,*) 'TCHA = ',TCHA
  120. ** WRITE(6,*) 'K1 = ',K1
  121. ** WRITE(6,*) 'K2 = ',K2
  122.  
  123. CALL TEQUIV(K1,RO1000,TCHA,FPRG,K2,SIGP0,DSIGP0,TEQ)
  124.  
  125. * WRITE(6,*) 'TEQ = ',TEQ
  126. MU = (SIGP0+DSIGP0)/FPRG
  127.  
  128. * WRITE(6,*) 'MU = ',MU
  129. * WRITE(6,*) 'DT = ',DT
  130.  
  131. DEPIN0 = (K1*RO1000*((24.D0*(TEQ+DT-TCHA)/1.D3)
  132. & **(3.D0*(1.D0-MU)/4.D0))*
  133. & (EXP((10.D0*MU-7.5D0)/K2))*(SIGP0+DSIGP0) - DSIGP0)
  134. & /EACIER
  135.  
  136. DELAS = DEPST(1) - DEPIN0
  137. DEFP(1) = DEPIN0
  138. VARF(1) = VAR0(1) + DEFP(1)
  139. * SIGF(1) = (AIRE*SIGP0) + (EACIER*AIRE*DEFP(1))
  140. * SIGF(1) = AIRE*SIGP0
  141. SIGF(1) = SIG0(1) + AIRE*EACIER*DELAS
  142. SIGPF = SIGF(1)/AIRE
  143.  
  144. * WRITE(6,*) 'DEFP(1) = ',DEFP(1)
  145. * WRITE(6,*) 'VARF(1) = ',VARF(1)
  146. * WRITE(6,*) 'SIGF(1) = ',SIGF(1)
  147. * WRITE(6,*) 'SIGPF = ',SIGPF
  148.  
  149. RETURN
  150. END
  151.  
  152.  
  153.  
  154.  

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