1 : C DSTGR2 SOURCE CHAT 05/01/12 22:56:47 5004
2 : SUBROUTINE DSTGR2(XE,NGRA,NSTRS,DDHOMU,HS4,HS5,HS6,BGR)
3 : *
4 : * CALCUL LA MATRICE BGR RELATIVE
5 : * AUX CISAILLEMENTS TRANSVERSES
6 : *
7 : ************************************************************************
8 : *
9 : * ENTREES :
10 : *
11 : * XE : COORDONNEES DES NOEUDS
12 : * NSTRS :NOMBRE DE COMPOSANTES DE CONTRAINTES
13 : * NGRA :NOMBRE DE COMPOSANTES DE GRADIENT
14 : * DDHOMU :MATRICE DE HOOKE
15 : *
16 : * SORTIES :
17 : *
18 : * HS4,HS5,HS6 : LES FONCTIONS QUI RELIE LES TERMES (BETA)SK A {U}
19 : * BGR : MATRICE DE BGR
20 : *
21 : * AUTEUR : X.Z. SUO FEVRIER 96
22 : ************************************************************************
23 : IMPLICIT INTEGER(I-N)
24 : IMPLICIT REAL*8(A-H,O-Z)
25 : DIMENSION XE(3,*),DDHOMU(NSTRS,*),HS4(*),HS5(*),HS6(*)
26 : DIMENSION BGR(NGRA,*),SX(3),SY(3),RL(3),EXX(3),EYY(3)
27 : DIMENSION BCISA(2,9),HFX(9),HFY(9)
28 : DIMENSION ASS(6),BSS(6),AXX(3),BXX(3),AYY(3),BYY(3)
29 : *
30 : DO 10 K=4,6
31 : IF(K.EQ.4)THEN
32 : IJ=1
33 : I=2
34 : J=3
35 : ELSEIF(K.EQ.5)THEN
36 : IJ=2
37 : I=3
38 : J=1
39 : ELSE
40 : IJ=3
41 : I=1
42 : J=2
43 : ENDIF
44 : SX(IJ)=XE(1,I)-XE(1,J)
45 : SY(IJ)=XE(2,I)-XE(2,J)
46 : RL(IJ)=SQRT(SX(IJ)*SX(IJ)+SY(IJ)*SY(IJ))
47 : EXX(IJ)=-SX(IJ)/RL(IJ)
48 : EYY(IJ)=-SY(IJ)/RL(IJ)
49 : 10 CONTINUE
50 : AIR=ABS(0.5D0*(SX(1)*SY(2)-SX(2)*SY(1)))
51 : SOMM=SX(IJ)+SY(IJ)+RL(IJ)+EXX(IJ)+EYY(IJ)
52 : *
52 : CALL CISAF(XE,EXX,EYY,DDHOMU,AIR,ASS,AXX,
53 : 1 AYY,BSS,BXX,BYY)
54 : *
54 : CALL CISAR(ASS,BSS,AXX,AYY,BXX,BYY,RL,EXX,EYY,
55 : 1 HS4,HS5,HS6,HFX,HFY)
56 : *
57 : DO 20 IC=1,9
58 : BCISA(1,IC)=HFX(IC)
59 : BCISA(2,IC)=HFY(IC)
60 : 20 CONTINUE
61 : *
62 : * REMPLISSAGE DANS BGR
63 : *
64 : DO 40 IL=1,2
65 : IL1=IL+6
66 : IL2=IL*3
67 : DO 40 IA=1,3
68 : IC1=3*(2*IA-1)
69 : IC2=3*IA-2
70 : BGR(IL1,IC1)=BCISA(IL,IC2)
71 : BGR(IL2,IC1+1)=BCISA(IL,IC2+1)
72 : BGR(IL2,IC1+2)=BCISA(IL,IC2+2)
73 : 40 CONTINUE
74 : *
75 : RETURN
76 : END
77 :
78 :
© Cast3M 2003 - All rights reserved.
Disclaimer