Télécharger nloc2.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : nloc2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. graph='N';
  5. saut page;
  6. mess 'test CONN et NLOC';
  7. *********************************************************************
  8. *
  9. * construction de connectivites sur des domaines differents
  10. * mesh12=1/4 de cylindre, mesh1 et mesh2=1/2 cylindre et
  11. * mesh=cylindre complet, avec ou sans symetrie de tel facon
  12. * que la reduction sur mesh12 d'un calcul non local mene a l'aide
  13. * de ces 4 conectivites soit toujours identique
  14. *
  15. *********************************************************************
  16. *
  17. opti dime 3 mode plan defo;
  18. *
  19. * 1) maillage: mesh est un maillage bi-symetrique
  20. * mesh12 en est le 1/4 (compose de me1[pri6] et me2[cub8])
  21. * mesh1 et mesh2 2 en sont des moities contenant mesh12
  22. *
  23. p1=0 0 0; p2=1 0 0; p22=.01 0 0; p3=0 1 0; p33=0 .01 0;
  24. pz=0 0 1000;
  25. *
  26. opti elem pri6;
  27. c1=(c 3 p22 p1 p33);
  28. nnc1=nbno c1;
  29. j=1; pci=c1 poin j;
  30. repeter lab1 (nnc1-1);
  31. j=j+1;
  32. pcf=c1 poin j;
  33. mecj=surf ((p1 d 1 pci) et (pci d 1 pcf) et (pcf d 1 p1)) 'PLANE';
  34. si (j ega 2); me1=mecj;
  35. sinon ; me1=me1 et mecj;
  36. finsi;
  37. pci=pcf;
  38. fin lab1;
  39. me1=me1 volu 1 tran pz;
  40. *
  41. opti elem cub8;
  42. me2=dall c1 (p33 d 3 p3) (c 3 p3 p1 p2) (p2 d 3 p22) 'PLAN';
  43. me2=me2 volu 1 tran pz;
  44. *
  45. *
  46. mesh12=me1 et me2; elim mesh12 1.d-4;
  47. *
  48. mesh2=mesh12 syme 'PLAN' p1 p2 pz;
  49. mesh3=mesh12 syme 'PLAN' p1 p3 pz;
  50. mesh4=mesh2 syme 'PLAN' p1 p3 pz;
  51. *
  52. mesh=mesh12 et mesh2 et mesh3 et mesh4;
  53. elim mesh 1.d-5;
  54. *
  55. mmesh1=mesh12 et mesh2; mesh1=mmesh1;
  56. mmesh2=mesh12 et mesh3; mesh2=mmesh2;
  57. *
  58. * 2) modele (tous les mesh ont mesh12 en commun)
  59. *
  60. mo=mode mesh mecanique consti constia;
  61. mo12=mode mesh12 mecanique consti constia;
  62. mo1=mode mesh1 mecanique consti constia;
  63. mo2=mode mesh2 mecanique consti constia;
  64. *
  65. * 3) listmot des composantes a moyenner et longueur caracteristique
  66. *
  67. cc=mots 'SCAL';
  68. lc=.4;
  69. *
  70. * 4) calcul sur le cylindre complet
  71. *
  72. mess 'calcul sur le cylindre complet';
  73. *
  74. * 4.1) connectivite (normale)
  75. *
  76. co =conn mo lc 'NORMAL' 'NO-MESH';
  77. *
  78. * 4.2) chamelem a moyenner
  79. * (construit par projection d'un chpoint nul partout sauf au centre)
  80. *
  81. ssca=(manu 'CHPO' p1 1 'SCAL' 1.)
  82. + (manu 'CHPO' mesh 1 'SCAL' 0.);
  83. csca=chan 'CHAM' ssca mo;
  84. ssca=chan 'STRESSES' mo csca;
  85. *
  86. * 4.3) moyenne non-locale
  87. *
  88. moy=nloc ssca co cc;
  89. *
  90. * 4.4) reduction sur le 1/4
  91. *
  92. mmoy=(redu moy me1) et (redu moy me2);
  93. *
  94. * 5) calcul sur le 1/4 de cylindre
  95. *
  96. mess 'calcul sur le 1/4 de cylindre';
  97. *
  98. * 5.1) connectivite avec symetries par rapport a 2 axes obtenues
  99. * de facon "orthogonale" a partir d'une symetrie point et de
  100. * 2 symetries droite
  101. *
  102. co12ecn=conn mo12 lc 'NORMAL' 'NO-MES12';
  103. co12ecv=conn mo12 lc 'PLAN' p1 p2 pz 'HO-MES12';
  104. co12ech=conn mo12 lc 'PLAN' p1 p3 pz 'VE-MES12';
  105. co12ecp=conn mo12 lc 'DROITE' p1 pz 'PT-MES12';
  106. co12vhp=co12ecn et co12ecv et co12ech et co12ecp;
  107. *
  108. * 5.2) chamelem a moyenner
  109. * (construit par projection d'un chpoint nul partout sauf au centre)
  110. *
  111. ssca12=(manu 'CHPO' p1 1 'SCAL' 1.)
  112. + (manu 'CHPO' mesh12 1 'SCAL' 0.);
  113. csca12=chan 'CHAM' ssca12 mo12;
  114. ssca12=chan 'STRESSES' mo12 csca12;
  115. *
  116. * 5.3) moyenne non-locale
  117. *
  118. moy12=nloc ssca12 co12vhp cc;
  119. *
  120. * 6) calcul sur le 1/2 cylindre vertical
  121. *
  122. mess 'calcul sur le 1/2 cylindre vertical';
  123. *
  124. * 6.1) connectivite avec une symetrie droite
  125. *
  126. co1ecn=conn mo1 lc 'NORMAL' 'NO-MESH1';
  127. co1ecv=conn mo1 lc 'PLAN' p1 p3 pz 'VE-MESH1';
  128. co1v=co1ecn et co1ecv;
  129. *
  130. * 6.2) chamelem a moyenner
  131. * (construit par projection d'un chpoint nul partout sauf au centre)
  132. *
  133. ssca1=(manu 'CHPO' p1 1 'SCAL' 1.)
  134. + (manu 'CHPO' mesh1 1 'SCAL' 0.);
  135. csca1=chan 'CHAM' ssca1 mo1;
  136. ssca1=chan 'STRESSES' mo1 csca1;
  137. *
  138. * 6.3) moyenne non-locale
  139. *
  140. moy1=nloc ssca1 co1v cc;
  141. *
  142. * 6.4) reduction sur le 1/4
  143. *
  144. mmoy1=(redu moy1 me1) et (redu moy1 me2);
  145. *
  146. *
  147. * 7) calcul sur le 1/2 cylindre horizontal
  148. *
  149. mess 'calcul sur le 1/2 cylindre horizontal';
  150. *
  151. * 7.1) connectivite avec une symetrie droite
  152. *
  153. co2ecn=conn mo2 lc 'NORMAL' 'NO-MESH2';
  154. co2ech=conn mo2 lc 'PLAN' p1 p2 pz 'VE-MESH2';
  155. co2h=co2ecn et co2ech;
  156. *
  157. * 7.2) chamelem a moyenner
  158. * (construit par projection d'un chpoint nul partout sauf au centre)
  159. *
  160. ssca2=(manu 'CHPO' p1 1 'SCAL' 1.)
  161. + (manu 'CHPO' mesh2 1 'SCAL' 0.);
  162. csca2=chan 'CHAM' ssca2 mo2;
  163. ssca2=chan 'STRESSES' mo2 csca2;
  164. *
  165. * 7.3) moyenne non-locale
  166. *
  167. moy2=nloc ssca2 co2h cc;
  168. *
  169. * 7.4) reduction sur le 1/4
  170. *
  171. mmoy2=(redu moy2 me1) et (redu moy2 me2);
  172. *
  173. * 8) erreur
  174. *
  175. m12m =(abs (moy12 - mmoy )) masque 'SUPERIEUR' 'SOMME' 1.d-10;
  176. m12m1=(abs (moy12 - mmoy1)) masque 'SUPERIEUR' 'SOMME' 1.d-10;
  177. m12m2=(abs (moy12 - mmoy2)) masque 'SUPERIEUR' 'SOMME' 1.d-10;
  178. lerr=m12m + m12m1 + m12m2;
  179. si (lerr > 0); erre 5;
  180. sinon; erre 0; finsi;
  181. *
  182. fin;
  183.  
  184.  
  185.  
  186.  
  187.  
  188.  
  189.  
  190.  

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