gmpFemNLSolver.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_NLSOLVER_H_
#define _GEMA_PLUGIN_FEM_NLSOLVER_H_

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

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

  //destructor
  virtual ~GmpFemNLSolver();
  
  void enableFastUpdate(bool mode) { _fastUpdate = mode; }

  // initializing Fem NL Solver
  virtual bool init();
 

  //virtual bool step(double dt, double* newt);
  virtual bool runStep(double dt, double* newt, bool* conv, double* err, int* niter, bool eval=false);
  // Geostatic procedure
  virtual bool geostatic(QString svId = "u", QString gaId = "Em");

protected:
  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 stateControlMapDataLoadedStaticSolver(QVariantMap* map);
  virtual bool fillStateControlMapDataTransientSolver(QVariantMap* map);
  virtual bool stateControlMapDataLoadedTransientSolver(QVariantMap* map);

//private:
  // general
  virtual bool initElementSets(GmNumSolver*);
  virtual void setElementSets(bool, bool, bool, bool, bool);
  virtual bool applyFixedBoundaryConditions(int*);
  virtual bool IntialFixedBoundaryConditions(int*);
  virtual bool applyNonzeroFixedBCsbyModification(int*);
  virtual bool isDynamic() const;
  virtual bool isNonlinear() const;

  virtual bool convergence() const;
  virtual bool convergenceN(GmVector& _erL, GmVector& _erD) const;
  virtual bool convergenceDC(GmVector& _er) const;
  virtual void update();
  virtual void restore();

  //Corrector 
  //virtual bool corrector(GmVector& _er, GmVector) const;

  //Iterative Error
  virtual bool LoadError(GmVector& _er)const;
  virtual bool DispError(GmVector& _er) const;
  virtual bool VelocityError(GmVector& _er)const;
  virtual bool FluxError(GmVector& _er) const;
  virtual bool FluxMagnitude(GmVector& _qa) const;

  virtual bool MaxErr(GmVector, double& ermax, double &tolmax);
  virtual bool MaxVal(GmVector& array,  GmVector) const;
  virtual bool PhygroupName(int & kgroup,  QString var) const;
  virtual bool divergingIncrement(double ermax, double& ermax_1, double& ermax_2);
  virtual bool adaptativeTimeAbqs(bool conv, int iter, double dt, double& newdt, double errmax, int inc, bool iter4);
  virtual bool adaptativeTime(bool conv, int iter, double dt, double& newdt, double err, double tol);
  virtual void saveInternalForceAttributes(GmVector& dv);

  virtual void setLoad();
  virtual void setTime();
  virtual void setState();
  virtual void setTransientSystem(double);
  virtual void setTransientSystem2(double);

  virtual GmpFemPhysics::FemResultType updateResidualForceAndStiffnessMatrix(double l);
  virtual GmpFemPhysics::FemResultType fillStiffnessMatrixBroyden();
  virtual GmpFemPhysics::FemResultType fillStiffnessMatrixBFGS();

  virtual double evalLocalError(GmVector, GmVector);

  virtual double evalTimeIncrementSize(int, double, double, GmVector&);

  // Strategies of continuation methods adopted by static solver 
  virtual double loadPredictor();
  virtual double loadCorrector();

  // Solvers
  virtual bool staticLinearSolver(double step, double* newt, bool* conv, double* err, int* niter, bool eval);
  virtual bool staticSolver(double step, double* newt, bool* conv, double* err, int* niter, bool eval);
  virtual bool transientSolver(double dt, double* newt, bool* conv, double* err, int* niter, bool eval);
  virtual bool transientLinearSolver(double);
  virtual bool iterativeSolver(double dt, double* newt, bool* conv, double* err, int* niter, bool eval);
  
  // To check
  virtual bool stepStandardDc(double dt, double* newt, bool* conv, double* err, int* niter, bool eval);
  virtual bool stepTransientNonlinearAut(double dt, double* newt);
  virtual bool stepNonlinearAut(double dt, double* newt);

  // atributtes
  double     _dl;
  double     _dl0;
  double     _ddl;
  double     _Lold;
  double     _tin;
  double     _attempt = 0.0;
  // double     _ts;                
  double     _dt;
  double     _oldIter;
  double     _hlast;

  // controls time incrementation
  int        _ninc = 0;          
  int        _incT = 0;          
  bool       _iter4old = false;  



  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   _erL;
  GmVector   _erD;
  GmVector   _r0;
  // controls fluxe
  GmVector     _qta;      
  GmVector     _qtamax;   


  QList<int> _fixedDof;

  bool _fastUpdate; 

  private:
    //GmVector _FeGs;

};

#endif