Esope source files

   1 : C HOTAN     SOURCE    CHAT      06/08/24    21:41:54     5529           
   2 :       SUBROUTINE HOTAN
   3 : *_____________________________________________________________________
   4 : *
   5 : *      creation d'un mchaml de matrice de hooke tangente
   6 : *
   7 : *
   8 : *
   9 : *     ho1=hotan mod1 si1 va1 ( ca1 ) (xprec) (flo1) (dt) (flo2) ;
  10 : *
  11 : *        mod1   modele de calcul, type mmodel
  12 : *        si1    champ par element de contraintes,type mchaml
  13 : *        va1    champ par element de variables internes,type mchaml
  14 : *        ca1    champ par element de caracteristiques ( materielle
  15 : *                 et/ou geometriques ),type mchaml
  16 : *        flo1   flottant (1.d-3 par defaut)
  17 : *        ho1    champ par element resultat, de type mchaml et de sous
  18 : *                 type matrice de hooke tangent
  19 : *
  20 : *
  21 : *_____________________________________________________________________
  22 : *
  23 :       IMPLICIT INTEGER(I-N)
  24 :       IMPLICIT REAL*8(A-H,O-Z)
  25 : -INC CCOPTIO
  26 : *
  27 :       CHARACTER*4 DEUXMOT(2)
  28 : *     on lit un ou deux flottants xprec et dt
  29 :       DATA DEUXMOT/'PREC','DT  '/                                               *
  30 :       XPREC=1.D-3
  31 :       DTPS=0.D0
  32 :  1    CONTINUE
  32 :       CALL LIRMOT(DEUXMOT,2,IVAL,0)
  33 :       IF ( IVAL .EQ. 1) THEN
  33 :         CALL LIRREE(XPREC,1,IRTFLO)
  34 :         GOTO 1
  35 :       ELSE IF ( IVAL .EQ. 2) THEN
  35 :         CALL LIRREE(DTPS,1,IRTFLO)
  36 :         GOTO 1
  37 :       ENDIF
  38 : *
  39 : *     lecture d'un model
  40 : *
  40 :       CALL LIROBJ('MMODEL',IPMODL,0,IRTM)
  41 : *
  42 : *        lecture du mchaml de contraintes
  43 : *
  43 :       CALL LIROBJ('MCHAML',IPCHE1,1,IRT1)
  44 :       IF(IERR.NE.0) GOTO 666
  45 : *
  46 : *        lecture du mchaml de variables internes
  47 : *
  47 :       CALL LIROBJ('MCHAML',IPCHE2,1,IRT2)
  48 :       IF(IERR.NE.0) GOTO 666
  49 : *
  50 : *        lecture du mchaml de caracteristiques
  51 : *
  51 :       CALL LIROBJ('MCHAML',IPCHE3,1,IRT3)
  52 :       IF(IERR.NE.0) GOTO 666
  53 : *
  53 :       CALL HOTANP(IPMODL,IPCHE1,IPCHE2,IPCHE3,XPREC,DTPS,IPCHOT,IRET)
  54 : *
  55 : *     ecriture du champs par element de matrice de hooke
  56 : *
  57 :       IF (IRET.EQ.1) THEN
  57 :             CALL ECROBJ('MCHAML',IPCHOT)
  58 :       ENDIF
  59 :   666 CONTINUE
  60 :       RETURN
  61 :       END
  62 : 
  63 :  
  64 :  
  65 :