Rate-independent non-associative Drucker Prager with hardening/softening. More...

#include <SolidMechanicsPlasticDruckerPragerHyperbolic.h>

Inheritance diagram for SolidMechanicsPlasticDruckerPragerHyperbolic:

Public Types
enum	FrictionDilation { friction = 0 , dilation = 1 }
	bbb (friction) and bbb_flow (dilation) are computed using the same function, onlyB, and this parameter tells that function whether to compute bbb or bbb_flow More...

typedef DataFileName	DataFileParameterType

Public Member Functions
	SolidMechanicsPlasticDruckerPragerHyperbolic (const InputParameters &parameters)

virtual std::string	modelName () const override

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.

void	bothAB (Real intnl, Real &aaa, Real &bbb) const
	Calculates aaa and bbb as a function of the internal parameter intnl.

void	dbothAB (Real intnl, Real &daaa, Real &dbbb) const
	Calculates d(aaa)/d(intnl) and d(bbb)/d(intnl) as a function of the internal parameter intnl.

void	onlyB (Real intnl, int fd, Real &bbb) const
	Calculate bbb or bbb_flow.

void	donlyB (Real intnl, int fd, Real &dbbb) const
	Calculate d(bbb)/d(intnl) or d(bbb_flow)/d(intnl)

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.

virtual void	activeConstraints (const std::vector< Real > &f, const RankTwoTensor &stress, Real intnl, const RankFourTensor &Eijkl, std::vector< bool > &act, RankTwoTensor &returned_stress) const
	The active yield surfaces, given a vector of yield functions.

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.

RankFourTensor	dflowPotential_dstress (const RankTwoTensor &stress, Real intnl) const override
	The derivative of the flow potential with respect to stress.

virtual RankTwoTensor	df_dsig (const RankTwoTensor &stress, Real bbb) const override
	Function that's used in dyieldFunction_dstress and flowPotential.

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 RankTwoTensor	dyieldFunction_dstress (const RankTwoTensor &stress, Real intnl) const override
	The derivative of yield function with respect to stress.

virtual Real	dyieldFunction_dintnl (const RankTwoTensor &stress, Real intnl) const override
	The derivative of yield function with respect to the internal parameter.

virtual RankTwoTensor	flowPotential (const RankTwoTensor &stress, Real intnl) const override
	The flow potential.

virtual RankTwoTensor	dflowPotential_dintnl (const RankTwoTensor &stress, Real intnl) const override
	The derivative of the flow potential with respect to the internal parameter.

virtual Real	hardPotential (const RankTwoTensor &stress, Real intnl) const
	The hardening potential.

virtual RankTwoTensor	dhardPotential_dstress (const RankTwoTensor &stress, Real intnl) const
	The derivative of the hardening potential with respect to stress.

virtual Real	dhardPotential_dintnl (const RankTwoTensor &stress, Real intnl) const
	The derivative of the hardening potential with respect to the internal parameter.

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 SolidMechanicsHardeningModel &	_mc_cohesion
	Hardening model for cohesion.

const SolidMechanicsHardeningModel &	_mc_phi
	Hardening model for tan(phi)

const SolidMechanicsHardeningModel &	_mc_psi
	Hardening model for tan(psi)

const MooseEnum	_mc_interpolation_scheme
	The parameters aaa and bbb are chosen to closely match the Mohr-Coulomb yield surface.

const bool	_zero_cohesion_hardening
	True if there is no hardening of cohesion.

const bool	_zero_phi_hardening
	True if there is no hardening of friction angle.

const bool	_zero_psi_hardening
	True if there is no hardening of dilation angle.

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
void	initializeAandB (Real intnl, Real &aaa, Real &bbb) const
	Returns the Drucker-Prager parameters A nice reference on the different relationships between Drucker-Prager and Mohr-Coulomb is H Jiang and X Yongli, A note on the Mohr-Coulomb and Drucker-Prager strength criteria, Mechanics Research Communications 38 (2011) 309-314.

void	initializeB (Real intnl, int fd, Real &bbb) const
	Returns the Drucker-Prager parameters A nice reference on the different relationships between Drucker-Prager and Mohr-Coulomb is H Jiang and X Yongli, A note on the Mohr-Coulomb and Drucker-Prager strength criteria, Mechanics Research Communications 38 (2011) 309-314.

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
const Real	_smoother2
	smoothing parameter for the cone's tip

const bool	_use_custom_returnMap
	whether to use the custom returnMap function

const bool	_use_custom_cto
	Whether to use the custom consistent tangent operator calculation.

const unsigned	_max_iters
	max iters for custom return map loop

Real	_aaa

Real	_bbb

Real	_bbb_flow

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 non-associative Drucker Prager with hardening/softening.

The cone's tip is smoothed in a hyperbolic fashion Most functions (eg flowPotential, etc) are simply inherited from SolidMechanicsPlasticDruckerPrager. Note df_dsig is over-ridden

Definition at line 21 of file SolidMechanicsPlasticDruckerPragerHyperbolic.h.

Member Enumeration Documentation

◆ FrictionDilation

enum SolidMechanicsPlasticDruckerPrager::FrictionDilation

inherited

bbb (friction) and bbb_flow (dilation) are computed using the same function, onlyB, and this parameter tells that function whether to compute bbb or bbb_flow

Enumerator
friction
dilation

Definition at line 41 of file SolidMechanicsPlasticDruckerPrager.h.

  {
    friction = 0,
    dilation = 1
  };

Constructor & Destructor Documentation

◆ SolidMechanicsPlasticDruckerPragerHyperbolic()

SolidMechanicsPlasticDruckerPragerHyperbolic::SolidMechanicsPlasticDruckerPragerHyperbolic ( const InputParameters & parameters )

Definition at line 53 of file SolidMechanicsPlasticDruckerPragerHyperbolic.C.

  : SolidMechanicsPlasticDruckerPrager(parameters),
    _smoother2(Utility::pow<2>(getParam<Real>("smoother"))),
    _use_custom_returnMap(getParam<bool>("use_custom_returnMap")),
    _use_custom_cto(getParam<bool>("use_custom_cto")),
    _max_iters(getParam<unsigned>("max_iterations"))
{
}

