Télécharger thsig1.eso

Retour à la liste

Numérotation des lignes :

  1. C THSIG1 SOURCE CHAT 05/01/13 03:38:00 5004
  2. SUBROUTINE THSIG1(SHP,DJAC,XSTRS,NBNN,LRE,REL,R)
  3. C=====================================================================
  4. C
  5. C CALCULE SOMME H,I * SIG I J *H ,J
  6. C MODE DE FOURIER 0
  7. C
  8. C ENTREE
  9. C SHP(6,NBNO) = H H,X H,Y H,Z
  10. C DJAC = JACOBIEN
  11. C XSTRS(NSTRS) = CONTRAINTES
  12. C NBNN = NOMBRE DE NOEUDS
  13. C LRE = TAILLE DE LA MATRICE DE RIGIDITE
  14. C REL(LRE,LRE) = MATRICE DE RIGIDITE (LRE=3*NBNN)
  15. C R = rayon
  16. C
  17. C SORTIES
  18. C REL(LRE,LRE)=MATRICE DE RIGIDITE
  19. C JEANVOINE AOUT 95
  20. C=======================================================================
  21.  
  22. IMPLICIT INTEGER(I-N)
  23. IMPLICIT REAL*8(A-H,O-Z)
  24. DIMENSION SHP(6,*),XSTRS(*),REL(LRE,*)
  25.  
  26. DO 100 IA=1,NBNN
  27. DO 200 IB=1,NBNN
  28. C
  29. CC1
  30. 1 = SHP(2,IA)*XSTRS(1)*SHP(2,IB) + SHP(3,IA)*XSTRS(2)*SHP(3,IB)
  31. 1 + SHP(2,IA)*XSTRS(4)*SHP(3,IB) + SHP(3,IA)*XSTRS(4)*SHP(2,IB)
  32. 1 + (1.D0/(R**2))*SHP(1,IA)*XSTRS(3)*SHP(1,IB)
  33. C
  34. CC2
  35. 1 = SHP(2,IA)*XSTRS(1)*SHP(2,IB) + SHP(3,IA)*XSTRS(2)*SHP(3,IB)
  36. 1 + SHP(2,IA)*XSTRS(4)*SHP(3,IB) + SHP(3,IA)*XSTRS(4)*SHP(2,IB)
  37. C
  38. CC3
  39. 1 = SHP(2,IA)*XSTRS(1)*SHP(2,IB) + SHP(3,IA)*XSTRS(2)*SHP(3,IB)
  40. 1 + SHP(2,IA)*XSTRS(4)*SHP(3,IB) + SHP(3,IA)*XSTRS(4)*SHP(2,IB)
  41. 1 + (1.D0/(R**2))*SHP(1,IA)*XSTRS(3)*SHP(1,IB)
  42. C
  43. REL(1+3*(IA - 1),1+3*(IB - 1)) = DJAC*CC1 +
  44. 1 REL(1+3*(IA - 1),1+3*(IB - 1))
  45. REL(2+3*(IA - 1),2+3*(IB - 1)) = DJAC*CC2 +
  46. 1 REL(2+3*(IA - 1),2+3*(IB - 1))
  47. REL(3+3*(IA - 1),3+3*(IB - 1)) = DJAC*CC3 +
  48. 1 REL(3+3*(IA - 1),3+3*(IB - 1))
  49. C
  50. 200 CONTINUE
  51. 100 CONTINUE
  52.  
  53. RETURN
  54. END
  55.  
  56.  
  57.  
  58.  

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