Télécharger la plaquette Cast3M


fr en

Télécharger thsig2.eso

Retour à la liste

Numérotation des lignes : 

thsig2

  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