閆 濱，王 闖

(沈陽(yáng)農(nóng)業(yè)大學(xué)水利學(xué)院，遼寧沈陽(yáng)110866)

基于子母模型聯(lián)合反饋修正算法的閘墩溫控防裂仿真分析

閆 濱，王 闖

(沈陽(yáng)農(nóng)業(yè)大學(xué)水利學(xué)院，遼寧沈陽(yáng)110866)

應(yīng)用ANSYS軟件，在運(yùn)用APDL語(yǔ)言編制仿真程序的過(guò)程中將子母模型聯(lián)合反饋修正算法應(yīng)用于計(jì)算閘墩混凝土隨時(shí)間變化的溫度和應(yīng)力，計(jì)算中在時(shí)間和空間上進(jìn)行加密，計(jì)算一次子模型，就將結(jié)果反饋至母模型。采用該算法得到的溫度場(chǎng)及應(yīng)力場(chǎng)計(jì)算結(jié)果的精度均有所提高。在雙臺(tái)子河閘工程實(shí)例分析中應(yīng)用此算法模擬出閘墩拆模后第105 d時(shí)受到寒潮作用導(dǎo)致閘墩表面開(kāi)裂，與實(shí)際情況相符合，證明了在應(yīng)用ANSYS進(jìn)行仿真計(jì)算時(shí)運(yùn)用該算法可以較為準(zhǔn)確的模擬寒潮作用下閘墩的溫度場(chǎng)及應(yīng)力場(chǎng)。

子母模型聯(lián)合反饋修正算法；裂縫；閘墩；溫度應(yīng)力；ANSYS；仿真分析

0 引 言

大體積混凝土閘墩外部熱量傳遞速度比內(nèi)部快，內(nèi)外溫差存在較大差異，從而使得混凝土內(nèi)外變形不一致。當(dāng)混凝土的拉應(yīng)變超過(guò)其極限拉應(yīng)變時(shí)，混凝土將會(huì)開(kāi)裂，寒潮作用會(huì)使裂縫變得更加嚴(yán)重。因此做好溫控防裂工作則尤為重要。對(duì)于直接受寒潮作用影響的大體積混凝土閘墩，為了更加精確地仿真計(jì)算其表層混凝土的溫度和應(yīng)力，就需要精細(xì)劃分其表層網(wǎng)格，這不僅使得劃分網(wǎng)格的難度增大，而且還將成倍的增加整體網(wǎng)格及單元數(shù)量，致使計(jì)算機(jī)計(jì)算時(shí)間冗長(zhǎng)。為此，本文應(yīng)用ANSYS軟件將子母模型聯(lián)合反饋修正算法[1]應(yīng)用于混凝土隨時(shí)間變化的溫度和應(yīng)力的求解，在不刻意加細(xì)劃分網(wǎng)格的情況下，提高仿真計(jì)算精度。

1 子母模型聯(lián)合反饋修正算法

子母模型聯(lián)合反饋修正算法是有限元子模型法的優(yōu)化算法，從時(shí)間τn到τn+1母模型計(jì)算的精確解為Fn+1，將下一步計(jì)算的近似解記為Fn→n+1，則在子模型計(jì)算時(shí)，就需要加密時(shí)間步，即離散時(shí)間區(qū)域(τn，τn+1)為若干個(gè)時(shí)間子步(τn，τn1，τn2，…，τn+1)。在子步的時(shí)間范圍內(nèi)，在計(jì)算每一步子模型的時(shí)候都要反饋并修正母模型，此時(shí)不加密母模型，之后再將下一步子模型的計(jì)算反饋修正。從圖1可見(jiàn)，從τn到τn1的計(jì)算結(jié)果為Fn1，從τn1到τn+1的計(jì)算結(jié)果為Fn1→n+1，之后反饋給下一步子模型的計(jì)算，即從τn1到τn2的計(jì)算步，按照此方式繼續(xù)計(jì)算，這樣使得母模型在子域范圍內(nèi)的誤差減小，從而使子模型計(jì)算精度得到提高。在空間范圍內(nèi)，經(jīng)過(guò)加密的子模型比母模型計(jì)算精度要高，因此計(jì)算一次子模型，就反饋給母模型一次結(jié)果[1]。

