material Namespace Reference

Data Structures
class	BaseInfluenceFn
	A base class for computing influence function. More...
class	ConstInfluenceFn
	A class to implement constant influence function. More...
class	GaussianInfluenceFn
	A class to implement Gaussian influence function. More...
class	LinearInfluenceFn
	A class to implement linear influence function. More...
class	Material
	Collection of methods and database related to peridynamic material. More...
class	PdElastic
	A class providing methods to compute energy density and force of peridynamic material. More...
class	PdState
	A class providing methods to compute energy density and force of peridynamic material. More...
class	PmbMaterial
	A class providing methods to compute energy density and force of peridynamic material. More...
class	RnpMaterial
	A class providing methods to compute energy density and force of peridynamic material. More...

Functions
void	computeStateMx (data::ModelData *model, bool compute_in_parallel=false)
	Computes the moment \( m_x \) term in state-based peridynamic formulation.
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)
Conversion methods
double	toNu (double lambda, double mu)
	Compute Poisson's ratio from Lame parameters.
double	toNuEG (double E, double G)
	Compute Poisson's ratio from Young's and Shear modulus.
double	toE (double K, double nu)
	Compute Young's modulus E from Bulk modulus K and Poisson's ratio nu.
double	toK (double E, double nu)
	Compute Bulk modulus K from Young's modulus K and Poisson's ratio nu.
double	toLambdaE (double E, double nu)
	Compute Lame first parameter lambda from Young's modulus E and Poisson's ratio nu.
double	toLambdaK (double K, double nu)
	Compute Lame first parameter lambda from Bulk modulus K and Poisson's ratio nu.
double	toGE (double E, double nu)
	Compute shear modulus from Young's modulus E and Poisson's ratio nu.
double	toGK (double K, double nu)
	Compute shear modulus from Bulk modulus K and Poisson's ratio nu.
double	toELambda (double lambda, double nu)
	Compute Young's modulus E from Lame first parameter lambda and Poisson's ratio nu.
double	toGc (double KIc, double nu, double E)
	Compute critical energy release rate Gc from critical stress-intensity factor KIc, Poisson's ratio nu, and Young's modulus E.
double	toKIc (double Gc, double nu, double E)
	Compute critical stress-intensity factor KIc from critical energy release rate Gc, Poisson's ratio \( nu\), and Young's modulus E.

Function Documentation

computeHydrostaticStrain()

void material::computeHydrostaticStrain	(	data::ModelData *	model,
		bool	compute_in_parallel = false
	)

Definition at line 296 of file materialUtil.cpp.

                                                                                      {
 
  model->d_thetaX.resize(model->d_x.size());
  if (!compute_in_parallel) {
    for (size_t i = 0; i < model->d_x.size(); i++) {
      const auto &pti = model->getPtId(i);
      const auto &particle = model->getParticleFromAllList(pti);
 
      auto thetax = computeHydrostaticStrainI(i, model->d_xRef, model->d_u,
                                        model->d_vol,
                                        model->d_neighPd, model->d_neighPdSqdDist,
                                        particle->getMeshSize(),
                                        particle->getMaterial(),
                                        model->d_fracture_p.get(),
                                        particle->getDimension());
 
      model->setThetax(i, thetax);
    }
  } else {
 
    tf::Executor executor(util::parallel::getNThreads());
    tf::Taskflow taskflow;
 
    taskflow.for_each_index(
      (std::size_t) 0, model->d_x.size(), (std::size_t) 1, [model](std::size_t i) {
        const auto &pti = model->getPtId(i);
        const auto &particle = model->getParticleFromAllList(pti);
 
        auto thetax = computeHydrostaticStrainI(i, 
          model->d_xRef, 
          model->d_u,
          model->d_vol,
          model->d_neighPd, model->d_neighPdSqdDist,
          particle->getMeshSize(),
          particle->getMaterial(),
          model->d_fracture_p.get(),
          particle->getDimension());
 
        model->setThetax(i, thetax);
      }
    ); // for_each
 
    executor.run(taskflow).get();
  }
 
}

