C DOUOTA SOURCE CHAT 05/01/12 22:54:33 5004 SUBROUTINE DOUOTA(VELA,MATE,IFOU,MAPL, SIG,VAR, LHOOK,DDHOOK,IRET) C C======================================================================= C C MATRICE DE HOOKE TANGENTE DES ELEMENTS JOINTS C C ENTREES C VELA() = materiau dans un tableau de travail C MATE = Nom du materiau C IFOU = num{ro d'harmonique de fourier: IFOUR de CCOPTIO C MAPL = numero du materiau plastique C LHOOK = taille de la matrice de hooke C SIG() = etat de contrainte C VAR() = variables internes C C SORTIES C DDHOOK(LHOOK,LHOOK) = matrice de hooke C IRET = 1 si option existante 0 SINON C C Pierre Pegon (ISPRA) OCTOBER 93 C======================================================================= C IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) CHARACTER*8 MATE C DIMENSION VELA(*),DDHOOK(LHOOK,*),SIG(*),VAR(*) C IRET=1 CALL ZERO(DDHOOK,LHOOK,LHOOK) C C MATERIAU ISOTROPE (SEULEMENT) C IF(MATE.EQ.'ISOTROPE') THEN C C MODELE DE JOINT DILATANT EN CONTRAINTE OU DEF. PLANE C IF(MAPL.EQ.35.AND. & (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1)) THEN C DDHOOK(1,1)=VAR(7) DDHOOK(1,2)=VAR(8) DDHOOK(2,1)=VAR(9) DDHOOK(2,2)=VAR(10) C C MODELE DE JOINT_SOFT C ELSEIF(MAPL.EQ.56.AND. & (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1)) THEN C DDHOOK(1,1)=VAR(3) DDHOOK(1,2)=VAR(4) DDHOOK(2,1)=0.D0 DDHOOK(2,2)=VAR(5) C C MODELE D'ANCRAGE_ELIGEHAUSEN C ELSEIF(MAPL.EQ.91.AND. & (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1)) THEN C DDHOOK(1,1)=VAR(4) DDHOOK(1,2)=0.D0 DDHOOK(2,1)=0.D0 DDHOOK(2,2)=VAR(5) C C OPTION NON DEFINIE C ELSE IRET=0 ENDIF C C OPTION NON DEFINIE C ELSE IRET=0 ENDIF C RETURN END