1 : C DOHPTR SOURCE CHAT 05/01/12 22:54:04 5004
2 : SUBROUTINE DOHPTR(VELA,MATE,IFOU,WORK,LHOOK,DDHOOK,IRET)
3 : C
4 : C=======================================================================
5 : C
6 : C MATRICE DE HOOK DES POUTRES
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 WORK = tabbleau contenant les caracteristiques de poutres
13 : C LHOOK = taille de la matrice de hooke
14 : C
15 : C SORTIES
16 : C DDHOOK(LHOOK,LHOOK) = matrice de hooke
17 : C IRET = 1 si option existante 0 SINON
18 : C
19 : C Passage en nouveau CHAMELEM par I.Monnier le 4 05 90
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,*),WORK(*)
26 : C
27 : IRET=1
27 : CALL ZERO(DDHOOK,LHOOK,LHOOK)
28 : C
29 : C ON RECUPERE LES CARACTERISTIQUES
30 : C
31 : XINX=WORK(1)
32 : XINY=WORK(2)
33 : XINZ=WORK(3)
34 : SD=WORK(4)
35 : SREDY=WORK(5)
36 : SREDZ=WORK(6)
37 : C
38 : C MATERIAU ISOTROPE
39 : C
40 : IF(MATE.EQ.'ISOTROPE') THEN
41 : C
42 : C LA MATRICE EST DIAGONALE
43 : C
44 : IF(IFOU.EQ.2) THEN
45 : YOU=VELA(1)
46 : XNU=VELA(2)
47 : CISA=YOU*0.5D0/(1.D0+XNU)
48 : C
49 : DDHOOK(1,1)=YOU*SD
50 : DDHOOK(2,2)=CISA*SREDY
51 : DDHOOK(3,3)=CISA*SREDZ
52 : DDHOOK(4,4)=CISA*XINX
53 : DDHOOK(5,5)=YOU*XINY
54 : DDHOOK(6,6)=YOU*XINZ
55 : C
56 : ELSE
57 : IRET=0
58 : ENDIF
59 : ELSE
60 : IRET=0
61 : ENDIF
62 : RETURN
63 : END
64 :
© Cast3M 2003 - All rights reserved.
Disclaimer