Télécharger effmarti.procedur

Retour à la liste

Numérotation des lignes :

  1. * EFFMARTI PROCEDUR AF221230 16/11/03 21:15:02 9139
  2. DEBPROC EFFMARTI SIG1*'MCHAML' MOD1*'MMODEL' MAT1*'MCHAML'
  3. VECT1*'POINT' VECT2*'POINT'
  4. H1*'FLOTTANT' ENRE1*'FLOTTANT' ENRI1*'FLOTTANT'
  5. COT1*'FLOTTANT';
  6. *
  7. * --------------------------------------------------------
  8. * Procedure EFF_MARTI
  9. * Definition des efforts de membrane pour la couche
  10. * superieure et inferieure pour un element coque
  11. *
  12. * Develeppeurs:
  13. * Alberto FRAU
  14. * DEN/DANS/DM2S/SEMT/EMSI
  15. * Nicolas ILE
  16. * DEN/DANS/DM2S/SEMT/EMSI
  17. *
  18. * ENTREE:
  19. * SIG1: MCHAML des contraintes
  20. * MOD1: MMODEL associé au SIG1
  21. * MAT1: MCHAML associé au SIG1
  22. * VEC1: direction 1 pour le calcul des contraintes
  23. * VEC2: direction 2 pour le calcul des contraintes
  24. * H1: epaisseur de l'element plaque
  25. * ENRE1: enrobage couche externe
  26. * ENRi1: enrobage couche interne
  27. * COT1: terme cotg(theta)
  28. *
  29. * SORTIE:
  30. * SIG2: MCHAML avec les composants des efforts de membrane
  31. * pour la couche inferieure et superieure et les
  32. * efforts de tranchants hors plan
  33. * --------------------------------------------------------
  34. *
  35. * controle dimension du calcul
  36. 'SI' (('VALE' 'DIME') 'NEG' 3);
  37. 'MESS' 'Calcul possible seulement en 3D';
  38. QUIT EFFMARTI;
  39. 'FINSI';
  40. *
  41. * Test sur l'orientation du maillage
  42. VERSSENS ('EXTR' SIG1 'MAIL');
  43. *
  44. * Changement des directions
  45. SIG1 = RTENS SIG1 MOD1 MAT1 VECT1 VECT2;
  46. *
  47. * changement à GRAVITE
  48. MAI1 = EXTR SIG1 MAIL;
  49. SIG1B = 'CHAN' 'GRAVITE' MOD1 SIG1;
  50. SIG1B = REDU SIG1B MAI1;
  51. *
  52. *
  53. * Creaction des MCHAML vides
  54. N11E = 'VIDE' 'MCHAML';
  55. N12E = 'VIDE' 'MCHAML';
  56. N22E = 'VIDE' 'MCHAML';
  57. *
  58. N11I = 'VIDE' 'MCHAML';
  59. N12I = 'VIDE' 'MCHAML';
  60. N22I = 'VIDE' 'MCHAML';
  61. *
  62. M11 = 'VIDE' 'MCHAML';
  63. M22 = 'VIDE' 'MCHAML';
  64. M12 = 'VIDE' 'MCHAML';
  65. *
  66. V1 = 'VIDE' 'MCHAML';
  67. V2 = 'VIDE' 'MCHAML';
  68. VR = 'VIDE' 'MCHAML';
  69. *
  70. * extraction des efforts globaux
  71. M11 = M11 'ET' (EXCO SIG1B 'M11' 'M11T');
  72. M22 = M22 'ET' (EXCO SIG1B 'M22' 'M22T');
  73. M12 = M12 'ET' (EXCO SIG1B 'M12' 'M12T');
  74. V1 = V1 'ET' (EXCO SIG1B 'V1' 'V1T');
  75. V2 = V2 'ET' (EXCO SIG1B 'V2' 'V2T');
  76. *
  77. CHP1 = 'EXCO' V1 'V1T' 'SCAL';
  78. CHP2 = 'EXCO' V2 'V2T' 'SCAL';
  79. CHP3 = ((CHP1)*(CHP1)) + ((CHP2)*(CHP2));
  80. CHP3 = ((CHP3)**(0.5));
  81. AAA = MASQU CHP3 'EGINFE' (1.E-12);
  82. CHP3 = CHP3 + (((AAA))*(1.E99));
  83. VR = VR 'ET' (EXCO CHP3 'SCAL' 'VR');
  84. 'OUBL' CHP1;
  85. 'OUBL' CHP2;
  86. 'OUBL' CHP3;
  87. *
  88. * Determination des efforts de membrane pour les couches
  89. * inferieure et superieure
  90. N11S = 'EXCO' SIG1B 'N11' 'SCAL';
  91. N12S = 'EXCO' SIG1B 'N12' 'SCAL';
  92. N22S = 'EXCO' SIG1B 'N22' 'SCAL';
  93. M11S = 'EXCO' M11 'M11T' 'SCAL';
  94. M22S = 'EXCO' M22 'M22T' 'SCAL';
  95. M12S = 'EXCO' M12 'M12T' 'SCAL';
  96. V1S = 'EXCO' V1 'V1T' 'SCAL';
  97. V2S = 'EXCO' V2 'V2T' 'SCAL';
  98. V0S = 'EXCO' VR 'VR' 'SCAL';
  99. *
  100. DD1 = H1 - ((ENRE1 + ENRI1)) ;
  101. *
  102. N11E = N11E 'ET' (((0.5)*(N11S)) +
  103. ((M11S)/(DD1)) +
  104. (((V1S)*(V1S))/(((2.0)*(COT1))*(V0S))));
  105. N22E = N22E 'ET' (((0.5)*(N22S)) +
  106. ((M22S)/(DD1)) +
  107. (((V2S)*(V2S))/(((2.0)*(COT1))*(V0S))));
  108. N12E = N12E 'ET' (((0.5)*(N12S)) +
  109. ((M12S)/(DD1)) +
  110. (((V1S)*(V2S))/(((2.0)*(COT1))*(V0S))));
  111. N11I = N11I 'ET' (((0.5)*(N11S)) -
  112. ((M11S)/(DD1)) +
  113. (((V1S)*(V1S))/(((2.0)*(COT1))*(V0S))));
  114. N22I = N22I 'ET' (((0.5)*(N22S)) -
  115. ((M22S)/(DD1)) +
  116. (((V2S)*(V2S))/(((2.0)*(COT1))*(V0S))));
  117. N12I = N12I 'ET' (((0.5)*(N12S)) -
  118. ((M12S)/(DD1)) +
  119. (((V1S)*(V2S))/(((2.0)*(COT1))*(V0S))));
  120. *
  121. N11E = 'EXCO' N11E 'SCAL' 'N11E';
  122. N22E = 'EXCO' N22E 'SCAL' 'N22E';
  123. N12E = 'EXCO' N12E 'SCAL' 'N12E';
  124. N11I = 'EXCO' N11I 'SCAL' 'N11I';
  125. N22I = 'EXCO' N22I 'SCAL' 'N22I';
  126. N12I = 'EXCO' N12I 'SCAL' 'N12I';
  127. *
  128. * Assemblage
  129. SIG2 = N11E ET N22E ET N12E ET
  130. N11I ET N22I ET N12I ET
  131. M11 ET M22 ET M12 ET
  132. V1 ET V2 ET VR;
  133. *
  134. FINPROC SIG2;
  135.  
  136.  

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