1 : C DSTGR1 SOURCE CHAT 05/01/12 22:56:39 5004
2 : SUBROUTINE DSTGR1(IGAU,NBNO,XE,NGRA,QSI,ETA,SHPTOT,
3 : 1 SHP,EXCEN,HS4,HS5,HS6,BGR,DJAC)
4 : *
5 : * CALCUL LA MATRICE BGR RELATIVE AUX EFFETS DE
6 : * MEMBRANE ET DE FLEXION
7 : *
8 : *
9 : * ENTREES :
10 : *
11 : * IGAU :NUMERO DU POINT DE GAUSS
12 : * XE :COORDONNEES DES NOEUDS
13 : * EXCEN :EXCENTEMENT DE COQUE
14 : * NBNO :NOMBRE DE NOEUDS
15 : * NGRA :NOMBRE DE COMPOSANTES DE GRADIENTS
16 : * QSI :COORDONNEES QSI DES POINTS DE GAUSS
17 : * ETA :COORDONNEES ETA DES POINTS DE GAUSS
18 : * SHPTOT :FONCTIONS DE FORME ET LEURS DERIVEE ( EFFET DE MEMBRANE)
19 : * SHP :TABLEAU DE TRAVAIL
20 : * HS4,HS5:FONCTIONS QUI RELIENT LES ROTATIONS (BETA)SK A {U}
21 : * HS6
22 : *
23 : * SORTIES :
24 : *
25 : * BGR :MATRICE BGR
26 : * DJAC :JACOBIEN
27 : *
28 : *
29 : * AUTEUR : X.Z. SUO FEVRIER 1996
30 : *
31 : ************************************************************************
32 : IMPLICIT INTEGER(I-N)
33 : IMPLICIT REAL*8 (A-H,O-Z)
34 : DIMENSION XE(3,*),QSI(*),ETA(*),SHPTOT(6,3,*),SHP(6,*)
35 : DIMENSION HS4(*),HS5(*),HS6(*),BGR(NGRA,*)
36 : DIMENSION SX(3),SY(3),RL(3),EXX(3),EYY(3)
37 : DIMENSION HXABS(9),HXORD(9),HYABS(9)
38 : DIMENSION HYORD(9),B(4,9)
39 : *
40 : * MATRICE BGR RELATIVE A L'EFFET DE MEMBRANE
41 : *
42 : DO 10 NPOI=1,3
43 : SHP(1,NPOI)=SHPTOT(1,NPOI,IGAU)
44 : SHP(2,NPOI)=SHPTOT(2,NPOI,IGAU)
45 : SHP(3,NPOI)=SHPTOT(3,NPOI,IGAU)
46 : 10 CONTINUE
46 : CALL JACOBI(XE,SHP,2,3,DJAC)
47 : K=1
48 : DO 30 NPOI=1,3
49 : BGR(1,K )=SHP(2,NPOI)
50 : BGR(5,K+1)=SHP(3,NPOI)
51 : BGR(2,K )=SHP(3,NPOI)
52 : BGR(4,K+1)=SHP(2,NPOI)
53 : K=K+6
54 : 30 CONTINUE
55 : *
56 : * MODIFICATION SI EXCENTREMENT DE COQUE
57 : *
58 : DO 40 K=4,6
59 : IF(K.EQ.4)THEN
60 : IJ=1
61 : I=2
62 : J=3
63 : ELSEIF(K.EQ.5)THEN
64 : IJ=2
65 : I=3
66 : J=1
67 : ELSE
68 : IJ=3
69 : I=1
70 : J=2
71 : ENDIF
72 : SX(IJ)=XE(1,I)-XE(1,J)
73 : SY(IJ)=XE(2,I)-XE(2,J)
74 : RL(IJ)=SQRT(SX(IJ)*SX(IJ)+SY(IJ)*SY(IJ))
75 : EXX(IJ)=-SX(IJ)/RL(IJ)
76 : EYY(IJ)=-SY(IJ)/RL(IJ)
77 : 40 CONTINUE
78 : AIR=ABS(0.5D0*(SX(1)*SY(2)-SX(2)*SY(1)))
79 : *
79 : CALL DERIVO(SX,SY,RL,QSI(IGAU),ETA(IGAU),HS4,HS5,
80 : 1 HS6,EXX,EYY,HXABS,HXORD,HYABS,HYORD)
80 : CALL BGRMAT(SX,SY,HXABS,HYABS,HXORD,HYORD,B)
81 : *
82 : K=2
83 : KK=0
84 : DO 1027 NPOI=1,NBNO
85 : DO 1127 IX =1,3
86 : BGR(1,K+IX)=EXCEN*B(1,IX+KK)
87 : BGR(2,K+IX)=EXCEN*B(2,IX+KK)
88 : BGR(4,K+IX)=EXCEN*B(3,IX+KK)
89 : BGR(5,K+IX)=EXCEN*B(4,IX+KK)
90 : 1127 CONTINUE
91 : K=K+6
92 : KK=KK+3
93 : 1027 CONTINUE
94 : *
95 : RETURN
96 : END
97 :
98 :
99 :
© Cast3M 2003 - All rights reserved.
Disclaimer