Télécharger thsigh.eso

Retour à la liste

Numérotation des lignes :

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

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