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