新论文 | 基于物理驱动降维概率密度演化方程的高层钢筋混凝土结构精细化随机响应与抗震可靠度分析

文摘   科学   2024-07-14 14:54   德国  

Refined probabilistic response and seismic reliability evaluation of high-rise reinforced concrete structures via physically driven dimension-reduced probability density evolution equation

基于物理驱动降维概率密度演化方程的高层钢筋混凝土结构精细化随机响应与抗震可靠性分析

摘要
大型钢筋混凝土结构的动力可靠性分析是工程实践中最具挑战性的问题之一。尽管在过去的几十年中,人们对混凝土结构的力学分析进行了大量的研究,但结构参数和激励的随机性对其动力性能 (复杂非线性、损伤、能量耗散和塑性) 的显著影响仍无法忽视,这使得评估高自由度实际复杂结构的非线性随机响应和动力可靠性存在很大困难。本研究提出了一种物理驱动方法,通过综合精细化力学分析基于物理的不确定性传,对实际钢筋混凝土结构进行精细的概率响应和地震可靠性评估。在该方法中,材料参数可视为通过蔓式连接函数 (vine copulas) 表征的概率相依随机变量,地震动通过非平稳 Clough-Penzien 谱进行建模。任意感兴趣响应量的不确定性传播由降维概率密度演化方 (dimension-reduced probability density evolution equation, DR-PDEE) 控制。降维概率密度演化方程中的本征漂移函数是不确定性传播的物理驱动力,可以通过结构的代表性动力分析数据来识别。通过求解物理驱动的降维概率密度演化方程,可以计算结构的时变可靠度,这无法通过由于计算成本过高的一般蒙特卡罗模拟实现。最后,本研究展示了一个实际工程应用,对具有近 28 万个自由度的 24 层钢筋混凝土剪力墙结构进行了随机响应和抗震可靠性分析。

研究背景

近年来提出了降维概率密度演化方程,为解决高维非线性随机动力系统的响应与可靠性分析提供了一类高效精细化的实现途径。不论系统维数多高,若仅关心系统的某一维感兴趣响应过程,且其仅需满足路径连续性条件,即可证明,这一过程的瞬时概率密度函数精确满足一个一维偏微分方程,即降维概率密度演化方 (dimension-reduced probability density evolution equation, DR-PDEE) (Lyu & Chen 2022)。方程中的本征漂移函数是待定的,对于某些特殊高维系统可以给出其解析表达 (Sun & Chen 2022, Luo et al. 2023),而对于更一般的系统,需要根据代表性动力分析数据进行数值构造 (Lyu & Chen 2021, Chen & Lyu 2022, Hao et al. 2024)。以往的研究工作展示了这一方法对不同随机响应与动力可靠性问题的广泛适用性与有效性 (Luo et al. 2022a, b, Sun et al. 2023, Chen & Lyu 2024, Chen et al. 2024, 律梦泽 & 陈建兵 2024)。在此基础上,本研究进一步将这一方法应用于复杂实际工程结构问题,结合大型工程结构的精细化确定性非线性有限元动力分析方法,实现其随机动力响应与可靠性分析。

方法概述

混凝土结构精细化抗震非线性分析的控制方程包括平衡方程变形协调方程本构方程:
其中本构方程中的塑性应变和损伤因子可由混凝土弹塑性损伤本构关系确定 (李杰 等 2014)。控制方程中考虑了混凝土材性的概率相依随机参数,包括弹性模量、抗拉/压强度、拉/压峰值应变、拉/压下降段参数:

