剪力墙结构智能化生成式设计方法:从数据驱动到物理增强

2024-04-11 13:02廖文杰陆新征黄羽立赵鹏举费一凡郑哲
土木与环境工程学报 2024年1期
关键词:剪力墙工程师结构设计

廖文杰,陆新征,黄羽立,赵鹏举,费一凡,郑哲

(清华大学 土木工程安全与耐久教育部重点实验室;土木工程系,北京 100084)

人工智能的应用为土木工程领域的规划、设计、建造、维护和防灾技术带来了重大变革与重塑[1-4],智能建造已成为土木领域的重点发展方向。建筑结构智能化设计则是智能建造的重要内容,不仅可以减少大量繁琐的人工设计流程,还能为设计人员提供更加多样的初始设计选择,最终达到更好的设计效果[5-14]。建筑结构设计主要包括方案设计、优化设计及施工图设计阶段,其中方案设计对后续设计影响关键,且对设计知识与经验的需求较高,并对设计速度有较高的要求。

智能化方案设计通常采用进化算法等对结构设计方案进行搜索和设计[15-17],但是,进化算法难以有效学习结构设计数据与经验,迭代优化效率有限,难以满足快速完成结构方案设计的需求。相比而言,具备大数据特征提取与学习能力的深度生成算法快速进步,例如,基于卷积的生成神经网络、图(Graph)生成神经网络与生成对抗网络等,为建筑和结构智能化设计提供了新的可能性[18-27]。

生成对抗网络(Generative Adversarial Networks,GAN)可通过对既有设计图像和文本的学习,掌握设计数据中潜在的设计规律,实现新设计图像的生成。剪力墙结构是一种广泛应用的典型住宅建筑结构形式[28],具备设计需求较大、结构相对规则、图纸表达有效等特点,可作为智能设计的重要研究对象[12-14,25-26]。笔者基于生成对抗网络开展了从数据驱动到物理增强数据驱动的剪力墙结构智能化生成式设计方法研究,基于生成对抗网络的剪力墙结构智能化设计方法包括StructGAN[12]、txtimg2img[13](后续改称为 StructGAN-TXT)、StructGANPHY[14]。但是,目前对于数据驱动及物理增强智能设计方法的对比研究相对缺乏,不同方法中,数据特征分析与数据集创建、智能设计算法开发、测试评估方法以及典型案例分析的研究相对缺乏且不明确。笔者针对提出的3 种智能化结构设计方法开展详细的对比分析研究,明确不同方法的性能及适用范围,并对相关评价指标进行合理性分析,为未来的智能化设计方法研究提供对比模型和评价方法及指标的参考。

1 剪力墙结构智能化设计方法

基于生成对抗网络的剪力墙结构智能化设计方法如图1 所示,主要包括语义化模块、设计模块、结构建模模块[12-14]。目前,该方法主要应用于基于像素图表达的建筑和结构设计中,对于可采用像素图表达的不同结构类型的智能设计,具备一定的可拓展性和通用性。

图1 基于生成对抗网络的剪力墙结构智能化设计方法Fig.1 Intelligent structural design method for shear wall building using generative adversarial network

1)语义化模块是指将复杂的原始CAD 建筑图,通过对关键元素的提取,并采用颜色填充生成为包含建筑设计初步特征的语义化图纸。语义化过程可有效排除复杂信息的干扰,有效提升深度生成网络的设计效果,更详细的语义化方法详见既有研究[12]。该过程可通过人工提取数据特征,精度高但效率较低,也可进一步研发基于BIM(Building Information Modeling)和CAD 的数字化自动提取方法,效率高且精度可控。

2)设计模块是指采用深度生成网络,根据输入的语义化建筑图智能地生成语义化结构设计图。深度生成网络是生成对抗网络模型的生成器,在采用生成网络进行设计前,需要采用构建的建筑-结构图纸数据集对生成对抗网络进行训练,使生成网络基本掌握剪力墙结构设计能力。

3)建模模块是指基于生成的结构设计像素图建立对应的结构计算模型。通过像素图到矢量图的自动转化算法提取剪力墙构件的矢量坐标,并基于结构计算软件的应用程序接口(Application Programming Interface,API)建立对应的结构计算模型,完成整体结构的方案设计。

