不知火舞的被虐|伊人天伊人天天综合网|博洛尼亚天气|任你懆这里只有精品4|久久美日韩精品久久|掌中之物漫画免费阅读观看|0丨d老妇

VUMAT子程序，只有彈性沒有塑性的原因？（等效應(yīng)變增量一直為0）？

使用CEL計(jì)算水下爆炸的過程中，沖擊波傳遞到結(jié)構(gòu)上之后，結(jié)構(gòu)開始變形并且與結(jié)構(gòu)后方的歐拉域接觸，然后計(jì)算的最小時(shí)間增量步開始減小，直到計(jì)算報(bào)錯(cuò)。懷疑是接觸設(shè)置不合理，CEL的計(jì)算中歐拉域與拉格朗日網(wǎng)格的接觸通常怎么設(shè)置呢？

對(duì)于法向應(yīng)變增量計(jì)算方法就是內(nèi)聚力單元上下表面的法向分離量除以厚度（默認(rèn)為1）；我還查看了幫助手冊(cè)，里面說(shuō)兩個(gè)切向的名義應(yīng)變是積分點(diǎn)的切向位移除以厚度（默認(rèn)為1），但我用這種方法計(jì)算出來(lái)的和abaqus傳入的并不相等。 所以請(qǐng)各位大佬解惑，這個(gè)VUMAT傳入的剪應(yīng)變是怎樣計(jì)算出來(lái)的呢？主要針對(duì)COH3D6零厚度內(nèi)聚力單元

有償

我的理解如下：初始增量步開始時(shí)（時(shí)刻為TotalTime=0），根據(jù)輸入的荷載首先得到一個(gè)初始應(yīng)力數(shù)組，根據(jù)輸入的材料屬性獲得剛度矩陣D，此時(shí)應(yīng)該strainInc=0；初始增量步結(jié)束第二個(gè)增量步開始時(shí)（時(shí)刻為TotalTime=dt），stressOld等于初始應(yīng)力數(shù)組，此時(shí)的strainInc是用初始應(yīng)力數(shù)組和剛度矩陣D計(jì)算得到的嗎？

1：編寫的Vumat已經(jīng)能正常計(jì)算，單單元模擬與abaqus內(nèi)置程序計(jì)算結(jié)果誤差可忽略，2：編寫的Vumat在計(jì)算一些結(jié)構(gòu)性模型，比如沖擊等，計(jì)算很慢。不設(shè)置質(zhì)量縮放的話，abaqus給的初始時(shí)間增量為1e-14，然后顯示總計(jì)算步數(shù)超過2000，0000，最終報(bào)錯(cuò)。設(shè)置質(zhì)量縮放，能計(jì)算，但是仍然很慢。

自己編寫材料vumat，利用平面應(yīng)力單元模擬，但通過輸出應(yīng)變增量發(fā)現(xiàn)strainInc(k,3)在過程中始終等于零，所以最后3方向的應(yīng)變也始終為零。但是利用abaqus自帶材料模型模擬平面應(yīng)力狀態(tài)會(huì)發(fā)現(xiàn)3的應(yīng)變不為零。自己的子程序用三維單元驗(yàn)證是沒問題。請(qǐng)問大佬這是哪里出問題？

微信截圖_20230428152858.png

請(qǐng)問一下，如果需要在VUMAT子程序構(gòu)建滯回曲線，如下圖所示。那么該如何判斷加卸載的轉(zhuǎn)折點(diǎn)1以及轉(zhuǎn)折點(diǎn)2？我看別人實(shí)現(xiàn)的時(shí)候是用應(yīng)變增量小于0時(shí)為卸載，應(yīng)變增量大于0時(shí)為加載。但是在計(jì)算過程中發(fā)現(xiàn)，顯式分析的應(yīng)變會(huì)存在震蕩現(xiàn)象，就算是純加載階段，應(yīng)變增量也會(huì)頻繁出現(xiàn)負(fù)值，好像應(yīng)變的增減是完全隨機(jī)的。即使添加了瑞利阻尼也只是減小了這種震蕩現(xiàn)象，不能消除。