
發(fā)布
注冊
/
登錄ansys瞬態(tài)求解
關注創(chuàng)建者:王靖雯 創(chuàng)建時間:2023-03-08

ansys瞬態(tài)求解的實例教程
在瞬態(tài)問題的求解中,導數項可以寫成前后時間變量差值與時間間隔的比值:
代入后得到如下形式:
求解思路
在求解過程中,把Tn+1當作未知量,Tn作為已知量。這樣在每個時間點,求解方法和結構有限元方法一致。
初始時候,可以指定一個溫度作為全域已知初始溫度,然后在迭代過程中,Tn和Tn+1會逐漸接近,達到收斂狀態(tài)。
案例效果
設計案例如下,同時包含對流換熱邊界條件和熱流,時間總長10000s,每步時間間隔50s。
自研求解器和商用軟件結果對比如下,從結果可以看出,自研求解器結果與商用軟件結果一致。
自研求解器結果:最終溫度分布
商用軟件結果:最終溫度分布
自研求解器結果:平均溫度時間曲線
商用軟件結果:平均溫度時間曲線
展開 關鍵詞:熱源,瞬態(tài),熱傳導,有限元求解器,三角形單元,自研
在《瞬態(tài)熱傳導有限元求解器開發(fā)》一文中,我們介紹了自研的二維瞬態(tài)熱傳導求解器。
當時那個控制方程沒有考慮熱源,邊界條件中只涉及溫度、熱流、對流。然而在很多問題中,熱源才是最關鍵的邊界條件,比如電發(fā)熱、化學反應生熱。
熱源的處理
熱源是體熱,相對應的熱流是面熱。兩者處理方式類似,都是根據單位熱功率值和幾何尺寸計算熱功率,然后加到控制方程矩陣的右側,承擔類似于結構力學中的“載荷”的功能。
區(qū)別在于,熱源是作用在體上的,單位是W/m3,熱流是作用在面上,單位是W/m2。具體到編程上,熱源要分配到單元的三個節(jié)點上,熱流要分配到單元某個邊的兩個節(jié)點上。
從求解器編程的角度來說,這些邊界條件的處理方式都是固定和通用的。考驗一般出現在實際工程項目中使用自研求解器的時候。
在CAE軟件的開發(fā)中,交互端和求解器端永遠要解決的問題是,如何讓所有單元始終知道:
(1)它是誰?(材料參數,幾何參數);
(2)它在哪?(和其他單元的相對位置);
(3)它怎么了?(邊界條件)。
以熱源為例,在交互界面上,我們通過視口選擇單元,指定其體熱功率。那么前端數據在生成求解器輸入的時候,就要告知求解器所有單元的編號和其對應的體熱功率。
當求解器拿到單元編號以后,就需要索引或者計算其面積,并根據單元三個節(jié)點編號,將功率加到載荷列陣對應的位置。
驗證
設計案例如下,區(qū)域外部為20℃空氣,對流換熱系數取5W/(m2K),時間總長18000s,每步時間間隔60s。
自研求解器得到模型中心最終溫度是84.6℃,與商用軟件結果完全一致。
展開 (6)求解設置,設置求解時間為1s,設置初始子步為10,最小子步為10,最大子步為1000。
(7)載荷設置。設置轉動副載荷,主動輪轉動副載荷設置為轉速;從動輪轉動副載荷設置為轉矩,分別如下圖所示。
(8)計算求解。
提示:如果不收斂,可以通過調試網格質量,調試接觸算法,或者增加一個時間較短的分析步,該分析步用于轉速從0rpm斜坡加速到600rpm,而不是階躍加載,這樣利于收斂。
(9)求解完成之后即可查看結果,等效應力如下圖所示。也可以通過接觸工具查看接觸壓力云圖等。
我們來看下基于<strong style="color: rgb(25, 27, 31);">神工坊?</strong><strong>SIMFORGE?高性能工業(yè)仿真平臺</strong>的CFD瞬態(tài)計算,和其他仿真云平臺進行效率對比如何。</p><h2 class="ql-align-justify"><strong>一、模型與網格</strong></h2><p class="ql-align-justify">采用某品牌空調室外機作為瞬態(tài)分析的仿真模型,左側與后側的進口流域,以及前側的出口流域都考慮到計算中,并對空調內部結構簡化后進行網格劃分,最終<strong>網格單元數868萬</strong>,如下圖所示。其中,風扇葉片的<strong>旋轉速度是850rpm</strong>。</p><p class="ql-align-center"><img src="https://img.jishulink.com/msimage/202301/b604dd3cd04128e9266451ab90624f30.png?image_process=/format,webp/quality,q_40" alt="【仿真平臺性能測試】Fluent旋轉機械瞬態(tài)分析的圖1"></p><h2 class="ql-align-justify"><strong>二、求解設置</strong></h2><p class="ql-align-justify">根據該款旋轉機械的相關參數,經過理論計算得到該旋轉機械的最大速度為25.6m/s,折合馬赫數為0.075,為不可壓縮流動,故選擇壓力基求解器,湍流模型選用了適用于旋轉機械的k-ε Realizable模型。對于動區(qū)域計算模型,本次瞬態(tài)計算選擇了網格區(qū)域移動的滑移網格法,仿真的模擬時間為10s,相關設置如下。
展開 如果不考慮慣性力和阻尼,則可以使用靜力學分析來代替瞬態(tài)動力學分析。對于線性結構,它的瞬態(tài)動力學平衡方程如下:</p><p><br></p><p>在Ansys有限元分析軟件中,式共有三種求解方法分別為:完全法、模態(tài)疊加法和縮減法。完全法和縮減法采用直接積分求解瞬態(tài)動力學平衡方程。而模態(tài)疊加法則使用坐標轉換解耦后開始求解。</p><p><br></p><p>5.1.1 模態(tài)疊加法</p><p>針對模態(tài)疊加法,式中的可寫為:</p><p><br></p><p>式中:</p><p>為節(jié)點力隨時間變化量;</p><p>為關于矢量載荷的比例因子;</p><p>是在模態(tài)分析中的矢量載荷。</p><p>利用模態(tài)坐標表示節(jié)點位移可通過下式得到:</p><p><br></p><p>式中,是第階模態(tài)振型;</p><p>是所要提取的模態(tài)數量。</p><p>根據式可得利用模態(tài)疊加法計算瞬態(tài)動力學問題首先需要進行模態(tài)分析,因為在節(jié)點位移中包含了模態(tài)振型。</p><p>將式帶入可得:</p><p><br></p><p>在式左乘模態(tài)振型可得:</p><p><br></p><p>模態(tài)的正交條件如下:</p><p>將正交條件帶到式則為:</p><p><br></p><p>利用質量矩陣進行歸一化處理,得到系數為:</p><p><br></p><p>的系數為:</p><p><br></p><p>的系數為:</p><p><br></p><p>同樣可以將其寫為:</p><p><br></p><p>利用式、、、帶到式中,得到:</p><p>縮減法用于模態(tài)分析時,涉及到主自由度的選取。因為縮減法通過減少模型的自由度數簡化計算,如果選用QR方法進行模態(tài)分析,模態(tài)坐標系下的運動微分方程寫為:</p><p>5.1.2 直接積分法</p><p>在ANSYS中,隱式方法Newmark和HHT用于求解瞬態(tài)動力學問題。
展開 
ansys瞬態(tài)求解的相關專題、標簽、搜索
ansys瞬態(tài)求解的最新內容
關鍵詞:熱源,瞬態(tài),熱傳導,有限元求解器,三角形單元,自研
在《瞬態(tài)熱傳導有限元求解器開發(fā)》一文中,我們介紹了自研的二維瞬態(tài)熱傳導求解器。
當時那個控制方程沒有考慮熱源,邊界條件中只涉及溫度、熱流、對流。然而在很多問題中,熱源才是最關鍵的邊界條件,比如電發(fā)熱、化學反應生熱。
熱源的處理
熱源是體熱,相對應的熱流是面熱。兩者處理方式類似,都是根據單位熱功率值和幾何尺寸計算熱功率,然后加到控制方程矩陣的右側
概述
這篇文章介紹了:
如何使用 RCWA 求解器分析周期性多層結構(如光子晶體、衍射光柵)的光學響應;
RCWA 求解器的原理:在傅里葉域中劃分均勻層,并通過 S 矩陣雙向傳播計算透射、反射及各個光柵階的功率;
如何設置入射平面波的傳播方向(X/Y/Z 軸)、角度(θ/?)和偏振(s/p),以及反向傳播的兩種模式(鏡像 k 矢量和反向 k 矢量);
對比 RCWA
在CAE領域,選擇Standard(隱式)還是Explicit(顯式)求解器,本質上是在平衡“計算精度”與“時間尺度”。
1?? 隱式求解 (Implicit/Standard)
核心是求解 $Ku=F$。每一步都需要進行矩陣求逆和牛頓迭代,以確保力平衡。
特點: 絕對收斂。步長可以很大,不受穩(wěn)定性限制。
擅長: 靜力學、線性振動
關鍵詞:瞬態(tài),熱傳導,有限元求解器,三角形單元
熱傳遞有三種方式:熱傳導、熱對流、熱輻射。就熱傳導問題而言,無論是結構力學還是流體力學都會涉及,兩邊都沒拿它當外人。
前面的文章提到過,結構力學的有限元發(fā)展地非常成熟,大部分的剛度矩陣在文獻里面都推導好了。而流體力學的很多單元類型的有限元方程,可能需要自行推導完成。在熱傳導問題中,我采用加權余量法進行處理,推導出了符合結構力學有限元文獻中給出的剛度矩陣
概述
PCB 組件在工作時產生的熱量會直接影響其電性能與長期可靠性。過高的溫度或頻繁的溫度波動會引發(fā)材料老化、信號失真,并因材料間熱膨脹系數不匹配而產生熱應力,最終導致焊點開裂、器件失效等故障。因此,評估 PCB 可靠性必須進行瞬態(tài)熱力耦合分析,即先分析動態(tài)溫度場,再計算由此產生的熱應力。
目標
通過高保真建模仿真,系統觀察并量化印刷電路板(PCB)上關鍵元器件在瞬態(tài)熱載荷作用下的力學響應與應力表現
AnsysWB-手機跌落瞬態(tài)仿真
汽水易拉罐壓碎仿真模擬
使用 ANSYS CFX 對離心泵內的流動進行瞬態(tài)仿真。湍流模型采用 SST。同時包含 CFX 定義文件。
附件下載
聯系工作人員獲取附件
概述
本文使用兩個示例演示了如何使用 ZPL 創(chuàng)建用戶自定義解。第一個示例介紹了如何創(chuàng)建 ZPL 解以確保序列文件中像面的曲率半徑等于系統的 Petzval 曲率。第二個示例介紹了如何在非序列元件編輯器(Non-Sequential Component Editor)中基于其他物體的參數來約束的物體位置。
簡介
求解 ( Solve ) 是可以在諸如鏡頭數據編輯器或非序列元件編輯器之類的編輯器中主動調整特定值的功能
本文原刊登于Ansys.com:《Race to Faster Fluent Results with Ansys Gateway Powered by AWS》
作者:Thomas Lejeune | Ansys產品營銷高級經理
編輯整理:郭曉東 | Ansys主任應用工程師
Ansys Fluent用戶需要出色的計算速度和功能來求解大規(guī)模的問題,而他們現在可以利用專用的云平臺


