非线性 | 弧长法(Arc-Length Methods)_非线性

非线性 | 弧长法(Arc-Length Methods)

图1所示为弧长法求解过程，若以下标表示第个荷载步，上标表示第个荷载步下的第次迭代，显然，当荷载增量 ,则迭代路径为一条水平直线，即为著名的牛顿-拉夫逊方法。对于图2所示的求解问题，牛顿-拉夫逊方法不能跨过极值点得到完整的荷载-位移曲线。因此，弧长法最重要的就是求荷载增量。

而弧长法的荷载增量是变化的，可自动控制荷载，这样在原方程组的基础之上又增加了一个未知数，因此需要额外补充一个方程。如图3所示，某一荷载步迭代至收敛时总有

考虑系统方程组

在迭代过程中，逐渐趋于0，如果这两个值都为0，则说明该荷载步的迭代已收敛。在上一个迭代收敛点(如图1中的 )将作一阶泰勒展开

即

令，则

弧长法通过自动建立适当的荷载增量进一步优化了牛顿-拉夫逊方法，使用弧长法，可以跟踪复杂的荷载-变形路径。

来源：数值分析与有限元编程

免费Fortran编译器oneAPI

Parallel Studio XE 是商业版的高性能计算工具。学生可以申请一年试用期，到期要重新申请许可，稍显麻烦。现在来试试Intel推出的免费Fortran编译器oneAPI。不要许可，开箱即用。oneAPI包含Base Toolkit 和HPC Toolkit 两部分，Base Toolkit 里面有MKL函数库，如果不用MKL，这部分可以不用安装。Fortran编译器则在HPC Toolkit 里面，这是必须安装的。配置方法：(一) 安装VS2017或者VS2019记得安装时要勾选“C++桌面开发”这个选项。(二) 安装Intel oneAPI Base Toolkit如果不用MKL，可以跳过这一步。https://software.intel.com/content/www/us/en/develop/tools/oneapi/all-toolkits.html打开这个网址，往下拉，找到如图所示的位置。可以看到Intel oneAPI Base Toolkit包含许多组件，由于我们只需要MKL，故选择在线安装更快。点击Download，会跳转到另一个页面选择在线安装方式，点击Download，得到一个安装文件w_BaseKit_p_2021.1.0.2664.exe，这和用Visual Studio installer在线安装VS2019一样，边下载边安装。注意安装时选择costomize，即自定义模式，只需将MKL那个选上即可。安装完成后将安装路径下的两个文件夹ia32和intel64添加到环境变量path中。(三) 安装Intel oneAPI HPC Toolkit同样是上面那个网址，Intel oneAPI Base Toolkit的右边就是HPC Toolkit了。这里跟前面一样，点击Download，跳转到另一个页面选择在线安装方式，点击Download，得到一个安装文件w_HPCKit_p_2021.1.0.2682.exe，可以一边下载一边安装。注意安装时选择costomize，即自定义模式，只需将Fortran编译器和MPI这两个选项选上即可。安装完成打开VS2019,新建项目里面就有Fortran了。来源：数值分析与有限元编程