Télécharger la plaquette Cast3M


fr en

Télécharger transk.eso

Retour à la liste

Numérotation des lignes : 

transk

  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/
  32.       IJ=1
  33.       JTROIS=0
  34.       DO 100 J=1,NB
  35.       ITROIS=0
  36.       DO 200 I=1,J
  37.       ICH=IMLT(IJ)
  38.       IF(ICOM.EQ.1) ICH=4
  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     SELECTION DU TYPE DE CHANGEMENT DE REPERE
  45. C     TRANSFORME A EN AA A L AIDE DE B
  46.       GOTO (310,320,330,340) ICH
  47.   310 CALL VMULT1(A,AA,B)
  48.       GOTO 350
  49.   320 CALL VMULT2(A,AA,B)
  50.       GOTO 350
  51.   330 CALL VMULT3(A,AA,B)
  52.       GOTO 350
  53.   340 CALL VMULT4(A,AA,B)
  54.   350 CONTINUE
  55. C     TRANSVASE LA MATRICE AA DANS LE BLOC 3 3 DE VKE
  56. C     ET SYMMETRISE LE BLOC TRANSPOSE DE VKE CORRESPONDANT
  57.       DO 400 IA=1,3
  58.       DO 400 IB=1,3
  59.       VKE(ITROIS+IA,JTROIS+IB)=AA(IA,IB)
  60.       VKE(JTROIS+IB,ITROIS+IA)=AA(IA,IB)
  61.   400 CONTINUE
  62.       IJ=IJ+1
  63.   200 ITROIS=ITROIS+3
  64.   100 JTROIS=JTROIS+3
  65.       RETURN
  66.       END
  67.  
  68.  

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