Esope source files

   1 : C WORK      SOURCE    FANDEUR   10/12/14    21:20:05     6812           
   2 :       SUBROUTINE WORK                                                           
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 : C=======================================================================        
   6 : C                                                                               
   7 : C                   DENSITE D'ENERGIE                                           
   8 : C     CALCULE LE MCHAML   PRODUIT   SIGMA * GRADIENT                            
   9 : C                                                                               
  10 : C              SYNTAXE  WORK1=WORK MODL1  SIG1  GRAD1( GRAF1)                   
  11 : C      MODL1= OBJET DE TYPE MMODEL                                              
  12 : C      SIG1 = CHAMP DE CONTRAINTE (TYPE MCHAML)                                 
  13 : C      GRAD1= CHAMP DE GRADIENT (TYPE MCHAML)                                   
  14 : C      GRAF1= CHAMP DE GRADIENT DE FLEXION (TYPE MCHAML)                        
  15 : C      WORK1=CHAMP DE DENSITE D'ENERGIE (TYPE MCHAML)                           
  16 : C     CODE DE SUO X.Z.                                                          
  17 : C     PASSAGE AUX NOUVEAUX CHAMELEMS PAR P. DOWLATYARI LE 25/4/91               
  18 : C=======================================================================        
  19 : C                                                                               
  20 : -INC CCOPTIO                                                                    
  21 : C                                                                               
  22 :        IPMODL=0                                                                 
  23 :        IPCHE1=0                                                                 
  24 :        IPCHE2=0                                                                 
  25 :        IPCHE3=0                                                                 
  26 :        IPCHE4=0                                                                 
  27 : C                                                                               
  28 : C    LECTURE D'UN OBJET MMODEL                                                  
  29 : C                                                                               
  29 :        CALL LIROBJ('MMODEL',IPMODL,1,IRET)                                      
  30 :        IF(IERR.NE.0)RETURN                                                      
  31 : C                                                                               
  32 : C  ****  LECTURE DU PREMIER CHAMP/ELEMENT                                       
  33 : C                                                                               
  33 :         CALL LIROBJ('MCHAML',IPCHE1,1,IRET)                                     
  34 :         IF(IERR.NE.0)RETURN                                                     
  35 : C                                                                               
  36 : C  ****  LECTURE DU DEUXIEME CHAMP/ELEMENT                                      
  37 : C                                                                               
  37 :         CALL LIROBJ('MCHAML',IPCHE2,1,IRET)                                     
  38 :         IF(IERR.NE.0) RETURN                                                    
  39 : C                                                                               
  40 : C  ****  LECTURE DU TROISIEME CHAMP SI COQUE , QUI DOIT ETRE GRAF               
  41 : C                                                                               
  41 :         CALL LIROBJ('MCHAML',IPCHE3,0,IRET)                                     
  42 :         IF(IERR.NE.0)RETURN                                                     
  43 : C                                                                               
  43 :       CALL WORKP(IPMODL,IPCHE1,IPCHE2,IPCHE3,IPCHE4,IRET)                       
  44 : C                                                                               
  44 :       IF(IRET.NE.0) CALL ECROBJ('MCHAML',IPCHE4)                                
  45 : C                                                                               
  46 :       RETURN                                                                    
  47 :       END                                                                       
  48 : 
  49 :