1 : C RCDST SOURCE CHAT 05/01/13 02:46:04 5004
2 : SUBROUTINE RCDST(XE,NSTRS,LRE,DDHOMU,HS4,HS5,HS6,
3 : 1 REL,BGENE,IOP)
4 : *
5 : * CALCUL LA MATRICE B ET LA MATRICE DE RIGIDITE
6 : * RELATIVE AUX CISAILLEMENTS TRANSVERSES
7 : *
8 : ************************************************************************
9 : *
10 : * ENTREES :
11 : *
12 : * XE : COORDONNEES DES NOEUDS
13 : * EPAIST :EPAISSEUR
14 : * NSTRS :NOMBRE DE COMPOSANTES DE CONTRAINTES
15 : * LRE :NOMBRE DE DEGRES DE LIBERTE
16 : * DDHOMU :MATRICE DE HOOKE
17 : * OP := 0 LA SORTIE EST LA MATRICE DE RIGIDITE
18 : * = 1 LA SORTIE EST LA MATRICE B
19 : *
20 : * SORTIES :
21 : *
22 : * HS4,HS5,HS6 : LES FONCTIONS QUI RELIE LES TERMES (BETA)SK A {U}
23 : * REL :MATRICE DE RIGIDITE
24 : * BGENE :MATRICE DE B
25 : *
26 : * AUTEUR :P.DOWLATYARI JUILLET 91
27 : ************************************************************************
28 : IMPLICIT INTEGER(I-N)
29 : IMPLICIT REAL*8(A-H,O-Z)
30 : DIMENSION XE(3,*),DDHOMU(NSTRS,*),HS4(*),HS5(*)
31 : DIMENSION HS6(*),REL(LRE,*),BGENE(NSTRS,*)
32 : DIMENSION SX(3),SY(3),RL(3),EXX(3),EYY(3)
33 : DIMENSION DCIS(2,2),EKCIS(9,9),BCISA(2,9),HFX(9),HFY(9)
34 : DIMENSION ASS(6),BSS(6),AXX(3),BXX(3),AYY(3),BYY(3)
35 : *
36 : *
37 : DO 10 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 : 10 CONTINUE
57 : AIR=ABS(0.5D0*(SX(1)*SY(2)-SX(2)*SY(1)))
58 : *
58 : CALL CISAF(XE,EXX,EYY,DDHOMU,AIR,ASS,AXX,
59 : 1 AYY,BSS,BXX,BYY)
60 : *
60 : CALL CISAR(ASS,BSS,AXX,AYY,BXX,BYY,RL,EXX,EYY,
61 : 1 HS4,HS5,HS6,HFX,HFY)
62 : *
63 : DO 20 IC=1,9
64 : BCISA(1,IC)=HFX(IC)
65 : BCISA(2,IC)=HFY(IC)
66 : 20 CONTINUE
67 : *
68 : IF(IOP.EQ.0)THEN
69 : DCIS(1,1)=DDHOMU(7,7)
70 : DCIS(1,2)=DDHOMU(7,8)
71 : DCIS(2,1)=DDHOMU(8,7)
72 : DCIS(2,2)=DDHOMU(8,8)
73 : *
73 : CALL ZERO(EKCIS,9,9)
73 : CALL BDBST(BCISA,AIR,DCIS,9,2,EKCIS)
74 : *
75 : * REMPLISSAGE DANS REL
76 : *
77 : DO 30 IL=1,9
78 : IL1=(IL-1)/3
79 : IL2=IL+3*IL1+2
80 : DO 30 IC=1,9
81 : IC1=(IC-1)/3
82 : IC2=IC+IC1*3+2
83 : REL(IL2,IC2)=EKCIS(IL,IC)
84 : 30 CONTINUE
85 : ELSE
86 : *
87 : * REMPLISSAGE DANS BGENE
88 : *
89 : DO 40 IL=1,2
90 : IL1=IL+6
91 : DO 40 IC=1,9
92 : IC1=(IC-1)/3
93 : IC2=IC1*3+IC+2
94 : BGENE(IL1,IC2)=BCISA(IL,IC)
95 : 40 CONTINUE
96 : ENDIF
97 : *
98 : RETURN
99 : END
100 :
© Cast3M 2003 - All rights reserved.
Disclaimer