Télécharger la plaquette Cast3M


fr en

Télécharger proi.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  test_proi.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ******************************************************************
  5. * Projection de champs                                           *
  6. *    Syntaxe :                                                   *
  7. *    1) Champ par point aux noeuds du maillage GEO2              *
  8. *    2) Champ par element aux noeuds ou pt d'integration de MOD2 *
  9. *    3) Projection au sens des moindres carrés                   *
  10. *    4) Projection d'un champ par point et calcul de derivees    *
  11. *   Note : ce test verifie la fonctionnalite, pas le resultat    *
  12. *                                                                *
  13. * G. Folzan, Octobre 2013                                        *
  14. *                                                                *
  15. ******************************************************************
  16.  
  17. 'OPTI' 'DIME' 2 'ELEM' 'QUA4' ;
  18. *                                                                    *
  19.  
  20. * Creation des points
  21. PT_A       = 0   0 ;
  22. PT_B       = 0   1 ;
  23. PT_G       = 0   0 ;
  24. PT_H       = 0   1 ;
  25. PT_C       = 1   1 ;
  26. PT_D       = 1   0 ;
  27. PT_E       = 0.5 0 ;
  28. PT_F       = 0.5 1 ;
  29.  
  30.  
  31. *Creation des cotes du maillage                                       *
  32. SEG_AB       = 'DROI' 3 PT_A PT_B ;
  33. SEG_GH       = 'DROI' 3 PT_G PT_H ;
  34. SEG_BF       = 'DROI' 3 PT_B PT_F ;
  35. SEG_FC       = 'DROI' 3 PT_F PT_C ;
  36. SEG_CD       = 'DROI' 3 PT_C PT_D ;
  37. SEG_DE       = 'DROI' 3 PT_D PT_E ;
  38. SEG_EA       = 'DROI' 3 PT_E PT_A ;
  39. SEG_FE       = 'DROI' 3 PT_F PT_E ;
  40.  
  41.  
  42. *Creation des contours
  43. CONT_1       = SEG_AB et SEG_BF et SEG_FE et SEG_EA;
  44. CONT_2       = SEG_FC et SEG_CD et SEG_DE et SEG_FE;
  45.  
  46.  
  47. *Creation  de maillages differents
  48. *     1) 10 couches d'elements
  49. *     2) 13 couches d'elements
  50. *     3) non structure: quadrangle et triangle
  51.  
  52. MAIL_1       = SEG_AB 'TRAN' 10 PT_D ;
  53.  
  54. MAIL_2       = (SEG_GH 'TRAN' 13 PT_D) 'COUL' 'VERT' ;
  55.  
  56. mail3a       = surf cont_1;
  57. 'OPTI' 'DIME' 2 'ELEM' 'TRI3' ;
  58. mail3b       = surf cont_2;
  59. MAIL_3       = (MAIL3A et MAIL3B) coul 'ROUGE';
  60.  
  61. *                                                                      *
  62. *'TRAC' (MAIL_1 'ET' MAIL_2 'ET' MAIL_3) ;
  63. *                                                                      *
  64.  
  65. *Creation de modeles associes a chaque maillage
  66.  
  67. MOD1     = 'MODE' MAIL_1 'THERMIQUE' ;
  68. MOD2     = 'MODE' MAIL_2 'MECANIQUE' ;
  69. MOD3     = 'MODE' MAIL_3 'MECANIQUE' ;
  70.  
  71.  
  72. *Creation du champ de depart: champ par element defini au noeud
  73. CHPO1    = 'MANU' 'CHPO' MAIL_1 1 'T' 100. ;
  74. CHAM1    = 'CHAN' 'CHAM' CHPO1 MOD1 ;
  75. *                                                              
  76.  
  77. *Test syntaxe 1
  78. CHPO2 = 'PROI' CHAM1 MAIL_2;
  79. CHPO2 = 'PROI' MAIL_2 CHAM1;
  80. CHPO3 = 'PROI' CHAM1 MAIL_3;
  81. 'MENAGE' 'OBLI';
  82.  
  83. *Test syntaxe 2
  84. CHAM2    = 'PROI' MOD2 CHAM1;
  85. CHAM2    = 'PROI' CHAM1 MOD2;
  86. CHAM2    = 'PROI' 'GRAVITE' MOD2 CHAM1 ;
  87. CHAM2    = 'PROI' MOD2 'RIGIDITE' CHAM1 ;
  88. CHAM2    = 'PROI' CHAM1 MOD2 'MASSES' ;
  89. CHAM2    = 'PROI' CHAM1 'STRESSES' MOD2;
  90. CHAM3    = 'PROI' MOD3 CHAM1;
  91. CHAM3    = 'PROI' MOD3 CHAM1 'GRAVITE' ;
  92. CHAM3    = 'PROI' MOD3 CHAM1 'RIGIDITE' ;
  93. CHAM3    = 'PROI' MOD3 CHAM1 'MASSES' ;
  94. CHAM3    = 'PROI' MOD3 CHAM1 'STRESSES' ;
  95. 'MENAGE' 'OBLI';
  96.  
  97. *Test syntaxe 3
  98. CHAM2    = 'PROI' 'MINI' MOD2 CHAM1;
  99. CHAM2    = 'PROI' 'MINI' CHAM1 MOD2;
  100. CHAM2    = 'PROI' 'MINI' MOD2 CHAM1 5;
  101. CHAM3    = 'PROI' MOD3 CHAM1 'MINI' ;
  102. CHAM3    = 'PROI' 5 MOD3 'MINI' CHAM1 ;
  103. 'MENAGE' 'OBLI';
  104.  
  105. *Test syntaxe 4
  106. CHPO2 = PROI 'POLY' MAIL_1 MAIL_2 CHPO1 1 'PLAN';
  107. CHPO2 = PROI 'POLY' CHPO1 MAIL_1 MAIL_2 1 'PLAN';
  108. CHPO2 = PROI 'POLY' MAIL_1 MAIL_2 CHPO1 1 'PLAN' 'POID' 5 9;
  109. CHPO2 = PROI 'POLY' MAIL_1 CHPO1 MAIL_2 1 'AXIS';
  110. 'MENAGE' 'OBLI';
  111.  
  112. *  TEST DE PROJECTION SUR DES ELEMENTS COQUES                          *
  113. OPTI DIME 3 elem cub8;
  114.  
  115. *Creation du maillage et champ de depart (element massif)
  116. vec1 = 0. 0. 0.1;
  117. vol1 = mail_1 VOLU 1 'TRAN' vec1;
  118. vec1 = 0. 0. -0.05;
  119. depl vol1 'PLUS' vec1;
  120. mod1b =  'MODE' vol1 'THERMIQUE' ;
  121. CHPO1    = 'MANU' 'CHPO' vol1 1 'T' 100. ;
  122. CHAM1    = 'CHAN' 'CHAM' CHPO1 MOD1b ;
  123.  
  124. *Creation du maillage de reception (element coque)
  125. MOD2b    = 'MODE' MAIL_2 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 'COQ4';
  126. CAR2 = CARA MOD2b EPAI 0.1 ;
  127.  
  128. *Test syntaxe 2
  129. CHAM2    = 'PROI' CAR2 MOD2b CHAM1 ;
  130. CHAM2    = 'PROI' MOD2b CAR2 CHAM1 ;
  131. CHAM2    = 'PROI' CHAM1 MOD2b CAR2 ;
  132. CHAM2    = 'PROI' CAR2 MOD2b CHAM1 'GRAVITE' ;
  133. CHAM2    = 'PROI' CAR2 MOD2b CHAM1 'RIGIDITE' ;
  134. CHAM2    = 'PROI' CAR2 MOD2b CHAM1 'STRESSES' ;
  135. CHAM2    = 'PROI' CAR2 MOD2b CHAM1 'STRESSES' ;
  136. 'MENAGE' 'OBLI';
  137.  
  138. *Test syntaxe 3
  139. CHAM2    = 'PROI' 'MINI' CAR2 MOD2b CHAM1 ;
  140. CHAM2    = 'PROI' 'MINI' MOD2b CHAM1 CAR2 ;
  141. CHAM2    = 'PROI' 'MINI' CHAM1 CAR2 MOD2b ;
  142. CHAM2    = 'PROI' 'MINI' CAR2 MOD2b CHAM1 5; 
  143. 'MENAGE' 'OBLI';
  144.  
  145. fin;
  146.  
  147.  
  148.  
  149.  
  150.  
  151.  
  152.  

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