Modeling, analysis and computation are three pillars of computational science. We discuss them within the context of liquid crystal networks (LCNs). These materials couple a nematic liquid crystal with a rubbery material. When actuated with heat or light, the interaction of the liquid crystal with the rubber creates complex shapes. Thin bodies of LCNs are natural candidates for soft robotics applications. We start from the classical 3D trace energy formula and derive a reduced 2D membrane energy as the formal asymptotic limit of vanishing thickness, including both stretching and bending energies, and characterize the zero energy deformations. We design a sound numerical method and discuss its Gamma convergence. We present computations showing the geometric effects that arise from liquid crystal defects as well as computations of nonisometric origami within and beyond theory. This work is joint with L. Bouck, G. Benavides, and S. Yang.