“…Here, we adopt the neoclassical strain-energy function, and assume the isotropic phase at high temperature as the reference configuration [8,[11][12][13][14], instead of the cross-linking nematic phase [2,5,56,[62][63][64]74]. We exploit theoretically the multiplicative decomposition of the deformation gradient from the reference configuration to the current configuration into an elastic distortion followed by a natural shape change [20,[41][42][43][44]. This multiplicative decomposition is similar to those found in the constitutive theories of thermoelasticity, elastoplasticity, and morphoelasticity [19,37] (see [18,53] as well), but it is also different in the sense that the stress-free geometric change is superposed on the elastic deformation, which is directly applied to the reference state.…”