Télécharger la plaquette Cast3M


fr en

Télécharger dragon.eso

Retour à la liste

Numérotation des lignes : 

dragon

  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