静态自洽计算