Member Function Documentation

◆ activeConstraints()

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

virtualinherited

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

	f	values of the yield functions
	stress	stress tensor
	intnl	internal parameter
	Eijkl	elasticity tensor (stress = Eijkl*strain)
[out]	act	act[i] = true if the i_th yield function is active
[out]	returned_stress	Approximate value of the returned stress

Reimplemented in SolidMechanicsPlasticMeanCapTC, SolidMechanicsPlasticMohrCoulombMulti, SolidMechanicsPlasticTensileMulti, SolidMechanicsPlasticWeakPlaneShear, and SolidMechanicsPlasticWeakPlaneTensile.

Definition at line 186 of file SolidMechanicsPlasticModel.C.

{
  mooseAssert(f.size() == numberSurfaces(),
              "f incorrectly sized at " << f.size() << " in activeConstraints");
  act.resize(numberSurfaces());
  for (unsigned surface = 0; surface < numberSurfaces(); ++surface)
    act[surface] = (f[surface] > _f_tol);
}

◆ bothAB()

void SolidMechanicsPlasticDruckerPrager::bothAB	(	Real	intnl,
		Real &	aaa,
		Real &	bbb
	)		const

inherited

Calculates aaa and bbb as a function of the internal parameter intnl.

Definition at line 151 of file SolidMechanicsPlasticDruckerPrager.C.

{
  if (_zero_cohesion_hardening && _zero_phi_hardening)
  {
    aaa = _aaa;
    bbb = _bbb;
    return;
  }
  initializeAandB(intnl, aaa, bbb);
}

Referenced by CappedDruckerPragerStressUpdate::computeAllQ(), CappedDruckerPragerStressUpdate::initializeVars(), returnMap(), SolidMechanicsPlasticDruckerPrager::yieldFunction(), yieldFunction(), and CappedDruckerPragerStressUpdate::yieldFunctionValues().

◆ consistentTangentOperator()

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

overrideprotectedvirtual

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_old	trial stress before returning
intnl_old	internal parameter before returning
stress	current returned stress state
intnl	internal parameter
E_ijkl	elasticity tensor
cumulative_pm	the cumulative plastic multipliers

Returns: the consistent tangent operator: E_ijkl if not over-ridden

Reimplemented from SolidMechanicsPlasticModel.

Definition at line 203 of file SolidMechanicsPlasticDruckerPragerHyperbolic.C.

{
  if (!_use_custom_cto)
    return SolidMechanicsPlasticModel::consistentTangentOperator(
        trial_stress, intnl_old, stress, intnl, E_ijkl, cumulative_pm);
 
  const Real mu = E_ijkl(0, 1, 0, 1);
  const Real la = E_ijkl(0, 0, 0, 0) - 2.0 * mu;
  const Real bulky = 3.0 * la + 2.0 * mu;
  Real bbb;
  onlyB(intnl, friction, bbb);
  Real bbb_flow;
  onlyB(intnl, dilation, bbb_flow);
  Real dbbb_flow;
  donlyB(intnl, dilation, dbbb_flow);
  const Real bbb_flow_mod = bbb_flow + cumulative_pm[0] * dbbb_flow;
  const Real J2 = stress.secondInvariant();
  const RankTwoTensor sij = stress.deviatoric();
  const Real sq = std::sqrt(J2 + _smoother2);
 
  const Real one_over_h =
      1.0 / (-dyieldFunction_dintnl(stress, intnl) + mu * J2 / Utility::pow<2>(sq) +
             3.0 * bbb * bbb_flow_mod * bulky); // correct up to hard
 
  const RankTwoTensor df_dsig_timesE =
      mu * sij / sq + bulky * bbb * RankTwoTensor(RankTwoTensor::initIdentity); // correct
  const RankTwoTensor rmod_timesE =
      mu * sij / sq +
      bulky * bbb_flow_mod * RankTwoTensor(RankTwoTensor::initIdentity); // correct up to hard
 
  const RankFourTensor coeff_ep =
      E_ijkl - one_over_h * rmod_timesE.outerProduct(df_dsig_timesE); // correct
 
  const RankFourTensor dr_dsig_timesE = -0.5 * mu * sij.outerProduct(sij) / Utility::pow<3>(sq) +
                                        mu * stress.d2secondInvariant() / sq; // correct
  const RankTwoTensor df_dsig_E_dr_dsig =
      0.5 * mu * _smoother2 * sij / Utility::pow<4>(sq); // correct
 
  const RankFourTensor coeff_si =
      RankFourTensor(RankFourTensor::initIdentitySymmetricFour) +
      cumulative_pm[0] *
          (dr_dsig_timesE - one_over_h * rmod_timesE.outerProduct(df_dsig_E_dr_dsig));
 
  RankFourTensor s_inv;
  try
  {
    s_inv = coeff_si.invSymm();
  }
  catch (MooseException & e)
  {
    return coeff_ep; // when coeff_si is singular return the "linear" tangent operator
  }
 
  return s_inv * coeff_ep;
}

◆ dbothAB()

void SolidMechanicsPlasticDruckerPrager::dbothAB	(	Real	intnl,
		Real &	daaa,
		Real &	dbbb
	)		const

inherited

Calculates d(aaa)/d(intnl) and d(bbb)/d(intnl) as a function of the internal parameter intnl.

Definition at line 221 of file SolidMechanicsPlasticDruckerPrager.C.

