Télécharger la plaquette Cast3M


fr en

Télécharger rigjoi.eso

Retour à la liste

Numérotation des lignes : 

rigjoi

  1. C RIGJOI    SOURCE    AM        11/05/26    21:16:38     6982
  2. C RIGLI2    SOURCE    CHAT      05/01/13    03:00:02     5004
  3. C RIGLI2    SOURCE    DJER      94/12/19    21:15:07     654
  4.       SUBROUTINE RIGJOI(NMATT,REL,LRE,WORK,IDIM,CMATE)
  5. C====================================================================
  6. C
  7. C     MODIFIE REL (LRE,LRE) MATRICE DE RIGIDITE EN INTRODUISANT DES
  8. C     RAIDEUR LOCALES EQUIVALENTES
  9. C====================================================================
  10.       IMPLICIT INTEGER(I-N)
  11.       IMPLICIT REAL*8(A-H,O-Z)
  12.       CHARACTER*8 CMATE
  13.       DIMENSION REL(LRE,LRE),WORK(NMATT)
  14.       DIMENSION RE2(3,3),RE3(6,6)
  15. C
  16.       CALL ZERO(REL,LRE,LRE)
  17.       XZER = 0.D0
  18. C
  19.       IF (IDIM.EQ.3) THEN
  20.       CALL ZERO(RE3,6,6)
  21.        DO I =1,6
  22.          RE3(I,I)=WORK(I+6)
  23.        END DO
  24. * SI anisotrope
  25.        IF (CMATE.EQ.'ANISOTRO') THEN
  26.         K=0
  27.         DO I =2,6
  28.          DO J =1,I-1
  29.            K=K+1
  30.            RE3(I,J)=WORK(K+12)
  31.            RE3(J,I)=RE3(I,J)
  32.          END DO
  33.         END DO
  34.        ENDIF
  35. C
  36.        DO I =1,6
  37.         DO J =1,6
  38.           REL(I,J)=RE3(I,J)
  39.           REL(I+6,J)=-1*RE3(I,J)
  40.           REL(I,J+6)=-1*RE3(I,J)
  41.           REL(I+6,J+6)=RE3(I,J)
  42.         END DO
  43.        END DO
  44. *
  45. *
  46. *
  47.       ELSE IF (IDIM.EQ.2) THEN
  48.       CALL ZERO(RE2,3,3)
  49.        DO I =1,3
  50.          RE2(I,I)=WORK(I+2)
  51.        END DO
  52. * SI anisotrope
  53.        IF (CMATE.EQ.'ANISOTRO') THEN
  54.         K=0
  55.         DO I =2,3
  56.          DO J =1,I-1
  57.            K=K+1
  58.            RE2(I,J)=WORK(K+5)
  59.            RE2(J,I)=RE2(I,J)
  60.          END DO
  61.         END DO
  62.        ENDIF
  63. C
  64.        DO I =1,3
  65.         DO J =1,3
  66.           REL(I,J)=RE2(I,J)
  67.           REL(I,J+3)=-1*RE2(I,J)
  68.           REL(I+3,J)=-1*RE2(I,J)
  69.           REL(I+3,J+3)=RE2(I,J)
  70.         END DO
  71.        END DO
  72.       ENDIF
  73.     5 CONTINUE
  74. *
  75. *-------------------------------------------------------------
  76. *      WRITE(*,*) '**  REL APRES ( dans RIGJOI) ** '
  77. *      DO 112 IA=1,LRE
  78. *       WRITE(*,1991) (REL(IA,IC),IC=1,LRE)
  79. *  112 CONTINUE
  80. *1991  FORMAT(18(1X,E12.5))
  81. *-------------------------------------------------------------
  82.       RETURN
  83.       END
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  

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