傳熱學(xué)課件第四章

1、第四章導(dǎo)熱問題的數(shù)值解法第四章導(dǎo)熱問題的數(shù)值解法4-14-1導(dǎo)熱問題數(shù)值求解基本思想導(dǎo)熱問題數(shù)值求解基本思想4-24-2內(nèi)節(jié)點(diǎn)離散方程的建立內(nèi)節(jié)點(diǎn)離散方程的建立 4-3 4-3 邊界節(jié)點(diǎn)離散方程的建立及代數(shù)邊界節(jié)點(diǎn)離散方程的建立及代數(shù) 方程的求解方程的求解1 1 、重點(diǎn)內(nèi)容：、重點(diǎn)內(nèi)容： 掌握導(dǎo)熱問題數(shù)值解法的基本思路；掌握導(dǎo)熱問題數(shù)值解法的基本思路； 利用熱平衡法和泰勒級(jí)數(shù)展開法建立利用熱平衡法和泰勒級(jí)數(shù)展開法建立節(jié)點(diǎn)的離散方程。節(jié)點(diǎn)的離散方程。2 2 、掌握內(nèi)容：、掌握內(nèi)容：數(shù)值解法的實(shí)質(zhì)。數(shù)值解法的實(shí)質(zhì)。 求解導(dǎo)熱問題的三種基本方法求解導(dǎo)熱問題的三種基本方法：(1)(1)實(shí)驗(yàn)法實(shí)驗(yàn)法;

2、(2); (2)理論分析法；理論分析法；(3)(3)數(shù)值計(jì)算法數(shù)值計(jì)算法三種方法的特點(diǎn)三種方法的特點(diǎn)實(shí)驗(yàn)法實(shí)驗(yàn)法: : 是傳熱學(xué)的基本研究方法。是傳熱學(xué)的基本研究方法。 a a 適應(yīng)性不好；適應(yīng)性不好； b b 費(fèi)用昂貴費(fèi)用昂貴分析法分析法: : a a 能獲得所研究問題的精確解，可以為實(shí)驗(yàn)和數(shù)能獲得所研究問題的精確解，可以為實(shí)驗(yàn)和數(shù)值計(jì)算提供比較依據(jù)；值計(jì)算提供比較依據(jù)；b b 局限性很大，對(duì)復(fù)雜的局限性很大，對(duì)復(fù)雜的問題無法求解；問題無法求解；c c 分析解具有普遍性，各種情況分析解具有普遍性，各種情況的影響清晰可見的影響清晰可見 數(shù)值計(jì)算法數(shù)值計(jì)算法 有效解決復(fù)雜問題的方法；是具有一定精

3、度的近有效解決復(fù)雜問題的方法；是具有一定精度的近似方法。在很大程度上彌補(bǔ)了分析法的缺點(diǎn)，適似方法。在很大程度上彌補(bǔ)了分析法的缺點(diǎn)，適應(yīng)性強(qiáng)，特別對(duì)于復(fù)雜問題更顯其優(yōu)越性；與實(shí)應(yīng)性強(qiáng)，特別對(duì)于復(fù)雜問題更顯其優(yōu)越性；與實(shí)驗(yàn)法相比成本低。驗(yàn)法相比成本低。數(shù)值解法：數(shù)值解法： 有限差分法（有限差分法（finite-differencefinite-difference） 有限元法（有限元法（finite-lementfinite-lement） 邊界元法（邊界元法（boundary-elementboundary-element） 分子動(dòng)力學(xué)模擬（分子動(dòng)力學(xué)模擬（MDMD） 分析解法與數(shù)值解法的異同點(diǎn)

