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