1 : C HOOINC SOURCE CHAT 05/01/13 00:25:03 5004 2 : SUBROUTINE HOOINC(D,YUNG,XNU,ITYP) 3 : C-------------------------------------------- 4 : C INSPIRE DE INCA 5 : C-------------------------------------------- 6 : IMPLICIT INTEGER(I-N) 7 : IMPLICIT REAL*8(A-H,O-Z) 8 : DIMENSION D(*) 9 : C*** 10 : C*** MATERIAU ISOTROPE 11 : C*** 12 : DO 1 IA=1,36 13 : 1 D(IA)=0.D0 14 : IF(ITYP.EQ.6) GO TO 3 15 : IF(ITYP.EQ.4) GO TO 4 16 : C CALCUL DE LAMDA ET MU 17 : A2MU=YUNG/(1.0D0+XNU) 18 : ALAMDA=A2MU*XNU/(1.0D0-2.0D0*XNU) 19 : AMU=0.5D0*A2MU 20 : C MISE ENPLACE DANS D DES LAMDA ET MU 21 : 202 CONTINUE 22 : JB=0 23 : JD=21 24 : DO 2 IA=1,3 25 : JA=JB 26 : DO 20 IB=1,3 27 : JA=JA+1 28 : 20 D(JA)=ALAMDA 29 : JC=JB+IA 30 : D(JC)=D(JC)+A2MU 31 : JE=JD+IA 32 : D(JE)=AMU 33 : JD=JD+6 34 : JB=JB+6 35 : 2 CONTINUE 36 : RETURN 37 : C CAS CONTRAINTES PLANES 38 : 3 CONTINUE 39 : ALAMDA=YUNG/(1.D0-XNU*XNU) 40 : AMU=0.5D0*(1.D0-XNU)*ALAMDA 41 : A2MU=XNU*ALAMDA 42 : D(1)=ALAMDA 43 : D(8)=ALAMDA 44 : D(22)=AMU 45 : D(2)=A2MU 46 : D(7)=A2MU 47 : RETURN 48 : C CAS UNIAXIAL 49 : 4 CONTINUE 50 : D(15)=YUNG 51 : RETURN 52 : END 53 :
© Cast3M 2003 - All rights reserved.
Disclaimer