《ABAQUS输出轴力和弯矩.docx》由会员分享，可在线阅读，更多相关《ABAQUS输出轴力和弯矩.docx（5页珍藏版）》请在三一文库上搜索。
1、ABAQUS中如何通过 cutting surface 和 section print 输由桩的 轴力经过两个星期的摸索与学习，今天终于学会了桩轴力的输出。现总结如下：1 .主要步骤是先定义截面cutting surface ,然后用 section print输出轴力 sof。2 .所有操作均是在inp文件中进行修改的，而不是ABAQUS/CAE 中的编辑关键词(edit keywords ) 原因：在CAE中编辑关键词是可以修改inp文件，但CAE并不能识别所有的inp文件关键词，下面将举例说明。3 .最后提交的inp文件也不是在 CAE中导入模型文件(import model ),然后提
2、交job进行运算 的，而是在 ABAQUS命令窗口(小黑屏)进行的。原因同2中的一样，CAE并不能识别关键词*section print。好了，下面开始详细的步骤讲解吧！第一步：定义截面(cutting surface ),具体的关键语句为：*surface,type=cutting surface,name=cutsurf-10.6,25,0,0,1,0Set-pile解读：第一行，定义surface surface类型以及名称。第二行，定义截面上的一点(0.6,25,0)以及截面的法向量 (0,1,0)。法向量不一定是单位向量。第三行，截面所在的单元或集合。这个集合可以是事先在CAE里定义
3、好的。需要注意的是，这个 cutting surface是垂直于桩径方向的一个桩截面，而不是桩的侧表面。我一开始理解错了！此关键句在inp文件中的位置是在Assembly, name=Assembly这一行之后，我试了下在放在材 料定义之后，运算不成功。估计是因为我的Set-pile是在assembly里定义的，而不是在part里定义的单元集。若是在 part里定义了桩的集合，是不是可以将此关键句移动到材料定义之后，这 个我倒没试。如果是在CAE中通过编辑关键词来添加上述语句，将会有下面的错误提示：第二步：定义输出(section print ),具体的关键语句为：section print,
4、name=secprint-1,surface=cutsurf-1sof,som解读：第一行，定义输出的名称及截面。第二行，定义输出选项，sof为截面合力，som为截面弯矩。注意，此处第一行 surface=*必须是上一步中 cutting surface的名称(name=* ),因为要输出的就是第一步定义的截面轴力及弯矩。如果此处乱写surface= ?的话，估计到最后什么也输出不了 了。有些人还在第一句后面加上了axes=local,frequency=1,update=yes ,由于我的模型较简单，也就没管这几个东西。除了命名之外，inp文件中的大小写英文字符，是没有区别的。所以大家不要
5、再犹豫写*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 i npu
6、t file reader:*SECTIONPRINTThe 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文件进
7、行运算了。具体操作是：开始-所有程序-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后缀为好。
8、到此，计算完毕。黑屏窗口并不会像 CAE那样提示job completed或ANSYS中的solution is done !（害得我还以为没反应-二汗！）默认工作路径下已经生产了结果文件 .com .dat .log .msg .odb .prt .sta 等一系列后缀文件。第四步：提取数据用记事本打开刚才生产的.dat文件，在最下面会找到想要的输出结果THE FOLLOWING TABLES ARE PRINTEDSOFMSOF1SOF2SOF31.0556E+064.0515E+041.0548E+060.000三个方向的sof,以及合力SOFM。需要注意的是，此 sof值是合应力，单位
9、为 No无需乘以截面面积至此，大功告成！后话：（1）本人查看了仿真论坛里的一些帖子，发现输出轴力的方法有三：后处理模块中运用free body cut此法较为方便简单，可直接得到截面合力,无奈我的模型是二维轴对称的，不能进行自由体切割Displaying a free body cutYou can define a free body cut to view the resultant forces and moments transmitted across a selected surface of a model. Force vectors are displayed with a
10、single arrowhead and moment vectors with a double arrowhead.To create a free body cut:1. To display the entire model in the viewport, select Tools I-D Display Group -Plot _*All from the main menu bar.2. From the main menu bar, select Tools I* Free Body Cut _* Manager .3. Click Create in the Free Bod
11、y Cut Manager .4. From the dialog box that appears, select 3D element faces as the Selection method and click Continue .5. In the Free Body Cross-Section dialog box, select Surfaces as the Item and Pick from viewport as the Method .6. In the prompt area, set the selection method to by angle and acce
12、pt the default angle.a.7. Select the surface, highlighted in Figure 4 -33, to define the free body cut cross-section.From the Selection toolbar, toggle off the Select the Entity Closestto the Screen tool and ensure that the Select From AllEntities toolis selected.b. As you move the cursor in the vie
13、wport, Abaqus/CAE highlights all of the potential selections and adds ellipsis marks (.) next to the cursor arrow to indicate an ambiguous selection. Position the cursor so that one of the faces of the desired surface is highlighted, and click to display the first surface selection.Figure 4 -33 Sele
14、cted faces for the free body cross-section.c. Use the Next and Previous buttons to cycle through the possible selections until the appropriate vertical surface is highlighted, and click OK.8. Click Done in the prompt area to indicate your selection is complete. Click OK in the Free Body Cross-Sectio
15、n dialog box.9. In the Edit Free Body Cut dialog box, accept the default settings forthe Summation Point and the Component Resolution . Click OK to close the dialog box.10. Click Options in the Free Body Cut Manager .11. From the Free Body Plot Options dialog box, select the Force tab intoOK inthe C
16、olor & Style tabbed page. Click the resultant color sample change the color of the resultant force arrow.12. Once you have selected a new color for the resultant force arrow, clickin the Free BodyFigure 4 T4.the Free Body Plot Options dialog box and click Dismiss Cut Manager .The free body cut is di
17、splayed in the viewport, as shown inFigure 4 -34 Free body cut displayed on the connecting lug.输出截面上的应力值进行积分运算，主要方法是创建路径(path)或者拾取节点应力值，桩身的网格越密集，最后的积分值越精确，此法我已学会。cutting surface ft section print联合使用，由于手动修改inp文件，此法最 复杂，而且不知道最终的sof,软件是如何算出来的。如果需要得到桩深方向N个截面的轴力，就需要在cutting surface里定义N个截面。(2)特别感谢仿真论坛的 zs
18、q-w总斑竹、shanhuimin923 斑竹、firce 、laoyao239以及 ABAQUS- 岩土群(84418357)里的羽出惊人 (274550165)。ABAQUS6.10支持cae定义截面，这样不需要在inp里定义截面三要素(截面上任意一个点、截面法线、截面所在的集 合Set),从而方便很多。同时，cae支持定义多个截面，在inp里添加多条输出语句，最后都能在 dat文件中找到sof和som。值得一提的是，后处理free body cut同样可以输出之前cae定义好的截面上的弯矩和轴力。经验证，和sof、som输出值一模一样。方法一：* *定义切平面*surface,type=
19、cutting surface,name=fifth-ring-up* *定义输出一一（SOM为截面总弯矩，SOF为截面合力）*section print,name=forcemoment,surface=fifth-ring-up,update=yesSOM, SOF注意：关键词 surface 和 section print 的用法可以参见 abaqus 帮助 /Abaqus Keywords Referrence Manual 以及 Abaqus/CAE Users Manual 等。方法二：首先需要导出节点应力，节点应力可以：（1）从Visualization/Report/Field output 里选择Unique Nodal导出；（2）由于节点应力不能直接在inp中关键词输出，但是可以在inp文件中将节点附近单元应力导出，然后取平均值即可；语法为：* el print, elset=fifth-ring-outside-els将截面内外节点的应力值求出后，用材料力学方法求解。