NewmarkMats.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#ifndef _NEWMARK_MATS_H
#define _NEWMARK_MATS_H
 
#include "ElmMats.h"
 
 
class NewmarkMats : public ElmMats
{
public:
  NewmarkMats(double a1, double a2, double b = 0.0, double c = 0.0,
              bool generalizedAlpha = false);
  virtual ~NewmarkMats() {}
 
  virtual void setStepSize(double dt, int it) { h = dt; isPredictor = it == 0; }
 
  virtual const Matrix& getNewtonMatrix() const;
  virtual const Vector& getRHSVector() const;
 
  const Vector& dis(bool prev = false) const { return vec[prev ? 3 : 0]; }
  const Vector& vel(bool prev = false) const { return vec[prev ? 4 : 1]; }
  const Vector& acc(bool prev = false) const { return vec[prev ? 5 : 2]; }
 
protected:
  bool  isPredictor; 
  double h;          
 
  double alpha1; 
  double alpha2; 
  double beta;   
  double gamma;  
 
private:
  double alpha_m; 
  double alpha_f; 
  bool   slvDisp; 
};
 
#endif