1 : C DOHTIM SOURCE CHAT 05/01/12 22:54:14 5004
2 : C DOHPTR SOURCE PLA2 90/07/11 21:19:35
3 : SUBROUTINE DOHTIM(VELA,MATE,IFOU,WORK,LHOOK,DDHOOK,IRET)
4 : C
5 : C=======================================================================
6 : C
7 : C MATRICE DE HOOK DES POUTRES DE TIMOSHENKO
8 : C
9 : C ENTREES
10 : C VELA() = materiau dans un tableau de travail
11 : C MATE = Nom du materiau
12 : C IFOU = num{ro d'harmonique de fourier: IFOUR de CCOPTIO
13 : C WORK = tabbleau contenant les caracteristiques de poutres
14 : C LHOOK = taille de la matrice de hooke
15 : C
16 : C SORTIES
17 : C DDHOOK(LHOOK,LHOOK) = matrice de hooke
18 : C IRET = 1 si option existante 0 SINON
19 : C
20 : C PP 4/9/92
21 : C=======================================================================
22 : IMPLICIT INTEGER(I-N)
23 : IMPLICIT REAL*8(A-H,O-Z)
24 : CHARACTER*8 MATE
25 : C
26 : DIMENSION VELA(*),DDHOOK(LHOOK,*),WORK(*)
27 : C
28 : IRET=1
28 : CALL ZERO(DDHOOK,LHOOK,LHOOK)
29 : C
30 : C ON RECUPERE LES CARACTERISTIQUES
31 : C
32 : XINX=WORK(1)
33 : XINY=WORK(2)
34 : XINZ=WORK(3)
35 : SD=WORK(4)
36 : SREDY=WORK(5)
37 : SREDZ=WORK(6)
38 : C
39 : C SECTION "REDUITE"
40 : C
41 : IF ( SREDZ. NE.0. ) THEN
42 : SRZ=SREDZ
43 : ELSE
44 : SRZ=SD
45 : ENDIF
46 : IF ( SREDY. NE.0. ) THEN
47 : SRY=SREDY
48 : ELSE
49 : SRY=SD
50 : ENDIF
51 : C
52 : C MATERIAU ISOTROPE
53 : C
54 : IF(MATE.EQ.'ISOTROPE') THEN
55 : C
56 : C LA MATRICE EST DIAGONALE
57 : C
58 : IF(IFOU.EQ.2) THEN
59 : YOU=VELA(1)
60 : XNU=VELA(2)
61 : CISA=YOU*0.5D0/(1.D0+XNU)
62 : C
63 : DDHOOK(1,1)=YOU*SD
64 : DDHOOK(2,2)=CISA*SRY
65 : DDHOOK(3,3)=CISA*SRZ
66 : DDHOOK(4,4)=CISA*XINX
67 : DDHOOK(5,5)=YOU*XINY
68 : DDHOOK(6,6)=YOU*XINZ
69 : C
70 : ELSE
71 : IRET=0
72 : ENDIF
73 : ELSE
74 : IRET=0
75 : ENDIF
76 : RETURN
77 : END
78 :
© Cast3M 2003 - All rights reserved.
Disclaimer