
發(fā)布
注冊(cè)
/
登錄VUSDFLD
關(guān)注創(chuàng)建者:復(fù)合材料力學(xué)-君莫 創(chuàng)建時(shí)間:2016-10-08
VUSDFLD的視頻教程
VUSDFLD在模擬斷裂中的簡(jiǎn)單應(yīng)用
本課程介紹了VUSDFLD在模擬斷裂中的簡(jiǎn)單應(yīng)用,附件提供了視頻中模型的inp文件和子程序,里面也對(duì)使用方法有文字說(shuō)明,歡迎大家交流討論!
¥9.9 3分鐘 58播放
查看
abaqus二次開(kāi)發(fā)-vusdfld單元?jiǎng)h除
采用Fortran語(yǔ)言編寫(xiě)vusdfld,并通過(guò)計(jì)算單元場(chǎng)輸出,控制單元?jiǎng)h除。
¥10 13分鐘 317播放
查看
VUSDFLD的實(shí)例教程
ABAQUS顯示動(dòng)力學(xué)中,通過(guò)子程序的方式實(shí)現(xiàn)土體強(qiáng)度隨深度線性增加的幾種方法:
第一種:VUFIELD子程序
第二種:VUSDFLD子程序
第三種:VUFIELD和VUSDFLD子程序結(jié)合。
需要說(shuō)明的是:第三種方法,首先通過(guò)VUFIELD子程序建立場(chǎng)變量與坐標(biāo)的關(guān)聯(lián),然后通過(guò)VUSDFLD將VUFIELD定義的場(chǎng)變量保存為狀態(tài)變量,然后更新VUSDFLD定義的場(chǎng)變量。VUSDFLD無(wú)需通過(guò)CALL的方式調(diào)用VUFIELD的場(chǎng)變量。
以上僅為個(gè)人的理解,如有錯(cuò)誤,請(qǐng)各位不吝指出。
vufield_vusdfld.for
vusdfld.for
vufield.for
展開(kāi) Abaqus運(yùn)用Vusdfld實(shí)現(xiàn)自定義單元?jiǎng)h除
———公眾號(hào)‘CAE仿真實(shí)驗(yàn)室’出品
Abaqus提供了幾種單元?jiǎng)h除的方法,如下表所示,單元?jiǎng)h除的本質(zhì)是讓單元喪失承受應(yīng)力的能力,使其對(duì)模型剛度沒(méi)有貢獻(xiàn)。在顯式分析中,可以利用*depvar,delete實(shí)現(xiàn)單元?jiǎng)h除,該方法是通過(guò)給單元設(shè)定一個(gè)狀態(tài)變量來(lái)控制單元是否刪除,當(dāng)狀態(tài)變量為0時(shí),單元?jiǎng)h除,狀態(tài)變量為1時(shí),單元不刪除。今天就來(lái)介紹一些通過(guò)VUSDFLD子程序來(lái)自定義單元?jiǎng)h除準(zhǔn)則的例子
1、例子
單軸拉伸樣條,左端固定,右端給定一個(gè)速度拉伸,假設(shè)當(dāng)應(yīng)變超過(guò)0.02時(shí),單元就失效,刪除單元。為了讓中間某些單元更快達(dá)到失效應(yīng)變,事先在式樣中間刪除幾個(gè)單元,使中間應(yīng)變更大。
定義材料模型時(shí),設(shè)置Depvar和控制element deletion的狀態(tài)變量并設(shè)置User Defined Field.Vusdfld關(guān)鍵代碼如下,關(guān)于vusdfld用法參考之前有關(guān)usdfld的帖子。
2、結(jié)果
應(yīng)變結(jié)果如下,開(kāi)始時(shí)應(yīng)變均勻增大,達(dá)到失效應(yīng)變后,中間單元應(yīng)變極劇增大,說(shuō)明該單元基本沒(méi)有‘剛度’。
應(yīng)變?cè)茍D和應(yīng)變曲線
單元?jiǎng)h除后,單元仍然在,但是單元的應(yīng)力為0,右端拉力在失效前均勻增大,失效后,急劇下降,之后結(jié)果波動(dòng),此時(shí)結(jié)果已經(jīng)沒(méi)有意義~
應(yīng)力云圖和應(yīng)力曲線
3、結(jié)論
由上面結(jié)果可以看出,基本實(shí)現(xiàn)了單元?jiǎng)h除的功能,但是該例子較為簡(jiǎn)單,用于復(fù)雜模型還沒(méi)有驗(yàn)證。
展開(kāi) 利用VUSDFLD子程序可以實(shí)現(xiàn)變形過(guò)程中單元狀態(tài)的控制
用戶手冊(cè)中關(guān)于VUSDFLD的表述為:
VUSDFLD:即USER DEFINED FIELD,用戶自定義場(chǎng)變量,可以實(shí)現(xiàn)單元的刪除。
通常由兩種實(shí)現(xiàn)方式,在CAE界面材料屬性后直接加入場(chǎng)變量,或者編寫(xiě)用戶子程序VUSDFLD,如下圖:
其中用戶子程序的接口結(jié)構(gòu)以及變量含義如下:
其中前面的變量為與CAE交互的只讀變量,用戶需要定義的變量為
stateNew,以及field。
通常適用VUSDFLD時(shí)需要應(yīng)用另一個(gè)子函數(shù),即vgetvrm,來(lái)獲取積分點(diǎn)信息。接口形式為:
其中VAR是待獲取變量。并將返回值給予對(duì)應(yīng)的變量,如下:
常用于獲取的變量如:
S:所有的應(yīng)力分量
MISES:mises等效應(yīng)力
TRIAX:應(yīng)力三軸度
LODE:洛德角參數(shù)
PEEQ:等效塑性應(yīng)變
SDEG:折減剛度系數(shù)
TEMP:積分點(diǎn)溫度
需要注意的是:
(1)該子函數(shù)不能應(yīng)用于獲取用戶定義的狀態(tài)變量。同時(shí)注意2D與3D獲取的分量順序:
(2)需要在材料界面打開(kāi)用戶自定義場(chǎng)
(3)做單元?jiǎng)h除時(shí)也要指明利用哪個(gè)狀態(tài)變量表示材料狀態(tài),哪個(gè)狀態(tài)變量用于判斷是否符合應(yīng)力刪除的特征。從而實(shí)現(xiàn)單元的刪除。
本案例介紹如下:
1,模型幾何尺寸20*20mm包含中心區(qū)域直徑為2.5mm的缺口
2,使用包含1514個(gè)CPE4R單元,每個(gè)單元包含代表一個(gè)單獨(dú)的晶粒
3,分別固定X0和Y0在x,y方向的自由度,并施加X(jué)1方向10%的工程應(yīng)變
4,設(shè)置單元的最大Mises等效應(yīng)力為100Mpa,超過(guò)100Mpa后單元失效刪除。
結(jié)果如下:
初始時(shí)刻的應(yīng)力狀態(tài)
單元?jiǎng)h除效果展示
展開(kāi) 在編寫(xiě)vusdfld\vumat等子程序的時(shí)候,有時(shí)需要知道單元編號(hào),以便對(duì)單元直接賦值,不同于usdfld,umat子程序,他們的單元編號(hào)信息會(huì)傳入子程序中。但是顯示動(dòng)力學(xué)對(duì)應(yīng)的子程序傳入的是積分點(diǎn)編號(hào),并沒(méi)有單元編號(hào)。可以有一種取巧的方法來(lái)實(shí)現(xiàn)積分點(diǎn)和單元編號(hào)的對(duì)應(yīng)。具體思路如下:
以邊長(zhǎng)為50mm立方體為例,網(wǎng)格種子密度為5,共1000個(gè)單元,單元類型C3R8R。而在顯示動(dòng)力學(xué)中,使用單核運(yùn)算時(shí),abaqus會(huì)一次性傳入136個(gè)積分點(diǎn),1000個(gè)單元會(huì)傳入8次,最后一次為48個(gè)積分點(diǎn),可以通過(guò)write(*,*) ‘nblock’,nblock輸出來(lái)看傳入積分點(diǎn)的變化情況。因此可以在程序中設(shè)置一個(gè)全局變量,每8次循環(huán)后就歸0,重新循環(huán)。利用狀態(tài)變量statenew(k,1)來(lái)存儲(chǔ)單元。
在后處理中顯示單元編號(hào),可以看出云圖從1到1000,正好對(duì)應(yīng)所有單元的個(gè)數(shù),利用查詢值功能,可以看到單元ID和狀態(tài)變量的值是相同的。
知道單元編號(hào)之后就可以對(duì)單元進(jìn)行賦值操作,比如讀取一個(gè)隨機(jī)材料場(chǎng)文件,然后將數(shù)據(jù)賦值給單元。還可以利用該功能,實(shí)現(xiàn)類似“生死單元”的功能。一開(kāi)始就根據(jù)單元編號(hào)抑制一部分,在之后的分析步在激活。比如第一個(gè)分析步殺死一些單元,然后下一個(gè)分析步在激活,比如刪除編號(hào)為501的單元。
可以看到編號(hào)501的狀態(tài)變量為值為0,單元應(yīng)力為0,說(shuō)明該單元被刪除,不參與運(yùn)算過(guò)程。如果在場(chǎng)輸出中把STATUS勾選上,可以視覺(jué)上顯示出來(lái)單元?jiǎng)h除的效果。結(jié)果如下圖所示:
需要注意的是如果一個(gè)單元有多個(gè)積分點(diǎn),則循環(huán)中需要跳過(guò)一部分積分點(diǎn)。該方法需要單核計(jì)算才能使用。
最后有子程序開(kāi)發(fā)等相關(guān)需求,歡迎聯(lián)系我們。
展開(kāi) <p><strong>【注意】本貼子只包含子程序文件</strong></p><p>基于<a href="https://www.yqgqt.org.cn/qa/6302" rel="noopener noreferrer" target="_blank">abaqus子程序</a>VUSDFLD編寫(xiě)的由Einav與Randolph提出的西澳模型,用于求解軟黏土體劇烈變形后的強(qiáng)度變化,可應(yīng)用于的大變形計(jì)算。</p><div contenteditable="false" width="100%"><img src="https://img.jishulink.com/upload/202306/f69f50d42a81489ea1cb5e7a03da5c14.png" title="8$U(VZ82]O{OEMQB}[P(ZMB.png" alt="8$U(VZ82]O{OEMQB}[P(ZMB.png" style="max-width:760px;" data-mobile-src="https://img.jishulink.com/upload/202306/f69f50d42a81489ea1cb5e7a03da5c14.png?image_process=/format,webp" data-pc-src="https://img.jishulink.com/upload/202306/f69f50d42a81489ea1cb5e7a03da5c14.png?
展開(kāi) 
VUSDFLD的相關(guān)專題、標(biāo)簽、搜索
VUSDFLD的最新內(nèi)容
ABAQUS顯示動(dòng)力學(xué)中,通過(guò)子程序的方式實(shí)現(xiàn)土體強(qiáng)度隨深度線性增加的幾種方法:
第一種:VUFIELD子程序
第二種:VUSDFLD子程序
第三種:VUFIELD和VUSDFLD子程序結(jié)合。
需要說(shuō)明的是:第三種方法,首先通過(guò)VUFIELD子程序建立場(chǎng)變量與坐標(biāo)的關(guān)聯(lián),然后通過(guò)VUSDFLD將VUFIELD定義的場(chǎng)變量保存為狀態(tài)變量,然后更新VUSDFLD定義的場(chǎng)變量。
關(guān)鍵詞:圍巖;裂隙;MC準(zhǔn)則;VUSDFLD;二次開(kāi)發(fā);單元?jiǎng)h除算法;
1 引言
水工隧洞地質(zhì)情況復(fù)制,部分隧洞圍巖內(nèi)部含有天然缺陷,這些缺陷以裂隙的方式存在。在開(kāi)挖擾動(dòng)作用下,裂隙擴(kuò)展成宏觀裂縫,導(dǎo)致圍巖失穩(wěn)引發(fā)安全事故。因此,對(duì)隧洞圍巖裂隙擴(kuò)展規(guī)律展開(kāi)研究具有重要意義。本文以某引水隧洞為例,結(jié)合隧洞圍巖實(shí)際力學(xué)屬性,將MC準(zhǔn)則嵌入ABAQUS子程序VUSDFLD進(jìn)行二次開(kāi)發(fā)。
法的單樁基礎(chǔ)貫入過(guò)程模擬:考慮應(yīng)變軟化與應(yīng)變率效應(yīng)_abaqus cel實(shí)例 ABAQUS二次開(kāi)發(fā)-技術(shù)鄰</a></p><p>在abaqus軟件中基于<a href="https://www.yqgqt.org.cn/service/abaqus_cel" rel="noopener noreferrer" target="_blank">CEL</a>法的分層地基單樁基礎(chǔ)貫入過(guò)程模擬,通過(guò)編寫(xiě)VUSDFLD
USDFLD(Abaqus/標(biāo)準(zhǔn))和 VUSDFLD(Abaqus/Explicit)
USDFLD子程序允許用戶定義自定義場(chǎng)變量以用于邊界條件、初始化和輸出等。這些場(chǎng)變量可以是位移、應(yīng)力、應(yīng)變、溫度等,并且可以隨時(shí)間和位置而變化。這些場(chǎng)變量可以在模擬過(guò)程中更新,以便更好地模擬分析中的復(fù)雜現(xiàn)象。
4.
圖3 READ函數(shù)使用方法
隨后將數(shù)組A中數(shù)據(jù)存入狀態(tài)變量?jī)?nèi),即可實(shí)現(xiàn)晶粒度的繼承,如圖4所示為第二道次初始時(shí)刻的晶粒度分布,與圖2所示第一道次結(jié)束時(shí)刻晶粒度分布完全一致:
圖4 第二道次初始時(shí)刻晶粒度分布云圖
四、小結(jié)
本文所使用的方法可以方便地實(shí)現(xiàn)任意狀態(tài)變量在不同變形道次之間的相互繼承,對(duì)于VUMAT,VUSDFLD
<p><strong>【注意】本貼子只包含子程序文件</strong></p><p>基于<a href="https://www.yqgqt.org.cn/qa/6302" rel="noopener noreferrer" target="_blank">abaqus子程序</a>VUSDFLD編寫(xiě)的由Einav與Randolph提出的西澳模型,用于求解軟黏土體劇烈變形后的強(qiáng)度變化,可應(yīng)用于的大變形計(jì)算
此文檔為VUSDFLD子程序如何在CAE中激活的步驟詳解,希望可以為有需要的朋友帶來(lái)幫助!如果有不正之處也請(qǐng)大家批評(píng)指正(新手小白的瑟瑟發(fā)抖)。
發(fā)現(xiàn)了一些問(wèn)題,請(qǐng)查看最新版的文件!!
案例以Abaqus/Explicit內(nèi)置的JC本構(gòu)模型,以及編寫(xiě)的自定義硬化VUSDFLD和自定義材料VUMAT三種方式模擬taylor桿撞擊,材料使用鋼材,考慮熱膨脹,熱傳導(dǎo),和塑性產(chǎn)熱以及JC本構(gòu)模型,模擬使用顯式動(dòng)態(tài)溫度-位移分析步,使用C3D8RT單元,初始溫度為20°,沖擊速度為287000mm/s,沖擊時(shí)間為8e-5s,其中使用內(nèi)置本構(gòu),VUSDFLD以及VUMAT的結(jié)果如下圖所示:
