液晶弹性自由能计算