4、：分析解法與數(shù)值解法的異同點(diǎn)： 相同點(diǎn)：相同點(diǎn)：根本目的是相同的，即確定：根本目的是相同的，即確定： t=f(x,y,z) t=f(x,y,z) ； 熱流量。熱流量。 不同點(diǎn)：不同點(diǎn)：數(shù)值解法求解的是區(qū)域或時(shí)間空間坐標(biāo)數(shù)值解法求解的是區(qū)域或時(shí)間空間坐標(biāo)系中離散點(diǎn)的溫度分布代替連續(xù)的溫度場；分析解系中離散點(diǎn)的溫度分布代替連續(xù)的溫度場；分析解法求解的是連續(xù)的溫度場的分布特征，而不是分散法求解的是連續(xù)的溫度場的分布特征，而不是分散點(diǎn)的數(shù)值。點(diǎn)的數(shù)值。 對(duì)物理問題進(jìn)行數(shù)值解法的基本思路可以概對(duì)物理問題進(jìn)行數(shù)值解法的基本思路可以概括為：把原來在時(shí)間、空間坐標(biāo)系中連續(xù)的物理量括為：把原來在時(shí)間、空間坐標(biāo)系

5、中連續(xù)的物理量的場，如導(dǎo)熱物體的溫度場等，用有限個(gè)離散點(diǎn)上的場，如導(dǎo)熱物體的溫度場等，用有限個(gè)離散點(diǎn)上的值的集合來代替，通過求解按一定方法建立起來的值的集合來代替，通過求解按一定方法建立起來的關(guān)于這些值的代數(shù)方程，來獲得離散點(diǎn)上被求物的關(guān)于這些值的代數(shù)方程，來獲得離散點(diǎn)上被求物理量的值，該方法稱為理量的值，該方法稱為數(shù)值解法數(shù)值解法。 這些離散點(diǎn)上被求物理量值的集合稱為該物理量這些離散點(diǎn)上被求物理量值的集合稱為該物理量的的數(shù)值解數(shù)值解。 4-1 4-1 導(dǎo)熱問題數(shù)值求解的基本思想導(dǎo)熱問題數(shù)值求解的基本思想4.1.1 4.1.1 基本思想基本思想建立控制方程及定解條件建立控制方程及定解條件確定節(jié)

6、點(diǎn)（區(qū)域離散化）確定節(jié)點(diǎn)（區(qū)域離散化）建立節(jié)點(diǎn)物理量的代數(shù)方程建立節(jié)點(diǎn)物理量的代數(shù)方程設(shè)立溫度場的迭代初值設(shè)立溫度場的迭代初值求解代數(shù)方程求解代數(shù)方程是否收斂是否收斂解的分析解的分析改進(jìn)初場改進(jìn)初場是是否否4.1.2 4.1.2 物理問題的數(shù)值求解過程物理問題的數(shù)值求解過程0tyf3thf2thf1thx二維矩形域內(nèi)穩(wěn)態(tài)無內(nèi)熱二維矩形域內(nèi)穩(wěn)態(tài)無內(nèi)熱源，常物性的導(dǎo)熱問題源，常物性的導(dǎo)熱問題2 2 例題條件例題條件（a）（1 1）建立控制方程及定解條件）建立控制方程及定解條件 控制方程（即控制方程（即導(dǎo)熱微分方程導(dǎo)熱微分方程） 22220ttxy0tyf3thf2thf1thx二維矩形域內(nèi)無內(nèi)熱源

7、、穩(wěn)態(tài)、常物性的導(dǎo)熱問二維矩形域內(nèi)無內(nèi)熱源、穩(wěn)態(tài)、常物性的導(dǎo)熱問題采用數(shù)值解法的步驟：題采用數(shù)值解法的步驟：（2 2）區(qū)域離散化（確立節(jié)點(diǎn)）區(qū)域離散化（確立節(jié)點(diǎn)） 用一系列與坐標(biāo)軸平行的網(wǎng)格線把求解區(qū)域用一系列與坐標(biāo)軸平行的網(wǎng)格線把求解區(qū)域劃分成若干個(gè)子區(qū)域，用網(wǎng)格線的交點(diǎn)作為劃分成若干個(gè)子區(qū)域，用網(wǎng)格線的交點(diǎn)作為需要確定溫度值的空間位置，稱為需要確定溫度值的空間位置，稱為節(jié)點(diǎn)節(jié)點(diǎn) ( ( 結(jié)結(jié)點(diǎn)點(diǎn) ) ) ，節(jié)點(diǎn)的位置用該節(jié)點(diǎn)在兩個(gè)方向上的，節(jié)點(diǎn)的位置用該節(jié)點(diǎn)在兩個(gè)方向上的標(biāo)號(hào)標(biāo)號(hào) m m ， n n 表示。表示。 相鄰兩節(jié)點(diǎn)間的距離相鄰兩節(jié)點(diǎn)間的距離稱稱步長步長。xyxynm(m,n)MN

