Télécharger dragon.eso

Retour à la liste

Numérotation des lignes :

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

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