#include <abaqususermaterial.h>

Inheritance diagram for oofem::AbaqusUserMaterial:

Collaboration diagram for oofem::AbaqusUserMaterial:

Public Member Functions
	AbaqusUserMaterial (int n, Domain *d)
	Constructor.
virtual	~AbaqusUserMaterial ()
	Destructor.
void	initializeFrom (InputRecord &ir) override
void	giveInputRecord (DynamicInputRecord &input) override
std::unique_ptr< MaterialStatus >	CreateStatus (GaussPoint *gp) const override
FloatMatrixF< 6, 6 >	give3dMaterialStiffnessMatrix (MatResponseMode mode, GaussPoint gp, TimeStep tStep) const override
FloatMatrixF< 9, 9 >	give3dMaterialStiffnessMatrix_dPdF (MatResponseMode mode, GaussPoint gp, TimeStep tStep) const override
FloatMatrixF< 5, 5 >	givePlaneStrainStiffnessMatrix_dPdF (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const override
FloatArrayF< 6 >	giveRealStressVector_3d (const FloatArrayF< 6 > &strain, GaussPoint gp, TimeStep tStep) const override
	Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.
FloatArrayF< 9 >	giveFirstPKStressVector_3d (const FloatArrayF< 9 > &vF, GaussPoint gp, TimeStep tStep) const override
	Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.
int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep) override
const char *	giveClassName () const override
const char *	giveInputRecordName () const override
Public Member Functions inherited from oofem::StructuralMaterial
	StructuralMaterial (int n, Domain *d)
bool	hasMaterialModeCapability (MaterialMode mode) const override
const char *	giveClassName () const override
void	initializeFrom (InputRecord &ir) override
void	giveInputRecord (DynamicInputRecord &input) override
void	giveCharacteristicMatrix (FloatMatrix &answer, MatResponseMode type, GaussPoint gp, TimeStep tStep) const override
	Returns characteristic matrix of the receiver.
void	giveCharacteristicVector (FloatArray &answer, FloatArray &flux, MatResponseMode type, GaussPoint gp, TimeStep tStep) const override
	Returns characteristic vector of the receiver.
