small_elastic_tri_states.C File Reference

Go to the source code of this file.

Functions
void	umat_ (double *stress, double *statev, double *ddsdde, double *sse, double *spd, double *scd, double *rpl, double *ddsddt, double *drplde, double *drpldt, double *stran, double *dstran, double *time, double *dtime, double *temp, double *dtemp, double *predef, double *dpred, char *cmname, int *ndi, int *nshr, int *ntens, int *nstatv, double *props, int *nprops, double *coords, double *drot, double *pnewdt, double *celent, double *dfgrd0, double *dfgrd1, int *noel, int *npt, int *layer, int *kspt, int *kstep, int *kinc, short cmname_len)

Function Documentation

umat_()

void umat_	(	double *	stress,
		double *	statev,
		double *	ddsdde,
		double *	sse,
		double *	spd,
		double *	scd,
		double *	rpl,
		double *	ddsddt,
		double *	drplde,
		double *	drpldt,
		double *	stran,
		double *	dstran,
		double *	time,
		double *	dtime,
		double *	temp,
		double *	dtemp,
		double *	predef,
		double *	dpred,
		char *	cmname,
		int *	ndi,
		int *	nshr,
		int *	ntens,
		int *	nstatv,
		double *	props,
		int *	nprops,
		double *	coords,
		double *	drot,
		double *	pnewdt,
		double *	celent,
		double *	dfgrd0,
		double *	dfgrd1,
		int *	noel,
		int *	npt,
		int *	layer,
		int *	kspt,
		int *	kstep,
		int *	kinc,
		short	cmname_len
	)

Definition at line 17 of file small_elastic_tri_states.C.

{
  auto emod = props[0];

  const double temperature = temp[0];
  const double voltage = predef[0];

  // External field-dependent contribution
  emod *= (temperature + voltage) / 1000;

  // State-dependent contribution of elasticity modulus.
  emod *= (1 + statev[0] * 10.0);

  auto enu = props[1];
  auto ebulk3 = emod / (1.0 - 2.0 * enu);
  auto eg2 = emod / (1.0 + enu);
  auto eg = eg2 / 2.0;
  auto elam = (ebulk3 - eg2) / 3.0;

  // elastic stiffness
  for (int k1 = 0; k1 < *ndi; ++k1)
  {
    for (int k2 = 0; k2 < *ndi; ++k2)
      ddsdde[k1 * *ntens + k2] = elam;
    ddsdde[k1 * *ntens + k1] += eg2;
  }
  for (int k1 = *ndi; k1 < *ntens; ++k1)
    ddsdde[k1 * *ntens + k1] = eg;

  // calculate stress
  for (int k1 = 0; k1 < *ntens; ++k1)
    for (int k2 = 0; k2 < *ntens; ++k2)
      stress[k1] += ddsdde[k1 * *ntens + k2] * dstran[k2];

  statev[0] += std::fabs(dstran[0]) + std::fabs(dstran[1]);
}