各不确定性参数之间的概率相依性由蔓式连接函数进行表征 (Tao et al. 2020);结构遭受的随机地震动采用非平稳调制的 Clough-Penzien 谱进行建模,并采用线性滤波法模拟。
在结构参数与外部激励的不确定性下,结构响应需视为随机过程。认为结构某一关键响应过程在首次超越准则下超过某一既定安全阈值,则结构发生破坏,则可定义结构的首次超越动力可靠度:
结构的感兴趣响应过程的瞬时概率密度函数满足降维概率密度演化方程:
从结构参数与激励的不确定性输入到结构响应随机性输出的不确定性传播过程如下图所示:
联合结构的确定性物理方程与随机性传播的降维概率密度演化方程,则可建立复杂结构系统中物理驱动概率密度演化的综合控制方程:
对上述方程组采用适当数值方法求解,即可获得复杂结构系统中任意感兴趣响应过程的瞬时概率密度函数及其在既定安全阈值下发生首次超越破坏的时变失效概率。具体数值步骤包括:
(1) 随机结构参数与外部激励的代表性数据生成: 采用基于蔓式连接函数模型的条件采样方法给出结构随机参数的代表性数据;同时采用线性滤波法生成非平稳随机地震动的代表性时程。
(2) 结构的代表性确定性动力分析: 将随机结构参数与外部激励的代表性数据作为结构有限元数值模型的输入参数,对结构进行代表性有限元动力分析,获得结构感兴趣响应过程的代表性时程数据。
(3) 本征漂移函数数值构造: 根据结构的代表性响应时程分析数据,采用局部加权回归方法数值估计结构响应过程 (针对随机响应分析) 或其吸收边界过程 (针对动力可靠性分析) 的降维概率密度演化方程中的本征漂移函数数值解。

(4) 降维概率密度演化方程求解: 将本征漂移函数的数值解代入降维概率密度演化方程,采用路径积分法数值求解,获得结构响应的瞬时概率密度函数或时变失效概率的数值解。

数值结果

工程案例: 24 层钢筋混凝土剪力墙结构

结构 ABAQUS 有限元模型 (共 277404 自由度)

结构材性

有限元建模中的纤维梁模型与分层壳模型

结构前六阶模态与特征周期

混凝土材性参数的边缘概率分布

混凝土参数间概率相依性表征的蔓式连接函数

结构随机参数的概率分布与相依性蔓式连接函数模型

随机材性参数下混凝土单轴应力—应变曲线

Clough-Penzien 谱与非平稳地震动代表性时程

代表性地震动激励下构的非线性损伤特征

不同强度地震动激励下结构的非线性动力特征

结构各层层间位移角的代表性分析数据

(采用天河 2 号超级计算机进行400 次结构分析获得)

结构第 5 层层间位移响应的瞬时概率密度函数数值解

结构第 5 层层间位移响应的标准差与峰度时程数值解

不同安全阈值下结构的时变失效概率

不同安全阈值下仅考虑激励随机性与同时考虑结构参数和激励随机性的结构时变失效概率对比

研究结论
本研究提出了一种通过综合精细化的力学分析基于物理的不确定性传播,来分析实际钢筋混凝土结构在随机地震作用下的概率响应和抗震可靠性的物理驱动方法。在所提方法中,混凝土参数和地震激励的不确定性分别通过概率相依性量化和线性滤波技术建模。然后,真实物理随机系统中的不确定性传播可以由降维概率密度演化方程控制,该方程是关于感兴趣响应量瞬时概率密度函数的二维偏微分方程。降维概率密度演化方程中的本征漂移函数是驱动随机传播的物理 “驱动力”,可以通过结构精细有限元模型的代表性动力分析数据来识别。求解降维概率密度演化方可以获得不同阈值下的瞬时概率密度函数和动力可靠性。本文通过一个近 28 万自由度的高层钢筋混凝土剪力墙结构实际工程案例,展示了所提方法的工程适用性。结论包括:
  1. 基于概率相依性量化和物理驱动降维概率密度演化方程,提出了实际工程结构的 “不确定性表征—传播—可靠性分析完整框架。可直接应用于工程实践。
  2. 求解物理驱动的降维概率密度演化方程,可以获得工程结构在不同阈值下关键响应的时变失效概率。这可以用于指导不同地震水平下的结构安全性评估。
  3. 考虑结构参数不确定性对于工程实践中罕遇事件下小失效概率基于可靠性的结构设计十分重要。

在未来工作中,可以基于物理驱动的降维概率密度演化方程在建筑、桥梁和岩土工程中进行更多的工程应用。在所提框架中,还可以进一步考虑结构参数和激励的认知不确定性和空间变异性。

论文引用格式 (GB):
Lyu MZ, Chen JB, Shen JX, 2024. Refined probabilistic response and seismic reliability evaluation of high-rise reinforced concrete structures via physically driven dimension-reduced probability density evolution equation [J]. Acta Mechanica, 235: 1535-1561.

DOI:
10.1007/s00707-023-03666-4

论文作者介绍:
律梦泽 Meng-Ze Lyu
同济大学土木工程学院博士后, 合作导师是陈建兵教授
邮箱:lyumz@tongji.edu.cn
主页:
https://www.researchgate.net/profile/Meng-Ze-Lyu

