phaseFieldLinearizedElasticity.hpp

Go to the documentation of this file.

/*
Brian Staber (brian.staber@gmail.com)
*/

#ifndef PHASEFIELDLINEARIZEDELASTICITY_HPP
#define PHASEFIELDLINEARIZEDELASTICITY_HPP

#include "linearizedElasticity.hpp"
#include "damageField.hpp"

class phaseFieldLinearizedElasticity : public linearizedElasticity{

public:
  Epetra_SerialDenseMatrix elasticity;
  Teuchos::RCP<damageField> damageInterface;

  double gc, lc;
  double E, nu, lambda, mu;

  Epetra_Vector * damageSolutionOverlaped;

  phaseFieldLinearizedElasticity();
  ~phaseFieldLinearizedElasticity();

  Epetra_Map constructGaussMap();

  void initialize(Epetra_Comm & comm, Teuchos::ParameterList & Parameters);
  void computeDisplacement(Teuchos::ParameterList & ParameterList,
                           Epetra_FECrsMatrix & matrix, Epetra_Vector & lhs, Epetra_FEVector & rhs,
                           double & bc_disp);
  void updateDamageHistory(Epetra_Vector & damageHistory,
                           Epetra_Vector & displacement,
                           Epetra_Map & GaussMap);
  void staggeredAlgorithmDirichletBC(Teuchos::ParameterList & ParametersList,
                                     bool print);

  void get_elasticity_tensor(unsigned int & e_lid, unsigned int & gp, Epetra_SerialDenseMatrix & tangent_matrix);
  void get_elasticity_tensor_for_recovery(unsigned int & e_lid, Epetra_SerialDenseMatrix & tangent_matrix);

};
#endif