gmpFemNonLinearSolverOld.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_NONLINEARSOLVEROLD_H_
#define _GEMA_PLUGIN_FEM_NONLINEARSOLVEROLD_H_
 
#include "gmpFemSolver.h"
#include "gmpFemSolverOptions.h"
 
class GMP_FEM_PROCESS_API_EXPORT GmpFemNonLinearSolverOld : public GmpFemSolver
{
public:
  GmpFemNonLinearSolverOld(GmElementMesh*, GmSimulationData*, const QVector<GmpFemPhysics*>&, GmNumSolver*,
                        const GmpFemSolverOptions&, const GmLogCategory&);
 
  virtual ~GmpFemNonLinearSolverOld();
 
  void enableFastUpdate(bool mode) { _fastUpdate = mode; }
 
  virtual bool init();
 
  virtual bool initResultAttributes(int fatherItemId);
 
  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();
  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 GmpFemPhysics::FemResultType updateResidualForceAndStiffnessMatrix(double l);
 
  virtual GmpFemPhysics::FemResultType fillStiffnessMatrixBroyden();
  
  virtual GmpFemPhysics::FemResultType fillStiffnessMatrixBFGS();
 
  bool solveLinearSystem(GmSolverMatrix* Keq, GmVector& Feq, GmVector& x, int iter=0);
 
  bool evaluateSolverTolerance(int iter);
 
  virtual bool staticSolver(double, double* newt);
 
  virtual bool transientSolver(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;
  int        _It;
 
  GmVector   _dx;
  GmVector   _xold;
  GmVector   _dxt;
  GmVector   _ddx;
  GmVector   _ddxt;
  GmVector   _ddxr;
  GmVector   _FeGs;
  GmVector   _r0;      
  GmVector   _normR;   
 
  QVector<int> _fixedDof;
 
  bool _fastUpdate; 
private:
  bool     _nonlinear; 
 
  int _femBroydenItem = -1;
 
};
 
#endif