LagrangianObjectiveRateTest.C File Reference

Go to the source code of this file.

Functions
	TEST (LagrangianObjectiveRateTest, zeroIncrement)
	TEST (LagrangianObjectiveRateTest, truesdellJacobian)
	TEST (LagrangianObjectiveRateTest, jaumannJacobian)
	TEST (LagrangianObjectiveRateTest, rashidJacobian)

Function Documentation

TEST() [1/4]

TEST	(	LagrangianObjectiveRateTest	,
		jaumannJacobian
	)

Definition at line 103 of file LagrangianObjectiveRateTest.C.

{
  checkJacobianFD(&jaumann, isotropicC(), RankTwoTensor(0.03, -0.02, 0.01, 0.015, -0.01, 0.02));
}

TEST() [2/4]

TEST	(	LagrangianObjectiveRateTest	,
		rashidJacobian
	)

Definition at line 108 of file LagrangianObjectiveRateTest.C.

{
  // Add an antisymmetric part so r_hat = exp(skew(dL)) != I and the Rodrigues chain is exercised.
  RankTwoTensor dL(0.03, -0.02, 0.01, 0.015, -0.01, 0.02);
  dL(0, 1) += 0.02;
  dL(1, 0) -= 0.02;
  checkJacobianFD(&rashid, isotropicC(), dL);
}

TEST() [3/4]

TEST	(	LagrangianObjectiveRateTest	,
		truesdellJacobian
	)

Definition at line 98 of file LagrangianObjectiveRateTest.C.

{
  checkJacobianFD(&truesdell, isotropicC(), RankTwoTensor(0.03, -0.02, 0.01, 0.015, -0.01, 0.02));
}

TEST() [4/4]

TEST	(	LagrangianObjectiveRateTest	,
		zeroIncrement
	)

Definition at line 83 of file LagrangianObjectiveRateTest.C.

{
  Inputs in = makeInputs(isotropicC(), RankTwoTensor()); // dL = 0 -> dS = 0, dW = 0
  in.rotation = RankTwoTensor::Identity();
  in.rotation_old = RankTwoTensor::Identity();
  in.inv_df = RankTwoTensor::Identity();
  in.inv_def_grad = RankTwoTensor::Identity();
  const RankTwoTensor expected = in.cauchy_stress_old + in.dS;
  EXPECT_LT((truesdell(in, false).cauchy_stress - expected).L2norm(), 1.0e-12);
  EXPECT_LT((jaumann(in, false).cauchy_stress - expected).L2norm(), 1.0e-12);
  EXPECT_LT((greenNaghdi(in, false).cauchy_stress - expected).L2norm(), 1.0e-12);
  EXPECT_LT((rashid(in, false).cauchy_stress - expected).L2norm(), 1.0e-12);
}