abaqus/explici显式分析中t加速计算方法——子循环Subcycling

taishanbuzuo

在abaqus/explicit显式求解器中，也有子循环的方法，能够提高显示分析的效率。据说可以提高20%的计算速度，没比较过，但是确实快了。
原理是：单元尺度小的单元采用较小的时间增量，而单元尺度大的单元则采用较大的时间增量。如果不采用子循环，则左右单元都按照最小的时间增量进行计算。
详见abaqus analysis user's manual  11.8.1 Selective subcycling.
使用方法，在*step 关键字之前最后一行，增加如下关键字：

*Elset, elset=subcycle, instance=plate-mesh-1
   833,   834,   835,   836,   837,   838,   839,   840,   841,   842,   843,   844,   845,   846,   847,   848


*Subcycling, ELSET=subcycle
**
** STEP: Step-1
**
*Step, name=Step-1

其中，subcycle是自己建好的子循环单元集，里面存储的是模型中直接影响稳定时间增量的单元。


开始计算之后，在status file里面会提示，子循环内单元的稳定时间增量和子循环外的稳定时间增量：
Most critical elements in subcycling zones:
    Element number   Rank    Time increment   Increment ratio
    (Instance name)
   ----------------------------------------------------------
        2340          1       1.593618E-08      1.000000E+00
    PLATE-MESH-1
        6372          2       1.593618E-08      1.000000E+00
    PLATE-MESH-1
        9060          3       1.593618E-08      1.000000E+00
    PLATE-MESH-1
       13092          4       1.593618E-08      1.000000E+00
    PLATE-MESH-1
       15108          5       1.593618E-08      1.000000E+00
    PLATE-MESH-1
       19140          6       1.593618E-08      1.000000E+00
    PLATE-MESH-1
       21828          7       1.593618E-08      1.000000E+00
    PLATE-MESH-1
       25860          8       1.593618E-08      1.000000E+00
    PLATE-MESH-1
        2341          9       1.593939E-08      9.997984E-01
    PLATE-MESH-1
        6373         10       1.593939E-08      9.997984E-01
    PLATE-MESH-1
   Most critical elements in non-subcycling zone:
    Element number   Rank    Time increment   Increment ratio
    (Instance name)
   ----------------------------------------------------------
         229          1       3.837396E-07      4.152862E-02
    TANK1-1
        2745          2       4.098517E-07      3.888280E-02
    TANK1-1
        2413          3       4.123782E-07      3.864457E-02
    TANK1-1
        2744          4       4.179116E-07      3.813290E-02
    TANK1-1
        2389          5       4.207705E-07      3.787380E-02
    TANK1-1
        3682          6       4.298889E-07      3.707046E-02
    TANK1-1
           1          7       4.301034E-07      3.705197E-02
    JIAOHE-1-1
         545          8       4.301039E-07      3.705193E-02
    JIAOHE-1-1
          67          9       4.301041E-07      3.705192E-02
    JIAOHE-3-1
           3         10       4.301041E-07      3.705192E-02
    JIAOHE-4-1

taishanbuzuo

额   发现一个问题  算到一定时间就出错，还米找到原因。

AsAlways

露珠啊，有没有找到问题的说？话说我觉着这个功能很实用啊。目前正在学习。能实用CAE方法来实现吗？

AsAlways

似乎是不能用CAE的方法，不过这个添加起来也还算方便。最烦那种要用编译软件打开INP修改的了。修好了还要另建一个job（或者用黑乎乎的DOS窗口提交）,蛋疼。

lttommy

请问出现如图所示的错误是因为什么？