河海大學(xué)土木工程學(xué)院力學(xué)系　
   我國(guó)幅員遼闊，江海湖泊繁多，水利資源豐富，同時(shí)也是爆發(fā)洪澇災(zāi)害的根源。 如何充分利用水利資源、抗洪減災(zāi)、造福人類(lèi)是水利工作者首要的任務(wù)。近年來(lái)在黨中央和各地政府的關(guān)心和支持下， 各地紛紛都在興修水利，如水電站、攔洪大壩，其中三峽工程的建設(shè)是一個(gè)最典型的例子。一座大壩的建成，不是一兩個(gè)月就能完成， 而是好幾個(gè)月乃至一兩年的浩大工程，因此經(jīng)濟(jì)效益和安全可靠是兩個(gè)十分重要的問(wèn)題。如何保證大壩的安全可靠，一方面， 工程質(zhì)量是一個(gè)重要因素，另一方面，大壩本身結(jié)構(gòu)強(qiáng)度的可靠性是關(guān)鍵。如何保證結(jié)構(gòu)物的強(qiáng)度，一是需要有一個(gè)合理的結(jié)構(gòu)形狀， 另外要合理地安排施工期。對(duì)于這樣龐大的建筑物，采用試驗(yàn)?zāi)M是不現(xiàn)實(shí)的，而采用常規(guī)的分析辦法也不能得到滿意的結(jié)果。現(xiàn)在普遍采用的方法是數(shù)值模擬技術(shù)即計(jì)算機(jī)仿真，其中以有限元分析為最普遍。

    我們對(duì)龍首電站大壩整個(gè)施工過(guò)程的溫度場(chǎng)與溫度應(yīng)力的仿真分析，包括擋水、泄水、引水系統(tǒng)和水電站廠房等建筑物。攔河大壩由碾壓砼雙曲拱壩為主體及左岸重力壩和右岸推力墩組成，壩身中有泄水中孔及表孔、左岸沖沙孔和電站引水系統(tǒng)等，形體和結(jié)構(gòu)部極為復(fù)雜，壩址區(qū)河谷呈V型，邊坡陡峻，兩岸山體高出河面300m以上。壩區(qū)地質(zhì)條件良好。該地區(qū)屬大陸性氣候，夏季酷熱、雨量稀少、蒸發(fā)強(qiáng)盛， 冬季嚴(yán)寒、冰期長(zhǎng)達(dá)四個(gè)月，年平均氣溫為8.5℃，多年最高氣溫為37.2℃、最低氣溫為-33.0℃，日溫差較大。 由于大壩高度不是太大，但壩區(qū)氣溫變幅較大，故溫度荷載將是大壩最主要的荷載， 因而對(duì)大壩的溫度和溫度應(yīng)力的分析計(jì)算應(yīng)予重視，對(duì)大壩施工期的溫度、溫度應(yīng)力、總體應(yīng)力進(jìn)行了仿真計(jì)算分析。
    計(jì)算模型如圖所示，共計(jì)11,408個(gè)單元，節(jié)點(diǎn)數(shù)為14,550。采用ANSYS的熱分析和結(jié)構(gòu)分析模塊， 先由熱分析計(jì)算出整個(gè)施工期的溫度場(chǎng)分布，然后將各個(gè)時(shí)刻的溫度場(chǎng)作為結(jié)構(gòu)熱應(yīng)力分析的溫度條件計(jì)算出整個(gè)施工期的結(jié)構(gòu)應(yīng)力分布。 采用ANSYS軟件的單元“死活”功能來(lái)模擬大壩的混凝土澆注過(guò)程，即把整個(gè)模型建好并進(jìn)行全部網(wǎng)格剖分，在計(jì)算過(guò)程中， 可以按要求逐步激活模型中的某些單元，以達(dá)到仿真整個(gè)結(jié)構(gòu)分層、分塊施工過(guò)程的效果。在對(duì)大壩施工期溫度和溫度應(yīng)力仿真計(jì)算中， 我們進(jìn)行了3年多（1160天）逐日變化的溫度場(chǎng)和應(yīng)力場(chǎng)仿真計(jì)算。