Télécharger nloc1.dgibi

Retour à la liste

Numérotation des lignes :

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

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