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. * defoinition of the conformity relations due to the refinement
  31. rel = RELA mo;
  32. * loads
  33. psupe = su2 point droite ( 0 12) ( 10 12) 0.1;
  34. lisupe = elem ( contour su2) appu stric psupe;
  35. ff = pres ( redu mo sureal) massif lisupe -1.;
  36. * displacements conditions
  37. py0= point su2 droit pa pc 0.01;
  38. liy0= elem ( contou su2) appu strict py0;
  39. li2bc = liy0 elem compris pb pc;
  40. cl1= bloqu li2bc UY;
  41. cl2= bloq UX pb;
  42. cltot= cl1 et cl2;
  43. * compute elastic solution
  44. ri = rigi mo ma;
  45. displa = reso ( ri et cltot et rel) ff;
  46. stre = sigma displa mo ma;
  47. vm = vmis stre mo ma;
  48. si ( ega igraph 'O'); trac su2 vm mo ma; finsi;
  49. * compute stress intensity factor
  50. gt = table;
  51. gt.'OBJECTIF' = MOT 'J';
  52. lifis = liy0 elem compris pa pb;
  53. gt.'LEVRE_SUPERIEURE' = lifis;
  54. gt.'FRONT_FISSURE' = Pb;
  55. gt.'CARACTERISTIQUES' = ma;
  56. gt.'MODELE' = mo;
  57. gt.'SOLUTION_RESO' = displa;
  58. rea = reaction (ri et rel) displa;
  59. gt.'CHARGEMENTS_MECANIQUES'=rea;
  60. * gt.'BLOCAGES_MECANIQUES' = (cltot et rel);
  61. naa = 5;
  62. gmoy=0.;
  63. repe no naa;
  64. gt.'COUCHE' = &no;
  65. G_THETA gt;
  66. si ( &no ega 1) ; g2=prog gt.resultats;sinon;
  67. g2 = g2 et ( prog gt.resultats); finsi;
  68. gmoy=gmoy + gt.resultats;
  69. fin no;
  70. gmoy = gmoy / naa;
  71. irate=0;
  72. repe noo naa;
  73. g1= extr g2 &noo;
  74. si (abs ( ( g1 - gmoy) / gmoy) > 0.03); IRATE=1;FINSI;
  75. fin noo;
  76. xx = prog 1 pas 1 naa;
  77. ev= evol manu 'nb of rows' xx 'G ' g2;
  78. ttt=table;
  79. ttt.1 = mot 'MARQ CROI';
  80. tt2=table;
  81. tt2.1= ' G ';
  82. ttt.'TITRE'=tt2;
  83. si ( ega igraph 'O'); dess ev lege ttt; finsi;
  84. * use of the refined mesh for thermal computation
  85. mmt = model su2 thermique isotrope;
  86. mmm= mate mmt k 2.5 'C' 2.7;
  87. cod = condu mmt mmm;
  88. pc3= poin su2 droite pc ( pc plus ( 0 1));
  89. pc1= poin su2 droit pa ( pa plus ( 0 1 ));
  90. bl1 = bloqu T pc1; bl2 = bloq T pc3;ft=depi bl2 12.;
  91. te= resou ( cod et bl1 et bl2 ) ft;
  92. si ( ega igraph 'O'); trac te su2;finsi;
  93. tpb = redu te pb;tpbv= maxi tpb;
  94. si(abs ( ( tpbv - 6. ) / 6. ) > 0.001); irate=1;finsi;
  95. *
  96. * resultats
  97. *
  98. mess ' valeur moyenne de G ' gmoy;
  99. mess ' valeurs trouvées ' ; list g2;
  100. mess ' valeur de la temperature au point PB ' tpbv;
  101. mess ' valeur attendue au point pb est 6.';
  102. ** trois points non conforme reliés par une relation 42 53
  103. ** et 95 au milieu
  104.  
  105. p42 = su2 point proc (16 6); t42= redu te p42; t42 = maxi t42;
  106. p53 = su2 point proc (16 8); t53 = redu te p53 ; t53 = maxi t53;
  107. p95 = su2 point proc (16 7); t95 = redu te p95 ; t95 = maxi t95;
  108. tmo= t42 + t53 / 2.; eca = tmo - t95 ;
  109. er_rela= abs (eca / t95);
  110. si (er_rela > 1.e-8);
  111. mess ' t42 t53 t95 ' t4é t53 t95 ; erreur 5;
  112. finsi;
  113.  
  114.  
  115. si (irate ega 1);
  116. erreur 5;
  117. finsi;
  118. fin;
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  

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