節點和點的案例
ABAQUS 輸出節點坐標和積分點坐標
總結inp中添加關鍵字
輸出單元的積分點坐標：*EL FILE
COORD
輸出節點坐標：*NODE FILE
COORD
原貼出處：https://www.researchgate.net/post/How-to-find-integration-point-coordinates-in-Abaqus-CAE
這是帖子討論的，但是我的嘗試是兩個COORD生成的結果文件是一樣的，都是節點坐標
展開 你不知道的CAE小常識（二十一）（二十二）
你不知道的CAE小常識（二十一）（二十二）
關于有限單元法中節（結）點與積分點的幾點釋疑
節點和積分點是有限單元法（FEM）的兩個基本概念，初涉有限元計算的同志往往在這點上產生混淆，假設導師面試的時候，問單元應力是什么，若回答不慎，將貽笑大方，得不償失。本文試圖以簡略易懂的說法來闡述節點和積分點的區別。
1.節點位移是有限元法的基本未知量。節點構筑了問題域的幾何離散化形狀，節點是形函數的零點，通常形函數是以節點為依據進行假設的。形函數決定了單元內部各點運動的位移模式（常用帕斯卡三角形來選擇單元位移模式），這樣就形成了數學上所說的插值。
有限元法的原理就是將問題域分割成N多小單元，在每個單元內采用簡單的函數來近似表達單元的真實位移，將各單元再連接起來，就可以近似描述整個問題域的運動。因此，有限元法從根本上就是精確的，而不是準確的。
2.積分點是單元進行數值積分的已知量。有限元法中一般采用高斯積分，但是積分方法不限于高斯積分，如果有人用了Irons積分或者Hammer積分，請不要驚訝。在形成單元剛度矩陣和進行節點應力磨平的時候，需要高斯積分。
以等參單元為例，其剛度矩陣 ，這個就需要數值積分來快速計算，高斯點坐標及權系數如表4.2[王勖成]所示。
老師授課時一般對常應力單元進行推導，而常應力單元只有一個積分點，被積函數是常數，因此體現不出高斯積分來。很多老師對高斯積分在單元剛度矩陣的應用不予細述，導致部分同學對單元積分點認識不足。
3.單元應力指的是高斯積分點的應力，而非節點上的應力。有了位移模式，再通過虛功原理得到單元剛度矩陣，然后聚合總剛，求解平衡方程，就會把基本未知量——節點位移求出來了。通過節點位移得到單元應變結果，利用物理方程求得單元應力結果。
展開 關于有限單元法中節（結）點與積分點的幾點釋疑
節點和積分點是有限單元法（FEM）的兩個基本概念，初涉有限元計算的同志往往在這點上產生混淆，假設導師面試的時候，問單元應力是什么，若回答不慎，將貽笑大方，得不償失。本文試圖以簡略易懂的說法來闡述節點和積分點的區別。
1.節點位移是有限元法的基本未知量。節點構筑了問題域的幾何離散化形狀，節點是形函數的零點，通常形函數是以節點為依據進行假設的。形函數決定了單元內部各點運動的位移模式（常用帕斯卡三角形來選擇單元位移模式），這樣就形成了數學上所說的插值。
有限元法的原理就是將問題域分割成N多小單元，在每個單元內采用簡單的函數來近似表達單元的真實位移，將各單元再連接起來，就可以近似描述整個問題域的運動。因此，有限元法從根本上就是精確的，而不是準確的。
2.積分點是單元進行數值積分的已知量。有限元法中一般采用高斯積分，但是積分方法不限于高斯積分，如果有人用了Irons積分或者Hammer積分，請不要驚訝。在形成單元剛度矩陣和進行節點應力磨平的時候，需要高斯積分。
以等參單元為例，其剛度矩陣
，這個就需要數值積分來快速計算，高斯點坐標及權系數如表4.2[王勖成]所示。
老師授課時一般對常應力單元進行推導，而常應力單元只有一個積分點，被積函數是常數，因此體現不出高斯積分來。很多老師對高斯積分在單元剛度矩陣的應用不予細述，導致部分同學對單元積分點認識不足。
3.單元應力指的是高斯積分點的應力，而非節點上的應力。有了位移模式，再通過虛功原理得到單元剛度矩陣，然后聚合總剛，求解平衡方程，就會把基本未知量——節點位移求出來了。通過節點位移得到單元應變結果，利用物理方程求得單元應力結果。
在等參元中，單元中n+1階(n=p-m)高斯積分點上的應變或應力近似解比其它部位具有較高的精度，因此我們稱(n+1)階高斯積分點是等參元中的最佳應力點。
展開 ansys的取值函數
AREAND(N1,N2,N3)
AREAKP(K1,K2,K3)
由節點或關鍵點3點為頂點圍成的三角形的面積
ARNODE(N)
從選擇且與節點N的相關的單元分配給節點N的面積，對于二維平面，返回與節點相關的邊面積，對于軸對稱體，返回與節點相關的邊緣表面積，對于三維實體，返回與節點相關的端面積
有關節點和關鍵點的方向余弦
NORMNX(N1,N2,N3)
NORMNY(N1,N2,N3)
NORMNZ(N1,N2,N3)
NORMKX(K1,K2,K3)
NORMKY(K1,K2,K3)
NORMKZ(K1,K2,K3)
由節點或關鍵點三點確定的平面法線在X或Y或Z方向的方向余弦
有關節點與單元連接性的取值函數
ENEXTN(N,LOC)
與節點N相連接的單元。
展開 
有限元計算過程中積分點應力如何外插至節點處？【公式推導篇】
注：由于技術鄰排版風格有限，故部分內容顯示不全，感興趣的小伙伴可點擊原文進行閱覽：
有限元計算過程中積分點應力如何外插至節點處？【公式推導篇】
https://mp.weixin.qq.com/s/47byQ3b3e5UpbUp7Krs2mQ
本次分享的是：有限元計算過程中，單元積分點應力如何外推至節點？
有關積分點與節點的概念可點擊跳轉閱讀歷史推文：有限元基本概念-【節點和積分點】，現科普一下Q4單元、Q8單元、Q9單元的形函數和高斯積分方案。
Q4單元
Q8/9單元
應力外插
核心理念：坐標系的轉換。
假設是母單元的自然坐標系，是由高斯積分點控制的坐標系（術語可能不專業），假設高斯積分方案為。坐標系轉換關系：
單元內任一點的應力，由4個高斯積分點應力進行插值時，可表示為
其中，是基于高斯積分點的形函數，第一個積分點的坐標在母單元坐標系下為（-1，-1），根據上述的坐標系轉換的方式，在高斯積分點的坐標系下，第一個單元節點在高斯積分點坐標系下坐標為，將此坐標值代入第一個形函數，得，相同的道理，可推導至四個節點在4個形函數下的外插矩陣：
對于Q8、Q9單元，依然可采用高斯積分方案（減縮積分）。
展開 LMS Virtual.Lab Motion_教程59之網格移動
今天給大家介紹個網格處理小技巧，在剛柔耦合過程中，可能會出現網格和模型不重合的現象，這種情況的一種解決辦法就是通過平移和旋轉網格來實現重合，那么問題來了，在VL中怎么移動網格呢？移動網格的操作在哪個軟件里面操作更方便呢？
在VL中移動網格需要進入Structure模塊。
如果需要往某個方向平移網格，使用該功能。
如果需要往幾個方向組合平移網格，使用該功能。
如果需要旋轉網格，使用該功能。
如果需要使網格上的節點和幾何的點重合，使用該功能。
更多資料請關注百度網盤LMS_VL_Motion，Moiton交流群：324201728；Motion汽車模塊交流群：264418240；Durability交流群：83853780
展開 Moldex3D模流分析之簡易的XY繪圖來完整了解成型過程的重要性質變化
步驟2：點選結果頁簽的歷程功能，在彈出精靈中選擇欲觀察的分析組別和成型階段，然后點開數據字段的下拉選單，選擇欲觀察的數據來源。
注：除了成型特性結果之外，使用者可以選擇感測節點或探針節點；感測節點需要在分析之前設定完成，求解器會在整個分析過程紀錄節點的數據值并提供用戶觀察；探針節點可以在任何時間點透過 功能實時建立，程序會讀取多段輸出的紀錄繪制歷程曲線供用戶觀察，故探針節點的準確度會隨使用者多段輸出設定而變，若沒有多段輸出結果就無法觀察。
步驟3：點擊結果項下拉選單，選擇欲觀察的結果和節點，點選繪制歷程曲線結果。
步驟4：點擊圖表右上方叫出XY曲線設定精靈，提供用戶在不同頁簽更改圖表標題、X軸向/Y軸向、圖例及其他設定。
步驟5 : 在歷程曲線精靈中，利用回收站將結果清空，點輸出繪制圖片/數據。
注：使用儲存后自動打開文件選項來直接從輸出路徑觀看結果。
展開 SAP2000面荷載傳遞剖析
博文【膜、板、殼單元應用簡述】中提到“……在SAP2000中設置膜、板、殼單元時，封邊梁彎矩會出現不連續的情況，這是因為荷載傳遞時通過節點和剖分點傳遞的，板、殼單元由于考慮面外剛度，荷載按剛度由節點傳遞，梁端較大，膜單元由于忽略面外剛度，各點傳遞荷載相等，且無扭矩產生……”，換句話說，殼單元與膜單元即使在同樣導荷方式下，荷載傳遞也是不一樣的。
下面取3mx6m板為例就上述現象說明，板四邊封梁，自重乘數設置為0，恒載定義為100kN/m2，導荷方式選取“均布殼荷載/Uniform Shell Load”，計算結果如圖1、圖2。
（a）殼單元
（b）膜單元
圖 1 支座反力圖
圖 2a 殼單元梁內力
圖 2b 膜單元梁內力
由結果可發現，當板定義為殼單元時，板上荷載合力 (3x6x100) kN與支座反力 (450x4) kN是相等的，梁上傳來的節點力大小與剛度成正比，越靠近支座處數值越大，板剖分越細越明顯（側面說明殼單元考慮板面外剛度），越細越接近于真實力狀態，其中每個支座借助于板面外剛度傳遞過來的荷載為 (450-110.6-265.7=73.7) kN；當板定義為膜單元時，板上荷載合力 (3x6x100) kN與支座反力 (200x4) kN不平衡，每個單元向四個角點傳遞荷載，每個角點分擔每個單元的 (1x1x100/4=25) kN，梁上的節點分擔荷載相等（無扭矩產生），支座反力為 (125+75=200) kN，即由于不考慮面外剛度，支座反力僅包括梁傳來的荷載，有 (450-200=150) kN荷載沒有傳遞，保留在板內單元節點上。若不進行剖分，則支座反力為 (450) kN。
展開 ANSYS apdl中如何使齒輪饒定軸轉動
問題描述：在ANSYS中，實體單元和平板單元只有平動自由度，無旋轉自由度，網上提供的方法，在柱坐標系下通過固定徑向位移，在周向施加小位移，但這種方法只適合小位移的轉動，無法實現大位移，本文提出采用MPC184-銷軸單元和MPC-184剛性梁單元完成。本文就GUI的方式來介紹如何來創建齒輪的繞定軸旋轉。達到下圖的效果
step1 定義單元類型
a plat182 單元 模擬齒輪
b mpc184-剛性梁單元
c mpc184-銷軸單元 （本文繞Z軸旋轉，如定義為x軸，需要定義局部坐標系，繞y軸旋轉90°）
step2 定義局部坐標系默認，本文定義12號
step3 定義銷軸截面以及單元坐標系
step3 創建銷軸連接單元
在齒輪的中心點分配3號銷軸單元
step4 創建剛性梁單元
單元屬性旋旋轉2號剛性梁單元，去內徑的節點和圓中心點創建剛性梁單元
step5 施加載荷使齒輪旋轉2圈
setp6 求解設置
step7 時間后處理選擇內徑上的一個節點觀察ux，uy，rotz，可以看出齒輪旋轉2圈
展開 Abaqus的歷史數據輸出的問題
Abaqus的歷史數據輸出的問題
Abaqus結果輸出時有場輸出（Field Output）和歷史輸出(History Output)，場變量輸出用于描述某個量隨空間位置的變化，歷史變量用于描述某個量隨時間的變化，區別是場變量輸出大量的單元或節點上的計算結果，寫入odb文件的頻率低，用與生成后出的各個圖。歷史變量輸出少量單元或節點上的計算結果，寫入Odb的頻率高，用于生成X-Y圖。歷史變量允許單獨輸出某個獨立分量，經常用到的就是這個獨立分量的輸出，比如輸出結構最大點的位移變化曲線，或者盈利最大點的應力變化過程等。
在輸出歷史數據時，大多數的輸出結果比如位移和應力等，不能實現整個模型全部輸出，也就是需要指定set（參考集），需要設置相關的set后，結果輸出才會有這些數據。
定義Set的對話框如圖1所示，定義Set的類型時有三個，按照以往的理解，一般經常用到的應該是Node，也即節點集，想要輸出這個節點在求解過程中的Mises應力變化過程，選擇某個節點，通過Step模塊的History Output Requests，如圖2所示。
圖1
圖2
選擇剛才設置的Set，選擇需要輸出的結果項，此處為Mises。之后計算，計算后提取結果時，發現并沒有輸出該節點的Mises歷史數據。
筆者很快想到了，有可能是不能輸出節點的應力，因為涉及到節點解和單元積分點解的問題，這個之前已經討論過。基于此考慮，那么輸出該節點的位移場肯定沒問題，因為節點位移解是最先得到的節點解。
將Mises換成UT，如圖3所示。
同時，新建一個Element Set，因為Node Set無法輸出歷史數據，那么采用Element Set也肯定能夠輸出單元積分點的數據，如圖4所示。
展開 案例14-基于Nelson-Vaugh轉子代表模型的軸組件轉子動力學
定義曲面上的引導點
4. 在該引導點上創建節點質量
在三維實體模型中，粘接的面-面接觸對用于定義橫截面單元和引導點之間的接觸。引導點位于橫截面的中心。如上圖所示，使用CONTA174單元（3-D 8節點面-面接觸）對橫截面單元和/或接觸表面進行網格劃分。
類似地，對于二維軸對稱模型，粘接的點-面接觸對用于定義橫截面單元和引導點之間的接觸。引導點位于橫截面的中心。如下圖所示，使用CONTA174單元（3-D 8節點點-面接觸）對橫截面節點和/或接觸表面進行網格劃分。
當對兩個接觸對進行建模時，多點約束（MPC）算法用于KEYOPT（2）=2的粘接接觸。
TARGE170單元用于對兩個模型中的引導點進行網格劃分。該引導點的所有自由度（DOF）通過KEYOPT（4）=111111來約束。
以下輸入片段顯示了為二維軸對稱模型創建接觸對所涉及的步驟：
軸承建模
軸承用于橫向支撐轉子。使用COMBI214單元對兩個相同的無阻尼和線性正交各向異性軸承進行建模，如下圖所示：
KEYOPT（2）=1用于定義平行于YZ平面的平面中的單元。這些單元的自由度則是UY和UZ。下圖中應該添加X、Y和Z軸。
以下是軸承建模的一般流程：
1. 在軸承位置創建接觸對（幾何圖形在這些位置預先切片）。
2. 對三維實體模型使用CONTA174單元，對二維軸對稱模型使用CONTA175單元對接觸面進行網格劃分。
3. 在每個軸承位置（引導點位置）的轉子橫截面中心創建節點。它代表軸承的第二個節點。
4. 使用這些節點和引導點創建COMBI214單元。
必須創建下圖的接觸對來建模軸承，接觸對的建模方式與前一節圓盤的建模方式類似。
展開 
仿真筆記——關于CAE有限元分析網格，需要從整體上理解這幾點
通常說的高次單元和線性單元的區別主要體現在積分點的區別上，高次單元理論上具有更高的精度，但對塑性問題不合適。工程上常常通過增加網格數量和節點數量來提高線性單元的計算結果精度。
文章來源：CAE仿真學社
三大沖壓仿真軟件AutoForm/DynaForm/PAM-STAMP綜合對比
，不知道誰是原創：））；三個軟件都支持在多工步計算中自動修邊，但PAM-STAMP沒有單獨的修邊工序；下面為結果截圖
點擊圖片查看大圖
測試5：回彈測試
三個軟件都可以對修邊前后的工件進行回彈的模擬，結果如下：
點擊圖片查看大圖
設置過程：三個軟件的回彈設置都相對不是很復雜；
計算速度：AutoForm最快；PAM-STAMP次之，DynaForm最慢（LS971R6之后會有改善）
約束：三個軟件都支持無約束回彈，也都支持節點約束和物理點坐標約束（DynaForm5.8之后才支持）
計算結果：使用同樣的點約束設置，3個軟件的結果差異性較大，這個也是回彈計算的特點，這也是筆者不再進行回彈補償測試的原因，如果連回彈都不準確，補償就更談不上了，回彈補償的流程3個軟件都基本一致，只是PAM-STAMP需要單獨的軟件映射曲面罷了；
點擊圖片查看大圖
通過以上的測試結果，三個軟件都可以非常方便穩定的對全工序的沖壓過程進行模擬分析，從計算效率上講，AutoForm無愧于速度第一的美譽，目前確實是其他2款軟件沒法比擬的，從設置過程看，AutoForm和DynaForm的自動化程度要高一些，PAM-STAMP的設置方法已經好多年沒有變化了，設置比較繁復，且很多時候都需要手動修訂，從計算結果看，AutoForm的計算結果相對其他2款軟件，其應變偏小，結果看起來更理想，而PAM-STAMP則容易出現破裂，失效風險更大（同等設置條件），所以在日常應用時，3個軟件的設置參數需要進行積累和微調，不能以一個標準進行分析；
沖壓過程的分析也需要積累和標準化的，不同零件及工況需要不同的設置參數，材料數據的積累也很重要，這些都需要應用部門進行大量的工作；軟件只是工具，分析結果的準確度不光僅僅靠軟件，應用人員的經驗也很重要。
展開 有限元理論基礎及Abaqus內部實現方式研究系列47：約束關系（3）-船舶規范約束導致的Max Ratio問題
3.3 原理分析和猜測
3.3.1 約束關系的原理分析
即兩個主節點的UR2，UR3與船體端面的節點的位移無關。
同時，主節點只有U2，U3、UR1約束。導致主節點（即獨立點）剩下的UR2，UR3自由度沒有約束或者綁定在艙段，因為主節點這兩個自由度度沒有和任何其它節點組成單元或者約束關系，按我們前面的討論，剛度陣都是兩個節點之間的關系，所以導致剛度陣UR2和UR3相關的列或者行都是0，一個矩陣只要有一列或者一排都是0，那么理論上剛度陣K就無法求逆，所以正常來說，按船舶CSR規范的約束理論上都會導致兩個MPC主節點的UR2和UR3都約束不夠。
3.3.2 實際編程的分析
實際編程和理論相比存在差異：
3.3.2.1 模型或者精度誤差
實際編程不會直接按上面少了簡化自由度的（2）式來編寫，依然按全部自由度的（1）來編寫，這樣UR3前的rx1依然要寫到程序代碼中，由于創建模型、計算機存儲浮點數存在截斷誤差或者rx1可能也是某些向量叉乘得到的誤差，那么rx1不一定是絕對的0，有限元求解器不一定會報這個錯誤了。
譬如上面正確的1號模型，查看前端面的所有節點和獨立節點的x坐標，正常來說都應是121500這個x坐標，但發現70648號和72048號節點卻是121504，和主節點差了一個小量，相對誤差0.003%。這個小偏差導致主節點的UR2和UR3和艙段的U2、U3有約束關系，間接的綁定了這個主節點的UR2和UR3自由度。后端面也有兩個類似存在小誤差的節點。
將這四個點的X值人為修改為所在MPC同一平面后，iSolver計算也出現了2號模型的錯誤結果，此時明顯位移極大不正常。同時，Nastran也是計算失敗，報告PIVOT RATIOS過大的致命錯誤。
展開 基于Abaqus的混凝土箱梁熱力耦合分析
_iz=58558&from=article.pc_detail&lk3s=953192f4&x-expires=1745730756&x-signature=v6d3dad73HuWnY9OqhxK3uG%2BKW4%3D" style="margin: 0px auto 8px; padding: 0px; max-width: 100%; display: block; border-radius: 4px; cursor: zoom-in; height: auto;"><p class="pgc-img-caption" style="margin-top: 20px; border: 0px;"></p></div><p class="ql-align-center">圖1 箱梁日照熱邊界示意圖</p><p class="ql-align-justify">根據《Abaqus幫助文檔》第16.1.2和16.1.3節，Abaqus可以進行順序耦合傳熱分析，即先進行傳熱分析，再進行應力/變形分析。此時節點溫度作為時間的函數存儲在傳熱結果(. fi1)文件或輸出數據庫. odb)文件中，要在傳熱、溫度－位移耦合或應力／位移分析中定義載荷，可以從 abaqus emload實用程序創建的輸出數據庫. odb)文件中讀取節點集中熱通量和點載荷。</p><p class="ql-align-justify"><strong>【模型信息】</strong>箱梁截面尺寸信息如下圖所示，混凝土節段取3m，混凝土強度等級為C60。熱膨脹系數取10^-5C^-1oi，比熱容取960J/(kg*C)，導熱系數取1.62W/(m*C)，導熱系數取1.62W/(m*C)。
展開 