J:/Tecgraf/prj/femoop/trunk/path/path.h

// -------------------------------------------------------------------------
// Path.h - This header file contains the public data structure definitions
//          for the path control problem class within FEMOOP.
// -------------------------------------------------------------------------
//
// $Author: evandro $
// $Revision: 1.4 $
// $Date: 2000/04/05 18:46:38 $
// $State: Exp $
//
// -------------------------------------------------------------------------

#ifndef _PATH_H
#define _PATH_H

#include "ctrl.h"

// -------------------------------------------------------------------------
// Forward declarations:
//
class cVector;
class cSysMatrix;

// -------------------------------------------------------------------------
// Path Control Algorithm types:
//
typedef enum _pathtype
{
 NRAPHSON,     // Newton-Raphson (load control)
 DCM,          // Displacement control method
 GDCM,         // Generalized displacement control method
 INIORTAL,     // Initialy orthogonal arc-length
 UPDORTAL,     // Updated orthogonal arc-length
 CILAL,        // Cilindrical (quadratic) arc-length
 CONLINAL      // Cons. linearized arc-length

} ePathType;

// -------------------------------------------------------------------------
// Path class:
//
class cPath : public cCtrl
{
 protected:
  static int       _DetStabPnt;   // Flag for stability point evaluation
  static int       _BranchSwitch; // Flag for branch-switching
  static int       _TargetIte;    // Target number of iterations
  static int       _PrevIte;      // Number of iterations of previous step
  static double    _IncFactor;    // Incremental load factor of current step
  static double    _SwitchFactor; // Size of 1st step in secondary path
  static ePathType _PathType;     // Path algorithm
         double    _CSP0,_CSP1;   // Current stiffness parameter

 public:
  static  void   ReadTargetIte(void);
  static  void   ReadStabPnt  (void);
                 cPath        (void);
  virtual       ~cPath        (void) { }
  virtual void   Solver       (void);

 protected:
          int    StabSolver   (int, int, double, double, cSysMatrix *,
                               cVector &, cVector &, cVector &, cVector &, 
                               cVector &, cVector &, cVector &, cVector &,
                               cVector &, cVector &);
          void   CompStabPoint(int, int, double, double, cSysMatrix *,
                               cVector &, cVector &, cVector &, cVector &,
                               cVector &, cVector &, cVector &, cVector &,
                               cVector &, cVector &, int *);
  virtual int    Convergence  (cVector &, cVector &);
  virtual double IncLoadFac   (cVector &, cVector &, cVector &)
                              { return 0.0; }
};

#endif

---