pzelastoplasticanalysis.h

Go to the documentation of this file.

#ifndef ELASTOPLASTICANALYSIS_H
#define ELASTOPLASTICANALYSIS_H

#include "pznonlinanalysis.h"
#include "pzcompel.h"
#include "TPZGeoElement.h"
#include "pzfmatrix.h"
#include "pzvec.h"
#include "pzpostprocanalysis.h"
#include <iostream>


class TPZElastoPlasticAnalysis : public TPZNonLinearAnalysis {

public:
                TPZElastoPlasticAnalysis(TPZCompMesh *mesh,std::ostream &out);
                TPZElastoPlasticAnalysis();
                virtual ~TPZElastoPlasticAnalysis();
                
    virtual void IterativeProcessPrecomputedMatrix(std::ostream &out, REAL tol, int numiter, bool linesearch);
    
    virtual void IterativeProcess(std::ostream &out, REAL tol, int numiter, int niter_update_jac, bool linesearch);
    
                virtual void IterativeProcess(std::ostream &out,REAL tol,int numiter, bool linesearch, bool checkconv,bool &ConvOrDiverg);

                virtual void IterativeProcess(std::ostream &out,REAL tol,int numiter, bool linesearch, bool checkconv);
    
                //virtual REAL LineSearch(const TPZFMatrix<REAL> &Wn, const TPZFMatrix<REAL> &DeltaW, TPZFMatrix<REAL> &NextW, REAL RhsNormPrev, REAL &RhsNormResult, int niter);
    
    //Improved:A verification is made in order to check convergence.
    virtual REAL LineSearch(const TPZFMatrix<REAL> &Wn, const TPZFMatrix<REAL> &DeltaW, TPZFMatrix<REAL> &NextW, REAL RhsNormPrev, REAL &RhsNormResult, int niter, bool &converging );
                
                virtual void ManageIterativeProcess(std::ostream &out,REAL tol,int numiter,
                                                                                                                                                                int BCId, int nsteps, REAL PGRatio,
                                                                                                                                                                TPZFMatrix<REAL> & val1Begin, TPZFMatrix<REAL> & val1End,
                                                                                                                                                                TPZFMatrix<REAL> & val2Begin, TPZFMatrix<REAL> & val2End,
                                                                                                                                                                TPZPostProcAnalysis * ppAnalysis = NULL, int res = 0);
                
                virtual REAL AcceptSolution(const int ResetOutputDisplacements = 0);
    
    virtual void LoadSolution();

    
    TPZFMatrix<REAL> &CumulativeSolution()
    {
        return fCumSol;
    }
                
                void SetPrecond(TPZMatrixSolver<REAL> &precond);
                
                void SetBiCGStab(int numiter, REAL tol);
                
                void SetBiCGStab_Jacobi(int numiter, REAL tol);
                
                void SetLU();
                
                void TransferSolution(TPZPostProcAnalysis & ppanalysis);
                
    void AddNoPenetration(int matid, int direction)
    {
        fMaterialIds.insert(std::pair<int,int>(matid, direction));
    }
    
    void IdentifyEquationsToZero();
    
    void GetActiveEquations(TPZVec<int64_t> &activeEquations);
    
protected:      
                
                void SetUpdateMem(int update);
                
                void UpdatePrecond();
                
                
public:
                
                void CheckConv(std::ostream &out, REAL range);
                
                virtual void ComputeTangent(TPZFMatrix<REAL> &tangent, TPZVec<REAL> &coefs, int icase);
                
                virtual int NumCases();
                
                virtual void Residual(TPZFMatrix<REAL> &residual, int icase);
                
                static void SetAllCreateFunctionsWithMem(TPZCompMesh *cmesh);
                
    bool IsMultiPhysicsConfiguration()
    {
        return fMultiPhysics != NULL;
    }
    
    void SetMultiPhysics(TPZCompMesh *mphysics, TPZVec<TPZCompMesh *> &meshvec)
    {
        fMultiPhysics = mphysics;
        fMeshVec = meshvec;
    }
    
    void ResetMultiPhysics()
    {
        fMultiPhysics = 0;
        fMeshVec.Resize(0);
    }
    
                
protected:
                
                /* @brief Cumulative solution vector*/
                TPZFMatrix<REAL> fCumSol;
                
                TPZMatrixSolver<REAL> * fPrecond;
    
    std::set<int64_t> fEquationstoZero;
    
    // the second value of the map indicates x (0) or y (1) restraint
    std::multimap<int,int> fMaterialIds;
    
    TPZCompMesh *fMultiPhysics;
    
    TPZManVector<TPZCompMesh *,2> fMeshVec;
                
                static TPZCompEl * CreateCubeElWithMem(  TPZGeoEl *gel, TPZCompMesh &mesh, int64_t &index);
                static TPZCompEl * CreateLinearElWithMem(TPZGeoEl *gel, TPZCompMesh &mesh, int64_t &index);
                static TPZCompEl * CreatePointElWithMem( TPZGeoEl *gel, TPZCompMesh &mesh, int64_t &index);
                static TPZCompEl * CreatePrismElWithMem( TPZGeoEl *gel, TPZCompMesh &mesh, int64_t &index);
                static TPZCompEl * CreatePyramElWithMem( TPZGeoEl *gel, TPZCompMesh &mesh, int64_t &index);
                static TPZCompEl * CreateQuadElWithMem(  TPZGeoEl *gel, TPZCompMesh &mesh, int64_t &index);
                static TPZCompEl * CreateTetraElWithMem( TPZGeoEl *gel, TPZCompMesh &mesh, int64_t &index);
                static TPZCompEl * CreateTriangElWithMem(TPZGeoEl *gel, TPZCompMesh &mesh, int64_t &index);
                
};

#endif