深覆盖层地基高面板堆石坝应力变形动力有限元分析
|
摘要
随着我国国民经济建设的发展和西部大开发战略的实施,建在深覆盖层地基上的混凝土面板堆石坝工程越来越多。特别是位于西南、西北等地区的面板堆石坝大多建在砂土、砂砾石或碎石土等组成的深覆盖层地基上。因此研究这些位于强震区面板堆石坝石坝在地震频率高、强度大的情况下的应力变形问题,特别是抗震安全问题显得尤为重要。
本文采用有限元数值分析方法,结合堆石料与接触面的静动力本构模型,对面板堆石坝进行三维有限元静动力计算分析。首先根据面板堆石坝的结构特点,研究模拟覆盖层、堆石坝体、面板、堆石与面板接触面等各构成部分工作状态的有限单元形式。在分析面板堆石坝抗震研究现状的基础上,引入堆石料的静、动力本构模型及接触界面的本构模型。研究了面板堆石坝应力变形计算中诸如接触、分级加载、面板裂缝与接缝等特殊边界模拟的问题,在ADINA软件中实现了上述特殊边界的模拟。最后,结合实际工程,利用大型有限元软件ADINA建立了面板堆石坝三维静、动力分析的有限元模型；进行了施工期与运行期混凝土面板堆石坝应力变形的三维有限元静力仿真分析。动力分析中堆石料的动力本构模型采用Hardin—Drnevich双曲线本构模型,由于堆石料本构模型的复杂性和有限元软件的局限性,通过将非线性问题线性化,计算分析了大坝在地震时的加速度、剪应力和剪应变等。
本文所建立或提出的上述模型及方法,为强震区深覆盖层上混凝土面板堆石坝在地震情况下的应力变形特性研究和预测进行了一些理论探索,所获得的分析计算成果对类似工程也具有一定的参考和应用价值。
With the development of China's national economy and the western development strategy, more and more earth-rock dam engineering are designed in deep alluvial depoist, and the scale of the project is increasing. Especially in southwest and northwest and other regions, most of the water conservancy and hydropower projects are built on deep overburden foundation which composed of sand, gravel or crushed gravel. Study deformation and stability of these earth-rock dams located in high earthquake zone, in high frequency and strength earthquake situation, the issue of earthquake resistance safety is particularly prominent.
In this paper, use finite element numerical analysis method, combined with static and dynamic constitutive model and non-linear algorithm of rock fill and interface, carry out static and dynamic three-dimensional calculation and analysis. First of all, combined with CFRD structural characteristics, research and simulate finite element model coverage layer, rock fill dam body, face slab, contact face between face slab and rock fill and other components in the work state. On the basis of analyzing the CFRD earthquake resistance research, introduce static, dynamic constitutive model of rock fill and constitutive model of contact interface. Based on the research and analysis of nonlinear finite element calculation, some special boundary analogy issues such as contacts, grading load, face slab cracks and joints and so on are discussed, the above-mentioned special boundary analogy is achieved. Finally, in light of the actual engineering, using large-scale finite element software ADINA, established a three-dimensional static and dynamic finite element analysis model to CFRD. On basis of this, carried out three-dimensional finite element analysis in static simulation of CFRD stress and deformation during construction and operation periods. In dynamic analysis, dynamic constitutive model used Hardin-Reich hyperbolic constitutive model, as a result of specificity of rock fill constitutive model and limitations of finite element software, draw up external procedures made nonlinear problems lineated, lay the foundation for using equivalent linear model. By external procedures call ADINA software program iterative calculations, reasonably determine earthquake acceleration, shear stress and shear strain and so on in earthquake time.
In this paper, the above mentioned models and methods lay a necessary theoretical basis for the research and forecasting deformation and stability of concrete faced rock fill dam located in high earthquake zone and deep alluvial depoist, and the analysis calculation achievement also has an important reference and application value to the similar projects.
本文采用有限元数值分析方法,结合堆石料与接触面的静动力本构模型,对面板堆石坝进行三维有限元静动力计算分析。首先根据面板堆石坝的结构特点,研究模拟覆盖层、堆石坝体、面板、堆石与面板接触面等各构成部分工作状态的有限单元形式。在分析面板堆石坝抗震研究现状的基础上,引入堆石料的静、动力本构模型及接触界面的本构模型。研究了面板堆石坝应力变形计算中诸如接触、分级加载、面板裂缝与接缝等特殊边界模拟的问题,在ADINA软件中实现了上述特殊边界的模拟。最后,结合实际工程,利用大型有限元软件ADINA建立了面板堆石坝三维静、动力分析的有限元模型；进行了施工期与运行期混凝土面板堆石坝应力变形的三维有限元静力仿真分析。动力分析中堆石料的动力本构模型采用Hardin—Drnevich双曲线本构模型,由于堆石料本构模型的复杂性和有限元软件的局限性,通过将非线性问题线性化,计算分析了大坝在地震时的加速度、剪应力和剪应变等。
本文所建立或提出的上述模型及方法,为强震区深覆盖层上混凝土面板堆石坝在地震情况下的应力变形特性研究和预测进行了一些理论探索,所获得的分析计算成果对类似工程也具有一定的参考和应用价值。
With the development of China's national economy and the western development strategy, more and more earth-rock dam engineering are designed in deep alluvial depoist, and the scale of the project is increasing. Especially in southwest and northwest and other regions, most of the water conservancy and hydropower projects are built on deep overburden foundation which composed of sand, gravel or crushed gravel. Study deformation and stability of these earth-rock dams located in high earthquake zone, in high frequency and strength earthquake situation, the issue of earthquake resistance safety is particularly prominent.
In this paper, use finite element numerical analysis method, combined with static and dynamic constitutive model and non-linear algorithm of rock fill and interface, carry out static and dynamic three-dimensional calculation and analysis. First of all, combined with CFRD structural characteristics, research and simulate finite element model coverage layer, rock fill dam body, face slab, contact face between face slab and rock fill and other components in the work state. On the basis of analyzing the CFRD earthquake resistance research, introduce static, dynamic constitutive model of rock fill and constitutive model of contact interface. Based on the research and analysis of nonlinear finite element calculation, some special boundary analogy issues such as contacts, grading load, face slab cracks and joints and so on are discussed, the above-mentioned special boundary analogy is achieved. Finally, in light of the actual engineering, using large-scale finite element software ADINA, established a three-dimensional static and dynamic finite element analysis model to CFRD. On basis of this, carried out three-dimensional finite element analysis in static simulation of CFRD stress and deformation during construction and operation periods. In dynamic analysis, dynamic constitutive model used Hardin-Reich hyperbolic constitutive model, as a result of specificity of rock fill constitutive model and limitations of finite element software, draw up external procedures made nonlinear problems lineated, lay the foundation for using equivalent linear model. By external procedures call ADINA software program iterative calculations, reasonably determine earthquake acceleration, shear stress and shear strain and so on in earthquake time.
In this paper, the above mentioned models and methods lay a necessary theoretical basis for the research and forecasting deformation and stability of concrete faced rock fill dam located in high earthquake zone and deep alluvial depoist, and the analysis calculation achievement also has an important reference and application value to the similar projects.
引文
[1]赵增凯.我国深厚覆盖层上修建面板堆石坝的发展[J].面板堆石坝工程,2006,3
[2]郦能惠.高混凝土面板堆石坝新技术[M].北京：中国水利水电出版社,2007,06
[3]王世夏.水工设计的理论和方法[M].北京：中国水利水电出版社,2000,07
[4]赵剑明,常亚屏,陈宁.加强高土石坝抗震研究的现实意义及工作展望[J].世界地震工程,2004,01
[5]蒋国澄,傅志安,凤家骥.混凝土面板堆石坝工程[M].湖北：湖北科学技术出版社,1997
[6]郦能惠,米占宽等.深覆盖层地基上混凝土面板堆石坝有限元静动力应力变形分析报告[R].江苏：南京水利科学研究院研究,2003
[7]郦能惠,孙大伟,米占宽,张发中,李学强.覆盖层上高面板坝连接板长度优化分析[J].水利水电技术,2005,36
[8]安盛勋.察汗乌苏水电站工程趾板建在深覆盖层上混凝土面板堆石坝设计综述[C],蒋国澄,中国混凝土面板堆石坝20年.北京：中国水利水电出版社,2005：151-159
[9]徐泽平,邵宇,梁建辉.深覆盖层上面板堆石坝应力变形特性研究[c],蒋国澄,中国混凝土面板堆石坝20年.北京：中国水利水电出版社,2005：582-588
[10]郦能惠.李国英,形性状,赵魁芝,陈铁林.强震区高面板堆石坝静力和动力应力和变形[J].岩土工程学报,2004,2
[11]汪小刚.西部水工程中的岩土工程问题[J].岩土工程学报,2007,29(8)：1129-1134
[12]金丰年,浦奎英.关于粘弹性模型的讨论[J].岩土力学与工程学报,2005,04
[13]赵剑明,常亚屏,陈宁,黄丽清.面板堆石坝地震分析的若干研究现状和进展及工作展望[C],蒋国澄,中国混凝土面板堆石坝20年.北京：中国水利水电出版社,2005：158-168
[14]李万红,汪闻韶.土石坝非线性动力反应分析[J].水利学报.1994,02
[15]赵剑明等.面板坝三维非线性动力分析研究[R].北京：中国水利水电科学研究院,2000
[16]赵剑明,常亚屏,陈宁.深厚覆盖层上土石坝抗震分析的三维真非线性方法及应用[C].岩土力学与工程学会.第一届中国水利水电岩土力学与工程学术讨论会论文集.北京：中国水利水电出版社,2007
[17]罗守成.对深覆盖层地质问题的认识[J].水力发电,1995,4：21-24
[18]陈海军.河谷深覆盖层工程地质特性及评价方法[J].地质灾害与环境保护,1996,7(4)
[19]金辉.西南地区河谷深覆盖层基本特征及成因机理研究[D].四川：成都理工大学,2008.
[20]钟诚昌.深厚覆盖层地区工程物探方法技术[J].水力发电,1996,5：34-38.
[21]孟永旭.下坂地水库坝址深厚覆盖层工程特性及主要地质问题初步评价[J].陕西水利水电技术,2000,2：48-51
[22]徐泽平.混凝土面板堆石坝应力变形特性研究[M].郑州：黄河水利出版社,2005
[23]水电水利规划设计总院.深覆盖层地基上混凝土面板堆石坝关键技术研究[R].北京：水电水利规划设计总院,2004
[24]司洪祥.论土石坝应力应变计算中的邓肯模型[J].南京：水利水运科学研究,1985,4：89-97
[25]Duncan J.M, Chang C Y. Nonlinear Analysis of Stress and Strain [J].Journal of Soil Mechanics and Foundation Division.1970,96(5):1629-1653
[26]陈慧远.土石坝有限元分析[M].南京：河海大学出版社,1988
[27]钱家欢,殷宗泽土.工原理与计算第二版[M].北京：中国水利水电出版社,2006
[28]刘信声,徐秉业.工程塑性力学选讲[M].北京：清华大学出版社,1986
[29]郑雨天.岩石力学的弹塑粘性理论基础[M].北京：煤炭工业出版社,1988
[30]张学言.岩土塑性力学[M].北京：人民交通出版社,1993
[31]麦家煊,薛继龙,王津龙.三维接触单元在堆石坝面板应力位移分析中的应用研究[J].水力发电学报[J],2002,14(1)：39-44
[32]汪召华,高莲士,宋文晶.三维摩擦接触的有限元分析[J].清华大学学报(自然科学版),2002,42(S1)：93-96
[33]沈珠江.面板堆石坝应力应变分析的若干问题[C],中国水力发电工程学会,中国混凝土面板堆石坝十年学术研讨会论文集.北京：中国水力发电工程学会,1995：43-50
[34]岑威钧,李星.面板坝数值分析中接触面模型与接缝模型述评[J].水力发电,2007,33(2)：38-41
[35]潘家军,费胜.混凝土面板堆石坝三维非线性有限元应力变形分析[J].水电能源科学,2007,25
[36]赵剑明,汪闻韶,混凝土面板堆石坝面板地震反应分析[J].岩石力学与工程学报,2001,20(s2)
[37]顾淦臣.土石坝地震工程[M].南京：河海大学出版社,1989
[38]徐植信,胡再龙.结构地震反应分析[M].北京：高等教育出版社,1993
[39]李桂青.抗震结构计算理论和方法[M].北京：地震出版社,1985
[40]史良.黄土隧道抗震设计研究[D].西安：长安大学,2005：30
[41]朱伯芳.有限单元法原理与应用[M].北京：中国水利水电出版社,1998：311-357
[42]R.W.Clough J.Penzien-Dynamics of structures [M].McGraw Inc.1975
[43]李炎隆,混凝土面板极端破坏情况下堆石坝渗流与应力变形特性研究[D].西安：西安理工大学,2008：50
[44]岳戈,陈权.ADINA应用基础与实例详解[M].北京：人民交通出版社,2008
[45]Birbraer.ANSeismic analysis of structure[M].Nauka:Sanket-petersburg,1998
[46]徐汉忠,吴旭光,马贞信.计算动水压力附加质量的韦氏公式的扩充.河海大学学报,1997,25(5)：120-122
[47]中华人民共和国行业标准,混凝土面板堆石坝设计规范SL228-98[S].北京：中国水利水电出版社,1999
[48]Lysmer. Response to buried Structures to Travelling Waves[J].Geotech. Eng. IJiv. ASCE,1981,20(5):183-200
[49]张燎军,张慧星,王大胜,王伟华.黏弹性人工边界在ADINA中的应用[J].世界地震工程,2008,03
[50]Deeks A J.Randolph M F.Axisymmetric time-domain transmitting Boundaries[J]. Journal of Engineering Mechanics,1994,120 (1):25-42
[51]刘晶波,李彬.三维黏弹性静-动力统一人工边界[J].中国科学E辑工程科学材料科学2005,35(9)：966-980
[52]Lysmer J. Kulemeyer R L.Finite dynamic model for infinite media.Journal of Engineering Mechanics, ASCE,1969,95:759-877
[53]Liu J B, Lu Y D. A direct method for analysis of dynamic soil-structure interaction based on interface ideaIn: Zhang C H, Wolf J P, eds.Dynamic Soil Structure Interaction[M].Beijing:International Academia Publishers,1997:258-273
[54]Wolf J P, Song C.Finite-Element Modeling of Unbounded Media[M].Newyork:Wiley,1996
[55]董红健甘肃洮河九甸峡水利枢纽工程设计[J].水利规划与设计,2006.05
[56]梁宗仁.甘肃九甸峡水利枢纽深厚覆盖层工程特性[J].水利规划与设计,2007,03
[57]蒋国澄.中国混凝土面板堆石坝20年[M].北京：中国水利水电出版社,2005：1-25
[58]吕生玺,沈振中,温续余.九甸峡混凝土面板堆石坝地震特性研究[J].中国农村水利水电,2008,02
[2]郦能惠.高混凝土面板堆石坝新技术[M].北京：中国水利水电出版社,2007,06
[3]王世夏.水工设计的理论和方法[M].北京：中国水利水电出版社,2000,07
[4]赵剑明,常亚屏,陈宁.加强高土石坝抗震研究的现实意义及工作展望[J].世界地震工程,2004,01
[5]蒋国澄,傅志安,凤家骥.混凝土面板堆石坝工程[M].湖北：湖北科学技术出版社,1997
[6]郦能惠,米占宽等.深覆盖层地基上混凝土面板堆石坝有限元静动力应力变形分析报告[R].江苏：南京水利科学研究院研究,2003
[7]郦能惠,孙大伟,米占宽,张发中,李学强.覆盖层上高面板坝连接板长度优化分析[J].水利水电技术,2005,36
[8]安盛勋.察汗乌苏水电站工程趾板建在深覆盖层上混凝土面板堆石坝设计综述[C],蒋国澄,中国混凝土面板堆石坝20年.北京：中国水利水电出版社,2005：151-159
[9]徐泽平,邵宇,梁建辉.深覆盖层上面板堆石坝应力变形特性研究[c],蒋国澄,中国混凝土面板堆石坝20年.北京：中国水利水电出版社,2005：582-588
[10]郦能惠.李国英,形性状,赵魁芝,陈铁林.强震区高面板堆石坝静力和动力应力和变形[J].岩土工程学报,2004,2
[11]汪小刚.西部水工程中的岩土工程问题[J].岩土工程学报,2007,29(8)：1129-1134
[12]金丰年,浦奎英.关于粘弹性模型的讨论[J].岩土力学与工程学报,2005,04
[13]赵剑明,常亚屏,陈宁,黄丽清.面板堆石坝地震分析的若干研究现状和进展及工作展望[C],蒋国澄,中国混凝土面板堆石坝20年.北京：中国水利水电出版社,2005：158-168
[14]李万红,汪闻韶.土石坝非线性动力反应分析[J].水利学报.1994,02
[15]赵剑明等.面板坝三维非线性动力分析研究[R].北京：中国水利水电科学研究院,2000
[16]赵剑明,常亚屏,陈宁.深厚覆盖层上土石坝抗震分析的三维真非线性方法及应用[C].岩土力学与工程学会.第一届中国水利水电岩土力学与工程学术讨论会论文集.北京：中国水利水电出版社,2007
[17]罗守成.对深覆盖层地质问题的认识[J].水力发电,1995,4：21-24
[18]陈海军.河谷深覆盖层工程地质特性及评价方法[J].地质灾害与环境保护,1996,7(4)
[19]金辉.西南地区河谷深覆盖层基本特征及成因机理研究[D].四川：成都理工大学,2008.
[20]钟诚昌.深厚覆盖层地区工程物探方法技术[J].水力发电,1996,5：34-38.
[21]孟永旭.下坂地水库坝址深厚覆盖层工程特性及主要地质问题初步评价[J].陕西水利水电技术,2000,2：48-51
[22]徐泽平.混凝土面板堆石坝应力变形特性研究[M].郑州：黄河水利出版社,2005
[23]水电水利规划设计总院.深覆盖层地基上混凝土面板堆石坝关键技术研究[R].北京：水电水利规划设计总院,2004
[24]司洪祥.论土石坝应力应变计算中的邓肯模型[J].南京：水利水运科学研究,1985,4：89-97
[25]Duncan J.M, Chang C Y. Nonlinear Analysis of Stress and Strain [J].Journal of Soil Mechanics and Foundation Division.1970,96(5):1629-1653
[26]陈慧远.土石坝有限元分析[M].南京：河海大学出版社,1988
[27]钱家欢,殷宗泽土.工原理与计算第二版[M].北京：中国水利水电出版社,2006
[28]刘信声,徐秉业.工程塑性力学选讲[M].北京：清华大学出版社,1986
[29]郑雨天.岩石力学的弹塑粘性理论基础[M].北京：煤炭工业出版社,1988
[30]张学言.岩土塑性力学[M].北京：人民交通出版社,1993
[31]麦家煊,薛继龙,王津龙.三维接触单元在堆石坝面板应力位移分析中的应用研究[J].水力发电学报[J],2002,14(1)：39-44
[32]汪召华,高莲士,宋文晶.三维摩擦接触的有限元分析[J].清华大学学报(自然科学版),2002,42(S1)：93-96
[33]沈珠江.面板堆石坝应力应变分析的若干问题[C],中国水力发电工程学会,中国混凝土面板堆石坝十年学术研讨会论文集.北京：中国水力发电工程学会,1995：43-50
[34]岑威钧,李星.面板坝数值分析中接触面模型与接缝模型述评[J].水力发电,2007,33(2)：38-41
[35]潘家军,费胜.混凝土面板堆石坝三维非线性有限元应力变形分析[J].水电能源科学,2007,25
[36]赵剑明,汪闻韶,混凝土面板堆石坝面板地震反应分析[J].岩石力学与工程学报,2001,20(s2)
[37]顾淦臣.土石坝地震工程[M].南京：河海大学出版社,1989
[38]徐植信,胡再龙.结构地震反应分析[M].北京：高等教育出版社,1993
[39]李桂青.抗震结构计算理论和方法[M].北京：地震出版社,1985
[40]史良.黄土隧道抗震设计研究[D].西安：长安大学,2005：30
[41]朱伯芳.有限单元法原理与应用[M].北京：中国水利水电出版社,1998：311-357
[42]R.W.Clough J.Penzien-Dynamics of structures [M].McGraw Inc.1975
[43]李炎隆,混凝土面板极端破坏情况下堆石坝渗流与应力变形特性研究[D].西安：西安理工大学,2008：50
[44]岳戈,陈权.ADINA应用基础与实例详解[M].北京：人民交通出版社,2008
[45]Birbraer.ANSeismic analysis of structure[M].Nauka:Sanket-petersburg,1998
[46]徐汉忠,吴旭光,马贞信.计算动水压力附加质量的韦氏公式的扩充.河海大学学报,1997,25(5)：120-122
[47]中华人民共和国行业标准,混凝土面板堆石坝设计规范SL228-98[S].北京：中国水利水电出版社,1999
[48]Lysmer. Response to buried Structures to Travelling Waves[J].Geotech. Eng. IJiv. ASCE,1981,20(5):183-200
[49]张燎军,张慧星,王大胜,王伟华.黏弹性人工边界在ADINA中的应用[J].世界地震工程,2008,03
[50]Deeks A J.Randolph M F.Axisymmetric time-domain transmitting Boundaries[J]. Journal of Engineering Mechanics,1994,120 (1):25-42
[51]刘晶波,李彬.三维黏弹性静-动力统一人工边界[J].中国科学E辑工程科学材料科学2005,35(9)：966-980
[52]Lysmer J. Kulemeyer R L.Finite dynamic model for infinite media.Journal of Engineering Mechanics, ASCE,1969,95:759-877
[53]Liu J B, Lu Y D. A direct method for analysis of dynamic soil-structure interaction based on interface ideaIn: Zhang C H, Wolf J P, eds.Dynamic Soil Structure Interaction[M].Beijing:International Academia Publishers,1997:258-273
[54]Wolf J P, Song C.Finite-Element Modeling of Unbounded Media[M].Newyork:Wiley,1996
[55]董红健甘肃洮河九甸峡水利枢纽工程设计[J].水利规划与设计,2006.05
[56]梁宗仁.甘肃九甸峡水利枢纽深厚覆盖层工程特性[J].水利规划与设计,2007,03
[57]蒋国澄.中国混凝土面板堆石坝20年[M].北京：中国水利水电出版社,2005：1-25
[58]吕生玺,沈振中,温续余.九甸峡混凝土面板堆石坝地震特性研究[J].中国农村水利水电,2008,02