abaqus輸入應(yīng)變的案例
workbench里面輸入材料應(yīng)力應(yīng)變曲線
spm_id_from=333.337.search-card.all.click&vd_source=14f46fb3e57dd52fee0767943481ee71
輸入的曲線要去掉彈性應(yīng)變，為什么最后又增加了一行，是為了保持曲線為水平嗎？
Ls-Dyna MAT24號材料卡片應(yīng)力應(yīng)變曲線輸入以及驗(yàn)證
3、有效應(yīng)力應(yīng)變曲線
從真實(shí)應(yīng)力應(yīng)變曲線中移除彈性段，就得到有效應(yīng)力應(yīng)變曲線，所以要區(qū)分兩條曲線就看第一個(gè)數(shù)據(jù)點(diǎn)的應(yīng)力是否為0。一般而言，有限元軟件要求輸入的都是有效應(yīng)力應(yīng)變曲線，ls-dyna大部分情況下是，abaqus在考慮材料塑性時(shí)，第一行輸入的應(yīng)力應(yīng)變數(shù)據(jù)，應(yīng)變等于0，應(yīng)力不等于0，所以應(yīng)該也是有效應(yīng)力應(yīng)變曲線（屈服應(yīng)力-有效塑性應(yīng)變曲線）。值得注意的是，有效應(yīng)變不等于有效塑性應(yīng)變，兩者的計(jì)算方式是不一樣的，詳情可見https://www.dynasupport.com/howtos/general/effective-plastic-strain。本人才疏學(xué)淺，也并非固體專業(yè)畢業(yè)，很難將這個(gè)問題講清楚，期待有大佬能夠在回答中直觀解釋這個(gè)問題。
【狗骨拉伸試驗(yàn)仿真】
左端固支，右端約束除拉伸方向的所有自由度（由于模型是實(shí)體單元，因此只需要約束平動(dòng)自由度）。
1、應(yīng)力應(yīng)變曲線輸入
真實(shí)應(yīng)力應(yīng)變曲線在有效應(yīng)力應(yīng)變曲線基礎(chǔ)上增加了（0,0）數(shù)據(jù)點(diǎn)，屈服強(qiáng)度均為0.474GPa。
2、有效應(yīng)力-有效塑性應(yīng)變曲線輸出
咦，這個(gè)結(jié)果為啥這么奇怪，各位小伙伴有沒有疑惑，風(fēng)流倜儻、英俊瀟灑、才華橫溢、陽光帥氣的我發(fā)現(xiàn)了一點(diǎn)小問題，有效塑性應(yīng)變接近0的這一段好奇怪，于是我打開了有效應(yīng)力應(yīng)變曲線數(shù)據(jù)，發(fā)現(xiàn)了一點(diǎn)點(diǎn)貓膩。
查看有效應(yīng)力有效塑性應(yīng)變曲線數(shù)據(jù)，有效應(yīng)力一直在增加，可是為啥有效塑性應(yīng)變一直為0呢？經(jīng)過我聰明的小腦袋快速運(yùn)算，很快我想明白了原因。
展開 abaqus子程序VUSDFLD——考慮應(yīng)變率與應(yīng)變軟化效應(yīng)的軟土模型 ￥25
<p><strong>【注意】本貼子只包含子程序文件</strong></p><p>基于<a href="https://www.yqgqt.org.cn/qa/6302" rel="noopener noreferrer" target="_blank">abaqus子程序</a>VUSDFLD編寫的由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?
展開 ABAQUS中求解某部分單元的平均應(yīng)力或平均應(yīng)變 ￥10
1、參考模型：單向纖維的RVE模型；
2、腳本功能：針對指定的單元集合，在后處理中求解平均應(yīng)力和平均應(yīng)變。
3、應(yīng)用的公式：一階均勻化計(jì)算方法。對于 RVE 模型的平均真應(yīng)力和平均真應(yīng)變，可通過對 RVE 內(nèi)每一個(gè)單元的真應(yīng)力 （真應(yīng)變）取均值獲得。使用一階均勻化計(jì)算方法輸出的應(yīng)力和應(yīng)變適用于各種邊界條件，但需要對每個(gè)單元進(jìn)行應(yīng)力（應(yīng)變）的輸出和計(jì)算。
abaqus用戶材料參數(shù)輸入
abaqus安裝目錄下有一個(gè)文件夾，其中有一行代碼是控制用戶材料模塊中參數(shù)輸入的問題。
有誰知道這行代碼在哪里嗎？
有限元理論基礎(chǔ)及Abaqus內(nèi)部實(shí)現(xiàn)方式研究系列27： Abaqus內(nèi)部計(jì)算和顯示的應(yīng)變
Abaqus真正計(jì)算的變量度量可以通過它的子程序的輸入參數(shù)獲取。在Abaqus中，增量步即代表時(shí)刻點(diǎn)，可以查看增量點(diǎn)時(shí)刻的子程序輸入來猜測Abaqus的內(nèi)部量描述方式。UMAT子程序中，在材料本構(gòu)函數(shù)中要利用應(yīng)變增量和當(dāng)前應(yīng)力等物理量更新應(yīng)力，查看UMAT等子程序的接口：
可知其中STRAN和DSTRAN分別表示當(dāng)前增量步最后時(shí)刻的應(yīng)變全量和增量。具體的介紹也可參考下面視頻講解：
http://www.yqgqt.org.cn/college/video/c13034
變形率D在一維上代表對數(shù)應(yīng)變的導(dǎo)數(shù)，但三維上并不是對數(shù)應(yīng)變的導(dǎo)數(shù)，這是有很大區(qū)別的，同時(shí)，可以利用iSolver分別采用上述兩種應(yīng)變度量和Abaqus子程序接口的結(jié)果比對來確認(rèn)Abaqus計(jì)算的應(yīng)變是哪種度量。所以下面我們將找一個(gè)體單元和一個(gè)殼單元的例子來驗(yàn)證到底Abaqus計(jì)算和顯示的應(yīng)變是什么。
1.4 體單元的例子
1.4.1 算例介紹
體單元算例參數(shù)如下：
尺寸：5X1X0.1。
材料：Young’s Modulus 1e8， Poisson Ratio 0.3。
左側(cè)四個(gè)節(jié)點(diǎn)固支。
右側(cè)四個(gè)節(jié)點(diǎn)約束位移為5,1,1。
劃分為一個(gè)殼單元C3D8R。
幾何非線性開關(guān)NLGeom=On，且控制只迭代一次。
1.4.2 Abaqus的應(yīng)變
Abaqus中采用殼的UMAT子程序進(jìn)行計(jì)算。
展開 基于python編程操作ABAQUS輸入文件生成PD3D單元顆粒
在空間中生成剛性顆粒（注意是剛性顆粒）有下列幾種方法：
1.修改關(guān)鍵字，構(gòu)建粒子生成器模型生成隨機(jī)分布剛性顆粒
2.使用python語言直接在ABAQUS中生成顆粒，并進(jìn)行剛體綁定，使其成為剛性顆粒，或者直接生成解析剛體或離散剛體。
方法1生成顆粒的隨機(jī)性較好，操作簡單。方法2直接在ABAQUS界面生成顆粒，當(dāng)所需顆粒數(shù)量以萬為計(jì)量單位時(shí)，在前處理界面時(shí)就會(huì)卡死，對顯卡要求極高。因此，在僅考慮到這些弊端情況下，就已經(jīng)使研究人員頭皮發(fā)麻，無從下手。
在一些特定應(yīng)用場合下，比如所需顆粒數(shù)量數(shù)以萬計(jì)，我們只能采用方法1生成顆粒，但我們不僅僅是需要顆粒，還需將這些顆粒與其它模型進(jìn)行耦合求解計(jì)算，這個(gè)時(shí)候粒子生成器就會(huì)有局限性。此外，考慮到顆粒在空間中排布的多樣性，比如最典型的高斯分布，那么粒子生成器很難做到一步到位生成所需分布特征的顆粒。
我們今天介紹的通過python編程操作ABAQUS輸入文件生成PD3D單元顆粒，其可操作性更強(qiáng)，我們可以不采用粒子生成器內(nèi)部定義的隨機(jī)算法生成顆粒，用戶可以根據(jù)需求自定義顆粒分布算法，以契合實(shí)際工況。此外，可省去粒子生成顆粒的分析步，直接進(jìn)行工況建模求解計(jì)算。
本貼只是個(gè)人興趣，只提供思路，不提供源碼，用戶需了解ABAQUS的inp文件的書寫規(guī)則、python操作文件語法和生成顆粒的底層邏輯(分布模型)。感興趣的可以私信，提供編寫思路。
下面我們采用這一方法生成直徑2mm、3mm、4mm和5mm的混合顆粒，數(shù)量為1000。具體生成結(jié)果如下圖所示。
展開 c30~c45塑性損傷本構(gòu)ABAQUS輸入 ￥1
快速查詢混凝土彈模和拉壓強(qiáng)度標(biāo)準(zhǔn)值.xlsx
以下內(nèi)容也可在我發(fā)布的文檔中免費(fèi)下載
ABAQUS批量提交Job與Python讀取ODB結(jié)果應(yīng)力應(yīng)變數(shù)據(jù)
批量提交的核心是需要等待當(dāng)前計(jì)算任務(wù)的結(jié)束，上圖中若干個(gè)job的提交代碼如下：
# coding: utf-8
#微信公眾號：ABAQUS二次開發(fā)
#作者：阿信老師CAE
#email：axin_cae@163.com
#2022.03.17
from abaqus import *
from abaqusConstants import *
jobList = mdb.jobs.keys() #獲取所有計(jì)算任務(wù)的jobName
for jobName in jobList:
mdb.jobs[jobName].submit(consistencyChecking=OFF) #提交計(jì)算
mdb.jobs[jobName].waitForCompletion() #等待計(jì)算完成
print jobName , "is completed"
不過在實(shí)際的操作過程中，我們可能會(huì)需要避開一些job不提交，或者job窗口事先并沒有建立job，等等，總之實(shí)際問題永遠(yuǎn)比任何教程都復(fù)雜，不過只需要靈活面對就可以了，處理起來并不難。
展開 abaqus鋼管混凝土的本構(gòu)關(guān)系輸入
各位大佬，請問abaqus里面關(guān)于鋼管混凝土的核心混凝土的本構(gòu)關(guān)系怎么輸入？是按照一般的分成彈性和塑性分別輸入嗎，但是彈性模量不知道怎么取。本構(gòu)關(guān)系是根據(jù)韓林海的鋼管混凝土結(jié)構(gòu)計(jì)算的。
混凝土CDP模型，直接輸入ABAQUS中，建議方便 ￥20
本excle簡捷易懂，只需在excle表中更改彈模以及軸心抗壓強(qiáng)度自動(dòng)生成數(shù)據(jù)，表中列出了公式以及只需要輸入ABAQUS中的數(shù)據(jù)，十分容易上手
混凝土CDP模型.xlsx
abaqus Python二次開發(fā)之 交互輸入和提示框
#單輸入框 提示框 from abaqus import getInput from math import sqrt number = float(getInput('Enter a number:')) print sqrt(number) #多輸入提示框 from abaqus import getInputs fields = (('Width:','10'), ('Length:', '20'), ('Height:', '30')) length, width, height = getInputs(fields=fields, label='Specify block dimensions:', dialogTitle='Create Block', ) print length, width, height #警告提示框 from abaqus import getWarningReply, YES, NO reply = getWarningReply(message='Okay to continue?', buttons=(YES,NO)) if reply == YES: print 'YES clicked' elif reply == NO: print 'NO clicked'
展開 Abaqus通過USDFLD子程序進(jìn)行泥巖的應(yīng)變軟化模擬
本文基于泥巖的三軸壓縮試驗(yàn)曲線，建立考慮應(yīng)變軟化特性的泥巖彈塑性本構(gòu)模型，使用Abaqus及其子程序?qū)δ鄮r的三軸壓縮試驗(yàn)進(jìn)行了數(shù)值模擬。泥巖在受壓過程中主要經(jīng)過了5個(gè)階段，即壓密階段、彈性變 形、應(yīng)變硬化、應(yīng)變軟化、殘余階段。泥巖應(yīng)變軟化模型如下所示。
式中，ξ為強(qiáng)度參數(shù)，ξp為峰值強(qiáng)度參數(shù)，ξr為殘余階段強(qiáng)度參數(shù)，η為應(yīng)變軟化參數(shù)，η*為殘余階段的應(yīng)變軟化參數(shù)初始值。對于三軸壓縮試驗(yàn)，η用塑性剪切應(yīng)變來表示
塑性屈服準(zhǔn)則采用Mohr-Coulomb準(zhǔn)則，則粘聚力和內(nèi)摩擦角的參數(shù)演化可以用下式表示
在巖石的塑性變形過程中會(huì)產(chǎn)生比較明顯的剪脹現(xiàn)象，而用來描述這一現(xiàn)象的較常用的力學(xué)參數(shù)就是剪脹角 Ψ，
上述模型可以通過USDFLD子程序進(jìn)行實(shí)現(xiàn)，流程圖如下
有限元模型如下圖所示
計(jì)算得到的應(yīng)力云圖及不同圍壓下的載荷位移響應(yīng)如下所示
參考文獻(xiàn)：張力偉,賈善坡,鄒江濤,舒婧曦.泥巖的峰后軟化力學(xué)模型.中國科技論文,2016,11(21):2456-2461
有關(guān)于abaqus子程序開發(fā)的相關(guān)問題可以通過公眾號聯(lián)系我們.
公眾號: 320科技工作室
展開 hypermesh與abaqus聯(lián)合仿真：模型導(dǎo)入（單component多次輸入）
一、hypermesh模型導(dǎo)出
單個(gè)component（hypermesh）對應(yīng)part（abaqus），利用inp文件做中間文件。
在操作前注意：
將需要的網(wǎng)格放入hypermesh中的一個(gè)component中，比如本例中只需要三維網(wǎng)格，但有兩個(gè)component，因此將對應(yīng)的三維網(wǎng)格放入對應(yīng)的component中，刪除多余的網(wǎng)格和幾何信息。
Abaqus的命名規(guī)則有特殊要求，因此在文件導(dǎo)出hypermesh前需要檢查component的命名，注意命名不能有小數(shù)點(diǎn)、開頭必須是英文字符、不能是abaqus關(guān)鍵詞。
做好設(shè)置后按如下步驟操作：
二、abaqus導(dǎo)入
按component順序依次導(dǎo)入模型，成為abaqus中的part。
File—import-model 對應(yīng)的inp文件。
依次導(dǎo)入后，導(dǎo)入幾次就產(chǎn)生多少個(gè)model，然后按以下操作將模型進(jìn)行規(guī)整。完成操作后所有parts就匯集到了最終的一個(gè)model中
展開 ABAQUS混凝土損傷塑性模型損傷因子對本構(gòu)關(guān)系影響 附c40~c45混凝土損傷因子ABAQUS輸入
模擬云圖結(jié)果如下：
混凝土棱柱體模擬結(jié)果的應(yīng)力-應(yīng)變曲線如下：
從圖中我們可以看出：
1，在單調(diào)荷載下，設(shè)置損傷因子后混凝土峰值應(yīng)力和峰值應(yīng)變比無損傷因子的要低，但差別只有3%左右，可以忽略不計(jì)。
2，循環(huán)荷載設(shè)置為3圈，綠色曲線為無損傷因子的，藍(lán)色曲線為有損傷因子的。可以看到第一圈上升段兩曲線重合，峰值應(yīng)力相同，因?yàn)榇藭r(shí)混凝土還沒有損傷。第一圈卸荷時(shí)，無損傷因子的曲線峰值應(yīng)變更大。第二圈和第三圈的兩曲線混凝土峰值應(yīng)力均大大降低。但無損傷因子的混凝土峰值應(yīng)力和退化后的彈性模量比有損傷因子的要大。
因此，損傷因子在單調(diào)荷載下設(shè)置與否差別不大，但非常影響循環(huán)荷載下的混凝土應(yīng)力-應(yīng)變關(guān)系，損傷因子會(huì)降低循環(huán)荷載下?lián)p傷后的混凝土強(qiáng)度和剛度。
此外
損傷因子是可以通過計(jì)算確定的，計(jì)算方法有多種，文獻(xiàn)【1】中給出的下圖計(jì)算方法是其中的一種。
下載地址：c40~c45混凝土損傷因子ABAQUS輸入
展開 