Télécharger sol-asym+rela-unil.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : sol-asym+rela-unil.dgibi
  2.  
  3. *** Options ...
  4.  
  5. opti dime 2 mode plan defo echo 1 ;
  6.  
  7. valrefx = 1. ;
  8. valrefy = 2. ;
  9.  
  10. *** Points ...
  11.  
  12. p1 = 0 0 ;
  13. ma1 = manu POI1 p1;
  14.  
  15. *** Champ test ...
  16.  
  17. chpotst = manu chpo ma1 2 UX valrefx UY valrefy NATURE
  18. DIFFUS ;
  19.  
  20. *** Rigidité ...
  21.  
  22. ri1 = manu RIGI TYPE rigidite ma1 (mots UX UY)
  23. (prog 2 0 1 3) QUEL ;
  24.  
  25. cl1 = bloq mini uy p1 ;
  26. fodi1 = depi cl1 (valrefy * 0.99) ;
  27. fodi2 = depi cl1 (valrefy * 1.01) ;
  28.  
  29. *** Force ...
  30.  
  31. * fo1 = forc FX 2 FY 7 p1 ;
  32. fo1 = ri1 * chpotst ;
  33.  
  34. *** Résolution ...
  35.  
  36. list ri1 ;
  37. de0 = reso ri1 fo1 ;
  38. valux = extr de0 UX p1 ;
  39. valuy = extr de0 UY p1 ;
  40. diffx = valux - valrefx ;
  41. diffy = valuy - valrefy ;
  42. list diffx ;
  43. list diffy ;
  44. si((abs diffx > 1.e-10) ou (abs diffy > 1.e-10)) ;
  45. erre 5 ;
  46. finsi ;
  47.  
  48. ritot = ri1 et cl1 ;
  49. fotot = fo1 et fodi1 ;
  50. de1 = reso ritot fotot ;
  51. valux = extr de1 UX p1 ;
  52. valuy = extr de1 UY p1 ;
  53. diffx = valux - valrefx ;
  54. diffy = valuy - valrefy ;
  55. list diffx ;
  56. list diffy ;
  57. si((abs diffx > 1.e-10) ou (abs diffy > 1.e-10)) ;
  58. erre 5 ;
  59. finsi ;
  60.  
  61. ritot = ri1 et cl1 ;
  62. fotot = fo1 et fodi2 ;
  63. de2 = reso ritot fotot ;
  64. valux = extr de2 UX p1 ;
  65. valuy = extr de2 UY p1 ;
  66. diffx = valux - valrefx ;
  67. diffy = valuy - (valrefy * 1.01) ;
  68. list diffx ;
  69. list diffy ;
  70. si((abs diffx > 1.e-10) ou (abs diffy > 1.e-7)) ;
  71. erre 5 ;
  72. finsi ;
  73.  
  74.  
  75. fin ;
  76.  
  77.  
  78.  
  79.  

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