Accueil | Cast3M

   1 : $$$$ RESO     NOTICE  PV090527  24/11/12    21:15:09     12068          
   2 :                                              DATE     24/11/12
   3 :                                              
   4 :    Operateur RESO                           Voir aussi :
   5 :     --------------  
   6 :  CHPO1 (CHPO2 ..) = RESO RIG1 (PREC) CHPO3 (CHPO4..)
   7 :  ('NOID')('NOUNIL')('STAB')('ELIM' NBPASSE)('NOSTAB')('SOUC')
   8 :  (CHPOF) ('INIB' BLO1 LENTI1);
   9 : 
  10 :  (NB_MOD_RIG CHPO_CONTR) ... = ...  ('ENSE')          ;
  11 : 
  12 :     LICHP2 = RESO RIG1 LICHP1 ;
  13 : 
  14 :     (TAB2 = ) RESO RIG1 TAB1 ;
  15 : 
  16 :     Objet :
  17 :     _______
  18 : 
  19 :     L'operateur RESO construit une solution, si elle existe, du systeme
  20 : lineaire :   CHPO3 = RIG1 CHPO1 .
  21 :     L'operateur RESO construit les solutions, si elles existent, 
  22 :   de chacun des systemes lineaire :   CHPO3 = RIG1 CHPO1 ,
  23 :      CHPO1 pris dans LICHP1, ou dans TAB1, 
  24 :      CHPO3 range dans LICHP2, resp. TAB1 
  25 : 
  26 :     Commentaire :
  27 :     _____________
  28 : 
  29 :     RIG1  : objet de type RIGIDITE.
  30 : 
  31 :     CHPO3 : objet de type CHPOINT.
  32 : 
  33 :     CHPO1 : objet de type CHPOINT dont les composantes sont les duales
  34 :             de celles de CHPO3 par rapport a RIG1.
  35 : 
  36 :     LICHP1, LICHP2 : objet de type LISTCHPO
  37 : 
  38 :     TAB1 : objet de type TABLE, de sous-type 'LIAISONS_STATIQUES'
  39 :            Les indices sont de type ENTIER, pointent sur des objets
  40 :            TABLE, comportant les entrees
  41 :            - 'BLOCAGE', type RIGIDITE
  42 :            - 'FORCE', type CHPOINT, 2nd membre du systeme
  43 :            et completes par
  44 :            - 'DEFORMEE' type CHPOINT, solution du systeme
  45 :            - 'POINT_REPERE', type POINT, associe au deplacement calcule
  46 :      
  47 :     TAB2 : type TABLE, optionnel
  48 : 
  49 :     Remarque :
  50 :     __________
  51 : 
  52 :  1- En presence d'une famille de seconds membres CHPO3, CHPO4 ...,
  53 : l'operateur RESO construit la famille de solutions CHPO1, CHPO2 ..
  54 : respectivement associee.
  55 : 
  56 :  2- Si RIG1 contient des matrices issues de conditions unilaterales,
  57 : RESO appelle la procedure UNILATER pour fournir une solution du systeme.
  58 : Si il y a des matrices de frottement, il faut fournir le champ CHPOF
  59 : de forces limite de frottement.
  60 : 
  61 :  3- Les mots-cle 'NOID' et 'NOUNIL' sont utiles quand on emploie RESO a
  62 :     l'interieur d'une procedure :
  63 : 
  64 :       - 'NOID' desactive la vÃ©rificatin du residu ce qui autorise la
  65 :       resolution du systeme avec comme second membre la restriction de 
  66 :       CHPO3, (CHPO4 ..) a l'espace cible de RIG1.
  67 : 
  68 : 
  69 :       - 'NOUNIL' permet de resoudre le systeme en ignorant le caractere
  70 :         eventuellement unilateral de RIG1.
  71 : 
  72 :       - 'INIB' BLO1 LENTI1 permet dans le cas de contact d'indiquer un
  73 :         etat de contact initial. Voir la procedure UNILATER.
  74 : 
  75 :  4- Le mot cle 'STAB' fait utiliser pour la resolution un operateur 
  76 :     rendu positif par augmentation diagonale. CHPOx peut donc ne pas
  77 :     etre solution du probleme initial.
  78 :     Le mot-clÃ© 'NOSTAB' fait gconserver l'opÃ©rateur fourni.
  79 : 
  80 :  5- Le mot cle 'SOUC' provoque l'emission d'un soucis au lieu d'une
  81 :     erreur en cas d'impossibilite de resoudre le systeme.
  82 : 
  83 :  6- Le mot cle 'ELIM' permet de regler le nombre NBPASSE de passes
  84 :     d'elimination des inconnues soumises a des conditions ou relations
  85 :     imposees.
  86 : 
  87 :  7- Par "OPTION RESO DIRECTE" ou "OPTION RESO ITERATIVE"
  88 :     on peut choisir soit une resolution par methode de CROUT 
  89 :     ( methode par defaut) soit une resolution par methode de gradients
  90 :     conjugues avec preconditionnement ILU0 stabilise.
  91 : 
  92 :  8- Les champs par points CHPO1(2..) obtenus sont de nature diffuse.
  93 : 
  94 :  9- Le mot-cle 'ENSE' indique que au cas ou le systeme est singulier
  95 :     et la singularite est excite, RESO fournira un vecteur du noyau.
  96 :     Il y a alors deux resultats supplementaires : NB_MOD_RIG qui est le
  97 :     nombre de vecteurs du noyau retournes et CHPO_CONTR qui est le 
  98 :     champ des degres de liberte sur lesquels on ajoute une raideur
  99 :     pour pouvoir resoudre le systeme.
 100 : 
 101 : 10- PREC est la precision de l'operation. Le defaut est 1D-18.
 102 : 
 103 :     Exemple :
 104 :     _________
 105 : 
 106 :     RIG1 etant la raideur d'une structure , FORC1 un champ de force
 107 : s'exerçant sur cette structure, on obtiendra le champ de deplacements
 108 : DEP1 en resultant, par l'instruction :
 109 : 
 110 :                   DEP1 = RESO  RIG1 FORC1 ;
 111 :