经过两个星期的摸索与学习，今天终于学会了桩轴力的输出。现总结如下：
1.主要步骤是先定义截面cutting surface，然后用section print输出轴力sof。
2.所有操作均是在inp文件中进行修改的，而不是ABAQUS/CAE中的编辑关键词（edit keywords）。
原因：在CAE中编辑关键词是可以修改inp文件，但CAE并不能识别所有的inp文件关键词，下面将举例说明。
3.最后提交的inp文件也不是在CAE中导入模型文件（import model），然后提交job进行运算的，而是在ABAQUS命令窗口（小黑屏）进行的。
原因同2中的一样，CAE并不能识别关键词*section print。
好了，下面开始详细的步骤讲解吧！
第一步：定义截面（cutting surface），具体的关键语句为：
*surface,type=cutting surface,name=cutsurf-1
0.6,25,0,0,1,0
Set-pile
解读：
第一行，定义surface、surface类型以及名称。
第二行，定义截面上的一点（0.6,25,0）以及截面的法向量（0,1,0）。法向量不一定是单位向量。
第三行，截面所在的单元或集合。这个集合可以是事先在CAE里定义好的。
需要注意的是，这个cutting surface是垂直于桩径方向的一个桩截面，而不是桩的侧表面。我一开始理解错了！
此关键句在inp文件中的位置是在*Assembly, name=Assembly这一行之后，我试了下在放在材料定义之后，运算不成功。估计是因为我的Set-pile是在assembly里定义的，而不是在part里定义的单元集。若是在part里定义了桩的集合，是不是可以将此关键句移动到材料定义之后，这个我倒没试。
如果是在CAE中通过编辑关键词来添加上述语句，将会有下面的错误提示：
第二步：定义输出（section print），具体的关键语句为：
*section print,name=secprint-1,surface=cutsurf-1
sof,som
解读：
第一行，定义输出的名称及截面。
第二行，定义输出选项，sof为截面合力，som为截面弯矩。
注意，此处第一行surface=**必须是上一步中cutting surface的名称（name=**），因为要输出的就是第一步定义的截面轴力及弯矩。如果此处乱写surface=？的话，估计到最后什么也输出不了了。
有些人还在第一句后面加上了axes=local,frequency=1,update=yes，由于我的模型较简单，也就没管这几个东西。除了命名之外，inp文件中的大小写英文字符，是没有区别的。所以大家不要再犹豫写*SECTION PRINT还是*section print……
此关键句在inp文件中的位置是最后一个** HISTORY OUTPUT: H-Output-3之后，*Output, history, variable=PRESELECT之前。也就是我的inp文件中的倒数第四和第三行。如果放错了位置，估计也会引起不识别等问题，最终导致无输出结果。
如果是在CAE中通过导入inp文件来进行运算，在导入时，CAE下方窗口将会有下述警告信息：
WARNING: The following keywords/parameters are not yet supported by the input file reader:
---------------------------------------------------------------------------------
*SECTIONPRINT
The model "Job-20101215" has been imported from an input file.
Please scroll up to check for error and warning messages.
这再次验证了，ABAQUS/CAE是不能识别关键词*section print的。强制提交job是可以运算的，但是在最后的dat结果文件里没用任何输出结果。
第三步：保存并提交inp文件。
在保存之前需要仔细检查上面关键语句的位置及拼写，以及是否含有中文字符。
然后就可以在ABAQUS命令窗口提交inp文件进行运算了。具体操作是：开始→所有程序→Abaqus 6.9.1→Abaqus Command
输入：
abaqus job=***.inp（就是你刚才修改的inp文件名）
回车即可。
然后黑屏窗口会弹出一个警告语句：
Abaqus Warning:The .inp extension has been removed from the job indentifier
大概意思是说，后缀.inp被省略。由此可判断，在上面输入的时候，完全可以不写后缀名字。我之所以写了，是因为我的默认工作路径D:\Temp下还有个跟inp文件完全同名的.cae文件，而如果这两个文件信息不是完全一样的话，建议还是加上.inp后缀为好。
到此，计算完毕。黑屏窗口并不会像CAE那样提示job completed或ANSYS中的solution is done！（害得我还以为没反应…… -_- 汗！）默认工作路径下已经生产了结果文件 .com .dat .log .msg .odb .prt .sta等一系列后缀文件。
第四步：提取数据
用记事本打开刚才生产的.dat文件，在最下面会找到想要的输出结果
THE FOLLOWING TABLES ARE PRINTED
SOFM SOF1 SOF2 SOF3
1.0556E+06 4.0515E+04 1.0548E+06 0.000
三个方向的sof，以及合力SOFM。
需要注意的是，此sof值是合应力，单位为N。无需乘以截面面积……
至此，大功告成！
后话：
（1）本人查看了仿真论坛里的一些帖子，发现输出轴力的方法有三：
?后处理模块中运用free body cut，此法较为方便简单，可直接得到截面合力，无奈我的模型是二维轴对称的，不能进行自由体切割。
?输出截面上的应力值进行积分运算，主要方法是创建路径（path）或者拾取节点应力值，桩身的网格越密集，最后的积分值越精确，此法我已学会。
?cutting surface和section print联合使用，由于手动修改inp文件，此法最复杂，而且不知道最终的sof值，软件是如何算出来的。如果需要得到桩深方向N个截面的轴力，就需要在cutting surface里定义N个截面。
（2）特别感谢仿真论坛的zsq-w总斑竹、shanhuimin923斑竹、firce、laoyao239以及ABAQUS-岩土群（84418357）里的◆羽出惊人◆(274550165)。
--------------------------------------------------------------------------
如大家发现我的这个总结有什么问题，请不吝赐教！
后来又发了一篇帖子
IP卡
狗仔卡
发表于 2016-3-28 19:34:31
提升卡
沉默卡
喧嚣卡
变色卡
抢沙发
千斤顶
显身卡
