Esope source files

   1 : C EPTH      SOURCE    CHAT      06/08/24    21:31:12     5529           
   2 : 
   3 : C=======================================================================
   4 : C=                              E P T H                                =
   5 : C=                              -------                                =
   6 : C=                                                                     =
   7 : C=  OPERATEUR CAST3M "EPTH" :                                          =
   8 : C=  -------------------------                                          =
   9 : C=   EPT1  =  'EPTH'  MODL1  |  CHP1   |  CARA1  ;                     =
  10 : C=                           |  CHEL1  |                               =
  11 : C=                                                                     =
  12 : C=   Cet operateur sert a calculer les deformations dues a un champ    =
  13 : C=   de temperatures.                                                  =
  14 : C=                                                                     =
  15 : C=  ARGUMENTS :                                                        =
  16 : C=  -----------                                                        =
  17 : C=   MODL1  (MMODEL)    Modele (global) associe a la structure         =
  18 : C=   CHP1   (CHPOINT)   Temperatures aux NOEUDS                        =
  19 : C=   CHEL1  (MCHAML)    Temperatures donnees par ELEMENT               =
  20 : C=                      Sous-type 'TEMPERATURES'                       =
  21 : C=   CARA1  (MCHAML)    Caracteristiques des materiaux                 =
  22 : C=                      Sous-type 'CARACTERISTIQUES'                   =
  23 : C=                                                                     =
  24 : C=  RESULTAT :                                                         =
  25 : C=  ----------                                                         =
  26 : C=   EPT1   (CHPOINT/MCHAML)   Deformations d'origine thermique        =
  27 : C=======================================================================
  28 : 
  29 :       SUBROUTINE EPTH
  30 : 
  31 :       IMPLICIT INTEGER(I-N)
  32 :       IMPLICIT REAL*8 (A-H,O-Z)
  33 : 
  34 : -INC CCOPTIO
  35 : 
  36 : C  1 - LECTURE DES ARGUMENTS DE L'OPERATEUR
  37 : C ==========================================
  38 : C  1.1 - Lecture OBLIGATOIRE du modele (IPMODL)
  39 : C =====
  40 :       MOTERR(1:8)='MODELE'
  40 :       CALL MESLIR(-137)
  40 :       CALL LIROBJ('MMODEL',IPMODL,1,IRet)
  41 :       IF (IERR.NE.0) RETURN
  42 : C =====
  43 : C  1.2 - Lecture OBLIGATOIRE du champ de caracteristiques (IPCHA1)
  44 : C =====
  45 :       IPCHA1=0
  45 :       CALL LIROBJ('MCHAML',IPCHA1,1,IRet)
  46 :       IF (IERR.NE.0) RETURN
  47 : C =====
  48 : C  1.3 - Lecture OBLIGATOIRE du champ de temperatures donne par :
  49 : C           1) un CHPOINT (IPCHPO) converti en MCHAML (IPCHE2)
  50 : C        ou 2) un MCHAML (IPCHA2) qui est duplique en IPCHE2
  51 : C =====
  51 :       CALL LIROBJ('CHPOINT',IPCHPO,0,IRet)
  52 :       IF (IERR.NE.0) RETURN
  53 :       IF (IRet.NE.0) THEN
  54 :         IPCHE1=IPCHA1
  54 :         CALL CHAME1(0,IPMODL,IPCHPO,' ',IPCHE2,1)
  55 :         IF (IERR.NE.0) GOTO 10
  56 :       ELSE
  57 : C= 1.3.1 - Mise en ordre des MCHAML lus
  58 : C          1 -> IPCHE1 = Caracteristiques, 2 -> IPCHE2 = Temperatures
  59 :         IPCHA2=0
  59 :         CALL LIROBJ('MCHAML', IPCHA2,1,IRet)
  60 :         IF (IERR.NE.0) RETURN
  61 :         IPCHE1=0
  62 :         IPCHE2=0
  62 :         CALL RNGCHA(IPCHA1,IPCHA2,'CARACTERISTIQUES','TEMPERATURES',
  63 :      .               IPCHE1,IPCHE2)
  64 : C= 1.3.2 - ERREUR si les sous-types ne sont pas corrects
  65 :         IF (IPCHE1.EQ.0.OR.IPCHE2.EQ.0) THEN
  65 :           CALL ERREUR(554)
  66 :           RETURN
  67 :         ENDIF
  68 : C= 1.3.3 - Copie du MCHAML de temperatures pour la linearisation
  68 :         CALL ECROBJ('MCHAML',IPCHE2)
  68 :         CALL COPIER
  68 :         CALL LIROBJ('MCHAML',IPCHE2,1,IRet)
  69 :         IF (IERR.NE.0) GOTO 10
  70 :       ENDIF
  71 : 
  72 : C  2 - LINEARISATION DES TEMPERATURES POUR LES NOEUDS MILIEUX
  73 : C ============================================================
  73 :       CALL LINEAT(IPCHE2)
  74 : 
  75 : C  3 - CALCUL DU MCHAML DE DEFORMATIONS THERMIQUES
  76 : C =================================================
  76 :       CALL EPTHP(IPMODL,IPCHE1,IPCHE2,IPEPTH,IRET)
  77 : 
  78 : C  4 - ECRITURE DU MCHAML RESULTAT
  79 : C =================================
  79 :       IF (IRET.EQ.1) CALL ECROBJ('MCHAML',IPEPTH)
  79 :  10   CALL DTCHAM(IPCHE2)
  80 : 
  81 :       RETURN
  82 :       END
  83 : 
  84 :  
  85 :  
  86 :