GursonModel.hpp

Go to the documentation of this file.

//*****************************************************************//
//    Albany 2.0:  Copyright 2012 Sandia Corporation               //
//    This Software is released under the BSD license detailed     //
//    in the file "license.txt" in the top-level Albany directory  //
//*****************************************************************//

#if !defined(GursonModel_hpp)
#define GursonModel_hpp

#include <Intrepid_MiniTensor.h>
#include "Phalanx_ConfigDefs.hpp"
#include "Phalanx_Evaluator_WithBaseImpl.hpp"
#include "Phalanx_Evaluator_Derived.hpp"
#include "Phalanx_MDField.hpp"
#include "Albany_Layouts.hpp"
#include "LCM/models/ConstitutiveModel.hpp"

namespace LCM
{

template<typename EvalT, typename Traits>
class GursonModel: public LCM::ConstitutiveModel<EvalT, Traits>
{
public:

  typedef typename EvalT::ScalarT ScalarT;
  typedef typename EvalT::MeshScalarT MeshScalarT;
  typedef typename Sacado::mpl::apply<FadType, ScalarT>::type DFadType;

  using ConstitutiveModel<EvalT, Traits>::num_dims_;
  using ConstitutiveModel<EvalT, Traits>::num_pts_;
  using ConstitutiveModel<EvalT, Traits>::field_name_map_;

  GursonModel(Teuchos::ParameterList* p,
      const Teuchos::RCP<Albany::Layouts>& dl);

  virtual
  ~GursonModel()
  {};

  virtual
  void
  computeState(typename Traits::EvalData workset,
      std::map<std::string, Teuchos::RCP<PHX::MDField<ScalarT> > > dep_fields,
      std::map<std::string, Teuchos::RCP<PHX::MDField<ScalarT> > > eval_fields);

private:

  GursonModel(const GursonModel&);

  GursonModel& operator=(const GursonModel&);

  RealType sat_mod_, sat_exp_;

  RealType f0_;

  RealType kw_;

  RealType eN_, sN_, fN_;

  RealType fc_, ff_;

  RealType q1_, q2_, q3_;

  ScalarT
  YieldFunction(Intrepid::Tensor<ScalarT> const & s, ScalarT const & p,
      ScalarT const & fvoid, ScalarT const & eq, ScalarT const & K,
      ScalarT const & Y, ScalarT const & jacobian, ScalarT const & E);

  void
  ResidualJacobian(std::vector<ScalarT> & X,
      std::vector<ScalarT> & R, std::vector<ScalarT> & dRdX, const ScalarT & p,
      const ScalarT & fvoid, const ScalarT & eq, Intrepid::Tensor<ScalarT> & s,
      const ScalarT & mu, const ScalarT & kappa, const ScalarT & K,
      const ScalarT & Y, const ScalarT & jacobian);

};
}

#endif