Télécharger la plaquette Cast3M


fr en

Télécharger diaxyz.eso

Retour à la liste

Numérotation des lignes : 

diaxyz

  1. C DIAXYZ    SOURCE    CHAT      05/01/12    22:50:02     5004
  2.       SUBROUTINE DIAXYZ(ITYP,SIGT,SIGEL,VECPRO,IBO,A,B,C,B1,B2)
  3. C
  4. C         EFFECTUE LA TRANSFORMATION DE SIGAMA DIAGONALISE EN
  5. C         SIGAM GLOBAL
  6. C
  7. C         ENTREE  SIGT
  8. C         SORTIE SIGEL
  9. C             ITYP TYPE DE L'ELEMENT
  10. C
  11.       IMPLICIT INTEGER(I-N)
  12.       IMPLICIT REAL*8(A-H,O-Z)
  13.       DIMENSION SIGT(*),SIGEL(*),VECPRO(*)
  14. C
  15.       GO TO (70,71,71,73,74,75,76,77,78,79,80,81,71),ITYP
  16. C
  17. C   MODIF DU 05-01-96
  18. C   POUR ITYP=13 (COQUES EPAISSES) ON PASSAIT ICI (AU 83)
  19. C   DESORMAIS ON VA AU 71 COMME LES CALCULS EN CONTRAINTE PLANE
  20. C
  21. C   83 CONTINUE
  22. C               COQUES AVEC CISAILLEMENTS TRANSVERSES
  23. C      SIGEL(1)=B*(SIGT(1)-SIGT(2))
  24. C      SIGEL(2)=B*(SIGT(1)+SIGT(2))
  25. C      SIGEL(3)=SIGT(3)
  26. C      SIGEL(4)=SIGT(4)
  27. C      SIGEL(5)=SIGT(5)
  28. C      SIGEL(6)=SIGT(6)
  29. C      GO TO 740
  30.   81  CONTINUE
  31. C                      TUYAUX
  32.   80  CONTINUE
  33. C                      POUTRES
  34.       DO 810 I=1,6
  35. 810   SIGEL(I)=SIGT(I)
  36.       GO TO 740
  37.    79 CONTINUE
  38.       GO TO 740
  39.    78 CONTINUE
  40.       GO TO 740
  41.    77 SIGEL(1)=SIGT(1)
  42.       GO TO 740
  43.    76 DO 760 IA=1,4
  44.   760 SIGEL(IA)=SIGT(IA)
  45.       GO TO 740
  46.    75 CONTINUE
  47.       SIGEL(1)=B*(SIGT(1)-SIGT(2))
  48.       SIGEL(2)=B*(SIGT(1)+SIGT(2))
  49.       SIGEL(4)=SIGT(4)
  50.       SIGEL(3)=SIGT(3)
  51.       GO TO 740
  52.    74 CONTINUE
  53.       DO 7401 IB=1,3
  54.       SIGEL(IB)=0.
  55.       DO 7402 IC=1,3
  56.       IP=3*(IC-1)+IB
  57.  7402 SIGEL(IB)=SIGEL(IB)+VECPRO(IP)*SIGT(IC)
  58.  7401 SIGEL(3+IB)=SIGT(3+IB)
  59.       GO TO 740
  60.    71 CONTINUE
  61.       B1=B
  62.       B2=B
  63.       DO 710 I=1,IBO
  64.       IP=(I-1)*3
  65.       SIGEL(IP+1)=B1*SIGT(IP+1)-B2*SIGT(IP+2)
  66.       SIGEL(IP+2)=B2*SIGT(IP+1)+B1*SIGT(IP+2)
  67.   710 SIGEL(3+IP)=SIGT(3+IP)
  68.       GO TO 740
  69.    73 CONTINUE
  70.       SIGEL(3)=SIGT(3)
  71.       GO TO 740
  72.    70 CONTINUE
  73.       SIGT(1)=SIGT(1)*A
  74.       SIGT(2)=SIGT(2)*B
  75.       SIGT(3)=SIGT(3)*C
  76.       SIGEL(1)=SIGT(1)+SIGT(2)+SIGT(3)
  77.       SIGEL(2)=SIGT(1)-SIGT(2)+SIGT(3)
  78.       SIGEL(3)=SIGT(1)-2.*SIGT(3)
  79.       SIGEL(4)=SIGT(4)
  80.       SIGEL(5)=SIGT(5)
  81.       SIGEL(6)=SIGT(6)
  82.   740 CONTINUE
  83.       RETURN
  84.       END
  85.  
  86.  
  87.  

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