gmpFemNonLinearSolver.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2014 by Carlos Augusto Teixera Mendes
** All rights reserved.
**
** This file is part of the "GeMA" software. It's use should respect
** the terms in the license agreement that can be found together
** with this source code.
** It is provided AS IS, with NO WARRANTY OF ANY KIND,
** INCLUDING THE WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR
** A PARTICULAR PURPOSE.
**
************************************************************************/

#ifndef _GEMA_PLUGIN_FEM_NONLINEARSOLVER_H_
#define _GEMA_PLUGIN_FEM_NONLINEARSOLVER_H_

#include "gmpFemSolver.h"
#include "gmpFemSolverOptions.h"

class GMP_FEM_PROCESS_API_EXPORT GmpFemNonLinearSolver : public GmpFemSolver
{
public:
  GmpFemNonLinearSolver(GmElementMesh*, GmSimulationData*, const QList<GmpFemPhysics*>&, GmNumSolver*,
                        const GmpFemSolverOptions&, const GmLogCategory&);

  virtual ~GmpFemNonLinearSolver();

  void enableFastUpdate(bool mode) { _fastUpdate = mode; }

  virtual bool init();
 
  virtual bool step(double dt, double* newt, bool* conv, double* err, int* niter, bool eval);

  virtual bool geostatic();

  virtual bool pos();

  virtual bool crackDetection();

  virtual bool crackPropagation();

  virtual bool updateBC(double dt, int iter, bool updateBc);

protected:
  virtual bool initElementSets(GmNumSolver* solver);
  
  virtual void setElementSets(bool, bool, bool, bool, bool);
  
  virtual bool isNonLinear() const;
  
  virtual bool isTransient() const
  {
    return _solverOptions._type == GmpFemSolverOptions::transientSolver;
  }

  virtual bool convergence() const;
  virtual bool convergenceNew() const;

  virtual void saveInternalForceAttributes();

  virtual void update();

  virtual void restore();

  virtual void setLoad();

  virtual void setTime();
  
  virtual void setState();

  virtual double loadPredictor();

  virtual double loadCorrector();

  virtual double calcLocalError(GmVector x, GmVector dx);

  virtual bool stepStaticNonlinear(double, double* newt);

  virtual GmpFemPhysics::FemResultType updateResidualForceAndStiffnessMatrix(double l);

  virtual GmpFemPhysics::FemResultType fillStiffnessMatrixBroyden();
  
  virtual GmpFemPhysics::FemResultType fillStiffnessMatrixBFGS();

  bool solveLinearSystem(GmSolverMatrix* Keq, GmVector& Feq, GmVector& x);

  virtual bool stepTransientNonLinear(double, double* newt);

  virtual bool stepTransientNonlinearAut(double, double* newt);

  //bool prepareTransientSystem(double dt);
  GmpFemPhysics::FemResultType prepareTransientSystem(double dt);
  GmpFemPhysics::FemResultType prepareTransientSystemPredict(double dt);

  void vectorTreatment(GmVector& fint);
  bool timeIncrementAdjustment(double dt, int iter, GmVector& X1, GmVector& X0);
  // Dump methods
  virtual bool addStateItemsToGroup(GmStateDump* state, int groupId);
  virtual bool fillStateControlMapData(QVariantMap* map);
  virtual bool stateControlMapDataLoaded(QVariantMap* map);
  virtual bool stateLoadedCalcDerivedResultsNonLinear();

  virtual bool fillStateControlMapDataStaticSolver(QVariantMap* map);

  virtual bool fillStateControlMapDataTransientSolver(QVariantMap* map);

  //atributtes
  int        _dofIndex;
  double     _dl;
  double     _dl0;
  double     _ddl;
  double     _arc;
  double     _Lold;
  double     _tin;  
  double     _dt;
  double     _oldIter;
  double     _hlast;
  int        _It;

  GmVector   _v;
  GmVector   _a;
  GmVector   _dx;
  GmVector   _dv;
  GmVector   _xold;
  GmVector   _vold;
  GmVector   _dxt;
  GmVector   _ddx;
  GmVector   _ddv;
  GmVector   _ddxt;
  GmVector   _ddxr;
  GmVector   _FeGs;
  GmVector   _feLcs;
  GmVector   _fiLcs;
  GmVector   _r0;

  QList<int> _fixedDof;

  bool _fastUpdate; 
private:
  bool     _nonlinear; 

};

#endif