圖1 時(shí)間步長(zhǎng)誤差修正

子母模型聯(lián)合反饋修正算法的示意見(jiàn)圖2。從溫度場(chǎng)來(lái)看，該算法可將每一步子模型更高精度的溫度場(chǎng)計(jì)算結(jié)果反饋給母模型，從而提高在子域內(nèi)母模型下一步計(jì)算的溫度場(chǎng)的精度。從應(yīng)力場(chǎng)來(lái)看，計(jì)算溫度應(yīng)力時(shí)利用反饋修正后的溫度增量的計(jì)算結(jié)果比時(shí)間步和網(wǎng)格都沒(méi)有加密的母模型的計(jì)算結(jié)果精度要高。與此同時(shí)，該算法中對(duì)溫度場(chǎng)或應(yīng)力場(chǎng)時(shí)間步的計(jì)算是否加密可以根據(jù)實(shí)際需要來(lái)決定，例如早期澆筑的閘墩混凝土，由于水化熱作用使得閘墩溫度變化劇烈，那么需要考慮加密時(shí)間步，當(dāng)進(jìn)入準(zhǔn)穩(wěn)定期時(shí)，則不需要對(duì)時(shí)間步進(jìn)行加密[1]。此算法可以極大的提高仿真計(jì)算結(jié)果的精度。

圖2 子母模型聯(lián)合反饋修正算法

2 實(shí)例分析

2.1 工程概況

盤錦雙臺(tái)子河閘除險(xiǎn)加固工程新建18孔淺孔閘，共有19個(gè)閘墩，其中有9個(gè)中墩，8個(gè)縫墩，2個(gè)邊墩。閘墩長(zhǎng)13.98m，高8.08m，中墩寬2m，縫墩寬3m。閘墩由C30混凝土澆筑而成，抗?jié)B標(biāo)號(hào)為W6，抗凍標(biāo)號(hào)為F200，三級(jí)配骨料。閘底板澆筑完成2個(gè)多月后開(kāi)始澆筑閘墩，2014年9月25日開(kāi)始澆筑閘墩，2014年11月13日澆筑完畢。澆筑期間的外界日最高氣溫為23～6 ℃，外界日最低氣溫為16～-4 ℃。澆筑完成之后12～18h進(jìn)行流水養(yǎng)護(hù)。當(dāng)混凝土強(qiáng)度達(dá)到3.5MPa時(shí)，即澆筑完成之后的第3d天拆除側(cè)面模板，并采用土工布覆蓋保濕。越冬期間混凝土未采取保溫措施。拆模后的第104～106d陸續(xù)在各個(gè)閘墩上出現(xiàn)裂縫。僅8號(hào)閘墩未見(jiàn)明顯裂縫。本文以10號(hào)中墩為例驗(yàn)證子母模型聯(lián)合反饋修正算法應(yīng)用在閘墩溫控防裂仿真分析中的合理性。

2.2 閘墩開(kāi)裂原因分析

本工程閘室底板地基采用振沖碎石樁進(jìn)行加固處理，設(shè)計(jì)要求經(jīng)處理后的復(fù)合地基承載力達(dá)到200kPa。設(shè)計(jì)樁徑1m，間距1.8m。施工時(shí)段為3月29日～5月2日，平均處理深度17.5m，共計(jì)1 317根樁。施工過(guò)程中各項(xiàng)參數(shù)嚴(yán)格按照規(guī)范上限執(zhí)行。按照《水利水電工程振沖法地基處理技術(shù)規(guī)范》要求，施工結(jié)束并恢復(fù)期達(dá)到28d后對(duì)已完工程進(jìn)行靜載試驗(yàn)，共檢測(cè)5組，檢測(cè)點(diǎn)數(shù)滿足規(guī)范要求，檢測(cè)結(jié)果全部滿足設(shè)計(jì)200kPa要求。混凝土閘墩坐落于閘底板上，因此閘墩不存在地基不均勻沉陷問(wèn)題，也不會(huì)產(chǎn)生不均勻沉陷裂縫。

