Esope source files

   1 : C SIGJOI    SOURCE    AM        11/05/26    21:16:41     6982           
   2 : C SIGLIA    SOURCE    CHAT      05/01/13    03:16:39     5004           
   3 : C SIGLIA    SOURCE    DJER     23/01/96
   4 :       SUBROUTINE SIGJOI(NMATT,XDDL,WORK,LRE,XSTRS,IDIM,NSTRS,CMATE)
   5 : C====================================================================           
   6 : C     CALCULE DES CONTRAINTES  POUR JOI1 :
   7 : C      element de liaison a 2 noeuds (6 ddl par noeuds)                                                        
   8 : C     (ON CALCUL en fait DES EFFORTs NODAUX ) 
   9 : C====================================================================           
  10 :       IMPLICIT INTEGER(I-N)
  11 :       IMPLICIT REAL*8(A-H,O-Z)
  12 :       CHARACTER*8 CMATE                                                       
  13 :       DIMENSION XDDL(LRE),WORK(NMATT),XSTRS(*),DEPL(12)
  14 :       DIMENSION HK3(6,6),HK2(3,3)
  15 : *-------------------------------------------------------------
  16 :        DO I =1,LRE
  17 :          DEPL(I) = XDDL(I)
  18 :        END DO
  19 : C
  20 :       IF (IDIM.EQ.3) THEN
  21 : C CAS D'UN MATERIAU ORTHOTROPE
  22 : C CALCUL DES EFFORTS NODAUX
  23 :       IF (CMATE.EQ.'ORTHOTRO') THEN     
  24 :        DO I =1,NSTRS
  25 :          XSTRS(I)=WORK(I+6)*(DEPL(I+NSTRS)-DEPL(I))
  26 :        END DO
  27 : C CAS D'UN MATERIAU ANISOTROPE
  28 :       ELSE IF (CMATE.EQ.'ANISOTRO') THEN 
  29 : C INSERTION DES RAIDEURS MATERIELLES DANS UNE MATRICE 6*6
  30 :         K = 1      
  31 :         DO I=1,NSTRS
  32 :          HK3(I,I)=WORK(I+6)
  33 :          DO J=1,I-1
  34 :           HK3(I,J)=WORK(K+12)
  35 :           HK3(J,I)=HK3(I,J)
  36 :           K = K+1
  37 :          ENDDO
  38 :         ENDDO
  39 : C CALCUL DES EFFORTS NODAUX
  40 :         DO I=1,NSTRS
  41 :          DO J=1,NSTRS
  42 :           XSTRS(I)=XSTRS(I)+HK3(I,J)*(DEPL(J+NSTRS)-DEPL(J))
  43 :          ENDDO
  44 :         ENDDO
  45 :       ENDIF
  46 : C
  47 : *
  48 : C
  49 :       ELSE IF (IDIM.EQ.2) THEN
  50 : C CAS D'UN MATERIAU ORTHOTROPE
  51 : C CALCUL DES EFFORTS NODAUX
  52 :       IF (CMATE.EQ.'ORTHOTRO') THEN
  53 :        DO I =1,NSTRS
  54 :          XSTRS(I)=WORK(I+2)*(DEPL(I+NSTRS)-DEPL(I))
  55 :        END DO
  56 : C CAS D'UN MATERIAU ANISOTROPE
  57 :       ELSE IF (CMATE.EQ.'ANISOTRO') THEN
  58 : C INSERTION DES RAIDEURS MATERIELLES DANS UNE MATRICE 3*3
  59 :         K = 1      
  60 :         DO I=1,NSTRS
  61 :          HK2(I,I)=WORK(I+2)
  62 :          DO J=1,I-1
  63 :           HK2(I,J)=WORK(K+5)
  64 :           HK2(J,I)=HK2(I,J)
  65 :           K = K+1
  66 :          ENDDO
  67 :         ENDDO
  68 : C CALCUL DES EFFORTS NODAUX
  69 :         DO I=1,NSTRS
  70 :          DO J=1,NSTRS
  71 :           XSTRS(I)=XSTRS(I)+HK2(I,J)*(DEPL(J+NSTRS)-DEPL(J))
  72 :          ENDDO
  73 :         ENDDO
  74 :       ENDIF
  75 :       ENDIF
  76 : *
  77 : *
  78 : *
  79 : *-------------------------------------------------------------
  80 : *      WRITE(*,*) '**  EFFORT  ( dans siglia) ** '
  81 : *      WRITE(*,1992) (XSTRS(IC),IC=1,12)
  82 : *1992  FORMAT(6(6(1X,E12.5),/))
  83 : *-------------------------------------------------------------
  84 :       RETURN                                                                    
  85 :       END                                                                       
  86 :  
  87 :  
  88 :  
  89 :  
  90 :