MercuryOS.h

Go to the documentation of this file.

// This file is part of the MercuryDPM project (https://www.mercurydpm.org).
// Copyright (c), The MercuryDPM Developers Team. All rights reserved.
// License: BSD 3-Clause License; see the LICENSE file in the root directory.
 
#ifndef MERCURYDPM_MERCURYOS_H
#define MERCURYDPM_MERCURYOS_H
 
#include <Mercury3D.h>
#include <Species/HertzianViscoelasticMindlinSpecies.h>
#include <Logger.h>
 
class MercuryOS : public Mercury3D
{
    // if true, additional output files are written
    bool writeOutput_ = false;
    // if true, the simulation time is shortened
    bool test_ = false;
    // if true, smooth walls are used instead of the stl walls
    bool useMercuryWalls_ = false;
    // if true, the simulation time is shortened
    bool soft_ = false;
    
public:
    
    // sets the variable writeOutput_
    void writeOutput(bool writeOutput)
    {
        writeOutput_ = writeOutput;
    }
    
    // returns the variable writeOutput_
    bool writeOutput() const
    {
        return writeOutput_;
    }
 
    // sets the variable test_
    void test(bool test)
    {
        test_ = test;
    }
    
    // returns the variable test_
    bool test() const
    {
        return test_;
    }
    
    // sets the variable soft_
    void soft(bool soft)
    {
        soft_ = soft;
    }
    
    // returns the variable soft_
    bool soft() const
    {
        return soft_;
    }
    
    // sets the variable useMercuryWalls_
    void useMercuryWalls(bool useMercuryWalls)
    {
        useMercuryWalls_ = useMercuryWalls;
    }
    
    // returns the variable useMercuryWalls_
    bool useMercuryWalls() const
    {
        return useMercuryWalls_;
    }
 
private:
    
    static double getEffectiveElasticModulus(double E1, double v1, double E2, double v2)
    {
        return 1. / ((1 - v1 * v1) / E1 + (1 - v2 * v2) / E2);
    }
    
    static double getEffectiveShearModulus(double E1, double v1, double E2, double v2)
    {
        double G1 = E1 / 2 / (1 + v1), G2 = E2 / 2 / (1 + v2);
        return 1. / ((2 - v1) / G1 + (2 - v2) / G2);
    }
    
protected:
    // pointers to the species
    HertzianViscoelasticMindlinSpecies *m1, *m2, *steel;
 
public:
    void setMaterialProperties()
    {
        // Young's modulus and Poisson ratios of M1, M2, steel
        double E1 = 1e9, E2 = 0.5e9, ES = 210e9;
        double v1 = 0.2, v2 = 0.2, vS = 0.2;
        
        if (soft()) {
            E1 = 2e6;
            E2 = 1e6;
        }
        
        // effective elastic and shear moduli
        double E11 = getEffectiveElasticModulus(E1, v1, E1, v1);
        double E12 = getEffectiveElasticModulus(E1, v1, E2, v2);
        double E1S = getEffectiveElasticModulus(E1, v1, ES, vS);
        double E22 = getEffectiveElasticModulus(E2, v2, E2, v2);
        double E2S = getEffectiveElasticModulus(E2, v2, ES, vS);
        double ESS = getEffectiveElasticModulus(ES, vS, ES, vS);
        double G11 = getEffectiveShearModulus(E1, v1, E1, v1);
        double G12 = getEffectiveShearModulus(E1, v1, E2, v2);
        double G1S = getEffectiveShearModulus(E1, v1, ES, vS);
        double G22 = getEffectiveShearModulus(E2, v2, E2, v2);
        double G2S = getEffectiveShearModulus(E2, v2, ES, vS);
        double GSS = getEffectiveShearModulus(ES, vS, ES, vS);
        
        // restitution coefficients
        double r11 = 0.5, r12 = 0.45, r1S = 0.4, r22 = 0.4, r2S = 0.4, rSS = 0.6;
        // sliding friction coefficients
        double mu11 = 0.3, mu12 = 0.2, mu1S = 0.2, mu22 = 0.4, mu2S = 0.2, muSS = 0.5;
        
        HertzianViscoelasticMindlinSpecies species;
        species.setEffectiveElasticModulusAndRestitutionCoefficient(E11, r11);
        species.setEffectiveShearModulus(G11);
        species.setSlidingFrictionCoefficient(mu11);
        species.setDensity(2500);
        m1 = speciesHandler.copyAndAddObject(species);
    
        species.setEffectiveElasticModulusAndRestitutionCoefficient(E22, r22);
        species.setEffectiveShearModulus(G22);
        species.setSlidingFrictionCoefficient(mu22);
        species.setDensity(2000);
        m2 = speciesHandler.copyAndAddObject(species);
    
        species.setEffectiveElasticModulusAndRestitutionCoefficient(ESS, rSS);
        species.setEffectiveShearModulus(GSS);
        species.setSlidingFrictionCoefficient(muSS);
        species.setDensity(7200);
        steel = speciesHandler.copyAndAddObject(species);
        
        auto m12 = speciesHandler.getMixedObject(m1, m2);
        m12->setEffectiveElasticModulusAndRestitutionCoefficient(E12, r12);
        m12->setEffectiveShearModulus(G12);
        m12->setSlidingFrictionCoefficient(mu12);
        
        auto m1S = speciesHandler.getMixedObject(m1, steel);
        m1S->setEffectiveElasticModulusAndRestitutionCoefficient(E1S, r1S);
        m1S->setEffectiveShearModulus(G1S);
        m1S->setSlidingFrictionCoefficient(mu1S);
        
        auto m2S = speciesHandler.getMixedObject(m2, steel);
        m2S->setEffectiveElasticModulusAndRestitutionCoefficient(E2S, r2S);
        m2S->setEffectiveShearModulus(G2S);
        m2S->setSlidingFrictionCoefficient(mu2S);
    }
    
    
};
 
 
#endif //MERCURYDPM_MERCURYOS_H