沉管隧道地震响应分析的并行计算
|
摘要
采用全三维非线性建模方法建立了沉管隧道总体有限元模型,模型规模上节点总数超过100万,单元总数超过120万。由于系统的非线性,以及地震激励是一种时间历程,沉管隧道地震响应数值模拟的实质是一个超大规模非线性系统的瞬态响应计算问题,因此导致目前在地震安全性评价领域普遍采用的串行算法、串行软件和普通的计算机无法胜任此工作。根据所用超级计算机的体系结构特点,采用并行区域分解算法,通过对沉管隧道总体模型合理分区来进行地震响应分析。地震波载荷选用调幅后的唐山地震波,应用通用的非线性动力分析有限元程序LS-DYNA进行求解,计算结果可为隧道的抗震设计提供有力依据。
Using the high performance computer SGI Onyx3800,this paper introduced finite element method in simulating the seismic response of the immersed tunnel in Shanghai. Consisting the immersed tunnel,soil and river,a large three-dimensional finite element model was finished in the CAE software;and the number of nodes and elements in final finite element model exceeds 1.2 million and 1 million respectively. The seismic response numerical simulation for immersed tunnel is a transient response analysis of large-scale system. Thus,it is impossible to finish this simulation by serial algorithm,serial software and PC. The ultimate calculation adopted explicit dynamic algorithm and parallel domain decomposition method. The recursive coordinate bisection method was used and improved considering the construction feature of immersed tunnel. The result reveals the weak joints in this immersed tunnel under Tangshan seismic waves. It can also provide data and references to the aseismatic design of immersed tunnel and its flexible joints.
Using the high performance computer SGI Onyx3800,this paper introduced finite element method in simulating the seismic response of the immersed tunnel in Shanghai. Consisting the immersed tunnel,soil and river,a large three-dimensional finite element model was finished in the CAE software;and the number of nodes and elements in final finite element model exceeds 1.2 million and 1 million respectively. The seismic response numerical simulation for immersed tunnel is a transient response analysis of large-scale system. Thus,it is impossible to finish this simulation by serial algorithm,serial software and PC. The ultimate calculation adopted explicit dynamic algorithm and parallel domain decomposition method. The recursive coordinate bisection method was used and improved considering the construction feature of immersed tunnel. The result reveals the weak joints in this immersed tunnel under Tangshan seismic waves. It can also provide data and references to the aseismatic design of immersed tunnel and its flexible joints.
引文
[1] Kiyomiya O. Earthquake-resistant design features of immersed tunnels in Japan[J]. Tunneling and Underground Space Technology,1995,20(4):463–475.
[2]严松宏,高峰,李德武,等.南京长江沉管隧道的地震安全性评价[J].岩石力学与工程学报,2003,22(增2):2800–2803.(Yan Songhong,Gao Feng,Li Dewu,et al. Estimation on seismic safety of Nanjing Changjiang submerged tunnel[J]. Chinese Journal of Rock Mechanics and Engineering,2003,22(Supp.2):2800–2803.(in Chinese))
[3]杨辉,陆建飞,王建华,等.黄浦江越江隧道的动力抗震分析[J].上海交通大学学报,2001,35(10):1512–1516.(Yang Hui,Lu Jianfei,Wang Jianhua,et al. Vibration analysis of a subaqueous tunnel of Huangpu River during earthquake[J]. Journal of Shanghai Jiaotong University [J],2001,35(10):1512–1516.(in Chinese))
[4]周树荃,邓绍忠.有限元结构分析并行计算的若干研究进展[J].南京航空航天大学学报,1995,27(1):27–32.(Zhou Shuquan,Deng Shaozhong. Some recent advances in the parallel computation of finite element structural analysis[J]. Journal of Nanjing University of Aeronautics and Astronautics,1995,27(1):27–32.(in Chinese))
[5] Bitzarakis B , Papadrakakis M , Kotopulos A. Parallel solution techniques in computational structural mechanics[J]. Computer Methods in Applied Mechanics and Engineering,1997,148(1):75–104.
[6]舒继武,郑纬民,汪东升. MPP 中区域分解法的临界子区域数的确定方法[J].清华大学学报(自然科学版),2000,40(7):105–108.(Shu Jiwu,Zheng Weimin,Wang Dongsheng,et al. Selection of the ideal number of sub-domains for domain decomposition in MPP [J]. Journal of Tsinghua University(Sci. and Tec. Edition.),2000,40(7):105–108.(in Chinese))
[7]郑哲明,程建钢,姚振汉.有限元并行计算中网格自动区域划分的研究[J].工程力学,2002,19(6):54–57.(Zheng Zheming,Cheng Jiangang,Yao Zhenhan. On automatic mesh generation in parallel finite element analysis[J]. Engineering Mechanics,2002,19(6):54–57.(in Chinese))
[8] Ach K,Alart P,Barboteu M,et al. Parallel frictional contact algorithms and industrial applications[J]. Computer Methods in Applied Mechanics and Engineering,1999,177:169–181.
[9] Steve P,Steve A,Bruce H,et al. Parallel tran-sient dynamics simulations:algorithms for contact detection and smoothed particle hydrodynamics[J]. Journal of Parallel and Distributed Computing,1998,50:104–122.
[10] Kirzhner F,Rosenhouse G . Numerical analysis of tunnel dynamic response to earth motions[J]. Tunneling and Underground Space Technology,2000,15(3):249–258.接头相对位移/mm
[2]严松宏,高峰,李德武,等.南京长江沉管隧道的地震安全性评价[J].岩石力学与工程学报,2003,22(增2):2800–2803.(Yan Songhong,Gao Feng,Li Dewu,et al. Estimation on seismic safety of Nanjing Changjiang submerged tunnel[J]. Chinese Journal of Rock Mechanics and Engineering,2003,22(Supp.2):2800–2803.(in Chinese))
[3]杨辉,陆建飞,王建华,等.黄浦江越江隧道的动力抗震分析[J].上海交通大学学报,2001,35(10):1512–1516.(Yang Hui,Lu Jianfei,Wang Jianhua,et al. Vibration analysis of a subaqueous tunnel of Huangpu River during earthquake[J]. Journal of Shanghai Jiaotong University [J],2001,35(10):1512–1516.(in Chinese))
[4]周树荃,邓绍忠.有限元结构分析并行计算的若干研究进展[J].南京航空航天大学学报,1995,27(1):27–32.(Zhou Shuquan,Deng Shaozhong. Some recent advances in the parallel computation of finite element structural analysis[J]. Journal of Nanjing University of Aeronautics and Astronautics,1995,27(1):27–32.(in Chinese))
[5] Bitzarakis B , Papadrakakis M , Kotopulos A. Parallel solution techniques in computational structural mechanics[J]. Computer Methods in Applied Mechanics and Engineering,1997,148(1):75–104.
[6]舒继武,郑纬民,汪东升. MPP 中区域分解法的临界子区域数的确定方法[J].清华大学学报(自然科学版),2000,40(7):105–108.(Shu Jiwu,Zheng Weimin,Wang Dongsheng,et al. Selection of the ideal number of sub-domains for domain decomposition in MPP [J]. Journal of Tsinghua University(Sci. and Tec. Edition.),2000,40(7):105–108.(in Chinese))
[7]郑哲明,程建钢,姚振汉.有限元并行计算中网格自动区域划分的研究[J].工程力学,2002,19(6):54–57.(Zheng Zheming,Cheng Jiangang,Yao Zhenhan. On automatic mesh generation in parallel finite element analysis[J]. Engineering Mechanics,2002,19(6):54–57.(in Chinese))
[8] Ach K,Alart P,Barboteu M,et al. Parallel frictional contact algorithms and industrial applications[J]. Computer Methods in Applied Mechanics and Engineering,1999,177:169–181.
[9] Steve P,Steve A,Bruce H,et al. Parallel tran-sient dynamics simulations:algorithms for contact detection and smoothed particle hydrodynamics[J]. Journal of Parallel and Distributed Computing,1998,50:104–122.
[10] Kirzhner F,Rosenhouse G . Numerical analysis of tunnel dynamic response to earth motions[J]. Tunneling and Underground Space Technology,2000,15(3):249–258.接头相对位移/mm
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心