大體積混凝土閘墩外部熱量傳遞速度比內(nèi)部快，內(nèi)外溫差存在較大差異，從而使得混凝土內(nèi)外變形不一致。當(dāng)混凝土的拉應(yīng)變超過(guò)其極限拉應(yīng)變時(shí)，混凝土將會(huì)開(kāi)裂，寒潮作用會(huì)使裂縫變得更加嚴(yán)重?？紤]在越冬期間混凝土未采取保溫措施，因此，分析認(rèn)為溫度應(yīng)力是導(dǎo)致閘墩開(kāi)裂的主要原因之一。

根據(jù)雙臺(tái)子河閘除險(xiǎn)加固工程閘墩裂縫檢測(cè)報(bào)告，37條裂縫均未到達(dá)墩頂，其中有24條裂縫未到達(dá)墩底，在距離底板以上0.6～4m范圍內(nèi)是裂縫寬度平均最大值段，裂縫長(zhǎng)度在1.72～5.68m范圍，主要分布在偏下位置，上述特征符合由于底板約束作用產(chǎn)生裂縫的特征，因此，底板約束是閘墩開(kāi)裂的另一個(gè)主要原因。

綜上所述，閘墩裂縫形成的主要原因是溫度應(yīng)力和底板約束。

2.3 閘墩溫度場(chǎng)和應(yīng)力場(chǎng)仿真模型建立

2.3.1 三維有限元計(jì)算模型

本文應(yīng)用ANSYS軟件，在運(yùn)用APDL語(yǔ)言編制仿真程序的過(guò)程中將子母模型聯(lián)合反饋修正算法應(yīng)用于求解混凝土隨時(shí)間變化的溫度場(chǎng)和應(yīng)力場(chǎng)[2-7]，由于閘墩拆模后128d之內(nèi)的資料較為齊全且包括閘墩出現(xiàn)裂縫的時(shí)間點(diǎn)，因此以10號(hào)中墩為例，對(duì)閘墩拆模后128d之內(nèi)的溫度場(chǎng)和應(yīng)力場(chǎng)進(jìn)行仿真計(jì)算。在2014年10月23日10號(hào)中墩開(kāi)始施工，在10月26日拆除側(cè)面模板，在2015年2月8日閘墩表面出現(xiàn)裂縫，即拆模后第105d，在此期間閘墩受寒潮作用。

10號(hào)中墩模型長(zhǎng)13.98m、高8.08m、厚2m。由于先澆筑的閘底板約2個(gè)月之后閘墩才開(kāi)始澆筑，底板可視為老混凝土，其對(duì)上部結(jié)構(gòu)的作用等同于基礎(chǔ)，因此將其看作基礎(chǔ)，同時(shí)將其簡(jiǎn)化為長(zhǎng)13.98m、寬12m、厚2.5m的長(zhǎng)方體結(jié)構(gòu)。仿真計(jì)算中考慮絕熱溫升變化、材料屬性變化、外界氣溫變化、重力荷載、邊界條件(溫度邊界條件，熱流邊界條件，對(duì)流換熱邊界條件)等影響因素。采用當(dāng)?shù)厝掌骄鶜鉁刈鳛榉抡嬗?jì)算的氣溫?cái)?shù)據(jù)，混凝土的彈性模量隨時(shí)間變化[8]。在熱分析中閘墩的四周和頂面施加對(duì)流邊界條件；底板的上表面施加對(duì)流邊界條件，底面和四周施加絕熱邊界條件。結(jié)構(gòu)分析中在閘墩下底面施加法向位移約束條件，刪除底板，這樣可以簡(jiǎn)化計(jì)算。有限元計(jì)算模型采用六面體8節(jié)點(diǎn)單元，在熱分析中采用solid70單元，結(jié)構(gòu)分析中采用solid65單元。模型的單元總數(shù)為8 835個(gè)，節(jié)點(diǎn)總數(shù)為10 878個(gè)。采用右手坐標(biāo)系建模，坐標(biāo)原點(diǎn)設(shè)定在閘墩與底板接觸矩形面的上游側(cè)靠右岸的一個(gè)頂點(diǎn)上，X軸的正方向取順?biāo)鞣较?，Y軸取橫河流方向，由右岸指向左岸，Z軸豎直向上[9]，有限元計(jì)算模型見(jiàn)圖3。

