Model::Solid::Finite::CrystalPlastic

src/Model/Solid/Finite/CrystalPlastic.H

A basic and relatively untested implementation of finite-deformation FCC crystal plasticity. Use with caution.

Inherits from Model::Solid::Finite::PseudoLinear::Cubic to provide elastic response.

Plastic flow modeled as:

\[\mathbf{F}^p\mathbf{F}^{p-1} = \sum \dot{\gamma}_n \mathbf{a}_n\otimes\mathbf{N}_n\]

where \(\mathbf{F}=\mathbf{F}^e\mathbf{F}^p\) and \(\gamma_n\) are slips on system \(n\).

Power law rate hardening is used in the integration of slips:

\[\dot{\gamma}_n = \dot{\gamma}_0 \Big(\frac{\tau}{\tau_{crss}}\Big)^m \operatorname{sign}(\tau)\]

\(\phantom{0}\)

tau_crss
conditional

Critical resolved shear stress \(\tau_{crss}\)

Names by executable 2 executables, 2 uses
mechanics-2d-coverage-g++
model0.tau_crss
mechanics-3d-coverage-g++
model0.tau_crss
gammadot0
1 conditional

Rate hardening coefficient \(\dot{\gamma}_0\)

Names by executable 2 executables, 2 uses
mechanics-2d-coverage-g++
model0.gammadot0
mechanics-3d-coverage-g++
model0.gammadot0
m_rate_inv
0.5 conditional

Inverse of the hardening exponent \(\frac{1}{m}\)

Names by executable 2 executables, 2 uses
mechanics-2d-coverage-g++
model0.m_rate_inv
mechanics-3d-coverage-g++
model0.m_rate_inv
tstart
0 conditional

Time to activate plastic slip

Names by executable 2 executables, 2 uses
mechanics-2d-coverage-g++
model0.tstart
mechanics-3d-coverage-g++
model0.tstart