This class defines a yield surface according to Von-Mises theory. More...

#include <mohr_coulomb_yield_surface.h>

Collaboration diagram for Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >:

Public Member Functions
Life Cycle
	MohrCoulombYieldSurface ()
	Initialization constructor. More...

	MohrCoulombYieldSurface (MohrCoulombYieldSurface const &rOther)
	Copy constructor. More...

MohrCoulombYieldSurface &	operator= (MohrCoulombYieldSurface const &rOther)
	Assignment operator. More...

virtual	~MohrCoulombYieldSurface ()
	Destructor. More...

Static Public Member Functions
Operations
static void	CalculateEquivalentStress (const array_1d< double, VoigtSize > &rPredictiveStressVector, const Vector &rStrainVector, double &rEquivalentStress, ConstitutiveLaw::Parameters &rValues)
	This method the uniaxial equivalent stress. More...

static void	GetInitialUniaxialThreshold (ConstitutiveLaw::Parameters &rValues, double &rThreshold)
	This method returns the initial uniaxial stress threshold. More...

static void	CalculateDamageParameter (ConstitutiveLaw::Parameters &rValues, double &rAParameter, const double CharacteristicLength)
	This method returns the damage parameter needed in the exp/linear expressions of damage. More...

static void	CalculatePlasticPotentialDerivative (const array_1d< double, VoigtSize > &rPredictiveStressVector, const array_1d< double, VoigtSize > &rDeviator, const double J2, array_1d< double, VoigtSize > &rDerivativePlasticPotential, ConstitutiveLaw::Parameters &rValues)
	This method calculates the derivative of the plastic potential DG/DS. More...

static void	CalculateYieldSurfaceDerivative (const array_1d< double, VoigtSize > &rPredictiveStressVector, const array_1d< double, VoigtSize > &rDeviator, const double J2, array_1d< double, VoigtSize > &rFFlux, ConstitutiveLaw::Parameters &rValues)
	This script calculates the derivatives of the Yield Surf according to NAYAK-ZIENKIEWICZ paper International journal for numerical methods in engineering vol 113-135 1972. As: DF/DS = c1V1 + c2V2 + c3*V3. More...

static int	Check (const Properties &rMaterialProperties)
	This method defines the check to be performed in the yield surface. More...

static bool	IsWorkingWithTensionThreshold ()
	This method returns true if the yield surfacecompares with the tension tield stress. More...

static double	GetScaleFactorTension (const Properties &rMaterialProperties)
	This method returns the scaling factor of the yield surface surfacecompares with the tension tield stress. More...

Type Definitions
typedef TPlasticPotentialType	PlasticPotentialType
	The type of potential plasticity. More...

static constexpr SizeType	Dimension = PlasticPotentialType::Dimension
	The Plastic potential already defines the working simension size. More...

static constexpr SizeType	VoigtSize = PlasticPotentialType::VoigtSize
	The Plastic potential already defines the Voigt size. More...

static constexpr double	tolerance = std::numeric_limits<double>::epsilon()
	The machine precision zero tolerance. More...

	KRATOS_CLASS_POINTER_DEFINITION (MohrCoulombYieldSurface)
	Counted pointer of MohrCoulombYieldSurface. More...

Detailed Description

template<class TPlasticPotentialType>
class Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >

This class defines a yield surface according to Von-Mises theory.

The Mohr–Coulomb failure surface is a cone with a hexagonal cross section in deviatoric stress space The yield surface requires the definition of the following properties:

FRACTURE_ENERGY: A fracture energy-based function is used to describe strength degradation in post-peak regime
YOUNG_MODULUS: It defines the relationship between stress (force per unit area) and strain (proportional deformation) in a material in the linear elasticity regime of a uniaxial deformation.
YIELD_STRESS: Yield stress is the amount of stress that an object needs to experience for it to be permanently deformed. Does not require to be defined simmetrically, one YIELD_STRESS_COMPRESSION and other YIELD_STRESS_TENSION can be defined for not symmetric cases
COHESION: Is the intercept of the failure envelope with the tau axis
See also
https://en.wikipedia.org/wiki/Mohr%E2%80%93Coulomb_theory

Template Parameters

TPlasticPotentialType The plastic potential considered

Author
Alejandro Cornejo & Lucia Barbu

Member Typedef Documentation

◆ PlasticPotentialType

template<class TPlasticPotentialType >

typedef TPlasticPotentialType Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >::PlasticPotentialType

The type of potential plasticity.

Constructor & Destructor Documentation

◆ MohrCoulombYieldSurface() [1/2]

template<class TPlasticPotentialType >

Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >::MohrCoulombYieldSurface ( )

inline

Initialization constructor.

◆ MohrCoulombYieldSurface() [2/2]

template<class TPlasticPotentialType >

Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >::MohrCoulombYieldSurface ( MohrCoulombYieldSurface< TPlasticPotentialType > const & rOther )

inline

Copy constructor.

◆ ~MohrCoulombYieldSurface()

template<class TPlasticPotentialType >

virtual Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >::~MohrCoulombYieldSurface ( )

inlinevirtual

Destructor.

Member Function Documentation

◆ CalculateDamageParameter()

template<class TPlasticPotentialType >

static void Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >::CalculateDamageParameter	(	ConstitutiveLaw::Parameters &	rValues,
		double &	rAParameter,
		const double	CharacteristicLength
	)

inlinestatic

This method returns the damage parameter needed in the exp/linear expressions of damage.

Parameters

rAParameter	The damage parameter
rValues	Parameters of the constitutive law
CharacteristicLength	The equivalent length of the FE

◆ CalculateEquivalentStress()

template<class TPlasticPotentialType >

static void Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >::CalculateEquivalentStress	(	const array_1d< double, VoigtSize > &	rPredictiveStressVector,
		const Vector &	rStrainVector,
		double &	rEquivalentStress,
		ConstitutiveLaw::Parameters &	rValues
	)

inlinestatic

This method the uniaxial equivalent stress.

Parameters

rPredictiveStressVector	The predictive stress vector S = C:(E-Ep)
rStrainVector	The StrainVector vector
rValues	Parameters of the constitutive law

◆ CalculatePlasticPotentialDerivative()

template<class TPlasticPotentialType >

static void Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >::CalculatePlasticPotentialDerivative	(	const array_1d< double, VoigtSize > &	rPredictiveStressVector,
		const array_1d< double, VoigtSize > &	rDeviator,
		const double	J2,
		array_1d< double, VoigtSize > &	rDerivativePlasticPotential,
		ConstitutiveLaw::Parameters &	rValues
	)

inlinestatic

This method calculates the derivative of the plastic potential DG/DS.

Parameters

StressVector	The stress vector
Deviator	The deviatoric part of the stress vector
J2	The second invariant of the Deviator
rDerivativePlasticPotential	The derivative of the plastic potential
rValues	Parameters of the constitutive law

◆ CalculateYieldSurfaceDerivative()

template<class TPlasticPotentialType >

static void Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >::CalculateYieldSurfaceDerivative	(	const array_1d< double, VoigtSize > &	rPredictiveStressVector,
		const array_1d< double, VoigtSize > &	rDeviator,
		const double	J2,
		array_1d< double, VoigtSize > &	rFFlux,
		ConstitutiveLaw::Parameters &	rValues
	)

inlinestatic

This script calculates the derivatives of the Yield Surf according to NAYAK-ZIENKIEWICZ paper International journal for numerical methods in engineering vol 113-135 1972. As: DF/DS = c1*V1 + c2*V2 + c3*V3.

Parameters

rPredictiveStressVector	The predictive stress vector S = C:(E-Ep)
rDeviator	The deviatoric part of the stress vector
J2	The second invariant of the Deviator
rFFlux	The derivative of the yield surface
rValues	Parameters of the constitutive law

◆ Check()

template<class TPlasticPotentialType >

static int Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >::Check ( const Properties & rMaterialProperties )

inlinestatic

This method defines the check to be performed in the yield surface.

Returns: 0 if OK, 1 otherwise

◆ GetInitialUniaxialThreshold()

template<class TPlasticPotentialType >

static void Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >::GetInitialUniaxialThreshold	(	ConstitutiveLaw::Parameters &	rValues,
		double &	rThreshold
	)

inlinestatic

This method returns the initial uniaxial stress threshold.

Parameters

rThreshold	The uniaxial stress threshold
rValues	Parameters of the constitutive law

◆ GetScaleFactorTension()

template<class TPlasticPotentialType >

static double Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >::GetScaleFactorTension ( const Properties & rMaterialProperties )

inlinestatic

This method returns the scaling factor of the yield surface surfacecompares with the tension tield stress.

◆ IsWorkingWithTensionThreshold()

template<class TPlasticPotentialType >

static bool Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >::IsWorkingWithTensionThreshold ( )

inlinestatic

This method returns true if the yield surfacecompares with the tension tield stress.

◆ KRATOS_CLASS_POINTER_DEFINITION()

template<class TPlasticPotentialType >

Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >::KRATOS_CLASS_POINTER_DEFINITION ( MohrCoulombYieldSurface< TPlasticPotentialType > )

Counted pointer of MohrCoulombYieldSurface.

◆ operator=()

template<class TPlasticPotentialType >

MohrCoulombYieldSurface& Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >::operator= ( MohrCoulombYieldSurface< TPlasticPotentialType > const & rOther )

inline

Assignment operator.

Member Data Documentation

◆ Dimension

template<class TPlasticPotentialType >

constexpr SizeType Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >::Dimension = PlasticPotentialType::Dimension

staticconstexpr

The Plastic potential already defines the working simension size.

◆ tolerance

template<class TPlasticPotentialType >

constexpr double Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >::tolerance = std::numeric_limits<double>::epsilon()

staticconstexpr

The machine precision zero tolerance.

◆ VoigtSize

template<class TPlasticPotentialType >

constexpr SizeType Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >::VoigtSize = PlasticPotentialType::VoigtSize

staticconstexpr

The Plastic potential already defines the Voigt size.

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

/home/runner/work/Documentation/Documentation/master/applications/ConstitutiveLawsApplication/custom_constitutive/auxiliary_files/yield_surfaces/mohr_coulomb_yield_surface.h

Public Member Functions

Static Public Member Functions

Type Definitions

Detailed Description

template<class TPlasticPotentialType> class Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >

Member Typedef Documentation

◆ PlasticPotentialType

Constructor & Destructor Documentation

◆ MohrCoulombYieldSurface() [1/2]

◆ MohrCoulombYieldSurface() [2/2]

◆ ~MohrCoulombYieldSurface()

Member Function Documentation

◆ CalculateDamageParameter()

◆ CalculateEquivalentStress()

◆ CalculatePlasticPotentialDerivative()

◆ CalculateYieldSurfaceDerivative()

◆ Check()

◆ GetInitialUniaxialThreshold()

◆ GetScaleFactorTension()

◆ IsWorkingWithTensionThreshold()

◆ KRATOS_CLASS_POINTER_DEFINITION()

◆ operator=()

Member Data Documentation

◆ Dimension

◆ tolerance

◆ VoigtSize

template<class TPlasticPotentialType>
class Kratos::MohrCoulombYieldSurface< TPlasticPotentialType >