Télécharger ntnst.eso

Retour à la liste

Numérotation des lignes :

  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