1 : C DOHMUR SOURCE CHAT 05/01/12 22:53:48 5004
2 : SUBROUTINE DOHMUR(VELA,MATE,IFOU,WORK,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 WORK = tabbleau contenant les caracteristiques du mur global
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 PP 27/04/04
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 : SECT=WORK(1)
32 : XINZ=WORK(2)
33 : C
34 : C MATERIAU ISOTROPE
35 : C
36 : IF(MATE.EQ.'ISOTROPE') THEN
37 : C
38 : C LA MATRICE EST DIAGONALE
39 : C
40 : IF (IFOU.EQ.-1.OR.IFOU.EQ.-2) THEN
41 : YOU=VELA(1)
42 : XNU=VELA(2)
43 : CISA=YOU*0.5D0/(1.D0+XNU)
44 : C
45 : DDHOOK(1,1)=YOU *SECT
46 : DDHOOK(2,2)=YOU *XINZ
47 : DDHOOK(3,3)=CISA*SECT
48 : DDHOOK(4,4)=YOU *SECT
49 : DDHOOK(5,5)=YOU *XINZ
50 : C
51 : ELSE
52 : IRET=0
53 : ENDIF
54 : ELSE
55 : IRET=0
56 : ENDIF
57 : C
58 : RETURN
59 : END
60 :
61 :
© Cast3M 2003 - All rights reserved.
Disclaimer