abaqus 高斯積分的案例
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>
展開 有限元中的高斯-勒讓德積分
當(dāng)函數(shù)表達(dá)式比較復(fù)雜時(shí)，f(x)的原函數(shù)可能難以求出，而采用高斯積分，其省去了求f(x)原函數(shù)，只需要將數(shù)值代入f(x)的表達(dá)式即可求解。</p><p>到目前為止，高斯積分的公式已經(jīng)介紹完成，那么有兩個(gè)最直接最現(xiàn)實(shí)的問題出現(xiàn)了：（1）f(x)的表達(dá)式是什么形式時(shí)適合采用高斯積分，精度怎么樣；（2）xi和wi的取值是多少。</p><p>關(guān)于（1），實(shí)踐表明，當(dāng)f（x）的表達(dá)式為多項(xiàng)式時(shí)，高斯積分是合適的，并且，n點(diǎn)高斯積分可以準(zhǔn)確積分2n-1次多項(xiàng)式。</p><p>關(guān)于（2），xi和wi的取值一般較多的有限元教科書中會(huì)給出數(shù)值，如果沒有給出數(shù)值，也可以用多項(xiàng)式手動(dòng)算出具體值，另外，scipy庫，PETSc庫也直接給出了高斯積分的值和權(quán)重。</p><p><img src="https://img.jishulink.com/msimage/202308/f410d6e3bb89c8459660277304592181.png"></p><p>以下是高斯積分點(diǎn)積分多項(xiàng)式的一個(gè)例子：</p><div contenteditable="false" width="100%">
<img src="https://img.jishulink.com/upload/202308/865d98a129374c668e89080010b652c9.jpg" title="圖片4.jpg" alt="圖片4.jpg" style="max-width:760px;" data-mobile-src="https://img.jishulink.com/upload/202308/865d98a129374c668e89080010b652c9.jpg?
展開 等參數(shù)公式描述-四節(jié)點(diǎn)四邊形單元-高斯積分法-python編程 ￥6.99
等參數(shù)公式描述-四節(jié)點(diǎn)四邊形單元-高斯積分法-python編程
有限元中的高斯點(diǎn)與積分詳解(下)_《數(shù)值計(jì)算與程序設(shè)計(jì)》系列課程之九 ￥599
本課從實(shí)際問題出發(fā)，帶著問題去講解有限元中的高斯點(diǎn)與數(shù)值積分。一開始拋出了以下3個(gè)關(guān)鍵問題：
1.對(duì)于一個(gè)任意函數(shù)怎么去得到它的積分？
2.數(shù)值積分的本質(zhì)是什么？為什么簡(jiǎn)單地取幾個(gè)點(diǎn)就可得到積分值？此種方法的立足點(diǎn)在哪？
3.很多資料上都說“有限元求解精度嚴(yán)重依賴于網(wǎng)格質(zhì)量，過度扭曲的單元會(huì)導(dǎo)致結(jié)果不收斂或者精度極度惡化”，這只是為什么呢？扭曲單元到底影響的是有限元方法中的哪一步？
圍繞這3個(gè)問題，本課分別講了一下三個(gè)內(nèi)容：
1. 數(shù)值積分基本方法。
2. 有限元單元積分。
3. 誤差分析。
希望有興趣的同學(xué)多多支持下，你們的支持是我更新的動(dòng)力
展開 
有限元中的高斯點(diǎn)與積分詳解(上)_《數(shù)值計(jì)算與程序設(shè)計(jì)》系列課程之五 ￥599
本課從實(shí)際問題出發(fā)，帶著問題去講解有限元中的高斯點(diǎn)與數(shù)值積分。一開始拋出了以下3個(gè)關(guān)鍵問題：
1.對(duì)于一個(gè)任意函數(shù)怎么去得到它的積分？
2.數(shù)值積分的本質(zhì)是什么？為什么簡(jiǎn)單地取幾個(gè)點(diǎn)就可得到積分值？此種方法的立足點(diǎn)在哪？
3.很多資料上都說“有限元求解精度嚴(yán)重依賴于網(wǎng)格質(zhì)量，過度扭曲的單元會(huì)導(dǎo)致結(jié)果不收斂或者精度極度惡化”，這只是為什么呢？扭曲單元到底影響的是有限元方法中的哪一步？
圍繞這3個(gè)問題，本課分別講了一下三個(gè)內(nèi)容：
1. 數(shù)值積分基本方法。
2. 有限元單元積分。
3. 誤差分析。
本次課程分為上下兩課，第一課講了第一和第二個(gè)內(nèi)容。關(guān)鍵詞是：數(shù)值積分的本質(zhì)，有限元高斯積分（附件中包含1個(gè)小時(shí)的詳細(xì)課程視頻以及PPT）。
在第二課中，再繼續(xù)展開第三部分內(nèi)容，誤差分析，解決問題“扭曲單元到底影響的是有限元方法中的哪一步”。
希望有興趣的同學(xué)多多支持下，你們的支持是我更新的動(dòng)力
展開 abaqus中焊接高斯面熱源和高斯體熱源程序 ￥19.89
abaqus中焊接中高斯面熱源和高斯體熱源程序
ABAQUS直裂紋、斜裂紋圍道積分計(jì)算裂紋尖端J積分
之前算過一個(gè)關(guān)于裂紋擴(kuò)展的問題，當(dāng)時(shí)創(chuàng)建裂紋選擇的是contour intergral，后來又有人咨詢我裂紋尖端J積分的計(jì)算問題。我才恍然大悟，其實(shí)圍道積分方法還是適用于計(jì)算裂紋尖端在某時(shí)刻的J積分，至于動(dòng)態(tài)擴(kuò)展問題，還是交給XFEM吧（雖然也不太好）。
計(jì)算了幾種情況下的裂紋尖端J積分，包括直裂紋、斜裂紋以及裂紋尖端傾斜等三種情況。
部分試件的應(yīng)力分布及J積分結(jié)果如圖所示：
Abaqus高斯體熱源
為什么高斯圓錐體熱源不能在z軸上移動(dòng)呢？
abaqus_子程序_高斯熱源 ￥15
前期準(zhǔn)備(這個(gè)可以百度，有更詳細(xì)的安裝教程)：
首先需要安裝兩個(gè)軟件（InterFortran和vs studio，注意ABAQUS與軟件版本號(hào)的兼容性，不然會(huì)出現(xiàn)配置不成功的情況）。鑒于題主用的是abauqs2016，此處以abauqs2016為例子簡(jiǎn)述步驟。
Abaqus2016的兼容版本為VS2012，InterFortran2013（需要lic許可證，要不然只能用一年_maybe）.
安裝過程中需要記錄各自的安裝地址，軟件安裝完成后，在abaqus中的一個(gè)launcher.bat文件中將上述兩個(gè)文件的路徑寫入該文件（注意路徑要正確，同時(shí)嚴(yán)格按照格式書寫，禁止使用中文字符，禁止多莫名其妙的空格）。
顯示上圖所示即配置成功。
建模（建模過程僅敘述關(guān)鍵點(diǎn)）：
材料屬性：必須輸入用于溫度場(chǎng)計(jì)算的熱物性參數(shù)。將絕對(duì)零度設(shè)定為-273.15，此時(shí)計(jì)算出來的溫度單位為℃。
如果只是單純的溫度場(chǎng)分析,將分析步類型設(shè)定為熱傳遞，分析步時(shí)間按照熱源移動(dòng)速度，模型大小計(jì)算得出。
邊界條件：根據(jù)需求設(shè)定輻射、散熱系數(shù)。荷載選擇表面熱流，作用區(qū)域定位在熱源照射面。預(yù)定義場(chǎng)設(shè)定初始溫度。
網(wǎng)格類型為DC3D8。
提交計(jì)算：
提交作業(yè)，在通用模塊用戶子程序一欄鍵入.for子程序路徑。
嗯，大概就這樣。
付費(fèi)內(nèi)容包括子程序文件，CAE模型，聯(lián)系方式。
展開 abaqus高斯熱源中心點(diǎn)的位置
摘要：abaqus的高斯熱源網(wǎng)上有很多例題，能夠運(yùn)行，但是并沒有講的太詳細(xì)。我用自己的模型，稍作修改就發(fā)現(xiàn)加載的位置不對(duì)了，所以來研究一下熱源的中心位置（x0,y0,z0）的定義方法。這里使用surface flux進(jìn)行研究。
test 1：加載面選擇XOY平面，x0=0，y0=0。建立part時(shí)，長(zhǎng)方形的一個(gè)角為坐標(biāo)原點(diǎn)。
test 2：加載面選擇XOY平面，x0=0，y0=0.07，y向總長(zhǎng)度為0.14
test 3：現(xiàn)在想要熱源從上往下移動(dòng)，也就是沿著y軸負(fù)方向。除了給定速度為負(fù)值以為，初始位置改為（x0=0.03，y0=0.14）。x0.03時(shí)為了查看結(jié)果方便，y向總長(zhǎng)度為0.14
test 4：之前的測(cè)試都有一個(gè)容易被忽略的前提，我們建模的時(shí)候長(zhǎng)方形的左下角為坐標(biāo)原點(diǎn)，重新建立一個(gè)模型，使得長(zhǎng)方形最下面一條邊的中點(diǎn)為草圖的坐標(biāo)原點(diǎn)。
結(jié)論：熱源施加的初始位置和兩個(gè)因素有關(guān)
1、建模的時(shí)候草圖的原點(diǎn)
2、子程序中的坐標(biāo)x0，y0。這個(gè)點(diǎn)是相對(duì)于草圖中的原點(diǎn)的位置。也就是說當(dāng)草圖坐標(biāo)原點(diǎn)在模型之外時(shí)，選擇（x0=0，y0=0）時(shí)看不到加載效果的。
展開 Abaqus高斯熱源3D打印熱力耦合模擬
Abaqus高斯熱源3D打印熱力耦合模擬
ABAQUS輸出單元積分點(diǎn)坐標(biāo)
方法
在ABAQUS CAE的場(chǎng)輸出中選擇的坐標(biāo)點(diǎn)是節(jié)點(diǎn)的坐標(biāo)，而節(jié)點(diǎn)是從積分點(diǎn)插值出來的，單元積分點(diǎn)的信息相對(duì)真實(shí)。所以最好是獲取積分點(diǎn)的信息，其中積分點(diǎn)的坐標(biāo)無法在CAE中獲取，需要在關(guān)鍵字中添加。具體在每個(gè)分析步的單元輸出下面添加COORD，如果需要輸出節(jié)點(diǎn)的坐標(biāo)也可以在節(jié)點(diǎn)場(chǎng)輸出下面添加COORD（這和CAE中場(chǎng)輸出選擇節(jié)點(diǎn)坐標(biāo)的效果是一致的）。具體如下圖：
2.注意
在ODB結(jié)果中創(chuàng)建場(chǎng)輸出時(shí)會(huì)附帶著一份XYZ坐標(biāo)，這個(gè)應(yīng)該也可以當(dāng)做單元的坐標(biāo)，，但是我比較過這個(gè)附帶的坐標(biāo)和單元的COORD輸出的坐標(biāo)，有時(shí)候有點(diǎn)差別，可能是數(shù)據(jù)精度的問題。
展開 Abaqus移動(dòng)加載高斯熱源子程序，適于焊接、3D打印領(lǐng)域 ￥2.5
Abaqus移動(dòng)加載高斯熱源子程序，適于焊接、3D打印領(lǐng)域，fortran編制的子程序見“付費(fèi)后”的附件中。
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
這是帖子討論的，但是我的嘗試是兩個(gè)COORD生成的結(jié)果文件是一樣的，都是節(jié)點(diǎn)坐標(biāo)
ABAQUS焊接模擬-Python編寫移動(dòng)高斯熱源子程序-不帶生死單元
ABAQUS平板對(duì)接-不帶生死單元。Python編寫移動(dòng)高斯熱源子程序（包括高斯面熱源、雙橢球熱源）
模型作如下假設(shè)：材料為各向同性材料，不考慮熔池流動(dòng)及相變影響。
考慮到過來學(xué)習(xí)的大多都是和我一樣的學(xué)生黨，因此設(shè)置了一個(gè)大家都能接受的價(jià)格。
如果視頻中有什么錯(cuò)誤或沒講清的大家可以留言！！
