The thermo-order-mechanical behaviors of liquid crystal elastomers （LCEs） under biaxial loading are studied in this paper. Inverse method for nonlinear elastic problems is utilized by imposing biaxial stretching to thin rectangular samples. Neo-classical elastic energy is used together with the Landau-de Gennes nematic free energy. Under plane stress assumptions, the constitutive equations are derived. Due to the possible reorientations of the liquid crystal molecules induced by the imposed biaxial loading, the in-plane nonlinear stress-strain relations can have different expressions depending on which loading axis will have the largest effective principal strain. And the free energy is a multi-well non-convex potential function. As shown by some typical loading paths, the LCE samples will exhibit an anisotropic nonlinear elastic behavior, as long as the loading has not induced a reorientation of the liquid crystal molecules. When this did occur, jumps of stresses could take place for dead loadings due to the losing of stability.
