大型地下洞室三维弹塑性损伤动力有限元分析
|
摘要
采用动力时程分析法研究地下洞室地震响应特性。根据运动学动力平衡方程,以Newmark直接积分法为基础,结合三维弹塑性损伤有限元理论,编写动力时程计算程序。采用塑性附加荷载、损伤附加荷载思想考虑单元塑性及损伤情况,每时步计算对单元损伤状态进行继承,反映损伤不可逆性。采用增量变塑性刚度迭代法,求解非线性动力方程。计算结果表明,该方法一般3步内收敛,计算快,精度高。阐述大型地下洞室动力响应时程分析中,模型边界条件的设置方法及动荷载的输入方式。通过对渔子溪大型水电站地下厂房洞室群的模拟计算表明,该程序计算速度较其他商业程序速度快,计算结果合理,准确,可以应用于地下工程动力时程分析。
The seismic response characteristics of underground caverns are studied by employing dynamic time-history method. Based on the dynamic equilibrium equations of kinemato-dynamic time-history method and Newmark direct integration method,the computational codes are programmed by adopting three-dimensional elastoplastic finite element theory. The concepts of additional plastic load and additional damage load are considered to simulate the plasticity and damage of elements. The damage state of rock mass within each time step is inherited by subsequent time steps. Therefore,the irreversible trend of rock damage can be reflected. The incremental variable plastic stiffness iteration method is employed to solve dynamic equations. The calculation results show that the proposed iteration method has good convergence and fast computational speed. The iteration of each step can attain convergence within 3 steps. The way to setup boundary condition and the way to input dynamic loads are expounded by performing dynamic time-history analysis of large-scale complex underground caverns of Yuzixi hydropower station. The calculation results show that the computational speed of the proposed method in the codes is faster than that of common software. The calculation results are reasonable as well,which indicate that the proposed methods are rational and can be used to carry out dynamic time-history analysis of underground engineering.
The seismic response characteristics of underground caverns are studied by employing dynamic time-history method. Based on the dynamic equilibrium equations of kinemato-dynamic time-history method and Newmark direct integration method,the computational codes are programmed by adopting three-dimensional elastoplastic finite element theory. The concepts of additional plastic load and additional damage load are considered to simulate the plasticity and damage of elements. The damage state of rock mass within each time step is inherited by subsequent time steps. Therefore,the irreversible trend of rock damage can be reflected. The incremental variable plastic stiffness iteration method is employed to solve dynamic equations. The calculation results show that the proposed iteration method has good convergence and fast computational speed. The iteration of each step can attain convergence within 3 steps. The way to setup boundary condition and the way to input dynamic loads are expounded by performing dynamic time-history analysis of large-scale complex underground caverns of Yuzixi hydropower station. The calculation results show that the computational speed of the proposed method in the codes is faster than that of common software. The calculation results are reasonable as well,which indicate that the proposed methods are rational and can be used to carry out dynamic time-history analysis of underground engineering.
引文
[1]林皋.地下结构抗震分析综述(上)[J].世界地震工程,1990,(2):l–10.(LIN Gao.Overview of anti-seismic analysis for underground structures(I)[J]World Earthquake Engineering,1990,(2):l–10.(in Chinese))
[2]周维垣,杨强.岩石力学数值计算方法[M].北京:中国电力出版社,2005.(ZHOU Weiyuan,YANG Qiang.Numerical computational methods for rock mechanics[M].Beijing:China Electric Power Press,2005.(in Chinese))
[3]朱伯芳.有限单元法原理与应用[M].北京:中国水利水电出版社,1998.(ZHU Bofang.Finite element method theory and applications[M].Beijing:China Water Power Press,1998.(in Chinese))
[4]钱胜国,陈玲玲,吴新霞,等.爆破地震工程结构动力分析[M].北京:中国水利水电出版社,2006.(QIAN Shengguo,CHEN Lingling,WU Xinxia,et al.Project structural dynamic response analysis ofblasting ground vibration[M].Beijing:China Water Power Press,2006.(in Chinese))
[5]俞裕泰,肖明.大型地下洞室围岩稳定三维弹塑性有限元分析[J].岩石力学与工程学报,1987,6(1):47–56.(YU Yutai,XIAO Ming.Three-dimensional elastoplastic finite element analysis for surrounding rock stability of large-scale underground openings[J].Chinese Journal of Rock Mechanics and Engineering,1987,6(1):47–56.(in Chinese))
[6]FRANTZISKONIS G,DESAI C S.Constitutive model with strain softening[J].International Journal of Solids and Structures,1987,23(6):733–750.
[7]肖明.地下高压钢筋混凝土岔管渗水开裂三维数值分析计算[J].岩石力学与工程学报,2002,21(7):1022–1026.(XIAO Ming.Three-dimensional numerical analysis of lining seepage crack of underground concrete branch pipe with high water pressure[J].Chinese Journal of Rock Mechanics and Engineering,2002,21(7):1022–1026.(in Chinese))
[8]王全凤,黄庆丰.损伤结构弹塑性抗震性能分析及评价[J].华侨大学学报(自然科学版),2009,30(1):6–11.(WANG Quanfeng,HUANG Qingfeng.Elastoplastic performance of damaged structures and its evaluation[J].Journal of Huaqiao University(Natural Science),2009,30(1):6–11.(in Chinese))
[9]黄胜,陈卫忠,杨建平,等.地下工程地震动力响应及抗震研究[J].岩石力学与工程学报,2009,28(3):483–490.(HUANG Sheng,CHEN Weizhong,YANG Jianping,et al.Research on earthquake-induced dynamic responses and aseismic measures for underground engineering[J].Chinese Journal of Rock Mechanics and Engineering,2009,28(3):483–490.(in Chinese))
[10]LYSMER J,KULEMEYER R L.Finite dynamic model for infinite media[J].Journal of the Engineering Mechanics Division,ASCE,1969,95(4):859–877.
[11]廖振鹏,黄孔亮,杨柏坡,等.暂态波透射边界[J].中国科学(A辑),1984,(6):556–564.(LIAO Zhenpeng,HUANG Kongliang,YANG Baipo,et al.Transmitting boundary of transient wave[J].Science in China(Series A),1984,(6):556–564.(in Chinese))
[12]吴鸿庆,任侠.结构有限元分析[M].北京:中国铁道出版社,2000.(WU Hongqing,REN Xia.Finite element analysis of structures[M].Beijing:China Railway Publishing House,2000.(in Chinese))
[13]Hibbitt,Karlsson and Sorensen,Inc..ABAQUS theory manual and analysis user′s manual[R].Pawtucket,USA:Hibbitt,Karlsson and Sorensen,Inc.,2002.
[14]Itasca Consulting Group,Inc..FLAC3D users manual(version3.0)[R].Minneapolis,USA:Itasca Consulting Group,Inc.,2005.
[15]中华人民共和国国家标准编写组.GB17741–2005工程场地地震安全性评价技术规范[S].北京:中国标准出版社,2005.(The National Standards Compilation Group of People′s Republic of China.GB17741–2005Evaluation of seismic safety for engineering sites[S].Beijing:Standards Press of China,2005.(in Chinese))
[16]华东水利学院土力学教研室.土工原理与计算[M].北京:水利电力出版社,1984.(Soil Mechanics Teaching and Research Group,East China Institute of Water Conservancy.Soil mechanics theory and calculation[M].Beijing:Water Resources and Electric Power Press,1984.(in Chinese))
[2]周维垣,杨强.岩石力学数值计算方法[M].北京:中国电力出版社,2005.(ZHOU Weiyuan,YANG Qiang.Numerical computational methods for rock mechanics[M].Beijing:China Electric Power Press,2005.(in Chinese))
[3]朱伯芳.有限单元法原理与应用[M].北京:中国水利水电出版社,1998.(ZHU Bofang.Finite element method theory and applications[M].Beijing:China Water Power Press,1998.(in Chinese))
[4]钱胜国,陈玲玲,吴新霞,等.爆破地震工程结构动力分析[M].北京:中国水利水电出版社,2006.(QIAN Shengguo,CHEN Lingling,WU Xinxia,et al.Project structural dynamic response analysis ofblasting ground vibration[M].Beijing:China Water Power Press,2006.(in Chinese))
[5]俞裕泰,肖明.大型地下洞室围岩稳定三维弹塑性有限元分析[J].岩石力学与工程学报,1987,6(1):47–56.(YU Yutai,XIAO Ming.Three-dimensional elastoplastic finite element analysis for surrounding rock stability of large-scale underground openings[J].Chinese Journal of Rock Mechanics and Engineering,1987,6(1):47–56.(in Chinese))
[6]FRANTZISKONIS G,DESAI C S.Constitutive model with strain softening[J].International Journal of Solids and Structures,1987,23(6):733–750.
[7]肖明.地下高压钢筋混凝土岔管渗水开裂三维数值分析计算[J].岩石力学与工程学报,2002,21(7):1022–1026.(XIAO Ming.Three-dimensional numerical analysis of lining seepage crack of underground concrete branch pipe with high water pressure[J].Chinese Journal of Rock Mechanics and Engineering,2002,21(7):1022–1026.(in Chinese))
[8]王全凤,黄庆丰.损伤结构弹塑性抗震性能分析及评价[J].华侨大学学报(自然科学版),2009,30(1):6–11.(WANG Quanfeng,HUANG Qingfeng.Elastoplastic performance of damaged structures and its evaluation[J].Journal of Huaqiao University(Natural Science),2009,30(1):6–11.(in Chinese))
[9]黄胜,陈卫忠,杨建平,等.地下工程地震动力响应及抗震研究[J].岩石力学与工程学报,2009,28(3):483–490.(HUANG Sheng,CHEN Weizhong,YANG Jianping,et al.Research on earthquake-induced dynamic responses and aseismic measures for underground engineering[J].Chinese Journal of Rock Mechanics and Engineering,2009,28(3):483–490.(in Chinese))
[10]LYSMER J,KULEMEYER R L.Finite dynamic model for infinite media[J].Journal of the Engineering Mechanics Division,ASCE,1969,95(4):859–877.
[11]廖振鹏,黄孔亮,杨柏坡,等.暂态波透射边界[J].中国科学(A辑),1984,(6):556–564.(LIAO Zhenpeng,HUANG Kongliang,YANG Baipo,et al.Transmitting boundary of transient wave[J].Science in China(Series A),1984,(6):556–564.(in Chinese))
[12]吴鸿庆,任侠.结构有限元分析[M].北京:中国铁道出版社,2000.(WU Hongqing,REN Xia.Finite element analysis of structures[M].Beijing:China Railway Publishing House,2000.(in Chinese))
[13]Hibbitt,Karlsson and Sorensen,Inc..ABAQUS theory manual and analysis user′s manual[R].Pawtucket,USA:Hibbitt,Karlsson and Sorensen,Inc.,2002.
[14]Itasca Consulting Group,Inc..FLAC3D users manual(version3.0)[R].Minneapolis,USA:Itasca Consulting Group,Inc.,2005.
[15]中华人民共和国国家标准编写组.GB17741–2005工程场地地震安全性评价技术规范[S].北京:中国标准出版社,2005.(The National Standards Compilation Group of People′s Republic of China.GB17741–2005Evaluation of seismic safety for engineering sites[S].Beijing:Standards Press of China,2005.(in Chinese))
[16]华东水利学院土力学教研室.土工原理与计算[M].北京:水利电力出版社,1984.(Soil Mechanics Teaching and Research Group,East China Institute of Water Conservancy.Soil mechanics theory and calculation[M].Beijing:Water Resources and Electric Power Press,1984.(in Chinese))
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心