重启动问题

zhao12

我想先用牛顿-拉夫孙发进行计算，将近荷载极限值后再中断，然后再重启用弧长法算出下阶段。我是按照下面步骤进行的：
中断前所用的牛顿-拉夫孙法的计算文件为job-a，其中step如下：
** STEP: Step-1
**
*Step, name=Step-1, nlgeom=YES
*Static
0.1, 1., 1e-05, 1.
**
** LOADS
**
** Name: Load-1   Type: Shell edge load
*Dsload, op=NEW, constant resultant=YES
Surf-1, EDSHR, 600000.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=2
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step


重启动分析的iNP文件如下（文件名为job-2）：
*Heading
** Job name: Job-2 Model name: Model-1
** Generated by: Abaqus/CAE Version 6.8-1
*Preprint, echo=NO, model=NO, history=NO, contact=NO
*Restart, read, step=1, inc=10
** STEP: Step-2
**
** STEP: Step-2
**
*Step, name=Step-2, nlgeom=YES
*Static, riks
0.1, 1., 1e-05, , 1.,
**
** LOADS
**
** Name: Load-2   Type: Shell edge load
*Dsload, op=NEW
Surf-1, EDSHR, 420000.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=2
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step

然后我从command中输入
abaqus job=job-2 oldjob=job-a
结果发现，job-2仅仅是将job-a继续算下去（继续按照牛顿-拉夫孙法算到极值），却没有采用弧长发算出下降段。
请问各位大侠这是什么原因吗？我哪里出错了？谢谢！！

cet6china

*Restart, read, step=1, inc=10
这句里面加个END STEP参数试试。

zhao12

感谢cet6china 兄弟，我先试试看

zhao12

感谢cet6china 兄弟，已经出现弧长发了，不过我对后处理结果方面还有点疑惑。
我的INP文件修改按照如下：
*Heading
** Job name: Job-2 Model name: Model-1
** Generated by: Abaqus/CAE Version 6.8-1
*Preprint, echo=NO, model=NO, history=NO, contact=NO
*Restart, read, step=1,inc=20,END STEP
** STEP: Step-2
**
** STEP: Step-2
**
*Step, name=Step-2, nlgeom=YES
*Static, riks
0.1, 1., 1e-05, , 1.,
**
** LOADS
**
** Name: Load-1   Type: Shell edge load
*Dsload, op=NEW, constant resultant=YES
Surf-1, EDSHR, 600000.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step



我先输出某一点的位移，其两列数据如下：
0.679605        0.225407
0.779605        0.226152
0.879605        0.226894
0.979605        0.22765
1.07961        0.228418
1.17961        0.229202
1.27961        0.230002
。。。。。。。
第二列应该是位移值，第一列是在前面老的计算结果最后一步荷载步长0.679605的基础上增加弧长步，如弧长0.1增量相当于0.779605。
但是我再提取LPF（荷载比例系数）——在纯粹用弧长发计算结果中，这个数据第一列是弧长增量，第二列是荷载因子（乘以总荷载即为每增量一步荷载大小）。在这个时候这两列数据如下：
0        0
0.1        0.000248262
0.2        0.000469881
0.3        0.000672768
0.4        0.000865116
0.5        0.00104275
0.6        0.00118997
0.7        0.00131604
0.8        0.00140403
0.9        0.00144549
1        0.00143023
1.1        0.00136855
1.2        0.00127206
1.35        0.00101177
1.575        0.000397554
1.9125        -0.00110396
2.41875        -0.0049196
3.17812        -0.0161467
4.31719        -0.0589503
5.45625        -0.142872
6.59531        -0.231752
7.73438        -0.305404
8.87344        -0.369149
10.0125        -0.427127
11.1516        -0.479031
12.8602        -0.550572
14.5688        -0.615306
16.2773        -0.672382
17.9859        -0.724315
19.6945        -0.771348
21.4031        -0.814039
23.1117        -0.854057
24.8203        -0.892279
26.5289        -0.929486
28.2375        -0.964732
29.9461        -0.997861
31.6547        -1.02966
33.3633        -1.06056
35.0719        -1.09082
36.7805        -1.12013
38.4891        -1.14852
第一列看上去像是弧长因子，第二列怎么都不像是荷载因子，难道第二列数据加上0.679605（此数为第一次用牛顿-拉夫孙法计算后，重启用弧长法计算时的荷载因子）。但好像也不像啊。因为荷载最后变成负的了？
请教各位大侠，怎么才能找出荷载因子系数？