Manual pages | Cast3M

$$$$ RESO     NOTICE  CHAT      11/09/12    21:17:58     7124           
                                             DATE     11/09/12
                                             
  Operateur RESO                           Voir aussi :
    --------------   
 CHPO1 (CHPO2 ..) = RESO RIG1 (PREC) CHPO3 (CHPO4..)('NOID')('NOUNIL') 
                         |      (CHPOF) ('INIB' BLO1 LENTI1)
 (NB_MOD_RIG MAIL_CONTR) |      ('ENSE')          ;

    LICHP2 = RESO RIG1 LICHP1 ;

    (TAB2 = ) RESO RIG1 TAB1 ;

    Description :
    ____________

    The RESO operator constructs a solution (if any) of the linear
 system :   RIG1 CHPO1 = CHPO3 .
    The RESO operator constructs the solutions (if any) of the linear
 systems :   RIG1 CHPO1 = CHPO3 as,
    CHPO1 is taken in LICHP1, or TAB1,
    CHPO3 is stored in LICHP2, resp. TAB1

    Contents :
    _________

    RIG1  : RIGIDITE type object

    CHPO3 : CHPOINT type object

    CHPO1 : CHPOINT type object whose components are the duals of
            those of CHPO3 with respect to RIG1.
  
    TAB1 : type TABLE object, subtype 'LIAISONS_STATIQUES'
           The indices belong to the ENTIER type, and point at 
           objects of TABLE type, which include 
           - 'BLOCAGE', RIGIDITE type
           - 'FORCE', CHPOINT type, 2nd member of the system
           and are completed with
           - 'DEFORMEE' CHPOINT type, solution of the system
           - 'POINT_REPERE', POINT type, associated with the
           computed displacement
     
    TAB2 : TABLE type, optional

    Note :
    ________

 1- With a second-member family CHPO3, CHPO4 ..., the RESO operator
 constructs the respectively associated family of solutions CHPO1,
 CHPO2 .. 
 
 2- If RIG1 contains matrices stemmineg from unilateral solutions,
 RESO calls on the UNILATER procedure to provide a solution of 
 the system.
 If there are friction stifness, the field CHPOF of limit friction
 strengths is required.
  
 3- The key words 'NOID' and 'NOUNIL' are useful when RESO is employed
 within the procedure :

      - 'NOID' enables the user to solve the system with the 
 restriction of CHPO3, (CHPO4 ..) from the subspace connected with
 RIG1 as second member 
 
      - 'NOUNIL' enables the user to solve the system by possibly
 ignoring the unilateral aspect of RIG1 

      - 'INIB' BLO1 LENTI1 enables, in case of contact, the
 specification of a state of initial contact. See UNILATER procedure

 4- Withh "OPTION RESO DIRECTE" or "OPTION RESO ITERATIVE" it is
    possible to choose between a CROUT method ( default method) and
    a method of conjugate gradients with preconditioning
    (incomplete Crout).

 5- The result CHPO1(2..) have a 'diffuse' nature.

 6- The keyword ('ENSE') tells RESO to return a vector of the kernel
    in the case where the system is singular and the singularity is
    excited. Then there are two new results: NB_MOD_RIG which is the 
    number of kernel-vectors and MAIL_CONTR which is the mesh of
    the constrained nodes needed to solve the system.

 7- PREC is the precision of the operation. By default it is 1D-15.    

   Example :
    _________

    RIG1 being the stiffness of a structure, and FORC1 a force
 field applying to this structure, the displacement field will
 be obtained by means of the following instruction :

                  DEP1 = RESO  RIG1 FORC1 ;