https://mooseframework.inl.gov
Loading...
Searching...
No Matches
Public Types | Public Member Functions | Static Public Member Functions | Public Attributes | Static Public Attributes | Protected Member Functions | Static Protected Member Functions | Protected Attributes | Static Protected Attributes | Private Member Functions | Static Private Member Functions | Private Attributes | List of all members
SolidMechanicsPlasticMeanCapTC Class Reference

Rate-independent associative mean-cap tensile AND compressive failure with hardening/softening of the tensile and compressive strength. More...

#include <SolidMechanicsPlasticMeanCapTC.h>

Inheritance diagram for SolidMechanicsPlasticMeanCapTC:
[legend]

Public Types

typedef DataFileName DataFileParameterType
 

Public Member Functions

 SolidMechanicsPlasticMeanCapTC (const InputParameters &parameters)
 
virtual void activeConstraints (const std::vector< Real > &f, const RankTwoTensor &stress, Real intnl, const RankFourTensor &Eijkl, std::vector< bool > &act, RankTwoTensor &returned_stress) const override
 The active yield surfaces, given a vector of yield functions.
 
virtual bool useCustomReturnMap () const override
 Returns false. You will want to override this in your derived class if you write a custom returnMap function.
 
virtual bool useCustomCTO () const override
 Returns false. You will want to override this in your derived class if you write a custom consistent tangent operator function.
 
virtual bool returnMap (const RankTwoTensor &trial_stress, Real intnl_old, const RankFourTensor &E_ijkl, Real ep_plastic_tolerance, RankTwoTensor &returned_stress, Real &returned_intnl, std::vector< Real > &dpm, RankTwoTensor &delta_dp, std::vector< Real > &yf, bool &trial_stress_inadmissible) const override
 Performs a custom return-map.
 
virtual RankFourTensor consistentTangentOperator (const RankTwoTensor &trial_stress, Real intnl_old, const RankTwoTensor &stress, Real intnl, const RankFourTensor &E_ijkl, const std::vector< Real > &cumulative_pm) const override
 Calculates a custom consistent tangent operator.
 
virtual std::string modelName () const override
 
void initialize ()
 
void execute ()
 
void finalize ()
 
virtual unsigned int numberSurfaces () const
 The number of yield surfaces for this plasticity model.
 
virtual void yieldFunctionV (const RankTwoTensor &stress, Real intnl, std::vector< Real > &f) const
 Calculates the yield functions.
 
virtual void dyieldFunction_dstressV (const RankTwoTensor &stress, Real intnl, std::vector< RankTwoTensor > &df_dstress) const
 The derivative of yield functions with respect to stress.
 
virtual void dyieldFunction_dintnlV (const RankTwoTensor &stress, Real intnl, std::vector< Real > &df_dintnl) const
 The derivative of yield functions with respect to the internal parameter.
 
virtual void flowPotentialV (const RankTwoTensor &stress, Real intnl, std::vector< RankTwoTensor > &r) const
 The flow potentials.
 
virtual void dflowPotential_dstressV (const RankTwoTensor &stress, Real intnl, std::vector< RankFourTensor > &dr_dstress) const
 The derivative of the flow potential with respect to stress.
 
virtual void dflowPotential_dintnlV (const RankTwoTensor &stress, Real intnl, std::vector< RankTwoTensor > &dr_dintnl) const
 The derivative of the flow potential with respect to the internal parameter.
 
virtual void hardPotentialV (const RankTwoTensor &stress, Real intnl, std::vector< Real > &h) const
 The hardening potential.
 
virtual void dhardPotential_dstressV (const RankTwoTensor &stress, Real intnl, std::vector< RankTwoTensor > &dh_dstress) const
 The derivative of the hardening potential with respect to stress.
 
virtual void dhardPotential_dintnlV (const RankTwoTensor &stress, Real intnl, std::vector< Real > &dh_dintnl) const
 The derivative of the hardening potential with respect to the internal parameter.
 
bool KuhnTuckerSingleSurface (Real yf, Real dpm, Real dpm_tol) const
 Returns true if the Kuhn-Tucker conditions for the single surface are satisfied.
 
virtual Real spatialValue (const Point &) const
 
virtual const std::vector< Point > spatialPoints () const
 
void setPrimaryThreadCopy (UserObject *primary)
 
UserObject * primaryThreadCopy ()
 
SubProblem & getSubProblem () const
 
bool shouldDuplicateInitialExecution () const
 
void gatherSum (T &value)
 
void gatherMax (T &value)
 
void gatherMin (T &value)
 
void gatherProxyValueMax (T1 &proxy, T2 &value)
 
void gatherProxyValueMin (T1 &proxy, T2 &value)
 
std::set< UserObjectName > getDependObjects () const
 
const std::set< std::string > & getRequestedItems () override
 
const std::set< std::string > & getSuppliedItems () override
 
unsigned int systemNumber () const
 
virtual bool needThreadedCopy () const
 
virtual bool enabled () const
 
std::shared_ptr< MooseObject > getSharedPtr ()
 
std::shared_ptr< const MooseObject > getSharedPtr () const
 
bool isKokkosObject () const
 
MooseApp & getMooseApp () const
 
const std::string & type () const
 
const std::string & name () const
 
std::string typeAndName () const
 
MooseObjectParameterName uniqueParameterName (const std::string &parameter_name) const
 
MooseObjectName uniqueName () const
 
const InputParameters & parameters () const
 
const hit::Node * getHitNode () const
 
bool hasBase () const
 
const std::string & getBase () const
 
const T & getParam (const std::string &name) const
 
std::vector< std::pair< T1, T2 > > getParam (const std::string &param1, const std::string &param2) const
 
const T * queryParam (const std::string &name) const
 
const T & getRenamedParam (const std::string &old_name, const std::string &new_name) const
 
T getCheckedPointerParam (const std::string &name, const std::string &error_string="") const
 
bool isParamValid (const std::string &name) const
 
bool isParamSetByUser (const std::string &name) const
 
void connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const
 
void paramError (const std::string &param, Args... args) const
 
void paramWarning (const std::string &param, Args... args) const
 
void paramWarning (const std::string &param, Args... args) const
 
void paramInfo (const std::string &param, Args... args) const
 
std::string messagePrefix (const bool hit_prefix=true) const
 
std::string errorPrefix (const std::string &) const
 
void mooseError (Args &&... args) const
 
void mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const
 
