1 : C DOU266 SOURCE AM 12/02/27 23:39:37 7289
2 : SUBROUTINE DOU266(VELA,MATE,IFOU,LHOOK,DDHOOK,IRET)
3 : C
4 : C=======================================================================
5 : C
6 : C MATRICE DE HOOKE DES ELEMENTS DE ZONES COHESIVES ZCO2, ZCO3, ZCO4
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 LHOOK = taille de la matrice de hooke
13 : C
14 : C SORTIES
15 : C DDHOOK(LHOOK,LHOOK) = matrice de hooke
16 : C IRET = 1 si option existante 0 SINON
17 : C
18 : C=======================================================================
19 : C
20 : IMPLICIT INTEGER(I-N)
21 : IMPLICIT REAL*8(A-H,O-Z)
22 : CHARACTER*8 MATE
23 : C
24 : DIMENSION VELA(*),DDHOOK(LHOOK,*)
25 : DATA DEUX/2.0D0/
26 : C
27 : IRET=1
27 : CALL ZERO(DDHOOK,LHOOK,LHOOK)
28 : C
29 : C MATERIAU ISOTROPE
30 : C
31 : IF(MATE.EQ.'ZONE_COH') THEN
32 : C
33 : C CAS TRIDIMENSIONNEL
34 : C
35 : IF(IFOU.EQ.2) THEN
36 : ZK1 =VELA(1)
37 : ZK2 =VELA(1)
38 : ZK3 =VELA(2)
39 : C
40 : DDHOOK(1,1)=ZK1
41 : DDHOOK(2,2)=ZK2
42 : DDHOOK(3,3)=ZK3
43 : C
44 : C CAS CONTRAINTES PLANES DANS LE PLAN (S1,N)
45 : C DEFORMATIONS PLANES DANS LE PLAN (S1,N)
46 : C AXISYMETRIE DANS LE PLAN (S1,N)
47 : C
48 : ELSE IF(IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.
49 : & IFOU.EQ.-1.OR.IFOU.EQ.0) THEN
50 : ZK1 =VELA(1)
51 : ZK2 =VELA(2)
52 : C
53 : DDHOOK(1,1)=ZK1
54 : DDHOOK(2,2)=ZK2
55 : C
56 : C OPTION NON DEFINIE
57 : C
58 : ELSE
59 : IRET=0
60 : ENDIF
61 : C
62 : C AUTRES CAS : NON ENCORE PREVUS ...!
63 : C
64 : ELSE
65 : IRET=0
66 : ENDIF
67 : C
68 : RETURN
69 : END
70 :
71 :
72 :
73 :
© Cast3M 2003 - All rights reserved.
Disclaimer