圖3 有限元計(jì)算模型

2.3.2 計(jì)算參數(shù)

(1)閘墩關(guān)鍵點(diǎn)位置。仿真所用閘墩關(guān)鍵點(diǎn)的位置如圖4、5所示。B0、C0、D0、E0、F0、G0、H0、I0、J0與A0在一條軸線上。

圖4 閘墩正立面及表面關(guān)鍵點(diǎn)位置(單位：m)

圖5 1-1截面關(guān)鍵點(diǎn)位置(單位：m)

(2)閘墩彈性模量。根據(jù)文獻(xiàn)[10]，按照式(1)計(jì)算混凝土的彈性模量，即

E(t)=βE0×(1-e-φt)

(1)

式中，E(t)為混凝土在齡期為t時(shí)的彈性模量；β為混凝土中摻合料對(duì)彈性模量的修正系數(shù)，應(yīng)以現(xiàn)場(chǎng)試驗(yàn)數(shù)據(jù)的取值為準(zhǔn)，在施工準(zhǔn)備階段和現(xiàn)場(chǎng)無(wú)試驗(yàn)數(shù)據(jù)時(shí)，可按表1計(jì)算E0，本算例中，閘墩及底板混凝土中粉煤灰和礦渣的摻量分別是26%和0，那么β1=0.98，β2=1，故β=β1·β2=0.98；E0為標(biāo)準(zhǔn)情況下混凝土養(yǎng)護(hù)28 d的彈性模量；φ為系數(shù)，應(yīng)通過(guò)混凝土試驗(yàn)確定，當(dāng)無(wú)試驗(yàn)數(shù)據(jù)時(shí)，可近似取0.09。

表1 不同摻量摻合料彈性模量調(diào)整系數(shù)

(3)混凝土絕熱溫升?；炷恋慕^熱溫升應(yīng)用在實(shí)際溫度場(chǎng)的計(jì)算中，而在ANSYS軟件中是通過(guò)生熱率Hgen來(lái)實(shí)現(xiàn)混凝土的絕熱溫升。生熱率是混凝土在單位時(shí)間內(nèi)產(chǎn)生的水化熱，水化熱的時(shí)間函數(shù)基本形式[11]為

Q(t)=Qc·Mc(1-e-mt)

(2)

式中，Q(t)為齡期t時(shí)混凝土產(chǎn)生的水化熱，kJ/m3；Qc為單位質(zhì)量水泥產(chǎn)生的水化熱，kJ/kg；Mc為單位體積混凝土的水泥用量，kg/m3；m為水泥水化速率系數(shù)，d-1；t為齡期，d。

(4)氣溫資料。采用當(dāng)?shù)厝掌骄鶜鉁刈鳛榉抡嬗?jì)算的氣溫?cái)?shù)據(jù)。

(5)混凝土熱力學(xué)參數(shù)。依據(jù)文獻(xiàn)[10-13]計(jì)算得到混凝土熱力學(xué)參數(shù)見(jiàn)表2。施工現(xiàn)場(chǎng)風(fēng)速為Va=4.5 m/s。

