1 : C BGFDST SOURCE CHAT 05/01/12 21:39:30 5004
2 : SUBROUTINE BGFDST(IGAU,XE,NGRAF,QSI,
3 : 1 ETA,HS4,HS5,HS6,BGF)
4 : *
5 : * CALCUL LA MATRICE BGF RELATIVE AUX
6 : * EFFETS DE FLEXION
7 : *
8 : *
9 : * ENTREES :
10 : *
11 : * IGAU :NUMERO DU POINT DE GAUSS
12 : * XE :COORDONNEES DES NOEUDS
13 : * NGRAF:NOMBRE DE COMPOSANTES DE GRADIENTS
14 : * QSI :COORDONNEES QSI DES POINTS DE GAUSS
15 : * ETA :COORDONNEES ETA DES POINTS DE GAUSS
16 : * HS4,HS5:FONCTIONS QUI RELIENT LES ROTATIONS (BETA)SK A {U}
17 : * HS6
18 : *
19 : * SORTIES :
20 : *
21 : * BGF :MATRICE BGF
22 : *
23 : *
24 : * AUTEUR : X.Z. SUO FEVRIER 1996
25 : *
26 : ************************************************************************
27 : IMPLICIT INTEGER(I-N)
28 : IMPLICIT REAL*8 (A-H,O-Z)
29 : DIMENSION XE(3,*),QSI(*),ETA(*)
30 : DIMENSION HS4(*),HS5(*),HS6(*),BGF(NGRAF,*)
31 : DIMENSION SX(3),SY(3),RL(3),EXX(3),EYY(3)
32 : DIMENSION HXABS(9),HXORD(9),HYABS(9)
33 : DIMENSION HYORD(9),B(4,9)
34 : *
35 : * MATRICE BGF RELATIVE A L'EFFET DE FLEXION
36 : *
37 : DO 40 K=4,6
38 : IF(K.EQ.4)THEN
39 : IJ=1
40 : I=2
41 : J=3
42 : ELSEIF(K.EQ.5)THEN
43 : IJ=2
44 : I=3
45 : J=1
46 : ELSE
47 : IJ=3
48 : I=1
49 : J=2
50 : ENDIF
51 : SX(IJ)=XE(1,I)-XE(1,J)
52 : SY(IJ)=XE(2,I)-XE(2,J)
53 : RL(IJ)=SQRT(SX(IJ)*SX(IJ)+SY(IJ)*SY(IJ))
54 : EXX(IJ)=-SX(IJ)/RL(IJ)
55 : EYY(IJ)=-SY(IJ)/RL(IJ)
56 : 40 CONTINUE
57 : AIR=ABS(0.5D0*(SX(1)*SY(2)-SX(2)*SY(1)))
58 : *
58 : CALL DERIVO(SX,SY,RL,QSI(IGAU),ETA(IGAU),HS4,HS5,
59 : 1 HS6,EXX,EYY,HXABS,HXORD,HYABS,HYORD)
59 : CALL BGRMAT(SX,SY,HXABS,HYABS,HXORD,HYORD,B)
60 : *
61 : K=2
62 : KK=0
63 : DO 327 NPOI=1,3
64 : DO 527 IY=1,3
65 : BGF(1,K+IY)=B(1,IY+KK)
66 : BGF(2,K+IY)=B(2,IY+KK)
67 : BGF(4,K+IY)=B(3,IY+KK)
68 : BGF(5,K+IY)=B(4,IY+KK)
69 : 527 CONTINUE
70 : KK=KK+3
71 : K=K+6
72 : 327 CONTINUE
73 : *
74 : RETURN
75 : END
76 :
77 :
78 :
© Cast3M 2003 - All rights reserved.
Disclaimer