用于训练和测试的数据集构建方法如图2 所示,即将CAD 图纸数据进行语义化处理。在该数据集中,7 度设防且结构高度50 m 以下(简称:7d-H1)数据63 组、7 度设防且结构高度50~150 m(7d-H2)数据80 组、8 度设防(8d)数据81 组[12]。

图2 数据集构建Fig.2 Datasets construction

2 数据驱动和物理增强设计方法对比

StructGAN 采用的核心算法为图像合成生成对抗网络[12],StructGAN-TXT 采用的核心算法为文本-图像特征融合生成对抗网络[13],二者均为数据驱动方法;StructGAN-PHY 采用的核心算法为物理增强生成对抗网络[14],是一种物理增强方法。3 种方法的核心差异在于设计模块,分别为3 种不同的条件生成对抗网络模型,如图3 所示。图3(a)所示为StructGAN 模型,由图像生成器与图像真假判别器构成,基于建筑图输入生成结构设计图纸。其对应数据集的输入为建筑设计图(图2(b)),标签为结构设计图(图2(c))。训练过程中,图像生成器不断生成结构设计图,判别器不断判断生成图的真假,判别器判断正确则提升生成器,判断错误则提升判别器;生成器与判别器对抗训练不断提升,直到二者性能均难以提升达到纳什均衡,完成训练获得结构设计能力,详细训练过程可见相关研究源代码[12,29]。训练时,7d-H1、7d-H2、8d 数据集分别训练3 个StructGAN 设计模型[12]。

图3 StructGAN 核心算法对比Fig.3 Comparison on the critical algorithms of StructGAN

图3(b)所示为StructGAN-TXT 模型,基于建筑图和对应结构设计条件同时输入,生成满足建筑图和设计条件约束的结构设计图纸。相比于StructGAN,StructGAN-TXT 增加了文本输入,其对应数据集的输入为建筑设计图(图2(b))和抗震设防烈度及结构高度文本,标签为结构设计图(图2(c))。因此,其图像生成器需要同时对图像和文本特征进行提取和融合,进而基于文本-图像融合特征,生成结构设计图纸;判别器同样需要对真假图像真假文本进行判别;生成器和判别器对抗训练至性能稳定,进而具备结构设计生成能力。训练时,7d-H1、7d-H2、8d 所有图像数据混合,并对每个图像给予对应文本,创建数据集训练一个StructGANTXT 设计模型[13]。

图3(b)所示为StructGAN-PHY 模型,由图像生成器、判别器、物理评估器构成,基于建筑图输入生成器,生成满足建筑图和相应设计条件的结构设计图纸。其对应数据集的输入为建筑设计图(图2(b)),标签为结构设计图(图2(c)),输入的抗震设防烈度和结构高度并不会输入生成器中,而是辅助物理性能评估器进行力学性能评估。与StructGAN和StructGAN-TXT 相比,其网络架构中新增了物理性能评估器。物理性能评估器为基于深度判别神经网络的代理模型,通过对结构平面图纸及其对应力学性能数据集的训练,物理性能掌握对结构设计对应力学性能的预测,进而有效指导生成器的结构设计。训练时,7d-H1、7d-H2、8d 所有图像数据混合,创建多个设计条件组合,设防烈度包括6 度(0.05g)、7 度(0.10g)、7 度(0.15g)、8 度(0.20g)、8 度(0.30g),结构高度包括“<40 m”“40~60 m”“60~80 m”“80~100 m”“>100 m”,采用物理增强提升StructGAN-PHY 设计性能,得到不同设计条件对应的智能设计模型。相比而言,当StructGAN-TXT和StructGAN 没有特定设计条件下的数据时(例如6 度50 m),便无法开展对应的设计模型训练[14]。所以StructGAN-PHY 是在数据驱动学习的基础上,由物理性能评估器进行增强的智能设计方法,能有效地保证设计结果尽可能满足设计规范要求,避免由于数据数量不足或质量不佳导致的设计效果降低。

整体而言,StructGAN 与StructGAN-TXT 为数据驱动的方法,其设计能力与训练数据的数量和质量密切相关,而StructGAN-PHY 则是物理增强的方法,能通过GAN 对结构力学性能学习,有效克服数据数量或质量的影响,并保证设计结果力学性能的稳定性。

3 典型案例分析与对比

