1 : C DOCO88 SOURCE CHAT 05/01/12 22:53:09 5004
2 : SUBROUTINE DOCO88(VELA,MATE,IFOU,LHOOK,DDHOOK,IRET)
3 : C
4 : C=======================================================================
5 : C
6 : C MATRICE DE HOOKE DES ELEMENTS JOINTS JCT3 ET JCI4 (2D)
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 REMARQUE : ATTENTION : DANS LES CAS CONTRAINTES PLANES,
19 : C DEFORMATIONS PLANES ET AXISYMETRIQUE,
20 : C LA MATRICE DE HOOK N'A PLUS LA DIMENSION
21 : C (3X3), MAIS (2X2). C'EST DONC LES DIMENSIONS DU CAS
22 : C BIDIMENSIONNEL.
23 : C
24 : C CECI EST UN DOUBLE DE DOHO88, MAIS AVEC LA THEORIE DE PARVIS
25 : C POUR LE CALCUL DE LA MATRICE DE TRANSFORMATION. LA MATRICE DE
26 : C ROTATION ROTHOO TRANSFORME UN VECTEUR (EN L'OCCURENCE LE VECTEUR
27 : C DEPLACEMENT) DU REPERE D'ORTHOTROPIE AU REPERE GLOBAL. IL FAUT
28 : C BIEN VOIR QUE CECI EST UNE TRANSFORMATION D'UN VECTEUR EN UN AUTRE
29 : C VECTEUR, ET NON PAS LA TRANSFORMATION D'UN TENSEUR (EX: TENSEUR
30 : C DES CONTRAINTES) EN UN AUTRE TENSEUR. C'EST DONC CETTE SUBROUTINE
31 : C QUI EST LA BONNE, CAR JE RAISONNE EN TERMES DE DEPLACEMENTS, ET
32 : C NON EN TERMES DE CONTRAINTES, DANS LE CAS DES JOINTS (cf. LA
33 : C RELATION FORCES/DEPLACEMENTS ET NON PAS CONTRAINTES/DEFORMATIONS)
34 : C
35 : C P. PEGON 2001 d'apres DOUO88
36 : C=======================================================================
37 : C
38 : IMPLICIT INTEGER(I-N)
39 : IMPLICIT REAL*8(A-H,O-Z)
40 : CHARACTER*8 MATE
41 : C
42 : DIMENSION VELA(*),DDHOOK(LHOOK,*)
43 : REAL*8 D3HOO1(3,3),D3HOO2(3,3),ROTHOO(3,3)
44 : DATA DEUX/2.0D0/
45 : C
46 : IRET=1
46 : CALL ZERO(DDHOOK,LHOOK,LHOOK)
46 : CALL ZERO(D3HOO1,LHOOK,LHOOK)
47 : C
48 : C MATERIAU ISOTROPE
49 : C
50 : IF(MATE.EQ.'ISOTROPE') THEN
51 : C
52 : C CAS CONTRAINTES PLANES DANS LE PLAN (S1,N)
53 : C DEFORMATIONS PLANES DANS LE PLAN (S1,N)
54 : C AXISYMETRIE DANS LE PLAN (S1,N)
55 : C
56 : IF(IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN
57 : C
58 : C et juste le cisaillement!
59 : C
60 : ZK1 =VELA(1)
61 : C
62 : DDHOOK(1,1)=ZK1
63 : DDHOOK(2,2)=ZK1
64 : C
65 : ELSE
66 : IRET=0
67 : ENDIF
68 : ELSE
69 : IRET=0
70 : ENDIF
71 : C
72 : RETURN
73 : END
74 :
75 :
76 :
77 :
© Cast3M 2003 - All rights reserved.
Disclaimer