1 : C DOUOTA SOURCE CHAT 05/01/12 22:54:33 5004
2 : SUBROUTINE DOUOTA(VELA,MATE,IFOU,MAPL, SIG,VAR, LHOOK,DDHOOK,IRET)
3 : C
4 : C=======================================================================
5 : C
6 : C MATRICE DE HOOKE TANGENTE DES ELEMENTS JOINTS
7 : C
8 : C ENTREES
9 : C VELA() = materiau dans un tableau de travail
10 : C MATE = Nom du materiau
11 : C IFOU = num{ro d'harmonique de fourier: IFOUR de CCOPTIO
12 : C MAPL = numero du materiau plastique
13 : C LHOOK = taille de la matrice de hooke
14 : C SIG() = etat de contrainte
15 : C VAR() = variables internes
16 : C
17 : C SORTIES
18 : C DDHOOK(LHOOK,LHOOK) = matrice de hooke
19 : C IRET = 1 si option existante 0 SINON
20 : C
21 : C Pierre Pegon (ISPRA) OCTOBER 93
22 : C=======================================================================
23 : C
24 : IMPLICIT INTEGER(I-N)
25 : IMPLICIT REAL*8(A-H,O-Z)
26 : CHARACTER*8 MATE
27 : C
28 : DIMENSION VELA(*),DDHOOK(LHOOK,*),SIG(*),VAR(*)
29 : C
30 : IRET=1
30 : CALL ZERO(DDHOOK,LHOOK,LHOOK)
31 : C
32 : C MATERIAU ISOTROPE (SEULEMENT)
33 : C
34 : IF(MATE.EQ.'ISOTROPE') THEN
35 : C
36 : C MODELE DE JOINT DILATANT EN CONTRAINTE OU DEF. PLANE
37 : C
38 : IF(MAPL.EQ.35.AND.
39 : & (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1)) THEN
40 : C
41 : DDHOOK(1,1)=VAR(7)
42 : DDHOOK(1,2)=VAR(8)
43 : DDHOOK(2,1)=VAR(9)
44 : DDHOOK(2,2)=VAR(10)
45 : C
46 : C MODELE DE JOINT_SOFT
47 : C
48 : ELSEIF(MAPL.EQ.56.AND.
49 : & (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1)) THEN
50 : C
51 : DDHOOK(1,1)=VAR(3)
52 : DDHOOK(1,2)=VAR(4)
53 : DDHOOK(2,1)=0.D0
54 : DDHOOK(2,2)=VAR(5)
55 : C
56 : C MODELE D'ANCRAGE_ELIGEHAUSEN
57 : C
58 : ELSEIF(MAPL.EQ.91.AND.
59 : & (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1)) THEN
60 : C
61 : DDHOOK(1,1)=VAR(4)
62 : DDHOOK(1,2)=0.D0
63 : DDHOOK(2,1)=0.D0
64 : DDHOOK(2,2)=VAR(5)
65 : C
66 : C OPTION NON DEFINIE
67 : C
68 : ELSE
69 : IRET=0
70 : ENDIF
71 : C
72 : C OPTION NON DEFINIE
73 : C
74 : ELSE
75 : IRET=0
76 : ENDIF
77 : C
78 : RETURN
79 : END
80 :
81 :
82 :
© Cast3M 2003 - All rights reserved.
Disclaimer