1 : C CQ4KSG SOURCE CHAT 05/01/12 22:27:01 5004
2 : SUBROUTINE CQ4KSG(DJAC,XSTRS,SHP,REL)
3 : C=====================================================================
4 : C MATICE KSIGMA POUR L ELEMENT COQ4
5 : C
6 : C ENTREES
7 : C DJAC = JACOBIEN
8 : C XSTRS(6) = CONTRAINTES AU POINT DE GAUSS
9 : * SHP = Fonctions de forme et d{riv{es dans l'espace
10 : * g{om{trique, au centre de l'{l{ment.
11 : * SHP(2,... = d{riv{e par rapport @ "X" local,
12 : * SHP(3,... = d{riv{e par rapport @ "Y" local.
13 : *
14 : C SORTIES
15 : C REL (24,24) = MATRICE K SIGMA
16 : C
17 : C DELHOMMEAU
18 : * P.M. 17/06/91: reprise pour nouveau COQ4, avec utilisation de
19 : * "SHP" au lieu de "BGENE".
20 : *
21 : C=======================================================================
22 : IMPLICIT INTEGER(I-N)
23 : IMPLICIT REAL*8(A-H,O-Z)
24 : DIMENSION REL(24,*),XSTRS(*),SHP(6,*)
25 : C
26 : C
27 : SIGXX=XSTRS(1)
28 : SIGYY=XSTRS(2)
29 : SIGXY=XSTRS(3)
30 : C
31 : DO 100 IA=1,4
32 : JJ=(IA-1)*6+1
33 : DO 200 IB=1,4
34 : II=(IB-1)*6+1
35 : CC= SHP(2,IA)*SHP(2,IB)*SIGXX + SHP(3,IA)*SHP(3,IB)*SIGYY
36 : DD= (SHP(2,IA)*SHP(3,IB) + SHP(2,IB)*SHP(3,IA)) * SIGXY
37 : S =(CC+DD)*DJAC*4.D0
38 : REL(II,JJ) = S
39 : REL(II+1,JJ+1)= S
40 : REL(II+2,JJ+2)= S
41 : 200 CONTINUE
42 : 100 CONTINUE
43 : C
44 : RETURN
45 : END
46 :
47 :
© Cast3M 2003 - All rights reserved.
Disclaimer