8、（b）xyxynm(m,n)MN基本概念：網(wǎng)格線、節(jié)點(diǎn)、界面線、步長、基本概念：網(wǎng)格線、節(jié)點(diǎn)、界面線、步長、控制容積控制容積二維矩形二維矩形域內(nèi)穩(wěn)態(tài)域內(nèi)穩(wěn)態(tài)無內(nèi)熱源，無內(nèi)熱源，常物性的常物性的導(dǎo)熱問題導(dǎo)熱問題 （3 3）建立節(jié)點(diǎn)物理量的代數(shù)方程（離散方程）建立節(jié)點(diǎn)物理量的代數(shù)方程（離散方程） 節(jié)點(diǎn)上物理量的代數(shù)方程稱離散方程。節(jié)點(diǎn)上物理量的代數(shù)方程稱離散方程。 首先劃分各節(jié)點(diǎn)的類型；首先劃分各節(jié)點(diǎn)的類型； 其次，建立節(jié)點(diǎn)離散方程；其次，建立節(jié)點(diǎn)離散方程； 最后，代數(shù)方程組的形成。最后，代數(shù)方程組的形成。 對(duì)節(jié)點(diǎn)對(duì)節(jié)點(diǎn) (m,n) (m,n) 的代數(shù)方程，當(dāng)?shù)拇鷶?shù)方程，當(dāng) x=x=y y 時(shí)，時(shí)

9、，有：有： ,1,1,1,11()4m nmnmnm nm nttttt（4 4） 設(shè)立迭代初場設(shè)立迭代初場 代數(shù)方程組的求解方法有直接解法與迭代數(shù)方程組的求解方法有直接解法與迭代解法，傳熱問題的有限差分法中主要采用代解法，傳熱問題的有限差分法中主要采用迭代法。采用迭代法求解時(shí)，需對(duì)被求的溫迭代法。采用迭代法求解時(shí)，需對(duì)被求的溫度場預(yù)先設(shè)定一個(gè)解，這個(gè)解稱為度場預(yù)先設(shè)定一個(gè)解，這個(gè)解稱為初場初場，并，并在求解過程中不斷改進(jìn)。在求解過程中不斷改進(jìn)。 （5 5）求解代數(shù)方程組）求解代數(shù)方程組 本例中除本例中除 m=1 m=1 的左邊界上的左邊界上各節(jié)點(diǎn)的溫度已知外，其余各節(jié)點(diǎn)的溫度已知外，其余 (

10、M-1)N (M-1)N 個(gè)節(jié)點(diǎn)均需建立離散個(gè)節(jié)點(diǎn)均需建立離散方程，共有方程，共有 (M-1)N (M-1)N 個(gè)方程，個(gè)方程，則構(gòu)成一個(gè)封閉的代數(shù)方程則構(gòu)成一個(gè)封閉的代數(shù)方程組。組。xyxynm(m,n)MN求解時(shí)遇到的問題：求解時(shí)遇到的問題： 線性；線性； 非線性；非線性； 收斂性等。收斂性等。 2 2 ）非線性代數(shù)方程組：）非線性代數(shù)方程組：代數(shù)方程一經(jīng)建立，其中代數(shù)方程一經(jīng)建立，其中各項(xiàng)系數(shù)在整個(gè)求解過程中不斷更新。各項(xiàng)系數(shù)在整個(gè)求解過程中不斷更新。 3 3 ）是否收斂判斷：）是否收斂判斷：是指用迭代法求解代數(shù)方程是是指用迭代法求解代數(shù)方程是否收斂，即本次迭代計(jì)算所得之解與上一次迭代計(jì)