virtual void	giveStiffnessMatrix (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep) const
virtual void	giveRealStressVector (FloatArray &answer, GaussPoint gp, const FloatArray &reducedStrain, TimeStep tStep) const
virtual FloatArrayF< 4 >	giveRealStressVector_PlaneStrain (const FloatArrayF< 4 > &strain, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_3d.
virtual FloatArray	giveRealStressVector_StressControl (const FloatArray &reducedE, const IntArray &strainControl, GaussPoint gp, TimeStep tStep) const
	Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero·
virtual FloatArray	giveRealStressVector_ShellStressControl (const FloatArray &reducedE, const IntArray &strainControl, GaussPoint gp, TimeStep tStep) const
virtual FloatArrayF< 3 >	giveRealStressVector_PlaneStress (const FloatArrayF< 3 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 1 >	giveRealStressVector_1d (const FloatArrayF< 1 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 2 >	giveRealStressVector_Warping (const FloatArrayF< 2 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 2 >	giveRealStressVector_2dBeamLayer (const FloatArrayF< 2 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 5 >	giveRealStressVector_PlateLayer (const FloatArrayF< 5 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 3 >	giveRealStressVector_Fiber (const FloatArrayF< 3 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveRealStressVector_StressControl.
virtual FloatArrayF< 3 >	giveRealStressVector_2dPlateSubSoil (const FloatArrayF< 3 > &reducedE, GaussPoint gp, TimeStep tStep) const
	Default implementation is not provided.
virtual FloatArrayF< 6 >	giveRealStressVector_3dBeamSubSoil (const FloatArrayF< 6 > &reducedE, GaussPoint gp, TimeStep tStep) const
virtual FloatArrayF< 5 >	giveFirstPKStressVector_PlaneStrain (const FloatArrayF< 5 > &vF, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveFirstPKStressVector_3d.
virtual FloatArray	giveFirstPKStressVector_StressControl (const FloatArray &reducedvF, const IntArray &FControl, GaussPoint gp, TimeStep tStep) const
	Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero·
virtual FloatArrayF< 4 >	giveFirstPKStressVector_PlaneStress (const FloatArrayF< 4 > &vF, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveFirstPKStressVector_StressControl.
virtual FloatArrayF< 1 >	giveFirstPKStressVector_1d (const FloatArrayF< 1 > &vF, GaussPoint gp, TimeStep tStep) const
	Default implementation relies on giveFirstPKStressVector_StressControl.
virtual void	giveCauchyStressVector_3d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
virtual void	giveCauchyStressVector_PlaneStrain (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
virtual void	giveCauchyStressVector_PlaneStress (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
virtual void	giveCauchyStressVector_1d (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
virtual void	giveEshelbyStressVector_PlaneStrain (FloatArray &answer, GaussPoint gp, const FloatArray &reducedF, TimeStep tStep)
virtual FloatArrayF< 6 >	giveThermalDilatationVector (GaussPoint gp, TimeStep tStep) const
double	giveReferenceTemperature ()
virtual FloatArray	computeStressIndependentStrainVector (GaussPoint gp, TimeStep tStep, ValueModeType mode) const
FloatArrayF< 6 >	computeStressIndependentStrainVector_3d (GaussPoint gp, TimeStep tStep, ValueModeType mode) const
virtual void	giveStiffnessMatrix_dPdF (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
virtual void	give3dMaterialStiffnessMatrix_dCde (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
void	giveStressDependentPartOfStrainVector (FloatArray &answer, GaussPoint gp, const FloatArray &reducedStrainVector, TimeStep tStep, ValueModeType mode) const
int	setIPValue (const FloatArray &value, GaussPoint *gp, InternalStateType type) override
int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep) override
virtual FloatMatrixF< 3, 3 >	givePlaneStressStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 4, 4 >	givePlaneStressStiffnessMatrix_dPdF (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual void	givePlaneStressStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
virtual FloatMatrixF< 4, 4 >	givePlaneStrainStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual void	givePlaneStrainStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
virtual FloatMatrixF< 1, 1 >	give1dStressStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 1, 1 >	give1dStressStiffnessMatrix_dPdF (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual void	give1dStressStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint gp, TimeStep tStep)
virtual FloatMatrixF< 2, 2 >	give2dBeamLayerStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 5, 5 >	givePlateLayerStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 3, 3 >	giveFiberStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 3, 3 >	give2dPlateSubSoilStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
virtual FloatMatrixF< 6, 6 >	give3dBeamSubSoilStiffMtrx (MatResponseMode mmode, GaussPoint gp, TimeStep tStep) const
Public Member Functions inherited from oofem::Material
	Material (int n, Domain *d)
virtual	~Material ()=default
	Destructor.
virtual bool	isCharacteristicMtrxSymmetric (MatResponseMode rMode) const
virtual double	giveCharacteristicValue (MatResponseMode type, GaussPoint gp, TimeStep tStep) const
	Returns characteristic value of the receiver.
virtual double	give (int aProperty, GaussPoint *gp) const
virtual bool	hasProperty (int aProperty, GaussPoint *gp) const
virtual void	modifyProperty (int aProperty, double value, GaussPoint *gp)
double	giveCastingTime () const
virtual bool	isActivated (TimeStep *tStep) const
virtual bool	hasCastingTimeSupport () const
void	printYourself () override
	Prints receiver state on stdout. Useful for debugging.
virtual void	saveIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp)
virtual void	restoreIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp)
int	checkConsistency () override
virtual void	restoreConsistency (GaussPoint *gp)
virtual int	initMaterial (Element *element)
virtual MaterialStatus *	giveStatus (GaussPoint *gp) const
virtual int	packUnknowns (DataStream &buff, TimeStep tStep, GaussPoint ip)
virtual int	unpackAndUpdateUnknowns (DataStream &buff, TimeStep tStep, GaussPoint ip)
virtual int	estimatePackSize (DataStream &buff, GaussPoint *ip)
virtual double	predictRelativeComputationalCost (GaussPoint *gp)
virtual double	predictRelativeRedistributionCost (GaussPoint *gp)
virtual void	initTempStatus (GaussPoint *gp) const
void	saveContext (DataStream &stream, ContextMode mode) override
void	restoreContext (DataStream &stream, ContextMode mode) override
Public Member Functions inherited from oofem::FEMComponent
	FEMComponent (int n, Domain *d)
virtual	~FEMComponent ()=default
	Virtual destructor.
Domain *	giveDomain () const
virtual void	setDomain (Domain *d)
int	giveNumber () const
void	setNumber (int num)
virtual void	updateLocalNumbering (EntityRenumberingFunctor &f)
virtual void	initializeFrom (InputRecord &ir, int priority)
virtual void	initializeFinish ()
virtual void	postInitialize ()
	Performs post initialization steps. Called after all components are created and initialized.
virtual void	printOutputAt (FILE file, TimeStep tStep)
virtual Interface *	giveInterface (InterfaceType t)
std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros).

Private Attributes
void *	umatobj = nullptr
	Dynamically loaded umat.
void(*	umat )(double stress, double statev, double ddsdde, double sse, double spd, double scd, double rpl, double ddsddt, double drplde, double drpldt, double stran, double dstran, double time[2], double dtime, double temp, double dtemp, double predef[1], double dpred[1], char cmname[80], int ndi, int nshr, int ntens, int nstatv, double props, int nprops, double coords[3], double drot, double pnewdt, double celent, double dfgrd0, double dfgrd1, int noel, int npt, int layer, int kspt, int kstep, int kinc) = nullptr
	Pointer to the dynamically loaded umat-function (translated to C).
char	cmname [80]
	Name for material routine.
int	numState = 0
	Size of the state vector.
FloatArray	properties
	Material properties.
FloatArrayF< 6 >	initialStress
	Initial stress.
int	mStressInterpretation = 0
bool	mUseNumericalTangent = false
double	mPerturbation = 1e-7
std::string	filename
	Name of the file that contains the umat function.

Static Private Attributes
static std::size_t const	abq2oo9 [9] = { 0, 1, 2, 5, 4, 3, 6, 8, 7 }
static std::size_t const	abq2oo6 [6] = { 0, 1, 2, 5, 4, 3 }

Additional Inherited Members
Static Public Member Functions inherited from oofem::StructuralMaterial
static int	giveSymVI (int ind1, int ind2)
static int	giveVI (int ind1, int ind2)
static FloatMatrixF< 9, 9 >	convert_dSdE_2_dPdF_3D (const FloatMatrixF< 6, 6 > &dSdE, const FloatArrayF< 6 > &S, const FloatArrayF< 9 > &F)
static FloatMatrixF< 5, 5 >	convert_dSdE_2_dPdF_PlaneStrain (const FloatMatrixF< 4, 4 > &dSdE, const FloatArrayF< 4 > &S, const FloatArrayF< 5 > &F)
static FloatMatrixF< 4, 4 >	convert_dSdE_2_dPdF_PlaneStress (const FloatMatrixF< 3, 3 > &dSdE, const FloatArrayF< 3 > &S, const FloatArrayF< 4 > &F)
static FloatMatrixF< 1, 1 >	convert_dSdE_2_dPdF_1D (const FloatMatrixF< 1, 1 > &dSdE, const FloatArrayF< 1 > &S, const FloatArrayF< 1 > &F)
static void	computePrincipalValues (FloatArray &answer, const FloatArray &s, stressStrainPrincMode mode)
	Common functions for convenience.
static FloatArrayF< 3 >	computePrincipalValues (const FloatMatrixF< 3, 3 > &s)
static FloatArrayF< 3 >	computePrincipalValues (double I1, double I2, double I3)
static void	computePrincipalValDir (FloatArray &answer, FloatMatrix &dir, const FloatArray &s, stressStrainPrincMode mode)
static std::pair< FloatArrayF< 3 >, FloatMatrixF< 3, 3 > >	computePrincipalValDir (const FloatMatrixF< 3, 3 > &s)
static FloatArrayF< 6 >	computeDeviator (const FloatArrayF< 6 > &s)
static std::pair< FloatArrayF< 6 >, double >	computeDeviatoricVolumetricSplit (const FloatArrayF< 6 > &s)
static FloatArrayF< 6 >	computeDeviatoricVolumetricSum (const FloatArrayF< 6 > &dev, double mean)
static FloatArrayF< 6 >	applyDeviatoricElasticCompliance (const FloatArrayF< 6 > &stress, double EModulus, double nu)
static FloatArrayF< 6 >	applyDeviatoricElasticCompliance (const FloatArrayF< 6 > &stress, double GModulus)
static FloatArrayF< 6 >	applyDeviatoricElasticStiffness (const FloatArrayF< 6 > &strain, double EModulus, double nu)
static FloatArrayF< 6 >	applyDeviatoricElasticStiffness (const FloatArrayF< 6 > &strain, double GModulus)
static FloatArrayF< 6 >	applyElasticStiffness (const FloatArrayF< 6 > &strain, double EModulus, double nu)
static FloatArrayF< 6 >	applyElasticCompliance (const FloatArrayF< 6 > &stress, double EModulus, double nu)
static double	computeStressNorm (const FloatArrayF< 6 > &stress)
static double	computeFirstInvariant (const FloatArrayF< 6 > &s)
static double	computeSecondStressInvariant (const FloatArrayF< 6 > &s)
static double	computeThirdStressInvariant (const FloatArrayF< 6 > &s)
static double	computeFirstCoordinate (const FloatArrayF< 6 > &s)
static double	computeSecondCoordinate (const FloatArrayF< 6 > &s)
static double	computeThirdCoordinate (const FloatArrayF< 6 > &s)
static int	giveVoigtVectorMask (IntArray &answer, MaterialMode mmode)
static int	giveVoigtSymVectorMask (IntArray &answer, MaterialMode mmode)
static void	giveInvertedVoigtVectorMask (IntArray &answer, MaterialMode mmode)
static int	giveSizeOfVoigtVector (MaterialMode mmode)
static int	giveSizeOfVoigtSymVector (MaterialMode mmode)
static void	giveFullVectorForm (FloatArray &answer, const FloatArray &strainVector, MaterialMode matMode)
	Converts the reduced symmetric Voigt vector (2nd order tensor) to full form.
static void	giveFullVectorFormF (FloatArray &answer, const FloatArray &strainVector, MaterialMode matMode)
	Converts the reduced deformation gradient Voigt vector (2nd order tensor).
static void	giveFullSymVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode)
	Converts the reduced unsymmetric Voigt vector (2nd order tensor) to full form.
static void	giveReducedVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode)
	Converts the full symmetric Voigt vector (2nd order tensor) to reduced form.
static void	giveReducedSymVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode)
	Converts the full unsymmetric Voigt vector (2nd order tensor) to reduced form.
static void	giveFullSymMatrixForm (FloatMatrix &answer, const FloatMatrix &red, MaterialMode matMode)
	Converts the full unsymmetric Voigt matrix (4th order tensor) to reduced form.
static void	giveReducedMatrixForm (FloatMatrix &answer, const FloatMatrix &full, MaterialMode matMode)
	Converts the full symmetric Voigt matrix (4th order tensor) to reduced form.
static void	giveReducedSymMatrixForm (FloatMatrix &answer, const FloatMatrix &full, MaterialMode matMode)
	Converts the full unsymmetric Voigt matrix (4th order tensor) to reduced form.
static FloatArrayF< 6 >	transformStrainVectorTo (const FloatMatrixF< 3, 3 > &base, const FloatArrayF< 6 > &strain, bool transpose=false)
static FloatArrayF< 6 >	transformStressVectorTo (const FloatMatrixF< 3, 3 > &base, const FloatArrayF< 6 > &stress, bool transpose=false)
static double	computeVonMisesStress (const FloatArray &currentStress)
static double	computeVonMisesStress_3D (const FloatArrayF< 6 > &stress)
static double	computeVonMisesStress_PlaneStress (const FloatArrayF< 3 > &stress)
static FloatMatrixF< 6, 6 >	giveStrainVectorTranformationMtrx (const FloatMatrixF< 3, 3 > &base, bool transpose=false)
static FloatMatrixF< 3, 3 >	give2DStrainVectorTranformationMtrx (const FloatMatrixF< 2, 2 > &base, bool transpose=false)
static FloatMatrixF< 6, 6 >	giveStressVectorTranformationMtrx (const FloatMatrixF< 3, 3 > &base, bool transpose=false)
static FloatMatrixF< 3, 3 >	givePlaneStressVectorTranformationMtrx (const FloatMatrixF< 2, 2 > &base, bool transpose=false)
static void	sortPrincDirAndValCloseTo (FloatArray &pVal, FloatMatrix &pDir, const FloatMatrix &toPDir)
Static Public Attributes inherited from oofem::StructuralMaterial
static std::array< std::array< int, 3 >, 3 >	vIindex
	Voigt index map.
static std::array< std::array< int, 3 >, 3 >	svIndex
	Symmetric Voigt index map.
Protected Attributes inherited from oofem::StructuralMaterial
double	referenceTemperature = 0.
	Reference temperature (temperature, when material has been built into structure).
MatResponseMode	SCStiffMode = TangentStiffness
	stifness mode used in stress control
double	SCRelTol = 1.e-3
	relative tolerance for stress control
double	SCAbsTol = 1.e-12
	absolute stress tolerance for stress control
int	SCMaxiter = 100000
	maximum iterations for stress-control
Protected Attributes inherited from oofem::Material
Dictionary	propertyDictionary
double	castingTime
int	preCastingTimeMat
	Material existing before casting time - optional parameter, zero by default.
Protected Attributes inherited from oofem::FEMComponent
int	number
	Component number.
Domain *	domain
	Link to domain object, useful for communicating with other FEM components.

Detailed Description

This class allows for custom user materials from Abaqus (UMAT).

Note: Experimental, subject to change. Many optional arguments haven't been dealt with properly (and few umat functions actually use all of them).; Nothing for large deformations have been tested.

Author: Mikael Öhman

The umat material should be compiled as a shared library. For example, compile like this;

* gfortran -fPIC -fdefault-real-8 -shared -Wl,-soname,umat.so -o umat.so my_umat_code.f
*

The interface for a user material file for Abaqus is as follows: SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD, SCD,RPL,DDSDDT,DRPLDE,DRPLDT, STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME, NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT, CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)

CHARACTER*80 CMNAME DOUBLE PRECISION STRESS(NTENS),STATEV(NSTATV), DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS), STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1), PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3)

Definition at line 81 of file abaqususermaterial.h.

Constructor & Destructor Documentation

◆ AbaqusUserMaterial()

oofem::AbaqusUserMaterial::AbaqusUserMaterial	(	int	n,
		Domain *	d )

Constructor.

Definition at line 56 of file abaqususermaterial.C.

References oofem::StructuralMaterial::StructuralMaterial().

Referenced by give3dMaterialStiffnessMatrix(), and giveRealStressVector_3d().

◆ ~AbaqusUserMaterial()

oofem::AbaqusUserMaterial::~AbaqusUserMaterial ( )

virtual

Destructor.

Definition at line 60 of file abaqususermaterial.C.

References umatobj.

Member Function Documentation

◆ CreateStatus()

std::unique_ptr< MaterialStatus > oofem::AbaqusUserMaterial::CreateStatus ( GaussPoint * gp ) const

overridevirtual

Creates new copy of associated status and inserts it into given integration point.

Parameters

gp	Integration point where newly created status will be stored.

Returns: Reference to new status.

Reimplemented from oofem::Material.

Definition at line 142 of file abaqususermaterial.C.

References numState.

◆ give3dMaterialStiffnessMatrix()

FloatMatrixF< 6, 6 > oofem::AbaqusUserMaterial::give3dMaterialStiffnessMatrix	(	MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

Computes full 3d material stiffness matrix at given integration point, time, respecting load history in integration point.

Parameters

answer	Computed results.
mode	Material response mode.
gp	Integration point.
tStep	Time step (most models are able to respond only when tStep is current time step).

Todo: Make this hack fit more nicely into OOFEM in general;

Reimplemented from oofem::StructuralMaterial.

Definition at line 149 of file abaqususermaterial.C.

References AbaqusUserMaterial(), oofem::FloatArray::at(), oofem::GaussPoint::giveMaterialStatus(), giveRealStressVector_3d(), oofem::FloatArray::giveSize(), oofem::Material::giveStatus(), oofem::StructuralMaterialStatus::giveTempStrainVector(), oofem::StructuralMaterialStatus::giveTempStressVector(), oofem::AbaqusUserMaterialStatus::giveTempTangent(), oofem::FloatMatrix::printYourself(), oofem::FloatMatrix::setColumn(), oofem::FloatArray::subtract(), oofem::FloatArray::times(), and oofem::zeros().

◆ give3dMaterialStiffnessMatrix_dPdF()

FloatMatrixF< 9, 9 > oofem::AbaqusUserMaterial::give3dMaterialStiffnessMatrix_dPdF	(	MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

Todo: Make this hack fit more nicely into OOFEM in general;

Reimplemented from oofem::StructuralMaterial.

Definition at line 184 of file abaqususermaterial.C.

References giveFirstPKStressVector_3d(), oofem::Material::giveStatus(), mPerturbation, mUseNumericalTangent, and oofem::FloatMatrixF< N, M >::setColumn().

Referenced by givePlaneStrainStiffnessMatrix_dPdF().

◆ giveClassName()

const char * oofem::AbaqusUserMaterial::giveClassName ( ) const

inlineoverridevirtual

Returns: Class name of the receiver.

Implements oofem::FEMComponent.

Definition at line 161 of file abaqususermaterial.h.

◆ giveFirstPKStressVector_3d()

FloatArrayF< 9 > oofem::AbaqusUserMaterial::giveFirstPKStressVector_3d	(	const FloatArrayF< 9 > &	vF,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.

Todo: Check this. I'm just guessing. Maybe intrinsic time instead?

Todo: Right default value? umat routines may change this (although we ignore it)

Todo: Large deformations?

Todo: Include the characteristic element length

Todo: No idea about these parameters

Reimplemented from oofem::StructuralMaterial.

Definition at line 395 of file abaqususermaterial.C.

References abq2oo6, abq2oo9, oofem::Element::computeGlobalCoordinates(), oofem::det(), oofem::dot(), oofem::dotT(), E, oofem::eye(), oofem::from_voigt_form_9(), oofem::from_voigt_stress_6(), oofem::GaussPoint::giveElement(), oofem::TimeStep::giveMetaStepNumber(), oofem::GaussPoint::giveNaturalCoordinates(), oofem::FEMComponent::giveNumber(), oofem::FloatArray::givePointer(), oofem::FloatArrayF< N >::givePointer(), oofem::FloatMatrixF< N, M >::givePointer(), oofem::Material::giveStatus(), oofem::TimeStep::giveTargetTime(), oofem::TimeStep::giveTimeIncrement(), oofem::inv(), mStressInterpretation, numState, OOFEM_LOG_DEBUG, properties, S, oofem::Tdot(), oofem::to_voigt_form_33(), oofem::to_voigt_strain_33(), oofem::to_voigt_stress_33(), and umat.

Referenced by give3dMaterialStiffnessMatrix_dPdF().

◆ giveInputRecord()

void oofem::AbaqusUserMaterial::giveInputRecord ( DynamicInputRecord & input )

overridevirtual

Setups the input record string of receiver.

Parameters

input Dynamic input record to be filled by receiver.

Reimplemented from oofem::FEMComponent.

Definition at line 132 of file abaqususermaterial.C.

References _IFT_AbaqusUserMaterial_name, _IFT_AbaqusUserMaterial_numState, _IFT_AbaqusUserMaterial_properties, _IFT_AbaqusUserMaterial_userMaterial, cmname, filename, oofem::StructuralMaterial::giveInputRecord(), numState, properties, and oofem::DynamicInputRecord::setField().

◆ giveInputRecordName()

const char * oofem::AbaqusUserMaterial::giveInputRecordName ( ) const

inlineoverridevirtual

Returns: Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 162 of file abaqususermaterial.h.

References _IFT_AbaqusUserMaterial_Name.

◆ giveIPValue()

int oofem::AbaqusUserMaterial::giveIPValue	(	FloatArray &	answer,
		GaussPoint *	gp,
		InternalStateType	type,
		TimeStep *	tStep )

overridevirtual

Returns the integration point corresponding value in Reduced form.

Parameters

answer	Contain corresponding ip value, zero sized if not available.
gp	Integration point to which the value refers.
type	Determines the type of internal variable.
tStep	Determines the time step.

Returns: Nonzero if the assignment can be done, zero if this type of variable is not supported.

Reimplemented from oofem::Material.

Definition at line 578 of file abaqususermaterial.C.

References oofem::FloatArray::add(), oofem::StructuralMaterial::giveIPValue(), oofem::Material::giveStatus(), and initialStress.

◆ givePlaneStrainStiffnessMatrix_dPdF()

FloatMatrixF< 5, 5 > oofem::AbaqusUserMaterial::givePlaneStrainStiffnessMatrix_dPdF	(	MatResponseMode	mmode,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

Reimplemented from oofem::StructuralMaterial.

Definition at line 240 of file abaqususermaterial.C.

References give3dMaterialStiffnessMatrix_dPdF().

◆ giveRealStressVector_3d()

FloatArrayF< 6 > oofem::AbaqusUserMaterial::giveRealStressVector_3d	(	const FloatArrayF< 6 > &	strain,
		GaussPoint *	gp,
		TimeStep *	tStep ) const

overridevirtual

Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.

Todo: Check this. I'm just guessing. Maybe intrinsic time instead?

Todo: Right default value? umat routines may change this (although we ignore it)

Todo: Include the characteristic element length

Todo: No idea about these parameters

Reimplemented from oofem::StructuralMaterial.

Definition at line 248 of file abaqususermaterial.C.

References AbaqusUserMaterial(), abq2oo6, oofem::Element::computeGlobalCoordinates(), oofem::eye(), oofem::GaussPoint::giveElement(), oofem::GaussPoint::giveNaturalCoordinates(), oofem::FEMComponent::giveNumber(), oofem::FloatArray::givePointer(), oofem::FloatArrayF< N >::givePointer(), oofem::FloatMatrixF< N, M >::givePointer(), oofem::Material::giveStatus(), oofem::TimeStep::giveTargetTime(), oofem::TimeStep::giveTimeIncrement(), initialStress, numState, OOFEM_LOG_DEBUG, properties, and umat.

Referenced by give3dMaterialStiffnessMatrix().

◆ initializeFrom()

void oofem::AbaqusUserMaterial::initializeFrom ( InputRecord & ir )

overridevirtual

Reads the following values;

numstate (required, integer): number of state variables.
properties (required, FloatArray): material property values.
umat (required, string): Filename of umat file dynamically library.
name (optional, string, default "umat"): Name of material model (used for input to umat routine).

Reimplemented from oofem::FEMComponent.

Definition at line 74 of file abaqususermaterial.C.

References _IFT_AbaqusUserMaterial_initialStress, _IFT_AbaqusUserMaterial_name, _IFT_AbaqusUserMaterial_numericalTangent, _IFT_AbaqusUserMaterial_numericalTangentPerturbation, _IFT_AbaqusUserMaterial_numState, _IFT_AbaqusUserMaterial_properties, _IFT_AbaqusUserMaterial_userMaterial, cmname, filename, oofem::InputRecord::hasField(), oofem::StructuralMaterial::initializeFrom(), initialStress, IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, mPerturbation, mUseNumericalTangent, numState, OOFEM_ERROR, properties, umat, and umatobj.

Member Data Documentation

◆ abq2oo6

std::size_t const oofem::AbaqusUserMaterial::abq2oo6 = { 0, 1, 2, 5, 4, 3 }

staticprivate

Definition at line 124 of file abaqususermaterial.h.

Referenced by giveFirstPKStressVector_3d(), and giveRealStressVector_3d().

◆ abq2oo9

std::size_t const oofem::AbaqusUserMaterial::abq2oo9 = { 0, 1, 2, 5, 4, 3, 6, 8, 7 }

staticprivate

Definition at line 123 of file abaqususermaterial.h.

Referenced by giveFirstPKStressVector_3d(), and oofem::REGISTER_Material().

◆ cmname

char oofem::AbaqusUserMaterial::cmname[80]

private

Name for material routine.

Definition at line 96 of file abaqususermaterial.h.

Referenced by giveInputRecord(), and initializeFrom().

◆ filename

std::string oofem::AbaqusUserMaterial::filename

private

Name of the file that contains the umat function.

Definition at line 121 of file abaqususermaterial.h.

Referenced by giveInputRecord(), and initializeFrom().

◆ initialStress

FloatArrayF< 6 > oofem::AbaqusUserMaterial::initialStress

private

Initial stress.

Definition at line 103 of file abaqususermaterial.h.

Referenced by giveIPValue(), giveRealStressVector_3d(), and initializeFrom().

◆ mPerturbation

double oofem::AbaqusUserMaterial::mPerturbation = 1e-7

private

Size of perturbation if numerical tangent is used.

Definition at line 118 of file abaqususermaterial.h.

Referenced by give3dMaterialStiffnessMatrix_dPdF(), and initializeFrom().

◆ mStressInterpretation

int oofem::AbaqusUserMaterial::mStressInterpretation = 0

private

Flag to determine how the stress and Jacobian are interpreted. 0 implies that the P and dPdF are returned from the umat routine.

Definition at line 108 of file abaqususermaterial.h.

Referenced by giveFirstPKStressVector_3d().

◆ mUseNumericalTangent

bool oofem::AbaqusUserMaterial::mUseNumericalTangent = false

private

Flag to determine if numerical tangent should be used.

Definition at line 113 of file abaqususermaterial.h.

Referenced by give3dMaterialStiffnessMatrix_dPdF(), and initializeFrom().

◆ numState

int oofem::AbaqusUserMaterial::numState = 0

private

Size of the state vector.

Definition at line 98 of file abaqususermaterial.h.

Referenced by CreateStatus(), giveFirstPKStressVector_3d(), giveInputRecord(), giveRealStressVector_3d(), and initializeFrom().

◆ properties

FloatArray oofem::AbaqusUserMaterial::properties

private

Material properties.

Definition at line 101 of file abaqususermaterial.h.

Referenced by giveFirstPKStressVector_3d(), giveInputRecord(), giveRealStressVector_3d(), and initializeFrom().

◆ umat

void(* oofem::AbaqusUserMaterial::umat) (double *stress, double *statev, double *ddsdde, double *sse, double *spd, double *scd, double *rpl, double *ddsddt, double *drplde, double *drpldt, double *stran, double *dstran, double time[2], double *dtime, double *temp, double *dtemp, double predef[1], double dpred[1], char cmname[80], int *ndi, int *nshr, int *ntens, int *nstatv, double *props, int *nprops, double coords[3], double *drot, double *pnewdt, double *celent, double *dfgrd0, double *dfgrd1, int *noel, int *npt, int *layer, int *kspt, int *kstep, int *kinc) = nullptr

private

Pointer to the dynamically loaded umat-function (translated to C).

Definition at line 88 of file abaqususermaterial.h.

Referenced by giveFirstPKStressVector_3d(), giveRealStressVector_3d(), and initializeFrom().

◆ umatobj

void* oofem::AbaqusUserMaterial::umatobj = nullptr

private

Dynamically loaded umat.

Definition at line 85 of file abaqususermaterial.h.

Referenced by initializeFrom(), and ~AbaqusUserMaterial().

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

Public Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ AbaqusUserMaterial()

◆ ~AbaqusUserMaterial()

Member Function Documentation

◆ CreateStatus()

◆ give3dMaterialStiffnessMatrix()

◆ give3dMaterialStiffnessMatrix_dPdF()

◆ giveClassName()

◆ giveFirstPKStressVector_3d()

◆ giveInputRecord()

◆ giveInputRecordName()

◆ giveIPValue()

◆ givePlaneStrainStiffnessMatrix_dPdF()

◆ giveRealStressVector_3d()

◆ initializeFrom()

Member Data Documentation

◆ abq2oo6

◆ abq2oo9

◆ cmname

◆ filename

◆ initialStress

◆ mPerturbation

◆ mStressInterpretation

◆ mUseNumericalTangent

◆ numState

◆ properties

◆ umat

◆ umatobj