C EQUIV     SOURCE    CHAT      11/03/16    21:21:25     6902           
      SUBROUTINE EQUIV
C=======================================================================C
C                                                                       C
C     OPERATEUR DE RIGIDITE POUR LES STRUCTURES EN BETON PRECONTRAINTES C
C                                                                       C
C     ENTREES:                                                          C 
C                                                                       C
C        MODL           MODELE POUR LE BETON                            C
C        MODL1          MODELES DES ACIERS                              C
C        MATR           CHARACTERISTIQUES MECANIQUES DES MATERIAUX      C
C                                                                       C
C     SORTIES :                                                         C
C                                                                       C
C       POINTEUR IPOI6  CREER UN OBJET DE RIGIDITE                      C
C                                                                       C
C-----------------------------------------------------------------------C
C             Fichier JIANG modifié FLEURET 05/96                       C
C=======================================================================C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)                                                  
-INC CCOPTIO                                                                    
-INC CCHAMP                                                                     
-INC SMCHAML
-INC SMELEME
-INC SMMODEL
-INC SMCHPOI
C
      SEGMENT IPNN
        INTEGER IPNNUM(NSOU1)
      ENDSEGMENT
C
      DIMENSION IPM(200)    
      CHARACTER*4 CMOT  
      CHARACTER*6 CMOT1
      IPCHE1=0                                                                  
      IPCHE2=0                                                                  
C
C --- ON RECUPERE EVENTUELLEMENT DES OBJETS MODELE
C
      DO 10 I=1,200
 10   IPM(I)=0
      CALL LIRCHA(CMOT1,0,IRETOU)
      NOMCA=1
      IMD=NOMCA+1
      DO 11 I=1,IMD
      IPMODL=0 
      CALL LIROBJ('MMODEL',IPMODL,0,IRET)
 11   IPM(I)=IPMODL
      MMODEL =IPM(2)
      SEGACT MMODEL                                                     
      NSOU1=KMODEL(/1) 
C      SEGINI IPNN
C
      IF(CMOT1.EQ.'RIGI') THEN                                   
C                                                                               
C --- ON LIT UN 1-ER CHAMELEM
C                                                                               
      IPCHA1 = 0                                                                
      CALL LIROBJ('MCHAML',IPCHA1,1,IRET1)
      IF (IERR.NE.0) RETURN                                                     
C                                                                               
C --- ON LIT EVENTUELLEMENT UN 2EME CHAMELEM                                   
C                                                                               
      IPCHA2 = 0                                                                
      CALL LIROBJ('MCHAML',IPCHA2,0,IRET1)  
C                                                                               
C --- RANGEMENT DES POINTEURS :  'CARACTERISTIQUES' DANS IPCHE1           
C                                'MATRICE DE HOOKE' DANS IPCHE2              
C                                                                               
        CALL RNGCHA(IPCHA1,IPCHA2,'CARACTERISTIQUES','MATRICE DE HOOKE',      
     1              IPCHE1,IPCHE2)                
        IF ( IERR.NE.0) RETURN                         
        IF (IPCHE2.NE.0) THEN                 
          IMAT = 2                            
        ELSE IF (IPCHE1.NE.0) THEN             
          IMAT = 1                         
        ELSE                             
          IPCHE1=IPCHA1                    
          IMAT = 1                      
        ENDIF
      ELSE IF(CMOT1.EQ.'FORCES') THEN
        IPCHEP=0 
        CALL LIROBJ('MCHAML',IPCHEP,0,IRET)
      END IF
C                                                                               
C --- APPEL AU CALCUL DES RIGIDITES ELEMENTAIRES


C                                                                               
 2000 CONTINUE                                                                  
       CALL SUBNNM(IPM,IPNN,IMD,IRET)
      IF(CMOT1.EQ.'RIGI') THEN
        CALL RIGIA(IMD,IPCHE1,IPCHE2,IMAT,IPOI6,IPM,IPNN,IRET)
        IF (IRET.EQ.1) CALL ECROBJ('RIGIDITE',IPOI6) 
      ELSE IF(CMOT1.EQ.'FORCES') THEN
        CALL LSIGMP(IMD,IPCHEP,IPCHE2,IPCHE3,IMAT,
     &     IPCHE4,IRET,IPM,IPNN,IPCHE)          
        IF (IRET.EQ.1) THEN                                                    
C 
C --- ATTRIBUTION D'UNE NATURE DISCRETE AU CHPO QUI SORT
C
          MCHPOI = IPCHE4
          SEGACT MCHPOI*MOD
          JATTRI(1) = 2
          SEGDES MCHPOI
          CALL ECROBJ('CHPOINT',IPCHE4)                                   
      END IF                              
      END IF
      SEGSUP IPNN
      END