子程序UANISOHYPER_INV有效性验证及问题

大勇

我写了一个模拟各向异性超弹性材料的UANISOHYPER_INV的子程序，有些问题，这里贴出来希望大家给参详参详将应变能方程写为如下形式方便编写子程序，这里我的材料方向为1个，为z方向。

子程序UANISOHYPER_INV需要给出W对各个不变量的1到3阶偏导数和应变能的表达式，利用abaqus自带的核心功能可以对这种材料进行模拟。
验证子程序的时候几何模型为一块板，分别进行理论计算和abaqus建模计算，比较两者的应变能值
1. 给定一个位移场，x和y方向的伸长从1到3，对应的应变能比较如下


2. 施加力载荷的时候就出问题了，我在x方向施加了一个拉伸的pressure载荷，载荷0.1Mpa，出现初始步不收敛的问题。
问题很奇怪，施加位移场，计算没问题，但是施加力载荷出现不收敛的情况。
有可能是什么地方的问题呢？

大勇

                    

as997788

我想问一下，就是自己写的程序需要把等于0的量写出来吗，因为我的能量密度仅和I4相关，和其他的无关，需要把W对其他不变量的等于0的导数写出来吗