剛度矩陣
關注創建者：小白菜瓜 創建時間：2018-05-22
剛度矩陣的視頻教程
Matlab自編C3D8、C3D20、C3D10單元計算動荷載問題
概述：matlab自編程序求解C3D8（一階六面體）、C3D20（二階六面體）、C3D10（二階四面體）單元的剛度矩陣、質量矩陣和阻尼矩陣，并采用newmark求解動荷載問題，計算結果與abaqus保持一致。 視頻首先采用abaqus計算懸臂梁受動荷載算例，與此同時修改inp文件，并導出abaqus計算的總體剛剛度矩陣、質量矩陣和阻尼矩陣。
免費 28分鐘 157播放
查看
ABAQUS復模態分析
如果剛度矩陣和阻尼矩陣是非對稱的，則不能用常規方法將方程解耦，這時必須用復模態解耦，稱這種方法為復模態理論。對于剛度矩陣或阻尼矩陣不對稱性，復特征值分析可用于識別系統的不穩定性。例如：摩擦導致不對稱的剛度矩陣，在摩擦誘導下振動的剎車尖叫。
￥50 47分鐘 370播放
查看
基于ANSYS的實體單元扭矩施加方法總結
由于ANSYS中不能直接對實體單元施加力矩，傳統方法采用若干對力偶來代替扭矩，該方法容易導致局部應力集中；改進的方法引入一些特殊單元如rbe3單元、mpc184單元、mass21單元等，通過引入這些特殊單元，能夠比較好的實現扭矩的施加，但是特殊單元的引入又改變了整體剛度矩陣。為了解決由于引入特殊單元而導致影響整體剛度矩陣的問題，有學者等提出采用接觸單元能夠很好的解決扭矩的施加問題。
￥10 23分鐘 384播放
查看
剛度矩陣的實例教程
其可以直接導出full文件中的矩陣數據，而通過我下面給大家提供的源代碼便可以直接將導出的質量矩陣與剛度矩陣直接轉化為matlab中可使用變量，實現我們的矩陣提取操作。
Hbmat法：
! 提取剛度矩陣
/AUX2
FILE,'file',full ! 將’file’改為自己路徑下的.full文件名
HBMAT, 'Stiffness_mat', dat, , ASCII, STIFF, YES, YES !剛度矩陣
HBMAT, 'Mass_mat', dat, ,ASCII, MASS, YES, YES ! 質量矩陣
FINISH
成功導出后你會在你的工作路徑中看到儲存在Ansys中的剛度矩陣與質量矩陣。
2.2GUI方法
目前大部分此類教程都忽略了最基本的GUI方法，可能是自動帶入了經驗豐富的工程師角色，但對于初學者而言，GUI方法十分有助于理解與學習有限元軟件，因此我也將介紹提取剛度質量矩陣的GUI方法。
從而我們便能導出所需要的剛度矩陣與質量矩陣源文件了。
3.源代碼的使用與轉換（使用方法）
我們打開所導出剛度矩陣會發現里面的數字比較混亂，一時間分不出哪些是我們需要的數據，也沒有辦法直接拿它用來計算，這時，我們便需要使用我們自己的矩陣轉換m文件，其可以自動幫我們提取出剛度矩陣與質量矩陣，并生成matlab中的矩陣變量形式，我們可以直接運用這段源代碼來進行操作。
展開 在有限元分析中，ANSYS 可以導出大規模稀疏矩陣（如剛度矩陣、質量矩陣），通常使用 Harwell-Boeing (HB) CCS 格式。這些矩陣對后續二次開發、動力學分析或自定義求解器非常重要，但由于其稀疏和壓縮存儲形式，直接在 MATLAB 中讀取和使用并不方便。
本文提供了 兩個 MATLAB 函數，可直接從 ANSYS 導出的 HB 矩陣文件中讀取并重構成 MATLAB 稀疏矩陣：
1.剛度矩陣提取函數
輸入：ANSYS 導出的剛度矩陣 HB 文件（stiff.txt）
輸出：MATLAB 稀疏矩陣 K，可直接用于動力學計算或驗證
支持自動對稱化，保證數值正確
2.質量矩陣提取函數
輸入：ANSYS 導出的質量矩陣 HB 文件（mass.txt）
輸出：MATLAB 稀疏矩陣 M
使用與剛度矩陣同樣的解析邏輯，無需額外修改
案例說明：
本文以高速鐵路接觸網結構為例，展示了如何將 ANSYS 中導出的稀疏剛度矩陣和質量矩陣，在 MATLAB 中完整展開，并進行后續動力學分析準備。
通過該方法，可以將大規模有限元矩陣快速轉化為 MATLAB 可操作形式，為自定義振動分析、模態分析及其他科研或工程應用提供基礎。
優勢與應用：
支持大規模稀疏矩陣解析
自動對稱化，保證數值精度
適用于剛度矩陣、質量矩陣、其他 HB 格式矩陣
可作為動力學求解器或后處理工具的基礎模塊
使用方法：
1.使用以下代碼對ansys中生成的質量及剛度矩陣進行提取，file,5,full（5為工作目錄下full文件的文件名，例如：filename.full）。
展開 1 概述
MSC Nastran模型的剛度矩陣和質量矩陣，可以輸出為文本文件。工程實際中，工程師可以校核、集成矩陣，進行第二次開發，完成商用軟件和自研程序的完美集成。例如：工程師有一個計算線性動力學方程組的瞬態python程序，可以集成MSC Nastran的剛度矩陣和質量矩陣。
2 剛度矩陣和質量矩陣的輸出方法
1) 剛度矩陣和質量矩陣輸出至punch（.pch）文件
如果需要在其他MSC Nastran計算中，重用MSC Nastran模型的矩陣，可以將MSC Nastran矩陣輸出至Punch文件，方法為：
l 在MSC Nastran卡片中，添加參數：PARAM,EXTOUT,DMIGPCH
注：Punch文件中的矩陣，Patran不支持
2) 剛度矩陣和質量矩陣輸出至f06（.f06）文件
如果想直接在f06中查看輸出的矩陣，可以使用如下方法：
l 在執行控制部分（CEND前），添加如下卡片：
COMPILE EXTOUT $
ALTER 'RETURN'(,-1) $
MATPRN KAA,,,,// $
MATPRN MAA,,,,// $
l 添加如下參數（BEGIN BULK），例如：PARAM,EXTOUT,DMIGPCH
實例：
輸入文件：
剛度矩陣
質量矩陣
3 參考信息
適用版本：MSC Nastran 2005及以后版本。
展開 知道含義以后即可借助MATLAB或者其它軟件，讀取pch中的剛度矩陣，并編寫代碼完成剛度矩陣的輸出。圖2是小翼做的一個結構的總剛度矩陣的局部，經過與自己編寫的結構剛度矩陣輸出代碼計算的剛度矩陣對比，發現一致性較好，部分剛度系數誤差在3~4%左右，當然還是以NASTRAN的精度為準。
圖2 某結構總剛度矩陣局部
由（27a）式定義的單元剛度矩陣，由于應變B對于3結點三角形單元是常量陣，因此有
代入彈性矩陣D和應變矩陣B后，它的任一分塊矩陣可表示成
由(34)式立即可以得到
由此可見單元剛度矩陣是對稱矩陣。
為了進一步理解單元剛度矩陣的物理意義，我們同樣可以利用最小位能原理建立一
個單元的平衡，這就得到
Pe是單元結點載荷，當然應當包括其它相鄰單元對該單元的作用力。現把ae、Pe順序表示為
這是單元結點平衡方程，每個結點在x和y方向上各有一個平衡方程，3個結點共有六個平衡方程。方程左端是通過單元結點位移表示的單元結點內力，方程右端是單元結點外載。
令a1=1(ui=1),a2=a3=…＝a6=0
由（38）式可以得到
（39）式表明，單元剛度矩陣第一列元素的物理意義是：a1＝1。其他結點位移都為零時，需要在單元各結點位移方向上施加結點力的大小。當然，單元在這些結點力作用下處于平衡，因此在x和y 方向上結點力之和為零：
對于單元剛度矩陣中其他列的元素也可用同樣的方法得到它們的物理解釋。因此單元剛度矩陣中任一元素Kij物理意義為：當單元的第j個結點位移為單位位移而其他結點位移為零時，需在單元第j個結點位移方向上施加的結點力的大小。單元剛性大，則使結點產生單位位移所需施加的結點力就大。因此單元剛度矩陣中的每個元素反映了單元剛性的大小，稱為剛度系數。
展開 
剛度矩陣的相關專題、標簽、搜索
剛度矩陣的最新內容
如果輸入的應變率曲線出現交叉（即高應變率下的應力低于低應變率下的應力），或者硬化曲線呈現負斜率（未激活損傷模塊時），求解器的材料剛度矩陣將出現非正定，導致不可控的網格畸變。此外，必須通過外推確保表格覆蓋到極高應變率（如10000 /s），以防求解器在局部高變形區發生錯誤的常數外推。
CPU
雙路 Intel Xeon Platinum 8592+（64核×2）或 雙路 AMD EPYC 9655（96核×2）
千萬級自由度細網格 GCI、高維 PCE 展開、數字孿生實時驗證
內存
512GB–1TB DDR5-4800 ECC RDIMM
超大規模剛度矩陣直接內存求解
網格尺寸同時影響整體剛度矩陣的數值特性：粗網格可能導致結構剛度偏軟，而在接觸、屈曲與動態分析中，網格不足會降低求解穩定性與收斂性。因此，合理的網格尺度選擇是控制離散誤差與確保數值穩健性的核心步驟。基于模型特征尺寸與多輪劃分測試，本研究采用最大網格尺寸 18?mm、接觸面 6?mm，最終獲得 844?549 個節點與 723?723 個單元。由此可見，對稱建模顯著降低了網格規模與計算成本。
將多個物理場的自由度放在同一個大型剛度矩陣中，在一個求解器里同步迭代求解。適用于物理場之間相互作用強、必須實時反饋的場景（如壓電效應）。精度極高，但極度消耗計算資源。
流固耦合 FSI (Fluid-Structure Interaction) 工程中最常見的一類耦合。流體的流動產生壓力使固體發生變形，而固體的變形又反過來改變了流體的流場（如風機葉片形變、橋梁風振）。
Nonlinear)
線性分析假設位移與載荷成正比，剛度矩陣
$$$$ 固定不變，計算一次即可。而非線性分析中，剛度矩陣隨計算過程變化，需要通過牛頓-拉夫遜法等算法進行多次迭代，計算量呈幾何倍數增長。
當雅可比值為負或過小時，意味著單元發生了自交或極度扭曲，會導致剛度矩陣奇異，計算直接崩潰（Singular Matrix）。
2?? 長寬比 (Aspect Ratio) 單元最長邊與最短邊的比值。在應力梯度大的區域，長寬比過大會人為增加局部剛度，導致應力集中結果偏低。通常建議關鍵區域保持在3:1以內。
瞬態熱傳導有限元求解器開發2個月前
前面的文章提到過，結構力學的有限元發展地非常成熟，大部分的剛度矩陣在文獻里面都推導好了。而流體力學的很多單元類型的有限元方程，可能需要自行推導完成。在熱傳導問題中，我采用加權余量法進行處理，推導出了符合結構力學有限元文獻中給出的剛度矩陣，殊途同歸。
實際上，傳統的結構力學有限元三大控制方程：幾何方程、物理方程、平衡方程。
有限元思路
搞結構力學有限元和其他方向有限元最大的區別是：結構力學有限元發展的太成熟了，桿梁板殼，各種模型的剛度矩陣前輩都給你推導好了。我在開發結構力學有限元求解器的時候，都是先去查資料，直接就把單元剛度矩陣拿過來用。
五、理論框架的完整圖景
5.1 兩尺度耦合機制
微觀尺度（RVE內部） 宏觀尺度（結構響應）
↓ ↓
應變/損傷非均勻分布 ←——h——→ 均勻化能量密度
↓ ↓
高階梯度項（?ε, ?2ε） 退化剛度矩陣
論文采用了增廣拉格朗日方法（ALM）的巧妙策略：
引入輔助變量：將應變 作為獨立變量，不再通過位移求導得到約束條件：通過拉格朗日乘子 λ 和懲罰項 r ，強制滿足 降階優勢：最高階導數降為二階，普通二次單元即可求解
最終形成15×15的分塊剛度矩陣（每個節點3個位移+6個應變+6個拉格朗日乘子），雖然自由度增加，但避免了復雜的 C1 單元和迭代處理高階邊界條件的麻煩。
