Télécharger nonconforme.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : nonconforme.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. * utilisation de maillage non conforme crée par raff
  6. *
  7. * mettre igraph='O' ; pour avoir les dessins
  8. *
  9. igraph = mot 'N';
  10. opti elem qua4 mode plan defo dime 2;
  11. dens 2.;
  12. * mesh 10x6
  13. pa= 0 0; pb= 10 0;pc= 20 0;
  14. liab= pa droi pb;libc= pb droi pc;
  15. su = (liab et libc) trans ( 0 12);
  16. si ( ega igraph 'O'); trac su; finsi;
  17. * definition of density
  18. x y = coor su;
  19. distance = ((x - 10 ) * ( x- 10) + ( y * y)) ** 0.5;
  20. den = 0.3 + (0.18*distance);
  21. si ( ega igraph 'O'); trac su den; finsi;
  22. * new mesh
  23. su2= raff su den;
  24. hh = elem su2 SURE ;
  25. sureal = su2 diff hh ;
  26. * use of this mesh
  27. * definition of model and caracteristic
  28. mo= model su2 mecanique elastique isotrope ;
  29. ma = MATE mo YOUN 2.e5 NU 0.3 ;
  30. * loads
  31. psupe = su2 point droite ( 0 12) ( 10 12) 0.1;
  32. lisupe = elem ( contour su2) appu stric psupe;
  33. ff = pres ( redu mo sureal) massif lisupe -1.;
  34. * displacements conditions
  35. py0= point su2 droit pa pc 0.01;
  36. liy0= elem ( contou su2) appu strict py0;
  37. li2bc = liy0 elem compris pb pc;
  38. cl1= bloqu li2bc UY;
  39. cl2= bloq UX pb;
  40. cltot= cl1 et cl2;
  41. * compute elastic solution
  42. ri = rigi mo ma;
  43. displa = reso ( ri et cltot) ff;
  44. stre = sigma displa mo ma;
  45. vm = vmis stre mo ma;
  46. si ( ega igraph 'O'); trac su2 vm mo ma; finsi;
  47. * compute stress intensity factor
  48. gt = table;
  49. gt.'OBJECTIF' = MOT 'J';
  50. lifis = liy0 elem compris pa pb;
  51. gt.'LEVRE_SUPERIEURE' = lifis;
  52. gt.'FRONT_FISSURE' = Pb;
  53. gt.'CARACTERISTIQUES' = ma;
  54. gt.'MODELE' = mo;
  55. gt.'SOLUTION_RESO' = displa;
  56. rea = reaction ri displa;
  57. gt.'CHARGEMENTS_MECANIQUES'=rea;
  58. naa = 5;
  59. gmoy=0.;
  60. repe no naa;
  61. gt.'COUCHE' = &no;
  62. G_THETA gt;
  63. si ( &no ega 1) ; g2=prog gt.resultats;sinon;
  64. g2 = g2 et ( prog gt.resultats); finsi;
  65. gmoy=gmoy + gt.resultats;
  66. fin no;
  67. gmoy = gmoy / naa;
  68. irate=0;
  69. repe noo naa;
  70. g1= extr g2 &noo;
  71. si (abs ( ( g1 - gmoy) / gmoy) > 0.03); IRATE=1;FINSI;
  72. fin noo;
  73. xx = prog 1 pas 1 naa;
  74. ev= evol manu 'nb of rows' xx 'G ' g2;
  75. ttt=table;
  76. ttt.1 = mot 'MARQ CROI';
  77. tt2=table;
  78. tt2.1= ' G ';
  79. ttt.'TITRE'=tt2;
  80. si ( ega igraph 'O'); dess ev lege ttt; finsi;
  81. * use of the refined mesh for thermal computation
  82. mmt = model su2 thermique isotrope;
  83. mmm= mate mmt k 2.5 'C' 2.7;
  84. cod = condu mmt mmm;
  85. pc3= poin su2 droite pc ( pc plus ( 0 1));
  86. pc1= poin su2 droit pa ( pa plus ( 0 1 ));
  87. bl1 = bloqu T pc1; bl2 = bloq T pc3;ft=depi bl2 12.;
  88. te= resou ( cod et bl1 et bl2 ) ft;
  89. si ( ega igraph 'O'); trac te su2;finsi;
  90. tpb = redu te pb;tpbv= maxi tpb;
  91. si(abs ( ( tpbv - 6. ) / 6. ) > 0.001); irate=1;finsi;
  92. *
  93. * resultats
  94. *
  95. mess ' valeur moyenne de G ' gmoy;
  96. mess ' valeurs trouvées ' ; list g2;
  97. mess ' valeur de la temperature au point PB ' tpbv;
  98. mess ' valeur attendue au point pb est 6.';
  99. ** trois points non conforme reliés par une relation 42 53
  100. ** et 95 au milieu
  101.  
  102. p42 = su2 point proc (16 6); t42= redu te p42; t42 = maxi t42;
  103. p53 = su2 point proc (16 8); t53 = redu te p53 ; t53 = maxi t53;
  104. p95 = su2 point proc (16 7); t95 = redu te p95 ; t95 = maxi t95;
  105. tmo= t42 + t53 / 2.; eca = tmo - t95 ;
  106. er_rela= abs (eca / t95);
  107. si (er_rela > 1.e-8);
  108. mess ' t42 t53 t95 ' t4é t53 t95 ; erreur 5;
  109. finsi;
  110.  
  111.  
  112. si (irate ega 1);
  113. erreur 5;
  114. finsi;
  115. fin;
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  

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