ansys 高斯積分點的案例
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中平面應力單元高斯積分點的順序為：
展開 有限元中的高斯點與積分詳解(下)_《數值計算與程序設計》系列課程之九 ￥599
本課從實際問題出發，帶著問題去講解有限元中的高斯點與數值積分。一開始拋出了以下3個關鍵問題：
1.對于一個任意函數怎么去得到它的積分？
2.數值積分的本質是什么？為什么簡單地取幾個點就可得到積分值？此種方法的立足點在哪？
3.很多資料上都說“有限元求解精度嚴重依賴于網格質量，過度扭曲的單元會導致結果不收斂或者精度極度惡化”，這只是為什么呢？扭曲單元到底影響的是有限元方法中的哪一步？
圍繞這3個問題，本課分別講了一下三個內容：
1. 數值積分基本方法。
2. 有限元單元積分。
3. 誤差分析。
希望有興趣的同學多多支持下，你們的支持是我更新的動力
展開 有限元中的高斯點與積分詳解(上)_《數值計算與程序設計》系列課程之五 ￥599
本課從實際問題出發，帶著問題去講解有限元中的高斯點與數值積分。一開始拋出了以下3個關鍵問題：
1.對于一個任意函數怎么去得到它的積分？
2.數值積分的本質是什么？為什么簡單地取幾個點就可得到積分值？此種方法的立足點在哪？
3.很多資料上都說“有限元求解精度嚴重依賴于網格質量，過度扭曲的單元會導致結果不收斂或者精度極度惡化”，這只是為什么呢？扭曲單元到底影響的是有限元方法中的哪一步？
圍繞這3個問題，本課分別講了一下三個內容：
1. 數值積分基本方法。
2. 有限元單元積分。
3. 誤差分析。
本次課程分為上下兩課，第一課講了第一和第二個內容。關鍵詞是：數值積分的本質，有限元高斯積分（附件中包含1個小時的詳細課程視頻以及PPT）。
在第二課中，再繼續展開第三部分內容，誤差分析，解決問題“扭曲單元到底影響的是有限元方法中的哪一步”。
希望有興趣的同學多多支持下，你們的支持是我更新的動力
展開 Ansys 查看高斯點上的應力
許多時候我們需要在ANSYS中查看高斯點上的應或者和應變，然而我們看到的節點上的應力或者應變通常是由高斯點上的應力或者應變外插而來，這時候我們就需要用到ERESX這個命令了。
ERESX命令使用格式：ERESX，Key(GUI: Main>solution > Load Step Opts > Output Ctrls > Integration Pt或Main Menu > Preprocessor > Loads > Load
Step Opts > Output
Ctrls > Integration Pt)
Key為外插法控制鍵，有DEFA,YES和NO三個選項，分別對應著三種情況：
DEFA(默認設置)：除了具有塑性、蠕變或膨脹等非線性特性的單元意外，將積分點的結果進行外插擴展到所有單元的節點上。
YES: 將積分點的結果進行外插擴展到所有單元的節點上，僅將線性結果數據通過外插法擴展到這些具有塑性、蠕變或膨脹非線性特性的單元上。
NO: 將積分點上的結果復制(不是外插)到所有單元的節點上。
顯然，當我們不確定ANSYS是如何外推的，想直接查看高斯點上的應力、應變或其它結果的時候，我們就可以直接使用ERESX,no這個命令來查看了。
注意：對于非線性的數據ANSYS總是采用復制的方式擴展到節點上，而不是外推法,當 然，你也可以用ERESX,yes來采用外推法；這個命令同樣可以在prep7中使用；
轉載來源于
http://blog.sina.com.cn/s/blog_934e096a0102wkyb.html
展開 