{
  if (_zero_cohesion_hardening && _zero_phi_hardening)
  {
    daaa = 0;
    dbbb = 0;
    return;
  }
 
  const Real C = _mc_cohesion.value(intnl);
  const Real dC = _mc_cohesion.derivative(intnl);
  const Real cosphi = std::cos(_mc_phi.value(intnl));
  const Real dcosphi = -std::sin(_mc_phi.value(intnl)) * _mc_phi.derivative(intnl);
  const Real sinphi = std::sin(_mc_phi.value(intnl));
  const Real dsinphi = std::cos(_mc_phi.value(intnl)) * _mc_phi.derivative(intnl);
  switch (_mc_interpolation_scheme)
  {
    case 0: // outer_tip
      daaa = 2.0 * std::sqrt(3.0) *
             (dC * cosphi / (3.0 - sinphi) + C * dcosphi / (3.0 - sinphi) +
              C * cosphi * dsinphi / Utility::pow<2>(3.0 - sinphi));
      dbbb = 2.0 / std::sqrt(3.0) *
             (dsinphi / (3.0 - sinphi) + sinphi * dsinphi / Utility::pow<2>(3.0 - sinphi));
      break;
    case 1: // inner_tip
      daaa = 2.0 * std::sqrt(3.0) *
             (dC * cosphi / (3.0 + sinphi) + C * dcosphi / (3.0 + sinphi) -
              C * cosphi * dsinphi / Utility::pow<2>(3.0 + sinphi));
      dbbb = 2.0 / std::sqrt(3.0) *
             (dsinphi / (3.0 + sinphi) - sinphi * dsinphi / Utility::pow<2>(3.0 + sinphi));
      break;
    case 2: // lode_zero
      daaa = dC * cosphi + C * dcosphi;
      dbbb = dsinphi / 3.0;
      break;
    case 3: // inner_edge
      daaa = 3.0 * dC * cosphi / std::sqrt(9.0 + 3.0 * Utility::pow<2>(sinphi)) +
             3.0 * C * dcosphi / std::sqrt(9.0 + 3.0 * Utility::pow<2>(sinphi)) -
             3.0 * C * cosphi * 3.0 * sinphi * dsinphi /
                 std::pow(9.0 + 3.0 * Utility::pow<2>(sinphi), 1.5);
      dbbb = dsinphi / std::sqrt(9.0 + 3.0 * Utility::pow<2>(sinphi)) -
             3.0 * sinphi * sinphi * dsinphi / std::pow(9.0 + 3.0 * Utility::pow<2>(sinphi), 1.5);
      break;
    case 4: // native
      daaa = dC;
      dbbb = dsinphi / cosphi - sinphi * dcosphi / Utility::pow<2>(cosphi);
      break;
  }
}

Referenced by CappedDruckerPragerStressUpdate::computeAllQ(), SolidMechanicsPlasticDruckerPrager::dyieldFunction_dintnl(), and returnMap().

◆ df_dsig()

RankTwoTensor SolidMechanicsPlasticDruckerPragerHyperbolic::df_dsig	(	const RankTwoTensor &	stress,
		Real	bbb
	)		const

overrideprotectedvirtual

Function that's used in dyieldFunction_dstress and flowPotential.

Reimplemented from SolidMechanicsPlasticDruckerPrager.

Definition at line 74 of file SolidMechanicsPlasticDruckerPragerHyperbolic.C.

{
  return 0.5 * stress.dsecondInvariant() / std::sqrt(stress.secondInvariant() + _smoother2) +
         stress.dtrace() * bbb;
}

◆ dflowPotential_dintnl()

RankTwoTensor SolidMechanicsPlasticDruckerPrager::dflowPotential_dintnl	(	const RankTwoTensor &	stress,
		Real	intnl
	)		const

overrideprotectedvirtualinherited

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

Parameters

stress	the stress at which to calculate the flow potential
intnl	internal parameter

Returns: dr_dintnl(i, j) = dr(i, j)/dintnl

Reimplemented from SolidMechanicsPlasticModel.

Definition at line 136 of file SolidMechanicsPlasticDruckerPrager.C.

{
  Real dbbb;
  donlyB(intnl, dilation, dbbb);
  return stress.dtrace() * dbbb;
}

◆ 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

	stress	the stress at which to calculate the flow potential
	intnl	internal parameter
[out]	dr_dintnl	dr_dintnl[alpha](i, j) = dr[alpha](i, j)/dintnl

Reimplemented in SolidMechanicsPlasticMohrCoulombMulti, and SolidMechanicsPlasticTensileMulti.

Definition at line 137 of file SolidMechanicsPlasticModel.C.

{
  return dr_dintnl.assign(1, dflowPotential_dintnl(stress, intnl));
}

◆ dflowPotential_dstress()

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

overrideprotectedvirtual

The derivative of the flow potential with respect to stress.

Parameters

stress	the stress at which to calculate the flow potential
intnl	internal parameter

Returns: dr_dstress(i, j, k, l) = dr(i, j)/dstress(k, l)

Reimplemented from SolidMechanicsPlasticDruckerPrager.

Definition at line 81 of file SolidMechanicsPlasticDruckerPragerHyperbolic.C.

{
  RankFourTensor dr_dstress;
  dr_dstress = 0.5 * stress.d2secondInvariant() / std::sqrt(stress.secondInvariant() + _smoother2);
  dr_dstress += -0.5 * 0.5 * stress.dsecondInvariant().outerProduct(stress.dsecondInvariant()) /
                std::pow(stress.secondInvariant() + _smoother2, 1.5);
  return dr_dstress;
}

◆ 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

	stress	the stress at which to calculate the flow potential
	intnl	internal parameter
[out]	dr_dstress	dr_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.

{
  return dr_dstress.assign(1, dflowPotential_dstress(stress, intnl));
}

◆ dhardPotential_dintnl()

Real SolidMechanicsPlasticModel::dhardPotential_dintnl	(	const RankTwoTensor &	stress,
		Real	intnl
	)		const

protectedvirtualinherited

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

Parameters

stress	the stress at which to calculate the hardening potentials
intnl	internal parameter

Returns: the derivative

Reimplemented in SolidMechanicsPlasticMeanCapTC.

Definition at line 172 of file SolidMechanicsPlasticModel.C.

{
  return 0.0;
}

Referenced by SolidMechanicsPlasticModel::dhardPotential_dintnlV().

◆ 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

	stress	the stress at which to calculate the hardening potentials
	intnl	internal parameter
[out]	dh_dintnl	dh_dintnl[alpha] = dh[alpha]/dintnl

Definition at line 178 of file SolidMechanicsPlasticModel.C.

{
  dh_dintnl.resize(numberSurfaces(), dhardPotential_dintnl(stress, intnl));
}

◆ dhardPotential_dstress()

RankTwoTensor SolidMechanicsPlasticModel::dhardPotential_dstress	(	const RankTwoTensor &	stress,
		Real	intnl
	)		const

protectedvirtualinherited

