![预处理共轭梯度法识别桥梁动荷载分析与应用](https://wfqqreader-1252317822.image.myqcloud.com/cover/618/37204618/b_37204618.jpg)
上QQ阅读APP看本书,新人免费读10天
设备和账号都新为新人
2.2 第二识别法(IMII)
![img](https://epubservercos.yuewen.com/7E8DE4/19720709701116506/epubprivate/OEBPS/Images/txt002_10.jpg?sign=1739288239-zz2tix6kSV7kSmieNqIRVNLBfbGwtatf-0-c3b5103977824cd0c68803d34479ba5e)
图2.2.1 移动荷载识别简支梁模型
如图2.2.1所示,将桥梁考虑为一简支梁,其跨长为L,抗弯刚度为EI,单位长度质量为ρ,黏性比例阻尼为C,忽略剪切变形和转动惯量(即伯努利-欧拉梁)。假设一动荷载f(t)以速度c自梁左端支承处向右移动,则其振动微分方程:
![img](https://epubservercos.yuewen.com/7E8DE4/19720709701116506/epubprivate/OEBPS/Images/txt002_11.jpg?sign=1739288239-GPVDCfUcFwukxkUquZknYakcLolIcdUv-0-3d6ae667372fbd00fc253faea79fdb7e)
这里v(x,t)是梁在时刻t、位置x处的变形,δ(x-ct)是狄拉克函数。
式(2.2.1)的边界条件为
![img](https://epubservercos.yuewen.com/7E8DE4/19720709701116506/epubprivate/OEBPS/Images/txt002_12.jpg?sign=1739288239-ZO7qr91vYTHNByvehoKwu3ihxtMAuAyi-0-b74642db7911f6924b157f4e8b5f67ad)
和
基于模态叠加原理,假设梁的第n阶模态振型函数为则式(2.2.1)的解可表示为
![img](https://epubservercos.yuewen.com/7E8DE4/19720709701116506/epubprivate/OEBPS/Images/txt002_15.jpg?sign=1739288239-q4z6MnZlWWhp5UZytSffWigI874cTQIO-0-a3459170879cdc2666ba2688940b5df0)
矩阵形式为
![img](https://epubservercos.yuewen.com/7E8DE4/19720709701116506/epubprivate/OEBPS/Images/txt002_16.jpg?sign=1739288239-cZMEvMlD0aenIq2NuiPxFnMpvvWFUAGT-0-03d0dc7f1ad9ccf18a14ab36c260d6e9)
这里n为模态数,qn(t)(n=1,2,…,∞)是第n阶模态位移。将式(2.2.2)代入式(2.2.1),并在[0,L]内对x进行积分,利用边界条件和狄拉克函数特性,系统振动微分方程可用模态位移qn(t)表示为
![img](https://epubservercos.yuewen.com/7E8DE4/19720709701116506/epubprivate/OEBPS/Images/txt002_17.jpg?sign=1739288239-M0ZRMNyiAM14PZBTUGRBVy1lliiNChPW-0-5a5d107b61bbd411167ae82b1dbc4d2e)
这里分别为桥梁第n阶模态频率、阻尼比和模态力。
如有k个荷载,且第k个荷载到第一个荷载的距离为则式(2.2.4)可写为
![img](https://epubservercos.yuewen.com/7E8DE4/19720709701116506/epubprivate/OEBPS/Images/txt002_20.jpg?sign=1739288239-VCp2nbpwuzU2fTFyyuXZoicse9oNsdQs-0-72830f09dad2e86247ff3da8f4b1edc6)
x1,x2,…,xl处的模态位移可通过式(2.2.1)求得
![img](https://epubservercos.yuewen.com/7E8DE4/19720709701116506/epubprivate/OEBPS/Images/txt002_21.jpg?sign=1739288239-ieSosV4WAWGHZx4jlksi7OZWLXhsqZZN-0-049de1dd533fa5d843939ba00d6326d0)
梁上x1,x2,…,xl处的速度可通过位移的一次微分求得
![img](https://epubservercos.yuewen.com/7E8DE4/19720709701116506/epubprivate/OEBPS/Images/txt002_22.jpg?sign=1739288239-XT1M5m47NVGDfq0h0icd8LNOosD85awg-0-f4e138889cb6444a92c29ae3c27997b8)
进一步,梁上x1,x2,…,xl处的加速度可通过位移的二次微分求得
![img](https://epubservercos.yuewen.com/7E8DE4/19720709701116506/epubprivate/OEBPS/Images/txt002_23.jpg?sign=1739288239-tZxRIL1pmLVCqO804aIuszbzgWxwjg3B-0-b5ade39a30fd7af41af588650ec97249)
类似地,相应位置的弯矩可利用关系式M=-EI(∂2v/∂x2)求得
![img](https://epubservercos.yuewen.com/7E8DE4/19720709701116506/epubprivate/OEBPS/Images/txt002_24.jpg?sign=1739288239-caPIBmHczBrnkWGkk86qXyqPmXmmfjJq-0-c238b1b55a4bc5839574b3361f984ac4)
若f1,f2,…,fk为已知常量移动荷载,忽略阻尼的影响,则式(2.2.1)的解为
![img](https://epubservercos.yuewen.com/7E8DE4/19720709701116506/epubprivate/OEBPS/Images/txt002_25.jpg?sign=1739288239-cTfuBE5BVfgOmk6SzTnuSqw8If1bQmAf-0-ae18677c1d6a04463c8d0974a10a65f2)
这里若在一组常量移动荷载作用下,x1,x2,…,xl处的位移已知,则每个常量移动荷载可通过解下式方程求得
![img](https://epubservercos.yuewen.com/7E8DE4/19720709701116506/epubprivate/OEBPS/Images/txt002_27.jpg?sign=1739288239-JX00CmmGhcTDxjPP7LYgyVLUbYMXtVmY-0-326c6e59d2c8f3ec857316b5a336295a)
其矩阵形式可表示为
![img](https://epubservercos.yuewen.com/7E8DE4/19720709701116506/epubprivate/OEBPS/Images/txt002_28.jpg?sign=1739288239-d64fNIaW2FvKDChjrHIHXCyeOqBGxoaz-0-958b1243c404b85afdb367899e21fdcb)
这里
![img](https://epubservercos.yuewen.com/7E8DE4/19720709701116506/epubprivate/OEBPS/Images/txt002_29.jpg?sign=1739288239-jysqVbUlSofghFr4eOsj5zgIsZqRqq0M-0-184590a55190f7ef76d4c3cc51b435f7)
若l≥k,即位移的测量点数大于或等于车轴轴数,f可用最小二乘法求解:
![img](https://epubservercos.yuewen.com/7E8DE4/19720709701116506/epubprivate/OEBPS/Images/txt002_30.jpg?sign=1739288239-RQFq3pVs7kZikOFFOAEr7MGRF8S91KFP-0-b4c1c476f4b9e535bbf170f6dc1eb84f)
若已知的不是桥梁位移响应,而是弯矩响应,则同样可以从弯矩响应求得式(2.2.1)的解:
![img](https://epubservercos.yuewen.com/7E8DE4/19720709701116506/epubprivate/OEBPS/Images/txt002_31.jpg?sign=1739288239-TD0hNeWgYZZuPGZfcVPv9o0KuumR36Fj-0-dbddb5ef392576ecd029a46732566096)