本发明涉及土木工程领域的桥梁碰撞间隙设置宽度的评估，特别涉及非平 稳地震作用下桥梁系统碰撞间隙宽度的概率评价方法，具体地说是基于非平稳 随机振动的虚拟激励法、精细积分法和矩阵运算方法的桥梁系统概率的计算方 法，可用于判断既有桥梁碰撞间隙宽度的评估和为新建桥梁提供碰撞间隙宽度 设置的理论依据和参数分析方法。
背景技术：
地震本身具有随机性，所以采用随机振动法进行结构抗震分析是未来抗震 领域发展的方向，但是由于理论复杂，计算效率低下使其在土木工程中的应用 受到了延缓。国内外学者运用随机振动理论对大跨度桥梁进行了相关研究，在 随机振动领域做了重大的贡献，并取得相应成果。但是，对于复杂的大跨度结 构，再加上考虑多点激励和非平稳性等条件，传统的随机振动方法计算结果用 CQC组合计算巨大，很难用于实际工程。很多学者采用简化的SRSS组合方式， 但是这却忽略了相关项的影响，对于频率密集的大跨度复杂结构，振型之间相 互耦合，SRSS组合方式不能满足工程要求。
虚拟激励法采用矩阵表达形式概念明确，计算效率高。虚拟激励法将非平 稳随机振动转化为时程分析，简化了非平稳随机振动的求解方式。但通用有限 元软件没有自带的虚拟激励法计算模块，如何在通用有限元软件中实现随机振 动的虚拟激励法是推动随机振动在工程领域应用的关键，且国内外研究甚少。 一般情况下，需要研究人员自编程序来解决虚拟激励法的计算以及软件的前后 处理，这样就会导致大量时间耗费在编程建模上。传统的虚拟激励法很难在通 用有限元中实现，即使能实现也仅限于自由度比较少的结构。因为传统的虚拟 激励法将绝对位移分解为拟静力位移项和动力相对位移项后求解运动方程。在 一维多点激励时，求解结构的拟静力位移项和动力相对位移项都要先得到静力 影响矩阵，然后再根据静力影响矩阵确定拟静力位移和自由度方向的虚拟激励 荷载来求动力相对位移项，尤其是对于复杂结构，自由度相当多，这使得传统 虚拟激励荷载加载和求解非常耗时，更不用说多点非平稳激励分析了。如果还 需自己编写程序计算，那么虚拟激励法很难在实际工程中进行推广应用。
众所周知，时程分析法也是非常耗时。虽然虚拟激励法将非平稳随机振动 转化为时程分析，但是在每一个离散频率点做一次时程分析，那么离散频率有 成千上万个，导致时程分析也有成千上万次，计算量还是很大。精细积分法的 引入正好解决了这个问题，只需要两个时程分析就能完成所有的计算，这大大 提高了计算效率，便于工程应用推广。
另外，地震灾后调查发现，地震作用下桥梁设计提供的梁体间隙宽度不足 时，不能满足梁体之间的相对位移需求，从而导致梁体发生碰撞破坏。在频遇 地震中，梁体之间的碰撞导致梁体局部破坏。在罕遇地震中，梁体之间的碰撞 可能导致落梁，落梁过程中梁体可能会对桥墩有巨大冲击，造成全桥的垮塌， 进而使得灾后桥梁修复困难。如何提供合理的间隙宽度来减轻或避免梁体在地 震中的碰撞已成为桥梁抗震领域研究的关键问题。
再者，基于概率理论来研究桥梁碰撞间隙宽度的设计能表征地震动的随机 性对间隙宽度的影响。国内外学者对碰撞间隙宽度的研究较少，研究方法基本 为确定性动力时程分析和非确定性的随机振动两种方法。动力时程分析方法要 体现随机性需通过大量的时程积分运算才能得到相对位移响应峰值均值。如要 得到桥梁单个构件破坏的条件超越概率，需要计算不同地震强度水平下的峰值 响应均值和方差，计算量巨大，更不用说考虑桥梁整个系统的条件超越概率； 随机振动法可方便得到相对位移响应峰值均值，但是因其理论复杂，计算效率 低下，在实际工程中应用较少。基于随机振动方法研究各种因素对桥梁碰撞间 隙宽度设计的影响，主要因素有土层分布、地震空间性、地震非平稳性、桩-土 -结构相互作用、支座参数(刚度和阻尼)和地震动强度等。但是，上述的研究均 未涉及到桥梁系统碰撞超越概率。桥梁结构是一个完整的系统，从系统可靠度 角度来专门研究桥梁碰撞间隙宽度是非常必要的。
因此，如果将随机振动的虚拟激励法、精细积分法与矩阵系统可靠度结合 起来研究桥梁系统在各设计间隙宽度下的碰撞系统概率分布，进而就可以评估 间隙宽度设置对桥梁系统安全性的影响。
技术实现要素：
本发明的目的是要解决现有技术问题，从而提供一种非平稳地震作用下桥 梁系统碰撞间隙宽度的概率评价方法。
为达到上述目的，本发明是按照以下技术方案实施的：
非平稳地震作用下桥梁系统碰撞间隙宽度的概率评价方法，包括以下步骤：
步骤一、建立非平稳空间地震激励下的桥梁碰撞间隙宽度需求概率模型；
步骤二、引入精细积分法提高非平稳随机振动计算效率；
步骤三、建立了非平稳地震作用下震级与梁体间隙宽度超越概率的关系；
步骤四、基于矩阵系统运算方法建立桥梁系统所有碰撞处各碰撞间隙宽度 和桥梁系统碰撞概率关系。
具体地，所述步骤四中的所有碰撞处为梁体与梁体和梁体与桥台之间的碰 撞处。
进一步，所述步骤二采用方法是基于APDL语言二次开发功能及Matlab数 据处理功能实现非平稳地震多点激励快速模拟和数据的快速处理，运用直接求 解的虚拟激励法与精细积分法来计算两个瞬态分析得到所有频率点下的响应。
进一步，所述步骤三利用Matlab矩阵运算，通过矩阵迭代计算系数矩阵C 和概率矩阵P，快速计算不同震级下超越桥梁设计的间隙宽度的条件概率。
进一步，所述步骤四将单个间隙宽度超越条件概率按其事件本身的内在联 系进行组合，快速计算出系统条件超越概率，从而识别地震作用下任意超过间 隙宽度后发展碰撞的桥梁系统破坏条件概率。
本发明的原理为：
运用虚拟激励法求解非平稳随机振动，将其转换为时程分析，提高了随机 振动计算效率；精细积分法将时程分析缩减为两次，大大提高计算效率；引入 矩阵运算法再次提高计算效率。系统可靠度计算效率的提高为其在实际工程中 应用奠定了基础；基于虚拟激励法得到不同震级下随时间变化的梁体响应位移 峰值均值和均方值，可得到不同震级下的概率分布，然后通过矩阵迭代计算得 到系统发生碰撞的条件概率分布，进而评估间隙宽度设置对桥梁系统安全性的 影响
1.精细积分法将时程分析缩减为两次，大大提高计算效率；
引入矩阵运算法再次提高计算效率。系统可靠度计算效率的提高为其在实 际工程中应用奠定了基础。基于虚拟激励法得到不同震级下随时间变化的梁体 响应位移峰值均值和均方值，可得到不同震级下的概率分布，然后通过矩阵迭 代计算得到系统发生碰撞的条件概率分布，进而评估间隙宽度设置对桥梁系统 安全性的影响。
与现有技术相比，本发明具有以下有益效果：
1、本发明选择虚拟激励法求解随机振动，将计算效率提高了2-4个数量级；
2、多点激励的虚拟激励法在通用有限元软件中的快速模拟，不需要编制专 门程序，为工程师节约了时间，赢得了效益，从而推动了随机振动法在实际工 程中的应用；
3、运用虚拟激励法和精细积分法在频域-时域内混合求解得出各强度水平 非平稳地震下结构响应的均值和均方差(均值和均方差随时间变化)；
4、引入Davenport理论计算结构响应峰值均值和均方差，迅速得到地震强 度与平均需求之间的关系；
5、引入极值分布函数获得随时间变化的单个间隙宽度需求的条件超越概 率分布；
6、运用矩阵运算方法计算各极限状态的系统条件超越概率。
综述，本发明运用虚拟激励法、精细积分法和矩阵运算法三次提高系统条 件概率的计算效率，克服了系统概率求解低下的缺点，便于在实际工程中推广 应用；该发明可以从桥梁系统层面。评估桥梁各间隙宽度设置对地震作用下桥 梁动力响应的影响
附图说明
图1为本发明的变化场地条件下的空间分布图。
图2为本发明的时程分析的三角加载形式图。
图3为本发明的流程图。
具体实施方式
下面结合具体实施例对本发明作进一步描述，在此发明的示意性实施例以 及说明用来解释本发明，但并不作为对本发明的限定。
如图3所示，本发明的非平稳地震作用下桥梁系统碰撞间隙宽度的概率评 价方法，包括以下步骤：
步骤一、建立非平稳空间地震激励下的桥梁碰撞间隙宽度需求概率模型；
步骤二、引入精细积分法提高非平稳随机振动计算效率；
步骤三、建立了非平稳地震作用下震级与梁体间隙宽度超越概率的关系；
步骤四、基于矩阵系统运算方法建立桥梁系统所有碰撞处各碰撞间隙宽度 和桥梁系统碰撞概率关系。
如下以下对本发明进行详细计算，如图3所示，具体步骤如下：
1)绝对位移直接求解结构动力方程的虚拟激励法
运用虚拟激励法求解式(1)可得：
式中，Peiωt就是虚拟激励荷载，P是由激励功率谱分解得到，为结 构绝对位移虚拟响应。求解方程(2)式，得到绝对位移响应功率谱密度如下
第i阶谱矩阵表示为
基于Davenport理论得到峰值响应的期望和均方差分别为
峰值期望：
峰值均方差：式(5)和式(6)中，ye和σy分别为线性响应X(t)的峰值和标准差，并定义无量纲参 数η＝ye/σy；平均频率表示为
E(η)为η的近似峰值期望值表示 为E(η)≈(2lnvτ)1/2+γ/(2lnvτ)1/2；τ为地震持时，γ＝0.5772为欧拉常数。
基于式(1-6)可得到不同地震强度水平作用下桥梁梁体之间和梁体与桥台 之间的相对位移(△)的响应峰值，从而得到地震作用下梁体间的间隙宽度需求均 值和均方差，为超越碰撞间隙宽度的条件概率需求模型奠定基础。
2)碰撞间隙宽度条件碰撞概率模型
据研究发现相对位移峰值响应Xpeak(t)为随机过程且服从极值I型分布。 Xpeak(t)概率分布函数可表示为
G(Xpeak(t))＝exp{-exp{-αn(Xpeak(t)-un)}} (7)
式中
不同震级下的相对位移响应峰值均值和均方差由式(5)和式(6)得到并表示为 各震级条件下峰值均值和均方差分别为和那么在特定某一震级 下超过间隙宽度并使得桥梁发生碰撞的概率被称为条件碰撞概率。条件碰撞超 越概率表示如下
Pp/m(t)＝1-exp{-exp{-αn(Xpk/m(t)-un)}}(8) 式中
从式(8)中可知，震级条件下峰值均值和均方差已知情况下， Xpk/m(t)的分布可以判断任意间隙宽度时的条件碰撞超越概率。
3)基于矩阵的系统可靠度
假设桥梁结构第i个碰撞点有di个设计间隙宽度，i＝1,...,n，则样本空间被分 为个穷尽互斥基本事件，且用ej表示，j＝1,...,m。c表示任意事件系数 向量且由0和1组成。ej发生的概率可表示为pj＝P(ej)，j＝1,...,m。由于事件ej互斥，那么非平稳地震作用下系统事件的概率Esys(t)可表示为ej事件之和，则系 统条件概率表示为两向量的内积
式中p为概率矩阵的列向量。
桥梁系统每一个碰撞点在多个设计间隙宽度的情况下，系统碰撞事件概率 可表示为Psys(t)＝CTP(t)，Psys(t)表示一个随时间变化的概率矩阵，其中i行j列的 元素表示在第j个条件下第i个事件在时间t的概率。只要确定C和P(t)中各元 素，则碰撞系统事件概率可方便得到。系数矩阵C通过以下迭代方式得到
C[n]的第i列是第i个事件向量，记为有如下的运算关系
同理，p也可由矩阵迭代计算得到
其中，Pi表示第i个事件概率，且
步骤1)中式(2)的虚拟力P由激励功率谱分解得到且表征地震动空 间性。
式中，子矩阵元素Si'j'(iω,t)为3×3维矩阵，参照附图1，相应两个水平向(x，y) 和竖向(z)的地震动分量组成的互功率谱矩阵表示为
三维地面运动的水平向(x,y)分量取相同的功率谱密度函数，而竖向(z)的功率谱 密度取其水平向的0.6倍。故式(15)的非对角线元素为
式中，Si'j',xx(iω,t)，Si'j',yy(iω,t)，Si'j',zz(iω,t)由场地传递函数从基岩功率谱密度矩 阵函数得到。
式中，HU,i(iω)，HV,i(iω)，HW,i(iω)和HU,j(iω)，HV,j(iω)，HW,j(iω)分别为空间 点i-i′、j-j′在U，V，W方向的场地传递函数。i，j点x，y，z方向的自功率 谱密度函数为
式中，Gi,x(ω,t)为第i个空间点x向在时域和频域内的调制函数，Sii,xx(iω)是
基岩x向自功率谱密度函数，γij(iω)为空间点i，j的相干函数，形式如下
γij(iω)＝|γij(iω)|exp[-iωdij/vapp(ω)] (19)
其中，dij为空间点i，j距离在地震波传播方向的投影，vapp(ω)表示视波速。
为了得到式(16)-(18)中的频率传递函数，基于一维波动理论，假设基岩中地 震波由面外SH波或者面内P-SV波组成，场地频率传递函数HU,i(iω)，HV,i(iω)， HW,i(iω)表示为
KSHXSH＝PSH，KP-SVXP-SV＝PP-SV (20)
其中，XSH，PSH和XP-SV，PP-SV是关于SH波和P-SV波的位移和荷载向量；刚度 矩阵XSH和KP-SV主要由土体性质、入射波类型、入射角和圆频率ω确定；动力荷 载向量PSH和PP-SV主要依赖于基岩的特性、类型、频率和入射波的幅值。通过式(20) 可以得到在每一个离散频率点面内或面外的频率传递函数。
为了构造虚拟力，三维非平稳地震激励功率谱矩阵可分解为
式中，P是3m×r维矩阵，r是矩阵的秩，上标*和T分别表示复数共 轭和转置。
式(20)中的W(iω)为3m×3m矩阵，表示了行波效应和场地土条件对相位的 变化，表达式如下
其中，Ti,x表示地震波沿x方向到第i个支撑处所需时间，θi,x，θm,y，θm,z分别表 示在第i个支撑处x，y，z方向的相位角，将由下面的公式得到
式中，θij,xx为空间点i-i′和j-j′在x方向的相位差。
式(11)中的Sg(iω)表达式如下
式中，Sgi,x代表空间点i的x方向上的功率谱密度函数，水平向(x,y)与竖向(z) 的功率谱比值为3:2。
式(11)中的G(ω,t)为三维非平稳地面运动非均匀调制函数，形式如下
一般认为各个空间每一个方向上调制函数一致，故可取
Gi,x(ω,t)＝Gi,y(ω,t)＝Gi,z(ω,t)＝G(ω,t) (27)
式(11)中，矩阵R表示多维多点地面运动的相干矩阵，一般是正定或半正定 矩阵，可通过LDLT分解如下
R＝[Q]3m×r[QT]r×3m＝[|γij|]3m×3m(当i＝j时，|γij|＝1；当i≠j时，0≤|γij|<1) (28)
式中，R被分解为非零特征值αj与其相应特征向量(j＝1,2,…,r,r≤3m)的和，其 形式如下
部分相干情况下，第阶虚拟激励表示为
从式(30)可知，相干效应、行波效应、局部场地效应和非平稳性分别通过虚拟激 励项中的W(iω)、Sg(iω)和G(ω,t)来体现。
完全相干和完全不相干虚拟力可表示为
根据式(20)和式(21)可得第j个特征值下的绝对位移响应为
式中
式中，Ij(iω,t)是杜哈梅尔积分，h(t-τ)单位脉冲函数，在部分相干情况的绝对 位移响应功率谱密度函数表示为
同样，可以得到完全相干和完全不相干情况下的响应为
完全相干
式中
完全不相干
式中
由公式(33)、(36)、(38)可以看出，需要在每一个频率点下做一次瞬态分 析计算，虽然相对于传统的随机振动理论来说，直接求解的虚拟激励法已经非 常节约时间了，但是如果离散的频率点太细，那么计算的瞬态分析也就越多， 这样也就非常耗时，故下面引入精细积分方法并在通用有限元软件实现，这样 只需要两个瞬态分析就能确定结构任意一个关注节点的响应，这样大大缩短了 计算时间。
在有限元计算中，式(2)可另写成如下形式
式中
为确定频率点下的非平稳随机激励时程，Mb为桥梁基础处附加的大质量。
将每一个确定频点下的动力方程式(39)写成状态方程形式
其中
状态方程(41)的一般解为
式中eHt为指数矩阵，对式(43)进行离散积分，设时间步长为Δt，由递推法 可知，在ti时刻的响应V(ti)＝Vi可用ti-1时刻的响应V(ti-1)＝Vi-1来表示如下
式中，T＝eH△t为指数矩阵，改写成如下形式
T＝eH△t＝[eH△t/m]m (45)
令τ＝△t/m，当取m＝2N很大的，△τ非常小，由泰勒级数展开式有
当N＝20时，与一般计算机的舍入误差相比较，泰勒级数的截断误差要小的多， 不会由于截断带来数值误差，故T(τ)的计算精度实际上已经给出了计算机上的精 确解。
[I+Tai]≡[I+Ta,i-1]2＝[I+2Ta,i-1+Ta,i-1×Ta,i-1](i＝1,2,…,N) (47)
因此依次类推
[I+TaN]＝[I+Ta,N-1]2＝[I+Ta,N-2]4＝…＝[I+Ta0]m＝T(τ) (48)
每次运算单元矩阵[I]不参与计算，因为Tai很小，当它与单位矩阵[I]相加时 就成为其尾数，在计算机的舍入操作中就会自动舍去。为了避免这种情况的发 生，采用下列递推方式求解
[Tai]＝2[Ta，i-1]+[Ta，i-1][Ta，i-1](i＝1,2,…,N) (49)
假定在每一个时间积分步(ti-1,ti)中，荷载为线性变换，则非平稳虚拟力 F(ω,t)可离散为t0,t1,t2,…,tk处的随机变量F0,F1,F2,…,Fk,非平稳虚拟力可表示为
F(ω,τ)＝r0+r1(τ-ti-1) (50)
式中，r0，r1是一个是不变向量。
将式(50)代入式(44)可得
式中H-1的计算可根据式(42)可得
将式(51)进一步写成下式
令S1＝(I-T)(H-2/△t)+TH-1,S2＝(T-I)(H-2/△t)-H-1 (54)
式中，H-2＝H-1·H-1，则式(53)简化为
Vi＝T△Vi-1+S1Fi-1+S2Fi(i＝1,2,…,k) (55)
每一固定频点时，V0＝0，则通过式(55)可得到固定频率点各时刻的响应表达式为
式中
S3＝TS2+S1 (57)
若用Ai,0，Ai,1，…，Ai,i分别表示F0，F1，…，Fi前面的系数，Ai,i中第一 个i表示ti时刻，第二个表示第i个荷载离散点位置，则式(56)可表示为
Vi＝Ai,0F0+Ai,1F1+…+Ai,i-1Fi-1+Ai,iFi(i＝1,2,…,k) (58)
式中，Ai,0，Ai,1，…，Ai,i只与结构本身有关，而与外荷载无关，体现了结构 本身在外荷载激励下所表现出来的固有属性。Ai,i的计算过程是一个递推过程， 用Vi-1的系数Ai-1,0，Ai-1,1，…，Ai-1,i-1来推导下一步的系数为
将式(59)代入式(58)写成矩阵的形式，则
由式(60)可知，系数矩阵中所有其他所有的元素都是第一列和第二列元素的 重复，也就是说只要能得到第一列和第二列的所有元素，那么整个系数矩阵就 能确定。令F0＝1，Fk＝0(k>0)，进行一个瞬态分析就可以得到系数矩阵的第一列 元素。同理，令F0＝0，F1＝1，Fk＝0(k>2)，进行一个瞬态分析就可以得到系数矩 阵的第二列元素，求解加载方式如附图2。
只要得到系数矩阵，那么任意频率点处的非平稳时程荷载作用下的响应就 能根据式(60)得到。当运用绝对位移直接求解时，对于多维多点激励的非平稳地 震动荷载，需要对每一个激励点三个方向依次进行激励，分别得到所需的结构 关键点响应的系数矩阵A，然后将每一个方向的激励时程与对应的系数矩阵A 相乘后叠加就得到多维多点在固定频率点的结构响应。最后，将所有的频率点 下的激励时程都对应代入后就得到整个频率范围的结构响应。
本发明的技术方案不限于上述具体实施例的限制，凡是根据本发明的技术方案 做出的技术变形，均落入本发明的保护范围之内。
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