11、否收斂，即本次迭代計(jì)算所得之解與上一次迭代計(jì)算所得之解的偏差是否小于允許值。算所得之解的偏差是否小于允許值。 1 1 ）線性代數(shù)方程組：）線性代數(shù)方程組：代數(shù)方程一經(jīng)建立，其中代數(shù)方程一經(jīng)建立，其中各項(xiàng)系數(shù)在整個(gè)求解過程中不再變化；各項(xiàng)系數(shù)在整個(gè)求解過程中不再變化；（6 6） 解的分析解的分析 通過求解代數(shù)方程，獲得物體中的溫度分布，通過求解代數(shù)方程，獲得物體中的溫度分布，根據(jù)溫度場應(yīng)進(jìn)一步計(jì)算通過的熱流量，熱根據(jù)溫度場應(yīng)進(jìn)一步計(jì)算通過的熱流量，熱應(yīng)力及熱變形等。應(yīng)力及熱變形等。因此，對(duì)于數(shù)值分析計(jì)算所得的溫度場及其因此，對(duì)于數(shù)值分析計(jì)算所得的溫度場及其它物理量應(yīng)作詳細(xì)分析，以獲得定性或定量它

12、物理量應(yīng)作詳細(xì)分析，以獲得定性或定量上的結(jié)論。上的結(jié)論。 4.2 4.2 內(nèi)節(jié)點(diǎn)離散方程的建立方法內(nèi)節(jié)點(diǎn)離散方程的建立方法(1) Taylor(1) Taylor（泰勒）級(jí)數(shù)展開法；（泰勒）級(jí)數(shù)展開法；(2) (2) 控制容積平衡法控制容積平衡法( (熱平衡法熱平衡法) )4.2.1 4.2.1 泰勒級(jí)數(shù)展開法泰勒級(jí)數(shù)展開法( )2( )( )( )()( )1!2!nnfxfxfxf xxf xxxxn 2233441,234,2624mnmnmnmnmntxtxtxtttxxxxx用節(jié)點(diǎn)用節(jié)點(diǎn)(m,n)(m,n)的溫度的溫度t tm,nm,n來表示節(jié)點(diǎn)來表示節(jié)點(diǎn)(m-1,n)(m-1,n)

