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

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

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

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

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

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

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

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

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