hydratingconcretemat.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 hydratingconcretemat_h
#define hydratingconcretemat_h
 
#include "tm/Materials/isoheatmat.h"
#include "tm/Materials/hydratingisoheatmat.h"


#define _IFT_HydratingConcreteMat_Name "hydratingconcretemat"
#define _IFT_HydratingConcreteMat_referenceTemperature "referencetemperature"
#define _IFT_HydratingConcreteMat_relMatAge "relmatage"
#define _IFT_HydratingConcreteMat_hydrationModelType "hydrationmodeltype"
#define _IFT_HydratingConcreteMat_maxModelIntegrationTime "maxmodelintegrationtime"
#define _IFT_HydratingConcreteMat_minModelTimeStepIntegrations "minmodeltimestepintegrations"
#define _IFT_HydratingConcreteMat_conductivitytype "conductivitytype"
#define _IFT_HydratingConcreteMat_capacitytype "capacitytype"
#define _IFT_HydratingConcreteMat_densitytype "densitytype"
#define _IFT_HydratingConcreteMat_activationEnergy "activationenergy"
#define _IFT_HydratingConcreteMat_massCement "masscement"
#define _IFT_HydratingConcreteMat_reinforcementDegree "reinforcementdegree"
#define _IFT_HydratingConcreteMat_tau "tau"
#define _IFT_HydratingConcreteMat_beta "beta"
#define _IFT_HydratingConcreteMat_B1 "b1"
#define _IFT_HydratingConcreteMat_B2 "b2"
#define _IFT_HydratingConcreteMat_eta "eta"
#define _IFT_HydratingConcreteMat_DoHInf "dohinf"
#define _IFT_HydratingConcreteMat_DoH1 "doh1"
#define _IFT_HydratingConcreteMat_P1 "p1"
#define _IFT_HydratingConcreteMat_qpot "qpot"
#define _IFT_HydratingConcreteMat_wc "w/c"
#define _IFT_HydratingConcreteMat_ac "a/c"
#define _IFT_HydratingConcreteMat_rhoCem "rhocem"
#define _IFT_HydratingConcreteMat_rhoAgg "rhoagg"
#define _IFT_HydratingConcreteMat_Blaine "blaine"
#define _IFT_HydratingConcreteMat_alphaSet0 "alphaset0"
#define _IFT_HydratingConcreteMat_timeSet "timeset"
#define _IFT_HydratingConcreteMat_alphaCrit0 "alphacrit0"
#define _IFT_HydratingConcreteMat_B0 "b0"
#define _IFT_HydratingConcreteMat_timeToSeconds "timetoseconds"
 

 
namespace oofem {
class HydratingConcreteMat : public IsotropicHeatTransferMaterial
{
public:
    HydratingConcreteMat(int n, Domain * d);
 
    bool hasInternalSource() const override { return true; }
    bool hasCastingTimeSupport() const override { return true; }
    void computeInternalSourceVector(FloatArray &val, GaussPoint *gp, TimeStep *tStep, ValueModeType mode) const override;
 
    double giveCharacteristicValue(MatResponseMode mode,
                                   GaussPoint *gp,
                                   TimeStep *tStep) const override;
 
    const char *giveClassName() const override { return "HydratingConcreteMat"; }
 
    void initializeFrom(InputRecord &ir) override;
 
    // post-processing
    int giveIPValue(FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep) override;
    double giveIsotropicConductivity(GaussPoint *gp, TimeStep *tStep) const override;
    virtual double giveConcreteCapacity(GaussPoint *gp, TimeStep *tStep) const;
    virtual double giveConcreteDensity(GaussPoint *gp, TimeStep *tStep) const;

    int hydrationModelType = 0;
    double maxModelIntegrationTime = 0.;
    double minModelTimeStepIntegrations = 0.;
    double Qpot = 0.;
    double massCement = 0.;
    double activationEnergy = 0.;
    double referenceTemperature = 0.;
    double relMatAge = 0.;
    
    double tau = 0., beta = 0.;

    double B1 = 0., B2 = 0., eta = 0., DoHInf = 0.;
 
    double DoH1 = 0., P1=0.;
    double wc, ac;
    double rhoCem, rhoAgg;
    double Vc0, Vw0;
    double VCemSet, VCHSet, VGelSet;
    double a0;
    double ng;
    double alphaSet;
    double timeSet;
    double aSet;
    double zSet;
    double alphaCrit;
    double B0;
    double timeToSeconds;
 
protected:
    double GivePower(TimeStep *tStep, GaussPoint *gp, ValueModeType mode) const;
    double scaleTemperature(GaussPoint *gp) const;
    double affinity25(double alpha) const;

    int conductivityType, capacityType, densityType;
    double reinforcementDegree;
    std::unique_ptr<MaterialStatus> CreateStatus(GaussPoint *gp) const override;
};

class HydratingConcreteMatStatus : public TransportMaterialStatus
{
public:
    HydratingConcreteMatStatus(GaussPoint * g);
    double giveDoHActual() const;
    void updateYourself(TimeStep *tStep) override;
    void printOutputAt(FILE *file, TimeStep *tStep) const override;
    double power = 0.;
    double lastEvalTime = -1.e20;
    double equivalentTime = 0.;
    double lastEquivalentTime = 0., degreeOfHydration = 0., lastDegreeOfHydration = 0.;
    double zShell = 0., lastZShell = 0.;
    // Radius of cement particle
    double aCement = 0., lastACement = 0.;
    // Volume fractions of cement, gel, CH;
    double VCem = 0., lastVCem = 0., VGel = 0., lastVGel = 0., VCH = 0., lastVCH = 0.;
    // When the hydration model is called the first time - imposing relMatAge;
    bool firstCall = true;
};
} // end namespace oofem
#endif // hydratingconcretemat_h