Esope source files

   1 : C NTNST     SOURCE    CHAT      05/01/13    02:02:34     5004           
   2 : 
   3 : C=======================================================================
   4 : C=                            N T N S T                                =
   5 : C=                            ---------                                =
   6 : C=                                                                     =
   7 : C=  Fonction :                                                         =
   8 : C=  ----------                                                         =
   9 : C=  Ajoute, a la matrice de rigidite globale, la contribution du point =
  10 : C=  de Gauss donnee par la matrice elementaire N et le jacobien.       =
  11 : C=  Contribution = transposee(N) * N * Jacobien                        =
  12 : C=                                                                     =
  13 : C=  Parametres :                                                       =
  14 : C=  ------------                                                       =
  15 : C=   XGENE(NST,LRE)  (E)   Matrice elementaire N au point de Gauss     =
  16 : C=   DJAC            (E)   Jacobien au point de Gauss considere        =
  17 : C=   LRE             (E)   Nombre de DDLs de la matrice de rigidite    =
  18 : C=   NST             (E)   Nombre de lignes de la matrice N (XGENE)    =
  19 : C=   REL            (E/S)  Matrice de rigidite GLOBALE                 =
  20 : C=                                                                     =
  21 : C=  Remarque : Si le parametre NST est superieur a 10, il faudra redi- =
  22 : C=  ---------- mensionner le tableau de travail Wrk !                  =
  23 : C=======================================================================
  24 : 
  25 :       SUBROUTINE NTNST (XGENE,DJAC,LRE,NST,REL)
  26 : 
  27 :       IMPLICIT INTEGER(I-N)
  28 :       IMPLICIT REAL*8 (A-H,O-Z)
  29 : 
  30 : -INC CCREEL
  31 : 
  32 :       DIMENSION XGENE(NST,*),REL(LRE,*)
  33 :       DIMENSION Wrk(10)
  34 : 
  35 : C= Boucles sur les DDLs (iX et iY)
  36 :       DO iX=1,LRE
  37 :         DO i=1,NST
  38 :           Wrk(i)=DJAC*XGENE(i,iX)
  39 :         ENDDO
  40 :         DO iY=1,iX
  41 : C== CC = Jacobien * transposee(N)(iX,i) * N(i,iY) (somme sur i)
  42 :           CC=XZero
  43 :           DO i=1,NST
  44 :             CC=CC+XGENE(i,iY)*Wrk(i)
  45 :           ENDDO
  46 : C== REL est une matrice symetrique
  47 :           REL(iX,iY)=REL(iX,iY)+CC
  48 :           REL(iY,iX)=REL(iX,iY)
  49 :         ENDDO
  50 :       ENDDO
  51 : 
  52 :       RETURN
  53 :       END
  54 : 
  55 :  
  56 :