  We explore reshaping of nematoelastic films upon imbibing an isotropic
solvent under conditions when isotropic and nematic phases coexist. The
structure of the interphase boundary is computed taken into account the optimal
nematic orientation governed by interaction of gradients of the nematic order
parameter and solvent concentration. This structure determines the effective
line tension of the boundary. We further compute equilibrium shapes of deformed
thin sheets and cylindrical and spherical shells with the rectilinear or
circular shape of the boundary between nematic and isotropic domains. A
differential expansion or contraction near this boundary generates a folding
pattern spreading out into the bulk of both phases. The hierarchical ordering
of this pattern is most pronounced on a cylindrical shell.
