基于计算可视化的土木工程可靠度课程教学方法_张洁.pdf

高等建筑教育 2023 年第 32 卷第 2 期JOURNAL OF ARCHITECTURAL EDUCATION IN INSTITUTIONS OF HIGHER LEARNING Vol.32 No.2 2023 doi:10.11835/j.issn.1005-2909.2023.02.017欢迎按以下格式引用:张洁,黄宏伟.基于计算可视化的土木工程可靠度课程教学方法J.高等建筑教育,2023,32(2):132-141.修回日期:2021-06-23基金项目:上海市教育发展基金会和上海市教育委员会“曙光计划”资助项目(19SG19)作者简介:张洁(1980),男,同济大学土木工程学院教授,博士,主要从事工程风险相关的教学和科研研究,(E-mail)。基于计算可视化的土木工程可靠度课程教学方法张 洁,黄宏伟(同济大学 土木工程学院,上海 200092)摘要:随着基于可靠度理论的结构设计方法的快速发展,可靠度理论已经成为土木工程专业学生需要掌握的重要知识。由于可靠度理论较为抽象、常涉及复杂的数值计算,学生常无法完成学习流程,导致学生学习兴趣匮乏,学习效果有待提高。针对上述难点,提出了一种基于计算可视化的教学方法。该教学方法基于 EXCEL 这一通用工具,可将各种可靠度算法的中间过程和结果进行直观分析和展示,消除可靠度理论学习中涉及的数值分析要求,由此精简学习流程。该方法可帮助学生从复杂的数值计算中解脱出来,将精力聚焦在可靠度分析原理上,方便和快速地理解和掌握相关知识点,并将所学知识灵活应用于自己的研究领域。基于计算可视化的教学方法可大幅提高可靠度课程的学习效率和学习获得感,为土木工程可靠度理论的教学提供一个行之有效的方法。关键词:计算可视化;土木工程;可靠度理论;教学方法中图分类号:TU 201.1;G642 文献标志码:A 文章编号:1005-2909(2023)02-0132-10土木工程的设计、施工和运营中存在大量的不确定性因素,如地震的不确定性、风荷载的不确定性、材料的不确定性、计算模型的不确定性等1。这些因素的存在导致土木工程设施的性能很难准确预测。可靠度理论采用概率方法定量模拟土木工程中的各类不确定性因素,采用失效概率刻画不确定环境中工程结构的服役性能,可定量衡量各类不确定性因素对土木工程设施服役性能的影响,是现代结构规范的理论基础2-4。自 1992 年我国颁布工程结构可靠度设计统一标准第一版(GB 5015392)以来5,土木工程的相关规范已经逐步从以经验为主的安全系数设计方法过渡到基于可靠度理论的极限状态设计方法6。目前我国实行的最新可靠度规范为 2008 年颁布的工程结构可靠性设计统一标准修订版(GB 501532008)7。该规范又被称为“规范的规范”或“规范之母”。在此背景下,可靠度理论已被多所国内外高校列为土木工程硕士研究生的专业基础课张 洁,等 基于计算可视化的土木工程可靠度课程教学方法程8。以同济大学为例,于 2009 年率先开设了研究生课程土木工程的概率分析,主要阐述土木工程结构可靠度的基本概念与原理、失效概率计算方法、基于可靠度的设计理论等内容,旨在使学生在深入理解可靠度原理的基础上,培养利用可靠度理论分析土木工程结构的安全性和适用性等问题的能力。应教学国际化需求,该课程于 2010 年转为全英文面向国内、国际学生授课,目前已开课 10 余年。学习土木工程可靠度分析课程除需要土木工程相关专业课知识外,还需要较深的概率与统计、优化原理、数值计算方法、编程等知识和技巧。由于该课程教学内容涉及的数学公式多、专业背景广、理论抽象,教师与学生往往遇到教学不易、学习困难的两难困境,不但影响学生的学习热情和兴趣,还容易使学生产生畏惧感,难以达到预期的教学目的。如何提高土木工程可靠度课程的教学效果一致是学界关心的热点问题。例如,杨世浩和齐甦从教学体系和教学实践两个方面讨论了改进措施,认为可靠度教学不仅要重视教学内容方法,也要提高教师素质9;潘文军等认为可靠度教学的难点在于内容抽象、枯燥,建议加强教学与实践应用的结合,提高学生兴趣10;李会军和李宗利建议加强学生的软件使用能力及编程能力,强化典型实际工程案例的研究,从而改善教学质量8。这些研究很好地分析了目前可靠度教学的挑战以及今后的主要改革方向,为土木工程可靠度课程建设积累了宝贵的经验。为提高土木工程可靠度课程的教学效果,同济大学土木工程可靠度课程教学团队经过长期教学实践,提出了一种基于计算可视化的土木工程可靠度教学方法。多年实践表明,该方法可有效缩短学习流程,大幅降低教师教学和学生学习的难度,在帮助学生更深刻理解可靠度原理的同时,快速将所学的理论知识应用到相关研究中。本文的目的是对这一教学方法进行介绍,为促进土木工程可靠度课程的高质量建设提供新的途径。一、土木工程可靠度课程教学中面临的关键挑战土木工程可靠度分析课程主要存在以下三点挑战。(1)教学内容抽象。可靠度课程主要讲授土木工程中的各类不确定性现象及其模拟方法、不确定性条件下工程结构失效概率计算方法、基于可靠性的设计方法等内容。上述内容涉及的数学知识多,很难通过实体进行展示,内容较为抽象,教师难于讲授、学生难于理解,增加了教学难度11。在此情况下,教学内容的讲授只能浅尝辄止,很少有时间安排实践课,导致学生们对可靠度理论掌握不牢,影响学生的学习积极性。(2)编程要求高。土木工程可靠度问题很少有解析解,各类可靠度分析方法常需采用不同数值计算算法,难以直接通过手算完成。这导致教师授课中无法直观展示计算过程,导致实例分析如同黑箱操作,示意性弱12。对于学生而言,完成相应的习题也需要借助 Matlab 等平台进行编程,需要较高的编程技巧。即使对于少部分具有较好编程基础或对编程具有浓厚兴趣的同学,采用编程进行学习也会耗费相当多的学习时间。对于缺乏编程基础或编程能力较弱的同学,只能对各类可靠度分析方法“望洋兴叹”。较高的编程门槛大大增加了学生的学习难度,增加了畏难情绪。(3)学习获得感不高。由于内容抽象、编程要求高,采用传统教学方法将不得不花大量时间在讲授各类算法及其实现过程上13,导致很少有时间讲解可靠度理论在工程中的工程实例应用,以及理解和解读可靠度分析的结果。上述应用和结果解读的缺乏导致学生在学习时感觉课程偏于数学331 高等建筑教育 2023 年第 32 卷第 2 期 和理论,与土木工程的联系较弱。由于对可靠度理论在土木工程应用方面缺乏认知,因此学习获得感低,这也违背了土木工程可靠度课程采用可靠度理论解决土木工程问题的初衷。二、基于计算可视化教学方法的核心思想以上分析表明,可靠度理论课程教学的关键挑战在于内容抽象,计算复杂,难以掌握。传统的教学方法,其学习流程共分为四个步骤:(1)学习可靠度分析原理;(2)掌握可靠度分析的算法;(3)根据算法采用手算或者编程进行计算获得结果;(4)通过敏感性分析等手段对结果进行理解。在传统教学方法中,步骤(3)的实施需要进行大量手算或编程,耗费学生大量的时间和精力。当采用手算时,步骤(4)常因为计算量过大而难以实施,导致学生对可靠度分析结果缺乏解读和理解。步骤(3)和步骤(4)的实现困难,使不少学生无法完成整个学习流程,或在学习过程中尝浅则止,或无法完成足够的练习,导致学习效果欠佳。相比而言,EXCEL 作为一种通用的办公软件,具有敏感性分析、优化分析、随机数生成等功能。针对土木工程可靠度分析中的教学难点,同济大学可靠度课程教学团队提出通过 EXCEL 这一通用工具,将各类可靠度分析算法在 EXCEL 中实现,便捷进行计算和参数敏感性分析,由此将可靠度学习流程精简为三步:(1)学习可靠度分析原理;(2)掌握可靠度分析算法;(3)通过 EXCEL 实现算法并对计算结果进行理解。由于基于 EXCEL 的可靠度分析极易实现,由此消除了传统教学方法中由于计算困难无法完成学习流程的困境。此外,在基于 EXCEL 的可靠度分析中,所有中间过程和结果都可以通过 EXCEL 单元格直观记录和表达,整个计算过程对教师和学生都是“透明”的,这对学生加深对可靠度原理和算法的理解极为有利,因此这一教学方法在本文中成为基于计算可视化的教学方法。不少学生还在基于 EXCEL 的计算可视化教学方法中发现了 EXCEL 强大的计算分析功能,并将其应用到其他课程或科研之中,极大地增加了学生学习的乐趣,将传统可靠度方法学习中的畏难情绪变成了学习兴趣。在土木工程可靠度课程中,中心点法、验算点法、蒙特卡罗抽样以及基于可靠度的分项系数标定是几个典型的知识点和教学难点。下面以上述知识点为例,介绍计算可视化的实现过程。三、计算可视化教学方法典型案例(一)中心点法的计算可视化土木工程可靠度分析中,当功能函数为非线性时,可对该功能函数在随机变量均值点处进行线性泰勒展开,即可利用线性方程的性质获得结构性能指标(如安全系数、位移)均值和标准差的近似值。假定结构性能指标服从某一分布(如正态分布),即可根据性能指标的均值和标准差计算其结构的失效概率。由于该方法中功能函数在随机变量的均值点进行线性展开,故常被称为中心点法。中心点法是可靠度理论中的早期方法,概念简单,但精度较差,一般用于结构的初步分析,通常是结构可靠度理论中首先介绍的可靠度分析方法。当对功能函数进行线性展开时,需要对不同随机变量求偏导数。在传统教学法中,偏导数常需通过手算完成。对于极限状态方程复杂的实例,这一过程极为繁琐。在 EXCEL 中,偏导数可采用有限差分法获得,由此大幅简化中心点法所需的计算工作量。为中心点法的计算可视化过程,考虑如下功能函数:431张 洁,等 基于计算可视化的土木工程可靠度课程教学方法g(x1,x2)=x1+x1x2(1)其中 x1、x2为随机变量。x1的均值为 1、标准差为 2;x2的均值和标准差均为 2。设 g(x1,x2)0表示结构失效。图1给出了基于EXCEL的中心点法实现模板。由有限差分分析可知,功能函数在中心点处的梯度矩阵为3,0.99。根据中心点法计算原理,功能函数的均值()和标准差()可基于下式计算得到14:=1+12(2)=GCxGT(3)其中 1、2分别为 x1、x2的均值,G 为梯度矩阵,Cx为协方差矩阵。留意到中心点法中标准差计算需进行矩阵运算。该矩阵运算可通过 EXCEL 中矩阵运算函数 MMULT 来实现。通过上述计算模板,即可获得该问题的可靠度指标为 0.47。在建立上述计算模板后,通过改变第一部分单元格中对应的数字,即可考察不同情况下可靠度指标的变化,由此分析不同因素对可靠度指标的影响,使用极为方便。图 1 基于 EXCEL 的中心点法分析教师授课中,在通过上述模板对中心点法的概念进行讲解后,学生即可在 510 分钟内自行通过 EXCEL 操作重现计算过程。在课后作业中,学生可方便地将其扩展到浅基础、隧道、边坡等问题的可靠度分析中。如前所述,中心点法通常为可靠度课程中介绍的第一个可靠度分析方法。学生在通过该方法与自己研究方向相关的问题进行分析后,可立刻体会到可靠度理论解决自身研究问531 高等建筑教育 2023 年第 32 卷第 2 期 题的潜力,极大增加对可靠度课程的兴趣和信心。(二)验算点法的计算可视化中心点法无法考虑随机变量分布对可靠度分析结果的影响,且对功能函数的形式具有依赖性。由于验算点法可有效克服上述方法的不足,该方法在结构可靠度分析中受到了广泛重视。以上述可靠度问题为例,其可靠度指标()按验算点法可按下式进行计算15:=ming(x)=0nR-1nT(4)其中 n 为随机变量(x)在标准空间中的转换变量,R 为随机变量的相关系数矩阵。对于正态分布变量,其转换关系为:n=x-xx(5)其中 x、x分别是 x 的均值和标准差。由式(4)可知,验算点法计算可靠度指标是一个有约束优化问题。由于优化问题的求解涉及大量的迭代计算,采用手算法求解有约束优化问题极为繁琐。在传统教学中,部分例题和练习需要学生通过手工迭代的方法求解上述有约束优化问题,耗费了教师和学生大量的时间。如要回避手算的困