表2 混凝土熱力學(xué)參數(shù)

2.4 閘墩溫度場(chǎng)和應(yīng)力場(chǎng)仿真計(jì)算結(jié)果與分析

閘墩出現(xiàn)豎直方向裂縫，分析是由于閘墩在沿長(zhǎng)度方向上的拉應(yīng)力超過(guò)其抗拉強(qiáng)度產(chǎn)生；且裂縫具有明顯的規(guī)律性，基本位于閘墩中間位置距離底板20 cm左右至閘墩高度5 m左右范圍內(nèi)，因此著重研究位于閘墩中間位置沿長(zhǎng)度方向上的溫度場(chǎng)及應(yīng)力場(chǎng)。

為了驗(yàn)證采用子母模型聯(lián)合反饋修正算法是否能夠提高溫度場(chǎng)和應(yīng)力場(chǎng)的計(jì)算精度，對(duì)未采用此方法的閘墩溫度場(chǎng)及應(yīng)力場(chǎng)進(jìn)行仿真計(jì)算與此對(duì)比，從而突出采用子母模型聯(lián)合反饋修正算法進(jìn)行仿真所達(dá)到的效果。由于采用此方法仿真計(jì)算的溫度場(chǎng)與未采用此方法仿真計(jì)算的結(jié)果并無(wú)明顯差異，因此只對(duì)應(yīng)力場(chǎng)的仿真結(jié)果進(jìn)行對(duì)比分析，在溫度場(chǎng)的仿真計(jì)算結(jié)果分析部分僅分析采用此方法的溫度場(chǎng)仿真計(jì)算結(jié)果。

2.4.1 閘墩溫度場(chǎng)仿真計(jì)算

采用子母模型聯(lián)合反饋修正算法計(jì)算的閘墩拆模后128 d之內(nèi)表面關(guān)鍵點(diǎn)及截面中心處關(guān)鍵點(diǎn)溫度隨時(shí)間變化的仿真計(jì)算結(jié)果分別見(jiàn)圖6和圖7。由圖6可見(jiàn)，仿真計(jì)算的表面關(guān)鍵點(diǎn)溫度隨時(shí)間變化的走勢(shì)基本相同，選取關(guān)鍵點(diǎn)E的實(shí)測(cè)溫度數(shù)據(jù)來(lái)作對(duì)照分析。閘墩拆模后128 d之內(nèi)表面關(guān)鍵點(diǎn)E實(shí)測(cè)溫度隨時(shí)間變化曲線見(jiàn)圖8。閘墩拆模后128 d之內(nèi)當(dāng)?shù)厝掌骄鶜鉁厝鐖D9所示。

圖6 閘墩拆模后128 d之內(nèi)表面關(guān)鍵點(diǎn)溫度隨時(shí)間變化曲線

圖7 閘墩拆模后128 d之內(nèi)截面中心處關(guān)鍵點(diǎn)溫度隨時(shí)間變化曲線

圖8 閘墩拆模后128 d之內(nèi)表面關(guān)鍵點(diǎn)E實(shí)測(cè)溫度隨時(shí)間變化曲線

圖9 閘墩拆模后128 d之內(nèi)當(dāng)?shù)厝掌骄鶜鉁?/p>

由圖6、7可見(jiàn)，閘墩表層及內(nèi)部混凝土的溫度在拆模后持續(xù)下降。閘墩表層混凝土的溫度變化幅度很大，大約在拆模后第60 d左右溫度開(kāi)始逐漸平穩(wěn)在-6 ℃左右，直至3個(gè)月以后(約2015年3月)隨著外界氣溫的升高閘墩表層混凝土溫度有所回升，表層混凝土的溫度變化趨勢(shì)與圖9所示當(dāng)?shù)仄骄鶜鉁氐淖兓厔?shì)基本一致，表明閘墩表層混凝土溫度主要受外界氣溫變化的影響。仿真計(jì)算的拆模后128 d之內(nèi)的日平均氣溫變化較大，這也導(dǎo)致閘墩表層混凝土的溫度變化較大。相比而言，閘墩內(nèi)部混凝土溫度雖持續(xù)下降，但波動(dòng)幅度相對(duì)較小，大約在拆模后第70 d左右溫度開(kāi)始逐漸趨于平穩(wěn)，時(shí)間上較閘墩表層混凝土略微滯后。由此可以看出閘墩內(nèi)部混凝土溫度變化滯后于外界氣溫變化。