13、的的溫度溫度t tm-1,nm-1,n根據(jù)泰勒級(jí)數(shù)展開式，用節(jié)點(diǎn)根據(jù)泰勒級(jí)數(shù)展開式，用節(jié)點(diǎn)( (m,nm,n) )的溫度的溫度t tm,nm,n來表示節(jié)點(diǎn)來表示節(jié)點(diǎn)( (m+1,nm+1,n) )的溫度的溫度t tm+1,nm+1,n2233441,234,2624mnmnmnmnmntxtxtxtttxxxxx將上兩式相加可得將上兩式相加可得24421,1,24,212mnmnm nm ntxttttxxx22,mntx將上式改寫成將上式改寫成 的表達(dá)式，有的表達(dá)式，有)(222, 1, 1,22xoxtttxtnmnmnmnm)(2221,1,22yoytttytnmnmnmnm同樣可得：

14、同樣可得：表示未明確寫出的表示未明確寫出的級(jí)數(shù)余項(xiàng)中的級(jí)數(shù)余項(xiàng)中的XX的最低階數(shù)為的最低階數(shù)為2 2 根據(jù)導(dǎo)熱問題的控制方程根據(jù)導(dǎo)熱問題的控制方程 ( ( 導(dǎo)熱微分方程導(dǎo)熱微分方程 ) )1,1,1,122220mnmnmnmnmnmnttttttxy若若 x=x=y y 則有則有 ,1,1,1,11()4m nmnmnm nm nttttt22220ttxy得得一階一階基本思想：基本思想：對(duì)每個(gè)有限大小的控制容積應(yīng)用能量守對(duì)每個(gè)有限大小的控制容積應(yīng)用能量守恒，從而獲得溫度場的代數(shù)方程組，它從基本物理恒，從而獲得溫度場的代數(shù)方程組，它從基本物理現(xiàn)象和基本定律出發(fā)，不必事先建立控制方程，依現(xiàn)象和

15、基本定律出發(fā)，不必事先建立控制方程，依據(jù)能量守恒和據(jù)能量守恒和FourierFourier導(dǎo)熱定律即可。導(dǎo)熱定律即可。能量守恒：能量守恒：流入控制體的總熱流量控制體內(nèi)熱源生成熱流入控制體的總熱流量控制體內(nèi)熱源生成熱 流出控制體的總熱流量控制體內(nèi)能的增量流出控制體的總熱流量控制體內(nèi)能的增量ovi4.2.2 4.2.2 控制容積平衡法控制容積平衡法( (熱平衡法熱平衡法) )voi)(ovi 從所有方向流入控制體的凈熱流量從所有方向流入控制體的凈熱流量 控制體內(nèi)熱源生成熱控制體內(nèi)能的增量控制體內(nèi)熱源生成熱控制體內(nèi)能的增量注意：上面的公式對(duì)內(nèi)部節(jié)點(diǎn)和邊界節(jié)點(diǎn)均適用注意：上面的公式對(duì)內(nèi)部節(jié)點(diǎn)和邊界節(jié)點(diǎn)

16、均適用穩(wěn)態(tài)、無內(nèi)熱源時(shí)：穩(wěn)態(tài)、無內(nèi)熱源時(shí)：從所有方向流入控制體的總熱流量從所有方向流入控制體的總熱流量0 01,mnm nettyx ,1,m nm nnttxy ,1,m nm nsttxy 從節(jié)點(diǎn)通過界面?zhèn)鲗?dǎo)到節(jié)點(diǎn)從節(jié)點(diǎn)通過界面?zhèn)鲗?dǎo)到節(jié)點(diǎn) (m,n) (m,n) 的熱流量：的熱流量：1,mnm nwttyx 對(duì)元體對(duì)元體 (m,n). (m,n). 根據(jù)能量守恒定律可知：根據(jù)能量守恒定律可知： 0ewns 1,m nmnwttyx 1,m nmnettyx , 1,mnmnnttxy , 1,mnmnsttxy +=0穩(wěn)態(tài)、無內(nèi)熱源時(shí)：穩(wěn)態(tài)、無內(nèi)熱源時(shí)：從所有方向流入控制體的總熱流量從所有

17、方向流入控制體的總熱流量0 01,1,1,12222 0mnm nmnm nm nm nttttttxy化簡得說明：說明： 上述分析與推導(dǎo)在笛卡兒坐標(biāo)系中進(jìn)行的；上述分析與推導(dǎo)在笛卡兒坐標(biāo)系中進(jìn)行的； 熱平衡法概念清晰，過程簡捷；熱平衡法概念清晰，過程簡捷； 熱平衡法與建立微分方程的思路與過程一致，熱平衡法與建立微分方程的思路與過程一致，但不同的是前者是有限大小的元體，后者是微但不同的是前者是有限大小的元體，后者是微元體。元體。 4.3 4.3 邊界節(jié)點(diǎn)離散方程的建立邊界節(jié)點(diǎn)離散方程的建立及代數(shù)方程的求解及代數(shù)方程的求解 對(duì)于對(duì)于第一類邊界條件第一類邊界條件的熱傳導(dǎo)問題，處理比較簡單，因的熱傳

18、導(dǎo)問題，處理比較簡單，因?yàn)橐阎吔绲臏囟?，可將其以?shù)值的形式加入到內(nèi)節(jié)點(diǎn)的離為已知邊界的溫度，可將其以數(shù)值的形式加入到內(nèi)節(jié)點(diǎn)的離散方程中，組成封閉的代數(shù)方程組，直接求解。散方程中，組成封閉的代數(shù)方程組，直接求解。 對(duì)于對(duì)于第二類或第三類邊界條件第二類或第三類邊界條件的導(dǎo)熱問題，所有內(nèi)節(jié)點(diǎn)的的導(dǎo)熱問題，所有內(nèi)節(jié)點(diǎn)的離散方程組成的代數(shù)方程組是不封閉的，因未知邊界溫度，離散方程組成的代數(shù)方程組是不封閉的，因未知邊界溫度，因此應(yīng)對(duì)邊界上的節(jié)點(diǎn)補(bǔ)充相應(yīng)的代數(shù)方程，才能使方程組因此應(yīng)對(duì)邊界上的節(jié)點(diǎn)補(bǔ)充相應(yīng)的代數(shù)方程，才能使方程組封閉，以便求解。封閉，以便求解。 為了求解方便，將第二類邊界條件及第三類邊界條

