1 : C CLENDO SOURCE CHAT 05/01/12 22:04:48 5004 2 : SUBROUTINE CLENDO (DDP,YUNG,XNU,ET,CISA,IENDO,KERRE) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8(A-H,O-Z) 5 : DIMENSION ET(*),CISA(*),DDP(6,6) 6 : C 7 : C*********************************************************************** 8 : C***************** CALCUL DU MODULE D'ENDOMMAGEMENT ******************** 9 : C*********************************************************************** 10 : C 11 : C INITIALISATION 12 : C 12 : CALL ZDANUL(DDP,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 : COEF=YUNG/(UPUN*UMDN) 23 : C 24 : C*********************************************************************** 25 : C************************ PAS DE FISSURATION *************************** 26 : C*********************************************************************** 27 : C 28 : DDP(1,1)=UMUN*COEF 29 : DDP(1,2)=XNU*COEF 30 : DDP(1,3)=DDP(1,2) 31 : DDP(2,1)=DDP(1,2) 32 : DDP(2,2)=DDP(1,1) 33 : DDP(2,3)=DDP(1,2) 34 : DDP(3,1)=DDP(1,2) 35 : DDP(3,2)=DDP(1,2) 36 : DDP(3,3)=DDP(1,1) 37 : DDP(4,4)=YUNG/(2.D0*UPUN) 38 : DDP(5,5)=DDP(4,4) 39 : DDP(6,6)=DDP(4,4) 40 : IF(IENDO.EQ.0) RETURN 41 : C 42 : GO TO (100,200,300,400,500,600,700),IENDO 43 : KERRE=468 44 : RETURN 45 : C 46 : C*********************************************************************** 47 : C******************** FISSURATION EN X SEULEMENT *********************** 48 : C*********************************************************************** 49 : C 50 : 100 DENO=(EMET1*UMUN)+(ET1*UPUN*UMDN) 51 : USDENO=1.D0/DENO 52 : DDP(1,1)=YUNG*ET1*UMUN*USDENO 53 : DDP(1,2)=XNU*YUNG*ET1*USDENO 54 : DDP(1,3)=DDP(1,2) 55 : DDP(2,1)=DDP(1,2) 56 : DDP(2,2)=YUNG*((EMET1/UPUN)+(ET1*UMUN))*USDENO 57 : DDP(2,3)=XNU*YUNG*((EMET1/UPUN)+ET1)*USDENO 58 : DDP(3,1)=DDP(1,3) 59 : DDP(3,2)=DDP(2,3) 60 : DDP(3,3)=DDP(2,2) 61 : DDP(4,4)=CISA(1) 62 : DDP(6,6)=CISA(3) 63 : RETURN 64 : C 65 : C*********************************************************************** 66 : C******************** FISSURATION EN Y SEULEMENT *********************** 67 : C*********************************************************************** 68 : C 69 : 200 DENO=(EMET2*UMUN)+(ET2*UPUN*UMDN) 70 : USDENO=1.D0/DENO 71 : DDP(1,1)=YUNG*((EMET2/UPUN)+(ET2*UMUN))*USDENO 72 : DDP(1,2)=XNU*YUNG*ET2*USDENO 73 : DDP(1,3)=XNU*YUNG*((EMET2/UPUN)+ET2)*USDENO 74 : DDP(2,1)=DDP(1,2) 75 : DDP(2,2)=YUNG*ET2*UMUN*USDENO 76 : DDP(2,3)=DDP(1,2) 77 : DDP(3,1)=DDP(1,3) 78 : DDP(3,2)=DDP(2,3) 79 : DDP(3,3)=DDP(1,1) 80 : DDP(4,4)=CISA(1) 81 : DDP(5,5)=CISA(2) 82 : RETURN 83 : C 84 : C*********************************************************************** 85 : C******************** FISSURATION EN Z SEULEMENT *********************** 86 : C*********************************************************************** 87 : C 88 : 400 DENO=(EMET3*UMUN)+(ET3*UPUN*UMDN) 89 : USDENO=1.D0/DENO 90 : DDP(1,1)=YUNG*((EMET3/UPUN)+(ET3*UMUN))*USDENO 91 : DDP(1,2)=XNU*YUNG*((EMET3/UPUN)+ET3)*USDENO 92 : DDP(1,3)=XNU*YUNG*ET3*USDENO 93 : DDP(2,1)=DDP(1,2) 94 : DDP(2,2)=DDP(1,1) 95 : DDP(2,3)=DDP(1,3) 96 : DDP(3,1)=DDP(1,3) 97 : DDP(3,2)=DDP(2,3) 98 : DDP(3,3)=YUNG*ET3*UMUN*USDENO 99 : DDP(5,5)=CISA(2) 100 : DDP(6,6)=CISA(3) 101 : RETURN 102 : C 103 : C*********************************************************************** 104 : C******************** FISSURATION EN XY ENSEMBLE *********************** 105 : C*********************************************************************** 106 : C 107 : 300 DENO=(EMET1*EMET2)+(ET1*EMET2*UPUN*UMUN)+(EMET1*ET2*UPUN*UMUN)+ 108 : . (ET1*ET2*UPUN*UPUN*UMDN) 109 : USDENO=1.D0/DENO 110 : DDP(1,1)=YUNG*ET1*(EMET2+(ET2*UPUN*UMUN))*USDENO 111 : DDP(1,2)=XNU*YUNG*ET1*ET2*UPUN*USDENO 112 : DDP(1,3)=XNU*YUNG*ET1*(EMET2+(ET2*UPUN))*USDENO 113 : DDP(2,1)=DDP(1,2) 114 : DDP(2,2)=YUNG*ET2*(EMET1+(ET1*UPUN*UMUN))*USDENO 115 : DDP(2,3)=XNU*YUNG*ET2*(EMET1+(ET1*UPUN))*USDENO 116 : DDP(3,1)=DDP(1,3) 117 : DDP(3,2)=DDP(2,3) 118 : DDP(3,3)=YUNG*((EMET1*EMET2)+(ET1*EMET2)+(EMET1*ET2)+ 119 : . (ET1*ET2*UPUN*UMUN))*USDENO 120 : DDP(4,4)=CISA(1) 121 : DDP(5,5)=CISA(2) 122 : DDP(6,6)=CISA(3) 123 : RETURN 124 : C 125 : C*********************************************************************** 126 : C******************** FISSURATION EN YZ ENSEMBLE *********************** 127 : C*********************************************************************** 128 : C 129 : 600 DENO=(EMET2*EMET3)+(ET2*EMET3*UPUN*UMUN)+(EMET2*ET3*UPUN*UMUN)+ 130 : . (ET2*ET3*UPUN*UPUN*UMDN) 131 : USDENO=1.D0/DENO 132 : DDP(1,1)=YUNG*((EMET2*EMET3)+(ET2*EMET3)+(EMET2*ET3)+ 133 : . (ET2*ET3*UPUN*UMUN))*USDENO 134 : DDP(1,2)=XNU*YUNG*ET2*(EMET3+(ET3*UPUN))*USDENO 135 : DDP(1,3)=XNU*YUNG*ET3*(EMET2+(ET2*UPUN))*USDENO 136 : DDP(2,1)=DDP(1,2) 137 : DDP(2,2)=YUNG*ET2*(EMET3+(ET3*UPUN*UMUN))*USDENO 138 : DDP(2,3)=XNU*YUNG*ET2*ET3*UPUN*USDENO 139 : DDP(3,1)=DDP(1,3) 140 : DDP(3,2)=DDP(2,3) 141 : DDP(3,3)=YUNG*ET3*(EMET2+(ET2*UPUN*UMUN))*USDENO 142 : DDP(4,4)=CISA(1) 143 : DDP(5,5)=CISA(2) 144 : DDP(6,6)=CISA(3) 145 : RETURN 146 : C 147 : C*********************************************************************** 148 : C******************** FISSURATION EN ZX ENSEMBLE *********************** 149 : C*********************************************************************** 150 : C 151 : 500 DENO=(EMET3*EMET1)+(ET3*EMET1*UPUN*UMUN)+(EMET3*ET1*UPUN*UMUN)+ 152 : . (ET3*ET1*UPUN*UPUN*UMDN) 153 : USDENO=1.D0/DENO 154 : DDP(1,1)=YUNG*ET1*(EMET3+(ET3*UPUN*UMUN))*USDENO 155 : DDP(1,2)=XNU*YUNG*ET1*(EMET3+(ET3*UPUN))*USDENO 156 : DDP(1,3)=XNU*YUNG*ET3*ET1*UPUN*USDENO 157 : DDP(2,1)=DDP(1,2) 158 : DDP(2,2)=YUNG*((EMET3*EMET1)+(ET3*EMET1)+(EMET3*ET1)+ 159 : . (ET3*ET1*UPUN*UMUN))*USDENO 160 : DDP(2,3)=XNU*YUNG*ET3*(EMET1+(ET1*UPUN))*USDENO 161 : DDP(3,1)=DDP(1,3) 162 : DDP(3,2)=DDP(2,3) 163 : DDP(3,3)=YUNG*ET3*(EMET1+(ET1*UPUN*UMUN))*USDENO 164 : DDP(4,4)=CISA(1) 165 : DDP(5,5)=CISA(2) 166 : DDP(6,6)=CISA(3) 167 : RETURN 168 : C 169 : C*********************************************************************** 170 : C******************** FISSURATION EN XYZ ENSEMBLE ********************** 171 : C*********************************************************************** 172 : C 173 : 700 DENO=(EMET1*EMET2*EMET3)+ 174 : . (ET1*EMET2*EMET3)+(EMET1*ET2*EMET3)+(EMET1*EMET2*ET3)+ 175 : . (UPUN*UMUN*((EMET1*ET2*ET3)+(ET1*EMET2*ET3)+(ET1*ET2*EMET3)))+ 176 : . (UPUN*UPUN*UMDN*(ET1*ET2*ET3)) 177 : USDENO=1.D0/DENO 178 : DDP(1,1)=((EMET2*EMET3)+(ET2*EMET3)+(EMET2*ET3)+ 179 : . (ET2*ET3*UPUN*UMUN))*YUNG*ET1*USDENO 180 : DDP(1,2)=XNU*YUNG*ET1*ET2*(EMET3+(ET3*UPUN))*USDENO 181 : DDP(1,3)=XNU*YUNG*ET1*ET3*(EMET2+(ET2*UPUN))*USDENO 182 : DDP(2,1)=DDP(1,2) 183 : DDP(2,2)=((EMET1*EMET3)+(ET1*EMET3)+(EMET1*ET3)+ 184 : . (ET1*ET3*UPUN*UMUN))*YUNG*ET2*USDENO 185 : DDP(2,3)=XNU*YUNG*ET2*ET3*(EMET1+(ET1*UPUN))*USDENO 186 : DDP(3,1)=DDP(1,3) 187 : DDP(3,2)=DDP(2,3) 188 : DDP(3,3)=((EMET1*EMET2)+(ET1*EMET2)+(EMET1*ET2)+ 189 : . (ET1*ET2*UPUN*UMUN))*YUNG*ET3*USDENO 190 : DDP(4,4)=CISA(1) 191 : DDP(5,5)=CISA(2) 192 : DDP(6,6)=CISA(3) 193 : RETURN 194 : C 195 : END 196 :
© Cast3M 2003 - All rights reserved.
Disclaimer