Datasets examples

* fichier :  sol-asym+rela-unil.dgibi

*** Options ...

    opti dime 2 mode plan defo echo 1 ;

    valrefx = 1. ;
    valrefy = 2. ;

*** Points ...

    p1 = 0 0 ;
    ma1 = manu POI1 p1;

*** Champ test ...

    chpotst = manu chpo ma1 2 UX valrefx UY valrefy NATURE
              DIFFUS ;

*** Rigidité ...

    ri1 = manu RIGI TYPE rigidite ma1 (mots UX UY)
          (prog 2 0 1 3)  QUEL ;

    cl1   = bloq mini uy p1 ;
    fodi1 = depi cl1 (valrefy * 0.99) ;
    fodi2 = depi cl1 (valrefy * 1.01) ;

*** Force ...

*    fo1 = forc FX 2 FY 7 p1 ;
    fo1 = ri1 * chpotst ;

*** Résolution ...

    list ri1 ;
    de0 = reso ri1 fo1 ;
    valux = extr de0 UX p1 ;
    valuy = extr de0 UY p1 ;
    diffx = valux - valrefx ;
    diffy = valuy - valrefy ;
list diffx ;
list diffy ;
    si((abs diffx > 1.e-10) ou (abs diffy > 1.e-10)) ;
       erre 5 ;
    finsi ;

    ritot = ri1 et cl1 ;
    fotot = fo1 et fodi1 ;
    de1 = reso ritot fotot ;
    valux = extr de1 UX p1 ;
    valuy = extr de1 UY p1 ;
    diffx = valux - valrefx ;
    diffy = valuy - valrefy ;
list diffx ;
list diffy ;
    si((abs diffx > 1.e-10) ou (abs diffy > 1.e-10)) ;
       erre 5 ;
    finsi ;

    ritot = ri1 et cl1 ;
    fotot = fo1 et fodi2 ;
    de2 = reso ritot fotot ;
    valux = extr de2 UX p1 ;
    valuy = extr de2 UY p1 ;
    diffx = valux - valrefx ;
    diffy = valuy - (valrefy * 1.01) ;
list diffx ;
list diffy ;
    si((abs diffx > 1.e-10) ou (abs diffy > 1.e-7)) ;
       erre 5 ;
    finsi ;


   fin ;