Télécharger rela.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : rela.dgibi
  2. 'OPTI' echo 0 ;
  3. ************************************************************************
  4. * NOM : RELA
  5. * DESCRIPTION : Cas-test elementaire RELA où un des maillages est
  6. * un point (nouvelle syntaxe du 2019/01)
  7. *
  8. *
  9. *
  10. * LANGAGE : GIBIANE-CAST3M
  11. * AUTEUR : Stephane GOUNAND (CEA/DEN/DM2S/SEMT/LTA)
  12. * mel : stephane.gounand@cea.fr
  13. **********************************************************************
  14. * VERSION : v1, 12/02/2019, version initiale
  15. * HISTORIQUE : v1, 12/02/2019, création
  16. * HISTORIQUE :
  17. * HISTORIQUE :
  18. ************************************************************************
  19. *
  20. interact = faux ;
  21. *
  22. 'OPTI' 'DIME' 2 ;
  23. lok = vrai ; tol = '*' ('VALE' 'PREC') 2. ;
  24. *
  25. p0 = 0. 0. ; p1 = 1. 0.; p2 = 1. 1. ;
  26. p3 = 2. 0. ; p4 = 2. 1. ;
  27. p5 = pi pi ;
  28. *
  29. m1 = 'MANU' 'SEG2' p2 p4 ;
  30. m2 = 'MANU' 'TRI3' p0 p1 p2 ;
  31. m3 = 'MANU' 'SEG2' p1 p3 ;
  32. *
  33. mt = m1 'ET' m2 'ET' m3 ;
  34. mtb = mt 'PLUS' (0. 2.) ;
  35. *
  36. * Test 1 : deux maillages complexes
  37. *
  38. alf = '**' 2. 0.5 ;
  39. bet = PI ;
  40. rel = 'RELA' 'T' mtb '-' 'T' mt ;
  41. ch = 'MANU' 'CHPO' mt 1 'T' alf ;
  42. matd = 'MANU' 'RIGI' ch 'DIAG' ;
  43. chf = 'MANU' 'CHPO' mt 1 'Q' bet ;
  44. solex = 'MANU' 'CHPO' mt 1 'T' ('/' bet alf) ;
  45. *
  46. mat = matd 'ET' rel ;
  47. smb = chf ;
  48. *
  49. sol = 'RESO' mat smb ;
  50. * On vérifie que les valeurs sur mt et mtb sont les memes egales à PI/2
  51. * et que les multiplicateurs de Lagrange sont nuls.
  52. chb = 'REDU' sol mtb ;
  53. err1 = 'MAXI' ('-' chb (bet '/' alf)) 'ABS' ;
  54. 'MESS' ('CHAI' 'err1=' err1 ' tol=' tol) ;
  55. 'SI' ('>' err1 tol) ;
  56. 'MESS' ('CHAI' '!!! err1=' err1 ' > tol=' tol) ;
  57. lok = lok 'ET' faux ;
  58. 'FINS' ;
  59. *
  60. chbb = 'DEDU' 'TRAN' chb mtb mt ;
  61. err2 = 'MAXI' ('-' chbb ('REDU' solex mt)) 'ABS' ;
  62. 'MESS' ('CHAI' 'err2=' err2 ' tol=' tol) ;
  63. 'SI' ('>' err2 tol) ;
  64. 'MESS' ('CHAI' '!!! err2=' err2 ' > tol=' tol) ;
  65. lok = lok 'ET' faux ;
  66. 'FINS' ;
  67. *
  68. err3 = 'MAXI' ('EXCO' 'LX' sol) 'ABS' ;
  69. 'MESS' ('CHAI' 'err3=' err3 ' tol=' tol) ;
  70. 'SI' ('>' err3 tol) ;
  71. 'MESS' ('CHAI' '!!! err3=' err3 ' > tol=' tol) ;
  72. lok = lok 'ET' faux ;
  73. 'FINS' ;
  74. *
  75. * Test 2 : un maillage complexe en relation avec un noeud
  76. *
  77. rel2 = 'RELA' 'T' mt '-' 'T' p5 ;
  78. smb2 = 'NOMC' 'Q' ('+' ('COOR' 1 mt) ('*' bet ('COOR' 2 mt))) ;
  79. mat2 = matd 'ET' rel2 ;
  80. sol2 = 'RESO' mat2 smb2 ;
  81. *
  82. vp5 = 'EXTR' sol2 'T' p5 ;
  83. vp5ex = ('MAXI' ('RESU' smb2)) '/' alf '/' ('NBNO' mt) ;
  84. err4 = 'ABS' ('-' vp5 vp5ex) ;
  85. 'MESS' ('CHAI' 'err4=' err4 ' tol=' tol) ;
  86. 'SI' ('>' err4 tol) ;
  87. 'MESS' ('CHAI' '!!! err4=' err4 ' > tol=' tol) ;
  88. lok = lok 'ET' faux ;
  89. 'FINS' ;
  90. *
  91. * Test 3 : un mix
  92. *
  93. rel3 = 'RELA' 'T' mt '-' 2. 'T' p5 '+' 'T' mtb ;
  94. mtt = mt 'ET' mtb ;
  95. smb3 = 'NOMC' 'Q' ('+' ('COOR' 1 mtt) ('*' bet ('COOR' 2 mtt))) ;
  96. cht = 'MANU' 'CHPO' mtt 1 'T' alf ;
  97. matdt = 'MANU' 'RIGI' cht 'DIAG' ;
  98. mat3 = matdt 'ET' rel3 ;
  99. sol3 = 'RESO' mat3 smb3 ;
  100. *
  101. vp5b = 'EXTR' sol3 'T' p5 ;
  102. vp5exb = ('MAXI' ('RESU' smb3)) '/' alf '/' ('NBNO' mtt) ;
  103. err5 = 'ABS' ('-' vp5b vp5exb) ;
  104. 'MESS' ('CHAI' 'err5=' err5 ' tol=' tol) ;
  105. 'SI' ('>' err5 tol) ;
  106. 'MESS' ('CHAI' '!!! err5=' err5 ' > tol=' tol) ;
  107. lok = lok 'ET' faux ;
  108. 'FINS' ;
  109. *
  110. 'SAUT' 1 'LIGNE' ;
  111. 'SI' lok ;
  112. 'MESSAGE' 'Tout sest bien passe' ;
  113. 'SINON' ;
  114. 'MESSAGE' '!!! Il y a eu des erreurs' ;
  115. 'FINSI' ;
  116. 'SAUT' 1 'LIGNE' ;
  117. *
  118. 'SI' interact ;
  119. 'OPTION' 'DONN' 5 'ECHO' 1 ;
  120. 'FINSI' ;
  121. 'SI' ('NON' lok) ;
  122. 'ERREUR' 5 ;
  123. 'FINSI' ;
  124. *
  125. * End of dgibi file RELA
  126. *
  127. 'FIN' ;
  128.  
  129.  
  130.  
  131.  

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