1 : C ELAST3 SOURCE CHAT 05/01/12 23:29:35 5004 2 : SUBROUTINE ELAST3(ISENS,IFOUR,XMATT,MFR,NSTRS,DD,KERRE,INDIC) 3 : C 4 : C---------------------------------------------------------------------------- 5 : C Objet: Calcul de la matrice de Hooke endommagee(ISENS=1) 6 : C ou son inverse(ISENS=2) 7 : C dans le cadre du modele CERAMIQUE (materiau 65) 8 : C---------------------------------------------------------------------------- 9 : C 10 : C---------------------------------------------------------------------------- 11 : C Entree: ISENS= 1 DD en sortie est la matrice de Hooke endommagee 12 : C d'ou SIG=DD*EPSELAS 13 : C = 2 DD en sortie est la matrice de Hooke inverse 14 : C endommagee d'ou EPSELAS=DD*SIG 15 : C avec SIG contraintes et EPSELAS deformations elastiques 16 : C IFOUR= -2 EN CONTR. PLANES 17 : C -1 EN DEFORM. PLANES 18 : C 0 EN AXISYMETRIE 19 : C 1 EN SERIE DE FOURIER 20 : C 2 EN TRIDIM 21 : C NSTRS nombre de composantes des contraintes ou 22 : C des deformations 23 : C XMATT(NCOMAT) tableau des parametres scalaires du materiau 24 : C a une temperature T donnee 25 : C dont XMATT(1) module d'Young et XMATT(2) coeff de Poisson 26 : C MFR indice de la formulation mecanique (seulement massif 27 : C pour les materiaux endommageables) 28 : C INDIC = 0 ====> le point de Gauss n'est pas cassé 29 : C INDIC = 1 ====> le point de Gauss est cassé 30 : C-------------------------------------------------------------------------------- 31 : C 32 : C-------------------------------------------------------------------------------- 33 : C Sortie: DD(NSTRS,NSTRS) matrice de Hooke endommagee si ISENS=1 34 : C son inverse si ISENS=2 35 : C 36 : C KERRE indice qui regit les erreurs 37 : C Il vaut 99 si la formulation mecanique n'est pas disponible ou 38 : C s'il y a incompatibilite entre MFR et IFOUR 39 : C _______________________________________________________________________________ 40 : C 41 : IMPLICIT INTEGER(I-N) 42 : IMPLICIT REAL*8 (A-H,O-Z) 43 : DIMENSION XMATT(*) 44 : DIMENSION DD(NSTRS,NSTRS) 45 : C 45 : CALL ZERO(DD,NSTRS,NSTRS) 46 : KERRE=0 47 : C 48 : YUNG=XMATT(1) 49 : XNU=XMATT(2) 50 : C Modifs Cmau 51 : C D1=1.D0-XD 52 : D1=1.D0 53 : C 54 : C ==================================================================== 55 : C SI LE PT. DE GAUSS EST ROMPU, ON DIV. SON MOD. D'YOUNG PAR 1 MILLION 56 : C ==================================================================== 57 : IF (INDIC.EQ.1) THEN 58 : D1=1.D-6 59 : ENDIF 60 : C 61 : COEFF1=(D1*YUNG)/(1.D0-XNU*XNU) 62 : COEFF2=1.D0/(D1*YUNG) 63 : COEFF3=(D1*YUNG)/((1.D0+XNU)*(1.D0-2.D0*XNU)) 64 : COEFF4=(1.D0-2.D0*XNU)/2.D0 65 : C 66 : IF (MFR.EQ.1) THEN 67 : C 68 : C------------ 69 : C MASSIFS 70 : C------------ 71 : IF (IFOUR.NE.-2) THEN 72 : C 73 : C ---------------------------------- 74 : C DEFORMATIONS PLANES ET AXISYMETRIE 75 : C ---------------------------------- 76 : IF (ISENS.EQ.1) THEN 77 : C 78 : C---------------- 79 : C SIG=DD*EPSELAS 80 : C---------------- 81 : DD(1,1)=COEFF3*(1.D0-XNU) 82 : DD(1,2)=COEFF3*XNU 83 : DD(1,3)=DD(1,2) 84 : DD(2,2)=DD(1,1) 85 : DD(2,3)=DD(1,2) 86 : DD(3,3)=DD(1,1) 87 : DD(2,1)=DD(1,2) 88 : DD(3,1)=DD(1,2) 89 : DD(3,2)=DD(1,2) 90 : C 91 : LA=4 92 : LB=4 93 : C 94 : C --------------------------- 95 : C LB=6 POUR TRIDIM ET FOURIER 96 : C --------------------------- 97 : IF (IFOUR.GT.0) LB=6 98 : DO 5 IA=LA,LB 99 : DD(IA,IA)=COEFF3*COEFF4 100 : 5 CONTINUE 101 : ELSE 102 : C 103 : C------------------- 104 : C EPSELAS=DD*SIG 105 : C------------------- 106 : DD(1,1)=COEFF2 107 : DD(1,2)=-XNU*COEFF2 108 : DD(1,3)=DD(1,2) 109 : DD(2,2)=DD(1,1) 110 : DD(2,3)=DD(1,2) 111 : DD(3,3)=DD(1,1) 112 : DD(2,1)=DD(1,2) 113 : DD(3,1)=DD(1,3) 114 : DD(3,2)=DD(2,3) 115 : LA=4 116 : LB=4 117 : C 118 : C --------------------------- 119 : C LB=6 POUR TRIDIM ET FOURIER 120 : C --------------------------- 121 : IF (IFOUR.GT.0) LB=6 122 : DO 10 IA=LA,LB 123 : DD(IA,IA)=(1.D0+XNU)*DD(1,1)*2.D0 124 : 10 CONTINUE 125 : C 126 : ENDIF 127 : C 128 : C ------------- 129 : C CONTR. PLANES 130 : C ------------- 131 : ELSE 132 : IF (ISENS.EQ.1) THEN 133 : C 134 : C----------------- 135 : C SIG=DD*EPSELAS 136 : C----------------- 137 : DD(1,1)=COEFF1 138 : DD(1,2)=COEFF1*XNU 139 : DD(2,2)=COEFF1 140 : DD(4,4)=COEFF1*((1.D0-XNU)/2.D0) 141 : DD(2,1)=DD(1,2) 142 : ELSE 143 : C 144 : C------------------ 145 : C EPSELAS=DD*SIG 146 : C------------------ 147 : DD(1,1)=COEFF2 148 : DD(1,2)=-COEFF2*XNU 149 : DD(2,2)=COEFF2 150 : DD(4,4)=COEFF2*2.D0*(1.D0+XNU) 151 : DD(2,1)=DD(1,2) 152 : DD(3,1)=DD(1,2) 153 : DD(3,2)=DD(1,2) 154 : C 155 : ENDIF 156 : ENDIF 157 : C 158 : ELSE 159 : IF (MFR.EQ.5) THEN 160 : C 161 : C------------------------------------------ 162 : C COQUES EPAISSES (NON ENCORE DISPONIBLE) 163 : C------------------------------------------ 164 : C 165 : C FORMULATION NON DISPONIBLE 166 : ENDIF 167 : KERRE=99 168 : ENDIF 169 : C 170 : RETURN 171 : END 172 : 173 : 174 : 175 :
© Cast3M 2003 - All rights reserved.
Disclaimer