1 : C THSIGH SOURCE CHAT 05/01/13 03:38:11 5004 2 : SUBROUTINE THSIGH(SHP,DJAC,XSTRS,NBNN,IDIM,LRE,REL) 3 : C===================================================================== 4 : C 5 : C CALCULE SOMME H,I * SIG I J *H ,J 6 : C 7 : C ENTREE 8 : C SHP(6,NBNO) = H H,X H,Y H,Z 9 : C DJAC =JACOBIEN 10 : C XSTRS(NSTRS) = CONTRAINTES 11 : C NBNN = NOMBRE DE NOEUDS 12 : C IDIM = DIMENSION 13 : C LRE =TAILLE DE LA MATRICE DE RIGIDITE 14 : C REL(LRE,LRE) =MATRICE DE RIGIDITE (LRE=IDIM*NBNN 15 : C 16 : C SORTIES 17 : C REL(LRE,LRE)=MATRICE DE RIGIDITE 18 : C EBERSOLT OCT 85 19 : C======================================================================= 20 : IMPLICIT INTEGER(I-N) 21 : IMPLICIT REAL*8(A-H,O-Z) 22 : DIMENSION SHP(6,*),XSTRS(*),REL(LRE,*) 23 : KA=0 24 : DO 100 IA=1,NBNN 25 : KB=0 26 : DO 200 IB=1,IA 27 : CC=SHP(2,IA)*XSTRS(1)*SHP(2,IB)+SHP(3,IA)*XSTRS(2)*SHP(3,IB) 28 : 1 +(SHP(2,IA)*SHP(3,IB)+SHP(3,IA)*SHP(2,IB))*XSTRS(4) 29 : IF (IDIM.EQ.2) GOTO 366 30 : CC=CC+SHP(4,IA)*XSTRS(3)*SHP(4,IB) 31 : 1 +(SHP(2,IA)*SHP(4,IB)+SHP(4,IA)*SHP(2,IB))*XSTRS(5) 32 : 1 +(SHP(3,IA)*SHP(4,IB)+SHP(4,IA)*SHP(3,IB))*XSTRS(6) 33 : 366 CONTINUE 34 : CC=CC*DJAC 35 : DO 300 IC=1,IDIM 36 : REL(KA+IC,KB+IC)=REL(KA+IC,KB+IC)+CC 37 : REL(KB+IC,KA+IC)=REL(KA+IC,KB+IC) 38 : 300 CONTINUE 39 : 200 KB=KB+IDIM 40 : 100 KA=KA+IDIM 41 : RETURN 42 : END 43 :
© Cast3M 2003 - All rights reserved.
Disclaimer