https://mooseframework.inl.gov
Functions
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  ,
zeroIncrement   
)

Definition at line 83 of file LagrangianObjectiveRateTest.C.

84 {
85  Inputs in = makeInputs(isotropicC(), RankTwoTensor()); // dL = 0 -> dS = 0, dW = 0
90  const RankTwoTensor expected = in.cauchy_stress_old + in.dS;
91  EXPECT_LT((truesdell(in, false).cauchy_stress - expected).L2norm(), 1.0e-12);
92  EXPECT_LT((jaumann(in, false).cauchy_stress - expected).L2norm(), 1.0e-12);
93  EXPECT_LT((greenNaghdi(in, false).cauchy_stress - expected).L2norm(), 1.0e-12);
94  EXPECT_LT((rashid(in, false).cauchy_stress - expected).L2norm(), 1.0e-12);
95 }
T L2norm(const RankTwoTensorTempl< T > &r2tensor)
Per-qp quantities a rate reads.
Outputs greenNaghdi(const Inputs &in, bool need_jacobian)
sigma_{n+1} = J(dO)^{-1} (sigma_n + Deltasigma), with dO built from the polar rotation increment...
Outputs jaumann(const Inputs &in, bool need_jacobian)
sigma_{n+1} = J(dW)^{-1} (sigma_n + Deltasigma). Jaumann rate.
static RankTwoTensorTempl Identity()
RankTwoTensor dS
Constitutive small-stress increment (_small_stress - _small_stress_old).
RankTwoTensor cauchy_stress_old
Cumulative Cauchy stress at step n.
RankTwoTensor rotation
Green-Naghdi only: polar-decomposition rotation R (n+1 and n), the inverse incremental deformation gr...
Outputs truesdell(const Inputs &in, bool need_jacobian)
sigma_{n+1} = J(dL)^{-1} (sigma_n + Deltasigma). Truesdell rate.
Outputs rashid(const Inputs &in, bool need_jacobian)
sigma_{n+1} = r_hat (sigma_n + Deltasigma) r_hat^T, r_hat = exp(Deltaw).

◆ TEST() [2/4]

TEST ( LagrangianObjectiveRateTest  ,
truesdellJacobian   
)

Definition at line 98 of file LagrangianObjectiveRateTest.C.

99 {
100  checkJacobianFD(&truesdell, isotropicC(), RankTwoTensor(0.03, -0.02, 0.01, 0.015, -0.01, 0.02));
101 }
Outputs truesdell(const Inputs &in, bool need_jacobian)
sigma_{n+1} = J(dL)^{-1} (sigma_n + Deltasigma). Truesdell rate.

◆ TEST() [3/4]

TEST ( LagrangianObjectiveRateTest  ,
jaumannJacobian   
)

Definition at line 103 of file LagrangianObjectiveRateTest.C.

104 {
105  checkJacobianFD(&jaumann, isotropicC(), RankTwoTensor(0.03, -0.02, 0.01, 0.015, -0.01, 0.02));
106 }
Outputs jaumann(const Inputs &in, bool need_jacobian)
sigma_{n+1} = J(dW)^{-1} (sigma_n + Deltasigma). Jaumann rate.

◆ TEST() [4/4]

TEST ( LagrangianObjectiveRateTest  ,
rashidJacobian   
)

Definition at line 108 of file LagrangianObjectiveRateTest.C.

109 {
110  // Add an antisymmetric part so r_hat = exp(skew(dL)) != I and the Rodrigues chain is exercised.
111  RankTwoTensor dL(0.03, -0.02, 0.01, 0.015, -0.01, 0.02);
112  dL(0, 1) += 0.02;
113  dL(1, 0) -= 0.02;
114  checkJacobianFD(&rashid, isotropicC(), dL);
115 }
Outputs rashid(const Inputs &in, bool need_jacobian)
sigma_{n+1} = r_hat (sigma_n + Deltasigma) r_hat^T, r_hat = exp(Deltaw).