drprag Class Reference

#include <drprag.h>

Public Member Functions
void	changeparam (atsel &atm, vector &val)
double	cohesion (vector &qtr)
void	der_q_gamma (vector &dqdg)
void	derpotsigma (matrix &sig, matrix &dgds, vector &q)
void	deryieldfq (vector &qtr, vector &dfq)
void	deryieldfsigma (matrix &sig, matrix &dfds, vector &q)
	drprag (void)
double	give_consparam (long ipp, long ido)
void	giveirrstrains (long ipp, long ido, vector &epsp)
void	matstiff (matrix &d, long ipp, long ido)
void	nlstresses (long ipp, long im, long ido)
void	nonloc_nlstresses (long ipp, long im, long ido)
double	plasmodscalar (vector &qtr)
void	read (XFILE *in)
void	tangentstiff (matrix &d, matrix &td, long ipp, long ido)
void	updateq (long ipp, vector &epsp, vector &q)
void	updateval (long ipp, long im, long ido)
double	yieldfunction (matrix &sig, vector &q)
	~drprag (void)
Public Attributes
double	alpha
	material constant alpha
double	alpha1
	material constant alpha1
double	beta
	material constant beta
double	c
	cohesion
double	clim
	limit cohesion
double	phi
	friction angle
double	psi
	dilatation
strretalg	sra
	stress return algorithm
double	theta
	angle of linear hardening/softening

Detailed Description

This class defines Drucker-Prager plasticity model

structure of stored data plastic strains, consistency parameter, hardening parametr

28.3.2002

Definition at line 20 of file drprag.h.

Constructor & Destructor Documentation

drprag ( void )

This constructor inializes attributes to zero values.

Definition at line 23 of file drprag.cpp.

References c, clim, phi, psi, and theta.

~drprag ( void )

This destructor is only for the formal purposes.

Definition at line 34 of file drprag.cpp.

Member Function Documentation

void changeparam	(	atsel &	atm,
		vector &	val
	)

This function changes required material parameters and updates material constants alpha, alpha1 and beta. It is used for stochastic caluclations.

Parameters:

	atm	- structure with description of selected parameters
	val	- vector of new values for the required material parameters.

Return values:

The

function returns value of consistency parameter.

4.4.2002 by Tomas Koudelka

Definition at line 613 of file drprag.cpp.

References alpha, alpha1, atsel::atrib, beta, c, clim, atsel::num, phi, psi, and theta.

double cohesion ( vector & qtr )

This function computes value of cohesion which depends on the hardening parameter.

Parameters:

qtr

- vector of hardening parameters

Returns:: This function returns value of the cohesion coresponding to the given hardening parameter.

4.4.2002 by Tomas Koudelka

Definition at line 209 of file drprag.cpp.

References c, clim, and theta.

Referenced by deryieldfsigma(), and yieldfunction().

void der_q_gamma ( vector & dqdg )

This function computes derivatives of hardening paramters with respect of consistency parameter gamma.

Parameters:

dqdg

- matrix where the resulting derivatives are stored

Returns:: The function returns resulting vector of derivatives in the parameter dqdg

4.4.2002 by Tomas Koudelka

Definition at line 258 of file drprag.cpp.

References alpha1.

Referenced by plasmodscalar().

void derpotsigma	(	matrix &	sig,
		matrix &	dgds,
		vector &	q
	)

This function computes derivatives of plastic potential function with respect of vector sigma.

Parameters:

	sig	- stress tensor
	dgds	- matrix where the resulting derivatives are stored
	q	- vector of hardening parameters

Returns:: The function returns resulting matrix of derivatives in the parameter dgds.

4.4.2002 by Tomas Koudelka

Definition at line 164 of file drprag.cpp.

References alpha1, c, cmulm(), deviator(), fillm(), first_invar(), h(), psi, and tensornorm().

void deryieldfq	(	vector &	qtr,
		vector &	dfq
	)

This function computes derivatives of as-th yield function with respect of vector of hradening parameters.

Parameters:

	qtr	- vector of the hardening parameters
	dfds	- matrix where the resulting derivatives are stored

Returns:: The function returns resulting vector of derivatives in the parameter dfq.

4.4.2002 by Tomas Koudelka

Definition at line 233 of file drprag.cpp.

References beta, c, clim, and theta.

Referenced by plasmodscalar().

void deryieldfsigma	(	matrix &	sig,
		matrix &	dfds,
		vector &	q
	)

This function computes derivatives of as-th yield function with respect of vector sigma.

Parameters:

	sig	- stress tensor
	dfds	- matrix where the resulting derivatives are stored
	q	- vector of hardening parameters

Returns:: The function returns resulting matrix of derivatives in the parameter dfds.

4.4.2002 by Tomas Koudelka

Definition at line 118 of file drprag.cpp.

References alpha, cmulm(), cohesion(), deviator(), fillm(), first_invar(), h(), phi, and tensornorm().

Referenced by mechmat::dfdsigma(), mechmat::dgdsigma(), and tangentstiff().

double give_consparam	(	long	ipp,
		long	ido
	)

This function extracts consistency parametr gamma for the attained equilibrium state from the integration point eqother array.

Parameters:

	ipp	- integration point number in the mechmat ip array.
	ido	- index of internal variables for given material in the ipp other array

Return values:

The

function returns value of consistency parameter.

4.4.2002 by Tomas Koudelka

Definition at line 588 of file drprag.cpp.

References intpoints::eqother, mechmat::ip, Mm, and intpoints::ncompstr.

Referenced by mechmat::give_consparam().

void giveirrstrains	(	long	ipp,
		long	ido,
		vector &	epsp
	)

Function returns irreversible plastic strains.