void mooseErrorNonPrefixed (Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseDeprecatedNoTrace (Args &&... args) const
 
void mooseInfo (Args &&... args) const
 
void callMooseError (std::string msg, const bool with_prefix, const hit::Node *node=nullptr, const bool show_trace=true) const
 
std::string getDataFileName (const std::string &param) const
 
std::string getDataFileNameByName (const std::string &relative_path) const
 
std::string getDataFilePath (const std::string &relative_path) const
 
virtual void initialSetup ()
 
virtual void timestepSetup ()
 
virtual void jacobianSetup ()
 
virtual void residualSetup ()
 
virtual void customSetup (const ExecFlagType &)
 
const ExecFlagEnum & getExecuteOnEnum () const
 
UserObjectName getUserObjectName (const std::string &param_name) const
 
const T & getUserObject (const std::string &param_name, bool is_dependency=true) const
 
const T & getUserObjectByName (const UserObjectName &object_name, bool is_dependency=true) const
 
const UserObjectBase & getUserObjectBase (const std::string &param_name, bool is_dependency=true) const
 
const UserObjectBase & getUserObjectBaseByName (const UserObjectName &object_name, bool is_dependency=true) const
 
const std::vector< MooseVariableScalar * > & getCoupledMooseScalarVars ()
 
const std::set< TagID > & getScalarVariableCoupleableVectorTags () const
 
const std::set< TagID > & getScalarVariableCoupleableMatrixTags () const
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialProperty (const std::string &name, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialProperty (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialProperty (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialProperty (const std::string &name, const unsigned int state=0)
 
const ADMaterialProperty< T > & getADMaterialProperty (const std::string &name, MaterialData &material_data)
 
const ADMaterialProperty< T > & getADMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > & getADMaterialProperty (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOld (const std::string &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOld (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOld (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOlder (const std::string &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOlder (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOlder (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const MaterialPropertyName &name)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const MaterialPropertyName &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const MaterialPropertyName &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const MaterialPropertyName &name)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyByName (const std::string &prop_name_in)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyOldByName (const std::string &prop_name)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyOlderByName (const std::string &prop_name)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialProperty (const std::string &name)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyOld (const std::string &name)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyOlder (const std::string &name)
 
std::pair< const MaterialProperty< T > *, std::set< SubdomainID > > getBlockMaterialProperty (const MaterialPropertyName &name)
 
std::pair< Moose::Kokkos::MaterialProperty< T, dimension >, std::set< SubdomainID > > getKokkosBlockMaterialProperty (const MaterialPropertyName &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialProperty (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialProperty ()
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialPropertyByName (const std::string &prop_name)
 
const MaterialProperty< T > & getZeroMaterialProperty (Ts... args)
 
std::set< SubdomainID > getMaterialPropertyBlocks (const std::string &name)
 
std::vector< SubdomainName > getMaterialPropertyBlockNames (const std::string &name)
 
std::set< BoundaryID > getMaterialPropertyBoundaryIDs (const std::string &name)
 
std::vector< BoundaryName > getMaterialPropertyBoundaryNames (const std::string &name)
 
void checkBlockAndBoundaryCompatibility (std::shared_ptr< MaterialBase > discrete)
 
std::unordered_map< SubdomainID, std::vector< MaterialBase * > > buildRequiredMaterials (bool allow_stateful=true)
 
void statefulPropertiesAllowed (bool)
 
virtual bool getMaterialPropertyCalled () const
 
virtual const std::unordered_set< unsigned int > & getMatPropDependencies () const
 
virtual void resolveOptionalProperties ()
 
const GenericMaterialProperty< T, is_ad > & getPossiblyConstantGenericMaterialPropertyByName (const MaterialPropertyName &prop_name, MaterialData &material_data, const unsigned int state)
 
bool isImplicit ()
 
Moose::StateArg determineState () const
 
virtual void threadJoin (const UserObject &) override
 
virtual void threadJoin (const UserObject &) override
 
virtual void subdomainSetup () override
 
virtual void subdomainSetup () override
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
const GenericOptionalMaterialProperty< T, is_ad > & getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const GenericOptionalMaterialProperty< T, is_ad > & getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const OptionalMaterialProperty< T > & getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const OptionalMaterialProperty< T > & getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const OptionalADMaterialProperty< T > & getOptionalADMaterialProperty (const std::string &name)
 
const OptionalADMaterialProperty< T > & getOptionalADMaterialProperty (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOld (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOld (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOlder (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOlder (const std::string &name)
 
MaterialBase & getMaterial (const std::string &name)
 
MaterialBase & getMaterial (const std::string &name)
 
MaterialBase & getMaterialByName (const std::string &name, bool no_warn=false)
 
MaterialBase & getMaterialByName (const std::string &name, bool no_warn=false)
 
bool hasMaterialProperty (const std::string &name)
 
bool hasMaterialProperty (const std::string &name)
 
bool hasMaterialPropertyByName (const std::string &name)
 
bool hasMaterialPropertyByName (const std::string &name)
 
bool hasADMaterialProperty (const std::string &name)
 
bool hasADMaterialProperty (const std::string &name)
 
bool hasADMaterialPropertyByName (const std::string &name)
 
bool hasADMaterialPropertyByName (const std::string &name)
 
bool hasKokkosMaterialProperty (const std::string &name)
 
bool hasKokkosMaterialProperty (const std::string &name)
 
bool hasKokkosMaterialPropertyByName (const std::string &name)
 
bool hasKokkosMaterialPropertyByName (const std::string &name)
 
bool hasGenericMaterialProperty (const std::string &name)
 
bool hasGenericMaterialProperty (const std::string &name)
 
bool hasGenericMaterialPropertyByName (const std::string &name)
 
bool hasGenericMaterialPropertyByName (const std::string &name)
 
const Function & getFunction (const std::string &name) const
 
const Function & getFunctionByName (const FunctionName &name) const
 
bool hasFunction (const std::string &param_name) const
 
bool hasFunctionByName (const FunctionName &name) const
 
Moose::Kokkos::Function getKokkosFunction (const std::string &name) const
 
const T & getKokkosFunction (const std::string &name) const
 
Moose::Kokkos::Function getKokkosFunctionByName (const FunctionName &name) const
 
const T & getKokkosFunctionByName (const FunctionName &name) const
 
bool hasKokkosFunction (const std::string &param_name) const
 
bool hasKokkosFunctionByName (const FunctionName &name) const
 
bool isDefaultPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
bool hasPostprocessor (const std::string &param_name, const unsigned int index=0) const
 
bool hasPostprocessorByName (const PostprocessorName &name) const
 
std::size_t coupledPostprocessors (const std::string &param_name) const
 
const PostprocessorName & getPostprocessorName (const std::string &param_name, const unsigned int index=0) const
 
const VectorPostprocessorValue & getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValue & getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValue & getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValue & getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValue & getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValue & getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValue & getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValue & getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const ScatterVectorPostprocessorValue & getScatterVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValue & getScatterVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValue & getScatterVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValue & getScatterVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
bool hasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const
 
bool hasVectorPostprocessor (const std::string &param_name) const
 
bool hasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
bool hasVectorPostprocessorByName (const VectorPostprocessorName &name) const
 
const VectorPostprocessorName & getVectorPostprocessorName (const std::string &param_name) const
 
T & getSampler (const std::string &name)
 
Sampler & getSampler (const std::string &name)
 
T & getSamplerByName (const SamplerName &name)
 
Sampler & getSamplerByName (const SamplerName &name)
 
virtual void meshChanged ()
 
virtual void meshDisplaced ()
 
PerfGraph & perfGraph ()
 
const PostprocessorValue & getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValue & getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValue & getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValue & getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValue & getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValue & getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
virtual const PostprocessorValue & getPostprocessorValueByName (const PostprocessorName &name) const
 
virtual const PostprocessorValue & getPostprocessorValueByName (const PostprocessorName &name) const
 
const PostprocessorValue & getPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValue & getPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValue & getPostprocessorValueOlderByName (const PostprocessorName &name) const
 
const PostprocessorValue & getPostprocessorValueOlderByName (const PostprocessorName &name) const
 
bool isVectorPostprocessorDistributed (const std::string &param_name) const
 
bool isVectorPostprocessorDistributed (const std::string &param_name) const
 
bool isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const
 
bool isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const
 
const Distribution & getDistribution (const std::string &name) const
 
const T & getDistribution (const std::string &name) const
 
const Distribution & getDistribution (const std::string &name) const
 
const T & getDistribution (const std::string &name) const
 
const Distribution & getDistributionByName (const DistributionName &name) const
 
const T & getDistributionByName (const std::string &name) const
 
const Distribution & getDistributionByName (const DistributionName &name) const
 
const T & getDistributionByName (const std::string &name) const
 
const Parallel::Communicator & comm () const
 
processor_id_type n_processors () const
 
processor_id_type processor_id () const
 

Static Public Member Functions

static InputParameters validParams ()
 
static void callMooseError (MooseApp *const app, const InputParameters &params, std::string msg, const bool with_prefix, const hit::Node *node, const bool show_trace=true)
 
static void sort (typename std::vector< T > &vector)
 
static void sortDFS (typename std::vector< T > &vector)
 
static void cyclicDependencyError (CyclicDependencyException< T2 > &e, const std::string &header, NameFunc &&name_func)
 
static void cyclicDependencyError (CyclicDependencyException< T2 > &e, const std::string &header)
 

Public Attributes

const Real _f_tol
 Tolerance on yield function.
 
const Real _ic_tol
 Tolerance on internal constraint.
 
 usingCombinedWarningSolutionWarnings
 
const ConsoleStream _console
 

Static Public Attributes

static const std::string type_param
 
static const std::string name_param
 
static const std::string unique_name_param
 
static const std::string app_param
 
static const std::string moose_base_param
 
static const std::string kokkos_object_param
 
static constexpr PropertyValue::id_type default_property_id
 
static constexpr PropertyValue::id_type zero_property_id
 
static constexpr auto SYSTEM
 
static constexpr auto NAME
 

Protected Member Functions

Real yieldFunction (const RankTwoTensor &stress, Real intnl) const override
 The following functions are what you should override when building single-plasticity models.
 
RankTwoTensor dyieldFunction_dstress (const RankTwoTensor &stress, Real intnl) const override
 The derivative of yield function with respect to stress.
 
Real dyieldFunction_dintnl (const RankTwoTensor &stress, Real intnl) const override
 The derivative of yield function with respect to the internal parameter.
 
RankTwoTensor flowPotential (const RankTwoTensor &stress, Real intnl) const override
 The flow potential.
 
RankFourTensor dflowPotential_dstress (const RankTwoTensor &stress, Real intnl) const override
 The derivative of the flow potential with respect to stress.
 
RankTwoTensor dflowPotential_dintnl (const RankTwoTensor &stress, Real intnl) const override
 The derivative of the flow potential with respect to the internal parameter.
 
Real hardPotential (const RankTwoTensor &stress, Real intnl) const override
 The hardening potential.
 
virtual RankTwoTensor dhardPotential_dstress (const RankTwoTensor &stress, Real intnl) const override
 The derivative of the hardening potential with respect to stress.
 
virtual Real dhardPotential_dintnl (const RankTwoTensor &stress, Real intnl) const override
 The derivative of the hardening potential with respect to the internal parameter.
 
RankTwoTensor df_dsig (const RankTwoTensor &stress, Real intnl) const
 Derivative of the yield function with respect to stress.
 
virtual Real tensile_strength (const Real internal_param) const
 tensile strength as a function of residual value, rate, and internal_param
 
virtual Real dtensile_strength (const Real internal_param) const
 d(tensile strength)/d(internal_param) as a function of residual value, rate, and internal_param
 
virtual Real compressive_strength (const Real internal_param) const
 compressive strength as a function of residual value, rate, and internal_param
 
virtual Real dcompressive_strength (const Real internal_param) const
 d(compressive strength)/d(internal_param) as a function of residual value, rate, and internal_param
 
virtual void addPostprocessorDependencyHelper (const PostprocessorName &name) const override
 
virtual void addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &name) const override
 
virtual void addUserObjectDependencyHelper (const UserObjectBase &uo) const override
 
void addReporterDependencyHelper (const ReporterName &reporter_name) override
 
void flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const
 
InvalidSolutionID registerInvalidSolutionInternal (const std::string &message, const bool warning) const
 
const ReporterContextBase & getReporterContextBaseByName (const ReporterName &reporter_name) const
 
const ReporterName & getReporterName (const std::string &param_name) const
 
T & declareRestartableData (const std::string &data_name, Args &&... args)
 
ManagedValue< T > declareManagedRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
const T & getRestartableData (const std::string &data_name) const
 
T & declareRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
T & declareRecoverableData (const std::string &data_name, Args &&... args)
 
T & declareRestartableDataWithObjectName (const std::string &data_name, const std::string &object_name, Args &&... args)
 
T & declareRestartableDataWithObjectNameWithContext (const std::string &data_name, const std::string &object_name, void *context, Args &&... args)
 
std::string restartableName (const std::string &data_name) const
 
const T & getMeshProperty (const std::string &data_name, const std::string &prefix)
 
const T & getMeshProperty (const std::string &data_name)
 
bool hasMeshProperty (const std::string &data_name, const std::string &prefix) const
 
bool hasMeshProperty (const std::string &data_name, const std::string &prefix) const
 
bool hasMeshProperty (const std::string &data_name) const
 
bool hasMeshProperty (const std::string &data_name) const
 
std::string meshPropertyName (const std::string &data_name) const
 
PerfID registerTimedSection (const std::string &section_name, const unsigned int level) const
 
PerfID registerTimedSection (const std::string &section_name, const unsigned int level, const std::string &live_message, const bool print_dots=true) const
 
std::string timedSectionName (const std::string &section_name) const
 
bool isCoupledScalar (const std::string &var_name, unsigned int i=0) const
 
unsigned int coupledScalarComponents (const std::string &var_name) const
 
unsigned int coupledScalar (const std::string &var_name, unsigned int comp=0) const
 
libMesh::Order coupledScalarOrder (const std::string &var_name, unsigned int comp=0) const
 
const VariableValue & coupledScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValue & adCoupledScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< is_ad > & coupledGenericScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< false > & coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const
 
const GenericVariableValue< true > & coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const
 
const VariableValue & coupledVectorTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValue & coupledMatrixTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValue & coupledScalarValueOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValue & coupledScalarValueOlder (const std::string &var_name, unsigned int comp=0) const
 
const VariableValue & coupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValue & adCoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValue & coupledScalarDotDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValue & coupledScalarDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValue & coupledScalarDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValue & coupledScalarDotDu (const std::string &var_name, unsigned int comp=0) const
 
const VariableValue & coupledScalarDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
const MooseVariableScalar * getScalarVar (const std::string &var_name, unsigned int comp) const
 
virtual void checkMaterialProperty (const std::string &name, const unsigned int state)
 
virtual void getKokkosMaterialPropertyHook (const std::string &, const unsigned int)
 
void markMatPropRequested (const std::string &)
 
MaterialPropertyName getMaterialPropertyName (const std::string &name) const
 
void checkExecutionStage ()
 
const T & getReporterValue (const std::string &param_name, const std::size_t time_index=0)
 
const T & getReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)
 
const T & getReporterValue (const std::string &param_name, const std::size_t time_index=0)
 
const T & getReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)
 
const T & getReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)
 
const T & getReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)
 
const T & getReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)
 
const T & getReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 
const GenericMaterialProperty< T, is_ad > * defaultGenericMaterialProperty (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > * defaultGenericMaterialProperty (const std::string &name)
 
const MaterialProperty< T > * defaultMaterialProperty (const std::string &name)
 
const MaterialProperty< T > * defaultMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > * defaultADMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > * defaultADMaterialProperty (const std::string &name)
 

Static Protected Member Functions

static std::string meshPropertyName (const std::string &data_name, const std::string &prefix)
 

Protected Attributes

const unsigned _max_iters
 max iters for custom return map loop
 
const bool _use_custom_returnMap
 Whether to use the custom return-map algorithm.
 
const bool _use_custom_cto
 Whether to use the custom consistent tangent operator algorithm.
 
const SolidMechanicsHardeningModel & _strength
 the tensile strength
 
const SolidMechanicsHardeningModel & _c_strength
 the compressive strength
 
const Moose::CoordinateSystemType & _coord_sys
 
const THREAD_ID _tid
 
SubProblem & _subproblem
 
FEProblemBase & _fe_problem
 
SystemBase & _sys
 
Assembly & _assembly
 
const bool _duplicate_initial_execution
 
std::set< std::string > _depend_uo
 
const bool & _enabled
 
MooseApp & _app
 
Factory & _factory
 
ActionFactory & _action_factory
 
const std::string & _type
 
const std::string & _name
 
const InputParameters & _pars
 
const ExecFlagEnum & _execute_enum
 
const ExecFlagType & _current_execute_flag
 
MooseApp & _restartable_app
 
const std::string _restartable_system_name
 
const THREAD_ID _restartable_tid
 
const bool _restartable_read_only
 
FEProblemBase & _mci_feproblem
 
FEProblemBase & _mdi_feproblem
 
MooseApp & _pg_moose_app
 
const std::string _prefix
 
FEProblemBase & _sc_fe_problem
 
const THREAD_ID _sc_tid
 
const Real & _real_zero
 
const VariableValue & _scalar_zero
 
const Point & _point_zero
 
const InputParameters & _mi_params
 
const std::string _mi_name
 
const MooseObjectName _mi_moose_object_name
 
FEProblemBase & _mi_feproblem
 
SubProblem & _mi_subproblem
 
const THREAD_ID _mi_tid
 
const bool _is_kokkos_object
 
const Moose::MaterialDataType _material_data_type
 
MaterialData & _material_data
 
bool _stateful_allowed
 
bool _get_material_property_called
 
std::vector< std::unique_ptr< PropertyValue > > _default_properties
 
std::unordered_set< unsigned int > _material_property_dependencies
 
const MaterialPropertyName _get_suffix
 
const bool _use_interpolated_state
 
const InputParameters & _ti_params
 
FEProblemBase & _ti_feproblem
 
bool _is_implicit
 
Real & _t
 
const Real & _t_old
 
int & _t_step
 
Real & _dt
 
Real & _dt_old
 
bool _is_transient
 
const Parallel::Communicator & _communicator
 

Static Protected Attributes

static const std::string _interpolated_old
 
static const std::string _interpolated_older
 

Private Member Functions

const Moose::FunctionBase * getKokkosFunctionByNameHelper (const FunctionName &name) const
 
const UserObjectBase & getUserObjectFromFEProblem (const UserObjectName &object_name, const THREAD_ID tid=0) const
 
const T & castUserObject (const UserObjectBase &uo_base, const std::string &param_name="") const
 
void mooseObjectError (const std::string &param_name, std::stringstream &oss) const
 
const std::string & userObjectType (const UserObjectBase &uo) const
 
const std::string & userObjectName (const UserObjectBase &uo) const
 
const PostprocessorName & getPostprocessorNameInternal (const std::string &param_name, const unsigned int index, const bool allow_default_value=true) const
 
bool isDefaultPostprocessorValueByName (const PostprocessorName &name) const
 
PostprocessorValue getDefaultPostprocessorValueByName (const PostprocessorName &name) const
 
void checkParam (const std::string &param_name, const unsigned int index=std::numeric_limits< unsigned int >::max()) const
 
bool postprocessorsAdded () const
 
const VectorPostprocessorValue & getVectorPostprocessorByNameHelper (const VectorPostprocessorName &name, const std::string &vector_name, bool broadcast, std::size_t t_index) const
 
const VectorPostprocessorContext< VectorPostprocessorValue > & getVectorPostprocessorContextByNameHelper (const VectorPostprocessorName &name, const std::string &vector_name) const
 
bool vectorPostprocessorsAdded () const
 
bool reportersAdded () const
 
void possiblyCheckHasReporter (const ReporterName &reporter_name, const std::string &param_name="") const
 
RestartableDataValue & registerRestartableDataOnApp (std::unique_ptr< RestartableDataValue > data, THREAD_ID tid) const
 
void registerRestartableNameWithFilterOnApp (const std::string &name, Moose::RESTARTABLE_FILTER filter)
 
RestartableData< T > & declareRestartableDataHelper (const std::string &data_name, void *context, Args &&... args) const
 
virtual std::string meshPropertyPrefix (const std::string &data_name) const
 
const RestartableDataValue & getMeshPropertyInternal (const std::string &data_name, const std::string &prefix) const
 
void mooseErrorInternal (Args &&... args) const
 
const VariableValue * getDefaultValue (const std::string &var_name) const
 
const ADVariableValue * getADDefaultValue (const std::string &var_name) const
 
void checkVar (const std::string &var_name) const
 
void validateExecutionerType (const std::string &name, const std::string &fn_name) const
 
Moose::MaterialDataType getMaterialDataType (const std::set< BoundaryID > &boundary_ids) const
 
unsigned int getMaxQps () const
 
void addConsumedPropertyName (const MooseObjectName &obj_name, const std::string &prop_name)
 
const PostprocessorValue & getPostprocessorValueInternal (const std::string &param_name, unsigned int index, std::size_t t_index) const
 
const PostprocessorValue & getPostprocessorValueInternal (const std::string &param_name, unsigned int index, std::size_t t_index) const
 
const PostprocessorValue & getPostprocessorValueByNameInternal (const PostprocessorName &name, std::size_t t_index) const
 
const PostprocessorValue & getPostprocessorValueByNameInternal (const PostprocessorName &name, std::size_t t_index) const
 
void possiblyCheckHasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const
 
void possiblyCheckHasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const
 
void possiblyCheckHasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
void possiblyCheckHasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 

Static Private Member Functions

static const hit::Node * getHitNode (const InputParameters &params)
 
static std::string messagePrefix (const InputParameters &params, const bool hit_prefix)
 

Private Attributes

UserObject * _primary_thread_copy
 
std::set< std::string > _supplied_uo
 
const ParallelParamObject & _parent
 
const MooseBase & _si_moose_base
 
const FEProblemBase * _si_problem
 
ExecFlagEnum _empty_execute_enum
 
const MooseObject & _fni_object
 
const InputParameters & _fni_params
 
FEProblemBase & _fni_feproblem
 
const THREAD_ID _fni_tid
 
const MooseObject & _uoi_moose_object
 
const FEProblemBase & _uoi_feproblem
 
const THREAD_ID _uoi_tid
 
const MooseObject & _ppi_moose_object
 
const InputParameters & _ppi_params
 
const FEProblemBase & _ppi_feproblem
 
std::map< PostprocessorName, std::unique_ptr< PostprocessorValue > > _default_values
 
const bool _broadcast_by_default
 
const MooseObject & _vpi_moose_object
 
const FEProblemBase & _vpi_feproblem
 
const THREAD_ID _vpi_tid
 
const InputParameters & _ri_params
 
FEProblemBase & _ri_fe_problem_base
 
const ReporterData & _ri_reporter_data
 
const MooseObject & _ri_moose_object
 
const InputParameters & _dni_params
 
FEProblemBase & _dni_feproblem
 
const MooseObject *const _dni_moose_object_ptr
 
const InputParameters & _si_params
 
FEProblemBase & _si_feproblem
 
THREAD_ID _si_tid
 
const RestartableDataMapName _metaname
 
std::string _restartable_name
 
MooseApp & _meta_data_app
 
const MooseObject *const _meta_data_object
 
const InputParameters & _sc_parameters
 
const std::string & _sc_name
 
const bool _sc_is_implicit
 
std::unordered_map< std::string, std::vector< MooseVariableScalar * > > _coupled_scalar_vars
 
std::unordered_map< std::string, std::unique_ptr< VariableValue > > _default_value
 
std::unordered_map< std::string, std::unique_ptr< ADVariableValue > > _dual_default_value
 
std::vector< MooseVariableScalar * > _coupled_moose_scalar_vars
 
std::unordered_map< std::string, std::vector< MooseVariableFieldBase * > > _sc_coupled_vars
 
std::set< TagID > _sc_coupleable_vector_tags
 
std::set< TagID > _sc_coupleable_matrix_tags
 
const MooseObject & _mi_moose_object
 
const bool _mi_boundary_restricted
 
const std::set< SubdomainID > & _mi_block_ids
 
const std::set< BoundaryID > & _mi_boundary_ids
 
std::vector< std::unique_ptr< OptionalMaterialPropertyProxyBase< MaterialPropertyInterface > > > _optional_property_proxies
 
const std::string _ti_name
 

Detailed Description

Rate-independent associative mean-cap tensile AND compressive failure with hardening/softening of the tensile and compressive strength.

The key point here is that the internal parameter is equal to the volumetric plastic strain. This means that upon tensile failure, the compressive strength can soften (using a SolidMechanicsHardening object) which physically means that a subsequent compressive stress can easily squash the material

Definition at line 24 of file SolidMechanicsPlasticMeanCapTC.h.

Constructor & Destructor Documentation

◆ SolidMechanicsPlasticMeanCapTC()

SolidMechanicsPlasticMeanCapTC::SolidMechanicsPlasticMeanCapTC ( const InputParameters &  parameters)

Definition at line 51 of file SolidMechanicsPlasticMeanCapTC.C.

52 : SolidMechanicsPlasticModel(parameters),
53 _max_iters(getParam<unsigned>("max_iterations")),
54 _use_custom_returnMap(getParam<bool>("use_custom_returnMap")),
55 _use_custom_cto(getParam<bool>("use_custom_cto")),
56 _strength(getUserObject<SolidMechanicsHardeningModel>("tensile_strength")),
57 _c_strength(getUserObject<SolidMechanicsHardeningModel>("compressive_strength"))
58{
59 // cannot check the following for all values of the internal parameter, but this will catch most
60 // errors
61 if (_strength.value(0) <= _c_strength.value(0))
62 mooseError("MeanCapTC: tensile strength (which is usually positive) must not be less than "
63 "compressive strength (which is usually negative)");
64}
GeneralUserObject::parameters
const InputParameters & parameters() const
GeneralUserObject::mooseError
void mooseError(Args &&... args) const
SolidMechanicsHardeningModel::value
virtual Real value(Real intnl) const
Definition SolidMechanicsHardeningModel.C:44
SolidMechanicsPlasticMeanCapTC::_use_custom_returnMap
const bool _use_custom_returnMap
Whether to use the custom return-map algorithm.
Definition SolidMechanicsPlasticMeanCapTC.h:68
SolidMechanicsPlasticMeanCapTC::_c_strength
const SolidMechanicsHardeningModel & _c_strength
the compressive strength
Definition SolidMechanicsPlasticMeanCapTC.h:77
SolidMechanicsPlasticMeanCapTC::_use_custom_cto
const bool _use_custom_cto
Whether to use the custom consistent tangent operator algorithm.
Definition SolidMechanicsPlasticMeanCapTC.h:71
SolidMechanicsPlasticMeanCapTC::_strength
const SolidMechanicsHardeningModel & _strength
the tensile strength
Definition SolidMechanicsPlasticMeanCapTC.h:74
SolidMechanicsPlasticMeanCapTC::_max_iters
const unsigned _max_iters
max iters for custom return map loop
Definition SolidMechanicsPlasticMeanCapTC.h:65
SolidMechanicsPlasticModel
Plastic Model base class The virtual functions written below must be over-ridden in derived classes t...
Definition SolidMechanicsPlasticModel.h:38

Member Function Documentation

◆ activeConstraints()

void SolidMechanicsPlasticMeanCapTC::activeConstraints ( const std::vector< Real > &  f,
const RankTwoTensor &  stress,
Real  intnl,
const RankFourTensor &  Eijkl,
std::vector< bool > &  act,
RankTwoTensor &  returned_stress 
) const
overridevirtual

The active yield surfaces, given a vector of yield functions.

This is used by FiniteStrainMultiPlasticity to determine the initial set of active constraints at the trial (stress, intnl) configuration. It is up to you (the coder) to determine how accurate you want the returned_stress to be. Currently it is only used by FiniteStrainMultiPlasticity to estimate a good starting value for the Newton-Rahson procedure, so currently it may not need to be super perfect.

Parameters
fvalues of the yield functions
stressstress tensor
intnlinternal parameter
Eijklelasticity tensor (stress = Eijkl*strain)
[out]actact[i] = true if the i_th yield function is active
[out]returned_stressApproximate value of the returned stress

Reimplemented from SolidMechanicsPlasticModel.

Definition at line 249 of file SolidMechanicsPlasticMeanCapTC.C.

255{
256 act.assign(1, false);
257
258 if (f[0] <= _f_tol)
259 {
260 returned_stress = stress;
261 return;
262 }
263
264 const Real tr = stress.trace();
265 const Real t_str = tensile_strength(intnl);
266 Real str;
267 Real dirn;
268 if (tr >= t_str)
269 {
270 str = t_str;
271 dirn = 1.0;
272 }
273 else
274 {
275 str = compressive_strength(intnl);
276 dirn = -1.0;
277 }
278
279 RankTwoTensor n; // flow direction
280 for (unsigned i = 0; i < 3; ++i)
281 for (unsigned j = 0; j < 3; ++j)
282 for (unsigned k = 0; k < 3; ++k)
283 n(i, j) += dirn * Eijkl(i, j, k, k);
284
285 // returned_stress = stress - gamma*n
286 // and taking the trace of this and using
287 // Tr(returned_stress) = str, gives
288 // gamma = (Tr(stress) - str)/Tr(n)
289 Real gamma = (stress.trace() - str) / n.trace();
290
291 for (unsigned i = 0; i < 3; ++i)
292 for (unsigned j = 0; j < 3; ++j)
293 returned_stress(i, j) = stress(i, j) - gamma * n(i, j);
294
295 act[0] = true;
296}
f
Real f(Real x)
Test function for Brents method.
Definition BrentsMethodTest.C:21
SolidMechanicsPlasticMeanCapTC::compressive_strength
virtual Real compressive_strength(const Real internal_param) const
compressive strength as a function of residual value, rate, and internal_param
Definition SolidMechanicsPlasticMeanCapTC.C:237
SolidMechanicsPlasticMeanCapTC::tensile_strength
virtual Real tensile_strength(const Real internal_param) const
tensile strength as a function of residual value, rate, and internal_param
Definition SolidMechanicsPlasticMeanCapTC.C:225
SolidMechanicsPlasticModel::_f_tol
const Real _f_tol
Tolerance on yield function.
Definition SolidMechanicsPlasticModel.h:170
libMesh::Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ compressive_strength()

Real SolidMechanicsPlasticMeanCapTC::compressive_strength ( const Real  internal_param) const
protectedvirtual

compressive strength as a function of residual value, rate, and internal_param

Definition at line 237 of file SolidMechanicsPlasticMeanCapTC.C.

238{
239 return _c_strength.value(internal_param);
240}

Referenced by activeConstraints(), df_dsig(), dflowPotential_dintnl(), dflowPotential_dstress(), dhardPotential_dintnl(), dhardPotential_dstress(), dyieldFunction_dintnl(), hardPotential(), returnMap(), and yieldFunction().

◆ consistentTangentOperator()

RankFourTensor SolidMechanicsPlasticMeanCapTC::consistentTangentOperator ( const RankTwoTensor &  trial_stress,
Real  intnl_old,
const RankTwoTensor &  stress,
Real  intnl,
const RankFourTensor &  E_ijkl,
const std::vector< Real > &  cumulative_pm 
) const
overridevirtual

Calculates a custom consistent tangent operator.

You may choose to over-ride this in your derived SolidMechanicsPlasticXXXX class.

(Note, if you over-ride returnMap, you will probably want to override consistentTangentOpertor too, otherwise it will default to E_ijkl.)

Parameters
stress_oldtrial stress before returning
intnl_oldinternal parameter before returning
stresscurrent returned stress state
intnlinternal parameter
E_ijklelasticity tensor
cumulative_pmthe cumulative plastic multipliers
Returns
the consistent tangent operator: E_ijkl if not over-ridden

Reimplemented from SolidMechanicsPlasticModel.

Definition at line 401 of file SolidMechanicsPlasticMeanCapTC.C.

408{
409 if (!_use_custom_cto)
410 return SolidMechanicsPlasticModel::consistentTangentOperator(
411 trial_stress, intnl_old, stress, intnl, E_ijkl, cumulative_pm);
412
413 Real df_dq;
414 Real alpha;
415 if (trial_stress.trace() >= tensile_strength(intnl_old))
416 {
417 df_dq = -dtensile_strength(intnl);
418 alpha = 1.0;
419 }
420 else
421 {
422 df_dq = dcompressive_strength(intnl);
423 alpha = -1.0;
424 }
425
426 RankTwoTensor elas;
427 for (unsigned int i = 0; i < 3; ++i)
428 for (unsigned int j = 0; j < 3; ++j)
429 for (unsigned int k = 0; k < 3; ++k)
430 elas(i, j) += E_ijkl(i, j, k, k);
431
432 const Real hw = -df_dq + alpha * elas.trace();
433
434 return E_ijkl - alpha / hw * elas.outerProduct(elas);
435}
RankTwoTensorTempl::outerProduct
RankFourTensorTempl< T > outerProduct(const RankTwoTensorTempl< T > &b) const
SolidMechanicsPlasticMeanCapTC::dtensile_strength
virtual Real dtensile_strength(const Real internal_param) const
d(tensile strength)/d(internal_param) as a function of residual value, rate, and internal_param
Definition SolidMechanicsPlasticMeanCapTC.C:231
SolidMechanicsPlasticMeanCapTC::dcompressive_strength
virtual Real dcompressive_strength(const Real internal_param) const
d(compressive strength)/d(internal_param) as a function of residual value, rate, and internal_param
Definition SolidMechanicsPlasticMeanCapTC.C:243
SolidMechanicsPlasticModel::consistentTangentOperator
virtual RankFourTensor consistentTangentOperator(const RankTwoTensor &trial_stress, Real intnl_old, const RankTwoTensor &stress, Real intnl, const RankFourTensor &E_ijkl, const std::vector< Real > &cumulative_pm) const
Calculates a custom consistent tangent operator.
Definition SolidMechanicsPlasticModel.C:252

◆ dcompressive_strength()

Real SolidMechanicsPlasticMeanCapTC::dcompressive_strength ( const Real  internal_param) const
protectedvirtual

d(compressive strength)/d(internal_param) as a function of residual value, rate, and internal_param

Definition at line 243 of file SolidMechanicsPlasticMeanCapTC.C.

244{
245 return _c_strength.derivative(internal_param);
246}
SolidMechanicsHardeningModel::derivative
virtual Real derivative(Real intnl) const
Definition SolidMechanicsHardeningModel.C:50

Referenced by consistentTangentOperator(), dflowPotential_dintnl(), dhardPotential_dintnl(), dyieldFunction_dintnl(), and returnMap().

◆ df_dsig()

RankTwoTensor SolidMechanicsPlasticMeanCapTC::df_dsig ( const RankTwoTensor &  stress,
Real  intnl 
) const
protected

Derivative of the yield function with respect to stress.

This is also the flow potential.

Parameters
stressthe stress at which to calculate the hardening potentials
intnlinternal parameter
Returns
the derivative

Definition at line 112 of file SolidMechanicsPlasticMeanCapTC.C.

113{
114 const Real tr = stress.trace();
115 const Real t_str = tensile_strength(intnl);
116 if (tr >= t_str)
117 return stress.dtrace();
118
119 const Real c_str = compressive_strength(intnl);
120 if (tr <= c_str)
121 return -stress.dtrace();
122
123 return -std::cos(libMesh::pi * (tr - c_str) / (t_str - c_str)) * stress.dtrace();
124}
RankTwoTensorTempl::dtrace
RankTwoTensorTempl< T > dtrace() const
libMesh::pi
const Real pi

Referenced by dyieldFunction_dstress(), and flowPotential().

◆ dflowPotential_dintnl()

RankTwoTensor SolidMechanicsPlasticMeanCapTC::dflowPotential_dintnl ( const RankTwoTensor &  stress,
Real  intnl 
) const
overrideprotectedvirtual

The derivative of the flow potential with respect to the internal parameter.

Parameters
stressthe stress at which to calculate the flow potential
intnlinternal parameter
Returns
dr_dintnl(i, j) = dr(i, j)/dintnl

Reimplemented from SolidMechanicsPlasticModel.

Definition at line 150 of file SolidMechanicsPlasticMeanCapTC.C.

152{
153 const Real tr = stress.trace();
154 const Real t_str = tensile_strength(intnl);
155 if (tr >= t_str)
156 return RankTwoTensor();
157
158 const Real c_str = compressive_strength(intnl);
159 if (tr <= c_str)
160 return RankTwoTensor();
161
162 const Real dt = dtensile_strength(intnl);
163 const Real dc = dcompressive_strength(intnl);
164 return std::sin(libMesh::pi * (tr - c_str) / (t_str - c_str)) * stress.dtrace() * libMesh::pi /
165 Utility::pow<2>(t_str - c_str) * ((tr - t_str) * dc - (tr - c_str) * dt);
166}

◆ dflowPotential_dintnlV()

void SolidMechanicsPlasticModel::dflowPotential_dintnlV ( const RankTwoTensor &  stress,
Real  intnl,
std::vector< RankTwoTensor > &  dr_dintnl 
) const
virtualinherited

The derivative of the flow potential with respect to the internal parameter.

Parameters
stressthe stress at which to calculate the flow potential
intnlinternal parameter
[out]dr_dintnldr_dintnl[alpha](i, j) = dr[alpha](i, j)/dintnl

Reimplemented in SolidMechanicsPlasticMohrCoulombMulti, and SolidMechanicsPlasticTensileMulti.

Definition at line 137 of file SolidMechanicsPlasticModel.C.

140{
141 return dr_dintnl.assign(1, dflowPotential_dintnl(stress, intnl));
142}
SolidMechanicsPlasticModel::dflowPotential_dintnl
virtual RankTwoTensor dflowPotential_dintnl(const RankTwoTensor &stress, Real intnl) const
The derivative of the flow potential with respect to the internal parameter.
Definition SolidMechanicsPlasticModel.C:131

◆ dflowPotential_dstress()

RankFourTensor SolidMechanicsPlasticMeanCapTC::dflowPotential_dstress ( const RankTwoTensor &  stress,
Real  intnl 
) const
overrideprotectedvirtual

The derivative of the flow potential with respect to stress.

Parameters
stressthe stress at which to calculate the flow potential
intnlinternal parameter
Returns
dr_dstress(i, j, k, l) = dr(i, j)/dstress(k, l)

Reimplemented from SolidMechanicsPlasticModel.

Definition at line 133 of file SolidMechanicsPlasticMeanCapTC.C.

135{
136 const Real tr = stress.trace();
137 const Real t_str = tensile_strength(intnl);
138 if (tr >= t_str)
139 return RankFourTensor();
140
141 const Real c_str = compressive_strength(intnl);
142 if (tr <= c_str)
143 return RankFourTensor();
144
145 return libMesh::pi / (t_str - c_str) * std::sin(libMesh::pi * (tr - c_str) / (t_str - c_str)) *
146 stress.dtrace().outerProduct(stress.dtrace());
147}

◆ dflowPotential_dstressV()

void SolidMechanicsPlasticModel::dflowPotential_dstressV ( const RankTwoTensor &  stress,
Real  intnl,
std::vector< RankFourTensor > &  dr_dstress 
) const
virtualinherited

The derivative of the flow potential with respect to stress.

Parameters
stressthe stress at which to calculate the flow potential
intnlinternal parameter
[out]dr_dstressdr_dstress[alpha](i, j, k, l) = dr[alpha](i, j)/dstress(k, l)

Reimplemented in SolidMechanicsPlasticMohrCoulombMulti, and SolidMechanicsPlasticTensileMulti.

Definition at line 123 of file SolidMechanicsPlasticModel.C.

126{
127 return dr_dstress.assign(1, dflowPotential_dstress(stress, intnl));
128}
SolidMechanicsPlasticModel::dflowPotential_dstress
virtual RankFourTensor dflowPotential_dstress(const RankTwoTensor &stress, Real intnl) const
The derivative of the flow potential with respect to stress.
Definition SolidMechanicsPlasticModel.C:117

◆ dhardPotential_dintnl()

Real SolidMechanicsPlasticMeanCapTC::dhardPotential_dintnl ( const RankTwoTensor &  stress,
Real  intnl 
) const
overrideprotectedvirtual

The derivative of the hardening potential with respect to the internal parameter.

Parameters
stressthe stress at which to calculate the hardening potentials
intnlinternal parameter
Returns
the derivative

Reimplemented from SolidMechanicsPlasticModel.

Definition at line 206 of file SolidMechanicsPlasticMeanCapTC.C.

208{
209 const Real tr = stress.trace();
210 const Real t_str = tensile_strength(intnl);
211 if (tr >= t_str)
212 return 0.0;
213
214 const Real c_str = compressive_strength(intnl);
215 if (tr <= c_str)
216 return 0.0;
217
218 const Real dt = dtensile_strength(intnl);
219 const Real dc = dcompressive_strength(intnl);
220 return -std::sin(libMesh::pi * (tr - c_str) / (t_str - c_str)) * libMesh::pi /
221 Utility::pow<2>(t_str - c_str) * ((tr - t_str) * dc - (tr - c_str) * dt);
222}

◆ dhardPotential_dintnlV()

void SolidMechanicsPlasticModel::dhardPotential_dintnlV ( const RankTwoTensor &  stress,
Real  intnl,
std::vector< Real > &  dh_dintnl 
) const
virtualinherited

The derivative of the hardening potential with respect to the internal parameter.

Parameters
stressthe stress at which to calculate the hardening potentials
intnlinternal parameter
[out]dh_dintnldh_dintnl[alpha] = dh[alpha]/dintnl

Definition at line 178 of file SolidMechanicsPlasticModel.C.

181{
182 dh_dintnl.resize(numberSurfaces(), dhardPotential_dintnl(stress, intnl));
183}
SolidMechanicsPlasticModel::numberSurfaces
virtual unsigned int numberSurfaces() const
The number of yield surfaces for this plasticity model.
Definition SolidMechanicsPlasticModel.C:55
SolidMechanicsPlasticModel::dhardPotential_dintnl
virtual Real dhardPotential_dintnl(const RankTwoTensor &stress, Real intnl) const
The derivative of the hardening potential with respect to the internal parameter.
Definition SolidMechanicsPlasticModel.C:172

◆ dhardPotential_dstress()

RankTwoTensor SolidMechanicsPlasticMeanCapTC::dhardPotential_dstress ( const RankTwoTensor &  stress,
Real  intnl 
) const
overrideprotectedvirtual

The derivative of the hardening potential with respect to stress.

Parameters
stressthe stress at which to calculate the hardening potentials
intnlinternal parameter
Returns
dh_dstress(i, j) = dh/dstress(i, j)

Reimplemented from SolidMechanicsPlasticModel.

Definition at line 189 of file SolidMechanicsPlasticMeanCapTC.C.

191{
192 const Real tr = stress.trace();
193 const Real t_str = tensile_strength(intnl);
194 if (tr >= t_str)
195 return RankTwoTensor();
196
197 const Real c_str = compressive_strength(intnl);
198 if (tr <= c_str)
199 return RankTwoTensor();
200
201 return -std::sin(libMesh::pi * (tr - c_str) / (t_str - c_str)) * libMesh::pi / (t_str - c_str) *
202 stress.dtrace();
203}

◆ dhardPotential_dstressV()

void SolidMechanicsPlasticModel::dhardPotential_dstressV ( const RankTwoTensor &  stress,
Real  intnl,
std::vector< RankTwoTensor > &  dh_dstress 
) const
virtualinherited

The derivative of the hardening potential with respect to stress.

Parameters
stressthe stress at which to calculate the hardening potentials
intnlinternal parameter
[out]dh_dstressdh_dstress[alpha](i, j) = dh[alpha]/dstress(i, j)

Definition at line 164 of file SolidMechanicsPlasticModel.C.

167{
168 dh_dstress.assign(numberSurfaces(), dhardPotential_dstress(stress, intnl));
169}
SolidMechanicsPlasticModel::dhardPotential_dstress
virtual RankTwoTensor dhardPotential_dstress(const RankTwoTensor &stress, Real intnl) const
The derivative of the hardening potential with respect to stress.
Definition SolidMechanicsPlasticModel.C:158

◆ dtensile_strength()

Real SolidMechanicsPlasticMeanCapTC::dtensile_strength ( const Real  internal_param) const
protectedvirtual

d(tensile strength)/d(internal_param) as a function of residual value, rate, and internal_param

Definition at line 231 of file SolidMechanicsPlasticMeanCapTC.C.

232{
233 return _strength.derivative(internal_param);
234}

Referenced by consistentTangentOperator(), dflowPotential_dintnl(), dhardPotential_dintnl(), dyieldFunction_dintnl(), and returnMap().

◆ dyieldFunction_dintnl()

Real SolidMechanicsPlasticMeanCapTC::dyieldFunction_dintnl ( const RankTwoTensor &  stress,
Real  intnl 
) const
overrideprotectedvirtual

The derivative of yield function with respect to the internal parameter.

Parameters
stressthe stress at which to calculate the yield function
intnlinternal parameter
Returns
the derivative

Reimplemented from SolidMechanicsPlasticModel.

Definition at line 92 of file SolidMechanicsPlasticMeanCapTC.C.

94{
95 const Real tr = stress.trace();
96 const Real t_str = tensile_strength(intnl);
97 if (tr >= t_str)
98 return -dtensile_strength(intnl);
99
100 const Real c_str = compressive_strength(intnl);
101 if (tr <= c_str)
102 return dcompressive_strength(intnl);
103
104 const Real dt = dtensile_strength(intnl);
105 const Real dc = dcompressive_strength(intnl);
106 return (dc - dt) / libMesh::pi * std::sin(libMesh::pi * (tr - c_str) / (t_str - c_str)) +
107 1.0 / (t_str - c_str) * std::cos(libMesh::pi * (tr - c_str) / (t_str - c_str)) *
108 ((tr - c_str) * dt - (tr - t_str) * dc);
109}

◆ dyieldFunction_dintnlV()

void SolidMechanicsPlasticModel::dyieldFunction_dintnlV ( const RankTwoTensor &  stress,
Real  intnl,
std::vector< Real > &  df_dintnl 
) const
virtualinherited

The derivative of yield functions with respect to the internal parameter.

Parameters
stressthe stress at which to calculate the yield function
intnlinternal parameter
[out]df_dintnldf_dintnl[alpha] = df[alpha]/dintnl

Reimplemented in SolidMechanicsPlasticMohrCoulombMulti, and SolidMechanicsPlasticTensileMulti.

Definition at line 96 of file SolidMechanicsPlasticModel.C.

99{
100 return df_dintnl.assign(1, dyieldFunction_dintnl(stress, intnl));
101}
SolidMechanicsPlasticModel::dyieldFunction_dintnl
virtual Real dyieldFunction_dintnl(const RankTwoTensor &stress, Real intnl) const
The derivative of yield function with respect to the internal parameter.
Definition SolidMechanicsPlasticModel.C:90

◆ dyieldFunction_dstress()

RankTwoTensor SolidMechanicsPlasticMeanCapTC::dyieldFunction_dstress ( const RankTwoTensor &  stress,
Real  intnl 
) const
overrideprotectedvirtual

The derivative of yield function with respect to stress.

Parameters
stressthe stress at which to calculate the yield function
intnlinternal parameter
Returns
df_dstress(i, j) = dyieldFunction/dstress(i, j)

Reimplemented from SolidMechanicsPlasticModel.

Definition at line 85 of file SolidMechanicsPlasticMeanCapTC.C.

87{
88 return df_dsig(stress, intnl);
89}
SolidMechanicsPlasticMeanCapTC::df_dsig
RankTwoTensor df_dsig(const RankTwoTensor &stress, Real intnl) const
Derivative of the yield function with respect to stress.
Definition SolidMechanicsPlasticMeanCapTC.C:112

◆ dyieldFunction_dstressV()

void SolidMechanicsPlasticModel::dyieldFunction_dstressV ( const RankTwoTensor &  stress,
Real  intnl,
std::vector< RankTwoTensor > &  df_dstress 
) const
virtualinherited

The derivative of yield functions with respect to stress.

Parameters
stressthe stress at which to calculate the yield function
intnlinternal parameter
[out]df_dstressdf_dstress[alpha](i, j) = dyieldFunction[alpha]/dstress(i, j)

Reimplemented in SolidMechanicsPlasticMohrCoulombMulti, and SolidMechanicsPlasticTensileMulti.

Definition at line 82 of file SolidMechanicsPlasticModel.C.

85{
86 df_dstress.assign(1, dyieldFunction_dstress(stress, intnl));
87}
SolidMechanicsPlasticModel::dyieldFunction_dstress
virtual RankTwoTensor dyieldFunction_dstress(const RankTwoTensor &stress, Real intnl) const
The derivative of yield function with respect to stress.
Definition SolidMechanicsPlasticModel.C:75

◆ execute()

void SolidMechanicsPlasticModel::execute ( )
virtualinherited

Implements GeneralUserObject.

Definition at line 45 of file SolidMechanicsPlasticModel.C.

46{
47}

◆ finalize()

void SolidMechanicsPlasticModel::finalize ( )
virtualinherited

Implements GeneralUserObject.

Definition at line 50 of file SolidMechanicsPlasticModel.C.

51{
52}

◆ flowPotential()

RankTwoTensor SolidMechanicsPlasticMeanCapTC::flowPotential ( const RankTwoTensor &  stress,
Real  intnl 
) const
overrideprotectedvirtual

The flow potential.

Parameters
stressthe stress at which to calculate the flow potential
intnlinternal parameter
Returns
the flow potential

Reimplemented from SolidMechanicsPlasticModel.

Definition at line 127 of file SolidMechanicsPlasticMeanCapTC.C.

128{
129 return df_dsig(stress, intnl);
130}

◆ flowPotentialV()

void SolidMechanicsPlasticModel::flowPotentialV ( const RankTwoTensor &  stress,
Real  intnl,
std::vector< RankTwoTensor > &  r 
) const
virtualinherited

The flow potentials.

Parameters
stressthe stress at which to calculate the flow potential
intnlinternal parameter
[out]rr[alpha] is the flow potential for the "alpha" yield function

Reimplemented in SolidMechanicsPlasticMohrCoulombMulti, and SolidMechanicsPlasticTensileMulti.

Definition at line 109 of file SolidMechanicsPlasticModel.C.

112{
113 return r.assign(1, flowPotential(stress, intnl));
114}
SolidMechanicsPlasticModel::flowPotential
virtual RankTwoTensor flowPotential(const RankTwoTensor &stress, Real intnl) const
The flow potential.
Definition SolidMechanicsPlasticModel.C:104

◆ hardPotential()

Real SolidMechanicsPlasticMeanCapTC::hardPotential ( const RankTwoTensor &  stress,
Real  intnl 
) const
overrideprotectedvirtual

The hardening potential.

Note that it is -1 for stress.trace() > _strength, and +1 for stress.trace() < _c_strength. This implements the idea that tensile failure will cause a massive reduction in compressive strength

Parameters
stressthe stress at which to calculate the hardening potential
intnlinternal parameter
Returns
the hardening potential

Reimplemented from SolidMechanicsPlasticModel.

Definition at line 169 of file SolidMechanicsPlasticMeanCapTC.C.

170{
171 // This is the key for this whole class!
172 const Real tr = stress.trace();
173 const Real t_str = tensile_strength(intnl);
174
175 if (tr >= t_str)
176 return -1.0; // this will serve to *increase* the internal parameter (so internal parameter will
177 // be a measure of volumetric strain)
178
179 const Real c_str = compressive_strength(intnl);
180 if (tr <= c_str)
181 return 1.0; // this will serve to *decrease* the internal parameter (so internal parameter will
182 // be a measure of volumetric strain)
183
184 return std::cos(libMesh::pi * (tr - c_str) /
185 (t_str - c_str)); // this interpolates C2 smoothly between 1 and -1
186}

◆ hardPotentialV()

void SolidMechanicsPlasticModel::hardPotentialV ( const RankTwoTensor &  stress,
Real  intnl,
std::vector< Real > &  h 
) const
virtualinherited

The hardening potential.

Parameters
stressthe stress at which to calculate the hardening potential
intnlinternal parameter
[out]hh[alpha] is the hardening potential for the "alpha" yield function

Definition at line 150 of file SolidMechanicsPlasticModel.C.

153{
154 h.assign(numberSurfaces(), hardPotential(stress, intnl));
155}
SolidMechanicsPlasticModel::hardPotential
virtual Real hardPotential(const RankTwoTensor &stress, Real intnl) const
The hardening potential.
Definition SolidMechanicsPlasticModel.C:145

◆ initialize()

void SolidMechanicsPlasticModel::initialize ( )
virtualinherited

Implements GeneralUserObject.

Definition at line 40 of file SolidMechanicsPlasticModel.C.

41{
42}

◆ KuhnTuckerSingleSurface()

bool SolidMechanicsPlasticModel::KuhnTuckerSingleSurface ( Real  yf,
Real  dpm,
Real  dpm_tol 
) const
inherited

Returns true if the Kuhn-Tucker conditions for the single surface are satisfied.

Parameters
yfYield function value
dpmplastic multiplier
dpm_toltolerance on plastic multiplier: viz dpm>-dpm_tol means "dpm is non-negative"

Definition at line 246 of file SolidMechanicsPlasticModel.C.

247{
248 return (dpm == 0 && yf <= _f_tol) || (dpm > -dpm_tol && yf <= _f_tol && yf >= -_f_tol);
249}

Referenced by SolidMechanicsPlasticTensileMulti::KuhnTuckerOK(), SolidMechanicsPlasticMohrCoulombMulti::KuhnTuckerOK(), and SolidMechanicsPlasticModel::returnMap().

◆ modelName()

std::string SolidMechanicsPlasticMeanCapTC::modelName ( ) const
overridevirtual

Implements SolidMechanicsPlasticModel.

Definition at line 450 of file SolidMechanicsPlasticMeanCapTC.C.

451{
452 return "MeanCapTC";
453}

◆ numberSurfaces()

unsigned SolidMechanicsPlasticModel::numberSurfaces ( ) const
virtualinherited

The number of yield surfaces for this plasticity model.

Reimplemented in SolidMechanicsPlasticMohrCoulombMulti, and SolidMechanicsPlasticTensileMulti.

Definition at line 55 of file SolidMechanicsPlasticModel.C.

56{
57 return 1;
58}

Referenced by SolidMechanicsPlasticModel::activeConstraints(), SolidMechanicsPlasticModel::dhardPotential_dintnlV(), SolidMechanicsPlasticModel::dhardPotential_dstressV(), SolidMechanicsPlasticModel::hardPotentialV(), and SolidMechanicsPlasticModel::returnMap().

◆ returnMap()

bool SolidMechanicsPlasticMeanCapTC::returnMap ( const RankTwoTensor &  trial_stress,
Real  intnl_old,
const RankFourTensor &  E_ijkl,
Real  ep_plastic_tolerance,
RankTwoTensor &  returned_stress,
Real &  returned_intnl,
std::vector< Real > &  dpm,
RankTwoTensor &  delta_dp,
std::vector< Real > &  yf,
bool &  trial_stress_inadmissible 
) const
overridevirtual

Performs a custom return-map.

You may choose to over-ride this in your derived SolidMechanicsPlasticXXXX class, and you may implement the return-map algorithm in any way that suits you. Eg, using a Newton-Raphson approach, or a radial-return, etc. This may also be used as a quick way of ascertaining whether (trial_stress, intnl_old) is in fact admissible.

For over-riding this function, please note the following.

(1) Denoting the return value of the function by "successful_return", the only possible output values should be: (A) trial_stress_inadmissible=false, successful_return=true. That is, (trial_stress, intnl_old) is in fact admissible (in the elastic domain). (B) trial_stress_inadmissible=true, successful_return=false. That is (trial_stress, intnl_old) is inadmissible (outside the yield surface), and you didn't return to the yield surface. (C) trial_stress_inadmissible=true, successful_return=true. That is (trial_stress, intnl_old) is inadmissible (outside the yield surface), but you did return to the yield surface. The default implementation only handles case (A) and (B): it does not attempt to do a return-map algorithm.

(2) you must correctly signal "successful_return" using the return value of this function. Don't assume the calling function will do Kuhn-Tucker checking and so forth!

(3) In cases (A) and (B) you needn't set returned_stress, returned_intnl, delta_dp, or dpm. This is for computational efficiency.

(4) In cases (A) and (B), you MUST place the yield function values at (trial_stress, intnl_old) into yf so the calling function can use this information optimally. You will have already calculated these yield function values, which can be quite expensive, and it's not very optimal for the calling function to have to re-calculate them.

(5) In case (C), you need to set: returned_stress (the returned value of stress) returned_intnl (the returned value of the internal variable) delta_dp (the change in plastic strain) dpm (the plastic multipliers needed to bring about the return) yf (yield function values at the returned configuration)

(Note, if you over-ride returnMap, you will probably want to override consistentTangentOpertor too, otherwise it will default to E_ijkl.)

Parameters
trial_stressThe trial stress
intnl_oldValue of the internal parameter
E_ijklElasticity tensor
ep_plastic_toleranceTolerance defined by the user for the plastic strain
[out]returned_stressIn case (C): lies on the yield surface after returning and produces the correct plastic strain (normality condition). Otherwise: not defined
[out]returned_intnlIn case (C): the value of the internal parameter after returning. Otherwise: not defined
[out]dpmIn case (C): the plastic multipliers needed to bring about the return. Otherwise: not defined
[out]delta_dpIn case (C): The change in plastic strain induced by the return process. Otherwise: not defined
[out]yfIn case (C): the yield function at (returned_stress, returned_intnl). Otherwise: the yield function at (trial_stress, intnl_old)
[out]trial_stress_inadmissibleShould be set to false if the trial_stress is admissible, and true if the trial_stress is inadmissible. This can be used by the calling prorgram
Returns
true if a successful return (or a return-map not needed), false if the trial_stress is inadmissible but the return process failed

Reimplemented from SolidMechanicsPlasticModel.

Definition at line 299 of file SolidMechanicsPlasticMeanCapTC.C.

309{
310 if (!(_use_custom_returnMap))
311 return SolidMechanicsPlasticModel::returnMap(trial_stress,
312 intnl_old,
313 E_ijkl,
314 ep_plastic_tolerance,
315 returned_stress,
316 returned_intnl,
317 dpm,
318 delta_dp,
319 yf,
320 trial_stress_inadmissible);
321
322 yf.resize(1);
323
324 Real yf_orig = yieldFunction(trial_stress, intnl_old);
325
326 yf[0] = yf_orig;
327
328 if (yf_orig < _f_tol)
329 {
330 // the trial_stress is admissible
331 trial_stress_inadmissible = false;
332 return true;
333 }
334
335 trial_stress_inadmissible = true;
336
337 // In the following we want to solve
338 // trial_stress - stress = E_ijkl * dpm * r ...... (1)
339 // and either
340 // stress.trace() = tensile_strength(intnl) ...... (2a)
341 // intnl = intnl_old + dpm ...... (3a)
342 // or
343 // stress.trace() = compressive_strength(intnl) ... (2b)
344 // intnl = intnl_old - dpm ...... (3b)
345 // The former (2a and 3a) are chosen if
346 // trial_stress.trace() > tensile_strength(intnl_old)
347 // while the latter (2b and 3b) are chosen if
348 // trial_stress.trace() < compressive_strength(intnl_old)
349 // The variables we want to solve for are stress, dpm
350 // and intnl. We do this using a Newton approach, starting
351 // with stress=trial_stress and intnl=intnl_old and dpm=0
352 const bool tensile_failure = (trial_stress.trace() >= tensile_strength(intnl_old));
353 const Real dirn = (tensile_failure ? 1.0 : -1.0);
354
355 RankTwoTensor n; // flow direction, which is E_ijkl * r
356 for (unsigned i = 0; i < 3; ++i)
357 for (unsigned j = 0; j < 3; ++j)
358 for (unsigned k = 0; k < 3; ++k)
359 n(i, j) += dirn * E_ijkl(i, j, k, k);
360 const Real n_trace = n.trace();
361
362 // Perform a Newton-Raphson to find dpm when
363 // residual = trial_stress.trace() - tensile_strength(intnl) - dpm * n.trace() [for
364 // tensile_failure=true]
365 // or
366 // residual = trial_stress.trace() - compressive_strength(intnl) - dpm * n.trace() [for
367 // tensile_failure=false]
368 Real trial_trace = trial_stress.trace();
369 Real residual;
370 Real jac;
371 dpm[0] = 0;
372 unsigned int iter = 0;
373 do
374 {
375 if (tensile_failure)
376 {
377 residual = trial_trace - tensile_strength(intnl_old + dpm[0]) - dpm[0] * n_trace;
378 jac = -dtensile_strength(intnl_old + dpm[0]) - n_trace;
379 }
380 else
381 {
382 residual = trial_trace - compressive_strength(intnl_old - dpm[0]) - dpm[0] * n_trace;
383 jac = -dcompressive_strength(intnl_old - dpm[0]) - n_trace;
384 }
385 dpm[0] += -residual / jac;
386 if (iter > _max_iters) // not converging
387 return false;
388 iter++;
389 } while (residual * residual > _f_tol * _f_tol);
390
391 // set the returned values
392 yf[0] = 0;
393 returned_intnl = intnl_old + dirn * dpm[0];
394 returned_stress = trial_stress - dpm[0] * n;
395 delta_dp = dpm[0] * dirn * returned_stress.dtrace();
396
397 return true;
398}
SolidMechanicsPlasticMeanCapTC::yieldFunction
Real yieldFunction(const RankTwoTensor &stress, Real intnl) const override
The following functions are what you should override when building single-plasticity models.
Definition SolidMechanicsPlasticMeanCapTC.C:67
SolidMechanicsPlasticModel::returnMap
virtual bool returnMap(const RankTwoTensor &trial_stress, Real intnl_old, const RankFourTensor &E_ijkl, Real ep_plastic_tolerance, RankTwoTensor &returned_stress, Real &returned_intnl, std::vector< Real > &dpm, RankTwoTensor &delta_dp, std::vector< Real > &yf, bool &trial_stress_inadmissible) const
Performs a custom return-map.
Definition SolidMechanicsPlasticModel.C:219

◆ tensile_strength()

Real SolidMechanicsPlasticMeanCapTC::tensile_strength ( const Real  internal_param) const
protectedvirtual

tensile strength as a function of residual value, rate, and internal_param

Definition at line 225 of file SolidMechanicsPlasticMeanCapTC.C.

226{
227 return _strength.value(internal_param);
228}

Referenced by activeConstraints(), consistentTangentOperator(), df_dsig(), dflowPotential_dintnl(), dflowPotential_dstress(), dhardPotential_dintnl(), dhardPotential_dstress(), dyieldFunction_dintnl(), hardPotential(), returnMap(), and yieldFunction().

◆ useCustomCTO()

bool SolidMechanicsPlasticMeanCapTC::useCustomCTO ( ) const
overridevirtual

Returns false. You will want to override this in your derived class if you write a custom consistent tangent operator function.

Reimplemented from SolidMechanicsPlasticModel.

Definition at line 444 of file SolidMechanicsPlasticMeanCapTC.C.

445{
446 return _use_custom_cto;
447}

◆ useCustomReturnMap()

bool SolidMechanicsPlasticMeanCapTC::useCustomReturnMap ( ) const
overridevirtual

Returns false. You will want to override this in your derived class if you write a custom returnMap function.

Reimplemented from SolidMechanicsPlasticModel.

Definition at line 438 of file SolidMechanicsPlasticMeanCapTC.C.

439{
440 return _use_custom_returnMap;
441}

◆ validParams()

InputParameters SolidMechanicsPlasticMeanCapTC::validParams ( )
static

Definition at line 21 of file SolidMechanicsPlasticMeanCapTC.C.

22{
23 InputParameters params = SolidMechanicsPlasticModel::validParams();
24 params.addRangeCheckedParam<unsigned>("max_iterations",
25 10,
26 "max_iterations>0",
27 "Maximum iterations for custom MeanCapTC return map");
28 params.addParam<bool>(
29 "use_custom_returnMap", true, "Whether to use the custom MeanCapTC returnMap algorithm.");
30 params.addParam<bool>("use_custom_cto",
31 true,
32 "Whether to use the custom consistent tangent operator computations.");
33 params.addRequiredParam<UserObjectName>("tensile_strength",
34 "A SolidMechanicsHardening UserObject that defines "
35 "hardening of the mean-cap tensile strength (it will "
36 "typically be positive). Yield function = trace(stress) "
37 "- tensile_strength for trace(stress)>tensile_strength.");
38 params.addRequiredParam<UserObjectName>(
39 "compressive_strength",
40 "A SolidMechanicsHardening UserObject that defines hardening of the "
41 "mean-cap compressive strength. This should always be less than "
42 "tensile_strength (it will typically be negative). Yield function = "
43 "- (trace(stress) - compressive_strength) for "
44 "trace(stress)<compressive_strength.");
45 params.addClassDescription(
46 "Associative mean-cap tensile and compressive plasticity with hardening/softening");
47
48 return params;
49}
InputParameters::addRequiredParam
void addRequiredParam(const std::string &name, const std::string &doc_string)
InputParameters::addParam
void addParam(const std::string &name, const std::initializer_list< typename T::value_type > &value, const std::string &doc_string)
InputParameters::addClassDescription
void addClassDescription(const std::string &doc_string)
InputParameters::addRangeCheckedParam
void addRangeCheckedParam(const std::string &name, const T &value, const std::string &parsed_function, const std::string &doc_string)
SolidMechanicsPlasticModel::validParams
static InputParameters validParams()
Definition SolidMechanicsPlasticModel.C:16

◆ yieldFunction()

Real SolidMechanicsPlasticMeanCapTC::yieldFunction ( const RankTwoTensor &  stress,
Real  intnl 
) const
overrideprotectedvirtual

The following functions are what you should override when building single-plasticity models.

The yield function

Parameters
stressthe stress at which to calculate the yield function
intnlinternal parameter
Returns
the yield function

Reimplemented from SolidMechanicsPlasticModel.

Definition at line 67 of file SolidMechanicsPlasticMeanCapTC.C.

68{
69 const Real tr = stress.trace();
70 const Real t_str = tensile_strength(intnl);
71 if (tr >= t_str)
72 return tr - t_str;
73
74 const Real c_str = compressive_strength(intnl);
75 if (tr <= c_str)
76 return -(tr - c_str);
77 // the following is zero at tr = t_str, and at tr = c_str
78 // it also has derivative = 1 at tr = t_str, and derivative = -1 at tr = c_str
79 // it also has second derivative = 0, at these points.
80 // This makes the complete yield function C2 continuous.
81 return (c_str - t_str) / libMesh::pi * std::sin(libMesh::pi * (tr - c_str) / (t_str - c_str));
82}

Referenced by returnMap().

◆ yieldFunctionV()

void SolidMechanicsPlasticModel::yieldFunctionV ( const RankTwoTensor &  stress,
Real  intnl,
std::vector< Real > &  f 
) const
virtualinherited

Calculates the yield functions.

Note that for single-surface plasticity you don't want to override this - override the private yieldFunction below

Parameters
stressthe stress at which to calculate the yield function
intnlinternal parameter
[out]fthe yield functions

Reimplemented in SolidMechanicsPlasticMohrCoulombMulti, and SolidMechanicsPlasticTensileMulti.

Definition at line 67 of file SolidMechanicsPlasticModel.C.

70{
71 f.assign(1, yieldFunction(stress, intnl));
72}
SolidMechanicsPlasticModel::yieldFunction
virtual Real yieldFunction(const RankTwoTensor &stress, Real intnl) const
The following functions are what you should override when building single-plasticity models.
Definition SolidMechanicsPlasticModel.C:61

Referenced by SolidMechanicsPlasticModel::returnMap().

Member Data Documentation

◆ _c_strength

const SolidMechanicsHardeningModel& SolidMechanicsPlasticMeanCapTC::_c_strength
protected

the compressive strength

Definition at line 77 of file SolidMechanicsPlasticMeanCapTC.h.

Referenced by compressive_strength(), dcompressive_strength(), and SolidMechanicsPlasticMeanCapTC().

◆ _f_tol

const Real SolidMechanicsPlasticModel::_f_tol
inherited

Tolerance on yield function.

Definition at line 170 of file SolidMechanicsPlasticModel.h.

Referenced by SolidMechanicsPlasticModel::activeConstraints(), activeConstraints(), SolidMechanicsPlasticMohrCoulombMulti::activeConstraints(), SolidMechanicsPlasticTensileMulti::activeConstraints(), SolidMechanicsPlasticWeakPlaneShear::activeConstraints(), SolidMechanicsPlasticWeakPlaneTensile::activeConstraints(), SolidMechanicsPlasticMohrCoulombMulti::doReturnMap(), SolidMechanicsPlasticTensileMulti::doReturnMap(), SolidMechanicsPlasticModel::KuhnTuckerSingleSurface(), SolidMechanicsPlasticTensileMulti::returnEdge(), SolidMechanicsPlasticMohrCoulombMulti::returnEdge000101(), SolidMechanicsPlasticMohrCoulombMulti::returnEdge010100(), SolidMechanicsPlasticModel::returnMap(), SolidMechanicsPlasticDruckerPragerHyperbolic::returnMap(), SolidMechanicsPlasticJ2::returnMap(), returnMap(), SolidMechanicsPlasticMohrCoulombMulti::returnPlane(), SolidMechanicsPlasticTensileMulti::returnPlane(), SolidMechanicsPlasticMohrCoulombMulti::returnTip(), SolidMechanicsPlasticTensileMulti::returnTip(), SolidMechanicsPlasticMohrCoulombMulti::SolidMechanicsPlasticMohrCoulombMulti(), and SolidMechanicsPlasticTensileMulti::SolidMechanicsPlasticTensileMulti().

◆ _ic_tol

const Real SolidMechanicsPlasticModel::_ic_tol
inherited

Tolerance on internal constraint.

Definition at line 173 of file SolidMechanicsPlasticModel.h.

◆ _max_iters

const unsigned SolidMechanicsPlasticMeanCapTC::_max_iters
protected

max iters for custom return map loop

Definition at line 65 of file SolidMechanicsPlasticMeanCapTC.h.

Referenced by returnMap().

◆ _strength

const SolidMechanicsHardeningModel& SolidMechanicsPlasticMeanCapTC::_strength
protected

the tensile strength

Definition at line 74 of file SolidMechanicsPlasticMeanCapTC.h.

Referenced by dtensile_strength(), SolidMechanicsPlasticMeanCapTC(), and tensile_strength().

◆ _use_custom_cto

const bool SolidMechanicsPlasticMeanCapTC::_use_custom_cto
protected

Whether to use the custom consistent tangent operator algorithm.

Definition at line 71 of file SolidMechanicsPlasticMeanCapTC.h.

Referenced by consistentTangentOperator(), and useCustomCTO().

◆ _use_custom_returnMap

const bool SolidMechanicsPlasticMeanCapTC::_use_custom_returnMap
protected

Whether to use the custom return-map algorithm.

Definition at line 68 of file SolidMechanicsPlasticMeanCapTC.h.

Referenced by returnMap(), and useCustomReturnMap().

The documentation for this class was generated from the following files:
Generated on Fri Aug 21 2026 13:43:53 for https://mooseframework.inl.gov by doxygen 1.9.8