OGS
List of all members | Public Types | Public Member Functions | Private Types | Private Attributes
ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim > Class Template Reference

Detailed Description

template<typename ShapeFunction, int DisplacementDim>
class ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >

Definition at line 37 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

#include <SmallDeformationLocalAssemblerMatrixNearFracture.h>

Inheritance diagram for ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >:
[legend]
Collaboration diagram for ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >:
[legend]

Public Types

using ShapeMatricesType
 
using NodalMatrixType = typename ShapeMatricesType::NodalMatrixType
 
using NodalVectorType = typename ShapeMatricesType::NodalVectorType
 
using ShapeMatrices = typename ShapeMatricesType::ShapeMatrices
 
using BMatricesType = BMatrixPolicyType<ShapeFunction, DisplacementDim>
 
using BMatrixType = typename BMatricesType::BMatrixType
 
using StiffnessMatrixType = typename BMatricesType::StiffnessMatrixType
 
using NodalForceVectorType = typename BMatricesType::NodalForceVectorType
 
using NodalDisplacementVectorType
 

Public Member Functions

 SmallDeformationLocalAssemblerMatrixNearFracture (SmallDeformationLocalAssemblerMatrixNearFracture const &)=delete
 
 SmallDeformationLocalAssemblerMatrixNearFracture (SmallDeformationLocalAssemblerMatrixNearFracture &&)=delete
 
 SmallDeformationLocalAssemblerMatrixNearFracture (MeshLib::Element const &e, std::size_t const n_variables, std::size_t const local_matrix_size, std::vector< unsigned > const &dofIndex_to_localIndex, NumLib::GenericIntegrationMethod const &integration_method, bool const is_axially_symmetric, SmallDeformationProcessData< DisplacementDim > &process_data)
 
void assemble (double const, double const, std::vector< double > const &, std::vector< double > const &, std::vector< double > &, std::vector< double > &, std::vector< double > &) override
 
void assembleWithJacobian (double const t, double const dt, Eigen::VectorXd const &local_u, Eigen::VectorXd &local_b, Eigen::MatrixXd &local_J) override
 
void preTimestepConcrete (std::vector< double > const &, double const, double const) override
 
void computeSecondaryVariableConcreteWithVector (double const t, Eigen::VectorXd const &local_x) override
 
