工程力学 ›› 2018, Vol. 35 ›› Issue (11): 79-85,91.doi: 10.6052/j.issn.1000-4750.2017.08.0656

• 土木工程学科 • 上一篇    下一篇

图形处理器加速算法在复杂高层结构非线性响应分析中的应用

李红豫1,2, 滕军3, 李祚华3, 张璐4   

  1. 1. 桂林理工大学土木与建筑工程学院, 桂林 541004;
    2. 桂林理工大学广西有色金属隐伏矿床勘查及材料开发协同创新中心, 桂林 541004;
    3. 哈尔滨工业大学深圳研究生院深圳市城市与土木工程防灾减灾重点实验室, 深圳 518055;
    4. 伊利诺伊大学芝加哥分校土木与材料工程系, 美国芝加哥 60607
  • 收稿日期:2017-01-11 修回日期:2017-08-30 出版日期:2018-11-07 发布日期:2018-11-07
  • 通讯作者: 李红豫(1985-),女,广西人,讲师,博士,主要从事高层建筑结构抗震和并行工程计算研究(E-mail:lihongyu@glut.edu.cn). E-mail:lihongyu@glut.edu.cn
  • 作者简介:滕军(1962-),男,辽宁人,教授,博士,博导,主要从事结构抗震抗风和结构振动控制方面研究(E-mail:tengj@hit.edu.cn);李祚华(1978-),男,吉林人,副教授,博士,主要从事高层建筑结构抗震和结构无损检测研究(E-mail:lizuohua@hitsz.edu.cn);张璐(1983-),男,吉林人,博士后,主要从事结构非线性分析和结构无损检测研究(E-mail:zhang889@uic.edu).
  • 基金资助:
    国家自然科学基金青年科学基金项目(51708147);国家自然科学基金重大国际(中美)合作研究项目(51261120374);广西自然科学基金青年科学基金项目(2017GXNSFBA198184);桂林理工大学科研启动基金项目(GLUTQD2016046);广西高等学校高水平创新团队及卓越学者计划项目(桂教人[2017]38号)

APPLICATION OF GRAPHICS PROCESSING UNIT BASED ALGORITHM IN NONLINEAR RESPONSE ANALYSIS TO COMPLEX HIGH-RISE BUILDING STRUCTURES

LI Hong-yu1,2, TENG Jun3, LI Zuo-hua3, ZHANG Lu4   

  1. 1. College of Civil Engineering and Architecture, Guilin University of Technology, Guilin 541004, China;
    2. Collaborative Innovation Center for Exploration of Hidden Nonferrous Metal Deposits and Development of New Materials in Guangxi, Guilin University of Technology, Guilin 541004, China;
    3. Shenzhen Key Laboratory of Urban & Civil Engineering Disaster Prevention & Reduction, Shenzhen Graduate School, Harbin Institute of Technology, Shenzhen 518055, China;
    4. Department of Civil and Materials Engineering, University of Illinois at Chicago, Chicago, IL 60607, USA
  • Received:2017-01-11 Revised:2017-08-30 Online:2018-11-07 Published:2018-11-07

摘要: 目前有限元分析软件多基于中央处理器的平台方式构建,在处理复杂高层结构非线性响应分析时暴露出计算耗时多、计算效率低以及对计算硬件要求高等问题。图形处理器由于其硬件构造的先天优势,可以提供十倍乃至上百倍于中央处理器的浮点运算和并行计算性能,因而为高层结构非线性计算所面临的瓶颈问题提供了一个切实可行的解决方法。该文在构建异构并行计算平台的基础上,提出一种适用于图形处理器加速的有限元并行数值计算方法。该方法利用精细化结构分析模型的自由度数据和图形处理器中的线程建立一一对应映射关系,对动力响应的隐式积分算法进行图形处理器线程级的并行化处理,并且结合EBE单元级的优化存储空间机制,降低系统方程组求解时对内存空间的需求。通过对比振动台试验结果对该方法进行验证,并对实际高层钢筋混凝土框筒结构工程进行弹塑性地震响应分析,结果显示该文所提方法在保证模型精度前提下能有效提高大型复杂高层结构非线性响应分析效率。

关键词: 非线性响应分析, 图形处理器, 高层结构, 并行计算, EBE

