1 : C CDRAGO SOURCE KICH 09/09/11 21:15:42 6505 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----------------------------------------------------------------------- 86 : * CALL ARED3M(SIG0,EPIN0,VAR0(2),XLIMT, 87 : * * SIGF,EPINF,VARF(2),XLIMD, 88 : * * DEPST,VDEP,VDEP1,NPA,VAL,VP,EPSPLS,VAL1,VP1, 89 : * * VAL2,VAL3,VP2,VP3,EPSTT,FDDD) 89 : CALL ARED3M(SIG0,EPIN0,VAR0(2),XLIMT, 90 : * SIGF,EPINF,VARF(2),XLIMD, 91 : * DEPST,AAA(25),AAA(61),NPA, 92 : * AAA(1),AAA(22),AAA(97),AAA(10),AAA(19), 93 : * AAA(103),AAA(112),AAA(121),AAA(124), 94 : * AAA(127),FDDD,AAA(133),AAA(142),AAA(148), 95 : * AAA(184),AAA(220),AAA(226),AAA(307), 96 : * AAA(388),AAA(394),AAA(400),KM,AAA) 96 : * CALL ARED3M(SIGT,EPST,ENDT,XLIMT, 97 : * * SIGD,EPSD,ENDD,XLIMD, 98 : * * DEPS,VDEP,VDEP1,NPA,VAL,VP,EPSPLS,VAL1,VP1, 99 : * * VAL2,VAL3,VP2,VP3,EPSTT,FDDD) 100 : VARF(1)=VARF(2)+VARF(3)+VARF(4) 101 : C----------------------------------------------------------------------- 101 : C CALL MATGE(VDEP,TRAD,NPA,IBIF0,IBIF1,IBIF2,IBIF3) 102 : C----------------------------------------------------------------------- 103 : 104 : DO 11 II=1,6 105 : DEFP(II)=DEPST(II) 106 : * SIGT(II)=SIGD(II) 107 : * EPST(II)=EPSD(II) 108 : * ENDT(II)=ENDD(II) 109 : 11 CONTINUE 110 : * XLIMT=XLIMD 111 : C----------------------------------------------------------------------- 112 : segsup DRA0 113 : RETURN 114 : END 115 : 116 : 117 : 118 : 119 : 120 : 121 : 122 : 123 : 124 :
© Cast3M 2003 - All rights reserved.
Disclaimer