26namespace HydroMechanics
30template <
typename ShapeFunctionDisplacement,
typename ShapeFunctionPressure,
43 std::size_t
const n_variables,
44 std::size_t
const local_matrix_size,
45 std::vector<unsigned>
const& dofIndex_to_localIndex,
47 bool const is_axially_symmetric,
52 double const )
override
62 std::vector<GlobalVector*>
const& x,
63 std::vector<NumLib::LocalToGlobalIndexMap const*>
const& dof_table,
64 std::vector<double>& cache)
const override;
68 std::vector<GlobalVector*>
const& x,
69 std::vector<NumLib::LocalToGlobalIndexMap const*>
const& dof_table,
70 std::vector<double>& cache)
const override;
74 std::vector<GlobalVector*>
const& x,
75 std::vector<NumLib::LocalToGlobalIndexMap const*>
const& dof_table,
76 std::vector<double>& cache)
const override;
80 std::vector<GlobalVector*>
const& ,
81 std::vector<NumLib::LocalToGlobalIndexMap const*>
const& ,
82 std::vector<double>& cache)
const override
90 std::vector<GlobalVector*>
const& ,
91 std::vector<NumLib::LocalToGlobalIndexMap const*>
const& ,
92 std::vector<double>& cache)
const override
100 std::vector<GlobalVector*>
const& ,
101 std::vector<NumLib::LocalToGlobalIndexMap const*>
const& ,
102 std::vector<double>& cache)
const override
110 std::vector<GlobalVector*>
const& ,
111 std::vector<NumLib::LocalToGlobalIndexMap const*>
const& ,
112 std::vector<double>& cache)
const override
119 const unsigned integration_point)
const override
124 return Eigen::Map<const Eigen::RowVectorXd>(N.data(), N.size());
129 Eigen::VectorXd
const& local_x,
130 Eigen::VectorXd
const& local_x_prev,
131 Eigen::VectorXd& local_rhs,
132 Eigen::MatrixXd& local_Jac)
override;
135 double const t,
double const dt,
136 Eigen::Ref<const Eigen::VectorXd>
const& p,
137 Eigen::Ref<const Eigen::VectorXd>
const& p_prev,
138 Eigen::Ref<const Eigen::VectorXd>
const& u,
139 Eigen::Ref<const Eigen::VectorXd>
const& u_prev,
140 Eigen::Ref<Eigen::VectorXd> rhs_p, Eigen::Ref<Eigen::VectorXd> rhs_u,
141 Eigen::Ref<Eigen::MatrixXd> J_pp, Eigen::Ref<Eigen::MatrixXd> J_pu,
142 Eigen::Ref<Eigen::MatrixXd> J_uu, Eigen::Ref<Eigen::MatrixXd> J_up);
145 double const t,
double const dt,
146 Eigen::VectorXd
const& local_x)
override;
149 double const t,
double const dt,
150 Eigen::Ref<const Eigen::VectorXd>
const& p,
151 Eigen::Ref<const Eigen::VectorXd>
const& u);
154 double const t, Eigen::Ref<Eigen::VectorXd> p);
169 ShapeFunctionDisplacement::NPOINTS>;
176 Eigen::aligned_allocator<IntegrationPointDataType>>
183 ShapeFunctionDisplacement::NPOINTS * GlobalDim;
188 typename ShapeMatricesTypeDisplacement::ShapeMatrices::ShapeType>
ShapeMatrixPolicyType< ShapeFunctionDisplacement, GlobalDim > ShapeMatricesTypeDisplacement
BMatrixPolicyType< ShapeFunctionDisplacement, GlobalDim > BMatricesType
std::vector< double > const & getIntPtFractureVelocity(const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
std::vector< double > const & getIntPtSigma(const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &cache) const override
std::vector< double > const & getIntPtFracturePermeability(const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
IntegrationPointDataMatrix< BMatricesType, ShapeMatricesTypeDisplacement, ShapeMatricesTypePressure, GlobalDim, ShapeFunctionDisplacement::NPOINTS > IntegrationPointDataType
std::vector< double > const & getIntPtFractureStress(const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
void postTimestepConcreteWithBlockVectors(double const t, double const dt, Eigen::Ref< const Eigen::VectorXd > const &p, Eigen::Ref< const Eigen::VectorXd > const &u)
static const int pressure_size
SecondaryData< typename ShapeMatricesTypeDisplacement::ShapeMatrices::ShapeType > _secondary_data
std::vector< IntegrationPointDataType, Eigen::aligned_allocator< IntegrationPointDataType > > _ip_data
void assembleWithJacobianConcrete(double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_x_prev, Eigen::VectorXd &local_rhs, Eigen::MatrixXd &local_Jac) override
std::vector< double > const & getIntPtEpsilon(const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &cache) const override
ShapeMatrixPolicyType< ShapeFunctionPressure, GlobalDim > ShapeMatricesTypePressure
void assembleBlockMatricesWithJacobian(double const t, double const dt, Eigen::Ref< const Eigen::VectorXd > const &p, Eigen::Ref< const Eigen::VectorXd > const &p_prev, Eigen::Ref< const Eigen::VectorXd > const &u, Eigen::Ref< const Eigen::VectorXd > const &u_prev, Eigen::Ref< Eigen::VectorXd > rhs_p, Eigen::Ref< Eigen::VectorXd > rhs_u, Eigen::Ref< Eigen::MatrixXd > J_pp, Eigen::Ref< Eigen::MatrixXd > J_pu, Eigen::Ref< Eigen::MatrixXd > J_uu, Eigen::Ref< Eigen::MatrixXd > J_up)
static const int displacement_size
static const int pressure_index
static const int displacement_index
std::vector< double > const & getIntPtDarcyVelocity(const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &cache) const override
HydroMechanicsLocalAssemblerMatrix(HydroMechanicsLocalAssemblerMatrix const &)=delete
void postTimestepConcreteWithVector(double const t, double const dt, Eigen::VectorXd const &local_x) override
void setPressureOfInactiveNodes(double const t, Eigen::Ref< Eigen::VectorXd > p)
HydroMechanicsLocalAssemblerMatrix(HydroMechanicsLocalAssemblerMatrix &&)=delete
static const int kelvin_vector_size
std::vector< double > const & getIntPtFractureAperture(const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
void preTimestepConcrete(std::vector< double > const &, double const, double const) override
HydroMechanicsProcessData< GlobalDim > & _process_data
Eigen::Map< const Eigen::RowVectorXd > getShapeMatrix(const unsigned integration_point) const override
Provides the shape matrix at the given integration point.
Eigen::Matrix< double, GlobalDim, 1 > GlobalDimVector
constexpr int kelvin_vector_dimensions(int const displacement_dim)
Kelvin vector dimensions for given displacement dimension.
Used for the extrapolation of the integration point values.