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