idm1.h

Go to the documentation of this file.

/*
 *
 *                 #####    #####   ######  ######  ###   ###
 *               ##   ##  ##   ##  ##      ##      ## ### ##
 *              ##   ##  ##   ##  ####    ####    ##  #  ##
 *             ##   ##  ##   ##  ##      ##      ##     ##
 *            ##   ##  ##   ##  ##      ##      ##     ##
 *            #####    #####   ##      ######  ##     ##
 *
 *
 *             OOFEM : Object Oriented Finite Element Code
 *
 *               Copyright (C) 1993 - 2025   Borek Patzak
 *
 *
 *
 *       Czech Technical University, Faculty of Civil Engineering,
 *   Department of Structural Mechanics, 166 29 Prague, Czech Republic
 *
 *  This library is free software; you can redistribute it and/or
 *  modify it under the terms of the GNU Lesser General Public
 *  License as published by the Free Software Foundation; either
 *  version 2.1 of the License, or (at your option) any later version.
 *
 *  This program is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 *  Lesser General Public License for more details.
 *
 *  You should have received a copy of the GNU Lesser General Public
 *  License along with this library; if not, write to the Free Software
 *  Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
 */
 
#ifndef idm1_h
#define idm1_h

//#define IDM_USE_MMAClosestIPTransfer
#define IDM_USE_MMAContainingElementProjection
//#define IDM_USE_MMAShapeFunctProjection
//#define IDM_USE_MMALeastSquareProjection
 
/*
 * Selects the use of mapped strain or projected strain from element.
 */
#define IDM_USE_MAPPEDSTRAIN
 
#include "material.h"
#include "sm/Materials/linearelasticmaterial.h"
#include "sm/Materials/isodamagemodel.h"
#include "sm/Materials/structuralms.h"
#include "randommaterialext.h"
#include "materialmapperinterface.h"
 
#ifdef IDM_USE_MMAClosestIPTransfer
 #include "mmaclosestiptransfer.h"
#endif
 
#ifdef IDM_USE_MMAContainingElementProjection
 #include "mmacontainingelementprojection.h"
#endif
 
#ifdef IDM_USE_MMAShapeFunctProjection
 #include "mmashapefunctprojection.h"
#endif
 
#ifdef IDM_USE_MMALeastSquareProjection
 #include "mmaleastsquareprojection.h"
#endif


#define _IFT_IsotropicDamageMaterial1_Name "idm1"
#define _IFT_IsotropicDamageMaterial1_e0 "e0"
#define _IFT_IsotropicDamageMaterial1_ef "ef"
#define _IFT_IsotropicDamageMaterial1_wf "wf"
#define _IFT_IsotropicDamageMaterial1_equivstraintype "equivstraintype"
#define _IFT_IsotropicDamageMaterial1_damageLaw "damlaw"
#define _IFT_IsotropicDamageMaterial1_k "k"
#define _IFT_IsotropicDamageMaterial1_md "md"
#define _IFT_IsotropicDamageMaterial1_ecsm "ecsm"
#define _IFT_IsotropicDamageMaterial1_At "at"
#define _IFT_IsotropicDamageMaterial1_Bt "bt"
#define _IFT_IsotropicDamageMaterial1_ft "ft"
#define _IFT_IsotropicDamageMaterial1_wkwf "wkwf"
#define _IFT_IsotropicDamageMaterial1_e1ef "e1ef"
#define _IFT_IsotropicDamageMaterial1_skft "skft"
#define _IFT_IsotropicDamageMaterial1_s1 "s1"
#define _IFT_IsotropicDamageMaterial1_sk "sk"
#define _IFT_IsotropicDamageMaterial1_wk "wk"
#define _IFT_IsotropicDamageMaterial1_e1 "e1"
#define _IFT_IsotropicDamageMaterial1_ek "ek"
#define _IFT_IsotropicDamageMaterial1_gf "gf"
#define _IFT_IsotropicDamageMaterial1_gft "gft"
#define _IFT_IsotropicDamageMaterial1_ep "ep"
#define _IFT_IsotropicDamageMaterial1_e2 "e2"
#define _IFT_IsotropicDamageMaterial1_nd "nd"
#define _IFT_IsotropicDamageMaterial1_checkSnapBack "checksnapback"
#define _IFT_IsotropicDamageMaterial1_n "griff_n"
#define _IFT_IsotropicDamageMaterial1_c1 "c1"
#define _IFT_IsotropicDamageMaterial1_c2 "c2"
#define _IFT_IsotropicDamageMaterial1_alphaps "alphaps"
#define _IFT_IsotropicDamageMaterial1_h "h"
#define _IFT_IsotropicDamageMaterial1_w_k "w_k"
#define _IFT_IsotropicDamageMaterial1_w_r "w_r"
#define _IFT_IsotropicDamageMaterial1_w_f "w_f"
#define _IFT_IsotropicDamageMaterial1_f_k "f_k"
#define _IFT_IsotropicDamageMaterial1_f_r "f_r"
 
namespace oofem {
#define IDM1_ITERATION_LIMIT 1.e-9

class IsotropicDamageMaterial1Status : public IsotropicDamageMaterialStatus, public RandomMaterialStatusExtensionInterface
{
public:
    IsotropicDamageMaterial1Status(GaussPoint *g);
 