References anonymous_namespace{materialUtil.cpp}::computeHydrostaticStrainI(), data::ModelData::d_fracture_p, data::ModelData::d_neighPd, data::ModelData::d_neighPdSqdDist, data::ModelData::d_thetaX, data::ModelData::d_u, data::ModelData::d_vol, data::ModelData::d_x, data::ModelData::d_xRef, util::parallel::getNThreads(), data::ModelData::getParticleFromAllList(), data::ModelData::getPtId(), and data::ModelData::setThetax().

Here is the call graph for this function:

computeStateMx()

void material::computeStateMx	(	data::ModelData *	model,
		bool	compute_in_parallel = false
	)

Computes the moment \( m_x \) term in state-based peridynamic formulation.

Definition at line 214 of file materialUtil.cpp.

                                                                            {
 
  model->d_mX.resize(model->d_x.size());
  if (!compute_in_parallel) {
    for (size_t i = 0; i < model->d_x.size(); i++) {
      const auto &pti = model->getPtId(i);
      const auto &particle = model->getParticleFromAllList(pti);
      auto mx = computeStateMxI(i, model->d_xRef, model->d_vol,
                                model->d_neighPd, model->d_neighPdSqdDist,
                                particle->getMeshSize(),
                                particle->getMaterial());
 
      model->setMx(i, mx);
    }
  } else {
 
    tf::Executor executor(util::parallel::getNThreads());
    tf::Taskflow taskflow;
 
    taskflow.for_each_index(
      (std::size_t) 0, model->d_x.size(), (std::size_t) 1, [model](std::size_t i) {
        const auto &pti = model->getPtId(i);
        const auto &particle = model->getParticleFromAllList(pti);
        auto mx = computeStateMxI(i, model->d_xRef, model->d_vol,
                                  model->d_neighPd, model->d_neighPdSqdDist,
                                  particle->getMeshSize(),
                                  particle->getMaterial());
 
        model->setMx(i, mx);
      }
    ); // for_each
 
    executor.run(taskflow).get();
  }
}

References anonymous_namespace{materialUtil.cpp}::computeStateMxI(), data::ModelData::d_mX, data::ModelData::d_neighPd, data::ModelData::d_neighPdSqdDist, data::ModelData::d_vol, data::ModelData::d_x, data::ModelData::d_xRef, util::parallel::getNThreads(), data::ModelData::getParticleFromAllList(), data::ModelData::getPtId(), and data::ModelData::setMx().

Referenced by PeriDEMModel::init().

Here is the call graph for this function:

Here is the caller graph for this function:

computeStateThetax()

void material::computeStateThetax	(	data::ModelData *	model,
		bool	compute_in_parallel = false
	)

Definition at line 250 of file materialUtil.cpp.

                                                                                {
 
  model->d_thetaX.resize(model->d_x.size());
  if (!compute_in_parallel) {
    for (size_t i = 0; i < model->d_x.size(); i++) {
      const auto &pti = model->getPtId(i);
      const auto &particle = model->getParticleFromAllList(pti);
 
      auto thetax = computeStateThetaxI(i, model->d_xRef, model->d_u,
                                        model->d_vol,
                                        model->d_neighPd, model->d_neighPdSqdDist,
                                        particle->getMeshSize(),
                                      particle->getMaterial(),
                                        model->d_fracture_p.get(),
                                        model->d_mX);
 
      model->setThetax(i, thetax);
 
    }
 
  } else {
 
    tf::Executor executor(util::parallel::getNThreads());
    tf::Taskflow taskflow;
 
    taskflow.for_each_index(
      (std::size_t) 0, model->d_x.size(), (std::size_t) 1, [model](std::size_t i) {
        const auto &pti = model->getPtId(i);
        const auto &particle = model->getParticleFromAllList(pti);
 
        auto thetax = computeStateThetaxI(i, model->d_xRef, model->d_u,
                                          model->d_vol,
                                          model->d_neighPd, model->d_neighPdSqdDist,
                                          particle->getMeshSize(),
                                          particle->getMaterial(),
                                          model->d_fracture_p.get(),
                                          model->d_mX);
 
        model->setThetax(i, thetax);
      }
    ); // for_each
 
    executor.run(taskflow).get();
  }
}

