桩土 剑桥模型 WARNING PROGRAM IS ASKED TO INVERT A SINGULAR MATRIX

yeoo111

算例是码头桩土结构，考虑桩土接触，porous elasticity材料收敛，但是渗流固结（Modified cam-clay model）不收敛，出现下面的问题：

INCREMENT     1 STARTS. ATTEMPT NUMBER  1, TIME INCREMENT  0.500   
  NUMBER OF EQUATIONS =   327736     NUMBER OF FLOATING PT. OPERATIONS =  2.20E+12
     CHECK POINT   START OF SOLVER   
     CHECK POINT  END OF SOLVER      
       ELAPSED USER TIME (SEC)      =   754.80   
       ELAPSED SYSTEM TIME (SEC)    =   39.900   
       ELAPSED TOTAL CPU TIME (SEC) =   794.70   
       ELAPSED WALLCLOCK TIME (SEC) =       1077
***WARNING: PROGRAM IS ASKED TO INVERT A SINGULAR MATRIX
***WARNING: PROGRAM IS ASKED TO INVERT A SINGULAR MATRIX
***WARNING: PROGRAM IS ASKED TO INVERT A SINGULAR MATRIX

***WARNING: THE PLASTICITY/CREEP/CONNECTOR FRICTION ALGORITHM DID NOT CONVERGE
             AT 817 POINTS

***NOTE: MATERIAL CALCULATIONS FAILED TO CONVERGE OR WERE NOT ATTEMPTED AT ONE
          OR MORE POINTS. CONVERGENCE IS JUDGED UNLIKELY.

  INCREMENT     1 STARTS. ATTEMPT NUMBER  2, TIME INCREMENT  0.125   
  NUMBER OF EQUATIONS =   327736     NUMBER OF FLOATING PT. OPERATIONS =  2.20E+12
     CHECK POINT   START OF SOLVER   
     CHECK POINT  END OF SOLVER      
       ELAPSED USER TIME (SEC)      =   687.80   
       ELAPSED SYSTEM TIME (SEC)    =   39.300   
       ELAPSED TOTAL CPU TIME (SEC) =   727.10   
       ELAPSED WALLCLOCK TIME (SEC) =        848
***WARNING: PROGRAM IS ASKED TO INVERT A SINGULAR MATRIX
***WARNING: PROGRAM IS ASKED TO INVERT A SINGULAR MATRIX
***WARNING: PROGRAM IS ASKED TO INVERT A SINGULAR MATRIX
***WARNING: PROGRAM IS ASKED TO INVERT A SINGULAR MATRIX
***WARNING: PROGRAM IS ASKED TO INVERT A SINGULAR MATRIX


接着就停止了运算，log文件中提示：

Run standard.exe
Abaqus License Manager checked out the following licenses:
Abaqus/Standard checked out 5 tokens.
9/24/2010 2:49:37 PM
Abaqus Error: The executable E:\SIMULIA\Abaqus\6.10-1\exec\standard.exe
aborted with system error code 5.
Please check the .dat, .msg, and .sta files for error messages if the files
exist.  If there are no error messages and you cannot resolve the problem,
please run the command "abaqus job=support information=support" to report and
save your system information.  Use the same command to run Abaqus that you
used when the problem occurred.  Please contact your local Abaqus support
office and send them the input file, the file support.log which you just
created, the executable name, and the error code.
Abaqus/Analysis exited with errors

单元数目较多，不好上传，请大家帮忙分析一下，***WARNING: PROGRAM IS ASKED TO INVERT A SINGULAR MATRIX
的解决办法，我已经搜索过，但是很少有满意的答案。

材料参数为：

*Material, name=NIANTU （注：粘土）
*Clay Plasticity
  0.2,     1.2,500000,     1.,     1.
*Density
2000.,
*Permeability, specific=1000.
1e-08, 1.2
*Porous Elastic
0.02, 0.3,   0.

msg文件的warning信息：

***WARNING: THERE ARE POINTS IN THE MODEL WHICH HAVE INITIAL STRESSES OUTSIDE
             THE INITIALLY DEFINED YIELD SURFACE.  THE VALUE DEFINING THE
             CENTER OF THE YIELD SURFACE WILL BE ADJUSTED TO REMOVE THE
             INCONSISTENCY.
THE FOLLOWING POINTS ARE OUTSIDE THE INITIALLY DEFINED YIELD SURFACE
             ELEMENT       POINT          OLD CENTER               NEW CENTER
                9906         1              5.0000E+05              2.8357E+06
                9906         2              5.0000E+05              2.8357E+06
                9906         3              5.0000E+05              2.8357E+06
                9906         4              5.0000E+05              2.8357E+06
                9906         5              5.0000E+05              2.8357E+06
                9906         6              5.0000E+05              2.8357E+06
                9906         7              5.0000E+05              2.8357E+06
                9906         8              5.0000E+05              2.8357E+06
                  :
                  :
                  :
               73162         5              5.0000E+05              4.5683E+06
               73162         6              5.0000E+05              4.5683E+06
               73162         7              5.0000E+05              4.5683E+06
               73162         8              5.0000E+05              4.5683E+06
               73163         1              5.0000E+05              2.1986E+06
               73163         2              5.0000E+05              2.1986E+06
               73163         3              5.0000E+05              2.1986E+06
               73163         4              5.0000E+05              2.1986E+06
               73163         5              5.0000E+05              2.1986E+06
（注：有些单元并没有用剑桥模型，但这儿也显示了改Pa，不知是为什么？奇怪！）
***WARNING: PROGRAM IS ASKED TO INVERT A SINGULAR MATRIX



先谢谢大家！

koquer

请问有人知道什么原因吗，我也遇到这问题了

yeoo111

不好意思，我没能解决这个，而且不再算这个项目了。