ansys剛度矩陣推導
關注創建者：王靖雯 創建時間：2023-03-07
ansys剛度矩陣推導的實例教程
材料的彈性特征可以使用剛度矩陣來進行完全描述。剛度矩陣中任意元素都可以根據由第一性原理計算得出的聲子色散曲線來確定。
在本次案例研究中，我們介紹了一個用立方體結構來分析硅、金剛石和銅單晶的剛度矩陣的方法。
如圖1，硅是具有金剛石結構的立方晶體，其剛度矩陣如下所示。
由于立方晶體具有對稱性，使用SIESTA模擬軟件計算可得該模型中有3個獨立元素：C11C11，C12C12，和C44C44。
材料的彈性特征可以使用剛度矩陣來進行完全描述。剛度矩陣中任意元素都可以根據由第一性原理計算得出的聲子色散曲線來確定。
在本次案例研究中，我們介紹了一個用立方體結構來分析硅、金剛石和銅單晶的剛度矩陣的方法。
如圖1，硅是具有金剛石結構的立方晶體，其剛度矩陣如下所示。
由于立方晶體具有對稱性，使用SIESTA模擬軟件計算可得該模型中有3個獨立元素：C11C11，C12C12，和C44C44。
本圖片是取自《有限元方法基礎教程》（國際單位制）（第五版）， 第六章 建立平面應力和平面應變剛度方程。
==》本博客是對于四節點四邊形單元的剛度矩陣的推導，沒有編程的實現。
根據取得的聲速，可以按以下公式計算剛度矩陣。
聲速與(1,0,0)/(1,1,1)方向上剛度矩陣的關系。
各個箭頭表示偏振方向，紅色箭頭表示縱波，綠色箭頭表示橫波。ρ表示硅的重量密度。
表1給出分析所得剛度矩陣中每個元素值和文獻中的值之對比。（C. Kittel 《固體物理學導論》，紅色括號）。
表1 剛度矩陣所得值與文獻中的值對比
材料的彈性特征可以使用剛度矩陣來進行完全描述。剛度矩陣中任意元素都可以根據由第一性原理計算得出的聲子色散曲線來確定。
在本次案例研究中，我們介紹了一個用立方體結構來分析硅、金剛石和銅單晶的剛度矩陣的方法。
如圖1，硅是具有金剛石結構的立方晶體，其剛度矩陣如下所示。
由于立方晶體具有對稱性，使用SIESTA模擬軟件計算可得該模型中有3個獨立元素：C11C11，C12C12，和C44C44。
圖1 硅的金剛石結構
立方晶體的剛度矩陣
硅的聲子色散曲線可以用SEISTA模擬軟件直接計算（如圖2）。
根據硅的晶胞中原子數為2這一事實依據，硅有三種聲振模式和三種光學模式。硅是具有金剛石結構的立方晶體，其色散曲線表現出不同晶體取向的色散特性。
這里我們重點關注從ΓX點指向的(1,0,0)方向，和從ΓL點指向的(1,1,1)方向擴展的聲振模式色散特性。聲振模式的色散特性在近場的長波區呈線性Γ，且該梯度曲線給出了聲速。
圖2 硅的聲子色散曲線
圖3 區域1和2的色散曲線
圖中虛線表示靠近該區域的每個傳播方向上的聲速Γ點，下標L/T代表縱波和橫波。
在（1,0,0）方向傳播的聲波由一個縱波和兩個簡并后的橫波組成，它們各自的聲速和剛度矩陣可以表示為下列關系式，由此結果可得C11C11 和 C44C44。
進一步使用該結果，根據（1,1,1）方向傳播的聲波可以得到c12。在本案例中，梯度（聲速）在Γ點可以通過使用圖中1和2區域的三階最小二乘法構造一條近似曲線獲得，這是SIESTA模塊的聲子分析功能。根據取得的聲速，可以按以下公式計算剛度矩陣。
聲速與(1,0,0)/(1,1,1)方向上剛度矩陣的關系。
展開 
ansys剛度矩陣推導的相關專題、標簽、搜索
ansys剛度矩陣推導的最新內容
在有限元分析中，ANSYS 可以導出大規模稀疏矩陣（如剛度矩陣、質量矩陣），通常使用 Harwell-Boeing (HB) CCS 格式。這些矩陣對后續二次開發、動力學分析或自定義求解器非常重要，但由于其稀疏和壓縮存儲形式，直接在 MATLAB 中讀取和使用并不方便。
本文提供了 兩個 MATLAB 函數，可直接從 ANSYS 導出的 HB 矩陣文件中讀取并重構成 MATLAB 稀疏矩陣：
材料的彈性特征可以使用剛度矩陣來進行完全描述。剛度矩陣中任意元素都可以根據由第一性原理計算得出的聲子色散曲線來確定。
在本次案例研究中，我們介紹了一個用立方體結構來分析硅、金剛石和銅單晶的剛度矩陣的方法。
如圖1，硅是具有金剛石結構的立方晶體，其剛度矩陣如下所示。
由于立方晶體具有對稱性，使用SIESTA模擬軟件計算可得該模型中有3個獨立元素：C11C11，C12C12，和C44C44
材料的彈性特征可以使用剛度矩陣來進行完全描述。剛度矩陣中任意元素都可以根據由第一性原理計算得出的聲子色散曲線來確定。
在本次案例研究中，我們介紹了一個用立方體結構來分析硅、金剛石和銅單晶的剛度矩陣的方法。
如圖1，硅是具有金剛石結構的立方晶體，其剛度矩陣如下所示。
由于立方晶體具有對稱性，使用SIESTA模擬軟件計算可得該模型中有3個獨立元素：C11C11，C12C12，和C44C44
1.引論
經常使用Ansys、Abaqus等一系列有限元分析軟件進行計算、學習的學生或工程師們都會知道在有限元分析建模與計算中剛度矩陣與質量矩陣的重要性。但是由于軟件的黑盒性質，大家往往在實際使用十分成熟的商業化軟件的過程中慢慢忽視了有限元及其衍生出的商業軟件背后的原理與方法。
這時，不管是在學習中還是在工程應用中往往都會遇到一個同樣的問題，那么就是如何將Ansys
材料的彈性特征可以使用剛度矩陣來進行完全描述。剛度矩陣中任意元素都可以根據由第一性原理計算得出的聲子色散曲線來確定。
在本次案例研究中，我們介紹了一個用立方體結構來分析硅、金剛石和銅單晶的剛度矩陣的方法。
如圖1，硅是具有金剛石結構的立方晶體，其剛度矩陣如下所示。
由于立方晶體具有對稱性，使用SIESTA模擬軟件計算可得該模型中有3個獨立元素：C11C11，C12C12，和C44C44
最近在考慮自己編寫的程序和商用軟件的驗證問題，有限元結構分析中最關鍵的一環就是剛度矩陣的獲得，如果涉及到模態分析，還有質量矩陣。考慮到商業軟件的成熟性，可以用ANSYS生成的剛度矩陣做參照來看自己編寫的程序是否正確，因此如何提取ANSYS中結構的剛度矩陣，并進行隨后的驗證或者二次開發是一個問題。
https://www.yqgqt.org.cn/content/post/1796144
材料的彈性特征可以使用剛度矩陣來進行完全描述。剛度矩陣中任意元素都可以根據由第一性原理計算得出的聲子色散曲線來確定。
在本次案例研究中，我們介紹了一個用立方體結構來分析硅、金剛石和銅單晶的剛度矩陣的方法。
如圖1，硅是具有金剛石結構的立方晶體，其剛度矩陣如下所示。
由于立方晶體具有對稱性，使用SIESTA模擬軟件計算可得該模型中有3個獨立元素：C11C11，C12C12，和C44C44
1. 背景
從事結構振動控制、車橋耦合振動、結構健康監測傳感器優化布置、結構動力性能分析等等一系列研究的同仁們應該都面臨過一個同樣的問題—“怎么把結構的剛度和質量矩陣建立出來？”。這對于那些數值分析高手和專家可能不是什么問題；但是對于科研剛入門的新手來說，這個難度還是相當大的。如果都靠自己寫程序來建立有限元模型，則對理論基礎、編程水平都有很高的要求，甚至程序做出來也未必能保證其正確性，是一個很讓人頭疼的問題
本圖片是取自《有限元方法基礎教程》（國際單位制）（第五版）， 第六章 建立平面應力和平面應變剛度方程。
==》本博客是對于四節點四邊形單元的剛度矩陣的推導，沒有編程的實現。
就ansys如何提取剛度矩陣、如何解讀提取的文檔以及利用Python進行解析。
在workbench中實現整個過程的參數化過程除了前幾次文章介紹的模型與網格，還應該包括材料參數的參數化定義。利用Python進行二次開發能夠實現材料參數的自由定義，比如來源于excel表格或者文檔的數據，通過Python代碼的自動讀取，參與到實際的有限元分析進程中。
結構有限元最后的求解過程總是歸結到求解一個大型矩陣方程