Eigen::Map< const Eigen::RowVectorXd > getShapeMatrix (const unsigned integration_point) const override
 Provides the shape matrix at the given integration point.
 
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 & getIntPtEpsilon (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 & getIntPtFractureStress (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
 
std::vector< double > const & getIntPtFractureAperture (const double, std::vector< GlobalVector * > const &, std::vector< NumLib::LocalToGlobalIndexMap const * > const &, std::vector< double > &cache) const override
 
- Public Member Functions inherited from ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface
 SmallDeformationLocalAssemblerInterface ()
 
 SmallDeformationLocalAssemblerInterface (std::size_t n_local_size, std::vector< unsigned > dofIndex_to_localIndex)
 
void assembleWithJacobian (double const t, double const dt, std::vector< double > const &local_x_, std::vector< double > const &, std::vector< double > &, std::vector< double > &, std::vector< double > &local_b_data, std::vector< double > &local_Jac_data) override
 
void computeSecondaryVariableConcrete (double const t, double const, Eigen::VectorXd const &local_x, Eigen::VectorXd const &) override
 
- Public Member Functions inherited from ProcessLib::LocalAssemblerInterface
virtual ~LocalAssemblerInterface ()=default
 
virtual void setInitialConditions (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, double const t, int const process_id)
 
virtual void initialize (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table)
 
virtual void preAssemble (double const, double const, std::vector< double > const &)
 
virtual void assembleForStaggeredScheme (double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_x_prev, int const process_id, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data)
 
virtual void assembleWithJacobianForStaggeredScheme (double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_x_prev, int const process_id, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data, std::vector< double > &local_Jac_data)
 
virtual void computeSecondaryVariable (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_prev, int const process_id)
 
virtual void preTimestep (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table, GlobalVector const &x, double const t, double const delta_t)
 
virtual void postTimestep (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, double const t, double const dt, int const process_id)
 
void postNonLinearSolver (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, double const t, double const dt, int const process_id)
 
virtual Eigen::Vector3d getFlux (MathLib::Point3d const &, double const, std::vector< double > const &) const
 
virtual Eigen::Vector3d getFlux (MathLib::Point3d const &, double const, std::vector< std::vector< double > > const &) const
 Fits to staggered scheme.
 
- Public Member Functions inherited from NumLib::ExtrapolatableElement
virtual ~ExtrapolatableElement ()=default
 

Private Types

using IntegrationPointDataType
 

Private Attributes

SmallDeformationProcessData< DisplacementDim > & _process_data
 
std::vector< FractureProperty * > _fracture_props
 
std::vector< JunctionProperty * > _junction_props
 
std::unordered_map< int, int > _fracID_to_local
 
std::vector< IntegrationPointDataType, Eigen::aligned_allocator< IntegrationPointDataType > > _ip_data
 
NumLib::GenericIntegrationMethod const & _integration_method
 
MeshLib::Element const & _element
 
bool const _is_axially_symmetric
 
SecondaryData< typename ShapeMatrices::ShapeType > _secondary_data
 

Additional Inherited Members

Member Typedef Documentation

◆ BMatricesType

template<typename ShapeFunction , int DisplacementDim>
using ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::BMatricesType = BMatrixPolicyType<ShapeFunction, DisplacementDim>

Definition at line 46 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

◆ BMatrixType

template<typename ShapeFunction , int DisplacementDim>
using ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::BMatrixType = typename BMatricesType::BMatrixType

Definition at line 48 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

◆ IntegrationPointDataType

template<typename ShapeFunction , int DisplacementDim>
using ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::IntegrationPointDataType
private
Initial value:
IntegrationPointDataMatrix<ShapeMatricesType, BMatricesType,
DisplacementDim>

Definition at line 148 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

◆ NodalDisplacementVectorType

template<typename ShapeFunction , int DisplacementDim>
using ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::NodalDisplacementVectorType
Initial value:
typename BMatricesType::NodalForceVectorType
ProcessLib::BMatrixPolicyType::NodalForceVectorType
VectorType< _number_of_dof > NodalForceVectorType
Rhs residual.
Definition BMatrixPolicy.h:44

Definition at line 51 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

◆ NodalForceVectorType

template<typename ShapeFunction , int DisplacementDim>
using ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::NodalForceVectorType = typename BMatricesType::NodalForceVectorType

Definition at line 50 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

◆ NodalMatrixType

template<typename ShapeFunction , int DisplacementDim>
using ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::NodalMatrixType = typename ShapeMatricesType::NodalMatrixType

Definition at line 43 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

◆ NodalVectorType

template<typename ShapeFunction , int DisplacementDim>
using ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::NodalVectorType = typename ShapeMatricesType::NodalVectorType

Definition at line 44 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

◆ ShapeMatrices

template<typename ShapeFunction , int DisplacementDim>
using ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::ShapeMatrices = typename ShapeMatricesType::ShapeMatrices

Definition at line 45 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

◆ ShapeMatricesType

template<typename ShapeFunction , int DisplacementDim>
using ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::ShapeMatricesType
Initial value:
ShapeMatrixPolicyType<ShapeFunction, DisplacementDim>

Definition at line 41 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

◆ StiffnessMatrixType

template<typename ShapeFunction , int DisplacementDim>
using ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::StiffnessMatrixType = typename BMatricesType::StiffnessMatrixType

Definition at line 49 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

Constructor & Destructor Documentation

◆ SmallDeformationLocalAssemblerMatrixNearFracture() [1/3]

template<typename ShapeFunction , int DisplacementDim>
ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::SmallDeformationLocalAssemblerMatrixNearFracture ( SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim > const & )
delete

◆ SmallDeformationLocalAssemblerMatrixNearFracture() [2/3]

template<typename ShapeFunction , int DisplacementDim>
ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::SmallDeformationLocalAssemblerMatrixNearFracture ( SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim > && )
delete

◆ SmallDeformationLocalAssemblerMatrixNearFracture() [3/3]

template<typename ShapeFunction , int DisplacementDim>
ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::SmallDeformationLocalAssemblerMatrixNearFracture ( MeshLib::Element const & e,
std::size_t const n_variables,
std::size_t const local_matrix_size,
std::vector< unsigned > const & dofIndex_to_localIndex,
NumLib::GenericIntegrationMethod const & integration_method,
bool const is_axially_symmetric,
SmallDeformationProcessData< DisplacementDim > & process_data )

Definition at line 47 of file SmallDeformationLocalAssemblerMatrixNearFracture-impl.h.

56 : SmallDeformationLocalAssemblerInterface(
57 n_variables * ShapeFunction::NPOINTS * DisplacementDim,
58 dofIndex_to_localIndex),
59 _process_data(process_data),
60 _integration_method(integration_method),
61 _element(e),
62 _is_axially_symmetric(is_axially_symmetric)
63{
64 std::vector<ShapeMatrices, Eigen::aligned_allocator<
65 typename ShapeMatricesType::ShapeMatrices>>
66 shape_matrices =
67 NumLib::initShapeMatrices<ShapeFunction, ShapeMatricesType,
68 DisplacementDim>(e, is_axially_symmetric,
69 _integration_method);
70
71 unsigned const n_integration_points =
72 _integration_method.getNumberOfPoints();
73
74 _ip_data.reserve(n_integration_points);
75 _secondary_data.N.resize(n_integration_points);
76
77 auto& solid_material = MaterialLib::Solids::selectSolidConstitutiveRelation(
78 _process_data.solid_materials, _process_data.material_ids, e.getID());
79
80 for (unsigned ip = 0; ip < n_integration_points; ip++)
81 {
82 _ip_data.emplace_back(solid_material);
83 auto& ip_data = _ip_data[ip];
84 auto const& sm = shape_matrices[ip];
85 ip_data.N = sm.N;
86 ip_data.dNdx = sm.dNdx;
87 ip_data.integration_weight =
88 _integration_method.getWeightedPoint(ip).getWeight() *
89 sm.integralMeasure * sm.detJ;
90
91 // Initialize current time step values
92 static const int kelvin_vector_size =
93 MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim);
94 ip_data._sigma.setZero(kelvin_vector_size);
95 ip_data._eps.setZero(kelvin_vector_size);
96
97 // Previous time step values are not initialized and are set later.
98 ip_data._sigma_prev.resize(kelvin_vector_size);
99 ip_data._eps_prev.resize(kelvin_vector_size);
100
101 ip_data._C.resize(kelvin_vector_size, kelvin_vector_size);
102
103 _secondary_data.N[ip] = sm.N;
104 }
105
106 for (auto fid : process_data._vec_ele_connected_fractureIDs[e.getID()])
107 {
108 _fracID_to_local.insert({fid, _fracture_props.size()});
109 _fracture_props.push_back(&_process_data.fracture_properties[fid]);
110 }
111
112 _junction_props = process_data._vec_ele_connected_junctionIDs[e.getID()] |
113 ranges::views::transform(
114 [&](auto const jid)
115 { return &_process_data.junction_properties[jid]; }) |
116 ranges::to<std::vector>;
117}
MathLib::WeightedPoint::getWeight
double getWeight() const
Definition WeightedPoint.h:80
NumLib::GenericIntegrationMethod::getWeightedPoint
MathLib::WeightedPoint const & getWeightedPoint(unsigned const igp) const
Definition GenericIntegrationMethod.h:36
ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture::_ip_data
std::vector< IntegrationPointDataType, Eigen::aligned_allocator< IntegrationPointDataType > > _ip_data
Definition SmallDeformationLocalAssemblerMatrixNearFracture.h:154
MaterialLib::Solids::selectSolidConstitutiveRelation
auto & selectSolidConstitutiveRelation(SolidMaterialsMap const &constitutive_relations, MeshLib::PropertyVector< int > const *const material_ids, std::size_t const element_id)
Definition SelectSolidConstitutiveRelation.h:36
MathLib::KelvinVector::kelvin_vector_dimensions
constexpr int kelvin_vector_dimensions(int const displacement_dim)
Kelvin vector dimensions for given displacement dimension.
Definition KelvinVector.h:24
NumLib::initShapeMatrices
std::vector< typename ShapeMatricesType::ShapeMatrices, Eigen::aligned_allocator< typename ShapeMatricesType::ShapeMatrices > > initShapeMatrices(MeshLib::Element const &e, bool const is_axially_symmetric, IntegrationMethod const &integration_method)
Definition InitShapeMatrices.h:53
EigenFixedShapeMatrixPolicy::ShapeMatrices
NumLib::ShapeMatrices< NodalRowVectorType, DimNodalMatrixType, DimMatrixType, GlobalDimNodalMatrixType > ShapeMatrices
Definition ShapeMatrixPolicy.h:87
ProcessLib::HeatTransportBHE::SecondaryData::N
std::vector< ShapeMatrixType, Eigen::aligned_allocator< ShapeMatrixType > > N
Definition SecondaryData.h:27

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::_fracID_to_local, ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::_fracture_props, ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::_integration_method, ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::_ip_data, ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::_junction_props, ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::_process_data, ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::_secondary_data, ProcessLib::LIE::SmallDeformation::SmallDeformationProcessData< DisplacementDim >::_vec_ele_connected_fractureIDs, ProcessLib::LIE::SmallDeformation::SmallDeformationProcessData< DisplacementDim >::_vec_ele_connected_junctionIDs, MeshLib::Element::getID(), NumLib::GenericIntegrationMethod::getNumberOfPoints(), MathLib::WeightedPoint::getWeight(), NumLib::GenericIntegrationMethod::getWeightedPoint(), NumLib::initShapeMatrices(), MathLib::KelvinVector::kelvin_vector_dimensions(), ProcessLib::HeatTransportBHE::SecondaryData< ShapeMatrixType >::N, and MaterialLib::Solids::selectSolidConstitutiveRelation().

Member Function Documentation

◆ assemble()

template<typename ShapeFunction , int DisplacementDim>
void ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::assemble ( double const ,
double const ,
std::vector< double > const & ,
std::vector< double > const & ,
std::vector< double > & ,
std::vector< double > & ,
std::vector< double > &  )
inlineoverridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 68 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

74 {
75 OGS_FATAL(
76 "SmallDeformationLocalAssembler: assembly without jacobian is not "
77 "implemented.");
78 }
OGS_FATAL
#define OGS_FATAL(...)
Definition Error.h:26

References OGS_FATAL.

◆ assembleWithJacobian()

template<typename ShapeFunction , int DisplacementDim>
void ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::assembleWithJacobian ( double const t,
double const dt,
Eigen::VectorXd const & local_u,
Eigen::VectorXd & local_b,
Eigen::MatrixXd & local_J )
overridevirtual

Reimplemented from ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 122 of file SmallDeformationLocalAssemblerMatrixNearFracture-impl.h.

126{
127 assert(_element.getDimension() == DisplacementDim);
128
129 auto const N_DOF_PER_VAR = ShapeFunction::NPOINTS * DisplacementDim;
130 auto const n_fractures = _fracture_props.size();
131 auto const n_junctions = _junction_props.size();
132 auto const n_enrich_var = n_fractures + n_junctions;
133
134 using BlockVectorType =
135 typename Eigen::VectorXd::FixedSegmentReturnType<N_DOF_PER_VAR>::Type;
136 using BlockMatrixType =
137 Eigen::Block<Eigen::MatrixXd, N_DOF_PER_VAR, N_DOF_PER_VAR>;
138
139 //--------------------------------------------------------------------------------------
140 // prepare sub vectors, matrices for regular displacement (u) and
141 // displacement jumps (g)
142 //
143 // example with two fractures with one intersection:
144 // |b(u)|
145 // b = |b(g1)|
146 // |b(g2)|
147 // |b(j1)|
148 //
149 // |J(u,u) J(u,g1) J(u,g2) J(u,j1) |
150 // J = |J(g1,u) J(g1,g1) J(g1,g2) J(g1,j1)|
151 // |J(g2,u) J(g2,g1) J(g2,g2) J(g2,j1)|
152 // |J(j1,u) J(j1,g1) J(j1,g2) J(j1,j1)|
153 //--------------------------------------------------------------------------------------
154 auto local_b_u = local_b.segment<N_DOF_PER_VAR>(0);
155 std::vector<BlockVectorType> vec_local_b_g;
156 for (unsigned i = 0; i < n_enrich_var; i++)
157 {
158 vec_local_b_g.push_back(
159 local_b.segment<N_DOF_PER_VAR>(N_DOF_PER_VAR * (i + 1)));
160 }
161
162 auto local_J_uu = local_J.block<N_DOF_PER_VAR, N_DOF_PER_VAR>(0, 0);
163 std::vector<BlockMatrixType> vec_local_J_ug;
164 std::vector<BlockMatrixType> vec_local_J_gu;
165 std::vector<std::vector<BlockMatrixType>> vec_local_J_gg(n_enrich_var);
166 for (unsigned i = 0; i < n_enrich_var; i++)
167 {
168 auto sub_ug = local_J.block<N_DOF_PER_VAR, N_DOF_PER_VAR>(
169 0, N_DOF_PER_VAR * (i + 1));
170 vec_local_J_ug.push_back(sub_ug);
171
172 auto sub_gu = local_J.block<N_DOF_PER_VAR, N_DOF_PER_VAR>(
173 N_DOF_PER_VAR * (i + 1), 0);
174 vec_local_J_gu.push_back(sub_gu);
175
176 for (unsigned j = 0; j < n_enrich_var; j++)
177 {
178 auto sub_gg = local_J.block<N_DOF_PER_VAR, N_DOF_PER_VAR>(
179 N_DOF_PER_VAR * (i + 1), N_DOF_PER_VAR * (j + 1));
180 vec_local_J_gg[i].push_back(sub_gg);
181 }
182 }
183
184 auto const nodal_u = local_u.segment<N_DOF_PER_VAR>(0);
185 std::vector<BlockVectorType> vec_nodal_g;
186 for (unsigned i = 0; i < n_enrich_var; i++)
187 {
188 auto sub = const_cast<Eigen::VectorXd&>(local_u).segment<N_DOF_PER_VAR>(
189 N_DOF_PER_VAR * (i + 1));
190 vec_nodal_g.push_back(sub);
191 }
192
193 //------------------------------------------------
194 // integration
195 //------------------------------------------------
196 unsigned const n_integration_points =
197 _integration_method.getNumberOfPoints();
198
199 MPL::VariableArray variables;
200 MPL::VariableArray variables_prev;
201 ParameterLib::SpatialPosition x_position;
202 x_position.setElementID(_element.getID());
203
204 for (unsigned ip = 0; ip < n_integration_points; ip++)
205 {
206 x_position.setIntegrationPoint(ip);
207
208 auto& ip_data = _ip_data[ip];
209 auto const& w = _ip_data[ip].integration_weight;
210
211 auto const& N = ip_data.N;
212 auto const& dNdx = ip_data.dNdx;
213
214 // levelset functions
215 Eigen::Vector3d const ip_physical_coords(
216 computePhysicalCoordinates(_element, N).data());
217 std::vector<double> const levelsets(
218 uGlobalEnrichments(_fracture_props, _junction_props,
219 _fracID_to_local, ip_physical_coords));
220
221 // u = u^hat + sum_i(enrich^br_i(x) * [u]_i) + sum_i(enrich^junc_i(x) *
222 // [u]_i)
223 NodalDisplacementVectorType nodal_total_u = nodal_u;
224 for (unsigned i = 0; i < n_enrich_var; i++)
225 {
226 nodal_total_u += levelsets[i] * vec_nodal_g[i];
227 }
228
229 auto const x_coord =
230 NumLib::interpolateXCoordinate<ShapeFunction, ShapeMatricesType>(
231 _element, N);
232 auto const B =
233 LinearBMatrix::computeBMatrix<DisplacementDim,
234 ShapeFunction::NPOINTS,
235 typename BMatricesType::BMatrixType>(
236 dNdx, N, x_coord, _is_axially_symmetric);
237
238 // strain, stress
239 auto const& eps_prev = ip_data._eps_prev;
240 auto const& sigma_prev = ip_data._sigma_prev;
241
242 auto& eps = ip_data._eps;
243 auto& sigma = ip_data._sigma;
244 auto& state = ip_data._material_state_variables;
245
246 eps.noalias() = B * nodal_total_u;
247
248 variables.mechanical_strain
249 .emplace<MathLib::KelvinVector::KelvinVectorType<DisplacementDim>>(
250 eps);
251
252 variables_prev.stress
253 .emplace<MathLib::KelvinVector::KelvinVectorType<DisplacementDim>>(
254 sigma_prev);
255 variables_prev.mechanical_strain
256 .emplace<MathLib::KelvinVector::KelvinVectorType<DisplacementDim>>(
257 eps_prev);
258 variables_prev.temperature = _process_data._reference_temperature;
259
260 auto&& solution = _ip_data[ip]._solid_material.integrateStress(
261 variables_prev, variables, t, x_position, dt, *state);
262
263 if (!solution)
264 {
265 OGS_FATAL("Computation of local constitutive relation failed.");
266 }
267
268 MathLib::KelvinVector::KelvinMatrixType<DisplacementDim> C;
269 std::tie(sigma, state, C) = std::move(*solution);
270
271 // r_u = B^T * Sigma = B^T * C * B * (u+phi*[u])
272 // r_[u] = (phi*B)^T * Sigma = (phi*B)^T * C * B * (u+phi*[u])
273 local_b_u.noalias() -= B.transpose() * sigma * w;
274 for (unsigned i = 0; i < n_enrich_var; i++)
275 {
276 vec_local_b_g[i].noalias() -=
277 levelsets[i] * B.transpose() * sigma * w;
278 }
279
280 // J_uu += B^T * C * B
281 local_J_uu.noalias() += B.transpose() * C * B * w;
282
283 for (unsigned i = 0; i < n_enrich_var; i++)
284 {
285 // J_u[u] += B^T * C * (levelset * B)
286 vec_local_J_ug[i].noalias() +=
287 B.transpose() * C * (levelsets[i] * B) * w;
288
289 // J_[u]u += (levelset * B)^T * C * B
290 vec_local_J_gu[i].noalias() +=
291 (levelsets[i] * B.transpose()) * C * B * w;
292
293 for (unsigned j = 0; j < n_enrich_var; j++)
294 {
295 // J_[u][u] += (levelset * B)^T * C * (levelset * B)
296 vec_local_J_gg[i][j].noalias() +=
297 (levelsets[i] * B.transpose()) * C * (levelsets[j] * B) * w;
298 }
299 }
300 }
301}
MaterialPropertyLib::VariableArray::mechanical_strain
std::variant< std::monostate, Eigen::Matrix< double, 4, 1 >, Eigen::Matrix< double, 6, 1 > > mechanical_strain
Definition VariableType.h:178
MaterialPropertyLib::VariableArray::stress
std::variant< std::monostate, Eigen::Matrix< double, 4, 1 >, Eigen::Matrix< double, 6, 1 > > stress
Definition VariableType.h:188
MeshLib::Element::getDimension
virtual constexpr unsigned getDimension() const =0
Get dimension of the mesh element.
MeshLib::Element::getID
std::size_t getID() const
Returns the ID of the element.
Definition Element.h:89
ParameterLib::SpatialPosition::setElementID
void setElementID(std::size_t element_id)
Definition SpatialPosition.h:59
ParameterLib::SpatialPosition::setIntegrationPoint
void setIntegrationPoint(unsigned integration_point)
Definition SpatialPosition.h:65
ProcessLib::BMatrixPolicyType::BMatrixType
MatrixType< _kelvin_vector_size, _number_of_dof > BMatrixType
Definition BMatrixPolicy.h:55
MathLib::KelvinVector::KelvinVectorType
Eigen::Matrix< double, kelvin_vector_dimensions(DisplacementDim), 1, Eigen::ColMajor > KelvinVectorType
Definition KelvinVector.h:47
MathLib::KelvinVector::KelvinMatrixType
Eigen::Matrix< double, kelvin_vector_dimensions(DisplacementDim), kelvin_vector_dimensions(DisplacementDim), Eigen::RowMajor > KelvinMatrixType
Definition KelvinVector.h:55
ProcessLib::LIE::computePhysicalCoordinates
Eigen::Vector3d computePhysicalCoordinates(MeshLib::Element const &e, Eigen::MatrixBase< Derived > const &shape)
Definition Utils.h:24
ProcessLib::LIE::uGlobalEnrichments
std::vector< double > uGlobalEnrichments(std::vector< FractureProperty * > const &frac_props, std::vector< JunctionProperty * > const &junction_props, std::unordered_map< int, int > const &fracID_to_local, Eigen::Vector3d const &x)
Definition LevelSetFunction.cpp:44
ProcessLib::LinearBMatrix::computeBMatrix
BMatrixType computeBMatrix(DNDX_Type const &dNdx, N_Type const &N, const double radius, const bool is_axially_symmetric)
Fills a B-matrix based on given shape function dN/dx values.
Definition LinearBMatrix.h:42

References ProcessLib::LinearBMatrix::computeBMatrix(), ProcessLib::LIE::computePhysicalCoordinates(), MaterialPropertyLib::VariableArray::mechanical_strain, OGS_FATAL, ParameterLib::SpatialPosition::setElementID(), ParameterLib::SpatialPosition::setIntegrationPoint(), MaterialPropertyLib::VariableArray::stress, MaterialPropertyLib::VariableArray::temperature, and ProcessLib::LIE::uGlobalEnrichments().

◆ computeSecondaryVariableConcreteWithVector()

template<typename ShapeFunction , int DisplacementDim>
void ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::computeSecondaryVariableConcreteWithVector ( double const t,
Eigen::VectorXd const & local_x )
overridevirtual

Implements ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 308 of file SmallDeformationLocalAssemblerMatrixNearFracture-impl.h.

311{
312 // Compute average value per element
313 using KV = MathLib::KelvinVector::KelvinVectorType<DisplacementDim>;
314 KV sigma_avg = KV::Zero();
315 auto const e_id = _element.getID();
316
317 unsigned const n_integration_points =
318 _integration_method.getNumberOfPoints();
319 for (unsigned ip = 0; ip < n_integration_points; ip++)
320 {
321 sigma_avg += _ip_data[ip]._sigma;
322 }
323 sigma_avg /= n_integration_points;
324
325 Eigen::Map<KV>(
326 &(*_process_data.element_stresses)[e_id * KV::RowsAtCompileTime]) =
327 MathLib::KelvinVector::kelvinVectorToSymmetricTensor(sigma_avg);
328}
MathLib::KelvinVector::kelvinVectorToSymmetricTensor
Eigen::Matrix< double, 4, 1 > kelvinVectorToSymmetricTensor(Eigen::Matrix< double, 4, 1, Eigen::ColMajor, 4, 1 > const &v)
Definition KelvinVector.cpp:142

References MathLib::KelvinVector::kelvinVectorToSymmetricTensor().

◆ getIntPtEpsilon()

template<typename ShapeFunction , int DisplacementDim>
std::vector< double > const & ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::getIntPtEpsilon ( const double t,
std::vector< GlobalVector * > const & x,
std::vector< NumLib::LocalToGlobalIndexMap const * > const & dof_table,
std::vector< double > & cache ) const
overridevirtual

Implements ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 344 of file SmallDeformationLocalAssemblerMatrixNearFracture-impl.h.

350{
351 return ProcessLib::getIntegrationPointKelvinVectorData<DisplacementDim>(
352 _ip_data, &IntegrationPointDataType::_eps, cache);
353}

◆ getIntPtFractureAperture()

template<typename ShapeFunction , int DisplacementDim>
std::vector< double > const & ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::getIntPtFractureAperture ( const double ,
std::vector< GlobalVector * > const & ,
std::vector< NumLib::LocalToGlobalIndexMap const * > const & ,
std::vector< double > & cache ) const
inlineoverridevirtual

Implements ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 132 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

137 {
138 cache.resize(0);
139 return cache;
140 }

◆ getIntPtFractureStress()

template<typename ShapeFunction , int DisplacementDim>
std::vector< double > const & ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::getIntPtFractureStress ( const double ,
std::vector< GlobalVector * > const & ,
std::vector< NumLib::LocalToGlobalIndexMap const * > const & ,
std::vector< double > & cache ) const
inlineoverridevirtual

Implements ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 122 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

127 {
128 cache.resize(0);
129 return cache;
130 }

◆ getIntPtSigma()

template<typename ShapeFunction , int DisplacementDim>
std::vector< double > const & ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::getIntPtSigma ( const double t,
std::vector< GlobalVector * > const & x,
std::vector< NumLib::LocalToGlobalIndexMap const * > const & dof_table,
std::vector< double > & cache ) const
overridevirtual

Implements ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerInterface.

Definition at line 332 of file SmallDeformationLocalAssemblerMatrixNearFracture-impl.h.

338{
339 return ProcessLib::getIntegrationPointKelvinVectorData<DisplacementDim>(
340 _ip_data, &IntegrationPointDataType::_sigma, cache);
341}

◆ getShapeMatrix()

template<typename ShapeFunction , int DisplacementDim>
Eigen::Map< const Eigen::RowVectorXd > ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::getShapeMatrix ( const unsigned integration_point) const
inlineoverridevirtual

Provides the shape matrix at the given integration point.

Implements NumLib::ExtrapolatableElement.

Definition at line 101 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

103 {
104 auto const& N = _secondary_data.N[integration_point];
105
106 // assumes N is stored contiguously in memory
107 return Eigen::Map<const Eigen::RowVectorXd>(N.data(), N.size());
108 }

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::_secondary_data, and ProcessLib::HeatTransportBHE::SecondaryData< ShapeMatrixType >::N.

◆ preTimestepConcrete()

template<typename ShapeFunction , int DisplacementDim>
void ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::preTimestepConcrete ( std::vector< double > const & ,
double const ,
double const  )
inlineoverridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 85 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

88 {
89 unsigned const n_integration_points =
90 _integration_method.getNumberOfPoints();
91
92 for (unsigned ip = 0; ip < n_integration_points; ip++)
93 {
94 _ip_data[ip].pushBackState();
95 }
96 }

References ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::_integration_method, ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::_ip_data, and NumLib::GenericIntegrationMethod::getNumberOfPoints().

Member Data Documentation

◆ _element

template<typename ShapeFunction , int DisplacementDim>
MeshLib::Element const& ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::_element
private

Definition at line 158 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

◆ _fracID_to_local

template<typename ShapeFunction , int DisplacementDim>
std::unordered_map<int, int> ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::_fracID_to_local
private

Definition at line 146 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

Referenced by ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::SmallDeformationLocalAssemblerMatrixNearFracture().

◆ _fracture_props

template<typename ShapeFunction , int DisplacementDim>
std::vector<FractureProperty*> ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::_fracture_props
private

Definition at line 144 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

Referenced by ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::SmallDeformationLocalAssemblerMatrixNearFracture().

◆ _integration_method

template<typename ShapeFunction , int DisplacementDim>
NumLib::GenericIntegrationMethod const& ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::_integration_method
private

Definition at line 156 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

Referenced by ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::SmallDeformationLocalAssemblerMatrixNearFracture(), and ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::preTimestepConcrete().

◆ _ip_data

template<typename ShapeFunction , int DisplacementDim>
std::vector<IntegrationPointDataType, Eigen::aligned_allocator<IntegrationPointDataType> > ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::_ip_data
private

Definition at line 154 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

Referenced by ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::SmallDeformationLocalAssemblerMatrixNearFracture(), and ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::preTimestepConcrete().

◆ _is_axially_symmetric

template<typename ShapeFunction , int DisplacementDim>
bool const ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::_is_axially_symmetric
private

Definition at line 159 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

◆ _junction_props

template<typename ShapeFunction , int DisplacementDim>
std::vector<JunctionProperty*> ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::_junction_props
private

Definition at line 145 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

Referenced by ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::SmallDeformationLocalAssemblerMatrixNearFracture().

◆ _process_data

template<typename ShapeFunction , int DisplacementDim>
SmallDeformationProcessData<DisplacementDim>& ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::_process_data
private

Definition at line 143 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

Referenced by ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::SmallDeformationLocalAssemblerMatrixNearFracture().

◆ _secondary_data

template<typename ShapeFunction , int DisplacementDim>
SecondaryData<typename ShapeMatrices::ShapeType> ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::_secondary_data
private

Definition at line 160 of file SmallDeformationLocalAssemblerMatrixNearFracture.h.

Referenced by ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::SmallDeformationLocalAssemblerMatrixNearFracture(), and ProcessLib::LIE::SmallDeformation::SmallDeformationLocalAssemblerMatrixNearFracture< ShapeFunction, DisplacementDim >::getShapeMatrix().

The documentation for this class was generated from the following files: