C ELASTI    SOURCE    FANDEUR   10/12/17    21:16:30     6815           

      SUBROUTINE ELASTI                                                         

*_______________________________________________________________________        
*                                                                               
*                     OPERATEUR ELASTICITE                                      
*                                                                               
*      SIG = ELAS MODL EPS MATER ( VA1 )                                     
*                OU                         EN ELASTICITE LINEAIRE              
*      EPS = ELAS MODL SIG MATER ( VA1 )                                     
*                                                                               
*            MODL   MODELE DE CALCUL,TYPE MMODEL                                
*            EPS    CHAMP PAR ELEMENT DE DEFORMATIONS,TYPE MCHAML              
*            SIG    CHAMP PAR ELEMENT DE CONTRAINTES,TYPE MCHAML               
*            MATER  CHAMP PAR ELEMENT DE CARACTERISTIQUE,TYPE MCHAML 
*            VA1    CHAMP PAR ELEMENT DE VARIABLES INTERNES,TYPE MCHAML                   
*_______________________________________________________________________        
*
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)                                                  
*
-INC CCOPTIO                                                                    
*
      IPCHE1=0                                                                  
      IPCHE2=0                                                                  
      IPCHE3=0                                                                  
      IPSTRS=0
      IRETOU=0
      IRET=0
*                                                                               
*     LECTURE DU MODELE                                                         
*                                                                               
      CALL LIROBJ('MMODEL',IPMODL,1,IRETOU)
      IF(IERR.NE.0) RETURN
*                                                                               
*     ON LIT LE 1-ER CHAMP/ELEMENT
*
      CALL LIROBJ('MCHAML',IPCHE1,1,IRETOU)
      IF(IERR.NE.0) RETURN
*
*     LECTURE DU MCHAML DE MATERIAU
*
      CALL LIROBJ('MCHAML',IPCHE2,1,IRETOU)
      IF(IERR.NE.0) RETURN
C                                                                               
C       LECTURE DU MCHAML DE VARIABLES INTERNES        
C                                                                               
      CALL LIROBJ('MCHAML',IPCHE3,0,IRETOU)                                     
      IF (IERR.NE.0) RETURN                                                   
*                                                                               
*     CALCUL DES CONTRAINTES OU DES DEFORMATIONS                                
* 
      CALL ELAS1 (IPMODL,IPCHE1,IPCHE2,IPCHE3, IPSTRS,IRET)  
* 
      IF (IRET.EQ.1) CALL ECROBJ('MCHAML',IPSTRS)

      RETURN                                                                    
      END