1 : C DOHPT2 SOURCE CHAT 05/01/12 22:54:01 5004
2 : C DOHPT2 SOURCE PLA2 90/07/11 21:19:35
3 : SUBROUTINE DOHPT2(VELA,MATE,IFOU,WORK,LHOOK,DDHOOK,IRET)
4 : C
5 : C=======================================================================
6 : C
7 : C MATRICE DE HOOK DES POUTRES 2D
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 Passage en nouveau CHAMELEM par I.Monnier le 4 05 90
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 : XINZ=WORK(2)
33 : SD=WORK(1)
34 : SREDY=WORK(3)
35 : C
36 : C MATERIAU ISOTROPE
37 : C
38 : IF(MATE.EQ.'ISOTROPE') THEN
39 : C
40 : C LA MATRICE EST DIAGONALE
41 : C
42 : IF(IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.-3) THEN
43 : YOU=VELA(1)
44 : XNU=VELA(2)
45 : CISA=YOU*0.5D0/(1.D0+XNU)
46 : C
47 : DDHOOK(1,1)=YOU*SD
48 : DDHOOK(2,2)=CISA*SREDY
49 : DDHOOK(3,3)=YOU*XINZ
50 : C
51 : ELSE
52 : IRET=0
53 : ENDIF
54 : ELSE
55 : IRET=0
56 : ENDIF
57 : RETURN
58 : END
59 :
60 :
© Cast3M 2003 - All rights reserved.
Disclaimer