板帶熱軋過程有限元快速求解方法及應(yīng)用.pdf

1、本文以國家自然科學(xué)基金和寶鋼聯(lián)合資助項(xiàng)目“板材軋制過程中有限元高速在線算法基礎(chǔ)”和韓國浦項(xiàng)國際合作項(xiàng)目“熱軋過程溫度高精度預(yù)測(cè)模型開發(fā)”為研究背景，通過理論分析有限元求解板材軋制過程計(jì)算時(shí)間和收斂穩(wěn)定性的各種影響因素，結(jié)合實(shí)際軋制過程，研究了加快有限元求解速度的各種方法。通過離線分析和在線調(diào)試，將有限元快速求解程序應(yīng)用于天津榮程750mm板帶精軋過程控制系統(tǒng)，實(shí)現(xiàn)了板帶熱軋過程有限元在線應(yīng)用。本文主要研究工作如下：
　　 (1)

2、基于剛塑性可壓縮有限元法，研究了能夠提高能量泛函和最優(yōu)化過程快速求解的措施，包括程序編制、矩陣存儲(chǔ)、線性方程組求解和阻尼因子搜索等。利用Compaq Visual Fortran語言開發(fā)了軋制參數(shù)計(jì)算程序FFE-RF2D，利用FFE-RF2D求解了某鋼廠板帶精軋過程軋制參數(shù)，并對(duì)剛塑性有限元求解過程各個(gè)模塊計(jì)算時(shí)間進(jìn)行了分析。此外，研究了編譯方式、單元?jiǎng)澐?、初始速度?chǎng)設(shè)定和相關(guān)參數(shù)對(duì)計(jì)算速度影響。
　　 (2)研究了Hessia

3、n矩陣強(qiáng)迫正定Cholesky分解法，從而避免了不正定或弱正定的Hessian矩陣降低求解速度和收斂穩(wěn)定性?；贖essian矩陣一維變帶寬壓縮存儲(chǔ)和線性方程組Cholesky求解方法，提出板材軋制過程有限元求解的NR-C算法(NR-Cholesky)，使單道次軋制過程計(jì)算時(shí)間少于200ms，提高了計(jì)算效率。在研究黃金分割法和Brent精確線性搜索方法基礎(chǔ)上，提出了NR-C-B算法(NR-Cholesky-Brent)，與NR法相比，N

4、R-C－B算法迭代步數(shù)明顯減少，單道次軋制過程計(jì)算時(shí)間少于150ms，進(jìn)一步提高了計(jì)算效率和收斂穩(wěn)定性。
　　 (3)在研究信賴域和擬牛頓方法基礎(chǔ)上，提出了混合信賴域和牛頓迭代的TR-NR算法以及NR混合BFGS的NR-BFGS算法。TR-NR相比NR算法迭代步數(shù)變化范圍較小，約為30—40之間，迭代求解過程較為穩(wěn)定。NR-BFGS算法不僅迭代步數(shù)明顯減少，而且平均每迭代步計(jì)算時(shí)間較短，在提高迭代求解速度同時(shí)，進(jìn)一步提高了收斂穩(wěn)

5、定性。利用各種算法對(duì)相同實(shí)際軋制過程進(jìn)行了求解與測(cè)試，迭代步數(shù)從多到少依次為NR、TR-NR、NR-C、NR-C-B和NR-BFGS，迭代求解穩(wěn)定性依次為NR、NR-C、NR-C-B、TR-NR和NR-BFGS; NR-BFGS算法計(jì)算時(shí)間最短，單元?jiǎng)澐譃?6×8條件下，單道次軋制過程迭代步數(shù)少于25，計(jì)算時(shí)間少于lOOms。
　　 (4)通過開發(fā)的FFE-TEMP2D程序分析了單元?jiǎng)澐趾蜁r(shí)間步長選取對(duì)溫度求解和計(jì)算時(shí)間的影響。

6、在研究集中熱容矩陣法基礎(chǔ)上，提出逐層細(xì)分網(wǎng)格法抑制溫度振蕩，提高計(jì)算精度，取得了較好效果。與均勻細(xì)分網(wǎng)格法相比，逐層細(xì)分網(wǎng)格法也提高了計(jì)算速度。利用開發(fā)程序求解了某鋼廠常規(guī)流程熱軋過程溫度變化，計(jì)算結(jié)果與實(shí)測(cè)值吻合良好，單元數(shù)目為200時(shí)，求解整個(gè)軋制過程計(jì)算時(shí)間約為700ms，計(jì)算效率較高。
　　 (5)提出了新的電磁感應(yīng)熱生成率模型作為內(nèi)熱源進(jìn)行溫度場(chǎng)有限元求解，利用ANSYS求解數(shù)據(jù)和實(shí)測(cè)值，在分析工作頻率、源電流密度、空

7、氣間隙和節(jié)點(diǎn)與表面距離的基礎(chǔ)上，建立了熱生成率模型。根據(jù)新的熱生成率模型對(duì)某鋼廠短流程熱軋過程進(jìn)行了溫度求解，計(jì)算所得溫度和實(shí)測(cè)值吻合良好。建立了軋制過程溫度、變形和相變耦合過程有限元求解模型，根據(jù)試驗(yàn)所得熱物性參數(shù)和W20變形抗力模型對(duì)某鋼廠1700mm無取向硅鋼W20粗軋過程某一道次軋制過程進(jìn)行了變形求解，獲得了溫度、應(yīng)力和相分布，為翹頭分析提供了理論依據(jù)。
　　 (6)在研究實(shí)際軋制過程控制系統(tǒng)基礎(chǔ)上，建立了適合FEM調(diào)用

8、的數(shù)據(jù)庫，主要包括熱物理性能模塊、層別數(shù)據(jù)模塊、系統(tǒng)常數(shù)模塊和其它控制參數(shù)模塊。以軋制參數(shù)求解的數(shù)學(xué)模型和有限元快速求解模型為基礎(chǔ)，開發(fā)了板帶熱軋過程離線模擬系統(tǒng)，為有限元在線應(yīng)用提供了軋制過程設(shè)定與計(jì)算、算法測(cè)試、結(jié)果比較和模型改進(jìn)的平臺(tái)。
　　 (7)以榮鋼750mm精軋過程控制系統(tǒng)為背景，研究了模型設(shè)定計(jì)算流程、數(shù)據(jù)傳輸和通訊方法，將軋制力和溫度有限元快速求解模型通過動(dòng)態(tài)鏈接庫形式嵌入到精軋過程控制程序中，實(shí)現(xiàn)了軋制過程有