面向桥梁抗震仿真的高性能并行数值方法

High Performance Parallel Numerical Method for Anti-Seismic Simulation of Bridges

为全面了解桥梁在地震环境中的动态结构响应问题,提出了一套高性能数值模拟方法.该方法以弹性模型为基础,在非结构化四面体网格上应用有限元技术,在时间上采用向后差分方法,构造一个全隐式的离散格式.为精确高效求解离散得到的大规模线性系统,采用Krylov子空间方法,同时通过重叠区域分解的思想来构造预条件算子以加快算法收敛速度.在数值实验部分,利用其对一个典型的钢管混凝土拱桥结构进行了数值仿真模拟,得到了该桥梁在地震作用下不同时刻各部分的位移、速度以及应变,发现桥梁在地震作用下最容易发生形变的是桥拱,最后还验证了该方法的并行可扩展性.

A set of high-performance numerical simulation method is proposed to fully find out about the dynamic structure response of bridges in seismic environment.Based on the elastic model,the finite element method is applied to unstructured tetrahedral mesh,and the backward difference method is used to construct a fully implicit discrete scheme in time.In order to solve the discretely obtained large-scale linear system with accuracy and efficiency,Krylov subspace method is adopted.At the same time,the preconditioner operator is constructed by the idea of overlapping domain decomposition to accelerate the convergence speed of the algorithm.In the part of numerical experiment,by means of a numerical simulation of a typical concrete filled steel tubular arch bridges structure,the displacement,velocity and response of each part of the bridge at different times under the action of earthquake are obtained.It is found that the bridge arch is the most prone to deformation under the action of earthquake.Finally,the parallel scalability of the method is verified.