Télécharger taille.eso

Retour à la liste

Numérotation des lignes :

taille
  1. C TAILLE SOURCE PV 20/03/31 14:34:08 10567
  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. segact mcoord
  37. IPMODL = 0
  38. IUNIF = 0
  39. c
  40. c Lecture du mot cle
  41. c
  42. CALL LIRMOT(MUNIF,3,IUNIF,0)
  43. IF (IERR.NE.0) RETURN
  44. c
  45. c Lecture du modele
  46. c
  47. CALL LIROBJ('MMODEL',IPMODL,1,IRET)
  48. IF (IERR.NE.0) RETURN
  49. C
  50. C Option : 'DIAMETRE_MIN'
  51. IF (IUNIF.EQ.1) THEN
  52. IPCHA1 = 0
  53. IPCHA2 = 0
  54. IPCHA3 = 0
  55. IPCHA4 = 0
  56. CALL LIROBJ('MCHAML',IPIN,0,IRET)
  57. IF (IERR.NE.0) RETURN
  58. IPCHA1=0
  59. IF (IRET .EQ. 1) THEN
  60. CALL REDUAF(IPIN,IPMODL,IPCHA1,0,IR,KER)
  61. IF(IR .NE. 1) CALL ERREUR(KER)
  62. IF(IERR .NE. 0) RETURN
  63. ENDIF
  64.  
  65. ICAS = 5
  66. CALL CFL1(IPMODL,IPCHA1,IPCHA2,IPCHA3,IPCHA4,ICAS)
  67. IF (IPCHA4 .NE. 0) CALL ECROBJ('MCHAML',IPCHA4)
  68.  
  69. C Option : 'DIRECTION' & Cas ou le mot cle UNIF est seul
  70. ELSE
  71. IF (IUNIF.EQ.2) THEN
  72. C Lecture eventuelle du mot cle 'UNIF' apres le mot cle 'DIRECTION'
  73. IUNIF = 0
  74. CALL LIRMOT(MUNIF(3),1,IUNIF,0)
  75. IF (IERR.NE.0) RETURN
  76. ENDIF
  77. c Calcul des parametres de taille par element a chaque point de Gauss
  78. c Si IUNIF non nul 0, toutes les composantes sont nulles.
  79. IPCHE = 0
  80. IRET = 0
  81. CALL TAILPO(IPMODL,IPCHE,IUNIF,IRET)
  82. IF (IRET.NE.0) CALL ECROBJ('MCHAML',IPCHE)
  83. ENDIF
  84.  
  85. RETURN
  86. END
  87.  
  88.  
  89.  
  90.  

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