漫湾电站坝体力学参数的反演分析
|
摘要
本文首先阐述了电站在社会经济中所起的重要作用以及建成后坝体具有危害性一面。接着介绍反分析和位移反分析的应用现状和发展动态，提出了应用位移反分析反演坝体力学参数的可行性。
介绍了有关漫湾电站的工程概况以及内部观测系统。由于本文是在观测资料的基础上进行的计算，观测资料的正确和精确性对坝体弹性模量反演结果有直接的影响。因此对观测值的处理方法进行了介绍，对漫湾电站的观测资料进行全面分析，确定出漫湾电站目前的运行状态是在弹性阶段，以及坝体与基岩之间的关系。
为了计算能够顺利进行，建立模型是很重要的。建模时进行了模型的简化以及单元类型的确立、材料基本参数的取得、边界条件和约束条件的确定。进行模型检验，确定了计算模型。
从有限元基本方程出发，采用最小二乘法作为误差控制准则，推出弹性模量的迭代式。在ANSYS中对其进行二次开发，通过利用ANSYS中的参数化设计语言APDL编写了关于漫湾电站的建模和计算程序，并且ANSYS中得以很好的运行，最终得到38个时刻的弹性模量值。并对计算得到的结果进行分类讨论分析。计算结果表明弹性模量经过最初一年的增长后就进入衰减期，但是衰减趋势逐渐平缓。
最后，在上述研究成果的基础上，总结全文。对需要改进之处及今后进一步开展的研究工作提出有益的建议。
At first, the important influence of the power plants and the danger of the built dams are expatiated in the social economy in this paper. The application and development of the back analysis and the back analysis of displacement are introduced and the research results of the field in our country and on the world are summarized. The feasibility of applying the back analysis of displacement to back analysis the mechanical parameters in the dams are proposed.
Secondly, the general situation and the inner observed systems of Manwan are introduced. The calculation of this paper is based on these observed data and the validity and preciseness of these data is very important to the result, so that the method dealing of these data is bring forward. The observed information is fully analyzed, therefore we can confirm that the present running state of Manwan hydroelectric power plant is elastic state.
The built model is very important to carry through the later calculation. The simplify of model, the selection of element type, the attainment of the basic parameters of material, the condition of boundary and restriction and so on is disposed. The model is verified after the employment.
Based on the basic function of FEA, the iterated formula of modulus of elasticity is deduced through mathematical method as deviation control rule. The second exploitation is made on ANSYS and the procedure of building model and calculation is complied using ANSYS Parameter Design Language-APDL. The procedure has finished thirty-eight times calculation. The results and attain some rule are analyzed. The results of the calculation indicate that the modulus of elasticity of Manwan dam has entered the attenuation moment after increase and the attenuation will gradually lessen.
At last, on the basis of analysis result the conclusion is drawn and the useful recommendation and perspectives are put forward.
介绍了有关漫湾电站的工程概况以及内部观测系统。由于本文是在观测资料的基础上进行的计算，观测资料的正确和精确性对坝体弹性模量反演结果有直接的影响。因此对观测值的处理方法进行了介绍，对漫湾电站的观测资料进行全面分析，确定出漫湾电站目前的运行状态是在弹性阶段，以及坝体与基岩之间的关系。
为了计算能够顺利进行，建立模型是很重要的。建模时进行了模型的简化以及单元类型的确立、材料基本参数的取得、边界条件和约束条件的确定。进行模型检验，确定了计算模型。
从有限元基本方程出发，采用最小二乘法作为误差控制准则，推出弹性模量的迭代式。在ANSYS中对其进行二次开发，通过利用ANSYS中的参数化设计语言APDL编写了关于漫湾电站的建模和计算程序，并且ANSYS中得以很好的运行，最终得到38个时刻的弹性模量值。并对计算得到的结果进行分类讨论分析。计算结果表明弹性模量经过最初一年的增长后就进入衰减期，但是衰减趋势逐渐平缓。
最后，在上述研究成果的基础上，总结全文。对需要改进之处及今后进一步开展的研究工作提出有益的建议。
At first, the important influence of the power plants and the danger of the built dams are expatiated in the social economy in this paper. The application and development of the back analysis and the back analysis of displacement are introduced and the research results of the field in our country and on the world are summarized. The feasibility of applying the back analysis of displacement to back analysis the mechanical parameters in the dams are proposed.
Secondly, the general situation and the inner observed systems of Manwan are introduced. The calculation of this paper is based on these observed data and the validity and preciseness of these data is very important to the result, so that the method dealing of these data is bring forward. The observed information is fully analyzed, therefore we can confirm that the present running state of Manwan hydroelectric power plant is elastic state.
The built model is very important to carry through the later calculation. The simplify of model, the selection of element type, the attainment of the basic parameters of material, the condition of boundary and restriction and so on is disposed. The model is verified after the employment.
Based on the basic function of FEA, the iterated formula of modulus of elasticity is deduced through mathematical method as deviation control rule. The second exploitation is made on ANSYS and the procedure of building model and calculation is complied using ANSYS Parameter Design Language-APDL. The procedure has finished thirty-eight times calculation. The results and attain some rule are analyzed. The results of the calculation indicate that the modulus of elasticity of Manwan dam has entered the attenuation moment after increase and the attenuation will gradually lessen.
At last, on the basis of analysis result the conclusion is drawn and the useful recommendation and perspectives are put forward.
引文
[1] 廖奇志，阮娅，廖远志，大坝在人类改造自然可持续发展中的作用．人民长江，Vol．32 No．12，2001．12
[2] 梁宜编译，关于大坝建设利与弊的研究综述．水利电力科技，Vol 28．1
[3] 程念高，我国坝工建设的发展、展望和挑战，大坝与安全，2001．2
[4] [印度]杰伊·纳拉杨·维亚斯，大坝和可持续发展．水利水电快报，Vol．23 No．7，2002．4
[5] 吴中如，老坝病变和机理探讨，见：庆祝中国水利杂志创刊50周年面向21世纪中国水利发展战略研讨会论文．2000．9
[6] 朱志强，抓住机遇，加快开发云南 水电资源，积极推进西电东送战略．中国电力新闻网，2002．3
[7] 新华网云南频道-小湾电站专题，抓住机遇，多方集资，加大培植力度，使电力建设跨上新台阶．2002．1
[8] 经济新闻，云南提出要把电力产业建成第一大支柱产业．2003．4
[9] 张俊芝，在役钢筋混凝土结构的耐久性与概率寿命，南昌：南昌水专学报，VOL．20 No．3，2001．9
[10] 结构设计安全度专题讨论综述，结构设计安全度专题讨论会．土木工程学报，Vol 32，6，75～77，1999．6
[11] 曹楚生，从大坝设计和风险分析看大坝安全，水利水电工程设计，Vol．19 No．1，2000
[12] 张俊芝，李桂青，服役重力坝系统可靠度及概率寿命探讨．水利学报，2000．4
[13] 李光伟，混凝土抗裂能力的评价．水利水电科技进展，VOL．21 No．2，2001．4
[14] 王成启，吴科如，不伺弹性模量的纤维对高强混凝土力学性能的影响．混凝土与水泥制品，No．3，2002．3
[15] [荷兰]C.C.安托诺夫等，老化混凝土坝运行能力评价准则．俄刊水工建设，No．2，1995
[16] 于庆荣，颜德姮，程文爙等主编，混凝土结构．北京：中国建筑工业出版社，1997．6
[17] 王时越，张立翔，徐人平，弹性模量对混凝上疲劳性能的影响．昆明：昆明理工大学学报，Vol．26 No．5，2001．10
[18] 尚伟，雷声隆，既有混凝土弹性模量的证据理论估计方法．武汉：武汉工业大学学报，Vol．20 No．1，1998．3
[19] 顾冲时，蔡新，吴中如，探讨混凝土坝空间位移的正反分析模型[J]．工程力学，1997，14(2)：138-144．
[20] 刘迎曦，王登刚，张家良等，材料物性参数识别的梯度正则化方法[J]．计算力学学报，2000，17(1)：69-75．
[21] 刘迎曦，王登刚，李守巨等，混凝土重力坝弹性模量识别的一种新方法[J]．大连：大连理工大学学报，2000，40(2)：144-147．
[22] 岳建平，华锡生．坝体综合模量反演中水位等因素的影响[J]．南京：河海大学学报，1994，22(1)：99-101
[23] 李守巨，刘迎曦，王登刚等，混凝土重力坝材料参数识别的正则化最小二乘法．计算物理，Vol．17 No．6，2000．11
[24] Gioda G, Maier G, Direct search soluton of an inverse problem in elasticplasticity: identifcation of cohesion, friction angle and in-situ stress by pressure tunnel tests, Int. J. Num, Meth. Engrg, 1980,15:1822～1848
[25] 吴中如，沈长松，阮焕祥，论混凝土坝变形统计模型的因子选择．南京：河海大学学报，1988(6)
[26] 吴中如，论混凝土坝安全监控的确定性模型和混合模型．水利学报，1989(5)
[27] 吴中如，陈继禹，范树平，用反演分析法推求连拱混凝土地力学参数和断裂韧度．大坝观测与土工测试，1986．1
[28] 顾冲时，吴中如，蔡新，探讨混凝土坝空间位移场的正反分析模型．工程力学，1997
[29] 杨林德著，岩土工程问题的反演理论与工程实践．北京：科学出版社，1999
[30] 李金光，初始地应力场的反演分析研究：[硕士论文]．武汉：武汉大学，2001
[31] 陈斌，岩土工程随机反演分析及工程应用：[博士论文]．南京：河海大学，2001
[32] Kavangh K and Clough R. Finite element Application in the characterizaton of Elastic solids, Int. J. solids structures, 1971. 7, 11～13
[33] Kirsten H.A. Determination of Rock Mass Elastic Moduli by Back Analysis od Deformation Measurements, Exporation for Rock Engineering. ed by Z,T. Bieniamshi Balkema, 1976, 1:165～172
[34] Giode G, Pandotfi A, Cividini A. A comparative evalution of some back analysis algorithms and their application to in-situ load tests, In:Proc 2nd Int Symp on Field Measurement in Geom, 1987, 1131～1144
[35] Giode G, Jurina G, Numerical Identification of soll structure interaction pressures. Int J for Num & Anal Meth in Geomech, 1981.5, 33～56
[36] Sakurai S, Abe S.A design approach to dimensioning underground openings. 3rd international conference in 1979, Numerical Methods in Geomechanics, Aachen, 649～661
[37] Sakurai S, Takeuchi K. Back analysis of measured displacement of tunnel. Rock Mech and Rock Engrg, 1983,16; 173～180
[38] 顾冲时，吴中如，蔡新，探讨混凝土坝空间位移场的正反分析模型．工程力学，1997
[39] Wang Zhiyin, Liu Huaiheng, Back analysis of measured rheologic displacements of underground openings. In:Proc 6th Conf on Num Meth in Geom, Austria, 1988.2291～2297
[40] Wang Zhiyin, Li Yunpeng, Back analysis of viscoparameters andstrata stress in underground openings. In:Proc Iht Symp on Underground Eng, New Delhi, 1988,181～186
[41] Wang Sijing, Yang Zhifa, The back analysis method from displacement for viscoelastic rock mass, In:Proc 2nd Int Symp on FMGM'87,1987,1059～1068
[42] 刘怀恒，地下工程位移反分析原理应用及其反分析．西安：西安矿业学院学报，1989，9(1)，17～24
[43] 李云鹏，王芝银，粘弹性位移反分析的边界元法．西安：西安矿业学
院学报，1988，8(3)：1～10
[44] 郭乙木，万力等，坝体力学参数的混合模型优化反演法及其应用，水电能源科学，2001，19(2)，36～38
[45] 杨代泉，林见，有限元法反演坝体混凝土平均弹性模量E_(cr)和热膨胀系数α_(cr)．14～19
[46] 吴中如，沈长松，阮焕祥，论混凝土坝变形同济模型的因子选．南京：河海大学学报，1988(6)
[47] 吴中如，论混凝土坝安全监控的确定性模型和混合模型．水利学报，1989(5)
[48] 吴中如，陈继禹，范树平，用反演分析法推求连拱混凝土地力学参数和断裂韧度．大坝观测与土工测试，1986．1
[49] 顾冲时，吴中如，蔡新，探讨混凝土坝空间位移场的正反分析模型．工程力学，1997
[50] 雷在林，漫湾水电站重力坝运行状态分析：[硕士论文]．昆明：昆明理工大学，2003
[51] 熊云红，漫湾电站大坝内部观测自动化简介．大坝观测与土工测试，1997，21(3)：29～31
[52] 王崇华，漫湾电站倒垂设施施工与竣工．葛洲坝水电，1996，(2)：91～92
[53] 漫湾电站观测班，坝基扬压力观测成果分析(报告)．1999
[54] 能源部，水利部．混凝土大坝安全监测技术规范(SDJ336-89)(试行)．北京：水利电力出版社，1990
[55] 湖北省水利发电学会，漫湾水电工程大坝安全监测资料反馈分析(报告)．1997
[56] 长江科学院，漫湾水电工程大坝内部监测资料综合分析(报告)．1995年
[57] 凌川，冯兴常．漫湾水电站工程大坝安全监测成果分析．云南水力发电，1996，(2)：81～91
[58] B.L.阿姆斯塔特，可靠性数学(M)．北京：科学出版社，1978
[59] P.R.贝文顿，数据处理和误差发现(M)．知识出版社，1986
[60] 赵从宝，实用材料试验数据处理(M)．北京：经济管理出版社，1992
[61] 龚曙光，ANSYS基础应用及范例解析．北京：机械工业出版社，2003
[62] 徐芝纶，弹性力学．北京：高等教育出版社，1990
[63] 顾冲时，郑东健，吴中如，大坝及岩基的仿真有限元模型的精度分析．水电能源科学，Vol．20 No．2，2002．6
[64] 刘北辰，工程计算力学——理论与应用．北京：机械工业出版社，1994
[65] 朱伯芳，有限单元法原理与应用(第二版)，中国水利水电出版社，1998
[66] 谭建国，使用ANSYS6.0进行有限元分析．北京：北京大学出版社，2002．5
[67] 长江葛洲坝工程局中心试验室．漫湾水电工砼性能试验与施工配合比设计报告(报告)．1990
[68] 葛洲坝漫湾工程施工局试验室．漫湾水电站混凝土施工资料汇总表(报告)．1992
[69] 漫湾水电站工程管理局，长江科技开发公司技术开发部．云南漫湾水电站工程大坝变形、渗流监测成果分析(1995年度)(报告)．1996
[70] 胡荣辉，张五禄，水工建筑物．水利电力出版社，1993．6
[71] 吴中如，阮焕祥，混凝土坝观测资料的反分析．南京：河海大学学报，Vol．17 No．2，1989．3
[72] 刘迎曦，李守巨，陈昌林等，丰满混凝土重力坝材料参数识别研究．水利学报，No．9，1999．9
[73] Maier G, Jurina L, Podolak K. On model identification problem in rock Complex Formations[C].Capri, 1977: 257～261
[2] 梁宜编译，关于大坝建设利与弊的研究综述．水利电力科技，Vol 28．1
[3] 程念高，我国坝工建设的发展、展望和挑战，大坝与安全，2001．2
[4] [印度]杰伊·纳拉杨·维亚斯，大坝和可持续发展．水利水电快报，Vol．23 No．7，2002．4
[5] 吴中如，老坝病变和机理探讨，见：庆祝中国水利杂志创刊50周年面向21世纪中国水利发展战略研讨会论文．2000．9
[6] 朱志强，抓住机遇，加快开发云南 水电资源，积极推进西电东送战略．中国电力新闻网，2002．3
[7] 新华网云南频道-小湾电站专题，抓住机遇，多方集资，加大培植力度，使电力建设跨上新台阶．2002．1
[8] 经济新闻，云南提出要把电力产业建成第一大支柱产业．2003．4
[9] 张俊芝，在役钢筋混凝土结构的耐久性与概率寿命，南昌：南昌水专学报，VOL．20 No．3，2001．9
[10] 结构设计安全度专题讨论综述，结构设计安全度专题讨论会．土木工程学报，Vol 32，6，75～77，1999．6
[11] 曹楚生，从大坝设计和风险分析看大坝安全，水利水电工程设计，Vol．19 No．1，2000
[12] 张俊芝，李桂青，服役重力坝系统可靠度及概率寿命探讨．水利学报，2000．4
[13] 李光伟，混凝土抗裂能力的评价．水利水电科技进展，VOL．21 No．2，2001．4
[14] 王成启，吴科如，不伺弹性模量的纤维对高强混凝土力学性能的影响．混凝土与水泥制品，No．3，2002．3
[15] [荷兰]C.C.安托诺夫等，老化混凝土坝运行能力评价准则．俄刊水工建设，No．2，1995
[16] 于庆荣，颜德姮，程文爙等主编，混凝土结构．北京：中国建筑工业出版社，1997．6
[17] 王时越，张立翔，徐人平，弹性模量对混凝上疲劳性能的影响．昆明：昆明理工大学学报，Vol．26 No．5，2001．10
[18] 尚伟，雷声隆，既有混凝土弹性模量的证据理论估计方法．武汉：武汉工业大学学报，Vol．20 No．1，1998．3
[19] 顾冲时，蔡新，吴中如，探讨混凝土坝空间位移的正反分析模型[J]．工程力学，1997，14(2)：138-144．
[20] 刘迎曦，王登刚，张家良等，材料物性参数识别的梯度正则化方法[J]．计算力学学报，2000，17(1)：69-75．
[21] 刘迎曦，王登刚，李守巨等，混凝土重力坝弹性模量识别的一种新方法[J]．大连：大连理工大学学报，2000，40(2)：144-147．
[22] 岳建平，华锡生．坝体综合模量反演中水位等因素的影响[J]．南京：河海大学学报，1994，22(1)：99-101
[23] 李守巨，刘迎曦，王登刚等，混凝土重力坝材料参数识别的正则化最小二乘法．计算物理，Vol．17 No．6，2000．11
[24] Gioda G, Maier G, Direct search soluton of an inverse problem in elasticplasticity: identifcation of cohesion, friction angle and in-situ stress by pressure tunnel tests, Int. J. Num, Meth. Engrg, 1980,15:1822～1848
[25] 吴中如，沈长松，阮焕祥，论混凝土坝变形统计模型的因子选择．南京：河海大学学报，1988(6)
[26] 吴中如，论混凝土坝安全监控的确定性模型和混合模型．水利学报，1989(5)
[27] 吴中如，陈继禹，范树平，用反演分析法推求连拱混凝土地力学参数和断裂韧度．大坝观测与土工测试，1986．1
[28] 顾冲时，吴中如，蔡新，探讨混凝土坝空间位移场的正反分析模型．工程力学，1997
[29] 杨林德著，岩土工程问题的反演理论与工程实践．北京：科学出版社，1999
[30] 李金光，初始地应力场的反演分析研究：[硕士论文]．武汉：武汉大学，2001
[31] 陈斌，岩土工程随机反演分析及工程应用：[博士论文]．南京：河海大学，2001
[32] Kavangh K and Clough R. Finite element Application in the characterizaton of Elastic solids, Int. J. solids structures, 1971. 7, 11～13
[33] Kirsten H.A. Determination of Rock Mass Elastic Moduli by Back Analysis od Deformation Measurements, Exporation for Rock Engineering. ed by Z,T. Bieniamshi Balkema, 1976, 1:165～172
[34] Giode G, Pandotfi A, Cividini A. A comparative evalution of some back analysis algorithms and their application to in-situ load tests, In:Proc 2nd Int Symp on Field Measurement in Geom, 1987, 1131～1144
[35] Giode G, Jurina G, Numerical Identification of soll structure interaction pressures. Int J for Num & Anal Meth in Geomech, 1981.5, 33～56
[36] Sakurai S, Abe S.A design approach to dimensioning underground openings. 3rd international conference in 1979, Numerical Methods in Geomechanics, Aachen, 649～661
[37] Sakurai S, Takeuchi K. Back analysis of measured displacement of tunnel. Rock Mech and Rock Engrg, 1983,16; 173～180
[38] 顾冲时，吴中如，蔡新，探讨混凝土坝空间位移场的正反分析模型．工程力学，1997
[39] Wang Zhiyin, Liu Huaiheng, Back analysis of measured rheologic displacements of underground openings. In:Proc 6th Conf on Num Meth in Geom, Austria, 1988.2291～2297
[40] Wang Zhiyin, Li Yunpeng, Back analysis of viscoparameters andstrata stress in underground openings. In:Proc Iht Symp on Underground Eng, New Delhi, 1988,181～186
[41] Wang Sijing, Yang Zhifa, The back analysis method from displacement for viscoelastic rock mass, In:Proc 2nd Int Symp on FMGM'87,1987,1059～1068
[42] 刘怀恒，地下工程位移反分析原理应用及其反分析．西安：西安矿业学院学报，1989，9(1)，17～24
[43] 李云鹏，王芝银，粘弹性位移反分析的边界元法．西安：西安矿业学
院学报，1988，8(3)：1～10
[44] 郭乙木，万力等，坝体力学参数的混合模型优化反演法及其应用，水电能源科学，2001，19(2)，36～38
[45] 杨代泉，林见，有限元法反演坝体混凝土平均弹性模量E_(cr)和热膨胀系数α_(cr)．14～19
[46] 吴中如，沈长松，阮焕祥，论混凝土坝变形同济模型的因子选．南京：河海大学学报，1988(6)
[47] 吴中如，论混凝土坝安全监控的确定性模型和混合模型．水利学报，1989(5)
[48] 吴中如，陈继禹，范树平，用反演分析法推求连拱混凝土地力学参数和断裂韧度．大坝观测与土工测试，1986．1
[49] 顾冲时，吴中如，蔡新，探讨混凝土坝空间位移场的正反分析模型．工程力学，1997
[50] 雷在林，漫湾水电站重力坝运行状态分析：[硕士论文]．昆明：昆明理工大学，2003
[51] 熊云红，漫湾电站大坝内部观测自动化简介．大坝观测与土工测试，1997，21(3)：29～31
[52] 王崇华，漫湾电站倒垂设施施工与竣工．葛洲坝水电，1996，(2)：91～92
[53] 漫湾电站观测班，坝基扬压力观测成果分析(报告)．1999
[54] 能源部，水利部．混凝土大坝安全监测技术规范(SDJ336-89)(试行)．北京：水利电力出版社，1990
[55] 湖北省水利发电学会，漫湾水电工程大坝安全监测资料反馈分析(报告)．1997
[56] 长江科学院，漫湾水电工程大坝内部监测资料综合分析(报告)．1995年
[57] 凌川，冯兴常．漫湾水电站工程大坝安全监测成果分析．云南水力发电，1996，(2)：81～91
[58] B.L.阿姆斯塔特，可靠性数学(M)．北京：科学出版社，1978
[59] P.R.贝文顿，数据处理和误差发现(M)．知识出版社，1986
[60] 赵从宝，实用材料试验数据处理(M)．北京：经济管理出版社，1992
[61] 龚曙光，ANSYS基础应用及范例解析．北京：机械工业出版社，2003
[62] 徐芝纶，弹性力学．北京：高等教育出版社，1990
[63] 顾冲时，郑东健，吴中如，大坝及岩基的仿真有限元模型的精度分析．水电能源科学，Vol．20 No．2，2002．6
[64] 刘北辰，工程计算力学——理论与应用．北京：机械工业出版社，1994
[65] 朱伯芳，有限单元法原理与应用(第二版)，中国水利水电出版社，1998
[66] 谭建国，使用ANSYS6.0进行有限元分析．北京：北京大学出版社，2002．5
[67] 长江葛洲坝工程局中心试验室．漫湾水电工砼性能试验与施工配合比设计报告(报告)．1990
[68] 葛洲坝漫湾工程施工局试验室．漫湾水电站混凝土施工资料汇总表(报告)．1992
[69] 漫湾水电站工程管理局，长江科技开发公司技术开发部．云南漫湾水电站工程大坝变形、渗流监测成果分析(1995年度)(报告)．1996
[70] 胡荣辉，张五禄，水工建筑物．水利电力出版社，1993．6
[71] 吴中如，阮焕祥，混凝土坝观测资料的反分析．南京：河海大学学报，Vol．17 No．2，1989．3
[72] 刘迎曦，李守巨，陈昌林等，丰满混凝土重力坝材料参数识别研究．水利学报，No．9，1999．9
[73] Maier G, Jurina L, Podolak K. On model identification problem in rock Complex Formations[C].Capri, 1977: 257～261