利用k-means算法实现二维平面点的聚类，包括了运行源代码和结果图(The k-means algorithm is used to realize the clustering of two-dimensional plane points, including the running source code and the result graph)

纳入了机械损伤和应变速率的断裂UMAT子程序，不能用于平面应力。(The fracture UMAT subroutine incorporating mechanical damage and strain rate can not be used for plane stress.)

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