三角形單元
關注創建者：素釵_XFEM 創建時間：2021-05-21
三角形單元的視頻教程
（未完結，謹慎下單）ABAQUS用戶單元子程序（UEL）理論推導及程序實現
，分別基于matlab和fortran語言進行編程，讓大家更加直觀的掌握UEL的工作原理和實現方法； （五）基于理論推導的結果，以平面應力問題的三角形單元為例演示程序的實現； （六）基于UEL實現擴展有限元方法（XFEM）模擬材料的夾雜； （七）待定..........
￥80 3小時7分鐘 2513播放
查看
ABAQUS UEL自定義單元子程序手把手實例研究（Fortran語言）
； (8) 3節點三角形單元的UEL子程序; (9) ABAQUS官方手冊上的UEL子程序; (10) 網上一些經典的UEL子程序; (11) XFEM（擴展有限元）的UEL子程序(學習咨詢)。
￥200 11小時42分鐘 20516播放
查看
Matlab有限元編程從入門到精通40講：快速獲得各典型有限元案例的Matlab代碼
涉及的單元有桿單元，梁單元，三角形單元，四邊形單元，板殼，四面體/六面體實體單元等等，非線性問題涉及材料非線性、幾何非線性、接觸非線性等非線性問題！總課時預計大于45課時，并不斷更新中（1.0版本10課時，2.0版本已更新至45課時）！ 隨課附贈學習包：【課程講義】 【案例源碼】 【理論文本】 【30余本經典有限元教材】?。?/p>
￥499 20小時26分鐘 3036播放
查看
三角形單元的實例教程
1 引言
有限元分析的單元類型以及單元的質量影響著數值模擬的結果，在論壇中有人問到了這個問題，因此本文僅對二維有限元實體單元類型作簡要的討論，討論的條件包括: (1) 巖土工程問題使用的主要單元類型；(2) 不包括結構元、界面元【板樁墻(Sheet Pile Wall)模擬---FEM中的界面元】以及一些特殊類型的單元，例如不連續元【使用Plaxis的不連續單元(Discontinuity Elements)】。
2 單元類型
(1) 3節點的三角形單元T3
三節點的三角形單元(3 Noded Trangles)是有限元分析中最簡單的實體單元，在實踐中已經很少使用，典型的應用包括RS2, ADONIS【二維FEM分析軟件ADONIS(V3.50.1)命令流】，EnFEM【EnFEM---一個小型的有限元分析教學軟件】等。
(2) 6節點的三角形單元T6
六節點的三角形單元(6 Noded Trangles)是有限元分析中經常使用的實體單元，6節點的三角形對位移進行了二階插值，數值積分涉及三個高斯點。如果劃分足夠多的單元，那么T6在標準的平面應變分析中可以得出非常精確的結果，不過在軸對稱應用或在使用強度折減法(phi-c )時應小心使用，它很可能高估其破壞載荷和安全系數。因此最好使用下面的T15單元。T6典型的應用包括Plaxis, RS2, ADONIS, EnFEM等。
(3) 15節點的三角形單元T15
15節點的三角形單元(15 Noded Trangles)能夠更精確地表示應力應變關系，特別是對于軸對稱問題，目前只有Plaxis提供這種單元。T5對位移進行了四階插值，對12個高斯點(應力點)進行數值積分。T15單元能夠精確計算一些巖土工程特殊問題的應力，例如不可壓縮土層中的塌陷計算。
展開 利用上面討論的三角形常應變單元解平面問題，其具體步驟可歸納如下:
1)將要計算的彈性體劃分成三角形單元。對結點進行編號，列出結點坐標作為輸入信息。
(2)對單元進行編號，列出單元三個結點的號碼作為輸入信息。
(3)計算載荷的等效結點力，把等效結點力作為輸入信息。
(4)按照(6)式計算各單元的常數bi、ci、bj、cj、bm、cm，再按照(4)計算2A。
(5)按照(35)式計算各單元的剛度矩陣。
(6)形成整體剛度矩陣。
(7)處理約束及消除剛體位移。
(8)解線性方程組(32)式，求結點位移。
(9)按照(20)式計算應力矩陣，再按(18)式計算單元應力。根據需要計算主應力和主方向。
通常步驟(4)至(9)均由計算機來完成，而步驟(1)至(3)可以用手工完成，也可由計算機來完成。在實現以上各步驟時，為了達到一定的計算精度，節約計算機存儲量，縮短計算機運行時間等目的，還需要注意下列事項。
1、利用對稱性
在劃分單元前要研究一下，計算對象是否有對稱變形或反對稱變形存在，從而確定是否需要取整個物體，還是取部分物體作為計算模型。例如圖8a所示受純彎曲的梁，它對于x,y軸都對稱，而載荷對于y軸對稱，對于x軸反對稱?？梢姡瑧蛻円鄬⒕哂型瑯拥膶ΨQ和反對稱特性，所以我們只需計算1/4梁就行了。分離體如圖8b所示。對于刪去部分結構的影響可以這樣考慮：對于處于y軸對稱面內各結點的x方向位移和y方向分布力都應等于零，而對于處在x軸反對稱面上的各結點的x方向位移和y方向分布力亦都應等于零。這些條件相當于安置如圖8b中的約束。圖中o點上安置y方向的約束是為了消除剛體位移而設置的。又例如在分析圖9中所承受均勻壓力的厚壁圓筒時，根據結構和載荷軸對稱的性質，我們可以取出一個小扇形(圖中陰影部分)進行計算。
展開 元素（單元）的形狀原則上是任意的。二維問題一般采用三角形單元或矩形單元，三維空間可采用四面體或多面體等，每個單元的頂點稱為節點（或結點）。
步驟2：單元分析
進行分片插值，即將分割單元中任意點的未知函數用該分割單元中形狀函數及離散網格點上的函數值展開，即建立一個線性插值函數。
步驟3：求解近似變分方程
用有限個單元將連續體離散化，通過對有限個單元作分片插值求解各種力學、物理問題的一種數值方法。有限元法把連續體離散成有限個單元：桿系結構的單元是每一個桿件；連續體的單元是各種形狀（如三角形、四邊形、六面體等）的單元體。每個單元的場函數是只包含有限個待定節點參量的簡單場函數，這些單元場函數的集合就能近似代表整個連續體的場函數。根據能量方程或加權殘量方程可建立有限個待定參量的代數方程組，求解此離散方程組就得到有限元法的數值解。
有限元法已被用于求解線性和非線性問題，并建立了各種有限元模型，如協調、不協調、混合、雜交、擬協調元等。有限元法十分有效、通用性強、應用廣泛，已有許多大型或專用程序系統供工程設計使用。結合計算機輔助設計技術，有限元法也被用于計算機輔助制造中。
展開 算例中采用的三角形單元為ABAQUS單元庫中的二維三角形單元CPS3。
算例一：平面三角形單元
一、三角形單元有限元計算參數
圖3 平面三角形單元示意圖
三角形單元的相關計算參數如表2所示：
表2 三角形單元計算參數
該三角形單元共有三節點（如圖1所示），每個節點有2個自由度
二、邊界條件及載荷設置
將節點1與節點3的位移自由度施加約束，在節點2的水平x方向施加100N的節點力。
單元網格的形成實際上屬于有限元計算中的前處理部分，即確定單元節點信息，當模型較為復雜時，用戶可在Abaqus、Ansys等大型商業有限元軟件中進行建模，導出網格信息。
當模型較為簡單時，如二維平面板模型，用戶可基于一些較為基礎的網格生成算法，在自己的程序中通過控制模型長、寬等信息，即可生成有限元網格。
看似應用有限，但是在一些比較復雜的領域內，往往需要先在簡單的模型中得到理論驗證，如此以來，有利于自編程代碼的完整性，即前處理、內核計算、后處理于一體。
本篇推文，木木就帶著大家，學習一下CST、LST單元網格的自動生成。
CST單元網格
單元自動網格劃分
如下圖所示，為3節點三角形單元網格生成示意圖，圖中NXE和NYE分別是模型橫向和縱向單元個數，dhx和dhy分別是單元的橫向、縱向長度。
展開 
三角形單元的相關專題、標簽、搜索
三角形單元的最新內容
考慮熱源的瞬態熱傳導有限元求解器2小時前
關鍵詞：熱源，瞬態，熱傳導，有限元求解器，三角形單元，自研
在《瞬態熱傳導有限元求解器開發》一文中，我們介紹了自研的二維瞬態熱傳導求解器。
當時那個控制方程沒有考慮熱源，邊界條件中只涉及溫度、熱流、對流。然而在很多問題中，熱源才是最關鍵的邊界條件，比如電發熱、化學反應生熱。
熱源的處理
熱源是體熱，相對應的熱流是面熱。
基于ansys apdl建立單元截面分層的材料參數
建立的截面，多少段，多少個自定義截面
也有使用CAX6M三角形單元類型來模擬坯料的，2種單元區別對比如下：
材料
假設用于坯料的材料模型依據 Lippmann (1979) 給出。其楊氏模量為 200 GPa，泊松比為 0.3，密度為 7833 kg/m3。采用了一個率無關的 von Mises 彈塑性材料模型，其中屈服應力為 700 MPa，硬化斜率為 0.3 GPa。命名為METAL。
使用子程序法定義任意單元刪除準則，不受算法模型限制。
在仿真界，Garbage In, Garbage Out 是鐵律。無論你用 Ansys 還是 Abaqus，網格（Mesh）質量直接決定了結果的生死。今天聊聊幾個繞不開的核心指標：
1?? 雅可比比率 (Jacobian Ratio) 衡量單元從理想形狀（如正方形）映射到實際形狀的變形程度。理想值為1.0。當雅可比值為負或過小時，意味著單元發生了自交或極度扭曲，會導致剛度矩陣奇異，計算直接崩潰
瞬態熱傳導有限元求解器開發2個月前
關鍵詞：瞬態，熱傳導，有限元求解器，三角形單元
熱傳遞有三種方式：熱傳導、熱對流、熱輻射。就熱傳導問題而言，無論是結構力學還是流體力學都會涉及，兩邊都沒拿它當外人。
前面的文章提到過，結構力學的有限元發展地非常成熟，大部分的剛度矩陣在文獻里面都推導好了。而流體力學的很多單元類型的有限元方程，可能需要自行推導完成。
使用高階單元
改用三角形二次單元作為速度單元，三角形線性單元作為壓力的單元，并且不使用罰函數，速度和壓力耦合求解，結果就對了。
改進后速度、壓力結果
文獻結果
利用Sim 3D 2206版本對赫姆霍茲諧振器進行聲學仿真求解時，出現“流體邊界條件 Visco-Thermal Treatment(1)部分或完全分布于非流體和非多孔彈性單元”的錯誤提示。麻煩問一下3個月前
[圖片]
Ansys Workbench ACT插件，由窗口選中體單元，提取體積和表面積，計算幾何特征尺寸
問題：
在FKM關于結構疲勞評估計算方法中指出：零部件特征尺寸，影響疲勞結果評估。原因是材料的應力壽命曲線是由標準試樣進行試驗測試獲得的。當零部件的特征尺寸與測試樣件不一致時，需要考慮零部件的特征尺寸這一因素。（一般而言，當零部件的尺寸大于材料標準測試樣件時，零部件的表面或內部缺陷發生的概率會增加
對于實際應用中承受非線性彈簧單元Combin39的實際應用。
在ANSYS Workbench里提供了兩種方法，一種是WB的雙向彈簧，輸入數據表格，其本質上采用是LINK8單元進行模擬，而不是非線性彈簧combin39。
而利用Combin39單元，需要建立彈簧單元后，插入命令流來實現，對于只承受壓縮載荷的力-位移曲線，輸入到最后，是需要稍等小的正位移和正力數值。
