Télécharger thsig2.eso

Retour à la liste

Numérotation des lignes :

  1. C THSIG2 SOURCE CHAT 05/01/13 03:38:04 5004
  2. SUBROUTINE THSIG2(SHP,DJAC,XSTRS,NBNN,LRE,REL,N,R)
  3. C=====================================================================
  4. C
  5. C CALCULE SOMME H,I * SIG I J *H ,J
  6. C MODE DE FOURIER
  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 N = Mode de Fourier
  16. C R = rayon
  17. C
  18. C SORTIES
  19. C REL(LRE,LRE)=MATRICE DE RIGIDITE
  20. C JEANVOINE AOUT 95
  21. C=======================================================================
  22.  
  23. IMPLICIT INTEGER(I-N)
  24. IMPLICIT REAL*8(A-H,O-Z)
  25. DIMENSION SHP(6,*),XSTRS(*),REL(LRE,*)
  26.  
  27. DO 100 IA=1,NBNN
  28. DO 200 IB=1,NBNN
  29. C
  30. CC1
  31. 1 = SHP(2,IA)*XSTRS(1)*SHP(2,IB) + SHP(3,IA)*XSTRS(2)*SHP(3,IB)
  32. 1 + SHP(2,IA)*XSTRS(4)*SHP(3,IB) + SHP(3,IA)*XSTRS(4)*SHP(2,IB)
  33. 1 + ((1+(N*N))/(R**2))*SHP(1,IA)*XSTRS(3)*SHP(1,IB)
  34. C
  35. CC2
  36. 1 = SHP(2,IA)*XSTRS(1)*SHP(2,IB) + SHP(3,IA)*XSTRS(2)*SHP(3,IB)
  37. 1 + SHP(2,IA)*XSTRS(4)*SHP(3,IB) + SHP(3,IA)*XSTRS(4)*SHP(2,IB)
  38. 1 + ((N*N)/(R**2))*SHP(1,IA)*XSTRS(3)*SHP(1,IB)
  39. C
  40. CC3
  41. 1 = SHP(2,IA)*XSTRS(1)*SHP(2,IB) + SHP(3,IA)*XSTRS(2)*SHP(3,IB)
  42. 1 + SHP(2,IA)*XSTRS(4)*SHP(3,IB) + SHP(3,IA)*XSTRS(4)*SHP(2,IB)
  43. 1 + ((1+(N*N))/(R**2))*SHP(1,IA)*XSTRS(3)*SHP(1,IB)
  44. C
  45. CC4 = 2*N/(R*R)*SHP(1,IA)*XSTRS(3)*SHP(1,IB)
  46. C
  47. REL(1+3*(IA - 1),1+3*(IB - 1)) = DJAC*CC1 +
  48. 1 REL(1+3*(IA - 1),1+3*(IB - 1))
  49. REL(2+3*(IA - 1),2+3*(IB - 1)) = DJAC*CC2 +
  50. 1 REL(2+3*(IA - 1),2+3*(IB - 1))
  51. REL(3+3*(IA - 1),3+3*(IB - 1)) = DJAC*CC3 +
  52. 1 REL(3+3*(IA - 1),3+3*(IB - 1))
  53. REL(3+3*(IA - 1),1+3*(IB - 1)) = DJAC*CC4 +
  54. 1 REL(3+3*(IA - 1),1+3*(IB - 1))
  55. REL(1+3*(IA - 1),3+3*(IB - 1)) = DJAC*CC4 +
  56. 1 REL(1+3*(IA - 1),3+3*(IB - 1))
  57. C
  58. 200 CONTINUE
  59. 100 CONTINUE
  60.  
  61. RETURN
  62. END
  63.  
  64.  
  65.  

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