The derivative of the hardening potential with respect to stress.

Parameters

stress	the stress at which to calculate the hardening potentials
intnl	internal parameter

Returns: dh_dstress(i, j) = dh/dstress(i, j)

Reimplemented in SolidMechanicsPlasticMeanCapTC.

Definition at line 158 of file SolidMechanicsPlasticModel.C.

{
  return RankTwoTensor();
}

Referenced by SolidMechanicsPlasticModel::dhardPotential_dstressV().

◆ 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

	stress	the stress at which to calculate the hardening potentials
	intnl	internal parameter
[out]	dh_dstress	dh_dstress[alpha](i, j) = dh[alpha]/dstress(i, j)

Definition at line 164 of file SolidMechanicsPlasticModel.C.

{
  dh_dstress.assign(numberSurfaces(), dhardPotential_dstress(stress, intnl));
}

◆ donlyB()

void SolidMechanicsPlasticDruckerPrager::donlyB	(	Real	intnl,
		int	fd,
		Real &	dbbb
	)		const

inherited

Calculate d(bbb)/d(intnl) or d(bbb_flow)/d(intnl)

Parameters

intnl	the internal parameter
fd	if fd==friction then bbb is calculated. if fd==dilation then bbb_flow is calculated
bbb	either bbb or bbb_flow, depending on fd

Definition at line 179 of file SolidMechanicsPlasticDruckerPrager.C.

{
  if (_zero_phi_hardening && (fd == friction))
  {
    dbbb = 0;
    return;
  }
  if (_zero_psi_hardening && (fd == dilation))
  {
    dbbb = 0;
    return;
  }
  const Real s = (fd == friction) ? std::sin(_mc_phi.value(intnl)) : std::sin(_mc_psi.value(intnl));
  const Real ds = (fd == friction) ? std::cos(_mc_phi.value(intnl)) * _mc_phi.derivative(intnl)
                                   : std::cos(_mc_psi.value(intnl)) * _mc_psi.derivative(intnl);
  switch (_mc_interpolation_scheme)
  {
    case 0: // outer_tip
      dbbb = 2.0 / std::sqrt(3.0) * (ds / (3.0 - s) + s * ds / Utility::pow<2>(3.0 - s));
      break;
    case 1: // inner_tip
      dbbb = 2.0 / std::sqrt(3.0) * (ds / (3.0 + s) - s * ds / Utility::pow<2>(3.0 + s));
      break;
    case 2: // lode_zero
      dbbb = ds / 3.0;
      break;
    case 3: // inner_edge
      dbbb = ds / std::sqrt(9.0 + 3.0 * Utility::pow<2>(s)) -
             3 * s * s * ds / std::pow(9.0 + 3.0 * Utility::pow<2>(s), 1.5);
      break;
    case 4: // native
      const Real c =
          (fd == friction) ? std::cos(_mc_phi.value(intnl)) : std::cos(_mc_psi.value(intnl));
      const Real dc = (fd == friction)
                          ? -std::sin(_mc_phi.value(intnl)) * _mc_phi.derivative(intnl)
                          : -std::sin(_mc_psi.value(intnl)) * _mc_psi.derivative(intnl);
      dbbb = ds / c - s * dc / Utility::pow<2>(c);
      break;
  }
}

Referenced by CappedDruckerPragerStressUpdate::computeAllQ(), consistentTangentOperator(), SolidMechanicsPlasticDruckerPrager::dflowPotential_dintnl(), returnMap(), and CappedDruckerPragerStressUpdate::setIntnlDerivatives().

◆ dyieldFunction_dintnl()

Real SolidMechanicsPlasticDruckerPrager::dyieldFunction_dintnl	(	const RankTwoTensor &	stress,
		Real	intnl
	)		const

overrideprotectedvirtualinherited

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

Parameters

stress	the stress at which to calculate the yield function
intnl	internal parameter

Returns: the derivative

Reimplemented from SolidMechanicsPlasticModel.

Definition at line 107 of file SolidMechanicsPlasticDruckerPrager.C.

{
  Real daaa;
  Real dbbb;
  dbothAB(intnl, daaa, dbbb);
  return stress.trace() * dbbb - daaa;
}

Referenced by consistentTangentOperator().

◆ 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

	stress	the stress at which to calculate the yield function
	intnl	internal parameter
[out]	df_dintnl	df_dintnl[alpha] = df[alpha]/dintnl

Reimplemented in SolidMechanicsPlasticMohrCoulombMulti, and SolidMechanicsPlasticTensileMulti.

Definition at line 96 of file SolidMechanicsPlasticModel.C.

{
  return df_dintnl.assign(1, dyieldFunction_dintnl(stress, intnl));
}

◆ dyieldFunction_dstress()

RankTwoTensor SolidMechanicsPlasticDruckerPrager::dyieldFunction_dstress	(	const RankTwoTensor &	stress,
		Real	intnl
	)		const

overrideprotectedvirtualinherited

The derivative of yield function with respect to stress.

Parameters

stress	the stress at which to calculate the yield function
intnl	internal parameter

Returns: df_dstress(i, j) = dyieldFunction/dstress(i, j)

Reimplemented from SolidMechanicsPlasticModel.

Definition at line 98 of file SolidMechanicsPlasticDruckerPrager.C.

{
  Real bbb;
  onlyB(intnl, friction, bbb);
  return df_dsig(stress, bbb);
}

◆ 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

	stress	the stress at which to calculate the yield function
	intnl	internal parameter
[out]	df_dstress	df_dstress[alpha](i, j) = dyieldFunction[alpha]/dstress(i, j)

Reimplemented in SolidMechanicsPlasticMohrCoulombMulti, and SolidMechanicsPlasticTensileMulti.

Definition at line 82 of file SolidMechanicsPlasticModel.C.

{
  df_dstress.assign(1, dyieldFunction_dstress(stress, intnl));
}

◆ 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 SolidMechanicsPlasticDruckerPrager::flowPotential	(	const RankTwoTensor &	stress,
		Real	intnl
	)		const

overrideprotectedvirtualinherited

The flow potential.

Parameters

stress	the stress at which to calculate the flow potential
intnl	internal parameter

Returns: the flow potential

