Télécharger la plaquette Cast3M


fr en

Télécharger cdrago.eso

Retour à la liste

Numérotation des lignes : 

cdrago

  1. C CDRAGO    SOURCE    PV        17/12/08    21:15:32     9660           
  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.  
  132.  

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