1 : C DEFSIG SOURCE PV 09/03/12 21:19:11 6325 2 : SUBROUTINE DEFSIG(MFR,NDEF, 3 : . INPLAS,IND,WRK1,WRK5,WTRAV, 4 : . IVASTF,IVARIF,IVADEP,COB,XMOB, 5 : . IB,IGAU,CMATE,MATE,MELE,KERRE) 6 : * 7 : ************************************************** 8 : * ENTREES 9 : ************************************************** 10 : * 11 : * MFR : formulation de l'élément 12 : * NSTRS : nombre de composantes des contraintes 13 : * NVARI : nombre des variables internes 14 : * NDEF : nombre des déformations 15 : * INPLAS : numéro du matériau inélastique 16 : * IND : 17 : * SIGF : contraintes à la fin du pas (WRK1) 18 : * SIG0 : contraintes au début du pas (WRK1) 19 : * DSIGT : incrément de contrainte au cours du pas (WTRAV) 20 : * EPINF : déformations inélastiques à la fin du pas (WRK5) 21 : * EPIN0 : déformations inélastiques au début du pas (WRK5) 22 : * EPST0 : déformations totales au début du pas (WRK5) 23 : * VARF : variables internes à la fin du pas (WRK1) 24 : * IVASTF : pointeur sur un segment mptval de contraintes 25 : * IVARIF : pointeur sur un segment mptval de variables internes 26 : * IVADEP : pointeur sur un segment mptval de deformations inelastiques 27 : * COB : porosité 28 : * IB : numéro de l'élément 29 : * IGAU : numéro du point de Gauss 30 : * DEFP : incrément de déformations inélastiques (WRK1) 31 : * CMATE : nom du matériau 32 : * MATE : numéro du matériau élastique 33 : * MELE : numéro élément fini 34 : * 35 : ************************************************** 36 : * SORTIES 37 : ************************************************** 38 : * 39 : * KERRE : indice d'erreur 40 : * SEGMENT COMPRENANT : 41 : * - les contraintes finales 42 : * - les variables internes finales 43 : * - les incréments de déformations inélastiques 44 : * 45 : ************************************************** 46 : * 47 : IMPLICIT INTEGER(I-N) 48 : IMPLICIT REAL*8(A-H,O-Z) 49 : * 50 : -INC CCOPTIO 51 : -INC SMCHAML 52 : * 53 : SEGMENT MPTVAL 54 : INTEGER IPOS(NS) ,NSOF(NS) 55 : INTEGER IVAL(NCOSOU) 56 : CHARACTER*16 TYVAL(NCOSOU) 57 : ENDSEGMENT 58 : * 59 : SEGMENT WRK1 60 : REAL*8 DDHOOK(LHOOK,LHOOK),SIG0(NSTRS),DEPST(NSTRS) 61 : REAL*8 SIGF(NSTRS),VAR0(NVARI),VARF(NVARI) 62 : REAL*8 DEFP(NSTRS),XCAR(ICARA) 63 : ENDSEGMENT 64 : * 65 : SEGMENT WRK5 66 : REAL*8 EPIN0(NSTRS),EPINF(NSTRS),EPST0(NSTRS) 67 : ENDSEGMENT 68 : * 69 : SEGMENT WTRAV 70 : REAL*8 DDAUX(LHOOK,LHOOK),VALMAT(NUMAT) 71 : REAL*8 VALCAR(NUCAR),DSIGT(NSTRS) 72 : REAL*8 TXR(IDIM,IDIM),DDHOMU(LHOOK,LHOOK) 73 : REAL*8 XLOC(IDIM,IDIM),XGLOB(IDIM,IDIM) 74 : REAL*8 D1HOOK(LHOOK,LHOOK),ROTHOO(LHOOK,LHOOK) 75 : ENDSEGMENT 76 : * 77 : CHARACTER*8 CMATE 78 : * 79 : * rearrangement pour milieu poreux 80 : * 81 : NSTRS=SIGF(/1) 82 : NVARI=VARF(/1) 83 : KERRE=0 84 : IF(MFR.EQ.33) THEN 85 : * 86 : *ZZZZZ CAS DES JOINTS NON LINEAIRES : A VOIR ! 87 : * 88 : * calcul des contraintes totales 89 : * 90 : SIGF(NSTRS) = SIG0(NSTRS)+DEPST(NSTRS)/XMOB 91 : DO 1994 IC=1,3 92 : IF(IFOUR.EQ.-2.AND.IC.EQ.3) GO TO 1994 93 : SIGF(IC) = SIGF(IC) 94 : & -COB*(EPST0(NSTRS)+DEPST(NSTRS)) 95 : * 96 : * ce qui suit a ete pompe plus bas 97 : * 98 : IF(IND.EQ.1.OR.INPLAS.EQ.31.OR.INPLAS.EQ.30 99 : & .OR.INPLAS.EQ.37.OR.INPLAS.EQ.66 100 : & .OR.INPLAS.EQ.141)THEN 101 : SIGF(NSTRS) = SIGF(NSTRS) 102 : & +COB*(DEPST(IC)- (EPINF(IC)-EPIN0(IC))) 103 : ELSE 104 : SIGF(NSTRS) = SIGF(NSTRS) 105 : & +COB*(DEPST(IC)- DEFP(IC)) 106 : ENDIF 107 : 1994 CONTINUE 108 : ENDIF 109 : * 110 : MPTVAL=IVASTF 111 : DO 1116 IC=1,NSTRS 112 : MELVAL=IVAL(IC) 113 : VELCHE(IGAU,IB)=SIGF(IC) 114 : 1116 continue 115 : c 116 : c et les variables internes finales 117 : c 118 : MPTVAL=IVARIF 119 : IF(MFR.EQ.7.AND.CMATE.EQ.'SECTION')THEN 120 : DO IC=1,NVARI 121 : MELVAL=IVAL(IC) 122 : IELCHE(IGAU,IB)=NINT(VARF(IC)) 123 : END DO 124 : * 125 : ELSE IF(INPLAS.EQ.74) THEN 126 : * 127 : * Modele MAXWELL - composante EPSE 128 : * 129 : MELVAL=IVAL(1) 130 : VELCHE(IGAU,IB)=VARF(1) 131 : * 132 : * les autres composantes 133 : * 134 : DO IC=2,NVARI 135 : MELVAL=IVAL(IC) 136 : IELCHE(IGAU,IB)=NINT(VARF(IC)) 137 : END DO 138 : * 139 : ELSE 140 : DO 1117 IC=1,NVARI 141 : MELVAL=IVAL(IC) 142 : VELCHE(IGAU,IB)=VARF(IC) 143 : 1117 continue 144 : ENDIF 145 : c 146 : c et les increments de deformations plastiques 147 : c 148 : MPTVAL=IVADEP 149 : DO 1118 IC=1,NDEF 150 : MELVAL=IVAL(IC) 151 : IF(IND.EQ.1.OR.INPLAS.EQ.31.OR.INPLAS.EQ.30 152 : & .OR.INPLAS.EQ.37.OR.INPLAS.EQ.66.OR.INPLAS.EQ.118 153 : & .OR. INPLAS.EQ.141)THEN 154 : VELCHE(IGAU,IB)=EPINF(IC)-EPIN0(IC) 155 : ELSE 156 : VELCHE(IGAU,IB)=DEFP(IC) 157 : ENDIF 158 : 1118 continue 159 : c 160 : 1000 RETURN 161 : END 162 : 163 : 164 : 165 : 166 : 167 : 168 : 169 : 170 : 171 :
© Cast3M 2003 - All rights reserved.
Disclaimer