由圖9可見(jiàn)，閘墩拆模后60 d之內(nèi)當(dāng)?shù)厝掌骄鶜鉁刈兓容^大，之后逐漸趨于平穩(wěn)，在拆模后第53 d和第105 d混凝土溫度達(dá)到128 d之內(nèi)的最低值-11 ℃，其次是在第94 d時(shí)溫度接近-11 ℃。第103～105 d溫降最快，溫差達(dá)到11 ℃。

由圖8可見(jiàn)，關(guān)鍵點(diǎn)E實(shí)測(cè)溫度隨時(shí)間變化曲線變化較為劇烈，這是因?yàn)橥饨鐨鉁仉S時(shí)間在不斷發(fā)生變化，E點(diǎn)作為閘墩混凝土表面的關(guān)鍵點(diǎn)，其溫度直接受外界氣溫的影響；由于仿真計(jì)算的外界氣溫?cái)?shù)據(jù)采用的是日平均氣溫，由圖6所示仿真計(jì)算的溫度隨時(shí)間變化過(guò)程與實(shí)測(cè)的溫度變化過(guò)程稍有偏差，但走勢(shì)大體相同。可見(jiàn)，基于子母模型聯(lián)合反饋修正算法仿真計(jì)算的溫度場(chǎng)與實(shí)測(cè)溫度場(chǎng)較為接近，仿真計(jì)算結(jié)果符合實(shí)際情況。

2.4.2 閘墩應(yīng)力場(chǎng)仿真計(jì)算及子母模型聯(lián)合反饋修正算法運(yùn)用的對(duì)比分析

未采用和采用子母模型聯(lián)合反饋修正算法計(jì)算的閘墩拆模后128 d之內(nèi)表面關(guān)鍵點(diǎn)沿X方向的應(yīng)力隨時(shí)間變化仿真曲線分別見(jiàn)圖10、11。

圖10 閘墩拆模后128 d之內(nèi)表面關(guān)鍵點(diǎn)沿X方向的應(yīng)力隨時(shí)間變化曲線

圖11 閘墩拆模后128 d之內(nèi)表面關(guān)鍵點(diǎn)沿X方向的應(yīng)力隨時(shí)間變化曲線

由于計(jì)算中采用子母模型聯(lián)合反饋修正算法的計(jì)算模型在時(shí)間和空間上進(jìn)行加密與未采用此方法的計(jì)算模型網(wǎng)格劃分略有不同，因此關(guān)鍵點(diǎn)的選取有所差異，但因?yàn)樗嘘P(guān)鍵點(diǎn)的應(yīng)力隨時(shí)間變化的走勢(shì)均是一致的，所以關(guān)鍵點(diǎn)的選取并不影響仿真計(jì)算結(jié)果。

由圖10與圖11比較可以看出兩幅圖的走勢(shì)基本相同，但極值存在差異，采用此方法較未采用此方法極大值要略微大些，極小值略微小些，這是由于在計(jì)算每一步子模型的時(shí)候都要反饋并修正母模型，仿真計(jì)算精度與未采用此算法的仿真計(jì)算精度有所差異。從圖10可以看出，最大拉應(yīng)力出現(xiàn)在大約105 d剛過(guò)時(shí)，拉應(yīng)力大約為1.9 MPa，并未達(dá)到C30混凝土的抗拉強(qiáng)度2.01 MPa，若用未采用此算法仿真計(jì)算的結(jié)果進(jìn)行分析閘，墩在這期間并未開(kāi)裂，顯然與事實(shí)不符，這就表明未采用此算法的仿真計(jì)算效果并不是很好。

