Télécharger la plaquette Cast3M


fr en

Télécharger taille.eso

Retour à la liste

Numérotation des lignes : 

taille

  1. C TAILLE    SOURCE    CB215821  23/06/14    21:15:05     11676          
  2. c     mr millard      93    20:00:16     354
  3.  
  4.       SUBROUTINE TAILLE
  5.  
  6. c=====================================================================
  7. c    operateur calculant les parametres de taille par element
  8. c        aux points de gauss (tenseur de taille)
  9. c
  10. c    ancienne syntaxe :
  11. c    ------------------
  12. c
  13. c    chel = tail modl
  14. c            ob1 : modele de calcul, type mmodel
  15. c            chel: chamelem de sous type scalaires contenant
  16. c                    les parametres de taille aux points de gauss
  17. c
  18. c=====================================================================
  19. c
  20.       IMPLICIT INTEGER(I-N)
  21.       IMPLICIT REAL*8(A-H,O-Z)
  22.  
  23.  
  24. -INC PPARAM
  25. -INC CCOPTIO
  26. -INC SMCOORD
  27. *
  28.       CHARACTER*4 MUNIF(3)
  29.       DATA MUNIF/'DIAM','DIRE','UNIF'/
  30. c
  31. c           ipmodl          pointeur sur un modele mmodel
  32. c           ipche           pointeur sur un chamelem  de
  33. c                           sous type scalaires contenant
  34. c                           les parametres de taille
  35. c
  36.       IPMODL = 0
  37.       IUNIF  = 0
  38. c
  39. c Lecture du mot cle
  40. c
  41.       CALL LIRMOT(MUNIF,3,IUNIF,0)
  42.       IF (IERR.NE.0) RETURN
  43. c
  44. c Lecture du modele
  45. c
  46.       CALL LIROBJ('MMODEL',IPMODL,1,IRET)
  47.       CALL ACTOBJ('MMODEL',IPMODL,1)
  48.       IF (IERR.NE.0) RETURN
  49.  
  50.       SEGACT,MCOORD
  51. C
  52. C Option : 'DIAMETRE_MIN'
  53.       IF (IUNIF.EQ.1) THEN
  54.         IPCHA1 = 0
  55.         IPCHA2 = 0
  56.         IPCHA3 = 0
  57.         IPCHA4 = 0
  58.         CALL LIROBJ('MCHAML',IPIN,0,IRET)
  59.         IF (IERR.NE.0) RETURN
  60.  
  61.         IPCHA1=0
  62.         IF (IRET .EQ. 1) THEN
  63.           CALL REDUAF(IPIN,IPMODL,IPCHA1,0,IR,KER)
  64.           IF(IR   .NE. 1) CALL ERREUR(KER)
  65.           IF(IERR .NE. 0) RETURN
  66.         ENDIF
  67.  
  68.         ICAS = 5
  69.         CALL CFL1(IPMODL,IPCHA1,IPCHA2,IPCHA3,IPCHA4,ICAS)
  70.         IF (IPCHA4 .NE. 0) THEN
  71.           CALL ACTOBJ('MCHAML',IPCHA4,1)
  72.           CALL ECROBJ('MCHAML',IPCHA4)
  73.         ENDIF
  74.  
  75. C Option : 'DIRECTION' & Cas ou le mot cle UNIF est seul
  76.       ELSE
  77.         IF (IUNIF.EQ.2) THEN
  78. C Lecture eventuelle du mot cle 'UNIF' apres le mot cle 'DIRECTION'
  79.           IUNIF = 0
  80.           CALL LIRMOT(MUNIF(3),1,IUNIF,0)
  81.           IF (IERR.NE.0) RETURN
  82.         ENDIF
  83.  
  84. c Calcul des parametres de taille par element a chaque point de Gauss
  85. c Si IUNIF non nul 0, toutes les composantes sont nulles.
  86.         IPCHE = 0
  87.         IRET  = 0
  88.         CALL TAILPO(IPMODL,IPCHE,IUNIF,IRET)
  89.         IF (IRET.NE.0) THEN
  90.           CALL ACTOBJ('MCHAML',IPCHE,1)
  91.           CALL ECROBJ('MCHAML',IPCHE)
  92.         ENDIF
  93.       ENDIF
  94.  
  95.       SEGDES,MCOORD
  96.  
  97.       RETURN
  98.       END
  99.  
  100.  

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