1 : C ECBPEL SOURCE CHAT 05/01/12 23:20:26 5004 2 : SUBROUTINE ECBPEL(SIG0,NSTRS,DEPST,VAR0,XMAT,NCOMAT,XCAR,ICARA, 3 : 1NVARI,SIGF,VARF,DEFP,MFR1,DDAUX,CMATE,VALMAT,VALCAR,N2EL,N2PTEL, 4 : 2 NBPGAU,IFOU,IB,IGAU,EPAIST,MELE,NPINT,NBGMAT,NELMAT,SECT, 5 : 3 LHOOK,TXR,XLOC,XGLOB,D1HOOK,ROTHOO,DDHOMU,CRIGI,DSIGT,KERRE,DT) 6 : C 7 : C ================================================================== 8 : C CE SOUS-PROGRAMME EST APPELE DANS "ECOUL2" VIA ECOU25. 9 : C IL PREPARE L'INTEGRATION DE LA FORMULE DE RELAXATION 10 : C DE CABLES DE PRECONTRAINTE SELON LA METHODE DU TEMPS 11 : C EQUIVALENT PRECONISEE PAR LE BPEL 91 12 : C 13 : C ENTREES: 14 : C ------- 15 : C NSTRS = NBR. DE COMPOSANTES DES CONTR. OU DES DEFORM. 16 : C SIG0(NSTRS) = CONTR. AU DEBUT DU PAS D'INTEGRATION 17 : C DEPST(NSTRS) = INCREMENT DES DEFORM. CALCULE ELASTIQUEMENT 18 : C A PARTIR DE L'INCREMENT DES DEFORM. TOTALES 19 : C NVARI = NBR. DE VARIABLES INTERNES 20 : C VAR0(NVARI) = VARIABLES INTERNES AU DEBUT DU PAS D'INTEGRATION 21 : C CE TABLEAU NE CONTIENT QU'UN INDICE CORRESPONDANT A 22 : C - LA DEFORM. INELASTIQUE CUMULEE 23 : C NCOMAT = NBR. DE CARACTERISTIQUES MECANIQUES DU MATERIAU 24 : C XMAT(NCOMAT) = CARACTERISTIQUES MECANIQUES DU MATERIAU 25 : C MFR = INDICE DE LA FORMULATION MECANIQUE 26 : C ICARA = NBR. DE CARACT. GEOMETRIQUES DES ELEMENTS FINIS 27 : C XCAR(ICARA) = CARACT. GEOMETRIQUES DES ELEMENTS FINIS 28 : * DDAUX = MATRICE DE HOOKE ELASTIQUE 29 : * NSTRS = NBRE DE COMPOSANTES DES DEFORMATIONS 30 : * CMATE = NOM DU MATERIAU 31 : * VALMAT= TABLEAU DE CARACTERISTIQUES DU MATERIAU 32 : * VALCAR= TABLEAU DE CARACTERISTIQUES GEOMETRIQUES 33 : * N2EL = NBRE D ELEMENTS DANS SEGMENT DE HOOKE 34 : * N2PTEL= NBRE DE POINTS DANS SEGMENT DE HOOKE 35 : * IFOU = OPTION DE CALCUL 36 : * IB = NUMERO DE L ELEMENT COURANT 37 : * IGAU = NUMERO DU POINT COURANT 38 : * EPAIST= EPAISSEUR 39 : * NBPGAU= NBRE DE POINTS DE GAUSS 40 : * MELE = NUMERO DE L ELEMENT FINI 41 : * NPINT = NBRE DE POINTS D INTEGRATION 42 : * NBGMAT= NBRE DE POINTS DANS SEGMENT DE CARACTERISTIQUES 43 : * NELMAT= NBRE D ELEMENTS DANS SEGMENT DE CARACTERISTIQUES 44 : * SECT = SECTION 45 : * LHOOK = TAILLE DE LA MATRICE DE HOOKE 46 : * TXR,XLOC,XGLOB,D1HOOK,ROTHOO,DDHOMU,CRIGI = TABLEAUX UTILISES 47 : * POUR LE CALCUL DE LA MATRICE DE HOOKE 48 : C 49 : C SORTIES: 50 : C ------- 51 : C SIGF(NSTRS)= CONTR. A LA FIN DU PAS D'INTEGRATION 52 : C VARF(NVARI)= VARIABLES INTERNES A LA FIN DU PAS D'INTEGRATION 53 : C DEFP(NSTRS)= INCREMENT DES DEFORM. PLASTIQUES A LA FIN DU PAS 54 : C D'INTEGRATION 55 : C KERRE = INDICE QUI REGIT LES ERREURS 56 : C = 77 SI LA DEFORM. PLAST. CUMULEE ENDOMMAGEE (2IEME VAR. 57 : C INT.) EST EN DEHORS DE LA COURBE DE TRACTION, DS. 58 : C LE CAS DE L'ECROUISSAGE ET DE L'ENDOMM. ISOTROPES. 59 : C CECI PEUT SE PRODUIRE SUITE A L'APPEL A "CRIDAM" 60 : C = 99 SI LA FORMULATION MECANIQUE N'EST PAS DISPONIBLE 61 : C POUR LE MODELE CONSIDERE OU S'IL Y A INCOMPATIBILITE 62 : C ENTRE MFR ET IFOUR 63 : C ================================================================== 64 : C ICI IL FAUT PROGRAMMER EN FORTRAN PUR 65 : C =================================================================== 66 : C 67 : IMPLICIT INTEGER(I-N) 68 : IMPLICIT REAL*8(A-H,O-Z) 69 : -INC CCOPTIO 70 : DIMENSION SIG0(*),DEPST(*),VAR0(*),XMAT(*),XCAR(*),SIGF(*), 71 : & VARF(*),DEFP(*) 72 : DIMENSION RSIG0(6),RDSIGT(6),RDEFP(6),RSIGF(6),DSIGT(*) 73 : DIMENSION VALMAT(*),VALCAR(*) 74 : DIMENSION TXR(IDIM,*),CRIGI(12) 75 : DIMENSION DDAUX(LHOOK,*),DDHOMU(LHOOK,*) 76 : DIMENSION XLOC(IDIM,*),XGLOB(IDIM,*) 77 : DIMENSION D1HOOK(LHOOK,*),ROTHOO(LHOOK,*) 78 : CHARACTER*8 CMATE 79 : REAL*8 K1,K2,MU 80 : C 81 : C CALCUL DU TEMPS EQUIVALENT 82 : C 83 : EACIER = XMAT(1) 84 : AIRE = XCAR(1) 85 : TENS0 = SIG0(1) 86 : SIGP0 = (TENS0/AIRE) + (EACIER*DEPST(1)) 87 : EPIN0 = VAR0(1) 88 : DSIGP0 = EPIN0*EACIER 89 : RO1000 = XMAT(5) 90 : FPRG = XMAT(6) 91 : TCHA = XMAT(7) 92 : 93 : IF (RO1000 .LE. 2.5D0) THEN 94 : K1 = 6.D-3 95 : K2 = 1.1D0 96 : ELSE 97 : K1 = 8.D-3 98 : K2 = 1.25D0 99 : IF (RO1000 .GE. 8.D0) THEN 100 : K2 = 1.5D0 101 : ENDIF 102 : ENDIF 103 : 104 : * WRITE(6,*) 'DT = ',DT 105 : * WRITE(6,*) 'DEPST(1) = ',DEPST(1) 106 : * WRITE(6,*) 'SIG0(1) = ',SIG0(1) 107 : * WRITE(6,*) 'VAR0(1) = ',VAR0(1) 108 : * WRITE(6,*) 'SIGP0 = ',SIGP0 109 : * WRITE(6,*) 'DSIGP0 = ',DSIGP0 110 : 111 : ** WRITE(6,*) 'EACIER = ',EACIER 112 : ** WRITE(6,*) 'AIRE = ',AIRE 113 : ** WRITE(6,*) 'TENS0 = ',TENS0 114 : ** WRITE(6,*) 'SIGP0 = ',SIGP0 115 : ** WRITE(6,*) 'EPIN0 = ',EPIN0 116 : ** WRITE(6,*) 'DSIGP0 = ',DSIGP0 117 : ** WRITE(6,*) 'RO1000 = ',RO1000 118 : ** WRITE(6,*) 'FPRG = ',FPRG 119 : ** WRITE(6,*) 'TCHA = ',TCHA 120 : ** WRITE(6,*) 'K1 = ',K1 121 : ** WRITE(6,*) 'K2 = ',K2 122 : 122 : CALL TEQUIV(K1,RO1000,TCHA,FPRG,K2,SIGP0,DSIGP0,TEQ) 123 : 124 : * WRITE(6,*) 'TEQ = ',TEQ 125 : MU = (SIGP0+DSIGP0)/FPRG 126 : 127 : * WRITE(6,*) 'MU = ',MU 128 : * WRITE(6,*) 'DT = ',DT 129 : 130 : DEPIN0 = (K1*RO1000*((24.D0*(TEQ+DT-TCHA)/1.D3) 131 : & **(3.D0*(1.D0-MU)/4.D0))* 132 : & (EXP((10.D0*MU-7.5D0)/K2))*(SIGP0+DSIGP0) - DSIGP0) 133 : & /EACIER 134 : 135 : DELAS = DEPST(1) - DEPIN0 136 : DEFP(1) = DEPIN0 137 : VARF(1) = VAR0(1) + DEFP(1) 138 : * SIGF(1) = (AIRE*SIGP0) + (EACIER*AIRE*DEFP(1)) 139 : * SIGF(1) = AIRE*SIGP0 140 : SIGF(1) = SIG0(1) + AIRE*EACIER*DELAS 141 : SIGPF = SIGF(1)/AIRE 142 : 143 : * WRITE(6,*) 'DEFP(1) = ',DEFP(1) 144 : * WRITE(6,*) 'VARF(1) = ',VARF(1) 145 : * WRITE(6,*) 'SIGF(1) = ',SIGF(1) 146 : * WRITE(6,*) 'SIGPF = ',SIGPF 147 : 148 : RETURN 149 : END 150 : 151 :
© Cast3M 2003 - All rights reserved.
Disclaimer