Télécharger la plaquette Cast3M


fr en

Télécharger normalisation-1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  normalisation-1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. *** Options ...
  6.  
  7.     opti dime 2 mode plan defo echo 1 ;
  8.  
  9.     critprec = 1.e-14 ;
  10.  
  11. *** Points ...
  12.  
  13.     dens 1. ;
  14.  
  15.     p1 = 0 0 ;
  16.  
  17.  
  18. *** Raideurs ...
  19.  
  20.     ri1 = manu rigi (manu poi1 p1) (mots UX UY) (prog 1 2
  21.                                                       2 3 ) ;
  22. *** Force ...
  23.  
  24.     fo1 = forc FX 1 FY 1 p1 ;
  25.  
  26. *** Première résolution ...
  27.  
  28.     de1 = reso ri1 fo1 ;
  29.  
  30. *** Deuxième, avec norm auto ...
  31.  
  32.     opti norm AUTO ;
  33.  
  34.     ri2 = 1 * ri1 ;
  35.     de2 = reso ri2 fo1 ;
  36.     dd = abs (de2 - de1) ;
  37.     si((maxi dd) > critprec) ;
  38.        erre 5 ;
  39.     finsi ;
  40.  
  41. *** Troisième, avec normalisation des primales ...
  42.  
  43.     opti norm (mots UX UY) (prog 2 3) ;
  44.     ri3 = 1 * ri1 ;
  45.     de3 = reso ri3 fo1 ;
  46.     dd = abs (de3 - de1) ;
  47.     si((maxi dd) > critprec) ;
  48.        erre 5 ;
  49.     finsi ;
  50.  
  51. *** Quatrième, avec normalisation des primales et des duales ...
  52.  
  53.     opti norm (mots UX UY) (prog 2 3) 
  54.               (mots FY FX) (prog 5 4) ;
  55.  
  56.     ri4 = 1 * ri1 ;
  57.     de4 = reso ri4 fo1 ;
  58.     dd = abs (de4 - de1) ;
  59.     si((maxi dd) > critprec) ;
  60.        erre 5 ;
  61.     finsi ;
  62.  
  63. *** Puis, le cas non-symétrique ...
  64.  
  65.     ra1 = manu rigi (manu poi1 p1) (mots UX UY)
  66.          (prog 1 2 3 4) QUEL ;
  67.  
  68.     fa1 = forc FX 3 FY 7 p1 ;
  69.  
  70. *** Première résolution ...
  71.  
  73.  
  74.     da1 = reso ra1 fa1 ;
  75.  
  76. *** On essaye avec la normalisation ...
  77.  
  78.     opti norm (mots UX UY) (prog 2 3) 
  79.               (mots FY FX) (prog 5 4) ;
  80.  
  81.     ra2 = 1 * ra1 ;
  82.     da2 = reso ra2 fa1 ;
  83.     dd = abs (da2 - da1) ;
  84.     si((maxi dd) > critprec) ;
  85.        erre 5 ;
  86.     finsi ;
  87.  
  88. *** Bye ...
  89.  
  90.     fin ;
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  

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