陈建兵 Jian-Bing Chen, 通讯作者
同济大学土木工程学院特聘教授, 国家杰出青年科学基金获得者

邮箱:chenjb@tongji.edu.cn

主页:

https://jcersm.tongji.edu.cn/1e/f9/c13092a139001/page.htm

申家旭 Jia-Xu Shen
北京工业大学建筑工程学院助理研究员

邮箱:shenjx@bjut.edu.cn
主页:

https://yanzhao.bjut.edu.cn/info/1488/15674.htm

相关参考文献:

[1] Chen JB, Lyu MZ, 2022. Globally-evolving-based generalized density evolution equation for nonlinear systems involving randomness from both system parameters and excitations [J]. Proceedings of the Royal Society A - Mathematical Physical & Engineering Sciences, 478 (2264): 20220356.

[2] Chen JB, Lyu JH, Spanos PD, Li J, 2024. Stochastic response analysis of a spar-type fowt subjected to extreme waves by a novel filter wave model and the DR-PDEE [J]. Journal of Engineering Mechanics, 150 (4): 04024009.

[3] Chen JB, Lyu MZ, 2024. Probabilistic response determination of high-dimensional nonlinear dynamical systems enforced by parametric multiple Poisson white noises [J]. Nonlinear Dynamics, 112: 11283-11298.

[4Hao TT, Yan WJ, Chen JB, Sun TT, Yuen KV, 2024. Multi-output multi-physics-informed neural network for learning dimension-reduced probability density evolution equation with unknown spatio-temporal-dependent coefficients [J]. Mechanical Systems & Signal Processing, 220: 111683.

[5] Luo Y, Chen JB, Spanos PD, 2022a. Determination of monopile offshore structure response to stochastic wave loads via analog filter approximation and GV-GDEE procedure [J]. Probabilistic Engineering Mechanics, 67: 103197.

[6] Luo Y, Spanos PD, Chen JB, 2022b. Stochastic response determination of multi-dimensional nonlinear systems endowed with fractional derivative elements by the GE-GDEE [J]. International Journal of Non-Linear Mechanics, 147: 104247.

[7] Luo Y, Lyu MZ, Chen JB, Spanos PD, 2023. Equation governing the probability density evolution of multi-dimensional linear fractional differential systems subject to Gaussian white noise [J]. Theoretical & Applied Mechanics Letters, 13: 100436.

[8] Lyu MZ, Chen JB, 2021. First-passage reliability of high-dimensional nonlinear systems under additive excitation by the ensemble-evolving-based generalized density evolution equation [J]. Probabilistic Engineering Mechanics, 63: 103119.

[9] Lyu MZ, Chen JB, 2022. A unified formalism of the GE-GDEE for generic continuous responses and first-passage reliability analysis of multi-dimensional nonlinear systems subjected to non-white-noise excitations [J]. Structural Safety, 98: 102233.

[10] Sun TT, Chen JB, 2022. Physically driven exact dimension-reduction of a class of nonlinear multi-dimensional systems subjected to additive white noise [J]. ASCE-ASME Journal of Risk & Uncertainty in Engineering Systems, Part A - Civil Engineering, 8 (2): 04022012.

[11] Sun TT, Lyu MZ, Chen JB, 2023. Property of intrinsic drift coefficients in globally-evolving-based generalized density evolution equation for the first-passage reliability assessment [J]. Acta Mechanica Sinica, 39: 722471.

[12] Tao JJ, Chen JB, Ren XD, 2020. Copula-based quantification of probabilistic dependence configurations of material parameters in damage constitutive modeling of concrete [J]. Journal of Structural Engineering, 146 (9): 04020194.

[13] 李杰, 吴建营, 陈建兵, 2014. 混凝土随机损伤力学 [M]. 北京: 科学出版社.

[14] 律梦泽, 陈建兵, 2024. 非平稳随机激励下高维非线性系统可靠度分析的概率密度全局演化方法 [J]. 振动工程学报 (在线发表).

工程可靠性与随机力学
同济大学工程可靠性与随机力学国际联合研究中心 (JCERSM) 成立于2016年。中心中方主任为中国科学院院士李杰教授,外方主任为美国工程院院士、中国科学院外籍院士 Spanos 教授。
 最新文章