1 : C CLOUVE SOURCE CHAT 05/01/12 22:08:47 5004 2 : SUBROUTINE CLOUVE (AAA,YUNG,XNU,ET,CISA,IENDO,KERRE) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8(A-H,O-Z) 5 : DIMENSION ET(*),CISA(*),AAA(6,6),UMBETA(3) 6 : C 7 : C*********************************************************************** 8 : C***************** CALCUL DU MODULE DE FISSURATION ********************* 9 : C*********************************************************************** 10 : C 11 : C INITIALISATION 12 : C 12 : CALL ZDANUL(AAA,36) 13 : ET1=ET(1) 14 : ET2=ET(2) 15 : ET3=ET(3) 16 : UPUN=1.D0+XNU 17 : UMUN=1.D0-XNU 18 : UMDN=1.D0-2.D0*XNU 19 : EMET1=YUNG-ET1 20 : EMET2=YUNG-ET2 21 : EMET3=YUNG-ET3 22 : CISEL=YUNG/(2.D0*UPUN) 23 : UMBETA(1)=1.D0-(CISA(1)/CISEL) 24 : UMBETA(2)=1.D0-(CISA(2)/CISEL) 25 : UMBETA(3)=1.D0-(CISA(3)/CISEL) 26 : C 27 : C*********************************************************************** 28 : C************************ PAS DE FISSURATION *************************** 29 : C*********************************************************************** 30 : C 31 : IF(IENDO.EQ.0) RETURN 32 : C 33 : GO TO (100,200,300,400,500,600,700),IENDO 34 : KERRE=468 35 : RETURN 36 : C 37 : C*********************************************************************** 38 : C******************** FISSURATION EN X SEULEMENT *********************** 39 : C*********************************************************************** 40 : C 41 : 100 DENO=(EMET1*UMUN)+(ET1*UPUN*UMDN) 42 : USDENO=1.D0/DENO 43 : AAA(1,1)=EMET1*UMUN*USDENO 44 : AAA(1,2)=EMET1*XNU*USDENO 45 : AAA(1,3)=EMET1*XNU*USDENO 46 : AAA(4,4)=UMBETA(1) 47 : AAA(6,6)=UMBETA(3) 48 : RETURN 49 : C 50 : C*********************************************************************** 51 : C******************** FISSURATION EN Y SEULEMENT *********************** 52 : C*********************************************************************** 53 : C 54 : 200 DENO=(EMET2*UMUN)+(ET2*UPUN*UMDN) 55 : USDENO=1.D0/DENO 56 : AAA(2,1)=EMET2*XNU*USDENO 57 : AAA(2,2)=EMET2*UMUN*USDENO 58 : AAA(2,3)=EMET2*XNU*USDENO 59 : AAA(4,4)=UMBETA(1) 60 : AAA(5,5)=UMBETA(2) 61 : RETURN 62 : C 63 : C*********************************************************************** 64 : C******************** FISSURATION EN Z SEULEMENT *********************** 65 : C*********************************************************************** 66 : C 67 : 400 DENO=(EMET3*UMUN)+(ET3*UPUN*UMDN) 68 : USDENO=1.D0/DENO 69 : AAA(3,1)=EMET3*XNU*USDENO 70 : AAA(3,2)=EMET3*XNU*USDENO 71 : AAA(3,3)=EMET3*UMUN*USDENO 72 : AAA(5,5)=UMBETA(2) 73 : AAA(6,6)=UMBETA(3) 74 : RETURN 75 : C 76 : C*********************************************************************** 77 : C******************** FISSURATION EN XY ENSEMBLE *********************** 78 : C*********************************************************************** 79 : C 80 : 300 DENO=(EMET1*EMET2)+(ET1*EMET2*UPUN*UMUN)+(EMET1*ET2*UPUN*UMUN)+ 81 : . (ET1*ET2*UPUN*UPUN*UMDN) 82 : USDENO=1.D0/DENO 83 : AAA(1,1)=EMET1*(EMET2+(ET2*UPUN*UMUN))*USDENO 84 : AAA(1,2)=EMET1*(XNU*ET2*UPUN)*USDENO 85 : AAA(1,3)=XNU*EMET1*(EMET2+(ET2*UPUN))*USDENO 86 : AAA(2,1)=EMET2*(XNU*ET1*UPUN)*USDENO 87 : AAA(2,2)=EMET2*(EMET1+(ET1*UPUN*UMUN))*USDENO 88 : AAA(2,3)=XNU*EMET2*(EMET1+(ET1*UPUN))*USDENO 89 : AAA(4,4)=UMBETA(1) 90 : AAA(5,5)=UMBETA(2) 91 : AAA(6,6)=UMBETA(3) 92 : RETURN 93 : C 94 : C*********************************************************************** 95 : C******************** FISSURATION EN YZ ENSEMBLE *********************** 96 : C*********************************************************************** 97 : C 98 : 600 DENO=(EMET2*EMET3)+(ET2*EMET3*UPUN*UMUN)+(EMET2*ET3*UPUN*UMUN)+ 99 : . (ET2*ET3*UPUN*UPUN*UMDN) 100 : USDENO=1.D0/DENO 101 : AAA(2,1)=XNU*EMET2*(EMET3+(ET3*UPUN))*USDENO 102 : AAA(2,2)=EMET2*(EMET3+(ET3*UPUN*UMUN))*USDENO 103 : AAA(2,3)=EMET2*(XNU*ET3*UPUN)*USDENO 104 : AAA(3,1)=XNU*EMET3*(EMET2+(ET2*UPUN))*USDENO 105 : AAA(3,2)=EMET3*(XNU*ET2*UPUN)*USDENO 106 : AAA(3,3)=EMET3*(EMET2+(ET2*UPUN*UMUN))*USDENO 107 : AAA(4,4)=UMBETA(1) 108 : AAA(5,5)=UMBETA(2) 109 : AAA(6,6)=UMBETA(3) 110 : RETURN 111 : C 112 : C*********************************************************************** 113 : C******************** FISSURATION EN ZX ENSEMBLE *********************** 114 : C*********************************************************************** 115 : C 116 : 500 DENO=(EMET3*EMET1)+(ET3*EMET1*UPUN*UMUN)+(EMET3*ET1*UPUN*UMUN)+ 117 : . (ET3*ET1*UPUN*UPUN*UMDN) 118 : USDENO=1.D0/DENO 119 : AAA(1,1)=EMET1*(EMET3+(ET3*UPUN*UMUN))*USDENO 120 : AAA(1,2)=XNU*EMET1*(EMET3+(ET3*UPUN))*USDENO 121 : AAA(1,3)=EMET1*(XNU*ET3*UPUN)*USDENO 122 : AAA(3,1)=EMET3*(XNU*ET1*UPUN)*USDENO 123 : AAA(3,2)=XNU*EMET3*(EMET1+(ET1*UPUN))*USDENO 124 : AAA(3,3)=EMET3*(EMET1+(ET1*UPUN*UMUN))*USDENO 125 : AAA(4,4)=UMBETA(1) 126 : AAA(5,5)=UMBETA(2) 127 : AAA(6,6)=UMBETA(3) 128 : RETURN 129 : C 130 : C*********************************************************************** 131 : C******************** FISSURATION EN XYZ ENSEMBLE ********************** 132 : C*********************************************************************** 133 : C 134 : 700 DENO=(EMET1*EMET2*EMET3)+ 135 : . (ET1*EMET2*EMET3)+(EMET1*ET2*EMET3)+(EMET1*EMET2*ET3)+ 136 : . (UPUN*UMUN*((EMET1*ET2*ET3)+(ET1*EMET2*ET3)+(ET1*ET2*EMET3)))+ 137 : . (UPUN*UPUN*UMDN*(ET1*ET2*ET3)) 138 : USDENO=1.D0/DENO 139 : AAA(1,1)=((EMET1*EMET2*EMET3)+(EMET1*ET2*EMET3)+(EMET1*EMET2*ET3)+ 140 : . (EMET1*ET2*ET3*UPUN*UMUN))*USDENO 141 : AAA(1,2)=XNU*EMET1*ET2*(EMET3+(ET3*UPUN))*USDENO 142 : AAA(1,3)=XNU*EMET1*ET3*(EMET2+(ET2*UPUN))*USDENO 143 : AAA(2,1)=XNU*EMET2*ET1*(EMET3+(ET3*UPUN))*USDENO 144 : AAA(2,2)=((EMET1*EMET2*EMET3)+(ET1*EMET2*EMET3)+(EMET1*EMET2*ET3)+ 145 : . (ET1*EMET2*ET3*UPUN*UMUN))*USDENO 146 : AAA(2,3)=XNU*EMET2*ET3*(EMET1+(ET1*UPUN))*USDENO 147 : AAA(3,1)=XNU*EMET3*ET1*(EMET2+(ET2*UPUN))*USDENO 148 : AAA(3,2)=XNU*EMET3*ET2*(EMET1+(ET1*UPUN))*USDENO 149 : AAA(3,3)=((EMET1*EMET2*EMET3)+(ET1*EMET2*EMET3)+(EMET1*ET2*EMET3)+ 150 : . (ET1*ET2*EMET3*UPUN*UMUN))*USDENO 151 : AAA(4,4)=UMBETA(1) 152 : AAA(5,5)=UMBETA(2) 153 : AAA(6,6)=UMBETA(3) 154 : RETURN 155 : C 156 : END 157 :
© Cast3M 2003 - All rights reserved.
Disclaimer