Reimplemented from SolidMechanicsPlasticModel.

Definition at line 117 of file SolidMechanicsPlasticDruckerPrager.C.

{
  Real bbb_flow;
  onlyB(intnl, dilation, bbb_flow);
  return df_dsig(stress, bbb_flow);
}

◆ flowPotentialV()

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

virtualinherited

The flow potentials.

Parameters

	stress	the stress at which to calculate the flow potential
	intnl	internal parameter
[out]	r	r[alpha] is the flow potential for the "alpha" yield function

Reimplemented in SolidMechanicsPlasticMohrCoulombMulti, and SolidMechanicsPlasticTensileMulti.

Definition at line 109 of file SolidMechanicsPlasticModel.C.

{
  return r.assign(1, flowPotential(stress, intnl));
}

◆ hardPotential()

Real SolidMechanicsPlasticModel::hardPotential	(	const RankTwoTensor &	stress,
		Real	intnl
	)		const

protectedvirtualinherited

The hardening potential.

Parameters

stress	the stress at which to calculate the hardening potential
intnl	internal parameter

Returns: the hardening potential

Reimplemented in SolidMechanicsPlasticMeanCapTC.

Definition at line 145 of file SolidMechanicsPlasticModel.C.

{
  return -1.0;
}

Referenced by SolidMechanicsPlasticModel::hardPotentialV().

◆ hardPotentialV()

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

virtualinherited

The hardening potential.

Parameters

	stress	the stress at which to calculate the hardening potential
	intnl	internal parameter
[out]	h	h[alpha] is the hardening potential for the "alpha" yield function

Definition at line 150 of file SolidMechanicsPlasticModel.C.

{
  h.assign(numberSurfaces(), hardPotential(stress, intnl));
}

◆ initialize()

void SolidMechanicsPlasticModel::initialize ( )

virtualinherited

Implements GeneralUserObject.

Definition at line 40 of file SolidMechanicsPlasticModel.C.

41{

42}

◆ initializeAandB()

void SolidMechanicsPlasticDruckerPrager::initializeAandB	(	Real	intnl,
		Real &	aaa,
		Real &	bbb
	)		const

privateinherited

Returns the Drucker-Prager parameters A nice reference on the different relationships between Drucker-Prager and Mohr-Coulomb is H Jiang and X Yongli, A note on the Mohr-Coulomb and Drucker-Prager strength criteria, Mechanics Research Communications 38 (2011) 309-314.

This function uses Table1 in that reference, with the following translations aaa = k bbb = -alpha/3

Parameters

	intnl	The internal parameter
[out]	aaa	The Drucker-Prager aaa quantity
[out]	bbb	The Drucker-Prager bbb quantity

Definition at line 272 of file SolidMechanicsPlasticDruckerPrager.C.

{
  const Real C = _mc_cohesion.value(intnl);
  const Real cosphi = std::cos(_mc_phi.value(intnl));
  const Real sinphi = std::sin(_mc_phi.value(intnl));
  switch (_mc_interpolation_scheme)
  {
    case 0: // outer_tip
      aaa = 2.0 * std::sqrt(3.0) * C * cosphi / (3.0 - sinphi);
      bbb = 2.0 * sinphi / std::sqrt(3.0) / (3.0 - sinphi);
      break;
    case 1: // inner_tip
      aaa = 2.0 * std::sqrt(3.0) * C * cosphi / (3.0 + sinphi);
      bbb = 2.0 * sinphi / std::sqrt(3.0) / (3.0 + sinphi);
      break;
    case 2: // lode_zero
      aaa = C * cosphi;
      bbb = sinphi / 3.0;
      break;
    case 3: // inner_edge
      aaa = 3.0 * C * cosphi / std::sqrt(9.0 + 3.0 * Utility::pow<2>(sinphi));
      bbb = sinphi / std::sqrt(9.0 + 3.0 * Utility::pow<2>(sinphi));
      break;
    case 4: // native
      aaa = C;
      bbb = sinphi / cosphi;
      break;
  }
}

Referenced by SolidMechanicsPlasticDruckerPrager::bothAB(), and SolidMechanicsPlasticDruckerPrager::SolidMechanicsPlasticDruckerPrager().

◆ initializeB()

void SolidMechanicsPlasticDruckerPrager::initializeB	(	Real	intnl,
		int	fd,
		Real &	bbb
	)		const

privateinherited

Returns the Drucker-Prager parameters A nice reference on the different relationships between Drucker-Prager and Mohr-Coulomb is H Jiang and X Yongli, A note on the Mohr-Coulomb and Drucker-Prager strength criteria, Mechanics Research Communications 38 (2011) 309-314.

This function uses Table1 in that reference, with the following translations aaa = k bbb = -alpha/3

Parameters

	intnl	The internal parameter
	fd	If fd == frction then the friction angle is used to set bbb, otherwise the dilation angle is used
[out]	bbb	The Drucker-Prager bbb quantity

Definition at line 303 of file SolidMechanicsPlasticDruckerPrager.C.

{
  const Real s = (fd == friction) ? std::sin(_mc_phi.value(intnl)) : std::sin(_mc_psi.value(intnl));
  switch (_mc_interpolation_scheme)
  {
    case 0: // outer_tip
      bbb = 2.0 * s / std::sqrt(3.0) / (3.0 - s);
      break;
    case 1: // inner_tip
      bbb = 2.0 * s / std::sqrt(3.0) / (3.0 + s);
      break;
    case 2: // lode_zero
      bbb = s / 3.0;
      break;
    case 3: // inner_edge
      bbb = s / std::sqrt(9.0 + 3.0 * Utility::pow<2>(s));
      break;
    case 4: // native
      const Real c =
          (fd == friction) ? std::cos(_mc_phi.value(intnl)) : std::cos(_mc_psi.value(intnl));
      bbb = s / c;
      break;
  }
}

Referenced by SolidMechanicsPlasticDruckerPrager::onlyB(), and SolidMechanicsPlasticDruckerPrager::SolidMechanicsPlasticDruckerPrager().

◆ 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

yf	Yield function value
dpm	plastic multiplier
dpm_tol	tolerance on plastic multiplier: viz dpm>-dpm_tol means "dpm is non-negative"

