|
| v0.14.0
|
|
| OpSpatialPhysical (const std::string &field_name, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const double alpha_u) |
|
MoFEMErrorCode | integrate (EntData &data) |
|
| OpAssembleVolume (const std::string &field, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const char type) |
|
| OpAssembleVolume (std::string row_field, std::string col_field, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const char type, const bool assemble_symmetry) |
|
| OpAssembleBasic (const std::string &field_name, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const char type) |
|
| OpAssembleBasic (std::string row_field, std::string col_field, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const char type, const bool assemble_symmetry, ScaleOff scale_off=[]() { return 1;}) |
|
virtual MoFEMErrorCode | integrate (int row_side, EntityType row_type, EntData &data) |
|
virtual MoFEMErrorCode | integrate (EntData &row_data, EntData &col_data) |
|
virtual MoFEMErrorCode | assemble (EntData &data) |
|
virtual MoFEMErrorCode | assemble (int row_side, EntityType row_type, EntData &data) |
|
virtual MoFEMErrorCode | assemble (int row_side, int col_side, EntityType row_type, EntityType col_type, EntData &row_data, EntData &col_data) |
|
MoFEMErrorCode | doWork (int side, EntityType type, EntData &data) |
|
MoFEMErrorCode | doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntData &row_data, EntData &col_data) |
|
Definition at line 57 of file HMHNeohookean.cpp.
◆ OpSpatialPhysical()
EshelbianPlasticity::HMHNeohookean::OpSpatialPhysical::OpSpatialPhysical |
( |
const std::string & |
field_name, |
|
|
boost::shared_ptr< DataAtIntegrationPts > |
data_ptr, |
|
|
const double |
alpha_u |
|
) |
| |
◆ integrate()
MoFEMErrorCode EshelbianPlasticity::HMHNeohookean::OpSpatialPhysical::integrate |
( |
EntData & |
data | ) |
|
|
virtual |
Reimplemented from EshelbianPlasticity::OpAssembleBasic< VolUserDataOperator >.
Definition at line 99 of file HMHNeohookean.cpp.
102 auto neohookean_ptr =
103 boost::dynamic_pointer_cast<HMHNeohookean>(
dataAtPts->physicsPtr);
104 if (!neohookean_ptr) {
106 "Pointer to HMHNeohookean is null");
112 "Stretch selector is not equal to LOG");
116 "Exponent base is not equal to exp(1)");
121 const auto c10 = neohookean_ptr->c10;
122 const auto K = neohookean_ptr->K;
131 int nb_integration_pts = data.
getN().size1();
132 auto v = getVolume();
133 auto t_w = getFTensor0IntegrationWeight();
134 auto t_approx_P_adjont_log_du =
135 getFTensor1FromMat<size_symm>(
dataAtPts->adjointPdUAtPts);
136 auto t_u = getFTensor2SymmetricFromMat<3>(
dataAtPts->stretchTensorAtPts);
137 auto t_grad_h1 = getFTensor2FromMat<3, 3>(
dataAtPts->wGradH1AtPts);
139 getFTensor2SymmetricFromMat<3>(
dataAtPts->logStretchTotalTensorAtPts);
141 getFTensor2SymmetricFromMat<3>(
dataAtPts->logStretchDotTensorAtPts);
143 getFTensor4DdgFromMat<3, 3, 1>(
dataAtPts->diffStretchTensorAtPts);
145 getFTensor2SymmetricFromMat<3>(
dataAtPts->logStretchTensorAtPts);
147 getFTensor2SymmetricFromMat<3>(
dataAtPts->logStretch2H1AtPts);
158 auto get_ftensor2 = [](
auto &
v) {
160 &
v[0], &
v[1], &
v[2], &
v[3], &
v[4], &
v[5]);
163 int nb_base_functions = data.
getN().size2();
165 for (
int gg = 0; gg != nb_integration_pts; ++gg) {
176 t_Ldiff_u(
i,
j,
L) = t_diff_u(
i,
j,
k,
l) * t_L(
k,
l,
L);
182 t_Ldiff_u(
i,
j,
L) = (t_diff_u(
i,
m,
k,
l) * t_h1(
m,
j)) * t_L(
k,
l,
L);
187 t_Ldiff_u(
i,
j,
L) = t_diff_u(
i,
j,
k,
l) * t_L(
k,
l,
L);
194 t_Sigma_u(
i,
j) = 2.0 *
c10 * (t_u(
i,
m) * t_h1(
m,
j));
195 const double tr = t_total_log_u(
i,
j) * t_kd_sym(
i,
j);
197 const double Simga_J = -2 *
c10 +
K * (
J - 1) *
J;
206 t_residual(
L) = t_approx_P_adjont_log_du(
L);
207 t_residual(
L) -= t_Ldiff_u(
i,
j,
L) * t_Sigma_u(
i,
j);
208 t_residual(
L) -= (t_L(
i,
j,
L) *
t_kd(
i,
j)) * Simga_J;
209 t_residual(
L) -= t_L(
i,
j,
L) * t_viscous_P(
i,
j);
212 ++t_approx_P_adjont_log_du;
217 auto t_nf = getFTensor1FromPtr<size_symm>(&*
nF.data().begin());
219 for (; bb != nb_dofs /
size_symm; ++bb) {
220 t_nf(
L) += t_row_base_fun * t_residual(
L);
224 for (; bb != nb_base_functions; ++bb)
◆ alphaU
const double EshelbianPlasticity::HMHNeohookean::OpSpatialPhysical::alphaU |
|
private |
The documentation for this struct was generated from the following file:
constexpr static auto size_symm
static double exponentBase
FTensor::Index< 'J', DIM1 > J
VectorDouble nF
local right hand side vector
#define FTENSOR_INDEX(DIM, I)
FTensor::Tensor0< FTensor::PackPtr< double *, 1 > > getFTensor0N(const FieldApproximationBase base)
Get base function as Tensor0.
static enum RotSelector gradApproximator
static enum StretchSelector stretchSelector
static boost::function< double(const double)> f
const VectorInt & getIndices() const
Get global indices of dofs on entity.
boost::shared_ptr< DataAtIntegrationPts > dataAtPts
data at integration pts
FTensor::Index< 'i', SPACE_DIM > i
constexpr auto field_name
const double v
phase velocity of light in medium (cm/ns)
MatrixDouble & getN(const FieldApproximationBase base)
get base functions this return matrix (nb. of rows is equal to nb. of Gauss pts, nb....
FTensor::Index< 'j', 3 > j
@ MOFEM_DATA_INCONSISTENCY
FTensor::Index< 'm', 3 > m
Kronecker Delta class symmetric.
OpAssembleVolume(const std::string &field, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const char type)
FTensor::Index< 'k', 3 > k
MatrixDouble K
local tangent matrix
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
FTensor::Index< 'l', 3 > l