Abstract: Currently, most of the commercial finite element (FE) softwares are based on the CPU architectures, which causes massively time consuming, low efficiency, and rigor of requirements of hardware during analyzing the nonlinear response of high-rise structures. Meanwhile, the emergence of GPU based algorithms presents significantly superior performance in floating-point operation and parallel computation due to its special configuration. Therefore, GPU based algorithms can provide a feasible solution for the perplexing issues of nonlinear computation of high-rise structures. Our work is to develop a parallel FE algorithm by introducing GPU and to construct a corresponding heterogeneous platform, ultimately leading to speed up the computation. Firstly, the mapping between the degrees of freedom (DOFs) of a refined model and the threads of GPU is formed. Then, the implicit integration algorithm for solving the dynamic response will be parallelized in threads; meanwhile, the strategies of storage are optimized in terms of element-by-element scale and the demand of memory was reduced while solving the equations. All of the GPU based algorithms have been validated by comparing with the experimental results of a shaking table. Moreover, the validated algorithms are extended to apply to the analysis of the elastic-plastic seismic response of a practical high-rise reinforced concrete frame tube structure. The results show that the proposed algorithm can not only guarantee the accuracy but also improve efficiency dramatically in the procedure of structural nonlinear response analyses.

Key words: nonlinear response analysis, graphics processing unit (GPU), high-rise building structures, parallel computing, EBE

中图分类号: 

  • TU973+.2
