HydroMechanics/LocalAssembler/IntegrationPointDataFracture.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include <Eigen/Core>
 
#include "MaterialLib/FractureModels/FractureModelBase.h"
 
namespace ProcessLib
{
namespace LIE
{
namespace HydroMechanics
{
template <typename HMatricesType,
          typename ShapeMatrixTypeDisplacement,
          typename ShapeMatrixTypePressure,
          int DisplacementDim>
struct IntegrationPointDataFracture final
{
    explicit IntegrationPointDataFracture(
        MaterialLib::Fracture::FractureModelBase<DisplacementDim>&
            fracture_material_)
        : fracture_material(fracture_material_),
          material_state_variables(
              fracture_material.createMaterialStateVariables()),
          darcy_velocity(GlobalDimVectorType::Zero())
    {
    }
 
    using GlobalDimVectorType = Eigen::Matrix<double, DisplacementDim, 1>;
 
    typename HMatricesType::HMatrixType H_u;
    typename HMatricesType::ForceVectorType sigma_eff, sigma_eff_prev;
    typename HMatricesType::ForceVectorType w, w_prev;
 
    typename ShapeMatrixTypePressure::NodalRowVectorType N_p;
    typename ShapeMatrixTypePressure::GlobalDimNodalMatrixType dNdx_p;
 
    double aperture = 0.0;
    double aperture0 = 0.0;
    double permeability = 0.0;
 
    MaterialLib::Fracture::FractureModelBase<DisplacementDim>&
        fracture_material;
    std::unique_ptr<typename MaterialLib::Fracture::FractureModelBase<
        DisplacementDim>::MaterialStateVariables>
        material_state_variables;
 
    Eigen::MatrixXd C;
    double integration_weight;
 
    GlobalDimVectorType darcy_velocity;
 
    void pushBackState()
    {
        w_prev = w;
        sigma_eff_prev = sigma_eff;
        material_state_variables->pushBackState();
    }
 
    EIGEN_MAKE_ALIGNED_OPERATOR_NEW;
};
 
}  // namespace HydroMechanics
}  // namespace LIE
}  // namespace ProcessLib