Esope source files

   1 : C FOFISS    SOURCE    FANDEUR   10/12/14    21:16:29     6812           
   2 : C                                                                               
   3 :       SUBROUTINE FOFISS                                                         
   4 : C                                                                               
   5 : C=======================================================================        
   6 : C                                                                               
   7 : C            OPERATEUR FORCES INTERNES                                          
   8 : C                                                                               
   9 : C          CET OPERATEUR NOUS PERMET DE CALCULER LE CHPOIN FOR1                 
  10 : C          DES FORCES NODALES RESULTANT DE L'INTEGRATION DU CHA-                
  11 : C          MELEM SIG1 DE CONTRAINTES ISSU PAR EXEMPLE DE L'OPE-                 
  12 : C          RATEUR  SIGMA  ET DU CHAMELEM GRAD1 DE GRADIENTS ( ET                
  13 : C          GRAF1 DE GRAFLEXIS ) ISSU PAR EXEMPLE DE L'OPERATEUR                 
  14 : C          GRAD ( ET GRAF ).                                                    
  15 : C                                                                               
  16 : C          SYNTAXE  : FOR1=FOFISS MOD1 SIG1 GRAD1 (GRAF1)(CAR1)  ;              
  17 : C                                                                               
  18 : C         ENTREES :                                                             
  19 : C                MOD1=OBJET DE TYPE MMODEL                                      
  20 : C                SIG1=CHAPM /ELMTS DE CONTRAINTES (TYPE MCHAML)                 
  21 : C                GRAD1=CHAPM /ELMTS DE GRADIENTS  (TYPE MCHAML)                 
  22 : C                GRAF1=CHAPM /ELMTS DE GRAFLEXIS  (TYPE MCHAML)                 
  23 : C                CAR1=CHAMP /ELMTS DE CARACTERISTIQUES (FACULTATIF)             
  24 : C                POUR LES POUTRES ET LES COQUES (TYPE MCHAML)                   
  25 : C         SORTIES:                                                              
  26 : C                  FOR1 =CHAMP PAR POINT DONNANT LES FORCES NODALES             
  27 : C    CODE JUILLET 87 SUO X.Z.                                                   
  28 : C    PASSAGE AUX NOUVEAUX MCHAMLS PAR P.DOWLATYARI SEP 91                       
  29 : C=======================================================================        
  30 : C                                                                               
  31 :       IMPLICIT INTEGER(I-N)
  32 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  33 : C                                                                               
  34 : -INC CCOPTIO                                                                    
  35 : C                                                                               
  36 :        IPMODL=0                                                                 
  37 :        IPCHE1=0                                                                 
  38 :        IPCHE2=0                                                                 
  39 :        IPCHE3=0                                                                 
  40 :        IPCHE4=0                                                                 
  41 :        IPCHA3=0                                                                 
  42 :        IPCHA4=0                                                                 
  43 : C                                                                               
  44 : C  ****  LECTURE D'UN OBJET MMODEL                                              
  45 : C                                                                               
  45 :        CALL LIROBJ('MMODEL',IPMODL,1,IRT1)                                      
  46 :        IF(IERR.NE.0)RETURN                                                      
  47 : C                                                                               
  48 : C  ****  LECTURE DU PREMIER CHAMP/ELEMENT                                       
  49 : C                                                                               
  49 :         CALL LIROBJ('MCHAML',IPCHE1,1,IRT1)                                     
  50 :         IF(IERR.NE.0)RETURN                                                     
  51 : C                                                                               
  52 : C  ****  LECTURE DU DEUXIEME CHAMP/ELEMENT                                      
  53 : C                                                                               
  53 :         CALL LIROBJ('MCHAML',IPCHE2,1,IRT1)                                     
  54 :         IF(IERR.NE.0) RETURN                                                    
  55 : C                                                                               
  56 : C  ****  LECTURE EVENTUELLE DU TROISIEME CHAMP                                  
  57 : C                                                                               
  57 :         CALL LIROBJ('MCHAML',IPCHA3,0,IRT1)                                     
  58 :         IF(IERR.NE.0)RETURN                                                     
  59 : C                                                                               
  60 : C  ****  LECTURE EVENTUELLE DU QUATRIEME CHAMP                                  
  61 : C                                                                               
  61 :         CALL LIROBJ('MCHAML',IPCHA4,0,IRT1)                                     
  62 :         IF(IERR.NE.0)RETURN                                                     
  63 : C
  64 :        IF (IPCHA3.NE.0.OR.IPCHA4.NE.0) THEN                                     
  64 :          CALL RNGCHA(IPCHA3,IPCHA4,'GRAFLEXI',                                  
  65 :      1                            'CARACTERISTIQUES',IPCHE3,IPCHE4)             
  66 :          IF(IERR.NE.0) RETURN
  67 :        ENDIF                                                                    
  68 : C                                                                               
  68 :       CALL FOFIS1(IPMODL,IPCHE1,IPCHE2,IPCHE3,IPCHE4,IPCHP1,IRET)               
  69 : C                                                                               
  69 :       IF(IRET.EQ.1) CALL ECROBJ('CHPOINT',IPCHP1)                               
  70 : C
  71 :       RETURN                                                                    
  72 :       END                                                                       
  73 : 
  74 :