Definition at line 246 of file SolidMechanicsPlasticModel.C.

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

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

◆ modelName()

std::string SolidMechanicsPlasticDruckerPragerHyperbolic::modelName ( ) const

overridevirtual

Reimplemented from SolidMechanicsPlasticDruckerPrager.

Definition at line 92 of file SolidMechanicsPlasticDruckerPragerHyperbolic.C.

{
  return "DruckerPragerHyperbolic";
}

◆ 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.

{
  return 1;
}

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

◆ onlyB()

void SolidMechanicsPlasticDruckerPrager::onlyB	(	Real	intnl,
		int	fd,
		Real &	bbb
	)		const

inherited

Calculate bbb or bbb_flow.

Parameters

intnl	the internal parameter
fd	if fd==friction then bbb is calculated. if fd==dilation then bbb_flow is calculated
bbb	either bbb or bbb_flow, depending on fd

Definition at line 163 of file SolidMechanicsPlasticDruckerPrager.C.

{
  if (_zero_phi_hardening && (fd == friction))
  {
    bbb = _bbb;
    return;
  }
  if (_zero_psi_hardening && (fd == dilation))
  {
    bbb = _bbb_flow;
    return;
  }
  initializeB(intnl, fd, bbb);
}

Referenced by CappedDruckerPragerStressUpdate::computeAllQ(), consistentTangentOperator(), SolidMechanicsPlasticDruckerPrager::dyieldFunction_dstress(), SolidMechanicsPlasticDruckerPrager::flowPotential(), CappedDruckerPragerStressUpdate::initializeVars(), returnMap(), CappedDruckerPragerStressUpdate::setIntnlDerivatives(), and CappedDruckerPragerStressUpdate::setIntnlValues().

◆ returnMap()

bool SolidMechanicsPlasticDruckerPragerHyperbolic::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

overrideprotectedvirtual

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_stress	The trial stress
	intnl_old	Value of the internal parameter
	E_ijkl	Elasticity tensor
	ep_plastic_tolerance	Tolerance defined by the user for the plastic strain
[out]	returned_stress	In case (C): lies on the yield surface after returning and produces the correct plastic strain (normality condition). Otherwise: not defined
[out]	returned_intnl	In case (C): the value of the internal parameter after returning. Otherwise: not defined
[out]	dpm	In case (C): the plastic multipliers needed to bring about the return. Otherwise: not defined
[out]	delta_dp	In case (C): The change in plastic strain induced by the return process. Otherwise: not defined
[out]	yf	In case (C): the yield function at (returned_stress, returned_intnl). Otherwise: the yield function at (trial_stress, intnl_old)
[out]	trial_stress_inadmissible	Should 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 98 of file SolidMechanicsPlasticDruckerPragerHyperbolic.C.

{
  if (!(_use_custom_returnMap))
    return SolidMechanicsPlasticModel::returnMap(trial_stress,
                                                 intnl_old,
                                                 E_ijkl,
                                                 ep_plastic_tolerance,
                                                 returned_stress,
                                                 returned_intnl,
                                                 dpm,
                                                 delta_dp,
                                                 yf,
                                                 trial_stress_inadmissible);
 
  yf.resize(1);
 
  yf[0] = yieldFunction(trial_stress, intnl_old);
 
  if (yf[0] < _f_tol)
  {
    // the trial_stress is admissible
    trial_stress_inadmissible = false;
    return true;
  }
 
  trial_stress_inadmissible = true;
  const Real mu = E_ijkl(0, 1, 0, 1);
  const Real lambda = E_ijkl(0, 0, 0, 0) - 2.0 * mu;
  const Real bulky = 3.0 * lambda + 2.0 * mu;
  const Real Tr_trial = trial_stress.trace();
  const Real J2trial = trial_stress.secondInvariant();
 
  // Perform a Newton-Raphson to find dpm when
  // residual = (1 + dpm*mu/ll)sqrt(ll^2 - s^2) - sqrt(J2trial) = 0, with s=smoother
  // with ll = sqrt(J2 + s^2) = aaa - bbb*Tr(stress) = aaa - bbb(Tr_trial - p*3*bulky*bbb_flow)
  Real aaa;
  Real daaa;
  Real bbb;
  Real dbbb;
  Real bbb_flow;
  Real dbbb_flow;
  Real ll;
  Real dll;
  Real residual;
  Real jac;
  dpm[0] = 0;
  unsigned int iter = 0;
  do
  {
    bothAB(intnl_old + dpm[0], aaa, bbb);
    dbothAB(intnl_old + dpm[0], daaa, dbbb);
    onlyB(intnl_old + dpm[0], dilation, bbb_flow);
    donlyB(intnl_old + dpm[0], dilation, dbbb_flow);
    ll = aaa - bbb * (Tr_trial - dpm[0] * bulky * 3 * bbb_flow);
    dll = daaa - dbbb * (Tr_trial - dpm[0] * bulky * 3 * bbb_flow) +
          bbb * bulky * 3 * (bbb_flow + dpm[0] * dbbb_flow);
    residual = bbb * (Tr_trial - dpm[0] * bulky * 3 * bbb_flow) - aaa +
               std::sqrt(J2trial / Utility::pow<2>(1 + dpm[0] * mu / ll) + _smoother2);
    jac = dbbb * (Tr_trial - dpm[0] * bulky * 3 * bbb_flow) -
          bbb * bulky * 3 * (bbb_flow + dpm[0] * dbbb_flow) - daaa +
          0.5 * J2trial * (-2.0) * (mu / ll - dpm[0] * mu * dll / ll / ll) /
              Utility::pow<3>(1 + dpm[0] * mu / ll) /
              std::sqrt(J2trial / Utility::pow<2>(1.0 + dpm[0] * mu / ll) + _smoother2);
    dpm[0] += -residual / jac;
    if (iter > _max_iters) // not converging
      return false;
    iter++;
  } while (residual * residual > _f_tol * _f_tol);
 
  // set the returned values
  yf[0] = 0;
  returned_intnl = intnl_old + dpm[0];
 
  bothAB(returned_intnl, aaa, bbb);
  onlyB(returned_intnl, dilation, bbb_flow);
  ll = aaa - bbb * (Tr_trial - dpm[0] * bulky * 3.0 * bbb_flow);
  returned_stress =
      trial_stress.deviatoric() / (1.0 + dpm[0] * mu / ll); // this is the deviatoric part only
 
  RankTwoTensor rij = 0.5 * returned_stress.deviatoric() /
                      ll; // this is the derivatoric part the flow potential only
 
  // form the returned stress and the full flow potential
  const Real returned_trace_over_3 = (aaa - ll) / bbb / 3.0;
  for (unsigned i = 0; i < 3; ++i)
  {
    returned_stress(i, i) += returned_trace_over_3;
    rij(i, i) += bbb_flow;
  }
 
  delta_dp = rij * dpm[0];
 
  return true;
}

