ABAQUS在inp里添加彈簧單元全過程

2017-04-20  by:CAE仿真在線  來源:互聯(lián)網(wǎng)

模擬鋼筋混凝土粘結(jié)滑移,以我現(xiàn)在的水平來看,大概有三個思路,一個是用cohesive,一個是用Spring2彈簧單元,再有就是Umat了。我先學(xué)習(xí)的是彈簧單元。

ABAQUS里有三種彈簧單元,Spring1,Spring2和SpringA。其中,ABAQUS文檔的解釋如下:

ABAQUS在inp里添加彈簧單元全過程abaqus有限元技術(shù)圖片1


Spring2有兩個節(jié)點(diǎn),每個節(jié)點(diǎn)只有一個自由度。鋼筋采用Wire建模,Truss單元,混凝土實(shí)體單元。mesh之后會出節(jié)點(diǎn),所以盡量mesh成相同的節(jié)點(diǎn),用鋼筋和混凝土重合的節(jié)點(diǎn),在三個方向分別建立彈簧。垂直于鋼筋縱向的兩個彈簧設(shè)置成剛度較大,彈性。沿鋼筋方向的彈簧,需要確定粘結(jié)滑移本構(gòu),定義力和滑移位移的關(guān)系。對比一下:Spring1是彈性的,只有一個節(jié)點(diǎn),而SpringA的作用是沿兩點(diǎn)連線的,每個點(diǎn)有三個自由度。

線形的彈簧單元可以在Interaction模塊里,找到Special里面添加,非線性的就只能的通過修改inp文件了。內(nèi)容添加Assembly里,我的習(xí)慣是寫在End Assembly前一行。

對于添加內(nèi)容,大概有兩部分,一部分是對節(jié)點(diǎn)的定義,另一部分是對力-位移曲線的定義。兩者是分開的,通過elset識別。elset是Element set 的意思,就是定義一個單元集合,名字根據(jù)在abaqus里添加線性彈簧單元所生成的inp來看,是集合名字后邊跟"-spring"type=Spring2是固定格式,

ABAQUS在inp里添加彈簧單元全過程abaqus有限元技術(shù)圖片2

這句后邊是節(jié)點(diǎn)編號,如下圖。


ABAQUS在inp里添加彈簧單元全過程abaqus有限元仿真圖片3


值得注意的是,Instance名是Part名后邊跟編號形成的(如 “-1”),如果Part的名字有空格,就需要加引號了。如下圖

ABAQUS在inp里添加彈簧單元全過程abaqus有限元仿真圖片4

最左邊一列的數(shù)字是單元的編號,從1開始,順著排下去就行。中間和右邊兩列是單元的節(jié)點(diǎn),Wire-1是Instance的名字。查看單元編號和節(jié)點(diǎn)編號,在View--Part Display Options里,在mesh一欄分別勾選Show node labels和Show element labels。在mesh模塊的視圖里就可以顯示了。實(shí)體名稱和對應(yīng)的節(jié)點(diǎn)編號之間用“.”連接。

ABAQUS在inp里添加彈簧單元全過程abaqus有限元仿真圖片5


線性彈簧,一般給個大剛度就好,用于定義型鋼和混凝土接觸面發(fā)現(xiàn)法向的接觸。

ABAQUS在inp里添加彈簧單元全過程abaqus有限元仿真圖片6

非線性彈簧,根據(jù)粘結(jié)滑移本構(gòu)得來。

關(guān)于粘結(jié)滑移本構(gòu),公式很多,相差也非常大,這里因?yàn)橐?guī)范給出的最好算,還是采取用規(guī)范的值。對于型鋼混凝土,建議看看楊勇的論文。

ABAQUS在inp里添加彈簧單元全過程abaqus有限元培訓(xùn)資料圖片7


ABAQUS在inp里添加彈簧單元全過程abaqus有限元培訓(xùn)資料圖片8

得到切向的應(yīng)力之后,再根據(jù)網(wǎng)格的劃分,確定了每個單元的鋼筋的長度之后,可以得到F。

ABAQUS在inp里添加彈簧單元全過程abaqus有限元培訓(xùn)資料圖片9

就是應(yīng)力乘以面積。

*Spring, elset=X-Spring(非線性的Spring2,多一個Nonlinear,這個是很多教程都是這么說的,我測試后發(fā)現(xiàn),去掉之后也沒有關(guān)系)。之后定義彈簧的方向,X軸方向?yàn)?,1;Y軸方向?yàn)?.2;Z軸方向?yàn)?,3。兩個數(shù)字中,第一個數(shù)字的含義是節(jié)點(diǎn)1沿著1軸(即X軸)運(yùn)動,第二個數(shù)字的含義是節(jié)點(diǎn)2沿著1軸(即X軸)方向運(yùn)動;F=kX,X就是兩點(diǎn)在1軸方向運(yùn)動的位移差。如果前后數(shù)字不同,我理解的,就是向量的加減了。

這些定義好之后,inp就可以計算了。我這邊有一個完整的inp文件,類似于型鋼混凝土的小練習(xí),有需要看這個inp的,在回復(fù)中留下郵箱,我會發(fā)送到郵箱里。

大多數(shù)教程講到這里就沒有然后了。后邊的一些東西,我感覺還是挺關(guān)鍵的。兩點(diǎn),一個是inp文件做好后,不要導(dǎo)入abaqus再計算,因?yàn)閍baqus識別不了;二是直接在job里提交inp的方法。

打開abaqus,不用新建任何項目,直接到j(luò)ob模塊,新建一個job,點(diǎn)集creat之后,Source選擇Input File,如下圖:

ABAQUS在inp里添加彈簧單元全過程abaqus有限元培訓(xùn)資料圖片10

如圖所示,此時就可以計算了。

總結(jié)起來,用inp添加彈簧單元這個過程非常繁瑣,工作量巨大,幸運(yùn)的是,很多單元節(jié)點(diǎn)是有規(guī)律的,這樣,這個過程以后可以靠程序來實(shí)現(xiàn)。


開放分享:優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才

相關(guān)標(biāo)簽搜索:ABAQUS在inp里添加彈簧單元全過程 abaqus分析培訓(xùn) abaqus技術(shù)教程 abaqus巖土分析 鋼筋混凝土仿真 abaqus分析理論 abaqus軟件下載 abaqus umat用戶子程序編程 Abaqus代做 Abaqus基礎(chǔ)知識 Fluent、CFX流體分析 HFSS電磁分析 Ansys培訓(xùn) 

編輯
在線報名:
  • 客服在線請直接聯(lián)系我們的客服,您也可以通過下面的方式進(jìn)行在線報名,我們會及時給您回復(fù)電話,謝謝!
驗(yàn)證碼

全國服務(wù)熱線

1358-032-9919

廣州公司:
廣州市環(huán)市中路306號金鷹大廈3800
電話:13580329919
          135-8032-9919
培訓(xùn)QQ咨詢:點(diǎn)擊咨詢 點(diǎn)擊咨詢
項目QQ咨詢:點(diǎn)擊咨詢
email:kf@1cae.com