采用4 个实际工程案例进行设计,并开展性能分析、对比。需指出,由于设计资料的知识产权限制,将隐去工程项目的真实信息,以及真实的图纸和对应的结构设计模型,仅展示研究提出方法的设计结果。典型案例命名为“案例1-7d83m”,其中,“1”代表案例编号、“7d”代表设防烈度为7 度(0.1g)、“83m”代表结构高度83 m。

3.1 评价方法

图4(a)为基于图像相似性的评价方法,将生成设计的剪力墙布置与工程师的设计逐像素对比以及逐构件轮廓对比,通过式(1)~式(4)所示的交并比计算相似性,以SCV作为综合相似性指标[12-14]。图4(b)为基于结构力学模型进行力学性能分析的评价,其中像素图剪力墙构件矢量提取方法见Lu等[14]的研究。

图4 设计评价方法Fig.4 Evaluation methods of design

式中:SCV为生成设计与目标设计的综合相似性指标,该指标越大,代表相似性越高;RSWratio为剪力墙占比率,即图像中剪力墙像素面积(Aswall)占总墙体面积(Aswall+Ainwall)的比例;ηSWratio为剪力墙占比率一致性指标,该指标越大,表明生成设计剪力墙率(RSWratiogenerate)与目标设计剪力墙率(RSWratiotarget)越一致;SSIoU为剪力墙轮廓交并比一致性指标,即生成与目标设计的剪力墙轮廓交集面积(Ainter)与并集面积(Aunion)的比例;SWIoU为剪力墙像素一致性指标,k是总的像素类别(类别0 是背景,类别1 是剪力墙,类别2 是填充墙,类别3 是窗户,类别4是户外门洞);pii是生成正确的像素点数量,pij和pji则是生成错误的像素点数量。

需指出,平面设计结果的评价采用图4(a)所示的基于计算机视觉的评价方法,快速评价生成设计结果与工程师设计结果的一致性。整体结构设计评价采用图4(b)所示的基于结构力学性能的评价方法。整体结构模型构建就以工程师设计模型的参数为基准,剪力墙布置由像素到矢量的自动化程序转化,剪力墙厚度和材料参数以及沿结构高度的变化与工程师设计一致,梁构件设计则根据墙体布置进行适应性调整。基于上述设计和评价方法,保证3 种方法对比条件尽可能一致,更加有效体现差异。

3.2 案例1-7d83m 对比分析

1)案例1-7d83m 基本信息

结构总高度83 m,28 层,平面尺寸为:49.4 m×18.7 m。抗震设防烈度7度(0.1g),特征周期0.55 s。工程师完成的建筑和结构平面设计CAD 及其对应的语义化图纸如图5 所示。

图5 案例1-7d83m 建筑图和结构图Fig.5 Architectural and structural design drawings of Case 1-7d83m

2)案例1-7d83m 平面设计结果对比

在该案例的平面设计中,StructGAN、-TXT、-PHY 三种方法的设计结果以及对应的工程师设计结果如图6 所示。同时,基于SCV指标评价的智能设计与工程师设计的相似性结果也如图6 所示。在该案例中,StructGAN 和StructGAN-PHY 两种方法的设计结果较好,与工程师设计相似性高(即量化SCV指标较高,均为0.47),而StructGAN-TXT 的设计则差别较大(量化SCV指标仅为0.39)。具体分析,StructGAN 和StructGAN-PHY 生成设计的ηSWratio指标接近1,意味着在该案例中二者生成的剪力墙总像素数量与工程师设计的剪力墙像素量基本一致,而StructGAN-TXT 生成设计的ηSWratio指标则偏低,代表剪力墙总量与工程差异较大;3 个案例的SWIoU和SSIoU指标则较为接近(均接近0.5),表明剪力墙布置位置与工程师布置位置差异较小。

图6 案例1-7d83m 的智能生成式结构方案设计图Fig.6 Intelligent designed structural drawings of Case 1-7d83m

3)案例1-7d83m 整体结构性能对比

整体结构案例分析的模型如图7 所示,不同设计结果的剪力墙布置不同,但建模方法与工程师设计均一致。

图7 案例1-7d83m 整体结构设计ETABS 模型Fig.7 Overall structural design of ETABS model of Case 1-7d83m

开展结构动力特性分析以及基于振型分解反应谱法的结构力学响应分析,对应的结构动力特性和层间位移角的对比结果如表1 和图8 所示。StructGAN-TXT 的设计结果与工程师设计差异较大,生成的剪力墙数量偏多,结构刚度偏大,导致动力特性差异20%左右,层间变形的最大差异达到了30%;而StructGAN 和StructGAN-PHY 的设 计与工程师差异则较小,最大性能差异约10%。