◆ useCustomCTO()

bool SolidMechanicsPlasticDruckerPragerHyperbolic::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 272 of file SolidMechanicsPlasticDruckerPragerHyperbolic.C.

{
  return _use_custom_cto;
}

◆ useCustomReturnMap()

bool SolidMechanicsPlasticDruckerPragerHyperbolic::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 266 of file SolidMechanicsPlasticDruckerPragerHyperbolic.C.

{
  return _use_custom_returnMap;
}

◆ validParams()

InputParameters SolidMechanicsPlasticDruckerPragerHyperbolic::validParams ( )

static

Definition at line 21 of file SolidMechanicsPlasticDruckerPragerHyperbolic.C.

{
  InputParameters params = SolidMechanicsPlasticDruckerPrager::validParams();
  params.addParam<bool>("use_custom_returnMap",
                        true,
                        "Whether to use the custom returnMap "
                        "algorithm.  Set to true if you are using "
                        "isotropic elasticity.");
  params.addParam<bool>("use_custom_cto",
                        true,
                        "Whether to use the custom consistent tangent "
                        "operator computations.  Set to true if you are "
                        "using isotropic elasticity.");
  params.addClassDescription("J2 plasticity, associative, with hardening");
  params.addRangeCheckedParam<Real>("smoother",
                                    0.0,
                                    "smoother>=0",
                                    "The cone vertex at J2=0 is smoothed.  The maximum mean "
                                    "stress possible, which is Cohesion*Cot(friction_angle) for "
                                    "smoother=0, becomes (Cohesion - "
                                    "smoother/3)*Cot(friction_angle).  This is a non-hardening "
                                    "parameter");
  params.addRangeCheckedParam<unsigned>(
      "max_iterations",
      40,
      "max_iterations>0",
      "Maximum iterations to use in the custom return map function");
  params.addClassDescription(
      "Non-associative Drucker Prager plasticity with hyperbolic smoothing of the cone tip.");
  return params;
}

◆ yieldFunction()

Real SolidMechanicsPlasticDruckerPragerHyperbolic::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

stress	the stress at which to calculate the yield function
intnl	internal parameter

Returns: the yield function

Reimplemented from SolidMechanicsPlasticDruckerPrager.

Definition at line 64 of file SolidMechanicsPlasticDruckerPragerHyperbolic.C.

{
  Real aaa;
  Real bbb;
  bothAB(intnl, aaa, bbb);
  return std::sqrt(stress.secondInvariant() + _smoother2) + stress.trace() * bbb - aaa;
}

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

	stress	the stress at which to calculate the yield function
	intnl	internal parameter
[out]	f	the yield functions

Reimplemented in SolidMechanicsPlasticMohrCoulombMulti, and SolidMechanicsPlasticTensileMulti.

Definition at line 67 of file SolidMechanicsPlasticModel.C.

{
  f.assign(1, yieldFunction(stress, intnl));
}

Referenced by SolidMechanicsPlasticModel::returnMap().

Member Data Documentation

◆ _aaa

Real SolidMechanicsPlasticDruckerPrager::_aaa

privateinherited

Definition at line 108 of file SolidMechanicsPlasticDruckerPrager.h.

Referenced by SolidMechanicsPlasticDruckerPrager::bothAB(), and SolidMechanicsPlasticDruckerPrager::SolidMechanicsPlasticDruckerPrager().

◆ _bbb

Real SolidMechanicsPlasticDruckerPrager::_bbb

privateinherited

Definition at line 109 of file SolidMechanicsPlasticDruckerPrager.h.

Referenced by SolidMechanicsPlasticDruckerPrager::bothAB(), SolidMechanicsPlasticDruckerPrager::onlyB(), and SolidMechanicsPlasticDruckerPrager::SolidMechanicsPlasticDruckerPrager().

◆ _bbb_flow

Real SolidMechanicsPlasticDruckerPrager::_bbb_flow

privateinherited

Definition at line 110 of file SolidMechanicsPlasticDruckerPrager.h.

Referenced by SolidMechanicsPlasticDruckerPrager::onlyB(), and SolidMechanicsPlasticDruckerPrager::SolidMechanicsPlasticDruckerPrager().

◆ _f_tol

const Real SolidMechanicsPlasticModel::_f_tol

inherited

Tolerance on yield function.

Definition at line 170 of file SolidMechanicsPlasticModel.h.

◆ _ic_tol

const Real SolidMechanicsPlasticModel::_ic_tol

inherited

Tolerance on internal constraint.

Definition at line 173 of file SolidMechanicsPlasticModel.h.

◆ _max_iters

const unsigned SolidMechanicsPlasticDruckerPragerHyperbolic::_max_iters

private

max iters for custom return map loop

Definition at line 72 of file SolidMechanicsPlasticDruckerPragerHyperbolic.h.

Referenced by returnMap().

◆ _mc_cohesion

const SolidMechanicsHardeningModel& SolidMechanicsPlasticDruckerPrager::_mc_cohesion

protectedinherited

Hardening model for cohesion.

Definition at line 65 of file SolidMechanicsPlasticDruckerPrager.h.

Referenced by SolidMechanicsPlasticDruckerPrager::dbothAB(), SolidMechanicsPlasticDruckerPrager::initializeAandB(), and SolidMechanicsPlasticDruckerPrager::SolidMechanicsPlasticDruckerPrager().

◆ _mc_interpolation_scheme

const MooseEnum SolidMechanicsPlasticDruckerPrager::_mc_interpolation_scheme

