Fortran UMAT子程序：为何给AA赋值时F_F矩阵会异常变动

问题现象 #

调试嵌入Abaqus UMAT的Fortran代码过程中，fabric_tensor_Mandel子程序执行时出现参数非预期变动，经观测，F_F矩阵在执行以下循环代码段时会发生值异常变更：

异常观测触发代码：

DO i=1,NTENS
    DO j=1,NTENS
        AA(i,j) =(1/(1+AFFA))*F_F(i,j)
    END DO
END DO

异常产生核心原因 #

• INVERT矩阵求逆子程序调用参数顺序错位，引发栈内存越界写入
  原代码中求逆调用为CALL INVERT(FF_MATRIX,FF_inv,IPIV,6,X,FSMAL,IOUT)，对照Abaqus配套通用INVERT子程序的标准形参定义，参数顺序应为「输入矩阵、矩阵阶数、输出逆矩阵、主元索引数组、临时工作数组、收敛精度阈值、错误输出码」。当前调用把硬编码的矩阵阶数6放置在第4个传参位置，导致子程序内部按错误的参数映射关系读写内存：
  1. 把本该作为主元数组的IPIV参数位置，解析为整型的矩阵阶数，后续所有索引计算完全混乱
  2. 对后续所有参数的内存位置进行错位写入，直接覆盖栈上相邻的F_F、AA等形参数组的存储值
     这种越界写入的实际触发点在INVERT调用阶段，只是因为后续代码没有频繁访问F_F元素，直到执行AA赋值循环、逐值读取F_F时，才观测到数值已经被篡改。
• 未声明临时变量的隐式类型分配风险
  原代码未开启IMPLICIT NONE类型检查，IPIV、X、AFFA、IOUT等变量未做显式维度声明，在老Fortran隐式类型规则下会被分配为单精度标量内存。INVERT子程序内部把这些标量内存当成数组做下标写入时，会直接越界写入相邻栈内存块，进一步加剧参数篡改概率。
• 硬编码矩阵维度不匹配风险
  原代码中INVERT调用硬编码矩阵维度为6，若求解场景为平面应力/平面应变状态（Abaqus环境下对应NTENS=4），INVERT内部按6×6矩阵尺寸读写数据时，会直接超出FF_MATRIX、FF_inv数组的实际分配长度，造成内存越界。

对应修复方案 #

• 开启隐式类型检查，补全所有变量维度声明
  在INCLUDE 'ABA_PARAM.INC'之后添加IMPLICIT NONE，编译阶段即可直接定位所有未声明的变量。参照INVERT子程序的参数需求，补全临时数组的维度声明：

  DIMENSION IPIV(NTENS), X(NTENS)
  ! 补全未声明的标量变量类型
  REAL(8) :: AFFA, FSMAL
  INTEGER :: IOUT
  

• 修正INVERT调用的参数顺序，替换硬编码维度
  按照INVERT的形参顺序调整传参，把硬编码的矩阵阶数替换为实际传入的NTENS变量，避免参数错位和维度不匹配：

  FSMAL = 1.0d-10
  ! 修正后的正确调用形式
  CALL INVERT(FF_MATRIX, NTENS, FF_inv, IPIV, X, FSMAL, IOUT)
  

• 增加独立本地数组副本，隔离形参操作风险
  提前将F_F的数值拷贝到本地独立数组，所有AA赋值计算均使用本地副本，避免其他子程序的越界写入干扰形参F_F的数值：

  DIMENSION F_F_LOCAL(NTENS, NTENS)
  F_F_LOCAL = F_F
  DO i=1,NTENS
      DO j=1,NTENS
          AA(i,j) = (1.0d0/(1.0d0 + AFFA)) * F_F_LOCAL(i,j)
      END DO
  END DO
  

• 增加分步校验输出
  在INVERT调用结束、DOTPROD调用结束、AA赋值循环前后，打印F_F矩阵的对角线元素值，确认内存篡改的精确发生位置，避免对异常触发点的误判。

内容的提问来源于stack exchange，提问作者faezeh iranmanesh