Télécharger transk.eso

Retour à la liste

Numérotation des lignes :

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

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