Télécharger taille.eso

Retour à la liste

Numérotation des lignes :

  1. C TAILLE SOURCE CB215821 16/12/05 22:04:37 9237
  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. -INC CCOPTIO
  24. *
  25. CHARACTER*4 MUNIF(3)
  26. DATA MUNIF/'DIAM','DIRE','UNIF'/
  27. c
  28. c ipmodl pointeur sur un modele mmodel
  29. c ipche pointeur sur un chamelem de
  30. c sous type scalaires contenant
  31. c les parametres de taille
  32. c
  33. IPMODL = 0
  34. IUNIF = 0
  35. c
  36. c Lecture du mot cle
  37. c
  38. CALL LIRMOT(MUNIF,3,IUNIF,0)
  39. IF (IERR.NE.0) RETURN
  40. c
  41. c Lecture du modele
  42. c
  43. CALL LIROBJ('MMODEL',IPMODL,1,IRET)
  44. IF (IERR.NE.0) RETURN
  45. C
  46. C Option : 'DIAMETRE_MIN'
  47. IF (IUNIF.EQ.1) THEN
  48. IPCHA1 = 0
  49. IPCHA2 = 0
  50. IPCHA3 = 0
  51. IPCHA4 = 0
  52. CALL LIROBJ('MCHAML',IPIN,0,IRET)
  53. IF (IERR.NE.0) RETURN
  54. IPCHA1=0
  55. IF (IRET .EQ. 1) THEN
  56. CALL REDUAF(IPIN,IPMODL,IPCHA1,0,IR,KER)
  57. IF(IR .NE. 1) CALL ERREUR(KER)
  58. IF(IERR .NE. 0) RETURN
  59. ENDIF
  60.  
  61. ICAS = 5
  62. CALL CFL1(IPMODL,IPCHA1,IPCHA2,IPCHA3,IPCHA4,ICAS)
  63. IF (IPCHA4 .NE. 0) CALL ECROBJ('MCHAML',IPCHA4)
  64.  
  65. C Option : 'DIRECTION' & Cas ou le mot cle UNIF est seul
  66. ELSE
  67. IF (IUNIF.EQ.2) THEN
  68. C Lecture eventuelle du mot cle 'UNIF' apres le mot cle 'DIRECTION'
  69. IUNIF = 0
  70. CALL LIRMOT(MUNIF(3),1,IUNIF,0)
  71. IF (IERR.NE.0) RETURN
  72. ENDIF
  73. c Calcul des parametres de taille par element a chaque point de Gauss
  74. c Si IUNIF non nul 0, toutes les composantes sont nulles.
  75. IPCHE = 0
  76. IRET = 0
  77. CALL TAILPO(IPMODL,IPCHE,IUNIF,IRET)
  78. IF (IRET.NE.0) CALL ECROBJ('MCHAML',IPCHE)
  79. ENDIF
  80.  
  81. RETURN
  82. END
  83.  
  84.  
  85.  

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