基于Matlab和GPU的BESO方法的全流程并行计算策略

　　针对传统并行计算方法实现结构拓扑优化快速计算的硬件成本高、程序开发效率低的问题，提出了一种基于Matlab和图形处理器（GPU）的双向渐进结构优化（BESO）方法的全流程并行计算策略。首先，探讨了Matlab编程环境中实现GPU并行计算的三种途径的优缺点和适用范围；其次，分别采用内置函数直接并行的方式实现了拓扑优化算法中向量和稠密矩阵的并行化计算，采用MEX函数调用CUSOLVER库的形式实现了稀疏格式有限元方程组的快速求解，采用并行线程执行（ PTX）代码的方式实现了拓扑优化中单元敏度分析等优化决策的并行化计算。数值算例表明，基于Matlab直接开发GPU并行计算程序不仅编程效率高，而且还可以避免不同编程语言间的计算精度差异，最终使CPU并行程序可以在保持计算结果不变的前提下取得可观的加速比。