Télécharger la plaquette Cast3M


fr en

Télécharger ntnst.eso

Retour à la liste

Numérotation des lignes : 

ntnst

  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.  
  57.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales