(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202210429658.9 (22)申请日 2022.04.22 (71)申请人 交通运输部公路科 学研究所 地址 100088 北京市海淀区西土城路8号 (72)发明人 白冰　王磊　 (74)专利代理 机构 北京市盛峰律师事务所 11337 专利代理师 于国强 (51)Int.Cl. G06F 30/13(2020.01) G06F 30/20(2020.01) G06F 111/08(2020.01) G06F 119/02(2020.01) (54)发明名称 一种桥梁结构多设计点可靠度分析方法 (57)摘要 本发明涉及工程结构可靠性分析领域， 尤其 涉及一种桥梁 结构多设计点可靠度分析方法， 所 述方法包括根据桥梁结构分析确定其原始结构 功能函数和设计点位置和数量； 选取一设计点， 并利用附加函数屏蔽其余设计点， 确定屏蔽后结 构功能函数； 基于改进的交叉熵方法， 对屏蔽后 结构功能函数进行MCMC预抽样并确定重要抽样 函数分布 参数； 通过对所述屏蔽后结构功能函数 进行重要抽样， 计算得到当前选取的设计点对应 的结构失效概率； 通过叠加所有设计点的结构失 效概率后得到结构整体的失效概率。 本发明方案 实现了针对结构多设计点可靠度问题的高效稳 健求解， 有利于对大型桥梁结构进行可靠性优 化、 结构设计和结构评估。 权利要求书3页 说明书12页 附图6页 CN 114692282 A 2022.07.01 CN 114692282 A 1.一种桥梁结构多设计点可靠度分析 方法， 其特 征在于， 所述方法包括以下步骤： 基于桥梁结构的结构参数、 功能参数和所承受荷载确定原 始结构功能函数； 根据所述原始结构功能函数确定设计点位置和设计点数量， 并针对不同设计点构造对 应的附加函数； 选取一设计点， 并利用所述附加函数屏蔽其 余设计点， 确定屏蔽后结构功能函数； 基于改进的交叉熵方法， 对所述屏蔽后结构功能函数进行MCMC预抽样并根据所述MCMC 预抽样的样本确定 重要抽样概 率密度函数的参数值； 采用所述重要抽 样概率密度函数对所述屏蔽后结构功能函数进行重要抽 样， 以确定当 前选取的设计点对应的结构失效概 率； 待获取所有设计点对应的结构失效概 率后， 通过叠加组合得到结构整体的失效概 率。 2.根据权利要求1所述的桥梁结构多设计点可靠度分析方法， 其特征在于， 所述附加函 数为： 其中， ri为屏蔽曲面半径， si为比例系数， xi*为所要屏蔽的设计点； 所述ri及si的数值分 别由下式确定： ri＝γβform,i 其中， βform,i为xi*对应的一阶可靠度指标， 为原始结构功能函数G(xi*)的梯度， γ和 δ 为经验参数。 3.根据权利要求1所述的桥梁结构多设计点可靠度分析方法， 其特征在于， 所述屏蔽后 结构功能函数为： Gj(x)＝G(x)+ΣBi(x) 其中， Gj(x)为当前选取设计点的屏蔽后结构功能函数， G(x)为原始结构功能函数， Bi (x)为被屏蔽设计点xi*所对应的附加函数。 4.根据权利要求2所述的桥梁结构多设计点可靠度分析方法， 其特征在于， 当所述附加 函数Bi(x)中||x‑xi*||＝ri时， Bi(x)的一阶导数及二阶导数分别为： Bi′(x|||x‑xi*||＝ri)＝0 Bi″(x|||x‑xi*||＝ri)＝8siri2>0。 5.根据权利要求1所述的桥梁结构多设计点可靠度分析方法， 其特征在于， 所述对所述 屏蔽后结构功能函数进行M CMC预抽样， 具体包括以下步骤： 步骤1： 随机生成一条马氏链的n维初始变量x1＝[x1(1),…,x1(i),…,x1(n)]∈F； 步骤2： 由第k轮抽样获得的第i维变量xk(i)构造正态提议分布φi(·； xk(i))， 自所述正 态提议分布中抽取样本xi′； 步骤3： 计算 转移核：权　利　要　求　书 1/3 页 2 CN 114692282 A 2步骤4： 由均匀分布[0,1]随机抽取样本Ui， 并令 步骤5： 重复步骤2 ～4， 直至 完成n维 抽样， 确定第k+1轮生成样本xk+1为： 步骤6： 重复所述 步骤2～5， 直至 完成所要求数量的样本抽取。 6.根据权利要求5所述的桥梁结构多设计点可靠度分析方法， 其特征在于， 所述重要抽 样概率密度函数的参数值包括均值和协方差阵； 所述均值基于所述M CMC预抽样的样本确定， 所述均值计算 为： 其中， 为所述均值， xk为所述MCMC预抽样的样本 。 7.根据权利要求6所述的桥梁结构多设计点可靠度分析方法， 其特征在于， 所述协方差 阵基于所述M CMC预抽样的样本和所述均值确定， 具体为： 以所述均值所指方向为坐标轴构建新坐标系， 并基于所述MCMC预抽样的样本计算所述 均值所指方向的方差， 利用转换矩阵获取原坐标系下的协方差阵； 所述协方差阵计算 为： 其中， 为所述均值所指方向的方差， 为所述协方差阵。 8.根据权利要求1所述的桥梁结构多设计点可靠度分析方法， 其特征在于， 所述确定当 前选取的设计点对应的结构失效概 率， 具体为：权　利　要　求　书 2/3 页 3 CN 114692282 A 3