施工期面板堆石坝临时过水度汛断面渗流场与应力场的耦合分析[摘要]施工期面板堆石坝采用临时断面过水度汛施工方案时，其临时断面存在的渗流场与应力场是相互作用、相互影响的，为了研究施工期面板堆石坝临时度汛断面渗流场与应力场的变化规律，研究提出二者相互作用的计算方法，本文结合某面板堆石坝工程实例，采用非线性有限元法，对施工期临时过水度汛断面进行了渗流场与应力场的耦合分析，并与不考虑耦合作用的计算结果进行了对比。[关键词]施工期；临时断面过水度汛；应力变形；渗流；耦合分析[abstract]therockfilldamduringtheusingprovisionalsectionwaterhighalertwhenconstructionscheme,thetemporarysectionoftheexistingseepagefieldandstressfieldisinteractionandinfluenceeachother,inordertostudytherockfilldamtemporaryhighalertduringtheseepagefieldandstressfieldofthecrosssectionchangerule,andputforwardthecalculationmethodofinteractionbetween,combinedwithafacerockfilldamprojectexample,byusingnonlinearfiniteelementmethod,thetemporarywaterduringthehighalerttheseepagefieldandthecrosssectionofstressfieldcouplinganalysis,andwithdon'tconsiderthecoupledactionofthecalculationresultsarecompared.[keywords]constructionstage;Temporarysectionwaterhighalert;Stressdeformation;Seepage;Couplinganalysis：TU74文献标识码：A：1计算原理1.1渗流场与应力场耦合分析的基本原理1.1.1渗流场与应力场耦合分析的数学模型[3][4]应力场影响下稳定渗流场的数学模型为：(1)式中:为渗流区域；H（x，y，z）为渗流区域内的水头分布函数；为各向同性的渗透系数，是应力的函数，即；是第一类边界，为其水头分布；为第二类边界，为法线方向，为上的流量分布；为第三类边界条件（即渗流自由面边界），为法线方向。渗流场影响下应力场的数学模型为：（2）式中，为包括渗流自由面以上坝体的整个区域；为应力张量场，为应变张量场；为体力，是渗流场水头分布的函数；为位移场；为常数；为体积应变；为Delta置换坐标；为已知面力边界，为其法线方向余弦，为其上已知面力分布，也是渗流场水头H（x，y，z）的函数，由式（5-1）确定；为已知位移边界，为其上已知位移分布。式（1）和式（2）联立即构成渗流场与应力场耦合分析的数学模型。1.1.2耦合分析的方法与步骤目前，常用的耦合分析主要有间接耦合法和直接耦合法两种方法。由于直接耦合法在进行应力变形非线性方程组求解时收敛性较差，因此对于涉及材料非线性问题研究时，目前大多采用间接耦合法[24]。根据本文所研究问题的特点，本文即采用间接耦合法，实现步骤如下：（1）假定坝体及坝基的初始应力场为零，计算得到单元初始孔隙水压力及水头分布；（2）将每一时段单元孔隙水压力变化转化为等效结点荷载；（3）计算由于每一时段等效结点荷载变化所引起的坝体应力变形的变化；（4）计算由于坝体应力变形的变化而引起的渗透系数的改变，根据各单元新的渗透系数计算出新的孔隙水压力，并得出新的水头分布；（5）判断迭代过程是否收敛，其收敛条件为：前后迭代得到的位移或渗透水头非常接近时，迭代结束，否则返回步骤（2）继续迭代。1.2汛期度汛过程的考虑施工期坝面过水堆石坝是在汛期未修筑面板的情况下由坝面过水度汛，由上游垫层及坝面堆石承受水压力。坝体过水时，坝体上部水流作用于护面的力有：时均动水压力、脉动动水压力以及拖曳力等【8】。拟静力法是一种用静力学方法近似解决动力学问题的简易方法，采用这种方法，计算方法较为简单，且参数也易于确定，故本文对其应力变形的研究，采用拟静力法简化坝体上部水体对坝体的作用，也就是坝面过流时，将坝面上部的水体对坝体产生的力等效于静水压力产生的力来处理。由于在过流前，需对临时断面加以防护，坝面设铅丝笼网格护面，考虑到铅丝笼结构简单、透水性大，且适应坝体沉降较好，其高度相比坝高较小，因此在汛期度汛过程中不考虑铅丝笼对坝体的作用。2工程实例2.1工程概况某混凝土面板堆石坝[9]坝顶高程423.6m，河床面高程为355.0m，趾板建基面高程3...