有限元中的裂纹尖端奇异元问题

在计算裂纹体的应力强度因子时，通常使用奇异元素划分裂纹尖端附近的区域，以满足裂纹尖端的奇异性要求。

用有限元方法解决裂纹尖端问题时，如何通过程序描述八节点四边形奇异元？