表1 案例1-7d83m 结构动力特性与最大层间变形对比Table 1 Comparison on structural dynamic characteristics and maximum inter-story drift ratio of Case 1-7d83m

图8 案例1-7d83m 层间位移角性能对比Fig.8 Comparison on inter-story drift ratio of Case 1-7d83m

3.3 案例2-8d96m、案例3-7d77m、案例4-7d41m对比分析

案例1-7d83m 的分析表明,StructGAN-PHY 与StructGAN 的设计效果较好,为进一步研究不同智能设计方法的通用性和泛化性,进一步开展了案例2-8d96m、案例3-7d77m、案例4-7d41m 的对比分析。案例2-8d96m 的抗震设防烈度8 度(0.2g)、结构总高度96 m(30 层);案例3-7d77m 的抗震设防烈度7 度(0.1g)、结构总高度77 m(26 层);案例4-7d41m 的抗震设防烈度7 度(0.1g)、结构总高度41 m(14 层)。

1)案例2-8d96m 平面与整体设计结果对比

在该案例的平面设计中,StructGAN、-TXT、-PHY 三种方法的设计结果以及对应的工程师设计结果如图9 所示。在该案例中,直观视觉判断3 种方法的设计均与工程师设计非常接近,且量化的SCV指标较高(>0.5)。其中,ηSWratio均较接近1,意味着生成设计的总剪力墙数量与工程设计的基本一致;且SWIoU和SSIoU均大于0.5,意味着剪力墙布置的位置一致性较高。

图9 案例2-8d96m 的智能生成式结构方案设计图Fig.9 Intelligent designed structural drawings of Case 2-8d96m

进一步开展整体结构设计的对比分析,结果如表2 所示。3 种设计与工程师设计的结果差异均较小,动力特性差异在3%左右,最大层间变形在15%左右,且整体变形模式非常接近。其中,StructGAN-TXT 的差异最小(动力特性差异1.00%,最大层间变形差异3.79%)。

表2 案例2-8d96m 结构动力特性与最大层间变形对比Table 2 Comparison on structural dynamic characteristics and maximum inter-story drift ratio of Case 2-8d96m

2)案例3-7d77m 平面设计结果对比

该案例的3 种方法设计结果以及对应SCV评价结果如图10 所示。可以看到,StructGAN 和StructGAN-PHY 两种方法的设计结果较好,与工程师设计相似性高,而StructGAN-TXT 的设计则差别较大。主要原因在于StructGAN-TXT 设计的剪力墙墙体过多,ηSWratio指标偏小,仅0.65。

图10 案例3-7d77m 的智能生成式结构方案设计图Fig.10 Intelligent designed structural drawings of Case 3-7d77m

进一步分析整体结构案例对应的结构动力特性和层间位移角的对比结果如表3 所示。与平面设计结果的对比类似,StructGAN-TXT 的设计结果与工程师设计结果差异较大,剪力墙布置偏多,导致性能差异20% 左右,而StructGAN 和StructGAN-PHY的设计与工程师差异则为10%左右。

表3 案例3-7d77m 结构动力特性与最大层间变形对比Table 3 Comparison on structural dynamic characteristics and maximum inter-story drift ratio of Case 3-7d77m

3)案例4-7d41m 平面设计结果对比

该案例的3 种方法设计结果以及对应SCV评价结果如图11 所示。可以看出,StructGAN-TXT 和StructGAN-PHY 两种方法的设计结果较好,与工程师设计相似性高,而StructGAN 的设计则差别较大。主要原因在于StructGAN 设计的剪力墙墙体偏少,ηSWratio、SWIoU、SSIoU指标均偏小。

图11 案例4-7d41m 的智能生成式结构方案设计图Fig.11 Intelligent designed structural drawings of Case 4-7d41m

整体结构案例分析对应的结构动力特性和层间位移角的对比结果如表4 所示。与平面设计结果的对比类似,StructGAN 的设计结果与工程师差异较大,其剪力墙布置偏少,动力特性差异20%左右,层间变形的最大差异达到了40%;而StructGANTXT 和StructGAN-PHY 的设计与工程师差异同样偏大,动力特性差异10%左右,层间变形最大差异达到了20%左右。主要原因是7 度41 m 设计条件下对应的剪力墙需求较少,剪力墙布置的较小差异都容易引起较大的结构整体特性不同。

