基于分数阶模型的重力坝动力可靠性研究.docVIP

基于分数阶模型的重力坝动力可靠性研究 冯宇迪 周星德 ( 河海大学 ，江苏 南京 210098) 关键词： 中图分类号：TU 文献标识码：A Reliability analysis of Gravity dams with A Fractional Order Model Feng Yudi Zhou Xingde Yang Bo Li Yali (College of Civil and Transportation Engineering, Hohai University, Nanjing 210098, PR China) Abstract: Aiming at the flaws of pseudo-static method, a method of calculating dynamic reliability of gravity dam based on a fractional order model was developed. First, relative displacement angle was used as functional function to reflect dynamic time-history reliability. Second, a fractional order model was put forward in consideration of concrete aging phenomena and solved by Outaloup filter and limit state function. Moreover, correlative formula was deduced and Matlab was used to figure up the reliability. At last, the superiority of the proposed method was demonstrated by comparing the results calculated by the proposed method and pseudo-static method． Keyword: gravity dams, dynamic reliability, fractional order 近年来，我国西部地区建设了许多大坝，由于处于高地震区，且坝体存在几何、材料及荷载等不确定性[1]，这就需要对其进行可靠性分析。已经出现了Monte Carlo法[2]、一次二阶矩法[3]、基于copula样本抽样法[4]、基于广义坐标系[5]的可靠度算法等。为了更加接近实际，科研工作者已经开始研究土体正交各相异性产生的空间效应[6]、短期稳定时的动力可靠性[7]、滑动面等可靠性[8]、基于非线性剪切标准的可靠度评估[9]等。 目前有关动力可靠度计算，出于安全考虑，通常取峰值加速度，采用拟静力法计算，对于关键的阻尼项没有涉及，不能准确反映坝体实际的可靠性，为此，本文以重力坝为对象，研究其动力可靠度。首先，为了反映动力时程可靠性，功能函数采用了文献[10]提出的相对位移转角；其次，考虑混凝土存在老化现象[11]，引入了分数阶模型，分数阶求解采用了Outaloup滤波器和状态方程模式；再次，采用Monte-Carlo法计算时程响应，通过K-S法检验了可靠性指标中位移分布模式，进而推导了可靠性计算公式。最后，采用Matlab进行了编程，并与拟静力法计算结果进行了对比分析。 1 含分数阶的振动方程 对于堆石坝，假设系统阻尼为比例阻尼，其分数阶振动模型可以表示为 (1) 式中：、和分别为质量、阻尼和刚度矩阵，均为阶矩阵，，为n阶矩阵为地震波的加速度，为阶。 模态，设前r阶模态为，令 (2) 把式()代入到式()并左乘的转置，得 (3) 式中：，，，，，…分别为系统的阻尼比和圆频率。 式()包含r个单自由度系统，可表示为 (4) 对于初始处于静止状态，式()Laplace变换 ， (5) 对于式(5)中的分数阶可采用Oustalop滤波器近似，假设拟合频率段为，则可构造如下滤波器的传递函数，为 (6) 式中：N2N＋1为滤波器的阶次，，式()，则 (7) 式中：为正整数。把式()转换为状态方程 (8) 式中：和均为维列矢量，为维行矢量，为维矩阵。 把式()组装为状态方程，为 (9) 式中，，，。 图1 重力坝计算模型