protectedinherited

The parameters aaa and bbb are chosen to closely match the Mohr-Coulomb yield surface.

Various matching schemes may be used and this parameter holds the user's choice.

Definition at line 93 of file SolidMechanicsPlasticDruckerPrager.h.

Referenced by SolidMechanicsPlasticDruckerPrager::dbothAB(), SolidMechanicsPlasticDruckerPrager::donlyB(), SolidMechanicsPlasticDruckerPrager::initializeAandB(), and SolidMechanicsPlasticDruckerPrager::initializeB().

◆ _mc_phi

const SolidMechanicsHardeningModel& SolidMechanicsPlasticDruckerPrager::_mc_phi

protectedinherited

Hardening model for tan(phi)

Definition at line 68 of file SolidMechanicsPlasticDruckerPrager.h.

Referenced by SolidMechanicsPlasticDruckerPrager::dbothAB(), SolidMechanicsPlasticDruckerPrager::donlyB(), SolidMechanicsPlasticDruckerPrager::initializeAandB(), SolidMechanicsPlasticDruckerPrager::initializeB(), and SolidMechanicsPlasticDruckerPrager::SolidMechanicsPlasticDruckerPrager().

◆ _mc_psi

const SolidMechanicsHardeningModel& SolidMechanicsPlasticDruckerPrager::_mc_psi

protectedinherited

Hardening model for tan(psi)

Definition at line 71 of file SolidMechanicsPlasticDruckerPrager.h.

Referenced by SolidMechanicsPlasticDruckerPrager::donlyB(), SolidMechanicsPlasticDruckerPrager::initializeB(), and SolidMechanicsPlasticDruckerPrager::SolidMechanicsPlasticDruckerPrager().

◆ _smoother2

const Real SolidMechanicsPlasticDruckerPragerHyperbolic::_smoother2

private

smoothing parameter for the cone's tip

Definition at line 63 of file SolidMechanicsPlasticDruckerPragerHyperbolic.h.

Referenced by consistentTangentOperator(), df_dsig(), dflowPotential_dstress(), returnMap(), and yieldFunction().

◆ _use_custom_cto

const bool SolidMechanicsPlasticDruckerPragerHyperbolic::_use_custom_cto

private

Whether to use the custom consistent tangent operator calculation.

Definition at line 69 of file SolidMechanicsPlasticDruckerPragerHyperbolic.h.

Referenced by consistentTangentOperator(), and useCustomCTO().

◆ _use_custom_returnMap

const bool SolidMechanicsPlasticDruckerPragerHyperbolic::_use_custom_returnMap

private

whether to use the custom returnMap function

Definition at line 66 of file SolidMechanicsPlasticDruckerPragerHyperbolic.h.

Referenced by returnMap(), and useCustomReturnMap().

◆ _zero_cohesion_hardening

const bool SolidMechanicsPlasticDruckerPrager::_zero_cohesion_hardening

protectedinherited

True if there is no hardening of cohesion.

Definition at line 96 of file SolidMechanicsPlasticDruckerPrager.h.

Referenced by SolidMechanicsPlasticDruckerPrager::bothAB(), and SolidMechanicsPlasticDruckerPrager::dbothAB().

◆ _zero_phi_hardening

const bool SolidMechanicsPlasticDruckerPrager::_zero_phi_hardening

protectedinherited

True if there is no hardening of friction angle.

Definition at line 99 of file SolidMechanicsPlasticDruckerPrager.h.

Referenced by SolidMechanicsPlasticDruckerPrager::bothAB(), SolidMechanicsPlasticDruckerPrager::dbothAB(), SolidMechanicsPlasticDruckerPrager::donlyB(), and SolidMechanicsPlasticDruckerPrager::onlyB().

◆ _zero_psi_hardening

const bool SolidMechanicsPlasticDruckerPrager::_zero_psi_hardening

protectedinherited

True if there is no hardening of dilation angle.

Definition at line 102 of file SolidMechanicsPlasticDruckerPrager.h.

Referenced by SolidMechanicsPlasticDruckerPrager::donlyB(), and SolidMechanicsPlasticDruckerPrager::onlyB().

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

solid_mechanics/include/userobjects/SolidMechanicsPlasticDruckerPragerHyperbolic.h
solid_mechanics/src/userobjects/SolidMechanicsPlasticDruckerPragerHyperbolic.C

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

Detailed Description

Member Enumeration Documentation

◆ FrictionDilation

Constructor & Destructor Documentation

◆ SolidMechanicsPlasticDruckerPragerHyperbolic()

Member Function Documentation

◆ activeConstraints()

◆ bothAB()

◆ consistentTangentOperator()

◆ dbothAB()

◆ df_dsig()

◆ dflowPotential_dintnl()

◆ dflowPotential_dintnlV()

◆ dflowPotential_dstress()

◆ dflowPotential_dstressV()

◆ dhardPotential_dintnl()

◆ dhardPotential_dintnlV()

◆ dhardPotential_dstress()

◆ dhardPotential_dstressV()

◆ donlyB()

◆ dyieldFunction_dintnl()

◆ dyieldFunction_dintnlV()

◆ dyieldFunction_dstress()

◆ dyieldFunction_dstressV()

◆ execute()

◆ finalize()

◆ flowPotential()

◆ flowPotentialV()

◆ hardPotential()

◆ hardPotentialV()

◆ initialize()

◆ initializeAandB()

◆ initializeB()

◆ KuhnTuckerSingleSurface()

◆ modelName()

◆ numberSurfaces()

◆ onlyB()

◆ returnMap()

◆ useCustomCTO()

◆ useCustomReturnMap()

◆ validParams()

◆ yieldFunction()

◆ yieldFunctionV()

Member Data Documentation

◆ _aaa

◆ _bbb

◆ _bbb_flow

◆ _f_tol

◆ _ic_tol

◆ _max_iters

◆ _mc_cohesion

◆ _mc_interpolation_scheme

◆ _mc_phi

◆ _mc_psi

◆ _smoother2

◆ _use_custom_cto

◆ _use_custom_returnMap

◆ _zero_cohesion_hardening

◆ _zero_phi_hardening

◆ _zero_psi_hardening