abaqus積分位置的案例
ABAQUS直裂紋、斜裂紋圍道積分計算裂紋尖端J積分
之前算過一個關(guān)于裂紋擴(kuò)展的問題，當(dāng)時創(chuàng)建裂紋選擇的是contour intergral，后來又有人咨詢我裂紋尖端J積分的計算問題。我才恍然大悟，其實(shí)圍道積分方法還是適用于計算裂紋尖端在某時刻的J積分，至于動態(tài)擴(kuò)展問題，還是交給XFEM吧（雖然也不太好）。
計算了幾種情況下的裂紋尖端J積分，包括直裂紋、斜裂紋以及裂紋尖端傾斜等三種情況。
部分試件的應(yīng)力分布及J積分結(jié)果如圖所示：
ABAQUS輸出單元積分點(diǎn)坐標(biāo)
方法
在ABAQUS CAE的場輸出中選擇的坐標(biāo)點(diǎn)是節(jié)點(diǎn)的坐標(biāo)，而節(jié)點(diǎn)是從積分點(diǎn)插值出來的，單元積分點(diǎn)的信息相對真實(shí)。所以最好是獲取積分點(diǎn)的信息，其中積分點(diǎn)的坐標(biāo)無法在CAE中獲取，需要在關(guān)鍵字中添加。具體在每個分析步的單元輸出下面添加COORD，如果需要輸出節(jié)點(diǎn)的坐標(biāo)也可以在節(jié)點(diǎn)場輸出下面添加COORD（這和CAE中場輸出選擇節(jié)點(diǎn)坐標(biāo)的效果是一致的）。具體如下圖：
2.注意
在ODB結(jié)果中創(chuàng)建場輸出時會附帶著一份XYZ坐標(biāo)，這個應(yīng)該也可以當(dāng)做單元的坐標(biāo)，，但是我比較過這個附帶的坐標(biāo)和單元的COORD輸出的坐標(biāo)，有時候有點(diǎn)差別，可能是數(shù)據(jù)精度的問題。
展開 ABAQUS 輸出節(jié)點(diǎn)坐標(biāo)和積分點(diǎn)坐標(biāo)
總結(jié)inp中添加關(guān)鍵字
輸出單元的積分點(diǎn)坐標(biāo)：*EL FILE
COORD
輸出節(jié)點(diǎn)坐標(biāo)：*NODE FILE
COORD
原貼出處：https://www.researchgate.net/post/How-to-find-integration-point-coordinates-in-Abaqus-CAE
這是帖子討論的，但是我的嘗試是兩個COORD生成的結(jié)果文件是一樣的，都是節(jié)點(diǎn)坐標(biāo)
Abaqus中平面應(yīng)力單元高斯積分點(diǎn)的順序
可以輸出umat接口中的變量coords進(jìn)行查看
write(*,"(A,I4)") "npt = ", npt
write(*,"(A,3ES16.8)") "coords = ", coords
結(jié)果為：
npt = 1
coords = -5.77350269E-01 -5.77350269E-01 1.00000000E-02
npt = 2
coords = 5.77350269E-01 -5.77350269E-01 1.00000000E-02
npt = 3
coords = -5.77350269E-01 5.77350269E-01 1.00000000E-02
npt = 4
coords = 5.77350269E-01 5.77350269E-01 1.00000000E-02
因此Abaqus中平面應(yīng)力單元高斯積分點(diǎn)的順序?yàn)椋?/span>
展開 
Python提取Abaqus SPH結(jié)果粒子位置 ￥9.99
Abaqus軟件后處理中SPH粒子的渲染真的太差了，粒子大小、光照等都沒有有效的調(diào)整方法，于是想著從ODB文件中讀取出粒子位置信息，提取出來的信息還可以用于其他后處理及渲染等。
但是有很多問題：1. PC3D粒子的信息怎么輸出？與節(jié)點(diǎn)關(guān)系如何對應(yīng)？2. 可以在界面上選擇節(jié)點(diǎn)編號查詢節(jié)點(diǎn)，但是批量怎么讀取？3. fieldOutput中物理量只有位移等，沒有當(dāng)前坐標(biāo)。
想了一個笨辦法：
先在后處理模塊中用“Display Group”功能顯示所有流體粒子，然后利用查詢功能獲得所有當(dāng)前粒子的label，存到一個列表中；
在odb文件的assembly對象中獲得上面粒子label編號列表對應(yīng)的節(jié)點(diǎn)的坐標(biāo)，就是初始坐標(biāo)。
在steps.frames對象中獲得位移場變量，位移是一個容器，里面存儲了所有節(jié)點(diǎn)的位移，找出流體對應(yīng)節(jié)點(diǎn)的位移，位移與初始坐標(biāo)的和就是當(dāng)前坐標(biāo)。
將查詢到的結(jié)果（包括節(jié)點(diǎn)當(dāng)前坐標(biāo)、節(jié)點(diǎn)位移值等）按照節(jié)點(diǎn)順序?qū)懭雟tk文件，利用paraview進(jìn)行可視化。
如下為讀取及寫入vtk文件的Python代碼
展開 Abaqus中獲取積分點(diǎn)坐標(biāo)的三種方法
經(jīng)常有小伙伴問獲取積分點(diǎn)坐標(biāo)的方法，今天給大家介紹三種獲取積分點(diǎn)坐標(biāo)的方式，希望能給你們帶來幫助。
1 通過abaqus子程序獲取積分點(diǎn)坐標(biāo)
Abaqus一些子程序中可以直接獲取積分點(diǎn)坐標(biāo)，例如我們熟知的UMAT子程序中包含COORD參數(shù)，即為積分點(diǎn)坐標(biāo)。順帶一提的是，當(dāng)打開了幾何非線性時，該積分點(diǎn)是當(dāng)前構(gòu)形下的坐標(biāo)，如果未打開幾何非線性則為初始坐標(biāo)。
2通過history output輸出積分點(diǎn)坐標(biāo)
Abaqus可以直接在歷程變量history output中輸出積分點(diǎn)坐標(biāo)。直接在history output中勾選COORD選項，但是這里需要注意的是，Domain中的Set集合如果是node set，這里輸出來的是節(jié)點(diǎn)坐標(biāo)，當(dāng)這里是element set的時候，輸出來的才是積分點(diǎn)坐標(biāo)。
3通過等參單元映射函數(shù)計算
等參元中，為了方便計算，把整體坐標(biāo)映射到自然坐標(biāo)，然后在自然坐標(biāo)下進(jìn)行高斯積分。如果知道了自然坐標(biāo)下的高斯積分點(diǎn)，通過映射函數(shù)反算，便能得到整體坐標(biāo)下的高斯積分點(diǎn)坐標(biāo)。以四邊形等參單元為例，其以自然呢坐標(biāo)為變量的插值形函數(shù)如下
坐標(biāo)變換采取同樣的插值函數(shù)(叫做等參的原因)，整體坐標(biāo)和自然坐標(biāo)的關(guān)系式如下，如果知道自然坐標(biāo)下的高斯積分點(diǎn)，直接通過此公式計算其在整體坐標(biāo)下的坐標(biāo)。
展開 abaqus高斯熱源中心點(diǎn)的位置
摘要：abaqus的高斯熱源網(wǎng)上有很多例題，能夠運(yùn)行，但是并沒有講的太詳細(xì)。我用自己的模型，稍作修改就發(fā)現(xiàn)加載的位置不對了，所以來研究一下熱源的中心位置（x0,y0,z0）的定義方法。這里使用surface flux進(jìn)行研究。
test 1：加載面選擇XOY平面，x0=0，y0=0。建立part時，長方形的一個角為坐標(biāo)原點(diǎn)。
test 2：加載面選擇XOY平面，x0=0，y0=0.07，y向總長度為0.14
test 3：現(xiàn)在想要熱源從上往下移動，也就是沿著y軸負(fù)方向。除了給定速度為負(fù)值以為，初始位置改為（x0=0.03，y0=0.14）。x0.03時為了查看結(jié)果方便，y向總長度為0.14
test 4：之前的測試都有一個容易被忽略的前提，我們建模的時候長方形的左下角為坐標(biāo)原點(diǎn)，重新建立一個模型，使得長方形最下面一條邊的中點(diǎn)為草圖的坐標(biāo)原點(diǎn)。
結(jié)論：熱源施加的初始位置和兩個因素有關(guān)
1、建模的時候草圖的原點(diǎn)
2、子程序中的坐標(biāo)x0，y0。這個點(diǎn)是相對于草圖中的原點(diǎn)的位置。也就是說當(dāng)草圖坐標(biāo)原點(diǎn)在模型之外時，選擇（x0=0，y0=0）時看不到加載效果的。
展開 Abaqus降魔篇之圍道積分小探
1 分析目的
平面應(yīng)變的一塊板，在扭矩的作用下，分析裂紋端的應(yīng)力強(qiáng)度因子以及斷裂能
模擬采用CPE4R（使用CPE8R往往會得到更精確的結(jié)果）
2 模型建立
然后進(jìn)行切分，主要是切出預(yù)制裂紋面以及圍道外邊
在這里注意切那個圓的時候一定要將半徑點(diǎn)選到裂紋上，以便劃分網(wǎng)格
這里的圓半徑為0.5mm
3 建立材料以及相關(guān)屬性
4 建立模型，預(yù)制裂紋
選擇之前的切分線
選擇裂紋端點(diǎn)
這里的0.25是根據(jù)劃分網(wǎng)格的種子點(diǎn)確定的（4個點(diǎn)）
而collapsed的選項是要求對重疊單元的裂紋控制只對單一點(diǎn)有效
5 設(shè)置邊界條件
a首先是施加扭矩端的耦合
b然后是扭矩施加
另一端同理可得
c 限制位移
6 計算吧！
然后你可以看到應(yīng)力強(qiáng)度因子以及斷裂能
或者你比較喜歡看趨勢
在或者你對歷史輸出里面的5感到困惑
那就看下圖
展開 abaqus 中的 Johnson-Cook 模型如何控制損傷起始位置？
2 JC本構(gòu)——損傷演化段
The Johnson-Cook criterion (available only in Abaqus/Explicit) is a special case of the ductile criterion in which the equivalent plastic strain at the onset of damage, , is assumed to be of the form
上面的英文是幫助文檔中對于 Johnson-Cook 損傷準(zhǔn)則的解釋，具體意思：下面的公式是定義損傷起始/萌生時的等效塑性應(yīng)變，當(dāng)達(dá)到損傷起始等效塑性應(yīng)變，材料就會發(fā)生損傷。
式中，d1-d5是需要輸入的損傷參數(shù)，損傷演化段和塑性硬化段一樣，等號右側(cè)第二個括號與第三個括號分別是應(yīng)變率和溫度對于損傷的影響。
下面將討論修改參考應(yīng)變率對于損傷起始位置的影響：
把參考應(yīng)變率從4e-4修改成1，損傷的起始位置會從圖1右邊黃色框住的位置變成左邊，因?yàn)閰⒖紤?yīng)變率變大第二個括號變小，導(dǎo)致?lián)p傷起始等效塑性應(yīng)變變小，即損傷位置提前
圖1 修改參考應(yīng)變率對于損傷起始位置的影響
參考資料：
（1）TC4鈦合金動態(tài)力學(xué)性能及本構(gòu)模型研究_惠旭龍
（2）abaqus 幫助文檔
展開 abaqus怎么獲取單元積分點(diǎn)上的坐標(biāo)呢
如題 為什么我查詢的時候坐標(biāo)都是0
abaqus C3D8 單元 計算中采用了多少個積分點(diǎn)？
按照正常的理解，毫無.疑問，abaqus 全積分一定是采用了2x2x2=8個積分點(diǎn)。
從后處理結(jié)果來看，似乎也是如此，每個單元存在8個積分點(diǎn)。
然而，如果自己動手跑一遍程序，就會發(fā)現(xiàn)事實(shí)遠(yuǎn)非如此，采用全積分計算得到的結(jié)果與abaqus 存在差異，原因何在？
事實(shí)賞，abaqus C3D8 采用的選擇積分方式（selective intergation schema），即對于偏應(yīng)變，采用8個積分，對于球應(yīng)變，采用中心點(diǎn)積分。這樣計算得到的結(jié)果才能與abaqus 完全對標(biāo)，亦可從abaqus 幫助文檔得到答案。
展開 
有限元中單元積分點(diǎn)與節(jié)點(diǎn)應(yīng)力相互轉(zhuǎn)換(CPE4為例)（ABAQUS）
（注意：變量是a,b,c,d,而不是x,y.所以方程組是線性的）
第一個積分點(diǎn)的應(yīng)力和坐標(biāo)：S1，(X1，Y1);
第二個積分點(diǎn)的應(yīng)力和坐標(biāo)：S2，(X2，Y2);
第三個積分點(diǎn)的應(yīng)力和坐標(biāo)：S3，(X3，Y3);
第四個積分點(diǎn)的應(yīng)力和坐標(biāo)：S4，(X4，Y4);
現(xiàn)在的問題是：應(yīng)力分量S1，S2，S3，S4是已知的，我們需要知道真實(shí)的積分點(diǎn)的坐標(biāo)信息嗎？
答案：不需要，只需要知道積分點(diǎn)在整個單元相對位置即可。即等參元中的坐標(biāo)。（教材中有）
等參元的長和寬都為2.
而有限元中的積分是高斯積分，積分點(diǎn)的位置是固定的。由查表可知：
上表是一維的高斯積分點(diǎn)的坐標(biāo)，后面的加權(quán)系數(shù)不用管（我們不求積分）。由一維可以猜出二維（兩個一維）。二維有4個積分點(diǎn)，所以我們對應(yīng)一維選第二行的數(shù)據(jù)。
展開 Abaqus|智能優(yōu)化算法的反演加強(qiáng)筋位置，提升薄壁結(jié)構(gòu)穩(wěn)定性 ￥50
代碼是通過Python腳本來實(shí)現(xiàn)，其代碼主要包含三個模塊runAbaqus、main_DE_inverse、main_TS_inverse,分別代表執(zhí)行CAE計算、差分進(jìn)化算法反演和遍歷搜索算法反演。