由圖11可見(jiàn)，閘墩拆模后128 d之內(nèi)表面關(guān)鍵點(diǎn)應(yīng)力數(shù)值波動(dòng)很大，主要是由于表層混凝土直接受外界氣溫劇烈變化引起的。最大拉應(yīng)力出現(xiàn)在閘墩拆模后大約第105 d剛過(guò)時(shí)，拉應(yīng)力約為2.2 MPa，比C30混凝土的抗拉強(qiáng)度2.01 MPa要大。由圖9可見(jiàn)，從第103 d到第105 d時(shí)間內(nèi)氣溫驟降11 ℃，第105 d時(shí)的日平均氣溫為-11 ℃，是仿真計(jì)算的128 d之內(nèi)溫度最低的時(shí)候，并且閘墩沒(méi)有進(jìn)行保溫，由此判斷在此時(shí)閘墩開(kāi)裂。兩個(gè)影響溫度應(yīng)力的重要因素分別是溫降和約束。當(dāng)混凝土受到寒潮作用溫降過(guò)快且在比較大的約束作用下，此時(shí)會(huì)產(chǎn)生較大的溫度應(yīng)力。由于越冬期間混凝土未采取適當(dāng)?shù)谋卮胧饨鐨鉁貙?duì)閘墩表層混凝土造成直接影響，產(chǎn)生過(guò)大的溫度應(yīng)力致使閘墩開(kāi)裂。由此可見(jiàn)，在寒潮期間對(duì)閘墩采取保溫措施至關(guān)重要，可有效減小因寒潮作用產(chǎn)生的溫度應(yīng)力對(duì)閘墩混凝土的影響。

相比未采用子母模型聯(lián)合反饋修正算法的仿真過(guò)程而言，采用此方法的仿真過(guò)程很好的模擬出閘墩的開(kāi)裂時(shí)間，因此采用此方法的仿真計(jì)算的效果與精度要優(yōu)于未采用此方法的。同樣在計(jì)算效率上，采用此方法的計(jì)算效率要高于未采用此方法的。如采用同一臺(tái)電腦(Intel Core i5-4460 處理器，3.20 GHz，四核，8G內(nèi)存)，未采用此方法的計(jì)算分析時(shí)間約為580 s，而采用此方法的計(jì)算時(shí)間約為205 s，若模型更加復(fù)雜，那么將會(huì)節(jié)省更多的時(shí)間。

5 結(jié) 論

基于子母模型聯(lián)合反饋修正算法，應(yīng)用ANSYS軟件在網(wǎng)格劃分不刻意加細(xì)的情況下可以較為準(zhǔn)確、高效的仿真計(jì)算出在寒潮作用下閘墩混凝土隨時(shí)間變化的溫度和應(yīng)力，雙臺(tái)子河閘工程實(shí)例證明了溫度場(chǎng)及應(yīng)力場(chǎng)計(jì)算中子母模型聯(lián)合反饋修正算法與ANSYS軟件結(jié)合應(yīng)用的合理性。

雙臺(tái)子河閘工程實(shí)例中，應(yīng)用ANSYS基于子母模型聯(lián)合反饋修正算法仿真計(jì)算的閘墩混凝土溫度場(chǎng)和應(yīng)力場(chǎng)結(jié)果表明，寒潮主要作用于閘墩表層混凝土，而對(duì)內(nèi)部混凝土影響較小，這造成閘墩混凝土的內(nèi)外溫度相差很大，加上底板約束的作用，由此產(chǎn)生較大的溫度應(yīng)力致使閘墩出現(xiàn)裂縫。因此寒潮期間對(duì)閘墩采取適當(dāng)?shù)谋卮胧┲陵P(guān)重要，可有效減小溫度應(yīng)力，避免閘墩出現(xiàn)裂縫。

