Télécharger cdrago.eso

Retour à la liste

Numérotation des lignes :

  1. C CDRAGO SOURCE BP208322 17/03/01 21:15:25 9325
  2. SUBROUTINE CDRAGO(WRK52,WRK53,WRK54)
  3. C DRAGON SOURCE AM 00/12/13 21:32:32 4045
  4. c SUBROUTINE DRAGON(WRK0,WRK1,WRK5,DRA0,KERRE)
  5. C=======================================================================
  6. C
  7. C MODELE ELASTIQUE ENDOMMAGEABLE A ANISOTROPIE INDUITE
  8. C CHARGEMENT PROPORTIONNEL EN DEFORMATIONS
  9. C
  10. C - PROGRAMME PRINCIPAL -
  11. C VERSION 1.0
  12. C
  13. C=======================================================================
  14. C
  15. C CREATION :F.CORMERY
  16. C E.N.S.M.A - LMPM
  17. C JUIN 1993
  18. C
  19. C MODIFICATIONS: T.DESOYER, D.HALM
  20. C LMPM-ENSMA
  21. C AVRIL 1996
  22. C
  23. C=======================================================================
  24. IMPLICIT INTEGER(I-N)
  25. IMPLICIT REAL*8 (A-H,O-Z)
  26. -INC DECHE
  27. REAL*8 LAMB,MUT
  28. C-----------------------------------------------------------------------
  29. C-------- COMMONS SPECIFIQUES
  30. C-----------------------------------------------------------------------
  31. COMMON/APENTE/ LAMB,MUT,G,ALP1,BETA
  32. COMMON/MAED/B,C0,CPRIME
  33. COMMON/DSIGEP/ H
  34. C
  35. SEGMENT DRA0
  36. REAL*8 AAA(LDRA0)
  37. ENDSEGMENT
  38. *
  39. DIMENSION KM(6)
  40. *
  41. C-----------------------------------------------------------------------
  42. C------- DATA ET DIMENSION
  43. C-----------------------------------------------------------------------
  44. C N9 N18
  45. * DIMENSION SIGT(6),EPST(6),ENDT(6),VAL(3,3),VAL1(3,3),
  46. CC DIMENSION VAL(3,3),VAL1(3,3),
  47. C N24 N30 N36 N39 N42
  48. CC * SIGD(6),EPSD(6),ENDD(6),VP1(3),VP(3),
  49. C N78 N114 N120
  50. * * DEPS(6),VDEP(36),VDEP1(36),EPSPLS(6),
  51. CC * VDEP(36),VDEP1(36),EPSPLS(6),
  52. C N129 N138 N141 N144 N150
  53. CC * VAL2(3,3),VAL3(3,3),VP2(3),VP3(3),EPSTT(6)
  54. C
  55. C DATA ZERO/0.D0/
  56. C-----------------------------------------------------------------------
  57. C CHARGEMENT DES DONNEES DU PROBLEME (LECTURE SUR *.DAT)
  58. C-----------------------------------------------------------------------
  59. KERRE=0
  60. LDRA0=951
  61. SEGINI DRA0
  62. LAMB = XMAT(1) * XMAT(2) /(( 1.D0 + XMAT(2))*(1.D0-2.d0*XMAT(2)))
  63. MUT = XMAT(1) / ( 2.d0 *( 1.D0+XMAT(2)))
  64. G=XMAT(5)
  65. ALP1=XMAT(6)
  66. BETA=XMAT(7)
  67. C0=XMAT(8)
  68. CPRIME=XMAT(9)
  69. B=XMAT(10)
  70. * READ(97,303)LAMB,MUT,G,ALP1,BETA,B,C0,CPRIME
  71. * DO 5 I=1,6
  72.  
  73. *5 READ(97,304) SIGT(I)
  74. * DO 6 I=1,6
  75. *6 READ(97,304) ENDT(I)
  76. * DO 7 I=1,6
  77. *7 READ(97,304) EPST(I)
  78. * DO 8 I=1,6
  79. *8 READ(97,304) DEPS(I)
  80. * READ(97,305) MLL
  81. NPA=0
  82. DD= VAR0(1)
  83. XLIMT=C0+CPRIME*DD
  84. XLIMD=C0+CPRIME*DD
  85.  
  86. C-----------------------------------------------------------------------
  87. * CALL ARED3M(SIG0,EPIN0,VAR0(2),XLIMT,
  88. * * SIGF,EPINF,VARF(2),XLIMD,
  89. * * DEPST,VDEP,VDEP1,NPA,VAL,VP,EPSPLS,VAL1,VP1,
  90. * * VAL2,VAL3,VP2,VP3,EPSTT,FDDD)
  91. CALL ARED3M(SIG0,EPIN0,VAR0(2),XLIMT,
  92. * SIGF,EPINF,VARF(2),XLIMD,
  93. * DEPST,AAA(25),AAA(61),NPA,
  94. * AAA(1),AAA(22),AAA(97),AAA(10),AAA(19),
  95. * AAA(103),AAA(112),AAA(121),AAA(124),
  96. * AAA(127),FDDD,AAA(133),AAA(142),AAA(148),
  97. * AAA(184),AAA(220),AAA(226),AAA(307),
  98. * AAA(388),AAA(394),AAA(400),KM,AAA)
  99. * CALL ARED3M(SIGT,EPST,ENDT,XLIMT,
  100. * * SIGD,EPSD,ENDD,XLIMD,
  101. * * DEPS,VDEP,VDEP1,NPA,VAL,VP,EPSPLS,VAL1,VP1,
  102. * * VAL2,VAL3,VP2,VP3,EPSTT,FDDD)
  103. VARF(1)=VARF(2)+VARF(3)+VARF(4)
  104. C-----------------------------------------------------------------------
  105. C CALL MATGE(VDEP,TRAD,NPA,IBIF0,IBIF1,IBIF2,IBIF3)
  106. C-----------------------------------------------------------------------
  107.  
  108. DO 11 II=1,6
  109. DEFP(II)=DEPST(II)
  110. * SIGT(II)=SIGD(II)
  111. * EPST(II)=EPSD(II)
  112. * ENDT(II)=ENDD(II)
  113. 11 CONTINUE
  114. * XLIMT=XLIMD
  115. C-----------------------------------------------------------------------
  116. segsup DRA0
  117. RETURN
  118. END
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  

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