ansys中剛度因子的案例
如何從Ansys APDL中提取剛度矩陣與質量矩陣？ ￥69
1.引論
經常使用Ansys、Abaqus等一系列有限元分析軟件進行計算、學習的學生或工程師們都會知道在有限元分析建模與計算中剛度矩陣與質量矩陣的重要性。但是由于軟件的黑盒性質，大家往往在實際使用十分成熟的商業化軟件的過程中慢慢忽視了有限元及其衍生出的商業軟件背后的原理與方法。
這時，不管是在學習中還是在工程應用中往往都會遇到一個同樣的問題，那么就是如何將Ansys APDL運行中的產生的各種數據（例如：剛度矩陣、質量矩陣）導出成為我們熟悉的形式或文件格式，從而為我們所用，所分析。
因此我決定寫下此篇文章來幫助很多實際工作或學習中需要用到此類技能的同學、同事們，讓大家更了解Ansys APDL背后的工作原理與數據導出方式。
當然，在社區中早就有大佬回答過了這個問題，并給大家制作了相應的提取矩陣軟件，其軟件具備了簡單、便捷的操作方式，讓很多想要提取剛度矩陣與質量矩陣的同僚們受益，那么我為什么還要寫一篇這樣的文章重新提起這樣一個話題呢？這就又回到了我開頭所說的“原理與方法”，我在此更希望面對想要進一步學習了解軟件背后機理的群體，并在此基礎上保留教學的簡潔性，提供導出矩陣與轉換、列式、求解的源代碼，使其既兼顧基本原理，又可以讓大家直接上手使用，非常的便捷，也避免了很多因為優化不完全導致的運行bug。
2.有限元軟件導出剛度矩陣與質量矩陣的方法
在使用APDL進行求解時，每次在求解完成后都會在工作路徑下生成一個.full文件，而這個文件十分關鍵，其正是剛度矩陣與質量矩陣的所在之處。
展開 在ANSYS中計算裂縫應力強度因子的技巧
在ANSYS中計算裂縫應力強度因子的技巧
裂縫應力強度因子用ANSYS中怎么求呀。另外，建模時，裂紋應該怎么處理呀，難道只有畫出一條線嗎？
首先說一下裂紋怎么畫，其實裂紋很簡單啊。只要畫出裂紋的上下表面（線）就可以了，即使是兩個面（線）重合也一定要是兩個面（線）；如果考慮道對稱模型就更好辦了，裂紋尖點左面用一個面（線），右邊用另外一個面（線），加上對稱邊界約束。
再說一下裂尖點附近網格的劃分。ansys提供了一個kscon的命令，主要是使得crack
tip的第一層單元變成奇異單元，用來模擬斷裂奇異性(singularity)。當然這個步驟不是必須的，有的人說起用ansys算強度因子的時候就一定要用奇異單元，其實是誤區（原因下面解釋）
好了，回到強度因子的計算。其實只要學過一些斷裂力學都知道，K的求法很多。就拿Mode
I的KI來說吧，Ansys自己提供了一個辦法（displacement extrapolation）
，中文可能翻譯作“位移外推”法，其實就是根據解析解的位移公式來對計算數據進行fitting的。分3步走，如果你已經算完了：
第一步，先定義一個crack-tip的局部坐標系，這是ansys幫助文件中說的，其實如果你的裂紋尖端就是整體坐標原點的話，而且你的x-axis就順著裂紋，就沒有什么必要了。
第二步，定義一個始于crack-tip的path,什么什么？path怎么定義？？看看幫助吧，在索引里面查找fracture
mechanics，找到怎么計算斷裂強度因子。（my god,我這3步全是在copy幫助中的東東啊）。
第三步，Nodal
Calcs>Stress Int Factr ，別忘了，這是在后處理postproc中啊。
展開 提取ANSYS中的剛度矩陣和質量矩陣到MATLAB中進行二次開發 ￥88
最近在考慮自己編寫的程序和商用軟件的驗證問題，有限元結構分析中最關鍵的一環就是剛度矩陣的獲得，如果涉及到模態分析，還有質量矩陣。考慮到商業軟件的成熟性，可以用ANSYS生成的剛度矩陣做參照來看自己編寫的程序是否正確，因此如何提取ANSYS中結構的剛度矩陣，并進行隨后的驗證或者二次開發是一個問題。
https://www.yqgqt.org.cn/content/post/1796144
受上述帖子的啟發，使用MATLAB提取ANSYS中的剛度和質量矩陣，并進行模態分析驗證提取的矩陣的正確性。
首先，在ANSYS中使用HBMAT命令方法提取整體矩陣。
命令：HBMAT,fname,ext,--,form,matrx,rhs
其中：
Fname---輸出矩陣的路徑和文件名，缺省為當前工作路徑和當前工作文件名。
ext---輸出矩陣文件的擴展名，缺省為.matrix。
form---定義輸出矩陣文件的格式，其值可取：
=ASCII：ASCII碼格式；
=BIN：二進制格式。
matrix---定義輸出矩陣的類型，其值可取：
=STIFF：輸出剛度矩陣。可用于寫入了.FULL文件的任何類型的分析。
=MASS：輸出質量矩陣。可用于特征值屈曲、子結構分析、模態分析。
=DAMP：輸出阻尼矩陣。僅用于有阻尼的模態分析。
rhs---右邊項輸出控制（右邊項指用矩陣所表示方程的等號右端矢量，這里可為節點荷載向量），如rhs=YES則輸出，如rhs=NO則不輸出。
模態分析時，因僅LANB和QR法可生成完整的質量矩陣，因此也僅采用這兩種方法時才可使用HBMAT命令得到質量矩陣文件。
展開 ANSYS中整體、單元剛度和質量矩陣的提取
將輸出信息送到output windows中
! 這時用編輯器打開cp.out文件，可以看到按單元寫出的質量、剛度等矩
陣
ANSYS中整體、單元剛度和質量矩陣的提取.rar
『分享』ANSYS中整體、單元剛度和質量矩陣的提取
將輸出信息送到output windows中
! 這時用編輯器打開cp.out文件，可以看到按單元寫出的質量、剛度等矩陣
