供水系统震后水力分析算法并行化研究  

Research on parallel algorithm for post-earthquake hydraulic analysis of water supply system

作　　者：龙立[1,2] 贺金川 郑山锁 周炎[1,2] LONG Li;HE Jinchuan;ZHENG Shansuo;ZHOU Yan(School of Civil Engineering,Xi’an University of Architecture and Technology,Xi’an 710055,China;Key Laboratory of Structural Engineering and Earthquake Resistance,Ministry of Education,Xi’an University of Architecture and Technology,Xi’an 710055,China;Institute of Architectural Design and Research,Xi’an University of Architecture and Technology,Xi’an 710055,China)

机构地区：[1]西安建筑科技大学土木工程学院,陕西西安710055 [2]西安建筑科技大学结构工程与抗震教育部重点实验室,陕西西安710055 [3]西安建筑科技大学建筑设计研究院,陕西西安710055

出　　处：《华中科技大学学报（自然科学版）》2020年第12期121-126,共6页Journal of Huazhong University of Science and Technology(Natural Science Edition)

基　　金：国家重点研发计划资助项目(2019YFC1509302);国家自然科学基金资助项目(51678475);西安市科技计划项目(2019113813CXSF016SF026);陕西省教育厅产业化项目(18JC020)。

摘　　要：为了适应大规模供水系统震后水力分析的需求,提出了基于中央处理单元-图形处理单元(CPU-GPU)异构计算的供水系统震后水力分析算法.该算法结合对称逐次超松弛-近似逆预处理(SSOR-AI)和共轭梯度法(CG)求解线性方程组,其中CPU负责供水管网系统节点及管段的基本数据输入、迭代初值设定、收敛性判断等步骤,GPU负责迭代雅可比矩阵、节点流量闭合差矢量生成及线性方程组求解.算例表明该算法与传统牛顿-拉夫森方法计算所得的节点自由水压及节点流量最大误差分别为0.44%和0.51%,验证了该方法的精确性.该算法预处理效果良好,当供水系统规模较大时,所生成线性方程组的时间较串行方法最高可获得13倍加速效果;整体水力分析效率方面,与基于Cholesky分解的直接法相比可获得6倍加速效果,能够满足大规模供水系统震后水力分析的需求.In order to meet the demand of hydraulic analysis of large-scale water supply system after earthquake,a post-earthquake hydraulic analysis algorithm for water supply system based on central processing unit-graphics processing unit(CPU-GPU) heterogeneous computing was proposed.Symmetric successive over relaxation approximate inverse(SSOR-AI) preconditioner and conjugate gradient(CG) algorithm were combined by this algorithm to solve linear equations.The steps of basic data input of water supply network system nodes and pipe sections,iterative initial value set and convergence judgement were processed by CPU,while the generation of iterative Jacobian matrix and node flow closure difference vector and the solution of linear equations were processed by GPU.The case analysis shows that the maximum error of the node free water pressure and node flow calculated by the proposed algorithm and the traditional Newton-Raphson method is 0.44% and 0.51%,respectively,which verifies the accuracy of the proposed algorithm.And the proposed algorithm has good preprocessing effect.When the scale of the water supply system is large,the proposed algorithm can generate linear equations up to 13 times faster than the serial method.In terms of overall hydraulic analysis efficiency,it can obtain 6 times acceleration effect compared with the direct method based on Cholesky decomposition,which can meet the requirement of large-scale water supply system post-earthquake hydraulic analysis.

关 键 词：水管道系统 震后水力分析 并行计算 牛顿-拉夫森法 共轭梯度法 

分 类 号：TU990.3[建筑科学—市政工程]