Parameters:

	ipp	- integration point number in the mechmat ip array.
	ido	- index of the first internal variable for given material in the ipp other array
	epsp	- vector of irreversible strains

Returns:: The function returns the vector of irreversible strains in the parameter epsp

4.4.2002 by Tomas Koudelka

Definition at line 568 of file drprag.cpp.

References intpoints::eqother, mechmat::ip, Mm, and vector::n.

Referenced by mechmat::giveirrstrains().

void matstiff	(	matrix &	d,
		long	ipp,
		long	ido
	)

This function computes material stiffnes matrix.

Parameters:

	d	- allocated matrix structure for material stiffness matrix
	ipp	- integration point number
	ido	- index of internal variables for given material in the ipp other array

Returns:: The function returns calculated stiffness matrix in the parameter d.

4.4.2002 by Tomas Koudelka

Definition at line 324 of file drprag.cpp.

References mechmat::elmatstiff(), initial_stiff, matrix::m, Mm, Mp, matrix::n, probdesc::nlman, nonlinman::stmat, tangent_stiff, and tangentstiff().

Referenced by mechmat::matstiff().

void nlstresses	(	long	ipp,
		long	im,
		long	ido
	)

This function computes stresses at given integration point ipp, depending on the reached strains. The cutting plane algorithm is used. The stress and the other attribute of given integration point is actualized.

Parameters:

	ipp	- integration point number in the mechmat ip array.
	im	- index of material type for given ip
	ido	- index of internal variables for given material in the ipp other array

Returns:: The function updates the stress array of the given integration point.

4.4.2002 by Tomas Koudelka

Definition at line 421 of file drprag.cpp.

References cp, mechmat::cutting_plane(), deviator(), intpoints::eqother, first_invar(), strretalg::give_err(), strretalg::give_ni(), strretalg::give_tsra(), mechmat::ip, Mm, intpoints::ncompstr, intpoints::other, print_err(), second_invar(), sra, intpoints::ssst, intpoints::strain, stress, intpoints::stress, and vector_tensor().

Referenced by mechmat::compnonloc_nlstresses(), and mechmat::computenlstresses().

void nonloc_nlstresses	(	long	ipp,
		long	im,
		long	ido
	)

This function computes stresses at given integration point ipp, depending on the reached averaged nonlocal strains. The cutting plane algorithm is used. The stress and the other attribute of given integration point are actualized.

Parameters:

	ipp	- integration point number in the mechmat ip array.
	im	- index of material type for given ip
	ido	- index of internal variables for given material in the ipp other array

Returns:: The function updates the stress array of the given integration point.

4.4.2002 by Tomas Koudelka

Definition at line 485 of file drprag.cpp.

References cp, mechmat::cutting_plane(), deviator(), intpoints::eqother, first_invar(), strretalg::give_err(), strretalg::give_ni(), strretalg::give_tsra(), mechmat::ip, Mm, intpoints::ncompstr, intpoints::nonloc, intpoints::other, print_err(), second_invar(), sra, intpoints::ssst, intpoints::strain, stress, intpoints::stress, and vector_tensor().

Referenced by mechmat::compnonloc_nlstresses().

double plasmodscalar ( vector & qtr )

This function computes plastic modulus.

Parameters:

qtr

-vector of hardening parameters

Returns:: The function returns value of the plastic modulus.

4.4.2002 by Tomas Koudelka

Definition at line 274 of file drprag.cpp.

References der_q_gamma(), deryieldfq(), vector::n, and scprd().

Referenced by mechmat::plasmodscalar(), and tangentstiff().

void read ( XFILE * in )

This function reads material parameters from the opened text file given by the parameter in. Then it computes material constants alpha, alpha1 and beta

Parameters:

in

- pointer to the opned text file

25.3.2002 by Tomas Koudelka

Definition at line 50 of file drprag.cpp.

References alpha, alpha1, beta, c, clim, phi, psi, strretalg::read(), sra, theta, and xfscanf().

Referenced by mechmat::readmattype().

void tangentstiff	(	matrix &	d,
		matrix &	td,
		long	ipp,
		long	ido
	)

This function computes elastic-plastic material stiffnes matrix.

Parameters:

	d	- elastic stiffness matrix
	td	- allocated matrix structure for material tangent stiffness matrix
	ipp	- integration point number
	ido	- index of internal variables for given material in the ipp other array

Returns:: The function returns calculated tangent stiffness matrix in the parameter td.

4.4.2002 by Tomas Koudelka

Definition at line 354 of file drprag.cpp.

References allocv(), cmulm(), copym(), deryieldfsigma(), destrv(), intpoints::eqother, mechmat::givestress(), mechmat::ip, matrix::m, Mm, mxm(), mxv(), matrix::n, intpoints::ncompstr, planestrain, planestress, plasmodscalar(), scprd(), intpoints::ssst, stress, subm(), tensor_vector(), vector_tensor(), and vxv().

Referenced by matstiff().

void updateq	(	long	ipp,
		vector &	epsp,
		vector &	q
	)

This function computes new value of teh hardening parameter q.

Parameters:

	ipp	- integration point pointer
	epsp	- vector of the reached plastic strains
	q	- vector of the hardening parameters

Returns:: The function updates components of the vector of hardening parameters q.

4.4.2002 by Tomas Koudelka

Definition at line 302 of file drprag.cpp.

References cumulstrain(), mechmat::ip, Mm, intpoints::ssst, strain, and vector_tensor().

Referenced by mechmat::updateq().

void updateval	(	long	ipp,
		long	im,
		long	ido
	)

This function updates values in the other array reached in the previous equlibrium state to values reached in the new actual equilibrium state.

Parameters:

	ipp	- integration point number in the mechmat ip array.
	im	- index of material type for given ip
	ido	- index of internal variables for given material in the ipp other array