Télécharger thsig3.eso

Retour à la liste

Numérotation des lignes :

  1. C THSIG3 SOURCE CHAT 05/01/13 03:38:08 5004
  2. SUBROUTINE THSIG3(SHP,DJAC,XSTRS,NBNN,LRE,REL,R)
  3. C=====================================================================
  4. C
  5. C CALCULE SOMME H,I * SIG I J *H ,J
  6. C AXISYMETRIQUE
  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=2*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. C
  42. REL(1+2*(IA - 1),1+2*(IB - 1)) = DJAC*CC1 +
  43. 1 REL(1+2*(IA - 1),1+2*(IB - 1))
  44. REL(2+2*(IA - 1),2+2*(IB - 1)) = DJAC*CC2 +
  45. 1 REL(2+2*(IA - 1),2+2*(IB - 1))
  46. C
  47. 200 CONTINUE
  48. 100 CONTINUE
  49.  
  50. RETURN
  51. END
  52.  
  53.  
  54.  
  55.  

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