Esope source files

   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-----------------------------------------------------------------------
  94 : *      CALL ARED3M(SIG0,EPIN0,VAR0(2),XLIMT,
  95 : *     *            SIGF,EPINF,VARF(2),XLIMD,
  96 : *     *            DEPST,VDEP,VDEP1,NPA,VAL,VP,EPSPLS,VAL1,VP1,     
  97 : *     *            VAL2,VAL3,VP2,VP3,EPSTT,FDDD)
  97 :       CALL ARED3M(SIG0,EPIN0,VAR0(2),XLIMT,
  98 :      *            SIGF,EPINF,VARF(2),XLIMD,
  99 :      *            DEPST,AAA(25),AAA(61),NPA,
 100 :      *            AAA(1),AAA(22),AAA(97),AAA(10),AAA(19),
 101 :      *            AAA(103),AAA(112),AAA(121),AAA(124),
 102 :      *            AAA(127),FDDD,AAA(133),AAA(142),AAA(148), 
 103 :      *            AAA(184),AAA(220),AAA(226),AAA(307), 
 104 :      *            AAA(388),AAA(394),AAA(400),KM,AAA)
 104 : *      CALL ARED3M(SIGT,EPST,ENDT,XLIMT,
 105 : *     *            SIGD,EPSD,ENDD,XLIMD,
 106 : *     *            DEPS,VDEP,VDEP1,NPA,VAL,VP,EPSPLS,VAL1,VP1,
 107 : *     *            VAL2,VAL3,VP2,VP3,EPSTT,FDDD)
 108 :       VARF(1)=VARF(2)+VARF(3)+VARF(4)
 109 : C-----------------------------------------------------------------------
 109 : C      CALL MATGE(VDEP,TRAD,NPA,IBIF0,IBIF1,IBIF2,IBIF3)
 110 : C-----------------------------------------------------------------------
 111 :       
 112 :       DO 11 II=1,6
 113 :        DEFP(II)=DEPST(II)   
 114 : *      SIGT(II)=SIGD(II)
 115 : *      EPST(II)=EPSD(II)
 116 : *      ENDT(II)=ENDD(II)
 117 : 11    CONTINUE
 118 : *      XLIMT=XLIMD
 119 : C-----------------------------------------------------------------------
 120 :       RETURN
 121 :       END
 122 :  
 123 :  
 124 :  
 125 :  
 126 :