Télécharger vloc.eso

Retour à la liste

Numérotation des lignes :

  1. C VLOC SOURCE CB215821 19/08/01 21:16:43 10279
  2. C
  3. SUBROUTINE VLOC
  4. C=====================================================================
  5. C
  6. C Fonction : CALCULE LES VECTEURS D'ORTHOTROPIE
  7. C
  8. C Syntaxe : CHAM123 = VLOC MOD1 MAT1;
  9. C
  10. C Input : MOD1 : MODELE de calcul , type MMODEL
  11. C
  12. C Output : CHAM123 : CHAMELEM aux POINTS DE GAUSS de RIGIDITE
  13. C de sous type VECTEUR LOCAUX (de composantes
  14. C V1X V1Y V1Z V2X V2Y V3X V3Y V3Z par ex.)
  15. C contenant les vecteurs V1 V2 V3
  16. c de base du repere local d'orthotropie
  17. C
  18. C Creation : BP, 2017-01-12 (inspiré de VSUR, JACONO et RIGI3)
  19. C Modifs : ...merci de compléter...
  20. C
  21. C=====================================================================
  22. C
  23. IMPLICIT INTEGER(I-N)
  24. IMPLICIT REAL*8(A-H,O-Z)
  25.  
  26.  
  27. -INC PPARAM
  28. -INC CCOPTIO
  29.  
  30. IPMODL=0
  31. IPCHE1=0
  32. IPCHE2=0
  33. IPCHE3=0
  34. IRET=0
  35. C
  36. C=====================================================================
  37. C LECTURE DU MODELE
  38. C=====================================================================
  39. C
  40. CALL LIROBJ('MMODEL ',IPMODL,1,IRET1)
  41. CALL ACTOBJ('MMODEL ',IPMODL,1)
  42. IF (IERR.NE.0) RETURN
  43. C
  44. C=====================================================================
  45. C LECTURE DU MATERIAU
  46. C=====================================================================
  47. C
  48. CALL LIROBJ('MCHAML ',IPMATE,1,IRET2)
  49. CALL ACTOBJ('MCHAML ',IPMATE,1)
  50. IF (IERR.NE.0) RETURN
  51. C
  52. C=====================================================================
  53. C CALCUL DES VECTEURS :
  54. C=====================================================================
  55. C DU REPERE LOCAL D'ORTHOTROPIE
  56. CALL REDUAF(IPMATE,IPMODL,IPMAT2,0,IRET3,KERR)
  57. IF (IRET3.NE.1) CALL ERREUR(KERR)
  58. IF(IERR.NE.0) RETURN
  59. c CALL QUESUP(IPMODL,IPMAT2,3,0,ISUP,ISUP2)
  60. c IF(ISUP.GT.1) THEN
  61. c CALL ERREUR(???)
  62. c RETURN
  63. c ENDIF
  64. CALL VLOC2(IPMODL,IPMAT2,IPCHE,IRET)
  65. C
  66. C=====================================================================
  67. C ECRITURE DES MCHAML
  68. C=====================================================================
  69. IF(IRET.NE.0) THEN
  70. CALL ACTOBJ('MCHAML ',IPCHE,1)
  71. CALL ECROBJ('MCHAML ',IPCHE)
  72. ENDIF
  73.  
  74. END
  75.  
  76.  
  77.  

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