基于拉丁超立方抽样的蒙特卡罗数值模拟在桥梁结构可靠度分析中的应用基于拉丁超立方抽样的蒙特卡罗数值模拟在桥梁结构可靠度分析中的应用摘要：根据《公路工程结构可靠度设计统一标准》概率原理已经被引入并作为结构设计的一般方法。运用蒙特卡罗数值模拟对结构进行可靠度分析也已经被广泛接受，但是，蒙特卡罗法的抽样方法还存在优化的空间，拉丁超立方抽样法作为抽样研究领域的最新研究成果，运用到蒙特卡罗法中，可以提高抽样的效率、节省样本空间并减少运算时间。关键字：拉丁超立方抽样；桥梁可靠度；蒙特卡罗；标准误；线性功能函数中图分类号：[TU997]文献标识码：A引言梁结构的安全性是桥梁设计中的重大问题。桥梁工程的建设耗资巨大，一旦失效不仅会造成结构本身和人民生命财产的巨大损失，还往往产生难以估量的次生灾害和附加损失。桥梁结构在设计、施工、使用过程中具有种种影响结构安全、使用、耐久的不确定性因素。这些不确定性使我们在判断结构是否可靠时，不能用一个简单的“是”或“否”来衡量，而必须采用以概率形式为基础的可靠度指标来表达。在这样一个背景下，结构的可靠度理论与方法便应运而生。[1]目前，蒙特卡罗法在结构可靠度设计中已经成为一种普遍方法，得到了广泛的应用。作为一种基于概率论的数值模拟方法，抽样方法的好坏决定了其解的精确性，而现在常用的抽样方法在精确性、稳定性和可操作性方面均有改进的空间。[2]拉丁超立方抽样作为抽样领域的最新研究成果，运用到蒙特卡罗法中，可以提高抽样的效率，节省样本空间并减少运算时间。本文将拉丁超立方抽样应用到蒙特卡罗法中，并对北盘江大桥进行可靠度分析，通过实际抽样结果的对比证明了拉丁超立方抽样法较蒙特卡罗直接抽样的明显改进，得到了令人满意的结果。拉丁超立方抽样在北盘江大桥可靠度分析中的应用北盘江大桥工程介绍北盘江特大桥为预应力混凝土悬浇空腹（斜腿）式连续刚构桥（如图1所示）。图1主桥桥型布置图（单位：cm）主桥部分总长度度为895m，桥跨布置为82.5+220+290+220+82.5m。上部结构采用连续变截面现浇混凝土箱梁，宽度为10.5m，由钢筋混凝土桥墩支撑，主墩的高度分别为75m、90m、143m、68m，均采用双肢薄壁空心截面，其中中间两对桥墩设有空腹斜腿以提高根部刚度。桥台部分的支撑类型为简支。桥墩的纵向配筋主要采用直径32mm，28mm，20mm的HRB335钢筋，横向箍筋为直径为16mm和20mm的HRB335钢筋。有限元建模我们运用MIDAS2010对北盘江大桥进行有限元建模。有限元模型如图2所示。图2有限元模型三维视图为了便于后续分析，我们将梁截面从左至右编为1-286号，具体编号见图3所示。求出各截面的功能函数考虑到桥梁是变截面桥梁，且预应力钢筋形状复杂，直接得出各截面的荷载与荷载效应之间的关系有困难，所以我们采用的方法是运用MIDAS2010对各种荷载单独作用下的荷载效应值进行分析，再不断变换该荷载的荷载值，得出不同荷载值对应的荷载效应值，运用MATLAB的统计工具箱进行线性拟合，得出不同荷载与相应荷载效应之间的线性关系，进而可求出各截面的功能函数。求出桥梁各截面功能函数变量的概率分布根据《公路工程结构可靠度设计统一标准》的有关规定，并综合考虑北盘江大桥的实际情况及相关文献[9]，现取各荷载与荷载效应的概率分布如表1所示：取截面抗力概率分布如表2所示。表1北盘江大桥荷载及荷载效应概率分布表2北盘江大桥截面抗力概率分布运用拉丁超立方抽样求出给截面的失效概率本文通过MATLAB进行编程，进而求得各截面的失效概率。为了体现拉丁超立方抽样在结构功能函数为线性时相对于蒙特卡罗直接抽样的改进，本文用两种不同的方法进行抽样，并对结果进行对比。拉丁超立方抽样和蒙特卡罗直接抽样的比较抽样误差率的比较选定比较量为，考察两种抽样方法的样本均值相对于数学期望的误差。考虑到样本数量对样本均值的影响，本文对不同样本数量下的误差进行考察。表3给出了样本容量为N=100、1000、10000时的各截面抗弯功能函数误差率最大、最小值的比较。表4给出了样本容量为N=100、1000、10000时的各截面抗剪功能函数误差率最大、最小值的比较。表3各截面抗弯功能函数误差率...