19、件合并為了求解方便，將第二類邊界條件及第三類邊界條件合并起來考慮，用起來考慮，用q qww表示邊界上的熱流密度或熱流密度表達(dá)式。表示邊界上的熱流密度或熱流密度表達(dá)式。為使結(jié)果更具一般性，假設(shè)物體具有內(nèi)熱源為使結(jié)果更具一般性，假設(shè)物體具有內(nèi)熱源( ( 不必均勻分不必均勻分布布 ) ) 。xyqw邊界節(jié)點(diǎn)邊界節(jié)點(diǎn) (m,n) (m,n) 只代表半個(gè)元體，若邊界上有向該元只代表半個(gè)元體，若邊界上有向該元體傳遞的熱流密度為體傳遞的熱流密度為q qww ，據(jù)能量守恒定律：，據(jù)能量守恒定律： 4.3.1 4.3.1 邊界節(jié)點(diǎn)離散方程的建立邊界節(jié)點(diǎn)離散方程的建立( (1) 1) 平直邊界上的節(jié)點(diǎn)平直邊界上的

20、節(jié)點(diǎn)1,1,1,2022mnmnmnmnmnmnmnwttttxyxyttxxyyqy yx2,1,1,12124m nwm nmnm nm nx xqtttt(2) (2) 外部角點(diǎn)外部角點(diǎn)2,1,12122m nwm nmnm nx xqttt1,1,22042mnm nm nm nm nwttttyxxyx yxyq yx如圖所示，二維墻角計(jì)算區(qū)域中，該節(jié)點(diǎn)外角點(diǎn)僅如圖所示，二維墻角計(jì)算區(qū)域中，該節(jié)點(diǎn)外角點(diǎn)僅代表代表 1/4 1/4 個(gè)以個(gè)以 為邊長的元體。假設(shè)邊界上有為邊長的元體。假設(shè)邊界上有向該元體傳遞的熱流密度為向該元體傳遞的熱流密度為 ，則據(jù)能量守恒定律，則據(jù)能量守恒定律得其熱平

21、衡式為：得其熱平衡式為： x、 ywq(3) (3) 內(nèi)部角點(diǎn)內(nèi)部角點(diǎn)22,1,1,11,213(22)62wm nmnm nm nmnx qxttttt1,1,1,1,230242mnm nm nm nm nm nmnm nm nwttttttxyxxyyttyx yxyqx yx內(nèi)部角點(diǎn)代表了內(nèi)部角點(diǎn)代表了 3/4 3/4 個(gè)元體，在同樣的假設(shè)條件下個(gè)元體，在同樣的假設(shè)條件下xyqw討論關(guān)于邊界熱流密度的三種情況：討論關(guān)于邊界熱流密度的三種情況： （1 1）絕熱邊界）絕熱邊界即令上式即令上式 即可。即可。 0wq （2 2） 值不為零值不為零wq（3 3）對(duì)流邊界）對(duì)流邊界此時(shí)此時(shí) ，將此

