1 : C DOHTIF SOURCE CHAT 05/01/12 22:54:10 5004
2 : SUBROUTINE DOHTIF(VELA,MATE,IFOU,LHOOK,DDHOOK,IRET)
3 : C
4 : C=======================================================================
5 : C
6 : C MATRICE DE HOOK DES POUTRES DE TIMOSHENKO
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 LHOOK = taille de la matrice de hooke
13 : C
14 : C SORTIES
15 : C DDHOOK(LHOOK,LHOOK) = matrice de hooke
16 : C IRET = 1 si option existante 0 SINON
17 : C
18 : C=======================================================================
19 : C Pierre Pegon (ISPRA) Juillet/Aout 1993
20 : C=======================================================================
21 : IMPLICIT INTEGER(I-N)
22 : IMPLICIT REAL*8(A-H,O-Z)
23 : CHARACTER*8 MATE
24 : C
25 : DIMENSION VELA(*),DDHOOK(LHOOK,*)
26 : C
27 : IRET=1
27 : CALL ZERO(DDHOOK,LHOOK,LHOOK)
28 : C
29 : C MATERIAU ISOTROPE
30 : C
31 : IF(MATE.EQ.'SECTION') THEN
32 : C
33 : IF(IFOU.EQ.2) THEN
34 : C
35 : C LA MATRICE EST PLEINE
36 : C
37 : DDHOOK(1,1)=VELA(1)
38 : DDHOOK(1,5)=VELA(3)
39 : DDHOOK(5,1)=VELA(3)
40 : DDHOOK(1,6)=-VELA(2)
41 : DDHOOK(6,1)=-VELA(2)
42 : C
43 : DDHOOK(5,5)=VELA(6)
44 : DDHOOK(5,6)=-VELA(5)
45 : DDHOOK(6,5)=-VELA(5)
46 : C
47 : DDHOOK(6,6)=VELA(4)
48 : C
49 : DDHOOK(2,2)=VELA(8)
50 : DDHOOK(2,4)=-VELA(10)
51 : DDHOOK(4,2)=-VELA(10)
52 : C
53 : DDHOOK(3,3)=VELA(7)
54 : DDHOOK(3,4)=VELA(9)
55 : DDHOOK(4,3)=VELA(9)
56 : C
57 : DDHOOK(4,4)=VELA(11)+VELA(12)
58 : C
59 : ELSEIF(IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.-3) THEN
60 : C
61 : C LA MATRICE EST PLEINE
62 : C
63 : DDHOOK(1,1)=VELA(1)
64 : DDHOOK(1,3)=-VELA(2)
65 : DDHOOK(3,1)=-VELA(2)
66 : C
67 : DDHOOK(3,3)=VELA(3)
68 : C
69 : DDHOOK(2,2)=VELA(4)
70 : C
71 : ELSE
72 : IRET=0
73 : ENDIF
74 : ELSE
75 : IRET=0
76 : ENDIF
77 : C
78 : RETURN
79 : END
80 :
81 :
© Cast3M 2003 - All rights reserved.
Disclaimer