Télécharger @clpd.procedur

Retour à la liste

Numérotation des lignes :

  1. * @CLPD PROCEDUR BP208322 21/01/08 21:15:01 10844
  2. ****************************************************
  3. ****************************************************
  4. *
  5. * PROCEDURE POUR CONDITIONS AUX LIMITES
  6. *
  7. * PERIODIQUE EN DEFORMATION IMPOSEE
  8. *
  9. * L.GELEBART, G. TREGO 10/06
  10. ****************************************************
  11. *****************************************************
  12. DEBPROC @CLPD MAIL0*MAILLAGE DEF0*TABLE;
  13. EXX = DEF0 . 1; EYY = DEF0 . 2; EZZ = DEF0 . 3;
  14. EXY = DEF0 . 4; EXZ = DEF0 . 5; EYZ = DEF0 . 6;
  15. MOD0 = MODE MAIL0 MECANIQUE ELASTIQUE;
  16. *
  17. * Identification des points et surfaces exterieures
  18. * pour application des CL
  19. *=====================================================
  20. ENVE0 = ENVE MAIL0;
  21. elim0 = 1.e-4;
  22. ENVE_PT0 = CHANGE ENVE0 POI1;
  23. XMAX = MAXI (COOR 1 ENVE0);XMIN = MINI (COOR 1 ENVE0);
  24. YMAX = MAXI (COOR 2 ENVE0);YMIN = MINI (COOR 2 ENVE0);
  25. ZMAX = MAXI (COOR 3 ENVE0);ZMIN = MINI (COOR 3 ENVE0);
  26. LX0 = XMAX - XMIN; LY0 = YMAX - YMIN; LZ0 = ZMAX - ZMIN;
  27. PXMIN = ENVE_PT0 POIN 'PLAN'
  28. (XMIN 0. 0.) (XMIN 1. 0.) (XMIN 0. 1.) (LX0 * elim0);
  29. PYMIN = ENVE_PT0 POIN 'PLAN'
  30. (0. YMIN 0.) (1. YMIN 0.) (0. YMIN 1.) (LY0 * elim0);
  31. PZMIN = ENVE_PT0 POIN 'PLAN'
  32. (0. 0. ZMIN) (1. 0. ZMIN) (0. 1. ZMIN) (LZ0 * elim0);
  33. *
  34. SURFX0 = ENVE0 ELEM APPUYE STRICTEMENT PXMIN;
  35. SURFY0 = ENVE0 ELEM APPUYE STRICTEMENT PYMIN;
  36. SURFZ0 = ENVE0 ELEM APPUYE STRICTEMENT PZMIN;
  37. SURFX0 = SURFX0 CHAN POI1;
  38. SURFY0 = SURFY0 CHAN POI1;
  39. SURFZ0 = SURFZ0 CHAN POI1;
  40. SURFX1 = SURFX0 PLUS (LX0 0. 0.);
  41. **SURFY1 = SURFY0 PLUS (0. LY0 0.);
  42. **SURFZ1 = SURFZ0 PLUS (0. 0. LZ0);
  43. elim SURFX1 MAIL0 (LX0 * elim0);
  44. **elim SURFY1 MAIL0 (LY0 * elim0);
  45. **elim SURFZ1 MAIL0 (LZ0 * elim0);
  46. *trac (SURFX0 ET SURFY0 ET SURFZ0);
  47. *trac (SURFX1 ET SURFY1 ET SURFZ1);
  48. *
  49. A1 = ENVE_PT0 POIN PROC (XMIN YMIN ZMIN);
  50. *
  51. * Blocage du mvt de corps rigide
  52. *===========================================
  53. BL1 = BLOQUE 'UX' 'UY' 'UZ' A1;
  54. *BL1 = @CORIGI ENVE0;
  55. *
  56. * Application des conditions de periodicite
  57. *===========================================
  58. RIGCLP = BL1;
  59. *deplacements imposes*
  60. *faces perpendiculaires a X
  61. CLP1X = RELA 1. 'UX' SURFX1 - 1. 'UX' SURFX0;
  62. depCLP1X = DEPI CLP1X (EXX * LX0);
  63. CLP1Y = RELA 1. 'UY' SURFX1 - 1. 'UY' SURFX0;
  64. depCLP1Y = DEPI CLP1Y (EXY * LX0);
  65. CLP1Z = RELA 1. 'UZ' SURFX1 - 1. 'UZ' SURFX0;
  66. depCLP1Z = DEPI CLP1Z (EXZ * LX0);
  67. CLP1 = CLP1X ET (CLP1Y ET CLP1Z);
  68. depCLP1 = depCLP1X ET (depCLP1Y ET depCLP1Z);
  69.  
  70. *faces perpendiculaires a Y
  71. SURFY0 = SURFY0 DIFF (SURFY0 INTE SURFX0);
  72. SURFY1 = SURFY0 PLUS (0. LY0 0.);
  73. elim SURFY1 MAIL0 (LY0 * elim0);
  74.  
  75. CLP2X = RELA 1. 'UX' SURFY1 - 1. 'UX' SURFY0;
  76. depCLP2X = DEPI CLP2X (EXY * LY0);
  77. CLP2Y = RELA 1. 'UY' SURFY1 - 1. 'UY' SURFY0;
  78. depCLP2Y = DEPI CLP2Y (EYY * LY0);
  79. CLP2Z = RELA 1. 'UZ' SURFY1 - 1. 'UZ' SURFY0;
  80. depCLP2Z = DEPI CLP2Z (EYZ * LY0);
  81. CLP2 = CLP2X ET (CLP2Y ET CLP2Z);
  82. depCLP2 = depCLP2X ET (depCLP2Y ET depCLP2Z);
  83. *faces perpendiculaires Z
  84. SURFZ0 = SURFZ0 DIFF (SURFZ0 INTE SURFX0);
  85. SURFZ0 = SURFZ0 DIFF (SURFZ0 INTE SURFY0);
  86. SURFZ1 = SURFZ0 PLUS (0. 0. LZ0);
  87. elim SURFZ1 MAIL0 (LZ0 * elim0);
  88. CLP3X = RELA 1. 'UX' SURFZ1 - 1. 'UX' SURFZ0;
  89. depCLP3X = DEPI CLP3X (EXZ * LZ0);
  90. CLP3Y = RELA 1. 'UY' SURFZ1 - 1. 'UY' SURFZ0;
  91. depCLP3Y = DEPI CLP3Y (EYZ * LZ0);
  92. CLP3Z = RELA 1. 'UZ' SURFZ1 - 1. 'UZ' SURFZ0;
  93. depCLP3Z = DEPI CLP3Z (EZZ * LZ0);
  94. CLP3 = CLP3X ET (CLP3Y ET CLP3Z);
  95. depCLP3 = depCLP3X ET (depCLP3Y ET depCLP3Z);
  96. *
  97. DEPTOT = depCLP1 et (depCLP2 et depCLP3 );
  98. RIGCLP = RIGCLP et (CLP1 et (CLP2 et CLP3));
  99. *
  100. FINPROC RIGCLP DEPTOT;
  101. *********************************************
  102. * FIN DE LA PROCEDURE CLPD
  103. *********************************************
  104.  
  105.  
  106.  
  107.  

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