ansys單元剛度的案例
ANSYS中整體、單元剛度和質量矩陣的提取
指定輸出單元矩陣
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 將輸出信息送到output windows中
! 這時用編輯器打開cp.out文件，可以看到按單元寫出的質量、剛度等矩
陣
ANSYS中整體、單元剛度和質量矩陣的提取.rar
平面四邊形四節點單元計算程序與ANSYS結果對比
為什么要導出單元剛度矩陣
在學習有限元方法時，我們會需要編寫程序計算結構的單元剛度矩陣。此外，當我們需要做有限元軟件二次開發時，我們也需要驗證所做的開發是否正確。為了驗證程序正確性，我們可以從商業有限元軟件中導出單元剛度矩陣來驗證程序的計算結果。下面簡單介紹從ansys軟件中導出平面四邊形四節點單元的單元剛度矩陣。
平面四邊形四節點單元示例
如圖所示，計算這兩個單元組成單元剛度矩陣，并組裝成整體剛度矩陣，求解各個節點的位移。
展開 基于ABAQUS的UEL子程序定義4節點平面應變等參單元的剛度問題
摘要：
采用基于ABAQUS的UEL子程序開發4節點平面應變等參單元，采用雙線性形函數，4點高斯積分，本構關系為線彈性各向同性材料，得到的單元剛度矩陣和ABABUS自帶的CPE4單元的單元剛度矩陣（剛度矩陣輸出方式為*element matrix output, elset= ALLE, stiffness=yes, OUTPUT FILE=USER DEFINED）不同；對比ANSYS的單元剛度矩陣，結果顯示兩者也不相同。問題出在哪里呢？本文檔將對此問題進行回答。
本文可以作為ABAQUS高級子程序UEL的入門級教程，做UEL的應該關注下！
基于ABAQUS的UEL子程序定義4節點平面應變等參單元的剛度問題（技術鄰 藍牙）.pdf
展開 『分享』ANSYS中整體、單元剛度和質量矩陣的提取
指定輸出單元矩陣
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 將輸出信息送到output windows中
! 這時用編輯器打開cp.out文件，可以看到按單元寫出的質量、剛度等矩陣
一個單元也能干大事之單元剛度初探
本公眾號將推出“一個單元也能干大事”系列文章，試圖從最小的單元出發，來為初學者講解一些有限元基本理論。該系列文章力爭篇幅短小，簡單通俗，深入本質，一看就懂，并且希望能跟讀者討論互動。
一個單元能干什么？今天我們來研究一下四邊形等參單元，試圖帶大家探討一下它的單元剛度矩陣。來，先跟我一起玩個小游戲，放松一下。
1個單元，4個點，8個位移，4條邊，哇
2個單元，8個點，16個位移，8條變，哇哇
3個單元，12個點，24個位移，12條變，哇哇哇……
關于四邊形等參單元剛度矩陣的推導，一般有限元教材都有詳細的介紹，這里就不再討論。書看千遍，不如公式推一遍，公式推千遍，不如代碼擼一行，我們直接上代碼，進行計算。
單元剛度矩陣如下，這是一個8*8的矩陣，我們來研究一下，該單元剛度矩陣的特性。
一、單元剛度矩陣是對稱矩陣，即滿足Kij=Kji
二、單元剛陣主對角線元素恒為正值；因為主對角元素表示力的方向和位移方向一致，故總為正值。
三、我們計算一下剛度矩陣的行列式
可以看到，行列式|K|=0，即單元剛陣是奇異陣，從物理意義上來解釋，這是因為計算單元剛陣時沒有對單元的節點加以約束，雖然，單元處于平衡狀態，但容許單元產生剛體位移，故從單元剛度平衡方程不可能得到唯一位移解。
四、我們來計算一下剛度矩陣的特征值
可以看到，矩陣共有八個特征值，其中有三個零特征值。你覺得這是偶然嗎？不，冥冥之中自有天數。剛度陣的秩為5，說明剛度矩陣只有5行是線性無關的，需要約束其中的3個自由度，方程Ku=F才能求解。而約束3個自由度，就是為了消除3個剛體位移。
展開 單元剛度矩陣
由（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個結點位移方向上施加的結點力的大小。單元剛性大，則使結點產生單位位移所需施加的結點力就大。因此單元剛度矩陣中的每個元素反映了單元剛性的大小，稱為剛度系數。
展開 cohesive單元剛度選擇與參數解釋 ￥10
cohesive單元剛度選擇與參數解釋
平面四邊形單元的剛度矩陣的推導 ￥2
本圖片是取自《有限元方法基礎教程》（國際單位制）（第五版）， 第六章 建立平面應力和平面應變剛度方程。
==》本博客是對于四節點四邊形單元的剛度矩陣的推導，沒有編程的實現。
如何從Ansys APDL中提取剛度矩陣與質量矩陣？ ￥69
2.1命令流方法
（1）用戶程序法
（2）超單元法
（3）Hbmat法
這三種在網絡上都能很輕易的查到其使用的命令流與方法，但前兩種并不推薦大家使用，不過感興趣的朋友可以自行搜索盡可能多的了解一些。在這里我主要給大家介紹的方法是Hbmat法，也是我推薦大家使用的方法。其可以直接導出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方法。
從而我們便能導出所需要的剛度矩陣與質量矩陣源文件了。
展開 【ABAQUS模態動力學】Composite&abaqus 預應力模態分析&輸出單元剛度矩陣
劃分網格
定義邊界條件
5.2 預應力模態
預應力模態，按我的理解就是，假設t=0時，結構的剛度矩陣、質量矩陣為M0,K0；t=t1時，結構（分析對象）收到外部激勵的作用，使得結構的剛度矩陣，質量矩陣發生改變，結構的剛度矩陣、質量矩陣為M1,K1；t=t2時，開始進行特征值提取，此時求解的是t1狀態的結果。
從上面這個理解出發，ABAQUS預應力模態只要在frequency分析步之前進行General,Static分析步,打開NLGeom選項(分析過程中剛度矩陣會不斷變化)。
提取單元剛度矩陣：
【ABAQUS 二次開發筆記】輸出單元剛度矩陣 - hayden_william - 博客園
以上均為我的一點理解，不一定完全正確，本文僅作為個人學習記錄之用，其他概不負責。
展開 基于hypermesh和lsdyna的彈簧離散單元的建立及剛度K的計算驗證 ￥10
本貼為大家講解一下彈簧單元的建立和剛度K的一個計算驗證
模型如下：上下兩個鋼板，頂端rigid抓取的節點施加力，中間兩個節點創建spring，底面約束。
?abaqus提取單元的剛度矩陣和質量矩陣
abaqus提取單元的剛度矩陣和質量矩陣：
*Output, history, variable=PRESELECT,
*File Format,ASCII
*Element Matrix Output,Elset=Beam-1.Set-2,
File Name=shuchu,Frequency=1,Output File=User Defined,Stiffness=Yes
*End Step
*Step
*Matrix generate，stiffness
*End Step
ANSYS的生死單元模擬焊接過程
ANSYS的生死單元模擬焊接過程
1 概述
焊接模擬計算在CAE仿真是比較大的一塊內容，也是比較復雜的一個過程，幾個比較關鍵的問題是熱源函數的描述、單元的融覆、熱源的移動等等，通過單純的GUI操作，無論使ANSYS還是Abaqus都不大可能完成這個過程，通常需要借助軟件的內置語言。
本次主要介紹單元生死的應用，單元生死主要用于單元缺失的場合，比如凝固溶解過程，斷裂過程，焊接過程等等，這些過程都是非線性或者時間歷程過程，計算需要很多子步和迭代，為了在此過程中避免一遍一遍修改單元，便引入生死單元的概念，通俗的講就是通過一些方法讓單元失效，具體的改變是單元的彈性模量的改變，當單元死時，修改其彈性模量為非常小的值，讓其在求解過程中不起作用。
詳細地說，激活單元死這個狀態時，ANSYS程序將單元剛度矩陣乘以很小的因子，程序默認值為1E-6，死單元的單元載荷為0，從而不對載荷向量生效，同樣的，死單元的質量、阻尼、比熱等等參數也設置為0，單元的應力應變也因此為0。
2 前處理
前處理包括單元定義、材料定義和建模，單元定義是需要注意單元屬性，此次定義13號二維耦合單元，具有溫度和位移自由度。
材料屬性包括結構參數和熱參數，具體包含彈性模量，泊松比，屈服強度，塑性屬性，材料密度，熱膨脹系數，熱傳導系數，比熱容。焊接時溫度較高，定義材料通常需要定義多個溫度下的值。
展開 244 基于matlab的實現四節點板單元剛度矩陣求解 ￥24.9
基于matlab的實現四節點板單元剛度矩陣求解，振動模態分析。可自主輸入材料參數，板單元長寬厚尺寸。輸出振型結果，輸出多少階可自主設置。最后以可視化的形式展示，程序已調通，可直接運行。
ANSYS單元生死技術助力牛郎織女來相會
然后按載荷步順序定義單元生死，并結合網格大小給愛心加載位移邊界。
第五步：求解
總結
ANSYS單元的生死功能是通過修改單元剛度的方式實現的。
■ 為了達到讓單元“死”掉的效果，ANSYS程序會給單元剛度乘以一個很小的系數（默認為1.0E-6）。
■ 在載荷向量中，與被“殺死”的單元相聯系的單元荷載被設置為0，其質量、阻尼等一切對計算有影響的參數都會被設置為0。與之相似，當單元“活”的時候，也是通過修改剛度系數的方式實現的。
在較老的版本——R18之前，這些效果都需要通過命令流的方式來實現，但現在可以直接通過菜單來實現了。
以上便是今天鵲橋相會的仿真模擬過程了，希望大家可以通過這個案例發現在仿真中的樂趣~
最后祝大家七夕快樂~
Happy Chinese Valentine's Day
*本文版權歸上海安世亞太所有，如需轉載，請與我們聯系。
關注【上海安世亞太】，獲取更多原創文章、活動資訊
如果你覺得這篇文章對你有用，點個贊吧！
展開 