Numérotation des lignes :

C RIGJOI    SOURCE    AM        11/05/26    21:16:38     6982C RIGLI2    SOURCE    CHAT      05/01/13    03:00:02     5004C RIGLI2    SOURCE    DJER      94/12/19    21:15:07     654      SUBROUTINE RIGJOI(NMATT,REL,LRE,WORK,IDIM,CMATE)C====================================================================CC     MODIFIE REL (LRE,LRE) MATRICE DE RIGIDITE EN INTRODUISANT DESC     RAIDEUR LOCALES EQUIVALENTESC====================================================================      IMPLICIT INTEGER(I-N)      IMPLICIT REAL*8(A-H,O-Z)      CHARACTER*8 CMATE      DIMENSION REL(LRE,LRE),WORK(NMATT)      DIMENSION RE2(3,3),RE3(6,6)C      CALL ZERO(REL,LRE,LRE)      XZER = 0.D0C      IF (IDIM.EQ.3) THEN      CALL ZERO(RE3,6,6)       DO I =1,6         RE3(I,I)=WORK(I+6)       END DO* SI anisotrope       IF (CMATE.EQ.'ANISOTRO') THEN        K=0        DO I =2,6         DO J =1,I-1           K=K+1           RE3(I,J)=WORK(K+12)           RE3(J,I)=RE3(I,J)         END DO        END DO       ENDIFC       DO I =1,6        DO J =1,6          REL(I,J)=RE3(I,J)          REL(I+6,J)=-1*RE3(I,J)          REL(I,J+6)=-1*RE3(I,J)          REL(I+6,J+6)=RE3(I,J)        END DO       END DO***      ELSE IF (IDIM.EQ.2) THEN      CALL ZERO(RE2,3,3)       DO I =1,3         RE2(I,I)=WORK(I+2)       END DO* SI anisotrope       IF (CMATE.EQ.'ANISOTRO') THEN        K=0        DO I =2,3         DO J =1,I-1           K=K+1           RE2(I,J)=WORK(K+5)           RE2(J,I)=RE2(I,J)         END DO        END DO       ENDIFC       DO I =1,3        DO J =1,3          REL(I,J)=RE2(I,J)          REL(I,J+3)=-1*RE2(I,J)          REL(I+3,J)=-1*RE2(I,J)          REL(I+3,J+3)=RE2(I,J)        END DO       END DO      ENDIF    5 CONTINUE**-------------------------------------------------------------*      WRITE(*,*) '**  REL APRES ( dans RIGJOI) ** '*      DO 112 IA=1,LRE*       WRITE(*,1991) (REL(IA,IC),IC=1,LRE)*  112 CONTINUE*1991  FORMAT(18(1X,E12.5))*-------------------------------------------------------------      RETURN      END

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