    const char *giveClassName() const override { return "IsotropicDamageMaterial1Status"; }
 
    Interface *giveInterface(InterfaceType it) override;
};

class IsotropicDamageMaterial1 : public IsotropicDamageMaterial,
    public RandomMaterialExtensionInterface,
    public MaterialModelMapperInterface
{
protected:
    double e0 = 0.;
    double ef = 0.;
    double wf = 0.;

    double gf = 0.;

    double gft = 0.;

    double ek = 0.;

    double wk = 0.;

    double sk = 0.;

    double c1 = 3., c2 = 6.93;  // default value of Hordijk parameter

    double w_k = 0., w_r = 0., w_f = 0., f_k = 0., f_r = 0.;

    enum EquivStrainType {
        EST_Mazars=0,
        EST_Rankine_Smooth=1,
        EST_ElasticEnergy=2,
        EST_Mises=3,
        EST_Rankine_Standard=4,
        EST_ElasticEnergyPositiveStress=5,
        EST_ElasticEnergyPositiveStrain=6,
        EST_Griffith=7,
        EST_Unknown = 100
    };

    EquivStrainType equivStrainType = EST_Unknown;

    double k = 0.;

    double griff_n = 8.;

    int damageLaw = 0;

    enum SofteningType { ST_Unknown, ST_Exponential, ST_Linear, ST_Mazars, ST_Smooth, ST_SmoothExtended, ST_Exponential_Cohesive_Crack, ST_Linear_Cohesive_Crack, ST_BiLinear_Cohesive_Crack, ST_Disable_Damage, ST_PowerExponential, ST_DoubleExponential, ST_Hordijk_Cohesive_Crack, ST_ModPowerExponential, ST_Trilinear_Cohesive_Crack };

    SofteningType softType = ST_Unknown;

    double At = 0., Bt = 0.;
    double md = 1.;

    double e1 = 0., e2 = 0., s1 = 0., nd = 0.;
    int checkSnapBack = 0;

    double ep = 0., ft = 0.;

    double ps_alpha = 0., ps_H = 0.;

    ElementCharSizeMethod ecsMethod = ECSM_Unknown;

    Set *sourceElemSet = nullptr;
 
#ifdef IDM_USE_MMAClosestIPTransfer
    static MMAClosestIPTransfer mapper;
#endif
#ifdef IDM_USE_MMAContainingElementProjection
    static MMAContainingElementProjection mapper;
#endif
#ifdef IDM_USE_MMAShapeFunctProjection
    static MMAShapeFunctProjection mapper;
#endif
#ifdef IDM_USE_MMALeastSquareProjection
    static MMALeastSquareProjection mapper;
#endif
 
public:
    IsotropicDamageMaterial1(int n, Domain *d);
    virtual ~IsotropicDamageMaterial1();
 
    const char *giveClassName() const override { return "IsotropicDamageMaterial1"; }
    const char *giveInputRecordName() const override { return _IFT_IsotropicDamageMaterial1_Name; }
    void initializeFrom(InputRecord &ir) override;
    void giveInputRecord(DynamicInputRecord &input) override;
    static void computeStrainInvariants(const FloatArray &strainVector, double &I1e, double &J2e);
 
    bool isCrackBandApproachUsed() const { return ( this->softType == ST_Exponential_Cohesive_Crack || this->softType == ST_Linear_Cohesive_Crack || this->softType == ST_BiLinear_Cohesive_Crack || this->softType == ST_Trilinear_Cohesive_Crack || this->gf != 0. ); }
    double computeEquivalentStrain(const FloatArray &strain, GaussPoint *gp, TimeStep *tStep) const override;
 
    void computeEta(FloatArray &answer, const FloatArray &strain, GaussPoint *gp, TimeStep *tStep) const override;
    double computeDamageParam(double kappa, const FloatArray &strain, GaussPoint *gp) const override;
    double computeDamageParamForCohesiveCrack(double kappa, GaussPoint *gp) const;
    double damageFunction(double kappa, GaussPoint *gp) const;
    double damageFunctionPrime(double kappa, GaussPoint *gp) const override;
    double complianceFunction(double kappa, GaussPoint *gp) const;
 
    double evaluatePermanentStrain(double kappa, double omega) const override;
 
    Interface *giveInterface(InterfaceType it) override;
 
    int MMI_map(GaussPoint *gp, Domain *oldd, TimeStep *tStep) override;
    int MMI_update(GaussPoint *gp, TimeStep *tStep, FloatArray *estrain = nullptr) override;
    int MMI_finish(TimeStep *tStep) override;
 
    std::unique_ptr<MaterialStatus> CreateStatus(GaussPoint *gp) const override;
    MaterialStatus *giveStatus(GaussPoint *gp) const override;
 
    double give(int aProperty, GaussPoint *gp) const override;
 
    bool isCharacteristicMtrxSymmetric(MatResponseMode rMode) const override { return false; }
 
    void restoreContext(DataStream &stream, ContextMode mode) override;
    void saveContext(DataStream &stream, ContextMode mode) override;
protected:
    void initDamaged(double kappa, FloatArray &totalStrainVector, GaussPoint *gp) const override;
};
} // end namespace oofem
#endif // idm1_h