Abaqus熱分析的案例
Abaqus熱應(yīng)力分析、熱誘導(dǎo)振動(dòng)分析簡(jiǎn)單實(shí)例
1、一個(gè)金屬懸臂梁，一端固支，初始溫度20℃，溫度突變到120℃時(shí)由于膨脹及邊界約束而產(chǎn)生熱應(yīng)力，進(jìn)而引起振動(dòng)，這種振動(dòng)就是熱誘導(dǎo)振動(dòng)。
2、熱誘導(dǎo)振動(dòng)分析的成功應(yīng)用不多見(jiàn)，在哈勃太空望遠(yuǎn)鏡曾因熱誘導(dǎo)振動(dòng)問(wèn)題而發(fā)生故障。現(xiàn)在對(duì)航天器的分析中，熱誘導(dǎo)振動(dòng)屬于難點(diǎn)和重點(diǎn)。國(guó)內(nèi)曾有人對(duì)衛(wèi)星天線做過(guò)準(zhǔn)靜態(tài)熱誘導(dǎo)振動(dòng)分析，也有人對(duì)空間站太陽(yáng)能電池陣的桅桿做過(guò)基于模態(tài)的熱誘導(dǎo)振動(dòng)分析（可能類似Abaqus中的線性攝動(dòng)分析）。
3、熱應(yīng)力分析與熱誘導(dǎo)振動(dòng)分析進(jìn)行耦合分析，還有難度，問(wèn)題是多方面的。下面僅就準(zhǔn)靜態(tài)非耦合的熱誘導(dǎo)振動(dòng)分析為例，介紹由熱應(yīng)力引起的振動(dòng)。
4、懸臂梁材料屬性：
Conductity: 300W/(mK)
Density: 3000kg/m3
Elastic: E=3e10Pa, ν=0.3
Expansion: 3e-5 K-1
Specific Heat: 300J/(kgK)
5、分析結(jié)果
6、詳細(xì)步驟
見(jiàn)附件。
Abaqus熱應(yīng)力分析、熱誘導(dǎo)振動(dòng)分析簡(jiǎn)單實(shí)例-kxh.part4.rar
Abaqus熱應(yīng)力分析、熱誘導(dǎo)振動(dòng)分析簡(jiǎn)單實(shí)例-kxh.part1.rar
Abaqus熱應(yīng)力分析、熱誘導(dǎo)振動(dòng)分析簡(jiǎn)單實(shí)例-kxh.part2.rar
Abaqus熱應(yīng)力分析、熱誘導(dǎo)振動(dòng)分析簡(jiǎn)單實(shí)例-kxh.part3.rar
展開(kāi) Abaqus熱應(yīng)力分析、熱誘導(dǎo)振動(dòng)分析簡(jiǎn)單實(shí)例
1、一個(gè)金屬懸臂梁，一端固支，初始溫度20℃，溫度突變到120℃時(shí)由于膨脹及邊界約束而產(chǎn)生熱應(yīng)力，進(jìn)而引起振動(dòng)，這種振動(dòng)就是熱誘導(dǎo)振動(dòng)。
2、熱誘導(dǎo)振動(dòng)分析的成功應(yīng)用不多見(jiàn)，在哈勃太空望遠(yuǎn)鏡曾因熱誘導(dǎo)振動(dòng)問(wèn)題而發(fā)生故障。現(xiàn)在對(duì)航天器的分析中，熱誘導(dǎo)振動(dòng)屬于難點(diǎn)和重點(diǎn)。國(guó)內(nèi)曾有人對(duì)衛(wèi)星天線做過(guò)準(zhǔn)靜態(tài)熱誘導(dǎo)振動(dòng)分析，也有人對(duì)空間站太陽(yáng)能電池陣的桅桿做過(guò)基于模態(tài)的熱誘導(dǎo)振動(dòng)分析（可能類似Abaqus中的線性攝動(dòng)分析）。
3、熱應(yīng)力分析與熱誘導(dǎo)振動(dòng)分析進(jìn)行耦合分析，還有難度，問(wèn)題是多方面的。下面僅就準(zhǔn)靜態(tài)非耦合的熱誘導(dǎo)振動(dòng)分析為例，介紹由熱應(yīng)力引起的振動(dòng)。
4、懸臂梁材料屬性：
Conductity: 300W/(mK)
Density: 3000kg/m3
Elastic: E=3e10Pa, ν=0.3
Expansion: 3e-5 K-1
Specific Heat: 300J/(kgK)
詳細(xì)步驟
見(jiàn)附件。
Abaqus熱應(yīng)力分析、熱誘導(dǎo)振動(dòng)分析簡(jiǎn)單實(shí)例-kxh.part4.rar
Abaqus熱應(yīng)力分析、熱誘導(dǎo)振動(dòng)分析簡(jiǎn)單實(shí)例-kxh.part2.rar
Abaqus熱應(yīng)力分析、熱誘導(dǎo)振動(dòng)分析簡(jiǎn)單實(shí)例-kxh.part3.rar
Abaqus熱應(yīng)力分析、熱誘導(dǎo)振動(dòng)分析簡(jiǎn)單實(shí)例-kxh.part1.rar
展開(kāi) abaqus輻射熱分析（面與面之間的輻射） ￥20
e.gif
abaqus輻射熱分析，面與面之間的輻射。
abaqus橡膠熱仿真：減振橡膠疲勞黏滯生熱的仿真分析-源文件與子程序詳解
得到了生熱率與溫度、載荷頻率和應(yīng)變幅值的函數(shù)關(guān)系式。
利用依黏彈性理論得出的黏滯生熱率與溫度、載荷頻率和應(yīng)變幅值的函數(shù)關(guān)系式，編制了相應(yīng)的計(jì)算程序。建立了減振橡膠疲勞黏滯生熱的有限元分析方法。
通過(guò)將經(jīng)典疲勞模型中用作疲勞壽命預(yù)測(cè)指標(biāo)的最大主應(yīng)變替換為穩(wěn)態(tài)溫升，在冪律模型的基礎(chǔ)上開(kāi)發(fā)了一種方法來(lái)快速評(píng)估橡膠結(jié)構(gòu)的疲勞壽命。
08
—
源文件與操作步驟（沙漏試樣為例）
8.1分析流程
仿真分析主要包括三個(gè)環(huán)節(jié)：變形分析、熱源計(jì)算與熱分析。（1）在變形分析環(huán)節(jié)，對(duì)材料和減振元件施加設(shè)定的載荷歷史，采用超彈性本構(gòu)描述橡膠材料的力學(xué)行為，求解每個(gè)加載時(shí)刻有限元模型中各積分點(diǎn)的應(yīng)變狀態(tài)；（2）在熱源計(jì)算環(huán)節(jié)，對(duì)應(yīng)每一加載時(shí)刻，將變形分析中對(duì)應(yīng)的載荷頻率、應(yīng)變狀態(tài)（動(dòng)態(tài)應(yīng)變幅值）以及熱分析中得到的溫度作為輸入變量，通過(guò)自編的Fortran語(yǔ)言子程序，計(jì)算得到各積分點(diǎn)的黏滯生熱率；（3）依已知的材料參數(shù)和問(wèn)題的熱邊界條件進(jìn)行Abaqus熱分析，得出溫度分布后再將溫度場(chǎng)數(shù)據(jù)返回到自編子程序，對(duì)黏滯生熱強(qiáng)度和溫度場(chǎng)進(jìn)行迭代計(jì)算，從而得出橡膠材料和減振元件各位置的溫升歷程。
展開(kāi) 
基于abaqus的土壤熱硬化分析
基于abaqus的土壤熱硬化分析
一、問(wèn)題描述
此問(wèn)題為太沙基固結(jié)問(wèn)題考慮溫度的熱固結(jié)問(wèn)題，具體問(wèn)題描述內(nèi)容此處不做討論。
ANSYS與ABAQUS比較之實(shí)例8---帶孔平板的熱應(yīng)力分析1
本博文是關(guān)于ANSYS與ABAQUS比較之系列博文，本例子使用ABAQUS做熱應(yīng)力分析，后面會(huì)使用ANSYS對(duì)同一個(gè)問(wèn)題做熱應(yīng)力分析。
【問(wèn)題描述】
一個(gè)帶孔平板結(jié)構(gòu)如下圖
該平板上邊沿固定，左右兩邊是滾動(dòng)支座支撐。該板的初始溫度是25度，現(xiàn)在要求當(dāng)溫度升高到150度時(shí)，板中的應(yīng)力分布。
已知：材料的彈性模量是2e9pa, 泊松比是0.3,熱膨脹系數(shù)是1.35e-5/度。
【問(wèn)題分析】
1.
分析類型。這是一個(gè)平面應(yīng)力問(wèn)題，應(yīng)力的產(chǎn)生是因?yàn)闇囟鹊淖兓瘜?dǎo)致產(chǎn)生了熱應(yīng)變，而該熱應(yīng)變又被約束限制導(dǎo)致熱應(yīng)力的產(chǎn)生。
2.
非線性考慮。只有一個(gè)物體，不存在接觸非線性；沒(méi)有材料非線性；沒(méi)有幾何非線性。總之，這就是一個(gè)最簡(jiǎn)單的線彈性分析。
3. 幾何建模。由于該結(jié)構(gòu)左右對(duì)稱，只取一半研究。
4.
邊界條件。除了常規(guī)的位移邊界條件以外，對(duì)該板施加預(yù)定義溫度場(chǎng)25度，而在第一個(gè)分析步修改該溫度場(chǎng)的溫度為150度。
【求解過(guò)程】
1. 創(chuàng)建部件
只取一半建模，它是一個(gè)二維的可變形部件。
2. 定義材料屬性
只需要定義彈性模量，泊松比及線膨脹系數(shù)。
3. 定義截面屬性
創(chuàng)建均質(zhì)的實(shí)體截面，并將上述材料屬性賦予給它，然后將該截面屬性賦予給前面的部件。
4. 裝配部件
唯一的部件，導(dǎo)入到裝配即可。
5.設(shè)置分析步
兩個(gè)分析步。
新創(chuàng)建的分析步是最一般的靜力學(xué)通用分析步。
6.定義載荷和邊界條件
首先定義位移邊界條件，在初始分析步中，固定上邊，左右兩邊施加X(jué)方向的位移限制。
使用預(yù)定義場(chǎng)確定溫度。
對(duì)整塊板施加25度的初始溫度。
展開(kāi) Abaqus熱分析關(guān)于結(jié)果中溫度的后處理二次開(kāi)發(fā) ￥5
在abaqus中，進(jìn)行熱分析的時(shí)候，經(jīng)常要對(duì)得到的溫度結(jié)果進(jìn)行處理。如下圖所示：
我們要對(duì)這個(gè)面上的溫度數(shù)據(jù)，進(jìn)行處理，比如：
（1）我要找出這個(gè)面上溫度符合某一條件的區(qū)域。
（2）我要找出某一時(shí)刻，這一面上所有節(jié)點(diǎn)溫度的平均值。
要解決以上兩個(gè)問(wèn)題，最好的辦法就是用python對(duì)ABAQUS的ODB文件進(jìn)行訪問(wèn)處理。用python來(lái)處理ODB文件中的歷史輸出數(shù)據(jù)。我寫了一個(gè)函數(shù)，用以解決以上兩個(gè)方面的問(wèn)題。大家可以從代碼得到啟發(fā)，思考如何對(duì)abaqus的歷史輸出數(shù)據(jù)進(jìn)行訪問(wèn)。有相關(guān)ABAQUS二次開(kāi)發(fā)方面的問(wèn)題，我們可以探討探討。QQ:1917665644
具體代碼如下：
展開(kāi) Abaqus熱應(yīng)力分析和膨脹節(jié)應(yīng)用
本文將用Abaqus Standard演示一個(gè)熱分析案例。該分析將含有圓柱形殼體的耦合熱應(yīng)力問(wèn)題（例如工廠中使用的承裝高溫流體的管道）。該管道將與一個(gè)金屬伸縮接頭連接，這種接頭將用于承受管道的熱延伸。耦合分析的目的是通過(guò)預(yù)先設(shè)置的場(chǎng)來(lái)論證結(jié)果值的對(duì)應(yīng)關(guān)系。
熱和結(jié)構(gòu)域間的耦合計(jì)算
Abaqus為分析中的熱和結(jié)構(gòu)域間的耦合計(jì)算提供了許多建模方法。
通常會(huì)使用兩種方法：
1） 使用耦合溫度-位移元件（SAX2T，S8RT，C3D20RT等），在一次分析中結(jié)合熱負(fù)荷和位移。
2） 在第一次分析中對(duì)結(jié)構(gòu)的熱負(fù)荷建模，然后在隨后的應(yīng)力分析中映射該熱負(fù)荷（通過(guò)節(jié)點(diǎn)溫度，即NT）。
本文將演示第二種方法，在很多情況下，若考慮不同類型負(fù)載的所有影響（軸向壓力，環(huán)向壓力，熱負(fù)荷等），第一種方法可能更可行。
建模假定條件
1） 模型將包含一個(gè)圓柱形殼體結(jié)構(gòu)，一個(gè)承裝高溫流體的管道。對(duì)于這種分析，由于管道內(nèi)外兩側(cè)的溫度差異，所以只考慮熱負(fù)荷。管道將被認(rèn)為具有足夠的剛性，在流體壓力下不會(huì)周向變形（在靜態(tài)分析只有軸向平移允許的）。
2） 為了經(jīng)行分析時(shí)清晰，使用了圓柱形殼體模型（完整圓柱體）。然而，建議盡可能經(jīng)常使用軸對(duì)稱模型（或者具有適當(dāng)邊界條件的單個(gè)元件）以減少運(yùn)行時(shí)間。
3） 將考慮穩(wěn)態(tài)熱傳導(dǎo)分析。
4） 將考慮靜力通用分析。
5） 膨脹節(jié)和流體管的材料是線性彈性的。在現(xiàn)實(shí)生活中，材料不一定相同。每種的應(yīng)用都有所不同，但每當(dāng)非金屬膨脹節(jié)用于低壓應(yīng)用時(shí)，通常比其連接部件更加靈活（橡膠類材料）。
6） 壓力結(jié)果以兆帕（Mpa）為單位給出；位移結(jié)果以毫米（mm）為單位給出。
本研究中，將進(jìn)行兩項(xiàng)分析。首先是內(nèi)部溫度為200攝氏度，外部（環(huán)境）溫度為10攝氏度的管道上的熱傳導(dǎo)分析。這種溫差會(huì)導(dǎo)致管道沿周向和軸向膨脹。
展開(kāi) 針對(duì)某袋除塵器整體進(jìn)行ABAQUS有限元分析，考慮九項(xiàng)載荷工況，分析設(shè)備靜應(yīng)力、熱應(yīng)力、變形及熱膨脹數(shù)值 ￥15
某袋除塵殼體結(jié)構(gòu)選型如下：
箱體板厚5mm
箱體角柱：角鋼L90*56*8
箱體加強(qiáng)筋：角鋼L90*56*6
花板厚6mm
花板下加強(qiáng)筋：橫向?yàn)楸怃?0*6，縱向?yàn)楸怃?00*6
箱體中間支撐管：鋼管Φ60*5
圖1 袋除塵殼體結(jié)構(gòu)示意圖
2、 建立模型
按照殼體結(jié)構(gòu)示意圖建立幾何模型如圖2所示。
圖2 建立幾何模型
三、約束條件及載荷
立柱底部約束如圖3所示。
圖3 立柱底部邊界約束
載荷：
（1）自重（軟件考慮）；
（2） 頂部載荷：檢修載（按400kg/m2）；
（3） 花板處載荷：濾袋、濾籠、濾袋積灰（積灰厚度按5mm）共3.06t；
（4） 灰斗積灰重：滿灰9.6t；
（5） 保溫載荷：按25kg/m2；
（6） 負(fù)壓11000Pa或正壓8000Pa兩種工況分別施加；
（7） 煙道及檢修平臺(tái)載荷：上煙道（出氣端）900kg，下煙道（進(jìn)氣端）
400kg，上中下三層檢修平臺(tái)檢修載荷均為400×2.85×3.25=3705kg。
注：此項(xiàng)載荷殼體和鋼支架各占一半。
（8） 灰斗卸灰口載荷（方向按照幾何模型坐標(biāo)系）：FX=4700N，F(xiàn)Y=3500N，F(xiàn)Z=-4700N，MX=3690N.m，MY=4800N.m，MZ=5540N.m。
（9） 頂部牛腿處檢修荷載：?jiǎn)蝹€(gè)牛腿處載荷為1t，頂板為260×260，轉(zhuǎn)化為面壓添加，面壓為1×10×1000/260/260=0.148N/mm2。
下圖4所示為載荷添加圖示：
（a）負(fù)壓11000Pa （b）正壓8000Pa （c）花板處載荷
展開(kāi) ABAQUS培訓(xùn)案例之熱分析-熱輻射
圖1模型示意圖
今天給大家分享的是熱輻射分析。如圖1所示，模型由2個(gè)part組成，一個(gè)是fin，其周期對(duì)稱性在設(shè)置輻射換熱時(shí)可以設(shè)置，另一個(gè)代表周邊環(huán)境的ambient。模型先進(jìn)行了穩(wěn)態(tài)分析，然后建立2個(gè)瞬態(tài)分析步，實(shí)現(xiàn)環(huán)境溫度為800時(shí)對(duì)fin加熱的過(guò)程，和環(huán)境溫度38時(shí)fin部件的冷卻過(guò)程。當(dāng)然除了ambient和fin的輻射換熱，ambient和fin也分別建立了Surface film condition換熱。下面詳解每個(gè)步驟的設(shè)置。
目標(biāo)：輻射換熱與對(duì)流換熱設(shè)置，cavity radiation應(yīng)用。
材料：材料參數(shù)定義了Density為7800，Conductivity和SpecificHeat分別為50和500（SI單位制）。
分析步設(shè)置：本案例設(shè)置了三個(gè)分析步，step-1為熱穩(wěn)態(tài)分析步，step-2和step-3為瞬態(tài)分析步，如下圖所示。歷史輸出設(shè)置output三個(gè)節(jié)點(diǎn)的溫度輸出。
圖2 分析步設(shè)置
相互作用設(shè)置：定義了三個(gè)換熱條件，設(shè)置bot面換熱系數(shù)Surfacefilm condition為2500，熱沉溫度100，srfs面換熱系數(shù)Surface film condition為10，step-1熱沉溫度38，step-2時(shí)為800，step-3時(shí)改為38。設(shè)置srfs和samb之間的輻射換熱Cavity radiation，發(fā)射系數(shù)為0.8和1，并設(shè)置Symmetry對(duì)稱，如下圖所示。
圖3 換熱設(shè)置
邊界條件：設(shè)置所有區(qū)域初始溫度為77，ambient的溫度step-1時(shí)38，step-2時(shí)為800（加熱過(guò)程），step-3時(shí)改為38（冷卻過(guò)程）。
展開(kāi) Abaqus熱傳導(dǎo)與熱應(yīng)力分析基礎(chǔ)知識(shí)介紹
熱傳遞的分析目標(biāo)是研究熱量的傳遞過(guò)程。熱傳遞分析以熱變量或與熱相關(guān)的變量的形式來(lái)計(jì)算熱響應(yīng)，如溫度分布和溫度梯度以及熱通量。
熱傳遞分析包括兩種類型，第一種，非耦合的熱響應(yīng)，即純熱傳遞分析；第二種耦合的響應(yīng)（熱-應(yīng)力分析），分為順序耦合和完全耦合。純熱傳遞分析在Abaqus/Standard中完成，耦合響應(yīng)在Abaqus/Standard和Abaqus/Explicit中完成。
熱傳遞包括三種模式：
傳導(dǎo)，也被稱為“實(shí)體熱傳遞”，發(fā)生在物體內(nèi)的分子水平上，金屬是典型的熱的良導(dǎo)體，氣體則不是。
對(duì)流，是通過(guò)熱物質(zhì)（氣體或者流體）的流動(dòng)進(jìn)行熱量傳遞，包括自然對(duì)流和強(qiáng)制對(duì)流，如水泵、風(fēng)機(jī)或其他壓差作用引起的對(duì)流。
輻射，即電磁輻射，發(fā)生不需要介質(zhì)，真空中亦可。
熱傳遞可以上述一種或幾種模式的組合來(lái)進(jìn)行。在熱傳遞分析中用到的基本量有以下這些，如圖所示。
abaqus-復(fù)合材料仿真分析基礎(chǔ)篇.pdf
展開(kāi) 
希望對(duì)大家有幫助
觀看須知.txt
Abaqus熱應(yīng)力分析、熱誘導(dǎo)振動(dòng)分析簡(jiǎn)單實(shí)例-kxh.part1.rar
Abaqus熱應(yīng)力分析、熱誘導(dǎo)振動(dòng)分析簡(jiǎn)單實(shí)例-kxh.part1.rar
觀看須知.txt
Abaqus熱應(yīng)力分析、熱誘導(dǎo)振動(dòng)分析簡(jiǎn)單實(shí)例-kxh.part1.rar
ABAQUS熱應(yīng)力分析 附ABAQUS中初始地應(yīng)力的施加下載
熱應(yīng)力分析過(guò)程
ABAQUS 提供三種熱應(yīng)力分析程序：
1. 順序耦合熱應(yīng)力分析，最常用的方法
? 當(dāng)應(yīng)力是由熱量場(chǎng)存在造成的，并且熱求解過(guò)程與應(yīng)力狀態(tài)無(wú)關(guān)，也就是說(shuō)應(yīng)力依賴于熱產(chǎn)生，而熱并不依賴位移。
? 需要跑兩個(gè)分析: 先分析熱傳導(dǎo)，再將溫度結(jié)果導(dǎo)熱應(yīng)力分析
? 熱分析的結(jié)果，如溫度(位置，時(shí)間的函數(shù))被讀入應(yīng)力分析，作為一個(gè)預(yù)定義場(chǎng)。
2. 完全耦合熱應(yīng)力分析，最常用的方法
? 應(yīng)力依賴于溫度場(chǎng)并且溫度也依賴于應(yīng)力場(chǎng)。
? 只需要跑一個(gè)析。
3.
展開(kāi) abaqus梁?jiǎn)卧?em>熱分析？
各位大佬，請(qǐng)問(wèn)abaqus里三維梁?jiǎn)卧獰o(wú)法進(jìn)行heat transfer分析嗎，在mesh類型里，選擇heat transfer后單元類型的后面沒(méi)有增加T，計(jì)算的時(shí)候也會(huì)報(bào)錯(cuò)，顯示單元缺少屬性定義。
Abaqus穩(wěn)態(tài)熱分析實(shí)例
模型尺寸如下圖所示，
熱傳導(dǎo)系數(shù) = 52 W/m/°C比熱= 434 J/kg/°C
密度 = 7832 kg/m3
對(duì)流換熱系數(shù) = 750 W/m2/°C
邊界條件:
AB邊溫度 q = 100°C （固定溫度）熱流 DA= 0（絕熱邊）
BC 和 CD 邊與周圍介質(zhì)對(duì)流換熱周圍介質(zhì)溫度為0°C
本例采用國(guó)際單位制。E點(diǎn)理論計(jì)算結(jié)果為18.3度，下面使用Abaqus來(lái)計(jì)算并驗(yàn)證。
首先創(chuàng)建2D、shell的幾何模型，其次是材料參數(shù)的設(shè)置，與靜力分析不同，熱傳導(dǎo)需要設(shè)置熱傳導(dǎo)系數(shù)，在Mechanical>Conductivity里輸入52。本例是穩(wěn)態(tài)熱分析，因此只需要這一個(gè)參數(shù)，若為瞬態(tài)熱分析，則還需要比熱以及密度值。
其次進(jìn)行分析步設(shè)置，這一步與靜力分析也有所不同，選擇Gerneal>Heat Transfer作為分析步。默認(rèn)的為瞬態(tài)響應(yīng)，這里選擇穩(wěn)態(tài)分析，同靜力分析一樣，這里的時(shí)間1沒(méi)有真實(shí)含義，保持默認(rèn)。增量步設(shè)置與靜力分析一樣。
邊界條件由默認(rèn)的Mechanical改為Other，選擇溫度。選擇模型最下面的邊，這里定義為bottom集，給予100度的溫度。
下面設(shè)置模型與周圍空氣的對(duì)流。模型右面的邊（side）與上面的邊（top）與周圍環(huán)境發(fā)生熱交換，對(duì)流系數(shù)為750，Sink temperature為周圍環(huán)境的溫度，這里給0。
Mesh模塊中，需要將單元族改為Heat Transfer，確認(rèn)使用的是DC2D4單元。至此，熱分析的設(shè)置已經(jīng)完成。可以提交計(jì)算。在后處理中查詢右邊界從下網(wǎng)上0.2m處的溫度值為18.4151，與理論計(jì)算結(jié)果18.3相差不大。右圖為對(duì)模型網(wǎng)格加密的結(jié)果，顯示溫度值為18.29，接近理論解。
abaqus穩(wěn)態(tài)傳熱分析實(shí)例.pdf
展開(kāi) 