论文部分内容阅读
随着高层建筑结构向精细化建模发展,精细化的高层三维结构分析模型自由度数量巨大,计算工作量庞大,对软件分析平台的计算速度要求更高。然而目前国内外分析软件多基于单一CPU串行平台的构建方式,计算耗时成为突出问题,如何实现高性能低价格的有限元分析平台,成为迫切需要解决的问题。采用宏观有限元模型进行非线性分析能有效减小计算工作量、缩短计算时间,然而宏观有限元模型无法完整描述三维结构非线性特性,可能导致分析结果存在较大误差。为寻求能够兼顾计算精度与计算效率二者平衡的精细化模型分析软件平台,本文通过CPU串行计算与GPU高性能并行计算相结合,建立CPU-GPU异构平台,开发适用于GPU的非线性有限元并行算法程序和软件,实现高层结构弹塑性时程分析的高精度、高效率。论文主要研究工作如下:(1)CPU-GPU异构平台建立与验证。基于CPU和GPU的混合编程思路,让CPU处理数据的读写和逻辑控制等顺序型任务,GPU则处理计算密集的大规模数据操作,建立了与有限元“分而治之”并行思想统一的CPU-GPU异构平台。基于CUDATM并行编程模型分别开发了适用于GPU加速运算的共轭梯度法(CG)和预处理共轭梯度法(PCG)并行求解器,实现有限元大型稀疏线性方程组在CPU-GPU异构平台上的高效并行求解,采用多个适用于GPU计算的并行优化策略,进一步提高了CPU-GPU异构平台的计算效率,并通过框架结构算例初步验证了异构平台的可靠性和有效性。(2)基于CPU-GPU异构平台的纤维模型非线性分析方法研究。在CPU-GPU异构平台上建立了考虑轴力-弯矩-剪力-扭矩共同作用的纤维模型,对箍筋约束混凝土补充了材料损伤试验研究,以此建立了约束混凝土的材料模型。提出了适用于GPU计算,且能有效解决结构负刚度问题的位移增量迭代算法,从有限元离散化角度出发,建立了全局自由度计算数据与GPU线程一一对应映射关系,结合纤维模型、材料模型和位移控制算法开发了基于GPU的有限元并行程序,实现了高精度高效率的梁柱结构非线性模拟。(3)基于CPU-GPU异构平台的空间壳元非线性分析方法研究。通过分层式的平面应力单元和考虑剪切变形的Mindlin板单元组合,构造了基于CPU-GPU异构平台的高精度空间壳元模型,开发了适用于模拟剪力墙和楼板的二维混凝土滞回本构模型,给出了应力状态更新的混凝土本构积分方法,提出了适用于GPU计算的Newton-Raphson迭代算法,实现了高效的板壳结构非线性有限元分析。通过与试验结果对比,表明CPU-GPU异构平台建立的空间壳元能够较为准确地描述剪力墙复杂受力状态下的非线性行为,并且计算效率较高。(4)基于CPU-GPU异构平台的复杂高层钢筋混凝土结构非线性动力时程分析方法研究。在CPU-GPU异构平台上建立了高层建筑结构非线性地震反应分析的空间结构模型,提出了包含复杂材料本构计算的精细化空间纤维模型和空间分层式壳元模型的地震响应分析GPU并行计算方法。开发了基于GPU的弹塑性地震响应分析程序,结合EBE技术将结构的“整体”计算分解到“单元”上,有效减小系统方程组求解对内存空间的需求。通过与振动台试验结果进行对比分析,验证了程序的有效性。选取2个实际高层钢筋混凝土框筒结构工程进行大震弹塑性时程分析,实例分析表明本文所建立的CPU-GPU异构平台及所开发并行程序具有求解规模大、计算精度高、计算效率高、可移植性好、适应性强等优势,具有良好的工程应用前景。