Télécharger vloc.eso

Retour à la liste

Numérotation des lignes :

vloc
  1. C VLOC SOURCE CB215821 23/01/25 21:15:39 11573
  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. -INC SMCOORD
  30.  
  31. IPMODL=0
  32. IPCHE1=0
  33. IPCHE2=0
  34. IPCHE3=0
  35. IRET=0
  36. C
  37. C=====================================================================
  38. C LECTURE DU MODELE
  39. C=====================================================================
  40. C
  41. CALL LIROBJ('MMODEL ',IPMODL,1,IRET1)
  42. CALL ACTOBJ('MMODEL ',IPMODL,1)
  43. IF (IERR.NE.0) RETURN
  44. C
  45. C=====================================================================
  46. C LECTURE DU MATERIAU
  47. C=====================================================================
  48. C
  49. CALL LIROBJ('MCHAML ',IPMATE,1,IRET2)
  50. CALL ACTOBJ('MCHAML ',IPMATE,1)
  51. IF (IERR.NE.0) RETURN
  52. C
  53. C=====================================================================
  54. C CALCUL DES VECTEURS :
  55. C=====================================================================
  56. C DU REPERE LOCAL D'ORTHOTROPIE
  57. CALL REDUAF(IPMATE,IPMODL,IPMAT2,0,IRET3,KERR)
  58. IF (IRET3.NE.1) CALL ERREUR(KERR)
  59. IF(IERR.NE.0) RETURN
  60. c CALL QUESUP(IPMODL,IPMAT2,3,0,ISUP,ISUP2)
  61. c IF(ISUP.GT.1) THEN
  62. c CALL ERREUR(???)
  63. c RETURN
  64. c ENDIF
  65. SEGACT,MCOORD
  66. CALL VLOC2(IPMODL,IPMAT2,IPCHE,IRET)
  67. SEGDES,MCOORD
  68. C
  69. C=====================================================================
  70. C ECRITURE DES MCHAML
  71. C=====================================================================
  72. IF(IRET.NE.0) THEN
  73. CALL ACTOBJ('MCHAML ',IPCHE,1)
  74. CALL ECROBJ('MCHAML ',IPCHE)
  75. ENDIF
  76.  
  77. END
  78.  
  79.  
  80.  
  81.  

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