ABAQUS中t輸出樁的軸力講課稿

1、Good is good, but better carries it.精益求精，善益求善。ABAQUS中t輸出樁的軸力-經(jīng)過兩個星期的摸索與學(xué)習(xí)，今天終于學(xué)會了樁軸力的輸出?，F(xiàn)總結(jié)如下：1.主要步驟是先定義截面cuttingsurface，然后用sectionprint輸出軸力sof。2.所有操作均是在inp文件中進行修改的，而不是ABAQUS/CAE中的編輯關(guān)鍵詞（editkeywords）。原因：在CAE中編輯關(guān)鍵詞是可以修改inp文件，但CAE并不能識別所有的inp文件關(guān)鍵詞，下面將舉例說明。3.最后提交的inp文件也不是在CAE中導(dǎo)入模型文件（importmodel），然后提交job

2、進行運算的，而是在ABAQUS命令窗口（小黑屏）進行的。原因同2中的一樣，CAE并不能識別關(guān)鍵詞*sectionprint。好了，下面開始詳細的步驟講解吧！第一步：定義截面（cuttingsurface），具體的關(guān)鍵語句為：*surface,type=cuttingsurface,name=cutsurf-10.6,25,0,0,1,0Set-pile解讀：第一行，定義surface、surface類型以及名稱。第二行，定義截面上的一點（0.6,25,0）以及截面的法向量（0,1,0）。法向量不一定是單位向量。第三行，截面所在的單元或集合。這個集合可以是事先在CAE里定義好的。需要注意的是，這

3、個cuttingsurface是垂直于樁徑方向的一個樁截面，而不是樁的側(cè)表面。我一開始理解錯了！此關(guān)鍵句在inp文件中的位置是在*Assembly,name=Assembly這一行之后，我試了下在放在材料定義之后，運算不成功。估計是因為我的Set-pile是在assembly里定義的，而不是在part里定義的單元集。若是在part里定義了樁的集合，是不是可以將此關(guān)鍵句移動到材料定義之后，這個我倒沒試。如果是在CAE中通過編輯關(guān)鍵詞來添加上述語句，將會有下面的錯誤提示：第二步：定義輸出（sectionprint），具體的關(guān)鍵語句為：*sectionprint,name=secprint-1,su

4、rface=cutsurf-1sof,som解讀：第一行，定義輸出的名稱及截面。第二行，定義輸出選項，sof為截面合力，som為截面彎矩。注意，此處第一行surface=*必須是上一步中cuttingsurface的名稱（name=*），因為要輸出的就是第一步定義的截面軸力及彎矩。如果此處亂寫surface=？的話，估計到最后什么也輸出不了了。有些人還在第一句后面加上了axes=local,frequency=1,update=yes，由于我的模型較簡單，也就沒管這幾個東西。除了命名之外，inp文件中的大小寫英文字符，是沒有區(qū)別的。所以大家不要再猶豫寫*SECTIONPRINT還是*secti

5、onprint此關(guān)鍵句在inp文件中的位置是最后一個*HISTORYOUTPUT:H-Output-3之后，*Output,history,variable=PRESELECT之前。也就是我的inp文件中的倒數(shù)第四和第三行。如果放錯了位置，估計也會引起不識別等問題，最終導(dǎo)致無輸出結(jié)果。如果是在CAE中通過導(dǎo)入inp文件來進行運算，在導(dǎo)入時，CAE下方窗口將會有下述警告信息：WARNING:Thefollowingkeywords/parametersarenotyetsupportedbytheinputfilereader:-*SECTIONPRINTThemodelJob-20101215

6、hasbeenimportedfromaninputfile.Pleasescrolluptocheckforerrorandwarningmessages.這再次驗證了，ABAQUS/CAE是不能識別關(guān)鍵詞*sectionprint的。強制提交job是可以運算的，但是在最后的dat結(jié)果文件里沒用任何輸出結(jié)果。第三步：保存并提交inp文件。在保存之前需要仔細檢查上面關(guān)鍵語句的位置及拼寫，以及是否含有中文字符。然后就可以在ABAQUS命令窗口提交inp文件進行運算了。具體操作是：開始所有程序Abaqus6.9.1AbaqusCommand輸入：abaqusjob=*.inp（就是你剛才修改的in

7、p文件名）回車即可。然后黑屏窗口會彈出一個警告語句：AbaqusWarning:The.inpextensionhasbeenremovedfromthejobindentifier大概意思是說，后綴.inp被省略。由此可判斷，在上面輸入的時候，完全可以不寫后綴名字。我之所以寫了，是因為我的默認工作路徑D:Temp下還有個跟inp文件完全同名的.cae文件，而如果這兩個文件信息不是完全一樣的話，建議還是加上.inp后綴為好。到此，計算完畢。黑屏窗口并不會像CAE那樣提示jobcompleted或ANSYS中的solutionisdone！（害得我還以為沒反應(yīng)-_-汗?。┠J工作路徑下已經(jīng)生產(chǎn)了

8、結(jié)果文件.com.dat.log.msg.odb.prt.sta等一系列后綴文件。第四步：提取數(shù)據(jù)用記事本打開剛才生產(chǎn)的.dat文件，在最下面會找到想要的輸出結(jié)果THEFOLLOWINGTABLESAREPRINTEDSOFMSOF1SOF2SOF31.0556E+064.0515E+041.0548E+060.000三個方向的sof，以及合力SOFM。需要注意的是，此sof值是合應(yīng)力，單位為N。無需乘以截面面積至此，大功告成！后話：（1）本人查看了仿真論壇里的一些帖子，發(fā)現(xiàn)輸出軸力的方法有三：后處理模塊中運用freebodycut，此法較為方便簡單，可直接得到截面合力，無奈我的模型是二維軸對稱的，不能進行自由體切割。輸出截面上的應(yīng)力值進行積分運算，主要方法是創(chuàng)建路徑（path）或者拾取節(jié)點應(yīng)力值，樁身的網(wǎng)格越密集，最后的積分值