abaqus 輸出剛度陣的案例
對總體質(zhì)量針和剛度陣輸出問題的總結(jié)
因為前一陣遇到從NASTRAN中導(dǎo)出總體質(zhì)量陣和剛度陣的問題，在網(wǎng)上搜尋了一下，發(fā)現(xiàn)有很多種方法，我自己整理了個文件，希望對后面有此需要的人有所幫助。
注：我建的模型是22個節(jié)點，故總體質(zhì)量陣和剛度陣時132*132的
另外也想求助一下：導(dǎo)出質(zhì)量和剛度陣后大家一般用到什么地方？導(dǎo)出的矩陣因為是奇異陣不能用于求逆等運算，需要作些什么特殊處理？
同時NASTRAN中的矩陣有g(shù)-set,f-set,a-set,n-set等，有沒有高手知道在NASTRAN內(nèi)部計算特征值時用的是什么類型的矩陣呢，是直接對總體質(zhì)量和剛度陣分析，還是用的縮減后的a-set 矩陣？
輸出總體剛度陣和質(zhì)量陣.pdf
展開 【ABAQUS模態(tài)動力學(xué)】Composite&abaqus 預(yù)應(yīng)力模態(tài)分析&輸出單元剛度矩陣
劃分網(wǎng)格
定義邊界條件
5.2 預(yù)應(yīng)力模態(tài)
預(yù)應(yīng)力模態(tài)，按我的理解就是，假設(shè)t=0時，結(jié)構(gòu)的剛度矩陣、質(zhì)量矩陣為M0,K0；t=t1時，結(jié)構(gòu)（分析對象）收到外部激勵的作用，使得結(jié)構(gòu)的剛度矩陣，質(zhì)量矩陣發(fā)生改變，結(jié)構(gòu)的剛度矩陣、質(zhì)量矩陣為M1,K1；t=t2時，開始進(jìn)行特征值提取，此時求解的是t1狀態(tài)的結(jié)果。
從上面這個理解出發(fā)，ABAQUS預(yù)應(yīng)力模態(tài)只要在frequency分析步之前進(jìn)行General,Static分析步,打開NLGeom選項(分析過程中剛度矩陣會不斷變化)。
提取單元剛度矩陣：
【ABAQUS 二次開發(fā)筆記】輸出單元剛度矩陣 - hayden_william - 博客園
以上均為我的一點理解，不一定完全正確，本文僅作為個人學(xué)習(xí)記錄之用，其他概不負(fù)責(zé)。
展開 Abaqus一鍵輸出剛度矩陣和質(zhì)量矩陣的插件EMM ￥20
=======概述=======
EMM(Export Matlab Matrix)是集成在ABAQUS/CAE中的一個插件，能夠一鍵輸出Abaqus模型的單元及全局剛度、質(zhì)量、載荷矩陣，并自動轉(zhuǎn)換為MATLAB矩陣。
Abaqus可以輸出剛度、質(zhì)量等單元或者全局矩陣。但需要手動添加關(guān)鍵詞，較為麻煩，且輸出的.mtx矩陣格式比較亂，不易閱讀。
本工具能一鍵實現(xiàn)以下功能：
1. 輸出Abaqus模型的剛度、質(zhì)量、載荷矩陣到.mtx，包括單元和全局的；
2. 將上述矩陣轉(zhuǎn)換為MATLAB的.mat文件；
3. 用Matlab打開上述.mat文件。
如果用戶有任何問題或者需要合作，歡迎聯(lián)系我們。email: SnowWave02@qq.com
=======使用過程=======
打開Abaqus模型，切換到Job模塊，在Abaqus菜單欄的Plug-ins里看到iSolver插件的菜單。
點擊iSolver->Export Matlab Matrix…，彈出EMM界面。
勾選需要輸出的矩陣和Set，點擊EMM界面右側(cè)的Submit后將提交Abaqus計算，計算完畢得到如下類似矩陣。
=====演示視頻======
http://www.yqgqt.org.cn/college/video/c11343
=====下載======
用到的模型
Job-Cabin-T.rar
用戶手冊
Export Matlab Matrix工具說明(SnowWave02 20191217).pdf
EMM插件
展開 