一种非线性显式分层壳单元及其GPU并行计算实现  

Nonlinear explicit layered shell element and its GPU parallel computing implementation

作　　者：曹胜涛 李志山[2] 杨志勇[1] CAO Shengtao;LI Zhishan;YANG Zhiyong(China Academy of Building Research,Beijing 100013,China;Guangdong Provincial Key Laboratory of Earthquake Engineering and Applied Technology,Guangzhou University,Guangzhou 510405,China)

机构地区：[1]中国建筑科学研究院,北京100013 [2]广州大学广东省地震工程与应用技术重点实验室,广州510405

出　　处：《振动与冲击》2019年第22期60-69,84,共11页Journal of Vibration and Shock

基　　金：十三五高性能计算重点专项(2016YFB0200605)

摘　　要：通用有限元程序ABAQUS的钢筋混凝土显式分层壳单元被广泛应用于剪力墙抗震性能分析,但存在两个缺陷:①只能得到混凝土受压损伤和受拉损伤,无法反映混凝土剪切损坏,因此不易根据损伤类型对结构进行优化;②基于CPU并行计算,大规模计算效率较低。基于平面应力条件下的混凝土弹塑性损伤本构模型,根据混凝土损伤发展时的受力状态和工程实践需要将损伤分为受拉损伤、受压损伤和受剪损伤。结合可损伤分类的塑性损伤模型,给出了非线性壳元物理沙漏力和面内旋转力的构造方法,进而得到一种含面内旋转自由度的4节点24自由度四边形非线性显式分层壳单元。将该研究壳元在自主研发的基于CPU+GPU异构并行计算的非线性分析软件中完成开发实现;通过与ABAQUS benchmark算例结果对比,验证了开发内容的正确性;通过与剪力墙单调加载试验对比,验证了该研究壳元的合理性。分别采用自主研发软件和ABAQUS对上海地区某框架核心筒体系的超高层结构进行了罕遇地震非线性时程分析,结果表明:①自主研发软件与ABAQUS结果基本一致,而计算效率为ABAQUS计算效率的5.69倍;②自主研发软件得到的受拉损伤、受压损伤和受剪损伤损伤可更清晰地揭示核心筒在罕遇地震作用下的损坏演化规律和破坏模式。Based on an elastoplastic damage model of concrete under plane stress,the damage of concrete was classified into three types,tension,compression and shear,according to the concrete stress state and the demands of engineering practice,and a method for the computation of nonlinear physical hourglass force and drilling force was given.Then,a quadrilateral nonlinear explicit layered shell element with 4 nodes-24 degrees of freedom including in-plane drilling was proposed.The proposed shell element was implemented in a nonlinear analysis software which was independently-developed and based on CPU+GPU parallel computing.By comparing the results of benchmark example by the ABAQUS and the software,the correctness of the shell element was verified.The rationality of the shell element was also proved by the monotonic loading experimental results on a shear wall.In addition,the seismic nonlinear time history analysis of a high-rise frame-tube structure in Shanghai was carried out by using the software and ABAQUS respectively.The results show that:the results by the software and ABAQUS are basically the same,and the computational efficiency of the software is 5.69 times that of the ABAQUS,moreover,the tension damage,compression damage and shear damage shown by the software are more useful to reveal the damage evolution and failure of the tube.

关 键 词：弹塑性损伤模型 混凝土损伤类型 非线性显式分层壳单元 剪力墙损坏模式 CPU+GPU并行计算 

分 类 号：TU311[建筑科学—结构工程]