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    CB215821  26/08/24    21:18:42     12622          
  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 401 IA=1,3
  41.       DO 300 IB=1,3
  42.       A(IA,IB)=VKE(ITROIS+IA,JTROIS+IB)
  43.   300 CONTINUE
  44.   401 CONTINUE
  45. C     SELECTION DU TYPE DE CHANGEMENT DE REPERE
  46. C     TRANSFORME A EN AA A L AIDE DE B
  47.       GOTO (310,320,330,340) ICH
  48.   310 CALL VMULT1(A,AA,B)
  49.       GOTO 350
  50.   320 CALL VMULT2(A,AA,B)
  51.       GOTO 350
  52.   330 CALL VMULT3(A,AA,B)
  53.       GOTO 350
  54.   340 CALL VMULT4(A,AA,B)
  55.   350 CONTINUE
  56. C     TRANSVASE LA MATRICE AA DANS LE BLOC 3 3 DE VKE
  57. C     ET SYMMETRISE LE BLOC TRANSPOSE DE VKE CORRESPONDANT
  58.       DO 402 IA=1,3
  59.       DO 400 IB=1,3
  60.       VKE(ITROIS+IA,JTROIS+IB)=AA(IA,IB)
  61.       VKE(JTROIS+IB,ITROIS+IA)=AA(IA,IB)
  62.   400 CONTINUE
  63.   402 CONTINUE
  64.       IJ=IJ+1
  65.       ITROIS=ITROIS+3
  66.   200 CONTINUE
  67.       JTROIS=JTROIS+3
  68.   100 CONTINUE
  69.       RETURN
  70.       END
  71.  
  72.  
  73.  

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