Télécharger transg.eso

Retour à la liste

Numérotation des lignes :

  1. C TRANSG SOURCE CHAT 05/01/13 03:44:57 5004
  2. SUBROUTINE TRANSG(VKE,B,LRE,NBNO,ICOM)
  3. C=====================================================
  4. C
  5. C TRANSFORME VKE(LRE,LRE) DU REPERE LOCAL VKE ANTISYMETRIQUE
  6. C VKE = CORIOLIS
  7. C AU REPERE GENERAL A L AIDE DE LA MATRICE DE PASSAGE
  8. C B(3,3)
  9. C NBNO NOMBRE DE NOEUDS DE L ELEMENT 3 NOEUDS DKT
  10. C 4 NOEUDS DKQ
  11. C LRE NOMBRE DE DDL DE L ELEMENT
  12. C 6 DDL PAR NOEUD
  13. C ENTREES
  14. C VKE(LRE,LRE) = MATRICE DE RAIDEUR APRES CHGMNT DE REPERE
  15. C B(3,3) = MATRICE DE PASSAGE
  16. C LRE = NOMBRE DE D.D.L. DANS LA MATRICE DE RIGIDITE
  17. C NBNO = NOMBRE DE NOEUDS ( LRE = NBNO*6
  18. C ICOM = 0 SI LES D.D.L. DE FLEXION SONT DECOUPLES DES
  19. C D.D.L. DE MEMBRANE 1 SINON
  20. C SORTIES
  21. C VKE(LRE,LRE) = MATRICE DE CORIOLIS APRES CHGMNT DE REPERE
  22. C
  23. C D'APRES TRANSK COMBESCURE 2004
  24. C=====================================================
  25. IMPLICIT INTEGER(I-N)
  26. IMPLICIT REAL*8(A-H,O-Z)
  27. DIMENSION VKE(LRE,*),B(3,*)
  28. DIMENSION AA(3,3),A(3,3)
  29. DIMENSION IMLT(36)
  30. DATA IMLT/1,2,1,4,3,1,2,4,2,1,4,3,4,3,1,2,4,2,4,2,1,
  31. 1 4,3,4,3,4,3,1,2,4,2,4,2,4,2,1/
  32. IJ=1
  33. JTROIS=0
  34. DO 100 J=1,NB
  35. ITROIS=0
  36. DO 200 I=1,NB
  37. ITROIS =3*I - 3
  38. JTROIS =3*J - 3
  39. C TRANSVASE UN BLOC 3 3 DE VKE DANS LA MATRICE DE TRAVAIL A
  40. DO 300 IA=1,3
  41. DO 300 IB=1,3
  42. A(IA,IB)=VKE(ITROIS+IA,JTROIS+IB)
  43. 300 CONTINUE
  44. C MATRICE NON SYMETRIQUE
  45. CALL VMULT4(A,AA,B)
  46. 350 CONTINUE
  47. C TRANSVASE LA MATRICE AA DANS LE BLOC 3 3 DE VKE
  48. DO 400 IA=1,3
  49. DO 400 IB=1,3
  50. VKE(ITROIS+IA,JTROIS+IB)=AA(IA,IB)
  51. C VKE(JTROIS+IB,ITROIS+IA)=(-1.D0)*AA(IA,IB)
  52. 400 CONTINUE
  53. IJ=IJ+1
  54. 200 CONTINUE
  55. 100 CONTINUE
  56. C 200 ITROIS=ITROIS+3
  57. C 100 JTROIS=JTROIS+3
  58. RETURN
  59. END
  60.  
  61.  
  62.  

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