用comsol实现计算E^H

算例放在Homogenization.mph里面，主要分两块：
1.方程书写（Buv=Fv，Buv和Fv的确定见《基于均匀化方法的材料参数化优化设计研究》P22）；
2.周期边界条件（主要是采用线性拉伸算子，linear extrusion，先在变量里面定义，再加狄利克雷条件）。

几个测试：
1.不加其他边界条件无法收敛，加1个固定点约束，加4个端点为固定点约束都无法收敛；
2.加狄利克雷边界条件时，加1条边，2条边都能收敛，当不挖孔时，怎么加条件对结果没影响，此时的位移1e-15级别，本来就E^H=E，挖孔之后，加不同的狄利克雷条件会得到不同的E^H

更新：
由于未考虑泊松比的变化，已经重做，把算例放在comsol_c&s文件夹中