SnowConstitutiveModel.hpp

Go to the documentation of this file.

// This code is based on Jet framework.
// Copyright (c) 2018 Doyub Kim
// CubbyFlow is voxel-based fluid simulation engine for computer games.
// Copyright (c) 2020 CubbyFlow Team
// Core Part: Chris Ohk, Junwoo Hwang, Jihong Sin, Seungwoo Yoo
// AI Part: Dongheon Cho, Minseo Kim
// We are making my contributions/submissions to this project solely in our
// personal capacity and are not conveying any rights to any intellectual
// property of any third parties.
 
#ifndef CUBBYFLOW_SNOW_CONSTITUTIVE_MODEL_HPP
#define CUBBYFLOW_SNOW_CONSTITUTIVE_MODEL_HPP
 
#include <Core/Matrix/Matrix.hpp>
 
namespace CubbyFlow
{
template <size_t N>
struct SnowDeformationState
{
    using MatrixType = Matrix<double, N, N>;
 
    MatrixType elastic = MatrixType::MakeIdentity();
    MatrixType plastic = MatrixType::MakeIdentity();
};
 
template <size_t N>
class SnowConstitutiveModel final
{
 public:
    static_assert(N == 2 || N == 3, "Snow model supports only 2-D and 3-D.");
 
    using MatrixType = Matrix<double, N, N>;
    using State = SnowDeformationState<N>;
 
    explicit SnowConstitutiveModel(double youngsModulus = 1.4e5,
                                   double poissonRatio = 0.2,
                                   double criticalCompression = 2.5e-2,
                                   double criticalStretch = 7.5e-3,
                                   double hardeningCoefficient = 10.0);
 
    [[nodiscard]] State Update(const MatrixType& deformationGradientIncrement,
                               const State& state) const;
 
    [[nodiscard]] MatrixType ComputeKirchhoffStress(const State& state) const;
 
    [[nodiscard]] MatrixType ComputeFirstPiolaStressDifferential(
        const State& state, const MatrixType& differential) const;
 
    [[nodiscard]] double ComputeWaveSpeed(const State& state,
                                          double referenceDensity) const;
 
 private:
    [[nodiscard]] double ComputeHardening(const State& state) const;
 
    [[nodiscard]] static bool IsFinite(const MatrixType& matrix);
 
    static void ValidateDeformation(const MatrixType& matrix);
 
    double m_mu0 = 0.0;
    double m_lambda0 = 0.0;
    double m_criticalCompression;
    double m_criticalStretch;
    double m_hardeningCoefficient;
};
 
using SnowDeformationState2 = SnowDeformationState<2>;
using SnowDeformationState3 = SnowDeformationState<3>;
using SnowConstitutiveModel2 = SnowConstitutiveModel<2>;
using SnowConstitutiveModel3 = SnowConstitutiveModel<3>;
}  // namespace CubbyFlow
 
#include <Core/Particle/MPM/SnowConstitutiveModel-Impl.hpp>
 
#endif