Télécharger reso.notice

Retour à la liste thématique

Afficher cette notice en

Numérotation des lignes :
   1 : $$$$ RESO     NOTICE  PV        16/11/03    21:15:15     9158           
   2 :                                              DATE     16/11/03
   3 :                                              
   4 :    Operateur RESO                           Voir aussi :
   5 :     --------------  
   6 :  CHPO1 (CHPO2 ..) = RESO RIG1 (PREC) CHPO3 (CHPO4..)
   7 :  ('NOID')('NOUNIL')('STAB')('NOEL')('NOSTAB')
   8 :  (CHPOF) ('INIB' BLO1 LENTI1);
   9 : 
  10 :  (NB_MOD_RIG MAIL_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 :   RIG1 CHPO1 = CHPO3 .
  21 :     L'operateur RESO construit les solutions, si elles existent, 
  22 :   de chacun des systemes lineaire :   RIG1 CHPO1 = CHPO3 ,
  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' permet de resoudre le systeme avec comme second membre la
  65 :         restriction de CHPO3, (CHPO4 ..) a l'espace cible de RIG1.
  66 : 
  67 :       - 'NOUNIL' permet de resoudre le systeme en ignorant le caractere
  68 :         eventuellement unilateral de RIG1.
  69 : 
  70 :       - 'INIB' BLO1 LENTI1 permet dans le cas de contact d'indiquer un
  71 :         etat de contact initial. Voir la procedure UNILATER.
  72 : 
  73 :  4- Le mot cle 'STAB' fait utiliser pour la resolution un operateur 
  74 :     rendu positif par augmentation diagonale. CHPOx peut donc ne pas
  75 :     etre solution du probleme initial.
  76 :     Le mot-clé 'NOSTAB' fait gconserver l'opérateur fourni.
  77 : 
  78 :  5- Le mot cle 'NOEL' desactive l'ielimination des inconnues soumises
  79 :     a des conditions ou relations imposees.
  80 : 
  81 :  6- Par "OPTION RESO DIRECTE" ou "OPTION RESO ITERATIVE"
  82 :     on peut choisir soit une resolution par methode de CROUT 
  83 :     ( methode par defaut) soit une resolution par methode de gradients
  84 :     conjugues avec preconditionnement ILU0 stabilise.
  85 : 
  86 :  7- Les champs par points CHPO1(2..) obtenus sont de nature diffuse.
  87 : 
  88 :  8- Le mot-cle 'ENSE' indique que au cas ou le systeme est singulier
  89 :     et la singularite est excite, RESO fournira un vecteur du noyau.
  90 :     Il y a alors deux resultats supplementaires : NB_MOD_RIG qui est le
  91 :     nombre de vecteurs du noyau retournes et MAIL_CONTR qui est le 
  92 :     maillage des noeuds sur lesquels on a applique une contrainte
  93 :     pour pouvoir resoudre le systeme.
  94 : 
  95 :  9- PREC est la prevision de l'operation. Le defaut est 1D-15.
  96 : 
  97 :     Exemple :
  98 :     _________
  99 : 
 100 :     RIG1 etant la raideur d'une structure , FORC1 un champ de force
 101 : s'exerçant sur cette structure, on obtiendra le champ de deplacements
 102 : DEP1 en resultant, par l'instruction :
 103 : 
 104 :                   DEP1 = RESO  RIG1 FORC1 ;
 105 : 

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