Faults有限元程序在Windows平台上的移植

Transplant the Faults finite element code to Windows platform

Faults有限元软件综合高程、热流、断层等多个地球物理场参数,相较于Ansys、Adina等商用有限元软件能更好的处理断层的摩擦机制,因此Faults有限元软件在地学模拟中有独特优势。目前使用的Faults软件基于IBM科学与工程数学库(ESSL),使用环境较为苛刻,使用不便,需要移植源代码。利用Linpack免费数学库替代IBM服务器和科学工程数学库(ESSL)。文中阐述了利用Linpack替代原程序解带状矩阵数学库的方法,解决了刚度矩阵的存储和求解问题,对比理想边界条件下的计算结果,与原程序速度场计算结果进行对比。结果显示,对Faults有限单元程序在Windows平台上的移植是成功和可靠的。

The finite element software ＂Faults＂ considers a lot of geophysics parameters hke elevation, heat-flow, etc. To compare with other FEM software like Ansys, Adina, ＂Faults＂ is more convenient to handle the friction mechanism of fault. So ＂Faults＂ has advantages in FEM simulation of geophysics problem. But right now, ＂Faults＂ is based on IBM server and Engineering and Scientific Library which are very expensive. In this paper we demonstrate how to use free math libraries ＂linpack＂ to substitute the IBM server and Engineering and Scientific Library which is used by the original code. We have introduced the method to transplant and compared the computed result with the that published by original author. The result of comparing shows that the transplant is successful and reliable.