References anonymous_namespace{materialUtil.cpp}::computeStateThetaxI(), data::ModelData::d_fracture_p, data::ModelData::d_mX, data::ModelData::d_neighPd, data::ModelData::d_neighPdSqdDist, data::ModelData::d_thetaX, data::ModelData::d_u, data::ModelData::d_vol, data::ModelData::d_x, data::ModelData::d_xRef, util::parallel::getNThreads(), data::ModelData::getParticleFromAllList(), data::ModelData::getPtId(), and data::ModelData::setThetax().

Here is the call graph for this function:

toE()

double material::toE ( double  K, double  nu  )

inline

Compute Young's modulus E from Bulk modulus K and Poisson's ratio nu.

Parameters

K	Bulk modulus
nu	Poisson's ratio

Returns

E Young's modulus

Definition at line 51 of file materialUtil.h.

{ return K * (3. * (1. - 2. * nu)); }

Referenced by anonymous_namespace{main.cpp}::buildInputJson(), anonymous_namespace{main.cpp}::buildInputJson(), anonymous_namespace{main.cpp}::buildInputJson(), anonymous_namespace{main.cpp}::buildInputJson(), anonymous_namespace{main.cpp}::buildInputJson(), anonymous_namespace{main.cpp}::buildInputJson(), and getInputJson().

Here is the caller graph for this function:

toELambda()

double material::toELambda ( double  lambda, double  nu  )

inline

Compute Young's modulus E from Lame first parameter lambda and Poisson's ratio nu.

Parameters

lambda	Lame first parameter
nu	Poisson's ratio

Returns

E Young's modulus

Definition at line 106 of file materialUtil.h.

                                                    {
    return lambda * (1. + nu) * (1. - 2. * nu) / nu;
  }

toGc()

double material::toGc ( double  KIc, double  nu, double  E  )

inline

Compute critical energy release rate Gc from critical stress-intensity factor KIc, Poisson's ratio nu, and Young's modulus E.

Below conversion from KIc to Gc assumes plane-stress condition. For plane-stress condition, we need to modify the Young's modulus \( E\) to \( \frac{E}{1 - \nu^2} \) where \( \nu\) is the Poisson's ratio.

Parameters

KIc	Critical stress-intensity factor
nu	Poisson's ratio
E	Young's modulus

Returns

Gc Critical energy release rate

Definition at line 123 of file materialUtil.h.

{ return KIc * KIc / E; }

Referenced by anonymous_namespace{main.cpp}::buildInputJson().

Here is the caller graph for this function:

toGE()

double material::toGE ( double  E, double  nu  )

inline

Compute shear modulus from Young's modulus E and Poisson's ratio nu.

Parameters

E	Young's modulus
nu	Poisson's ratio

Returns

G Shear modulus

Definition at line 87 of file materialUtil.h.

{ return E / (2. * (1. + nu)); }

Referenced by anonymous_namespace{main.cpp}::buildInputJson(), anonymous_namespace{main.cpp}::buildInputJson(), anonymous_namespace{main.cpp}::buildInputJson(), anonymous_namespace{main.cpp}::buildInputJson(), anonymous_namespace{main.cpp}::buildInputJson(), anonymous_namespace{main.cpp}::buildInputJson(), and getInputJson().

Here is the caller graph for this function:

toGK()

double material::toGK ( double  K, double  nu  )

inline

Compute shear modulus from Bulk modulus K and Poisson's ratio nu.

Parameters

K	Bulk modulus
nu	Poisson's ratio

Returns

G Shear modulus

Definition at line 95 of file materialUtil.h.

                                          {
    return 3. * K * (1. - 2. * nu) / (2. * (1. + nu));
  }

toK()

double material::toK ( double  E, double  nu  )

inline

