GET_DF_PLANESTRESSHARDMIXVAF_FD - evaluates the yield function and its derivatives using finite differences

Comments

Assumptions: plane stress, isotropic hardening (Voce), kinematic hardening (Armstrong-Frederick).

Input Arguments

theta (double) - material parameters

H_isotropic (double) - isotropic hardening parameters

sigmaV (double) - Cauchy stress under plane stress conditions in Voigt notation (sigma_11, sigma_22, sigma_12)

gamma (double) - plastic multiplier

sigmaV_back (double) - back stress

FD (double) - finite differences stepsize

secondDerivative (double) - indicates whether second derivative should be computed

Output Arguments

derivatives of the yield function