摘要
为开发一个相比于商业有限元软件更加透明、操作性更高的有限元分析系统,本文编写了一个针对镦粗变形过程模拟的刚塑性有限元程序FEM_FORGE.此程序使用Fortran语言,以刚塑性有限元为理论基础,采用体积可压缩法,选用了比较简单的四边形等参单元,采用线性方法确定初始速度场,应用具有二次收敛性的Newton-Raphson法求解非线性方程组,然后通过不求导的Brent法优化计算结果,最后在MATLAB环境中进行后处理.目前,此程序仅能够模拟等温的镦粗过程,可方便地改变网格密度、变形条件,也可方便地显示应力、应变、应变速率以及任何用户自定义的场量.