Compute Bulk modulus K from Young's modulus K and Poisson's ratio nu.

Parameters

E	Young's modulus
nu	Poisson's ratio

Returns

K Bulk modulus

Definition at line 59 of file materialUtil.h.

{ return E / (3. * (1. - 2. * nu)); }

toKIc()

double material::toKIc ( double  Gc, double  nu, double  E  )

inline

Compute critical stress-intensity factor KIc from critical energy release rate Gc, Poisson's ratio \( nu\), and Young's modulus E.

Below conversion from Gc to KIc assumes plane-stress condition. For plane-stress condition, we need to modify the Young's modulus \( E\) to \( \frac{E}{1 - \nu^2} \) where \( \nu\) is the Poisson's ratio.

Parameters

Gc	Critical energy release rate
nu	Poisson's ratio
E	Young's modulus

Returns

KIc Critical stress-intensity factor

Definition at line 138 of file materialUtil.h.

{ return std::sqrt(Gc * E); }

toLambdaE()

double material::toLambdaE ( double  E, double  nu  )

inline

Compute Lame first parameter lambda from Young's modulus E and Poisson's ratio nu.

Parameters

E	Young's modulus
nu	Poisson's ratio

Returns

lambda Lame first parameter

Definition at line 68 of file materialUtil.h.

                                               {
    return E * nu / ((1. + nu) * (1. - 2. * nu));
  }

toLambdaK()

double material::toLambdaK ( double  K, double  nu  )

inline

Compute Lame first parameter lambda from Bulk modulus K and Poisson's ratio nu.

Parameters

K	Bulk modulus
nu	Poisson's ratio

Returns

lambda Lame first parameter

Definition at line 79 of file materialUtil.h.

{ return 3. * K * nu / (1. + nu); }

toNu()

double material::toNu ( double  lambda, double  mu  )

inline

Compute Poisson's ratio from Lame parameters.

Parameters

lambda	Lame first parameter
mu	Lame second parameter

Returns

nu Poisson's ratio

Definition at line 35 of file materialUtil.h.

{ return lambda * 0.5 / (lambda + mu); }

toNuEG()

double material::toNuEG ( double  E, double  G  )

inline

Compute Poisson's ratio from Young's and Shear modulus.

Parameters

E	Youngs modulus
G	Shear modulus

Returns

nu Poisson's ratio

Definition at line 43 of file materialUtil.h.

{ return E * 0.5 / G - 1.; }

updateBondFractureData()

void material::updateBondFractureData	(	data::ModelData *	model,
		bool	compute_in_parallel = false
	)

Definition at line 343 of file materialUtil.cpp.

                                                                                    {
 
  if (!compute_in_parallel) {
    for (size_t i = 0; i < model->d_x.size(); i++) {
      const auto &pti = model->getPtId(i);
      const auto &particle = model->getParticleFromAllList(pti);
 
      updateBondFractureDataI(i, model->d_xRef, model->d_neighPd, model->d_u,
                              particle->getMaterial(),
                                              model->d_fracture_p.get());
    }
  } else {
 
    tf::Executor executor(util::parallel::getNThreads());
    tf::Taskflow taskflow;
 
    taskflow.for_each_index(
      (std::size_t) 0, model->d_x.size(), (std::size_t) 1, [model](std::size_t i) {
        const auto &pti = model->getPtId(i);
        const auto &particle = model->getParticleFromAllList(pti);
 
        updateBondFractureDataI(i, 
          model->d_xRef, 
          model->d_neighPd, 
          model->d_u,
          particle->getMaterial(),
          model->d_fracture_p.get());
      }
    ); // for_each
 
    executor.run(taskflow).get();
  }
}

References data::ModelData::d_fracture_p, data::ModelData::d_neighPd, data::ModelData::d_u, data::ModelData::d_x, data::ModelData::d_xRef, util::parallel::getNThreads(), data::ModelData::getParticleFromAllList(), data::ModelData::getPtId(), and anonymous_namespace{materialUtil.cpp}::updateBondFractureDataI().

Here is the call graph for this function: