黄永刚单晶体塑性的UMAT原始文件，不是网上流传的版本(huang-crystal plasticity)

平面波展开法计算声子晶体，单包为正方形，散射体为圆形状(Calculation of phononic crystal by plane wave expansion method. A single packet is a square, and the scatterer is a circle shape.)