douota
C DOUOTA SOURCE CHAT 05/01/12 22:54:33 5004 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 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
© Cast3M 2003 - Tous droits réservés.
Mentions légales