1 : C T1 SOURCE CHAT 05/01/13 03:29:34 5004 2 : SUBROUTINE T1(SIG,DSIG,YOUN,ANU,RT1,RT2,H1,RDP,ADP, 3 : 1 XLAM1,XLAM2,IDAM,ANG) 4 : IMPLICIT INTEGER(I-N) 5 : IMPLICIT REAL*8(A-H,O-Z) 6 : -INC CCOPTIO 7 : DIMENSION SIG(3),DSIG(3),SFC(3),DSFC(3),TSFC(3),SFG(3),DSFG(3), 8 : 1 EPC(3),A(4),B(2),IDAM(3),SG(3),DSG(3) 9 : Y=YOUN/(1.D0-ANU*ANU) 10 : IF(RT1.LT.1.E-10) THEN 11 : H1=0.D0 12 : RT1=0.D0 13 : ENDIF 14 : C 15 : C-------------------------------------------- 16 : C BOITE D ECOULEMENT DE TRACTION 1 17 : C-------------------------------------------- 18 : C 19 : C----------------------------------------------- 20 : C ON SE PLACE DANS LE REPERE DE FISSURATION 21 : C----------------------------------------------- 22 : C 22 : CALL CHREP(SIG,SFG,ANG) 22 : CALL CHREP(DSIG,DSFG,ANG) 23 : IF(IIMPI.EQ.9) WRITE(IOIMP,9999) 24 : 9999 FORMAT(1X,'T1 ECOULEMENT SUIVANT TRACTION 1',/) 25 : 1 DL1=(SFG(1)+DSFG(1)-RT1)/(Y-H1) 26 : IF(DL1.LT.0.D0) THEN 27 : IDAM(1)=0 28 : RETURN 29 : ENDIF 30 : C 31 : C---------------------------------------------------------------- 32 : C CAS OU LA LIMITE EN TRACTION RESTE SUPERIEURE A 0 (1000) 33 : C---------------------------------------------------------------- 34 : C 35 : IF((RT1-H1*DL1).GE.-1.E-10) GOTO 1000 36 : C 37 : C-------------------------------------------------------- 38 : C CAS OU LA LIMITE EN TRACTION DEVIENT NULLE 39 : C-------------------------------------------------------- 40 : C 41 : DL1=RT1/H1 42 : EPC(1)=DL1 43 : EPC(2)=0.D0 44 : EPC(3)=0.D0 44 : CALL CPHOOB(EPC,SFC,YOUN,ANU) 45 : X=0.D0 46 : IF(DSFG(1).NE.0.D0) X=(SFC(1)-SFG(1))/DSFG(1) 47 : DO 10 ITYP=1,3 48 : SFC(ITYP)=SFG(ITYP)+X*DSFG(ITYP)-SFC(ITYP) 49 : 10 DSFC(ITYP)=(1.D0-X)*DSFG(ITYP) 49 : CALL CDP(SFC,ADP,RDP,VCDP) 49 : CALL CTRAF(SFC(2),RT2,VCTR2) 49 : CALL CCOAF(SFC(2),XLAM2,VCCO2) 50 : IF(VCTR2.GT.0.D0) GOTO 2000 51 : IF(VCCO2.GT.0.D0) GOTO 3000 52 : IF(VCDP.GT.0.D0) GOTO 4000 53 : RT1=0.D0 54 : H1=0.D0 55 : XLAM1=XLAM1+DL1 56 : DO 20 ITYP=1,3 57 : SFG(ITYP)=SFC(ITYP) 58 : 20 DSFG(ITYP)=DSFC(ITYP) 59 : GOTO 1 60 : 1000 EPC(1)=DL1 61 : EPC(2)=0.D0 62 : EPC(3)=0.D0 62 : CALL CPHOOB(EPC,SFC,YOUN,ANU) 63 : DO 1010 ITYP=1,3 64 : 1010 SFC(ITYP)=SFG(ITYP)+DSFG(ITYP)-SFC(ITYP) 65 : C 66 : C----------------------------------------------- 67 : C ON REGARDE SI LES AUTRES CRITERES NE 68 : C SONT PAS ENDOMMAGES PENDANT L ECOULEMENT 69 : C----------------------------------------------- 70 : C 70 : CALL CDP(SFC,ADP,RDP,VCDP) 70 : CALL CTRAF(SFC(2),RT2,VCTR2) 70 : CALL CCOAF(SFC(2),XLAM2,VCCO2) 71 : C 72 : C------------------------------------ 73 : C CRITERE DE TRACTION ENDOMMAGE 74 : C------------------------------------ 75 : C 76 : IF(VCTR2.GT.0.D0) GOTO 2000 77 : C 78 : C---------------------------------------- 79 : C CRITERE DE COMPRESSION ENDOMMAGE 80 : C---------------------------------------- 81 : C 82 : IF(VCCO2.GT.0.D0) GOTO 3000 83 : C 84 : C------------------------------------ 85 : C CRITERE DE DRUCKER ENDOMMAGE 86 : C------------------------------------ 87 : C 88 : IF(VCDP.GT.0.D0) GOTO 4000 89 : C 90 : C-------------------------------------------------------- 91 : C AUCUN AUTRE CRITERE ENDOMMAGE PENDANT L ECOULEMENT 92 : C-------------------------------------------------------- 93 : C 93 : CALL CHREP(SFC,SIG,-ANG) 94 : DO 1020 ITYP=1,3 95 : 1020 DSIG(ITYP)=0.D0 96 : RT1=RT1-H1*DL1 97 : IF(RT1.LE.1.E-10) THEN 98 : H1=0.D0 99 : RT1=0.D0 100 : ENDIF 101 : IDAM(1)=0 102 : IDAM(2)=0 103 : IDAM(3)=0 104 : XLAM1=XLAM1+DL1 105 : RETURN 106 : C 107 : C---------------------------------------- 108 : C CAS OU EN ENDOMMAGE LA TRACTION2 109 : C---------------------------------------- 110 : C 111 : 2000 A(1)=DSFG(1) 112 : A(3)=DSFG(2) 113 : A(2)=-Y+H1 114 : A(4)=-Y*ANU 115 : DET=A(1)*A(4)-A(2)*A(3) 116 : IF(DET.EQ.0.D0) THEN 117 : IDAM(1)=1 118 : IDAM(2)=1 119 : IDAM(3)=0 120 : RETURN 121 : ENDIF 122 : B(1)=RT1-SFG(1) 123 : B(2)=RT2-SFG(2) 123 : CALL SYLIN2(A,B,X,DL1) 124 : EPC(1)=DL1 125 : EPC(2)=0.D0 126 : EPC(3)=0.D0 126 : CALL CPHOOB(EPC,SFC,YOUN,ANU) 127 : DO 2010 ITYP=1,3 128 : DSFC(ITYP)=(1.D0-X)*DSFG(ITYP) 129 : 2010 SFC(ITYP)=SFG(ITYP)+X*DSFG(ITYP)-SFC(ITYP) 129 : CALL CDP(SFC,ADP,RDP,VCDP) 130 : IF(VCDP.GT.0.D0) GOTO 4000 131 : RT1=RT1-H1*DL1 132 : IF(RT1.LE.1.E-10) THEN 133 : H1=0.D0 134 : RT1=0.D0 135 : ENDIF 136 : IDAM(1)=1 137 : IDAM(2)=1 138 : IDAM(3)=0 139 : XLAM1=XLAM1+DL1 139 : CALL CHREP(DSFC,DSIG,-ANG) 139 : CALL CHREP(SFC,SIG,-ANG) 140 : RETURN 141 : C 142 : C---------------------------------------- 143 : C CAS OU EN ENDOMMAGE LA COMPRESSION 2 144 : C---------------------------------------- 145 : C 146 : 3000 A(1)=DSFG(1) 147 : A(3)=DSFG(2) 148 : A(2)=-Y+H1 149 : A(4)=-Y*ANU 150 : DET=A(1)*A(4)-A(2)*A(3) 151 : IF(DET.EQ.0.D0) THEN 152 : IDAM(1)=1 153 : IDAM(2)=-1 154 : IDAM(3)=0 155 : RETURN 156 : ENDIF 157 : B(1)=RT1-SFG(1) 158 : B(2)=-SFG(2) 158 : CALL SYLIN2(A,B,X,DL1) 159 : EPC(1)=DL1 160 : EPC(2)=0.D0 161 : EPC(3)=0.D0 161 : CALL CPHOOB(EPC,SFC,YOUN,ANU) 162 : DO 3010 ITYP=1,3 163 : DSFC(ITYP)=(1.D0-X)*DSFG(ITYP) 164 : 3010 SFC(ITYP)=SFG(ITYP)+X*DSFG(ITYP)-SFC(ITYP) 164 : CALL CDP(SFC,ADP,RDP,VCDP) 165 : IF(VCDP.GT.0.D0) GOTO 4000 166 : RT1=RT1-H1*DL1 167 : IF(RT1.LE.1.E-10) THEN 168 : H1=0.D0 169 : RT1=0.D0 170 : ENDIF 171 : IDAM(1)=1 172 : IDAM(2)=-1 173 : IDAM(3)=0 174 : XLAM1=XLAM1+DL1 174 : CALL CHREP(DSFC,DSIG,-ANG) 174 : CALL CHREP(SFC,SIG,-ANG) 175 : RETURN 176 : C 177 : C---------------------------------------- 178 : C CAS OU EN ENDOMMAGE LE DP 179 : C---------------------------------------- 180 : C 181 : 4000 CONTINUE 182 : EPC(1)=1.D0 183 : EPC(2)=0.D0 184 : EPC(3)=0.D0 184 : CALL CPHOOB(EPC,SFC,YOUN,ANU) 185 : E=Y-H1 186 : F=(SFG(1)-RT1)/E 187 : G=DSFG(1)/E 188 : DO 4010 ITYP=1,3 189 : SG(ITYP)=SFG(ITYP)-F*SFC(ITYP) 190 : 4010 DSG(ITYP)=DSFG(ITYP)-G*SFC(ITYP) 190 : CALL XDP(SG,DSG,RDP,ADP,X,ITEST) 191 : IF(ITEST.EQ.1) THEN 192 : IDAM(1)=1 193 : IDAM(2)=0 194 : IDAM(3)=1 195 : RETURN 196 : ENDIF 197 : DL1=(SFG(1)+X*DSFG(1)-RT1)/(Y-H1) 198 : EPC(1)=DL1 199 : EPC(2)=0.D0 200 : EPC(3)=0.D0 201 : RT1=RT1-H1*DL1 202 : IF(RT1.LE.1.E-10) THEN 203 : H1=0.D0 204 : RT1=0.D0 205 : ENDIF 205 : CALL CPHOOB(EPC,SFC,YOUN,ANU) 206 : DO 4020 ITYP=1,3 207 : SFC(ITYP)=X*DSFG(ITYP)-SFC(ITYP)+SFG(ITYP) 208 : 4020 DSFC(ITYP)=DSFG(ITYP)*(1.D0-X) 209 : XLAM1=XLAM1+DL1 209 : CALL CHREP(DSFC,DSIG,-ANG) 209 : CALL CHREP(SFC,SIG,-ANG) 210 : IDAM(1)=1 211 : IDAM(2)=0 212 : IDAM(3)=1 213 : RETURN 214 : END 215 :
© Cast3M 2003 - All rights reserved.
Disclaimer