|
OOFEM 3.0
|
#include <mplasticmaterial.h>
Public Member Functions | |
| MPlasticMaterial (int n, Domain *d) | |
| virtual | ~MPlasticMaterial () |
| bool | hasMaterialModeCapability (MaterialMode mode) const override |
| const char * | giveClassName () const override |
| LinearElasticMaterial * | giveLinearElasticMaterial () |
| Returns reference to undamaged (bulk) material. | |
| bool | isCharacteristicMtrxSymmetric (MatResponseMode rMode) const override |
| FloatMatrixF< 6, 6 > | give3dMaterialStiffnessMatrix (MatResponseMode, GaussPoint *gp, TimeStep *tStep) const override |
| void | giveRealStressVector (FloatArray &answer, GaussPoint *gp, const FloatArray &, 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< 4 > | giveRealStressVector_PlaneStrain (const FloatArrayF< 4 > &strain, GaussPoint *gp, TimeStep *tStep) const override |
| Default implementation relies on giveRealStressVector_3d. | |
| FloatArrayF< 3 > | giveRealStressVector_PlaneStress (const FloatArrayF< 3 > &strain, GaussPoint *gp, TimeStep *tStep) const override |
| Default implementation relies on giveRealStressVector_StressControl. | |
| FloatArrayF< 1 > | giveRealStressVector_1d (const FloatArrayF< 1 > &strain, GaussPoint *gp, TimeStep *tStep) const override |
| Default implementation relies on giveRealStressVector_StressControl. | |
| int | giveIPValue (FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) override |
| virtual int | giveSizeOfFullHardeningVarsVector () const |
| virtual int | giveSizeOfReducedHardeningVarsVector (GaussPoint *) const |
| std::unique_ptr< MaterialStatus > | CreateStatus (GaussPoint *gp) 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 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< 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< 9 > | giveFirstPKStressVector_3d (const FloatArrayF< 9 > &vF, GaussPoint *gp, TimeStep *tStep) const |
| Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress. | |
| 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 FloatMatrixF< 9, 9 > | give3dMaterialStiffnessMatrix_dPdF (MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) const |
| 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< 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< 5, 5 > | givePlaneStrainStiffnessMatrix_dPdF (MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) const |
| virtual void | givePlaneStrainStiffMtrx_dCde (FloatMatrix &answer, MatResponseMode mmode, GaussPoint *gp, TimeStep *tStep) |
| 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< 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 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. | |
| virtual const char * | giveInputRecordName () const =0 |
| 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). | |
Protected Types | |
| enum | ReturnMappingAlgoType { mpm_ClosestPoint , mpm_CuttingPlane } |
| Protected type to determine the return mapping algorithm. More... | |
| enum | functType { yieldFunction , loadFunction } |
| Type that allows to distinguish between yield function and loading function. More... | |
| enum | plastType { associatedPT , nonassociatedPT } |
Protected Member Functions | |
| void | closestPointReturn (FloatArray &answer, IntArray &activeConditionMap, FloatArray &gamma, GaussPoint *gp, const FloatArray &totalStrain, FloatArray &plasticStrainR, FloatArray &strainSpaceHardeningVariables, TimeStep *tStep) const |
| void | cuttingPlaneReturn (FloatArray &answer, IntArray &activeConditionMap, FloatArray &gamma, GaussPoint *gp, const FloatArray &totalStrain, FloatArray &plasticStrainR, FloatArray &strainSpaceHardeningVariables, TimeStep *tStep) const |
| void | computeGradientVector (FloatArray &answer, functType ftype, int isurf, GaussPoint *gp, const FloatArray &fullStressVector, const FloatArray &fullStressSpaceHardeningVars) const |
| void | computeResidualVector (FloatArray &answer, GaussPoint *gp, const FloatArray &gamma, const IntArray &activeConditionMap, const FloatArray &plasticStrainVectorR, const FloatArray &strainSpaceHardeningVariables, std ::vector< FloatArray > &gradVec) const |
| virtual FloatMatrix | giveConsistentStiffnessMatrix (MatResponseMode, GaussPoint *gp, TimeStep *tStep) const |
| virtual void | giveElastoPlasticStiffnessMatrix (FloatMatrix &answer, MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) const |
| void | computeAlgorithmicModuli (FloatMatrix &answer, GaussPoint *gp, const FloatMatrix &elasticModuliInverse, const FloatMatrix &hardeningModuliInverse, const FloatArray &gamma, const IntArray &activeConditionMap, const FloatArray &fullStressVector, const FloatArray &fullStressSpaceHardeningVars) const |
| void | computeDiagModuli (FloatMatrix &answer, GaussPoint *gp, FloatMatrix &elasticModuliInverse, FloatMatrix &hardeningModuliInverse) const |
| virtual void | computeStressSpaceHardeningVars (FloatArray &answer, GaussPoint *gp, const FloatArray &strainSpaceHardeningVariables) const =0 |
| virtual double | computeYieldValueAt (GaussPoint *gp, int isurf, const FloatArray &stressVector, const FloatArray &stressSpaceHardeningVars) const =0 |
| virtual void | computeHardeningReducedModuli (FloatMatrix &answer, GaussPoint *gp, const FloatArray &strainSpaceHardeningVariables, TimeStep *tStep) const =0 |
| virtual void | computeStressGradientVector (FloatArray &answer, functType ftype, int isurf, GaussPoint *gp, const FloatArray &stressVector, const FloatArray &stressSpaceHardeningVars) const =0 |
| virtual void | computeStressSpaceHardeningVarsReducedGradient (FloatArray &answer, functType ftype, int isurf, GaussPoint *gp, const FloatArray &stressVector, const FloatArray &stressSpaceHardeningVars) const =0 |
| virtual int | hasHardening () const |
| virtual void | computeReducedGradientMatrix (FloatMatrix &answer, int isurf, GaussPoint *gp, const FloatArray &stressVector, const FloatArray &stressSpaceHardeningVars) const =0 |
| virtual void | computeTrialStressIncrement (FloatArray &answer, GaussPoint *gp, const FloatArray &strainIncrement, TimeStep *tStep) const |
| virtual void | computeReducedElasticModuli (FloatMatrix &answer, GaussPoint *gp, TimeStep *tStep) const |
| FloatMatrixF< 3, 3 > | givePlaneStressStiffMtrx (MatResponseMode, GaussPoint *gp, TimeStep *tStep) const override |
| FloatMatrixF< 4, 4 > | givePlaneStrainStiffMtrx (MatResponseMode, GaussPoint *gp, TimeStep *tStep) const override |
| FloatMatrixF< 1, 1 > | give1dStressStiffMtrx (MatResponseMode, GaussPoint *gp, TimeStep *tStep) const override |
| FloatMatrixF< 2, 2 > | give2dBeamLayerStiffMtrx (MatResponseMode, GaussPoint *gp, TimeStep *tStep) const override |
| FloatMatrixF< 5, 5 > | givePlateLayerStiffMtrx (MatResponseMode, GaussPoint *gp, TimeStep *tStep) const override |
| FloatMatrixF< 3, 3 > | giveFiberStiffMtrx (MatResponseMode, GaussPoint *gp, TimeStep *tStep) const override |
Protected Attributes | |
| LinearElasticMaterial * | linearElasticMaterial = nullptr |
| Reference to bulk (undamaged) material. | |
| int | nsurf = 0 |
| Number of yield surfaces. | |
| enum oofem::MPlasticMaterial::ReturnMappingAlgoType | rmType = mpm_ClosestPoint |
| enum oofem::MPlasticMaterial::plastType | plType = associatedPT |
| 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. | |
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 ¤tStress) |
| 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. | |
This class implements a general plastic material. It is assumed to be a base class for many material models based on different yield conditions and hardening laws. Sress return mapping algorithm is based on general return mapping algorithm, with following assumptions
Definition at line 141 of file mplasticmaterial.h.
|
protected |
Type that allows to distinguish between yield function and loading function.
| Enumerator | |
|---|---|
| yieldFunction | |
| loadFunction | |
Definition at line 151 of file mplasticmaterial.h.
|
protected |
| Enumerator | |
|---|---|
| associatedPT | |
| nonassociatedPT | |
Definition at line 152 of file mplasticmaterial.h.
|
protected |
Protected type to determine the return mapping algorithm.
| Enumerator | |
|---|---|
| mpm_ClosestPoint | |
| mpm_CuttingPlane | |
Definition at line 149 of file mplasticmaterial.h.
| oofem::MPlasticMaterial::MPlasticMaterial | ( | int | n, |
| Domain * | d ) |
Definition at line 48 of file mplasticmaterial.C.
References oofem::StructuralMaterial::StructuralMaterial().
Referenced by give1dStressStiffMtrx(), give2dBeamLayerStiffMtrx(), giveFiberStiffMtrx(), givePlateLayerStiffMtrx(), giveRealStressVector_1d(), giveRealStressVector_3d(), giveRealStressVector_PlaneStrain(), giveRealStressVector_PlaneStress(), oofem::J2MPlasticMaterial::J2MPlasticMaterial(), and oofem::RankinePlasticMaterial::RankinePlasticMaterial().
|
virtual |
Definition at line 53 of file mplasticmaterial.C.
References linearElasticMaterial.
|
protected |
Definition at line 162 of file mplasticmaterial.C.
References oofem::FloatArray::add(), associatedPT, oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::IntArray::at(), oofem::FloatMatrix::beInverseOf(), oofem::FloatArray::beProductOf(), oofem::FloatArray::beTProductOf(), oofem::FloatMatrix::clear(), computeAlgorithmicModuli(), computeDiagModuli(), computeGradientVector(), computeHardeningReducedModuli(), oofem::FloatArray::computeNorm(), computeReducedElasticModuli(), computeResidualVector(), computeStressSpaceHardeningVars(), computeTrialStressIncrement(), computeYieldValueAt(), oofem::FloatArray::dotProduct(), oofem::GaussPoint::giveElement(), oofem::FEMComponent::giveNumber(), oofem::GaussPoint::giveNumber(), oofem::FloatMatrix::giveNumberOfRows(), oofem::MPlasticMaterialStatus::givePlasticStrainVector(), oofem::FloatArray::giveSize(), oofem::Material::giveStatus(), oofem::MPlasticMaterialStatus::giveStrainSpaceHardeningVars(), loadFunction, nonassociatedPT, nsurf, OOFEM_WARNING, PLASTIC_MATERIAL_MAX_ITERATIONS, plType, RES_TOL, oofem::FloatArray::resize(), oofem::FloatMatrix::resize(), oofem::FloatMatrix::solveForRhs(), oofem::FloatArray::subtract(), YIELD_TOL, yieldFunction, oofem::FloatArray::zero(), and oofem::IntArray::zero().
Referenced by giveRealStressVector().
|
protected |
Definition at line 751 of file mplasticmaterial.C.
References oofem::FloatMatrix::add(), oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::IntArray::at(), oofem::FloatMatrix::beInverseOf(), computeReducedGradientMatrix(), oofem::FloatMatrix::giveNumberOfRows(), oofem::FloatArray::giveSize(), hasHardening(), nsurf, oofem::FloatMatrix::resize(), oofem::FloatMatrix::times(), and oofem::FloatMatrix::zero().
Referenced by closestPointReturn(), and giveConsistentStiffnessMatrix().
|
protected |
Definition at line 1105 of file mplasticmaterial.C.
References oofem::FloatMatrix::at(), oofem::FloatMatrix::giveNumberOfRows(), oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().
Referenced by closestPointReturn().
|
protected |
Definition at line 652 of file mplasticmaterial.C.
References oofem::FloatArray::at(), computeStressGradientVector(), computeStressSpaceHardeningVarsReducedGradient(), oofem::GaussPoint::giveMaterialMode(), oofem::FloatArray::giveSize(), oofem::FloatArray::isNotEmpty(), and oofem::FloatArray::resize().
Referenced by closestPointReturn(), cuttingPlaneReturn(), giveConsistentStiffnessMatrix(), and giveElastoPlasticStiffnessMatrix().
|
protectedpure virtual |
Implemented in oofem::J2MPlasticMaterial, and oofem::RankinePlasticMaterial.
Referenced by closestPointReturn(), cuttingPlaneReturn(), giveConsistentStiffnessMatrix(), and giveElastoPlasticStiffnessMatrix().
|
protectedvirtual |
Definition at line 1139 of file mplasticmaterial.C.
References linearElasticMaterial.
Referenced by closestPointReturn(), computeTrialStressIncrement(), cuttingPlaneReturn(), giveConsistentStiffnessMatrix(), and giveElastoPlasticStiffnessMatrix().
|
protectedpure virtual |
Implemented in oofem::J2MPlasticMaterial, and oofem::RankinePlasticMaterial.
Referenced by computeAlgorithmicModuli().
|
protected |
Definition at line 697 of file mplasticmaterial.C.
References oofem::FloatArray::at(), oofem::IntArray::at(), oofem::MPlasticMaterialStatus::givePlasticStrainVector(), oofem::FloatArray::giveSize(), oofem::Material::giveStatus(), oofem::MPlasticMaterialStatus::giveStrainSpaceHardeningVars(), nsurf, and oofem::FloatArray::resize().
Referenced by closestPointReturn().
|
protectedpure virtual |
Implemented in oofem::J2MPlasticMaterial, and oofem::RankinePlasticMaterial.
Referenced by computeGradientVector().
|
protectedpure virtual |
Implemented in oofem::J2MPlasticMaterial, and oofem::RankinePlasticMaterial.
Referenced by closestPointReturn(), cuttingPlaneReturn(), giveConsistentStiffnessMatrix(), and giveElastoPlasticStiffnessMatrix().
|
protectedpure virtual |
Implemented in oofem::J2MPlasticMaterial, and oofem::RankinePlasticMaterial.
Referenced by computeGradientVector().
|
protectedvirtual |
Definition at line 734 of file mplasticmaterial.C.
References oofem::FloatArray::beProductOf(), computeReducedElasticModuli(), and oofem::GaussPoint::giveMaterialMode().
Referenced by closestPointReturn(), and cuttingPlaneReturn().
|
protectedpure virtual |
Implemented in oofem::J2MPlasticMaterial, and oofem::RankinePlasticMaterial.
Referenced by closestPointReturn(), and cuttingPlaneReturn().
|
overridevirtual |
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::RankinePlasticMaterial.
Definition at line 75 of file mplasticmaterial.C.
References giveSizeOfReducedHardeningVarsVector().
|
protected |
Definition at line 404 of file mplasticmaterial.C.
References oofem::FloatArray::add(), associatedPT, oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::IntArray::at(), oofem::FloatArray::beTProductOf(), computeGradientVector(), computeHardeningReducedModuli(), computeReducedElasticModuli(), computeStressSpaceHardeningVars(), computeTrialStressIncrement(), computeYieldValueAt(), oofem::GaussPoint::giveElement(), oofem::GaussPoint::giveMaterialMode(), oofem::FEMComponent::giveNumber(), oofem::GaussPoint::giveNumber(), oofem::FloatMatrix::giveNumberOfRows(), oofem::MPlasticMaterialStatus::givePlasticStrainVector(), oofem::Material::giveStatus(), oofem::MPlasticMaterialStatus::giveStrainSpaceHardeningVars(), oofem::StructuralMaterialStatus::giveStressVector(), loadFunction, nonassociatedPT, nsurf, OOFEM_WARNING, PLASTIC_MATERIAL_MAX_ITERATIONS, plType, oofem::FloatArray::resize(), oofem::FloatMatrix::resize(), oofem::FloatMatrix::solveForRhs(), oofem::FloatArray::subtract(), YIELD_TOL, yieldFunction, oofem::FloatArray::zero(), oofem::FloatMatrix::zero(), and oofem::IntArray::zero().
Referenced by giveRealStressVector().
|
overrideprotectedvirtual |
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.
| 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 1238 of file mplasticmaterial.C.
References giveConsistentStiffnessMatrix(), giveElastoPlasticStiffnessMatrix(), linearElasticMaterial, MPlasticMaterial(), mpm_ClosestPoint, and rmType.
|
overrideprotectedvirtual |
Method for computing 2d beam layer stiffness matrix of receiver. Default implementation computes 3d stiffness matrix using give3dMaterialStiffnessMatrix and reduces it to 2d beam layer 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 1259 of file mplasticmaterial.C.
References giveConsistentStiffnessMatrix(), giveElastoPlasticStiffnessMatrix(), linearElasticMaterial, MPlasticMaterial(), mpm_ClosestPoint, and rmType.
|
overridevirtual |
Computes full 3d material stiffness matrix at given integration point, time, respecting load history in integration point.
| 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). |
Reimplemented from oofem::StructuralMaterial.
Definition at line 1150 of file mplasticmaterial.C.
References giveConsistentStiffnessMatrix(), giveElastoPlasticStiffnessMatrix(), linearElasticMaterial, mpm_ClosestPoint, and rmType.
|
inlineoverridevirtual |
Implements oofem::FEMComponent.
Reimplemented in oofem::RankinePlasticMaterial.
Definition at line 160 of file mplasticmaterial.h.
|
protectedvirtual |
Definition at line 823 of file mplasticmaterial.C.
References associatedPT, oofem::FloatMatrix::at(), oofem::IntArray::at(), oofem::FloatMatrix::beInverseOf(), oofem::FloatMatrix::beProductOf(), oofem::FloatMatrix::beProductTOf(), oofem::FloatMatrix::beSubMatrixOf(), oofem::FloatMatrix::clear(), computeAlgorithmicModuli(), computeGradientVector(), computeHardeningReducedModuli(), computeReducedElasticModuli(), computeStressSpaceHardeningVars(), oofem::GaussPoint::giveMaterialMode(), oofem::FloatMatrix::giveNumberOfColumns(), oofem::FloatMatrix::giveNumberOfRows(), oofem::Material::giveStatus(), oofem::MPlasticMaterialStatus::giveStrainSpaceHardeningVars(), oofem::StructuralMaterialStatus::giveStressVector(), oofem::MPlasticMaterialStatus::giveTempActiveConditionMap(), oofem::MPlasticMaterialStatus::giveTempGamma(), oofem::MPlasticMaterialStatus::giveTempStateFlag(), linearElasticMaterial, loadFunction, nonassociatedPT, nsurf, plType, oofem::FloatMatrix::resize(), oofem::FloatMatrix::subtract(), and yieldFunction.
Referenced by give1dStressStiffMtrx(), give2dBeamLayerStiffMtrx(), give3dMaterialStiffnessMatrix(), giveFiberStiffMtrx(), givePlaneStrainStiffMtrx(), givePlaneStressStiffMtrx(), and givePlateLayerStiffMtrx().
|
protectedvirtual |
Definition at line 973 of file mplasticmaterial.C.
References oofem::FloatMatrix::add(), associatedPT, oofem::FloatMatrix::at(), oofem::IntArray::at(), oofem::FloatMatrix::beInverseOf(), oofem::FloatMatrix::beProductOf(), oofem::FloatMatrix::beProductTOf(), oofem::FloatMatrix::beTProductOf(), computeGradientVector(), computeHardeningReducedModuli(), computeReducedElasticModuli(), computeStressSpaceHardeningVars(), oofem::GaussPoint::giveMaterialMode(), oofem::FloatMatrix::giveNumberOfRows(), oofem::Material::giveStatus(), oofem::MPlasticMaterialStatus::giveStrainSpaceHardeningVars(), oofem::StructuralMaterialStatus::giveStressVector(), oofem::MPlasticMaterialStatus::giveTempActiveConditionMap(), oofem::MPlasticMaterialStatus::giveTempStateFlag(), linearElasticMaterial, loadFunction, oofem::FloatMatrix::negated(), nonassociatedPT, nsurf, plType, oofem::FloatMatrix::resize(), yieldFunction, and oofem::FloatMatrix::zero().
Referenced by give1dStressStiffMtrx(), give2dBeamLayerStiffMtrx(), give3dMaterialStiffnessMatrix(), giveFiberStiffMtrx(), givePlaneStrainStiffMtrx(), givePlaneStressStiffMtrx(), and givePlateLayerStiffMtrx().
|
overrideprotectedvirtual |
Method for computing 1d fiber stiffness matrix of receiver. Default implementation computes 3d stiffness matrix using give3dMaterialStiffnessMatrix and reduces it to 1d fiber 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 1307 of file mplasticmaterial.C.
References giveConsistentStiffnessMatrix(), giveElastoPlasticStiffnessMatrix(), linearElasticMaterial, MPlasticMaterial(), mpm_ClosestPoint, and rmType.
|
overridevirtual |
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::Material.
Definition at line 1331 of file mplasticmaterial.C.
References oofem::StructuralMaterial::computePrincipalValues(), oofem::GaussPoint::giveMaterialMode(), oofem::MPlasticMaterialStatus::givePlasticStrainVector(), oofem::Material::giveStatus(), and oofem::principal_strain.
|
inline |
Returns reference to undamaged (bulk) material.
Definition at line 163 of file mplasticmaterial.h.
References linearElasticMaterial.
|
overrideprotectedvirtual |
Method for computing plane strain stiffness matrix of receiver. Default implementation computes 3d stiffness matrix using give3dMaterialStiffnessMatrix and reduces it to plane strain 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. Note: as already described, if zero strain component is imposed (Plane strain, ..) this condition must be taken into account in geometrical relations, and corresponding component has to be included in reduced vector. (So plane strain conditions are \( \epsilon_z = \gamma_{xz} = \gamma_{yz} = 0 \), but relations for \( \epsilon_z\) and \(\sigma_z\) are included).
| 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 1215 of file mplasticmaterial.C.
References giveConsistentStiffnessMatrix(), giveElastoPlasticStiffnessMatrix(), linearElasticMaterial, mpm_ClosestPoint, and rmType.
|
overrideprotectedvirtual |
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.
Definition at line 1190 of file mplasticmaterial.C.
References giveConsistentStiffnessMatrix(), giveElastoPlasticStiffnessMatrix(), linearElasticMaterial, mpm_ClosestPoint, and rmType.
|
overrideprotectedvirtual |
Method for computing 2d plate layer stiffness matrix of receiver. Default implementation computes 3d stiffness matrix using give3dMaterialStiffnessMatrix and reduces it to 2d plate layer 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 1283 of file mplasticmaterial.C.
References giveConsistentStiffnessMatrix(), giveElastoPlasticStiffnessMatrix(), linearElasticMaterial, MPlasticMaterial(), mpm_ClosestPoint, and rmType.
|
overridevirtual |
Computes the real stress vector for given total strain and integration point. The total strain is defined as strain computed directly from displacement field at given time. The stress independent parts (temperature, eigenstrains) are subtracted in constitutive driver. The service should use previously reached equilibrium history variables. Also it should update temporary history variables in status according to newly reached state. The temporary history variables are moved into equilibrium ones after global structure equilibrium has been reached by iteration process.
| answer | Stress vector in reduced form. For large deformations it is treated as the second Piola-Kirchoff stress. |
| gp | Integration point. |
| reducedStrain | Strain vector in reduced form. For large deformations it is treated as the Green-Lagrange strain. |
| tStep | Current time step (most models are able to respond only when tStep is current time step). |
Reimplemented from oofem::StructuralMaterial.
Definition at line 82 of file mplasticmaterial.C.
References oofem::FloatArray::at(), closestPointReturn(), cuttingPlaneReturn(), oofem::MPlasticMaterialStatus::giveStateFlag(), oofem::Material::giveStatus(), oofem::StructuralMaterial::giveStressDependentPartOfStrainVector(), oofem::StructuralMaterialStatus::giveTempStressVector(), oofem::Material::initTempStatus(), oofem::MPlasticMaterialStatus::letTempPlasticStrainVectorBe(), oofem::MPlasticMaterialStatus::letTempStateFlagBe(), oofem::MPlasticMaterialStatus::letTempStrainSpaceHardeningVarsVectorBe(), oofem::StructuralMaterialStatus::letTempStrainVectorBe(), oofem::StructuralMaterialStatus::letTempStressVectorBe(), mpm_ClosestPoint, nsurf, rmType, oofem::MPlasticMaterialStatus::setTempActiveConditionMap(), and oofem::MPlasticMaterialStatus::setTempGamma().
Referenced by giveRealStressVector_1d(), giveRealStressVector_3d(), giveRealStressVector_PlaneStrain(), and giveRealStressVector_PlaneStress().
|
inlineoverridevirtual |
Default implementation relies on giveRealStressVector_StressControl.
Reimplemented from oofem::StructuralMaterial.
Definition at line 195 of file mplasticmaterial.h.
References giveRealStressVector(), and MPlasticMaterial().
|
inlineoverridevirtual |
Default implementation relies on giveRealStressVector for second Piola-Kirchoff stress.
Reimplemented from oofem::StructuralMaterial.
Definition at line 177 of file mplasticmaterial.h.
References giveRealStressVector(), and MPlasticMaterial().
|
inlineoverridevirtual |
Default implementation relies on giveRealStressVector_3d.
Reimplemented from oofem::StructuralMaterial.
Definition at line 183 of file mplasticmaterial.h.
References giveRealStressVector(), and MPlasticMaterial().
|
inlineoverridevirtual |
Default implementation relies on giveRealStressVector_StressControl.
Reimplemented from oofem::StructuralMaterial.
Definition at line 189 of file mplasticmaterial.h.
References giveRealStressVector(), and MPlasticMaterial().
|
inlinevirtual |
Reimplemented in oofem::J2MPlasticMaterial.
Definition at line 205 of file mplasticmaterial.h.
|
inlinevirtual |
Reimplemented in oofem::J2MPlasticMaterial.
Definition at line 206 of file mplasticmaterial.h.
Referenced by CreateStatus(), and oofem::RankinePlasticMaterial::CreateStatus().
|
inlineprotectedvirtual |
Reimplemented in oofem::J2MPlasticMaterial, and oofem::RankinePlasticMaterial.
Definition at line 256 of file mplasticmaterial.h.
Referenced by computeAlgorithmicModuli().
|
overridevirtual |
Tests if material supports material mode.
| mode | Required material mode. |
Reimplemented from oofem::Material.
Definition at line 60 of file mplasticmaterial.C.
|
inlineoverridevirtual |
Returns true if stiffness matrix of receiver is symmetric. Default implementation returns true.
Reimplemented from oofem::Material.
Definition at line 169 of file mplasticmaterial.h.
|
protected |
Reference to bulk (undamaged) material.
Definition at line 145 of file mplasticmaterial.h.
Referenced by oofem::J2MPlasticMaterial::compute3dElasticModuli(), computeReducedElasticModuli(), give1dStressStiffMtrx(), give2dBeamLayerStiffMtrx(), give3dMaterialStiffnessMatrix(), giveConsistentStiffnessMatrix(), giveElastoPlasticStiffnessMatrix(), giveFiberStiffMtrx(), giveLinearElasticMaterial(), givePlaneStrainStiffMtrx(), givePlaneStressStiffMtrx(), givePlateLayerStiffMtrx(), oofem::J2MPlasticMaterial::initializeFrom(), oofem::RankinePlasticMaterial::initializeFrom(), oofem::J2MPlasticMaterial::J2MPlasticMaterial(), oofem::RankinePlasticMaterial::RankinePlasticMaterial(), and ~MPlasticMaterial().
|
protected |
Number of yield surfaces.
Definition at line 147 of file mplasticmaterial.h.
Referenced by closestPointReturn(), computeAlgorithmicModuli(), computeResidualVector(), cuttingPlaneReturn(), giveConsistentStiffnessMatrix(), giveElastoPlasticStiffnessMatrix(), giveRealStressVector(), oofem::J2MPlasticMaterial::J2MPlasticMaterial(), and oofem::RankinePlasticMaterial::RankinePlasticMaterial().
|
protected |
|
protected |
Referenced by give1dStressStiffMtrx(), give2dBeamLayerStiffMtrx(), give3dMaterialStiffnessMatrix(), giveFiberStiffMtrx(), givePlaneStrainStiffMtrx(), givePlaneStressStiffMtrx(), givePlateLayerStiffMtrx(), giveRealStressVector(), oofem::J2MPlasticMaterial::initializeFrom(), and oofem::RankinePlasticMaterial::RankinePlasticMaterial().