表4 案例4-7d41m 结构动力特性与最大层间变形对比Table 4 Comparison on structural dynamic characteristics and maximum inter-story drift ratio of Case 4-7d41m

3.4 基于计算机视觉与力学分析的评价分析

通过对4 个典型案例的平面设计结果和整体结构分析和对比,分析结果汇总至表5 中。可以看出:1)在不同案例中,StructGAN-PHY 是3 种方法中最有效的,较少受到数据质量的限制,能更准确地匹配对应设计需求,且设计结果与工程师设计最接近;2)StructGAN 和StructGAN-TXT 的设计能力均受到训练数据的制约,仅有部分案例效果较好,意味着如果待设计的建筑与训练数据的特征域较为接近,则对应的设计结果较良好,反之则设计质量不佳;3)结构抗震性能需求较高的案例,设计结果通常较好,原因在于所需布置的剪力墙较多,设计变化较少;对于抗震性能需求较少的案例,对应布置的剪力墙较少,设计变化则较丰富,人工智能不能保证找到最合适的结果。

表5 不同案例设计结果分析对比Table 5 Comparison on different case studies

既有研究认为当生成设计与目标设计的交并比大于0.5 则相对较优秀[12-14],但目前尚未明确基于计算机视觉的平面设计评价指标与结构力学性能评价指标之间的关系。将所有案例分析结果绘于图12 中,可以看出:1)当SCV大于0.5 时,基本可以保证智能设计的结构力学性能与工程师设计的结构动力特性差异在5%以内,且层间位移角在10%以内;2)对于地震作用比较大的情况,例如8 度(0.2g)设防、96 m 案例中,即使SCV相差仅0.02 且动力特性差异在5%以内,层间位移角的差异仍旧会超过10%,主要原因在于较大地震作用会导致更大的结构变形,相对差异也变得更大。

图12 基于计算机视觉与结构力学分析的评价指标相关性Fig.12 Indexes relationship between computer visionbased and structural analysis-based evaluation

因此,可以认为基于计算机视觉的SCV指标与结构力学性能评价的指标具有高度的正相关性,即SCV越高则结构力学性能越优,更接近专家的优化设计。SCV=0.5 可作为基于视觉评价的合理性阈值,其对应的结构动力特性差异约为5%,层间位移角差异约为10%。

4 结论

针对建筑结构智能化方案设计方法进行了对比分析研究。从数据驱动到物理增强数据驱动方法,对比了核心算法构建和设计性能,并通过典型案例对比了相同设计条件下不同算法的实际表现,主要结论如下:

1)数据驱动与物理增强数据驱动是建筑结构智能化设计的两个阶段,数据驱动方法搭建了智能化设计的数据、算法、评价和应用方法的基础,物理增强数据驱动方法则进一步提升了算法性能。

2)数据驱动方法(StructGAN 和StructGANTXT)通常受限于特定设计条件下的数据质量与数量,而物理增强数据驱动(StructGAN-PHY)则可以有效降低对结构设计数据的依赖性。

3)对于抗震设防需求较低的设计,剪力墙布置的需求较少,其布置的位置更加灵活,导致智能化设计与工程师设计存在一定差异,未来将进一步提升相关设计条件下的智能化设计能力。

4)确定了基于计算机视觉与力学分析评价方法的正相关性,且SCV=0.5 可作为基于视觉评价的合理性阈值,可供未来智能化设计评价方法使用。

5)主要针对剪力墙结构开展了智能设计相关研究,框架、框架-剪力墙等多种结构类型的智能设计方法研究有待进一步从结构表达、智能算法和评估方法等方面开展。

采用的4 个典型案例分析已在GitHub 中开源(https://github.com/wenjie-liao/StructGAN-PHY/blob/main/StructGAN-TXT-PHY.zip)。

猜你喜欢
剪力墙工程师结构设计
高层建筑连体结构设计
《机械工程师》征订启事
Kenoteq的工程师研发环保砖块
结构设计优化在房屋建筑结构设计中的应用
民用建筑结构设计中短肢剪力墙技术应用
剪力墙结构设计应用浅析
青年工程师
非加劲钢板剪力墙显著屈服位移角的研究
BIM结构设计应用
某排水虹吸井结构设计