Télécharger la plaquette Cast3M


fr en

Télécharger hotan.eso

Retour à la liste

Numérotation des lignes : 

hotan

  1. C HOTAN     SOURCE    CB215821  19/08/01    21:16:06     10279          
  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.  
  26. -INC PPARAM
  27. -INC CCOPTIO
  28. *
  29.       CHARACTER*4 DEUXMOT(2)
  30. *     on lit un ou deux flottants xprec et dt
  31.       DATA DEUXMOT/'PREC','DT  '/
  32.       XPREC=1.D-3
  33.       DTPS=0.D0
  34.  1    CONTINUE
  35.       CALL LIRMOT(DEUXMOT,2,IVAL,0)
  36.       IF ( IVAL .EQ. 1) THEN
  37.         CALL LIRREE(XPREC,1,IRTFLO)
  38.         GOTO 1
  39.       ELSE IF ( IVAL .EQ. 2) THEN
  40.         CALL LIRREE(DTPS,1,IRTFLO)
  41.         GOTO 1
  42.       ENDIF
  43. *
  44. *     lecture d'un model
  45. *
  46.       CALL LIROBJ('MMODEL',IPMODL,0,IRTM)
  47.       IF(IRTM .EQ. 1)CALL ACTOBJ('MMODEL',IPMODL,1)
  48.       IF(IERR.NE.0) GOTO 666
  49. *
  50. *        lecture du mchaml de contraintes
  51. *
  52.       CALL LIROBJ('MCHAML',IPIN,1,IRT1)
  53.       CALL ACTOBJ('MCHAML',IPIN,1)
  54.       IF(IERR.NE.0) GOTO 666
  55.         CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  56.         IF(IR   .NE. 1) CALL ERREUR(KER)
  57.         IF(IERR .NE. 0) RETURN
  58. *
  59. *        lecture du mchaml de variables internes
  60. *
  61.       CALL LIROBJ('MCHAML',IPIN,1,IRT2)
  62.       CALL ACTOBJ('MCHAML',IPIN,1)
  63.       IF(IERR.NE.0) GOTO 666
  64.         CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
  65.         IF(IR   .NE. 1) CALL ERREUR(KER)
  66.         IF(IERR .NE. 0) RETURN
  67. *
  68. *        lecture du mchaml de caracteristiques
  69. *
  70.       CALL LIROBJ('MCHAML',IPIN,1,IRT3)
  71.       CALL ACTOBJ('MCHAML',IPIN,1)
  72.       IF(IERR.NE.0) GOTO 666
  73.         CALL REDUAF(IPIN,IPMODL,IPCHE3,0,IR,KER)
  74.         IF(IR   .NE. 1) CALL ERREUR(KER)
  75.         IF(IERR .NE. 0) RETURN
  76. *
  77.       CALL HOTANP(IPMODL,IPCHE1,IPCHE2,IPCHE3,XPREC,DTPS,IPCHOT,IRET)
  78.       IF(IERR .NE. 0) RETURN
  79. *
  80. *     ecriture du champs par element de matrice de hooke
  81. *
  82.       IF (IRET.EQ.1) THEN
  83.         CALL ACTOBJ('MCHAML  ',IPCHOT,1)
  84.         CALL ECROBJ('MCHAML  ',IPCHOT)
  85.       ENDIF
  86.   666 CONTINUE
  87.       END
  88.  
  89.  
  90.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales