This class defines a yield surface according to Drucker-Prager theory. More...

#include <drucker_prager_yield_surface.h>

Collaboration diagram for Kratos::DruckerPragerYieldSurface< TPlasticPotentialType >:

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

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

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

virtual	~DruckerPragerYieldSurface ()
	Destructor. More...

Static Public Member Functions
Operations
static void	CalculateEquivalentStress (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 > &rGFlux, 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 (DruckerPragerYieldSurface)
	Counted pointer of DruckerPragerYieldSurface. More...

Detailed Description

template<class TPlasticPotentialType>
class Kratos::DruckerPragerYieldSurface< TPlasticPotentialType >

This class defines a yield surface according to Drucker-Prager theory.

The Drucker–Prager yield criterion is similar to the von Mises yield criterion, with provisions for handling materials with differing tensile and compressive yield strengths. This criterion is most often used for concrete where both normal and shear stresses can determine failure. 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
FRICTION_ANGLE: Its definition is derived from the Mohr-Coulomb failure criterion and it is used to describe the friction shear resistance of soils together with the normal effective stress.
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
See also
https://en.wikipedia.org/wiki/Drucker%E2%80%93Prager_yield_criterion

Template Parameters

TPlasticPotentialType The plastic potential considered

Author
Alejandro Cornejo & Lucia Barbu

Member Typedef Documentation

◆ PlasticPotentialType

template<class TPlasticPotentialType >

typedef TPlasticPotentialType Kratos::DruckerPragerYieldSurface< TPlasticPotentialType >::PlasticPotentialType

The type of potential plasticity.

Constructor & Destructor Documentation

◆ DruckerPragerYieldSurface() [1/2]

template<class TPlasticPotentialType >

Kratos::DruckerPragerYieldSurface< TPlasticPotentialType >::DruckerPragerYieldSurface ( )

inline

Initialization constructor.

◆ DruckerPragerYieldSurface() [2/2]

template<class TPlasticPotentialType >

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

inline

Copy constructor.

◆ ~DruckerPragerYieldSurface()

template<class TPlasticPotentialType >

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

inlinevirtual

Destructor.

Member Function Documentation

◆ CalculateDamageParameter()

template<class TPlasticPotentialType >

static void Kratos::DruckerPragerYieldSurface< 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::DruckerPragerYieldSurface< TPlasticPotentialType >::CalculateEquivalentStress	(	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
rEquivalentStress	The effective stress or equivalent uniaxial stress is a scalar. It is an invariant value which measures the “intensity” of a 3D stress state.

◆ CalculatePlasticPotentialDerivative()

template<class TPlasticPotentialType >

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

inlinestatic

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

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
rGFlux	The derivative of the plastic potential
rValues	Parameters of the constitutive law

◆ CalculateYieldSurfaceDerivative()

template<class TPlasticPotentialType >

static void Kratos::DruckerPragerYieldSurface< 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::DruckerPragerYieldSurface< 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::DruckerPragerYieldSurface< 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::DruckerPragerYieldSurface< 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::DruckerPragerYieldSurface< TPlasticPotentialType >::IsWorkingWithTensionThreshold ( )

inlinestatic

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

◆ KRATOS_CLASS_POINTER_DEFINITION()

template<class TPlasticPotentialType >

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

Counted pointer of DruckerPragerYieldSurface.

◆ operator=()

template<class TPlasticPotentialType >

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

inline

Assignment operator.

Member Data Documentation

◆ Dimension

template<class TPlasticPotentialType >

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

staticconstexpr

The Plastic potential already defines the working simension size.

◆ tolerance

template<class TPlasticPotentialType >

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

staticconstexpr

The machine precision zero tolerance.

◆ VoigtSize

template<class TPlasticPotentialType >

constexpr SizeType Kratos::DruckerPragerYieldSurface< 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/drucker_prager_yield_surface.h

Public Member Functions

Static Public Member Functions

Type Definitions

Detailed Description

template<class TPlasticPotentialType> class Kratos::DruckerPragerYieldSurface< TPlasticPotentialType >

Member Typedef Documentation

◆ PlasticPotentialType

Constructor & Destructor Documentation

◆ DruckerPragerYieldSurface() [1/2]

◆ DruckerPragerYieldSurface() [2/2]

◆ ~DruckerPragerYieldSurface()

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::DruckerPragerYieldSurface< TPlasticPotentialType >