[1] 李云贵. 工程结构设计中的高性能计算[J]. 建筑结构学报, 2010, 31(6):89-95. Li Yungui. High-performance computing in structural design[J]. Journal of Building Structures, 2010, 31(6):89-95. (in Chinese)
[2] He Kai, Tan Sheldon X-D, Wang Hai, et al. GPU-Accelerated parallel sparse LU factorization method for fast circuit analysis[J]. IEEE Transactions on Very Large Scale Integration (VLSI) Systems, 2016, 24(3):1140-1150.
[3] Yang D, Peterson G D, Li H. Compressed sensing and Cholesky decomposition on FPGAs and GPUs[J]. Parallel Computing, 2012, 38(8):421-437.
[4] Naumov M. Parallel solution of sparse triangular linear systems in the preconditioned iterative methods on the GPU[R]//NVIDIA Corporation, Westford, MA, USA, Technology Report NVR-2011-001, 2011.
[5] Li R, Saad Y. GPU-accelerated preconditioned iterative linear solvers[J]. The Journal of Supercomputing, 2013, 63(2):443-466.
[6] 高家全, 王志超. 基于GPU的SSOR稀疏近似逆预条件研究[J]. 浙江工业大学学报, 2016, 44(2):140-145. Gao Jiaquan, Wang Zhichao. Research of the SSOR sparse approximate inverse preconditioner on GPUs[J]. Journal of Zhejiang University of Technology, 2016, 44(2):140-145. (in Chinese)
[7] 陈曦, 王冬勇, 任俊, 等. CPU-GPU混合计算构架在岩土工程有限元分析中的应用[J]. 土木工程学报, 2016, 49(6):105-112. Chen Xi, Wang Dongyong, Ren Jun, et al. Application of hybrid CPU-GPU computing platform in large-scale geotechnical finite element analysis[J]. China Civil Engineering Journal, 2016, 49(6):105-112. (in Chinese)
[8] Serban R, Melanz D, Li A, et al. A GPU-based preconditioned Newton-Krylov solver for flexible multibody dynamics[J]. International Journal for Numerical Methods in Engineering, 2015, 102(9):1585-1604.
[9] Fazanaro F I, Soriano D C, Suyama R, et al. Numerical characterization of nonlinear dynamical systems using parallel computing:The role of GPUs approach[J]. Communications in Nonlinear Science and Numerical Simulation, 2016, 37:143-162.
[10] Mafi R, Sirouspour S. GPU-based acceleration of computations in nonlinear finite element deformation analysis[J]. International Journal for Numerical Methods in Biomedical Engineering, 2014, 30(3):365-381.
[11] Yang Y, Yang C, Hsieh T. GPU parallelization of an object-oriented nonlinear dynamic structural analysis platform[J]. Simulation Modelling Practice and Theory, 2014, 40(1):112-121.
[12] Cai Y, Li G, Wang H, et al. Development of parallel explicit finite element sheet forming simulation system based on GPU architecture[J]. Advances in Engineering Software, 2012, 45(1):370-379.
[13] 刘小虎, 胡耀国, 符伟. 大规模有限元系统的GPU加速计算研究[J]. 计算力学学报, 2012, 29(1):146-152. Liu Xiaohu, Hu Yaoguo, Fu Wei. Solving large finite element system by GPU computation[J]. Chinese Journal of Computational Mechanics, 2012, 29(1):146-152. (in Chinese)
[14] 蔡勇, 李光耀, 王琥. 基于多重网格法和GPU并行计算的大规模壳结构快速计算方法[J]. 工程力学, 2014, 31(5):20-26. Cai Yong, Li Guangyao, Wang Hu. A fast calculation method for large-scale shell structure based on multigrid method and GPU parallel computing[J]. Engineering Mechanics, 2014, 31(5):20-26. (in Chinese)
[15] DS Simulia Corporation. Abaqus 6.12 analysis user's manual[Z]. 2012.
[16] 李红豫, 滕军, 李祚华. 钢筋混凝土框架结构非线性静、动力分析的高效计算平台HSNAS(GPU)——I程序开发[J]. 振动与冲击, 2016, 35(14):47-53. Li Hongyu, Teng Jun, Li Zuohua. An efficient platform HSNAS (GPU) for nonlinear static and dynamic analysis of reinforced concrete frames-I. Program development[J]. Journal of Vibration and Shock, 2016, 35(14):47-53. (in Chinese)
[17] Sanders J, Kandrot E. CUDA by example:an introduction to general-purpose GPU programming[M]. Addison-Wesley Professional, 2010.
[18] 孙宝印, 古泉, 张沛洲, 等. 钢筋混凝土框架结构弹塑性数值子结构分析方法[J]. 工程力学, 2016, 33(5):44-49. Sun Baoyin, Gu Quan, Zhang Peizhou, et al. Elastoplastic numerical substructure method of reinforced concrete frame structures[J]. Engineering Mechanics, 2016, 33(5):44-49. (in Chinese)
[19] 陶慕轩, 聂建国. 组合构件纤维模型的建模策略-单元划分和截面离散[J]. 工程力学, 2016, 33(2):96-103. Tao Muxuan, Nie Jianguo. Modeling strategies of fiber models for composite structural members:Element mesh and section discretization[J]. Engineering Mechanics, 2016, 33(2):96-103. (in Chinese)
[20] 吕伟荣, 罗雯, 蒋庆, 等. 基于分层壳单元的现浇混凝土空心楼板数值模拟[J]. 工程力学, 2015, 32(S1):172-183. Lv Weirong, Luo Wen, Jiang Qing, et al. Numerical simulation of cast-in-situ concrete hollow slabs based on layered shell elements[J]. Engineering Mechanics, 2015, 32(S1):172-183. (in Chinese)
[21] 李红豫, 滕军, 李祚华. 钢筋混凝土框架结构非线性静、动力分析的高效计算平台HSNAS(GPU)——Ⅱ 验证分析[J]. 振动与冲击, 2016, 35(14):54-60. Li Hongyu, Teng Jun, Li Zuohua. An efficient platform HSNAS(GPU) for nonlinear static and dynamic analysis of reinforced concrete frames-Ⅱ. Program verification and analysis[J]. Journal of Vibration and Shock, 2016, 35(14):54-60. (in Chinese)
[22] Teng Jun, Li Zuohua, Ou Jinping, et al. Fiber damage analysis model for RC beam-column based on EEP super-convergent computation[J]. Science China Technological Sciences, 2011, 54(10):2542-2548.
[1] 姜绍飞, 任晖, 骆剑彬. 基于云计算的框架结构参数并行辨识算法[J]. 工程力学, 2018, 35(4): 135-143.
[2] 曹胜涛, 李志山, 刘付钧. 一种非比例振型阻尼模型及在大规模非线性分析中的应用[J]. 工程力学, 2018, 35(11): 162-171.
[3] 王小庆, 金先龙, 杨志豪. 基于ALE的大型输水隧道地震动响应并行数值分析[J]. 工程力学, 2017, 34(3): 247-256.
[4] 唐逸豪,高振勋,蒋崇文,李椿萱. 基于LPT近似算法的CFD并行计算网格分配算法[J]. 工程力学, 2015, 32(5): 243-249,256.
[5] 陈尚鸿, 祁皑, 王素裹, 林伟. 多维多点输入下大跨度连体高层结构地震反应振动台阵试验研究[J]. 工程力学, 2014, 31(6): 212-217.
[6] 蔡勇, 李光耀, 王琥. 基于多重网格法和GPU并行计算的大规模壳结构快速计算方法[J]. 工程力学, 2014, 31(5): 20-26.
[7] 解琳琳,黄羽立,陆新征,林楷奇,叶列平. 基于OpenSees的RC框架-核心筒超高层建筑抗震弹塑性分析[J]. 工程力学, 2014, 31(1): 64-71.
[8] 董志君,李筱毅,阎培渝,欧阳蓉. 非荷载因素对超高层结构变形的影响研究[J]. 工程力学, 2013, 30(增刊): 165-168.
[9] 司林军,李国强,孙飞飞. 一种空间剪力墙非线性计算模型[J]. 工程力学, 2013, 30(7): 122-128,135.
[10] 张钰,潘鹏. 考虑子结构间相互作用的结构地震反应并行计算方法研究[J]. 工程力学, 2013, 30(5): 118-124.
[11] 刘铁林;石大川;刘 泓. 高层框架结构转动惯性对其地震响应的影响[J]. 工程力学, 2011, 28(增刊I): 110-114.
[12] 颜学渊;张永山;王焕定;魏陆顺. 三维隔震抗倾覆结构振动台试验[J]. 工程力学, 2010, 27(5): 91-096.
[13] 付朝江. 隐式非线性动力分析有限元并行求解格式[J]. 工程力学, 2010, 27(10): 27-033.
[14] 付俊峰;金 生. 用OpenMP实现三维复杂渗流场的并行计算[J]. 工程力学, 2009, 26(12): 216-221.
[15] 刘铁林;姜迎春;刘 伟. 地震作用下高层结构动力响应分析的一种数值方法[J]. 工程力学, 2008, 25(增刊Ⅰ): 0-178.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 王贞, 李强, 吴斌. 实时混合试验的自适应时滞补偿方法[J]. 工程力学, 2018, 35(9): 37 -43 .
[2] 潘天林, 吴斌. 基于桁架单元的能量一致积分方法[J]. 工程力学, 2018, 35(10): 1 -9,36 .
[3] 周大为, 陈飙松, 李云鹏, 张盛. 基于广义模态截断扩增方法的结构频响拓扑变量灵敏度分析[J]. 工程力学, 2018, 35(11): 1 -7,16 .
[4] 田建勃, 史庆轩, 刘云贺, 李慎, 马辉. PRC连梁-混合联肢剪力墙抗震性能分析[J]. 工程力学, 2018, 35(11): 53 -67 .
[5] 杨庆山, 田玉基, 陈波, 黄国庆. 行业标准《屋盖结构风荷载标准》的主要内容[J]. 工程力学, 2018, 35(7): 1 -6 .
[6] 王涛, 柳占立, 高岳, 庄茁. 基于给定参数的水力裂缝与天然裂缝相互作用结果的预测准则[J]. 工程力学, 2018, 35(11): 216 -222 .
[7] 张景, 何川, 耿萍, 陈枰良, 卢志楷. 盾构隧道环间接头弯曲状态非线性研究[J]. 工程力学, 2018, 35(11): 35 -44 .
[8] 朱世权, 李海元, 陈志华, 张辉. 攻角对空空导弹与载机分离过程的影响[J]. 工程力学, 2018, 35(9): 248 -256 .
[9] 邓明科, 杨铄, 王露. 高延性混凝土加固无筋砖墙抗震性能试验研究与承载力分析[J]. 工程力学, 2018, 35(10): 101 -111,123 .
[10] 夏阳, 廖科. 复杂三维曲梁结构的无闭锁等几何分析算法研究[J]. 工程力学, 2018, 35(11): 17 -25 .
X

近日,本刊多次接到来电,称有不法网站冒充《工程力学》杂志官网,并向投稿人收取高额费用。在此,我们郑重申明:

1.《工程力学》官方网站是本刊唯一的投稿渠道(原网站已停用),《工程力学》所有刊载论文必须经本刊官方网站的在线投稿审稿系统完成评审。我们不接受邮件投稿,也不通过任何中介或编辑收费组稿。

2.《工程力学》在稿件符合投稿条件并接收后会发出接收通知,请作者在接到版面费或审稿费通知时,仔细检查收款人是否为“《工程力学》杂志社”,千万不要汇款给任何的个人账号。请广大读者、作者相互转告,广为宣传!如有疑问,请来电咨询:010-62788648。

感谢大家多年来对《工程力学》的支持与厚爱,欢迎继续关注我们!

《工程力学》杂志社

2018年11月15日