FEM之求解器(12)---利用C++開發(fā)熱學(xué)求解器1
2017-01-15 by:CAE仿真在線 來源:互聯(lián)網(wǎng)
本文主要介紹用C++開發(fā)的
2D/3D熱學(xué)求解器FasThermal,一可以單獨(dú)求解熱學(xué)問題,二可以作為多物理場求解器的基礎(chǔ)。FasThermal目前有以下簡單功能:
1. 支持三角形和四面體單元
2. 支持穩(wěn)態(tài)的溫度場求解
3. 支持3種熱荷載方式
4.
求解溫度場的結(jié)果可以在FasView里顯示
熱學(xué)中的剛度矩陣數(shù)據(jù)比力學(xué)中少,熱的形函數(shù)主要是溫度場,溫度是標(biāo)量,1個節(jié)點(diǎn)1個自由度,不管是1D和3D單元,節(jié)點(diǎn)都有相同的自由度。因為簡單,所以一般開源力學(xué)求解器里都會有溫度場的計算。只不過溫度場的邊界多一些,常用的三類邊界條件都會用到。
熱傳遞的方式:
1. 傳導(dǎo):接觸物體之間的能量交換或者 一個物體內(nèi)部由于溫度梯度引起的內(nèi)部能量交換
2. 對流:在物體和周圍介質(zhì)之間發(fā)生的熱交換
3. 輻射:一個物體或者兩個物體之間通過電磁波進(jìn)行的能量交換
幾種常用的熱荷載方式:
1. DOF約束:指定的節(jié)點(diǎn)上的溫度數(shù)值
2. 集中荷載 :集中荷載(熱流)施加在點(diǎn)上
3. 面荷載:在面上的分布荷載(對流,熱流)
4. 體荷載:體積或區(qū)域荷載(導(dǎo)體自發(fā)熱)
為簡化,例子使用平面問題的熱傳遞問題。一個長方形,左邊邊上溫度保持100度,
右邊溫度保持0度,上下邊絕熱,求平面內(nèi)溫度的分布。
采用三角形單元。三角單元介紹和三類邊界條件參考FEM之單元(1)---三角單元介紹
具體熱矩陣公式的推導(dǎo)可以參考 孔祥謙的《有限單元法在傳熱學(xué)中的應(yīng)用》,目前絕版了,只能找到電子書。
1.
FaModel中建里長方形,設(shè)置邊界溫度和導(dǎo)熱系數(shù)。
2. 導(dǎo)出為 FasThermal
求解格式.
3. FasThermal求解
FasThermal求解結(jié)果
中間兩節(jié)點(diǎn)計算結(jié)果 T = 50
與理論值一樣
小結(jié):
1.
三角單元有三個點(diǎn),因此每個單元的熱剛度矩陣為3*3,與此對應(yīng)的溫度向量和溫度荷載為3*1向量。數(shù)據(jù)量比力學(xué)的要少很多。對于體單元,矩陣數(shù)據(jù)也很方便計算
2. 總剛的組裝方法與力學(xué)一樣。
3. 溫度隨時間變化的情形為瞬態(tài)分析,求解思想也和力學(xué)中的動力學(xué)一樣。
4. 例子中使用的最簡單的第一類邊界
5.
熱傳導(dǎo)FEM相對簡單,主要作為后面介紹的熱-結(jié)構(gòu),熱-固-流,熱-結(jié)構(gòu)-電磁 耦合的基礎(chǔ)。
我的更多文章:
- FEM之優(yōu)化算法(1)(2014-10-30 10:53:00)
-
FEM之單元(1)---三角單元介紹
- FEM之求解器加速(3)---線性方程組求解庫大全(最全介紹)(2014-10-04 09:47:04)
-
FEM之求解器加速(2)---顯示解法求解加速
-
FEM之求解器(10)---利用C++開發(fā)動力學(xué)求解器(后處理)
(2014-09-10 13:22:40) -
FEM之求解器(10)---利用C++開發(fā)動力學(xué)求解器
(2014-09-07 20:39:18) - FEM之求解器(11)---利用C++開發(fā)動力學(xué)求解器(2014-09-03 20:33:31)
-
FEM之求解器(9)---利用C++開發(fā)動力學(xué)求解器
-
FEM之求解器(6)---利用C++開發(fā)求解器1
開放分享:優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才
相關(guān)標(biāo)簽搜索:FEM之求解器(12)---利用C++開發(fā)熱學(xué)求解器1 Ansys有限元培訓(xùn) Ansys workbench培訓(xùn) ansys視頻教程 ansys workbench教程 ansys APDL經(jīng)典教程 ansys資料下載 ansys技術(shù)咨詢 ansys基礎(chǔ)知識 ansys代做 Fluent、CFX流體分析 HFSS電磁分析 Abaqus培訓(xùn)
編輯