To make things work with FCH we need a way to specify the active strain. I don't yet fully understand the mechanics but it amounts to, ultimately, computing an additional tensor FF_A s.t. we use PK1(FF) * FF_A^-1 to compute the load vector instead of just PK1(FF).
These should be supplied to Parts and then propagated down to compute_load_vector().
The original problem description is wrong - we need to modify FF and dependent quantities. Fortunately MechanicsValues can be easily modified to this end.
To make things work with FCH we need a way to specify the active strain. I don't yet fully understand the mechanics but it amounts to, ultimately, computing an additional tensor
FF_A
s.t. we usePK1(FF) * FF_A^-1
to compute the load vector instead of justPK1(FF)
.These should be supplied to
Part
s and then propagated down tocompute_load_vector()
.