materialUtil.h

Go to the documentation of this file.

/*
 * -------------------------------------------
 * Copyright (c) 2021 - 2026 Prashant K. Jha
 * -------------------------------------------
 * PeriDEM https://github.com/prashjha/PeriDEM
 *
 * Distributed under the Boost Software License, Version 1.0. (See accompanying
 * file LICENSE)
 */
 
#ifndef MATERIAL_UTIL_H
#define MATERIAL_UTIL_H
 
#include "util/point.h"
#include "mparticle/material.h"
#include "fracture/fracture.h"
#include "data/modelData.h"
#include <limits>
#include <string>
#include <vector>
 
namespace material {
 
  inline double toNu(double lambda, double mu) { return lambda * 0.5 / (lambda + mu); }
 
  inline double toNuEG(double E, double G) { return E * 0.5 / G - 1.; }
 
  inline double toE(double K, double nu) { return K * (3. * (1. - 2. * nu)); }
 
  inline double toK(double E, double nu) { return E / (3. * (1. - 2. * nu)); }
 
  inline double toLambdaE(double E, double nu) {
    return E * nu / ((1. + nu) * (1. - 2. * nu));
  }
 
  inline double toLambdaK(double K, double nu) { return 3. * K * nu / (1. + nu); }
 
  inline double toGE(double E, double nu) { return E / (2. * (1. + nu)); }
 
  inline double toGK(double K, double nu) {
    return 3. * K * (1. - 2. * nu) / (2. * (1. + nu));
  }
 
  inline double toELambda(double lambda, double nu) {
    return lambda * (1. + nu) * (1. - 2. * nu) / nu;
  }
 
  inline double toGc(double KIc, double nu, double E) { return KIc * KIc / E; }
 
  inline double toKIc(double Gc, double nu, double E) { return std::sqrt(Gc * E); }
 
void computeStateMx(data::ModelData *model, bool compute_in_parallel = false);
 
void computeStateThetax(data::ModelData *model, bool compute_in_parallel = false);
 
 
void computeHydrostaticStrain(data::ModelData *model, bool compute_in_parallel = false);
 
void updateBondFractureData(data::ModelData *model, bool compute_in_parallel = false);
 
 
} // namespace material
 
#endif // MATERIAL_UTIL_H