[1]由國(guó)文， 郭磊， 陳守開(kāi). 寒潮作用下大型水閘施工期溫控防裂仿真分析[J]. 水利水電科技進(jìn)展， 2015， 35(3)： 71-74．

[2]司政， 李守義， 陳培培， 等. 基于ANSYS的大體積混凝土溫度場(chǎng)計(jì)算程序開(kāi)發(fā)[J]. 長(zhǎng)江科學(xué)院院報(bào)， 2011， 29(9)： 53-56．

[3]李立峰. 基于ANSYS的混凝土水化熱溫度場(chǎng)讀取方法[J]. 山西建筑， 2007(1)： 74-75．

[4]王偉. ANSYS14. 0土木工程有限元分析從入門到精通[M]. 北京： 清華大學(xué)出版社， 2013．

[5]張朝暉. ANSYS12. 0熱分析工程應(yīng)用實(shí)戰(zhàn)手冊(cè)[M]. 北京： 中國(guó)鐵道出版社， 2010．

[6]王澤鵬. ANSYS12. 0熱力學(xué)有限元分析從入門到精通[M]. 北京： 機(jī)械工業(yè)出版社， 2010．

[7]王金龍. ANSYS12. 0有限元分析與范例解析[M]. 北京： 機(jī)械工業(yè)出版社， 2010．

[8]王暉， 劉超， 史志利. 鋼筋混凝土水池池壁施工溫度場(chǎng)仿真分析[J]. 特種結(jié)構(gòu)， 2007， 24(1)： 14-16．

[9]朱秋菊， 韓菊紅， 樂(lè)金朝. 閘墩施工期溫度應(yīng)力仿真分析[J]. 鄭州大學(xué)學(xué)報(bào)： 工學(xué)版， 2005， 26(1)： 47-49， 53．

[10]GB 50496—2009 大體積混凝土施工規(guī)范[s]．

[11]魏尊祥， 李志利. ANSYS模擬大體積混凝土澆筑過(guò)程的基本方法[J]. 城市建設(shè)理論研究： 電子版， 2013(14)， DOI： 10. 3969/j. issn. 2095—2104. 2013. 14. 192．

[12]趙英菊， 王社良， 康寧娟. ANSYS模擬大體積混凝土澆筑過(guò)程的參數(shù)分析[J]. 科技信息， 2007(14)： 96-97．

[13]朱伯芳. 大體積混凝土溫度應(yīng)力與溫度控制[M]. 北京： 中國(guó)電力出版社， 1999．

(責(zé)任編輯 王 琪)

Simulation Analysis on Temperature Control and Anti-cracking of Pier Based on Sub-master Model Correction Algorithm

YAN Bin, WANG Chuang

(College of Water Conservancy, Shenyang Agricultural University, Shenyang 110866, Liaoning, China)

The sub-master model correction algorithm is applied in APDL language simulation program to calculate the temperature and stress variation of concrete pier when using ANSYS software. The simulation calculation is encrypted in time and space, and each computation results of sub-model will be fed back to master pattern. The accuracies of temperature and stress fields calculated by the algorithm are improved. Applied this algorithm in the analysis of Shuangtaizi River sluice project, the simulation shows that the surface of sluice pier will cracked in Day 105 after slipform removal affected by cold, which consistent with actual situation. It is proved that the simulation calculation of temperature and stress fields of pier affected by cold is more accurate when using ANSYS combined with this algorithm．

sub-master model correction algorithm; crack; sluice pier; temperature stress; ANSYS; simulation analysis

2016-01-05

遼寧省重大科技計(jì)劃項(xiàng)目(2012212001)

閆濱(1972—)，女，遼寧沈陽(yáng)人，副教授，碩士生導(dǎo)師，博士，研究方向?yàn)樯鷳B(tài)水利及水工結(jié)構(gòu)健康監(jiān)控．

TV544

A

0559-9342(2016)11-0062-06