22、表達(dá)式代入上述方程，并，將此表達(dá)式代入上述方程，并將此項(xiàng)中的將此項(xiàng)中的 與等號(hào)前的與等號(hào)前的 合并。合并。對(duì)于對(duì)于 的情形有：的情形有：)(,nmfwtthq,m nt,m ntxy wq流入元體，流入元體， 取正，流出元體，取正，流出元體， 取負(fù)取負(fù)wq（a a）平直邊界）平直邊界（b b）外部角點(diǎn)）外部角點(diǎn)（c c）內(nèi)部角點(diǎn)）內(nèi)部角點(diǎn)2,1,1,12222m nm nmnm nm nfh xxh xttttt2,1,12212m nm nmnm nfh xxh xtttt2,1,11,1322322m nm nmnm nmnm nfh xxh xtttttt4.3.2 4.3.2 處理不規(guī)

23、則區(qū)域的階梯型逼進(jìn)法處理不規(guī)則區(qū)域的階梯型逼進(jìn)法當(dāng)計(jì)算區(qū)域出現(xiàn)曲線邊界或傾斜邊界時(shí)，常常采當(dāng)計(jì)算區(qū)域出現(xiàn)曲線邊界或傾斜邊界時(shí)，常常采用用階梯形的折線階梯形的折線來模擬真實(shí)邊界，然后用上述方來模擬真實(shí)邊界，然后用上述方法建立邊界節(jié)點(diǎn)的離散方程。法建立邊界節(jié)點(diǎn)的離散方程。4.3.3 4.3.3 代數(shù)方程的求解方法代數(shù)方程的求解方法 2 2）迭代法：）迭代法：先對(duì)要計(jì)算的場作出假設(shè)（設(shè)先對(duì)要計(jì)算的場作出假設(shè)（設(shè)定初場），在迭代計(jì)算中不斷予以改進(jìn)，直定初場），在迭代計(jì)算中不斷予以改進(jìn)，直到計(jì)算前的假定值與計(jì)算結(jié)果相差小于允許到計(jì)算前的假定值與計(jì)算結(jié)果相差小于允許值為止的方法，稱迭代計(jì)算收斂。值為止的方

24、法，稱迭代計(jì)算收斂。1 1）直接解法：）直接解法：通過有限次運(yùn)算獲得精確解通過有限次運(yùn)算獲得精確解的方法，如：矩陣求解，高斯消元法。的方法，如：矩陣求解，高斯消元法。 2 2 迭代法目前應(yīng)用較多的是：迭代法目前應(yīng)用較多的是： 1 1 ）雅可比迭代法（簡單迭代）：）雅可比迭代法（簡單迭代）：每次迭代每次迭代計(jì)算，均用上一次迭代計(jì)算出的值。計(jì)算，均用上一次迭代計(jì)算出的值。 2 2 ）高斯）高斯賽德爾迭代法：賽德爾迭代法：每次迭代計(jì)算，每次迭代計(jì)算，均是使用節(jié)點(diǎn)溫度的最新值。均是使用節(jié)點(diǎn)溫度的最新值。 在計(jì)算后面的節(jié)點(diǎn)溫度時(shí)應(yīng)按下式（采用最新值）在計(jì)算后面的節(jié)點(diǎn)溫度時(shí)應(yīng)按下式（采用最新值）例如：根據(jù)

25、第例如：根據(jù)第 k k 次迭代的數(shù)值次迭代的數(shù)值(k)n(k)2(k)1.ttt、可以求得節(jié)點(diǎn)溫度：可以求得節(jié)點(diǎn)溫度：)(1)(1)(212)(111) 1(1.kknnkkkbtatatat)()() 1(11) 1(22) 1(11) 1()(3)(3) 1(232) 1(131) 1(3)(2)(2)(222) 1(121) 1(2.knknnnknnnknknknkknnkkkkknnkkkbtatatatatbtatatatbtatatat設(shè)有一三元方程組設(shè)有一三元方程組： 11112 213 3121122 223 3231132 233 33a ta ta tba ta ta tba ta ta tb其中其中 （ i=1,2,3 i=1,2,3 ； j=1,2,3 j=1,2,3 ）及）及 是已是已知的系數(shù)（均不為零）及常數(shù)。知的系數(shù)（均不為零）及常數(shù)。, i jaib采用高斯采用高斯賽德爾迭代法的步驟：賽德爾迭代法的步驟： （1）將三元方程變形為迭式方程：）將三元方程變形為迭式方程： 1112 213 3112221 123 3223331 132 2331()1()1()tba ta tatb