Notices | Cast3M

$$$$ RESO     NOTICE  PV090527  24/11/12    21:15:09     12068          
                                             DATE     24/11/12
                                             
  Operateur RESO                           Voir aussi :
    --------------  
 CHPO1 (CHPO2 ..) = RESO RIG1 (PREC) CHPO3 (CHPO4..)
 ('NOID')('NOUNIL')('STAB')('ELIM' NBPASSE)('NOSTAB')('SOUC')
 (CHPOF) ('INIB' BLO1 LENTI1);

 (NB_MOD_RIG CHPO_CONTR) ... = ...  ('ENSE')          ;

    LICHP2 = RESO RIG1 LICHP1 ;

    (TAB2 = ) RESO RIG1 TAB1 ;

    Objet :
    _______

    L'operateur RESO construit une solution, si elle existe, du systeme
lineaire :   CHPO3 = RIG1 CHPO1 .
    L'operateur RESO construit les solutions, si elles existent, 
  de chacun des systemes lineaire :   CHPO3 = RIG1 CHPO1 ,
     CHPO1 pris dans LICHP1, ou dans TAB1, 
     CHPO3 range dans LICHP2, resp. TAB1 

    Commentaire :
    _____________

    RIG1  : objet de type RIGIDITE.

    CHPO3 : objet de type CHPOINT.

    CHPO1 : objet de type CHPOINT dont les composantes sont les duales
            de celles de CHPO3 par rapport a RIG1.

    LICHP1, LICHP2 : objet de type LISTCHPO

    TAB1 : objet de type TABLE, de sous-type 'LIAISONS_STATIQUES'
           Les indices sont de type ENTIER, pointent sur des objets
           TABLE, comportant les entrees
           - 'BLOCAGE', type RIGIDITE
           - 'FORCE', type CHPOINT, 2nd membre du systeme
           et completes par
           - 'DEFORMEE' type CHPOINT, solution du systeme
           - 'POINT_REPERE', type POINT, associe au deplacement calcule
     
    TAB2 : type TABLE, optionnel

    Remarque :
    __________

 1- En presence d'une famille de seconds membres CHPO3, CHPO4 ...,
l'operateur RESO construit la famille de solutions CHPO1, CHPO2 ..
respectivement associee.

 2- Si RIG1 contient des matrices issues de conditions unilaterales,
RESO appelle la procedure UNILATER pour fournir une solution du systeme.
Si il y a des matrices de frottement, il faut fournir le champ CHPOF
de forces limite de frottement.

 3- Les mots-cle 'NOID' et 'NOUNIL' sont utiles quand on emploie RESO a
    l'interieur d'une procedure :

      - 'NOID' desactive la vÃ©rificatin du residu ce qui autorise la
      resolution du systeme avec comme second membre la restriction de 
      CHPO3, (CHPO4 ..) a l'espace cible de RIG1.


      - 'NOUNIL' permet de resoudre le systeme en ignorant le caractere
        eventuellement unilateral de RIG1.

      - 'INIB' BLO1 LENTI1 permet dans le cas de contact d'indiquer un
        etat de contact initial. Voir la procedure UNILATER.

 4- Le mot cle 'STAB' fait utiliser pour la resolution un operateur 
    rendu positif par augmentation diagonale. CHPOx peut donc ne pas
    etre solution du probleme initial.
    Le mot-clÃ© 'NOSTAB' fait gconserver l'opÃ©rateur fourni.

 5- Le mot cle 'SOUC' provoque l'emission d'un soucis au lieu d'une
    erreur en cas d'impossibilite de resoudre le systeme.

 6- Le mot cle 'ELIM' permet de regler le nombre NBPASSE de passes
    d'elimination des inconnues soumises a des conditions ou relations
    imposees.

 7- Par "OPTION RESO DIRECTE" ou "OPTION RESO ITERATIVE"
    on peut choisir soit une resolution par methode de CROUT 
    ( methode par defaut) soit une resolution par methode de gradients
    conjugues avec preconditionnement ILU0 stabilise.

 8- Les champs par points CHPO1(2..) obtenus sont de nature diffuse.

 9- Le mot-cle 'ENSE' indique que au cas ou le systeme est singulier
    et la singularite est excite, RESO fournira un vecteur du noyau.
    Il y a alors deux resultats supplementaires : NB_MOD_RIG qui est le
    nombre de vecteurs du noyau retournes et CHPO_CONTR qui est le 
    champ des degres de liberte sur lesquels on ajoute une raideur
    pour pouvoir resoudre le systeme.

10- PREC est la precision de l'operation. Le defaut est 1D-18.

    Exemple :
    _________

    RIG1 etant la raideur d'une structure , FORC1 un champ de force
s'exerçant sur cette structure, on obtiendra le champ de deplacements
DEP1 en resultant, par l'instruction :

                  DEP1 = RESO  RIG1 FORC1 ;