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