1 : C DOHLIS SOURCE CHAT 05/01/12 22:53:39 5004
2 : C
3 : SUBROUTINE DOHLIS(VELA,MATE,IFOU,LHOOK,DDHOOK,IRET)
4 : C
5 : C=======================================================================
6 : C
7 : C MATRICE DE HOOK DU LINESPRING
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 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 : C
22 : IMPLICIT INTEGER(I-N)
23 : IMPLICIT REAL*8(A-H,O-Z)
24 : PARAMETER(UNDEMI=.5D0,XZER=0.0D0,XK=1.2D0)
25 : PARAMETER(UN=1.D0,DEUX=2.D0)
26 : CHARACTER*8 MATE
27 : C
28 : DIMENSION VELA(*),DDHOOK(LHOOK,*)
29 : C
30 : IRET=1
30 : CALL ZERO(DDHOOK,LHOOK,LHOOK)
31 : C
32 : C MATERIAU ISOTROPE
33 : C
34 : IF(MATE.EQ.'ISOTROPE') THEN
35 : C
36 : C TRIDIMENSIONNEL
37 : C
38 : IF(IFOU.EQ.2) THEN
39 : YOU=VELA(1)
40 : XNU=VELA(2)
41 : EE=YOU*UNDEMI/(UN-XNU*XNU)
42 : C
43 : DDHOOK(1,1)=EE
44 : DDHOOK(2,2)=EE
45 : ELSE
46 : IRET=0
47 : ENDIF
48 : ELSE
49 : IRET=0
50 : ENDIF
51 : RETURN
52 : END
53 :
© Cast3M 2003 - All rights reserved.
Disclaimer