OOFEM
2.4
OOFEM.org - Object Oriented Finite Element Solver
|
This class implements an isotropic elastoplastic material with Rankine yield condition, associated flow rule and linear isotropic softening, and with isotropic damage that leads to softening. More...
#include <rankinemat.h>
Public Member Functions | |
RankineMat (int n, Domain *d) | |
virtual | ~RankineMat () |
double | evalYieldFunction (const FloatArray &sigPrinc, const double kappa) |
double | evalYieldStress (const double kappa) |
double | evalPlasticModulus (const double kappa) |
void | performPlasticityReturn (GaussPoint *gp, const FloatArray &totalStrain) |
double | computeDamage (GaussPoint *gp, TimeStep *tStep) |
double | computeDamageParam (double tempKappa) |
double | computeDamageParamPrime (double tempKappa) |
virtual void | computeCumPlastStrain (double &kappa, GaussPoint *gp, TimeStep *tStep) |
virtual int | hasMaterialModeCapability (MaterialMode mode) |
Tests if material supports material mode. More... | |
virtual IRResultType | initializeFrom (InputRecord *ir) |
Initializes receiver according to object description stored in input record. More... | |
virtual int | hasNonLinearBehaviour () |
Returns nonzero if receiver is non linear. More... | |
virtual const char * | giveInputRecordName () const |
virtual const char * | giveClassName () const |
LinearElasticMaterial * | giveLinearElasticMaterial () |
Returns a reference to the basic elastic material. More... | |
virtual bool | isCharacteristicMtrxSymmetric (MatResponseMode rMode) |
Returns true if stiffness matrix of receiver is symmetric Default implementation returns true. More... | |
virtual MaterialStatus * | CreateStatus (GaussPoint *gp) const |
Creates new copy of associated status and inserts it into given integration point. More... | |
virtual void | giveRealStressVector_PlaneStress (FloatArray &answer, GaussPoint *gp, const FloatArray &reducesStrain, TimeStep *tStep) |
Default implementation relies on giveRealStressVector_StressControl. More... | |
virtual void | giveRealStressVector_1d (FloatArray &answer, GaussPoint *gp, const FloatArray &totalStrain, TimeStep *tStep) |
Default implementation relies on giveRealStressVector_StressControl. More... | |
Public Member Functions inherited from oofem::StructuralMaterial | |
StructuralMaterial (int n, Domain *d) | |
Constructor. More... | |
virtual | ~StructuralMaterial () |
Destructor. More... | |
virtual void | giveInputRecord (DynamicInputRecord &input) |
Setups the input record string of receiver. More... | |
virtual void | giveStiffnessMatrix (FloatMatrix &answer, MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) |
Computes the stiffness matrix for giveRealStressVector of receiver in given integration point, respecting its history. More... | |
virtual void | giveRealStressVector (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedStrain, TimeStep *tStep) |
Computes the real stress vector for given total strain and integration point. More... | |
virtual void | giveRealStressVector_3d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep) |
Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress. More... | |
virtual void | giveRealStressVector_PlaneStrain (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep) |
Default implementation relies on giveRealStressVector_3d. More... | |
virtual void | giveRealStressVector_StressControl (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, const IntArray &strainControl, TimeStep *tStep) |
Iteratively calls giveRealStressVector_3d to find the stress controlled equal to zero· More... | |
virtual void | giveRealStressVector_ShellStressControl (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, const IntArray &strainControl, TimeStep *tStep) |
virtual void | giveRealStressVector_Warping (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep) |
Default implementation relies on giveRealStressVector_StressControl. More... | |
virtual void | giveRealStressVector_2dBeamLayer (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep) |
Default implementation relies on giveRealStressVector_StressControl. More... | |
virtual void | giveRealStressVector_PlateLayer (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep) |
Default implementation relies on giveRealStressVector_StressControl. More... | |
virtual void | giveRealStressVector_Fiber (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep) |
Default implementation relies on giveRealStressVector_StressControl. More... | |
virtual void | giveRealStressVector_Lattice2d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep) |
virtual void | giveRealStressVector_Lattice3d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep) |
virtual void | giveRealStressVector_2dPlateSubSoil (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep) |
Default implementation is not provided. More... | |
virtual void | giveRealStressVector_3dBeamSubSoil (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedE, TimeStep *tStep) |
virtual void | giveEshelbyStressVector_PlaneStrain (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep) |
Prototype for computation of Eshelby stress. More... | |
void | give_dPdF_from (const FloatMatrix &dSdE, FloatMatrix &answer, GaussPoint *gp) |
void | convert_dSdE_2_dPdF (FloatMatrix &answer, const FloatMatrix &dSdE, const FloatArray &S, const FloatArray &F, MaterialMode matMode) |
virtual void | giveThermalDilatationVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep) |
Returns a vector of coefficients of thermal dilatation in direction of each material principal (local) axis. More... | |
double | giveReferenceTemperature () |
Returns the reference temperature of receiver. More... | |
virtual void | computeStressIndependentStrainVector (FloatArray &answer, GaussPoint *gp, TimeStep *tStep, ValueModeType mode) |
Computes reduced strain vector in given integration point, generated by internal processes in material, which are independent on loading in particular integration point. More... | |
virtual void | computeStressIndependentStrainVector_3d (FloatArray &answer, GaussPoint *gp, TimeStep *tStep, ValueModeType mode) |
virtual void | give3dMaterialStiffnessMatrix (FloatMatrix &answer, MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) |
Computes full 3d material stiffness matrix at given integration point, time, respecting load history in integration point. More... | |
virtual void | give3dMaterialStiffnessMatrix_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) |
Method for subtracting from reduced space strain vector its stress-independent parts (caused by temperature, shrinkage, creep and possibly by other phenomena). More... | |
void | giveStressDependentPartOfStrainVector_3d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedStrainVector, TimeStep *tStep, ValueModeType mode) |
virtual int | setIPValue (const FloatArray &value, GaussPoint *gp, InternalStateType type) |
Sets the value of a certain variable at a given integration point to the given value. More... | |
virtual void | give2dBeamLayerStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) |
Method for computing 2d beam layer stiffness matrix of receiver. More... | |
virtual void | givePlateLayerStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) |
Method for computing 2d plate layer stiffness matrix of receiver. More... | |
virtual void | giveFiberStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) |
Method for computing 1d fiber stiffness matrix of receiver. More... | |
virtual void | give2dLatticeStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) |
Method for computing 2d lattice stiffness matrix of receiver. More... | |
virtual void | give3dLatticeStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) |
Method for computing 3d lattice stiffness matrix of receiver. More... | |
virtual void | give2dPlateSubSoilStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) |
Method for computing stiffness matrix of plate subsoil model. More... | |
virtual void | give3dBeamSubSoilStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) |
Method for computing stiffness matrix of beam3d subsoil model. More... | |
virtual void | giveFirstPKStressVector_3d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep) |
Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress. More... | |
virtual void | giveFirstPKStressVector_PlaneStrain (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep) |
Default implementation relies on giveFirstPKStressVector_3d. More... | |
virtual void | giveFirstPKStressVector_PlaneStress (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep) |
Default implementation relies on giveFirstPKStressVector_3d. More... | |
virtual void | giveFirstPKStressVector_1d (FloatArray &answer, GaussPoint *gp, const FloatArray &reducedF, TimeStep *tStep) |
Default implementation relies on giveFirstPKStressVector_3d. More... | |
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 | givePlaneStressStiffMtrx_dPdF (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) |
virtual void | givePlaneStressStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) |
virtual void | givePlaneStrainStiffMtrx (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) |
Method for computing plane strain stiffness matrix of receiver. More... | |
virtual void | givePlaneStrainStiffMtrx_dPdF (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) |
virtual void | givePlaneStrainStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) |
virtual void | give1dStressStiffMtrx_dPdF (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) |
virtual void | give1dStressStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) |
Public Member Functions inherited from oofem::Material | |
Material (int n, Domain *d) | |
Constructor. More... | |
virtual | ~Material () |
Destructor. More... | |
virtual double | give (int aProperty, GaussPoint *gp) |
Returns the value of material property 'aProperty'. More... | |
virtual bool | hasProperty (int aProperty, GaussPoint *gp) |
Returns true if 'aProperty' exists on material. More... | |
virtual void | modifyProperty (int aProperty, double value, GaussPoint *gp) |
Modify 'aProperty', which already exists on material. More... | |
double | giveCastingTime () |
virtual bool | isActivated (TimeStep *tStep) |
virtual int | hasCastingTimeSupport () |
Tests if material supports casting time. More... | |
virtual void | printYourself () |
Prints receiver state on stdout. Useful for debugging. More... | |
virtual contextIOResultType | saveIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp) |
Stores integration point state to output stream. More... | |
virtual contextIOResultType | restoreIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp) |
Reads integration point state to output stream. More... | |
virtual int | checkConsistency () |
Allows programmer to test some internal data, before computation begins. More... | |
virtual int | initMaterial (Element *element) |
Optional function to call specific procedures when initializing a material. More... | |
virtual MaterialStatus * | giveStatus (GaussPoint *gp) const |
Returns material status of receiver in given integration point. More... | |
virtual int | packUnknowns (DataStream &buff, TimeStep *tStep, GaussPoint *ip) |
Pack all necessary data of integration point (according to element parallel_mode) into given communication buffer. More... | |
virtual int | unpackAndUpdateUnknowns (DataStream &buff, TimeStep *tStep, GaussPoint *ip) |
Unpack and updates all necessary data of given integration point (according to element parallel_mode) into given communication buffer. More... | |
virtual int | estimatePackSize (DataStream &buff, GaussPoint *ip) |
Estimates the necessary pack size to hold all packed data of receiver. More... | |
virtual double | predictRelativeComputationalCost (GaussPoint *gp) |
Returns the weight representing relative computational cost of receiver The reference material model is linear isotropic material - its weight is set to 1.0 The other material models should compare to this reference model. More... | |
virtual double | predictRelativeRedistributionCost (GaussPoint *gp) |
Returns the relative redistribution cost of the receiver. More... | |
virtual void | initTempStatus (GaussPoint *gp) |
Initializes temporary variables stored in integration point status at the beginning of new time step. More... | |
Public Member Functions inherited from oofem::FEMComponent | |
FEMComponent (int n, Domain *d) | |
Regular constructor, creates component with given number and belonging to given domain. More... | |
virtual | ~FEMComponent () |
Virtual destructor. More... | |
Domain * | giveDomain () const |
virtual void | setDomain (Domain *d) |
Sets associated Domain. More... | |
int | giveNumber () const |
void | setNumber (int num) |
Sets number of receiver. More... | |
virtual void | updateLocalNumbering (EntityRenumberingFunctor &f) |
Local renumbering support. More... | |
virtual contextIOResultType | saveContext (DataStream &stream, ContextMode mode, void *obj=NULL) |
Stores receiver state to output stream. More... | |
virtual contextIOResultType | restoreContext (DataStream &stream, ContextMode mode, void *obj=NULL) |
Restores the receiver state previously written in stream. More... | |
virtual void | printOutputAt (FILE *file, TimeStep *tStep) |
Prints output of receiver to stream, for given time step. More... | |
virtual Interface * | giveInterface (InterfaceType t) |
Interface requesting service. More... | |
std::string | errorInfo (const char *func) const |
Returns string for prepending output (used by error reporting macros). More... | |
Protected Member Functions | |
virtual void | give1dStressStiffMtrx (FloatMatrix &answer, MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) |
Method for computing 1d stiffness matrix of receiver. More... | |
virtual void | givePlaneStressStiffMtrx (FloatMatrix &answer, MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) |
Method for computing plane stress stiffness matrix of receiver. More... | |
void | evaluatePlaneStressStiffMtrx (FloatMatrix &answer, MatResponseMode mode, GaussPoint *gp, TimeStep *tStep, double gprime) |
Executive method used by local and gradient version. More... | |
void | computeEta (FloatArray &answer, RankineMatStatus *status) |
Computes derivatives of final kappa with respect to final strain. More... | |
virtual int | giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) |
Returns the integration point corresponding value in Reduced form. More... | |
Protected Attributes | |
LinearElasticMaterial * | linearElasticMaterial |
Reference to the basic elastic material. More... | |
double | E |
Young's modulus. More... | |
double | nu |
Poisson's ratio. More... | |
double | H0 |
Initial hardening modulus. More... | |
int | plasthardtype |
Type of plastic hardening (0=linear, 1=exponential) More... | |
double | delSigY |
Final increment of yield stress (at infinite cumulative plastic strain) More... | |
double | sig0 |
Initial (uniaxial) yield stress. More... | |
double | yieldtol |
Relative tolerance in yield condition. More... | |
double | a |
Parameter that controls damage evolution (a=0 turns damage off). More... | |
double | ep |
Total strain at peak stress sig0–Used only if plasthardtype=2. More... | |
double | md |
Exponent in hardening law–Used only if plasthardtype=2. More... | |
int | damlaw |
type of damage law (0=exponential, 1=exponential and damage starts after peak stress sig0) More... | |
double | param1 |
coefficient required when damlaw=1 or 2 More... | |
double | param2 |
coefficient required when damlaw=1 or 2 More... | |
double | param3 |
coefficient required when damlaw=2 More... | |
double | param4 |
coefficient required when damlaw=2 More... | |
double | param5 |
coefficient required when damlaw=2 More... | |
Protected Attributes inherited from oofem::StructuralMaterial | |
double | referenceTemperature |
Reference temperature (temperature, when material has been built into structure). More... | |
Protected Attributes inherited from oofem::Material | |
Dictionary | propertyDictionary |
Property dictionary. More... | |
double | castingTime |
Casting time. More... | |
Protected Attributes inherited from oofem::FEMComponent | |
int | number |
Component number. More... | |
Domain * | domain |
Link to domain object, useful for communicating with other FEM components. More... | |
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 int | giveVoigtVectorMask (IntArray &answer, MaterialMode mmode) |
Returns a mask of the vector indicies corresponding to components in a general (non-symmetric) second order tensor of some stress/strain/deformation measure that performes work. More... | |
static int | giveVoigtSymVectorMask (IntArray &answer, MaterialMode mmode) |
The same as giveVoigtVectorMask but returns a mask corresponding to a symmetric second order tensor. More... | |
static void | giveInvertedVoigtVectorMask (IntArray &answer, MaterialMode mmode) |
Gives the inverted version of giveVoigtVectorMask. More... | |
static int | giveSizeOfVoigtVector (MaterialMode mmode) |
Returns the size of reduced stress/strain vector according to given mode. More... | |
static int | giveSizeOfVoigtSymVector (MaterialMode mmode) |
Returns the size of symmetric part of a reduced stress/strain vector according to given mode. More... | |
static void | giveFullVectorForm (FloatArray &answer, const FloatArray &strainVector, MaterialMode matMode) |
Converts the reduced symmetric Voigt vector (2nd order tensor) to full form. More... | |
static void | giveFullVectorFormF (FloatArray &answer, const FloatArray &strainVector, MaterialMode matMode) |
Converts the reduced deformation gradient Voigt vector (2nd order tensor). More... | |
static void | giveFullSymVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode) |
Converts the reduced unsymmetric Voigt vector (2nd order tensor) to full form. More... | |
static void | giveReducedVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode) |
Converts the full symmetric Voigt vector (2nd order tensor) to reduced form. More... | |
static void | giveReducedSymVectorForm (FloatArray &answer, const FloatArray &vec, MaterialMode matMode) |
Converts the full unsymmetric Voigt vector (2nd order tensor) to reduced form. More... | |
static void | giveFullSymMatrixForm (FloatMatrix &answer, const FloatMatrix &red, MaterialMode matMode) |
Converts the full unsymmetric Voigt matrix (4th order tensor) to reduced form. More... | |
static void | giveReducedMatrixForm (FloatMatrix &answer, const FloatMatrix &full, MaterialMode matMode) |
Converts the full symmetric Voigt matrix (4th order tensor) to reduced form. More... | |
static void | giveReducedSymMatrixForm (FloatMatrix &answer, const FloatMatrix &full, MaterialMode matMode) |
Converts the full unsymmetric Voigt matrix (4th order tensor) to reduced form. More... | |
static void | transformStrainVectorTo (FloatArray &answer, const FloatMatrix &base, const FloatArray &strainVector, bool transpose=false) |
Transforms 3d strain vector into another coordinate system. More... | |
static void | transformStressVectorTo (FloatArray &answer, const FloatMatrix &base, const FloatArray &stressVector, bool transpose=false) |
Transforms 3d stress vector into another coordinate system. More... | |
static double | computeVonMisesStress (const FloatArray *currentStress) |
Computes equivalent of von Mises stress. More... | |
static void | giveStrainVectorTranformationMtrx (FloatMatrix &answer, const FloatMatrix &base, bool transpose=false) |
Computes 3d strain vector transformation matrix from standard vector transformation matrix. More... | |
static void | give2DStrainVectorTranformationMtrx (FloatMatrix &answer, const FloatMatrix &base, bool transpose=false) |
Computes 2d strain vector transformation matrix from standard vector transformation matrix. More... | |
static void | giveStressVectorTranformationMtrx (FloatMatrix &answer, const FloatMatrix &base, bool transpose=false) |
Computes 3d stress vector transformation matrix from standard vector transformation matrix. More... | |
static void | givePlaneStressVectorTranformationMtrx (FloatMatrix &answer, const FloatMatrix &base, bool transpose=false) |
Computes 2d stress vector transformation matrix from standard vector transformation matrix. More... | |
static void | sortPrincDirAndValCloseTo (FloatArray *pVal, FloatMatrix *pDir, FloatMatrix *toPDir) |
Method for sorting newly computed principal values (pVal) and corresponding principal directions (pDir) to be closed to some (often previous) principal directions (toPDir). More... | |
static void | computePrincipalValues (FloatArray &answer, const FloatArray &s, stressStrainPrincMode mode) |
Common functions for convenience. More... | |
static void | computePrincipalValDir (FloatArray &answer, FloatMatrix &dir, const FloatArray &s, stressStrainPrincMode mode) |
Computes principal values and directions of stress or strain vector. More... | |
static double | computeDeviatoricVolumetricSplit (FloatArray &dev, const FloatArray &s) |
Computes split of receiver into deviatoric and volumetric part. More... | |
static void | computeDeviatoricVolumetricSum (FloatArray &s, const FloatArray &dev, double mean) |
static void | applyDeviatoricElasticCompliance (FloatArray &strain, const FloatArray &stress, double EModulus, double nu) |
static void | applyDeviatoricElasticCompliance (FloatArray &strain, const FloatArray &stress, double GModulus) |
static void | applyDeviatoricElasticStiffness (FloatArray &stress, const FloatArray &strain, double EModulus, double nu) |
static void | applyDeviatoricElasticStiffness (FloatArray &stress, const FloatArray &strain, double GModulus) |
static void | applyElasticStiffness (FloatArray &stress, const FloatArray &strain, double EModulus, double nu) |
static void | applyElasticCompliance (FloatArray &strain, const FloatArray &stress, double EModulus, double nu) |
static double | computeStressNorm (const FloatArray &stress) |
static double | computeFirstInvariant (const FloatArray &s) |
static double | computeSecondStressInvariant (const FloatArray &s) |
static double | computeThirdStressInvariant (const FloatArray &s) |
static double | computeFirstCoordinate (const FloatArray &s) |
static double | computeSecondCoordinate (const FloatArray &s) |
static double | computeThirdCoordinate (const FloatArray &s) |
Static Public Attributes inherited from oofem::StructuralMaterial | |
static std::vector< std::vector< int > > | vIindex |
Voigt index map. More... | |
static std::vector< std::vector< int > > | svIndex |
Symmetric Voigt index map. More... | |
This class implements an isotropic elastoplastic material with Rankine yield condition, associated flow rule and linear isotropic softening, and with isotropic damage that leads to softening.
It differs from other similar materials (such as RankinePlasticMaterial) by implementation - here we use an efficient algorithm for this specific model and implement the plane stress case and the vertex return. Also, hardening and softening is incorporated.
Definition at line 85 of file rankinemat.h.
oofem::RankineMat::RankineMat | ( | int | n, |
Domain * | d | ||
) |
|
inlinevirtual |
Definition at line 142 of file rankinemat.h.
References computeCumPlastStrain(), computeDamage(), computeDamageParam(), computeDamageParamPrime(), evalPlasticModulus(), evalYieldFunction(), evalYieldStress(), hasMaterialModeCapability(), initializeFrom(), and performPlasticityReturn().
|
virtual |
Reimplemented in oofem::RankineMatGrad.
Definition at line 484 of file rankinemat.C.
References oofem::Material::giveStatus(), and oofem::RankineMatStatus::giveTempCumulativePlasticStrain().
Referenced by computeDamage(), oofem::RankineMatNl::computeLocalCumPlasticStrain(), and ~RankineMat().
double oofem::RankineMat::computeDamage | ( | GaussPoint * | gp, |
TimeStep * | tStep | ||
) |
Definition at line 470 of file rankinemat.C.
References computeCumPlastStrain(), computeDamageParam(), oofem::RankineMatStatus::giveDamage(), and oofem::Material::giveStatus().
Referenced by giveRealStressVector_1d(), giveRealStressVector_PlaneStress(), oofem::RankineMatGrad::giveRealStressVectorGrad(), and ~RankineMat().
double oofem::RankineMat::computeDamageParam | ( | double | tempKappa | ) |
Definition at line 436 of file rankinemat.C.
References a, damlaw, ep, param1, param2, param3, param4, and param5.
Referenced by oofem::RankineMatNl::computeDamage(), computeDamage(), and ~RankineMat().
double oofem::RankineMat::computeDamageParamPrime | ( | double | tempKappa | ) |
Definition at line 453 of file rankinemat.C.
References a, damlaw, ep, param1, param2, param3, param4, and param5.
Referenced by oofem::RankineMatNl::giveLocalNonlocalStiffnessContribution(), oofem::RankineMatGrad::givePlaneStressGprime(), oofem::RankineMatGrad::givePlaneStressStiffMtrx(), oofem::RankineMatNl::givePlaneStressStiffMtrx(), givePlaneStressStiffMtrx(), and ~RankineMat().
|
protected |
Computes derivatives of final kappa with respect to final strain.
Definition at line 612 of file rankinemat.C.
References oofem::FloatArray::at(), oofem::FloatArray::beProductOf(), oofem::StressVector::computePrincipalValDir(), E, evalPlasticModulus(), oofem::RankineMatStatus::giveDKappa(), oofem::StructuralMaterial::givePlaneStressVectorTranformationMtrx(), oofem::RankineMatStatus::giveTempCumulativePlasticStrain(), oofem::RankineMatStatus::giveTempEffectiveStress(), and nu.
Referenced by oofem::RankineMatNl::giveRemoteNonlocalStiffnessContribution(), and isCharacteristicMtrxSymmetric().
|
virtual |
Creates new copy of associated status and inserts it into given integration point.
gp | Integration point where newly created status will be stored. |
Reimplemented from oofem::Material.
Reimplemented in oofem::RankineMatNl, and oofem::RankineMatGrad.
Definition at line 166 of file rankinemat.C.
References oofem::FEMComponent::giveDomain().
Referenced by isCharacteristicMtrxSymmetric().
double oofem::RankineMat::evalPlasticModulus | ( | const double | kappa | ) |
Definition at line 261 of file rankinemat.C.
References delSigY, E, ep, H0, md, and plasthardtype.
Referenced by computeEta(), evaluatePlaneStressStiffMtrx(), oofem::RankineMatGrad::givePlaneStressKappaMatrix(), performPlasticityReturn(), and ~RankineMat().
|
protected |
Executive method used by local and gradient version.
(with different parameters gprime)
Definition at line 521 of file rankinemat.C.
References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::StressVector::computePrincipalValDir(), E, evalPlasticModulus(), oofem::RankineMatStatus::giveCumulativePlasticStrain(), oofem::RankineMatStatus::giveDKappa(), giveLinearElasticMaterial(), oofem::StructuralMaterial::givePlaneStressVectorTranformationMtrx(), oofem::Material::giveStatus(), oofem::StructuralMaterial::giveStiffnessMatrix(), oofem::RankineMatStatus::giveTangentShearStiffness(), oofem::RankineMatStatus::giveTempCumulativePlasticStrain(), oofem::RankineMatStatus::giveTempDamage(), oofem::RankineMatStatus::giveTempEffectiveStress(), nu, oofem::FloatMatrix::resize(), oofem::FloatMatrix::rotatedWith(), oofem::FloatMatrix::subtract(), oofem::FloatMatrix::times(), and oofem::FloatMatrix::zero().
Referenced by oofem::RankineMatGrad::givePlaneStressStiffMtrx(), oofem::RankineMatNl::givePlaneStressStiffMtrx(), givePlaneStressStiffMtrx(), and isCharacteristicMtrxSymmetric().
double oofem::RankineMat::evalYieldFunction | ( | const FloatArray & | sigPrinc, |
const double | kappa | ||
) |
Definition at line 231 of file rankinemat.C.
References oofem::FloatArray::at(), and evalYieldStress().
Referenced by performPlasticityReturn(), and ~RankineMat().
double oofem::RankineMat::evalYieldStress | ( | const double | kappa | ) |
Definition at line 238 of file rankinemat.C.
References delSigY, E, ep, H0, md, plasthardtype, and sig0.
Referenced by evalYieldFunction(), performPlasticityReturn(), and ~RankineMat().
|
protectedvirtual |
Method for computing 1d stiffness matrix of receiver.
Default implementation computes 3d stiffness matrix using give3dMaterialStiffnessMatrix and reduces it to 1d stiffness using reduce method described above. However, this reduction is quite time consuming and if it is possible, it is recommended to overload this method and provide direct method for computing particular stiffness matrix.
answer | Stiffness matrix. |
mmode | Material response mode. |
gp | Integration point, which load history is used. |
tStep | Time step (most models are able to respond only when tStep is current time step). |
Reimplemented from oofem::StructuralMaterial.
Definition at line 504 of file rankinemat.C.
References oofem::FloatMatrix::at(), E, oofem::Material::giveStatus(), oofem::RankineMatStatus::giveTempDamage(), OOFEM_ERROR, oofem::FloatMatrix::resize(), and oofem::FloatMatrix::times().
Referenced by isCharacteristicMtrxSymmetric().
|
inlinevirtual |
Reimplemented from oofem::StructuralMaterial.
Reimplemented in oofem::RankineMatNl, and oofem::RankineMatGrad.
Definition at line 160 of file rankinemat.h.
|
inlinevirtual |
Implements oofem::FEMComponent.
Reimplemented in oofem::RankineMatNl, and oofem::RankineMatGrad.
Definition at line 159 of file rankinemat.h.
References _IFT_RankineMat_Name.
|
protectedvirtual |
Returns the integration point corresponding value in Reduced form.
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. |
Reimplemented from oofem::StructuralMaterial.
Reimplemented in oofem::RankineMatNl, and oofem::RankineMatGrad.
Definition at line 650 of file rankinemat.C.
References oofem::FloatArray::at(), ep, oofem::RankineMatStatus::giveCumulativePlasticStrain(), oofem::RankineMatStatus::giveDamage(), oofem::RankineMatStatus::giveDissWork(), oofem::StructuralMaterial::giveIPValue(), oofem::RankineMatStatus::givePlasDef(), oofem::Material::giveStatus(), oofem::RankineMatStatus::giveStressWork(), oofem::FloatArray::resize(), and oofem::FloatArray::zero().
Referenced by oofem::RankineMatGrad::giveIPValue(), oofem::RankineMatNl::giveIPValue(), and isCharacteristicMtrxSymmetric().
|
inline |
Returns a reference to the basic elastic material.
Definition at line 163 of file rankinemat.h.
References linearElasticMaterial.
Referenced by evaluatePlaneStressStiffMtrx(), and oofem::RankineMatNl::givePlaneStressStiffMtrx().
|
protectedvirtual |
Method for computing plane stress stiffness matrix of receiver.
Default implementation computes 3d stiffness matrix using give3dMaterialStiffnessMatrix and reduces it to plane stress stiffness using reduce method described above. However, this reduction is quite time consuming and if it is possible, it is recommended to overload this method and provide direct method for computing particular stiffness matrix.
answer | Stiffness matrix. |
mmode | Material response mode. |
gp | Integration point, which load history is used. |
tStep | Time step (most models are able to respond only when tStep is current time step). |
Reimplemented from oofem::StructuralMaterial.
Reimplemented in oofem::RankineMatNl, and oofem::RankineMatGrad.
Definition at line 492 of file rankinemat.C.
References computeDamageParamPrime(), evaluatePlaneStressStiffMtrx(), oofem::Material::giveStatus(), and oofem::RankineMatStatus::giveTempCumulativePlasticStrain().
Referenced by isCharacteristicMtrxSymmetric().
|
virtual |
Default implementation relies on giveRealStressVector_StressControl.
Reimplemented from oofem::StructuralMaterial.
Reimplemented in oofem::RankineMatNl.
Definition at line 176 of file rankinemat.C.
References oofem::FloatArray::beScaled(), computeDamage(), oofem::RankineMatStatus::computeWork_1d(), E, oofem::Material::giveStatus(), oofem::RankineMatStatus::giveTempEffectiveStress(), oofem::Material::initTempStatus(), oofem::StructuralMaterialStatus::letTempStrainVectorBe(), oofem::StructuralMaterialStatus::letTempStressVectorBe(), performPlasticityReturn(), oofem::RankineMatStatus::setTempDamage(), and sig0.
Referenced by isCharacteristicMtrxSymmetric().
|
virtual |
Default implementation relies on giveRealStressVector_StressControl.
Reimplemented from oofem::StructuralMaterial.
Reimplemented in oofem::RankineMatNl.
Definition at line 202 of file rankinemat.C.
References oofem::FloatArray::beScaled(), computeDamage(), oofem::RankineMatStatus::computeWork_PlaneStress(), E, oofem::Material::giveStatus(), oofem::RankineMatStatus::giveTempEffectiveStress(), oofem::Material::initTempStatus(), oofem::StructuralMaterialStatus::letTempStrainVectorBe(), oofem::StructuralMaterialStatus::letTempStressVectorBe(), performPlasticityReturn(), oofem::RankineMatStatus::setTempDamage(), and sig0.
Referenced by isCharacteristicMtrxSymmetric().
|
virtual |
Tests if material supports material mode.
mode | Required material mode. |
Reimplemented from oofem::StructuralMaterial.
Reimplemented in oofem::RankineMatGrad.
Definition at line 73 of file rankinemat.C.
Referenced by ~RankineMat().
|
inlinevirtual |
Returns nonzero if receiver is non linear.
Reimplemented from oofem::Material.
Definition at line 158 of file rankinemat.h.
|
virtual |
Initializes receiver according to object description stored in input record.
This function is called immediately after creating object using constructor. Input record can be imagined as data record in component database belonging to receiver. Receiver may use value-name extracting functions to extract particular field from record.
ir | Input record to initialize from. |
Reimplemented from oofem::StructuralMaterial.
Reimplemented in oofem::RankineMatNl, and oofem::RankineMatGrad.
Definition at line 81 of file rankinemat.C.
References _IFT_RankineMat_a, _IFT_RankineMat_damlaw, _IFT_RankineMat_delsigy, _IFT_RankineMat_ep, _IFT_RankineMat_gf, _IFT_RankineMat_h, _IFT_RankineMat_param1, _IFT_RankineMat_param2, _IFT_RankineMat_param3, _IFT_RankineMat_param4, _IFT_RankineMat_param5, _IFT_RankineMat_plasthardtype, _IFT_RankineMat_sig0, _IFT_RankineMat_yieldtol, a, damlaw, delSigY, E, ep, H0, oofem::LinearElasticMaterial::initializeFrom(), oofem::StructuralMaterial::initializeFrom(), IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, oofem::IRRT_BAD_FORMAT, oofem::IRRT_OK, linearElasticMaterial, md, nu, OOFEM_ERROR, OOFEM_WARNING, param1, param2, param3, param4, param5, plasthardtype, sig0, and yieldtol.
Referenced by oofem::RankineMatGrad::initializeFrom(), oofem::RankineMatNl::initializeFrom(), and ~RankineMat().
|
inlinevirtual |
Returns true if stiffness matrix of receiver is symmetric Default implementation returns true.
Reimplemented from oofem::Material.
Definition at line 165 of file rankinemat.h.
References computeEta(), CreateStatus(), evaluatePlaneStressStiffMtrx(), give1dStressStiffMtrx(), giveIPValue(), givePlaneStressStiffMtrx(), giveRealStressVector_1d(), and giveRealStressVector_PlaneStress().
void oofem::RankineMat::performPlasticityReturn | ( | GaussPoint * | gp, |
const FloatArray & | totalStrain | ||
) |
Definition at line 283 of file rankinemat.C.
References oofem::StrainVector::applyElasticStiffness(), oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::StressVector::computePrincipalValDir(), E, evalPlasticModulus(), evalYieldFunction(), evalYieldStress(), oofem::RankineMatStatus::giveCumulativePlasticStrain(), oofem::RankineMatStatus::giveDKappa(), oofem::GaussPoint::giveMaterialMode(), oofem::RankineMatStatus::givePlasticStrain(), oofem::Material::giveStatus(), oofem::RankineMatStatus::letTempEffectiveStressBe(), oofem::RankineMatStatus::letTempPlasticStrainBe(), nu, OOFEM_ERROR, oofem::FloatArray::pY(), oofem::RankineMatStatus::setDKappa(), oofem::RankineMatStatus::setTangentShearStiffness(), oofem::RankineMatStatus::setTempCumulativePlasticStrain(), sig0, oofem::FloatArray::subtract(), and yieldtol.
Referenced by giveRealStressVector_1d(), giveRealStressVector_PlaneStress(), oofem::RankineMatGrad::giveRealStressVectorGrad(), oofem::RankineMatNl::updateBeforeNonlocAverage(), and ~RankineMat().
|
protected |
Parameter that controls damage evolution (a=0 turns damage off).
Definition at line 113 of file rankinemat.h.
Referenced by computeDamageParam(), computeDamageParamPrime(), and initializeFrom().
|
protected |
type of damage law (0=exponential, 1=exponential and damage starts after peak stress sig0)
Definition at line 122 of file rankinemat.h.
Referenced by computeDamageParam(), computeDamageParamPrime(), initializeFrom(), and RankineMat().
|
protected |
Final increment of yield stress (at infinite cumulative plastic strain)
Definition at line 104 of file rankinemat.h.
Referenced by evalPlasticModulus(), evalYieldStress(), initializeFrom(), and RankineMat().
|
protected |
Young's modulus.
Definition at line 92 of file rankinemat.h.
Referenced by computeEta(), evalPlasticModulus(), evaluatePlaneStressStiffMtrx(), evalYieldStress(), give1dStressStiffMtrx(), oofem::RankineMatGrad::givePlaneStressKappaMatrix(), oofem::RankineMatNl::giveRealStressVector_1d(), giveRealStressVector_1d(), oofem::RankineMatNl::giveRealStressVector_PlaneStress(), giveRealStressVector_PlaneStress(), oofem::RankineMatGrad::giveRealStressVectorGrad(), initializeFrom(), performPlasticityReturn(), and RankineMat().
|
protected |
Total strain at peak stress sig0–Used only if plasthardtype=2.
Definition at line 116 of file rankinemat.h.
Referenced by computeDamageParam(), computeDamageParamPrime(), evalPlasticModulus(), evalYieldStress(), giveIPValue(), initializeFrom(), and RankineMat().
|
protected |
Initial hardening modulus.
Definition at line 98 of file rankinemat.h.
Referenced by evalPlasticModulus(), evalYieldStress(), initializeFrom(), and RankineMat().
|
protected |
Reference to the basic elastic material.
Definition at line 89 of file rankinemat.h.
Referenced by giveLinearElasticMaterial(), initializeFrom(), and RankineMat().
|
protected |
Exponent in hardening law–Used only if plasthardtype=2.
Definition at line 119 of file rankinemat.h.
Referenced by evalPlasticModulus(), evalYieldStress(), initializeFrom(), and RankineMat().
|
protected |
Poisson's ratio.
Definition at line 95 of file rankinemat.h.
Referenced by computeEta(), evaluatePlaneStressStiffMtrx(), oofem::RankineMatGrad::givePlaneStressKappaMatrix(), initializeFrom(), performPlasticityReturn(), and RankineMat().
|
protected |
coefficient required when damlaw=1 or 2
Definition at line 125 of file rankinemat.h.
Referenced by computeDamageParam(), computeDamageParamPrime(), initializeFrom(), and RankineMat().
|
protected |
coefficient required when damlaw=1 or 2
Definition at line 128 of file rankinemat.h.
Referenced by computeDamageParam(), computeDamageParamPrime(), initializeFrom(), and RankineMat().
|
protected |
coefficient required when damlaw=2
Definition at line 131 of file rankinemat.h.
Referenced by computeDamageParam(), computeDamageParamPrime(), initializeFrom(), and RankineMat().
|
protected |
coefficient required when damlaw=2
Definition at line 134 of file rankinemat.h.
Referenced by computeDamageParam(), computeDamageParamPrime(), initializeFrom(), and RankineMat().
|
protected |
coefficient required when damlaw=2
Definition at line 137 of file rankinemat.h.
Referenced by computeDamageParam(), computeDamageParamPrime(), initializeFrom(), and RankineMat().
|
protected |
Type of plastic hardening (0=linear, 1=exponential)
Definition at line 101 of file rankinemat.h.
Referenced by evalPlasticModulus(), evalYieldStress(), and initializeFrom().
|
protected |
Initial (uniaxial) yield stress.
Definition at line 107 of file rankinemat.h.
Referenced by evalYieldStress(), oofem::RankineMatNl::giveRealStressVector_1d(), giveRealStressVector_1d(), oofem::RankineMatNl::giveRealStressVector_PlaneStress(), giveRealStressVector_PlaneStress(), oofem::RankineMatGrad::giveRealStressVectorGrad(), initializeFrom(), performPlasticityReturn(), and RankineMat().
|
protected |
Relative tolerance in yield condition.
Definition at line 110 of file rankinemat.h.
Referenced by initializeFrom(), and performPlasticityReturn().