基于状态空间法的圆弧曲梁面内动力学分析.pdf

文章编号：1000-4750(2023)Suppl-0019-06基于状态空间法的圆弧曲梁面内动力学分析刘兴喜1，杨博1，徐荣桥2(1.浙江理工大学土木工程系，浙江，杭州310018；2.浙江大学土木工程系，浙江，杭州310058)摘要：该文采用状态空间法研究了圆弧曲梁面内动力学特性。通过选择合适的状态变量，建立了相应的状态空间列式，并给出了固有频率和振动模态的求解过程。进一步通过引入辛内积的概念，建立了三种工程常见边界条件(简支、固支和自由)下圆弧曲梁面内振动模态关于质量和旋转惯量的正交关系式。在此基础上，运用模态叠加法给出了非齐次状态方程的解析解，并得到了圆弧曲梁在竖向移动集中荷载作用下的瞬态响应。数值算例结果表明该文方法是十分精确和可靠的。关键词：圆弧曲梁；自由振动；强迫振动；状态空间法；辛内积；模态正交性中图分类号：TB122文献标志码：Adoi:10.6052/j.issn.1000-4750.2022.06.S018IN-PLANEDYNAMICANALYSISOFCIRCULARLYCURVEDBEAMSBASEDONTHESTATESPACEMETHODLIUXing-xi1,YANGBo1,XURong-qiao2(1.DepartmentofCivilEngineering,ZhejiangSci-TechUniversity,Hangzhou,Zhejiang310018,China;2.DepartmentofCivilEngineering,ZhejiangUniversity,Hangzhou,Zhejiang310058,China)Abstract:Thestatespacemethodisadoptedtostudythein-planedynamicbehaviorofcircularlycurvedbeams.By selecting appropriate state variables,the state space formula is properly set up.The frequencies andcorrespondingmodalshapesofcircularlycurvedbeamsarethenobtainedonthebasisofstate-spaceformulae.Byutilizingtheconceptionofsymplecticinnerproduct,themodeorthogonalitywithrespecttothemassandrotaryinertiapropertiesforcurvedbeamsunderthreecommonboundaryconditionsinengineering(simply-supported,clampedandfree)isestablished.Basedontheorthogonalityrelation,themodesuperpositionmethodisutilizedtoobtainthesolutionoftheinhomogeneousequationforforcedvibrationanddynamicresponseofcurvedbeamssubjectedtoaverticalmovingconstantforce.Thenumericalresultsshowthattheproposedmethodisveryaccurateandreliable.Keywords:circularly curved beams;free vibration;forced vibration;state space method;sympletic innerproduct;modalorthogonality由于曲梁结构具有造型优美、施工简便及受力性能好等优点，已被广泛应用于土木工程、航空航天、工业机械等领域1。在实际工程中，曲梁不可避免地会受到动力荷载的作用(如地震作用、风荷载、移动荷载等)，动力学分析是这类结构计算分析的重要内容。与直梁相比，由于受到曲率的影响，曲梁的动力特性很复杂，目前缺乏对其进行较为精确分析的方法。因此，建立一套精确、可靠的计算方法，高效地获得曲梁动力学问题的精确解，具有十分重要的意义。DENHARTOG2最早使用 Rayleigh-Ritz 法得到 了 圆 弧 拱 面 内 振 动 的 频 率。BICKFORD 和收稿日期：2022-06-27；修改日期：2022-11-25基金项目：国家自然科学基金项目(12272335)通讯作者：徐荣桥(1972)，男，浙江人，教授，博士，从事组合桥梁结构、FRP加固分析理论等方面的研究(E-mail:).作者简介：刘兴喜(1989)，男，湖北人，讲师，博士，从事 FRP 加固分析理论、结构动力分析理论等方面的研究(E-mail:);杨博(1979)，男，宁夏人，教授，博士，主要从事复合材料力学与结构方面的研究(E-mail:).第40卷增刊Vol.40Suppl工程力学2023 年 6月June2023ENGINEERINGMECHANICS19STROM3运用传递矩阵法对曲梁面内和面外自由振动进行了研究。基于精确的动力刚度法，HOWSON和 JEMAH4、TSENG 等5分析了曲梁的面内自振特性。叶康生和殷振炜1、YANG 等6和袁驷等7采用有限元这一应用范围广泛的数值方法研究了曲梁的面内自由振动。虽然国内外学者对曲梁的自由振动特性做了大量的研究，关于其强迫振动的研究则相对不足。吴玉华和楼文娟8采用虚拟激励法与 Galerkin 法对水平随机地震作用下弹性圆拱的随机响应进行了研究。EFTEKHARI9运用基于微分求积法得到了移动集中荷载作用下曲梁的瞬态响应。但文献 8 和文献 9 都是基于 Euler曲梁理论，没有考虑剪切变形和转动惯量的影响。状态空间法通过选择能量对偶的两类物理量如位移和内力作为状态变量，能将描述结构力学特性的控制方程整理成简洁的矩阵形式10，即状态方程。这种方法具有计算效率高、易于编程等优点，已被广泛地应用于结构的静力和动力分析中1113。然而，在曲梁的强迫振动分析中尚鲜见报道。本文基于 Timoshenko 曲梁理论，通过选择合适的状态变量，建立了圆弧曲梁面内动力学问题的状态空间列式，并给出了固有频率和振动模态的求解过程。采用 SHEN 等11和 LIU 等13的方法，即通过引入辛内积，建立了普通边界条件下圆弧曲梁振动模态关于质量和旋转惯量的正交关系式。基于所建立的正交关系式，运用模态叠加法导出了非齐次状态方程即强迫振动方程的解析解，并得到了圆弧曲梁在竖向移动荷载作用下的瞬态响应。最后，给出了数值算例，结果表明采用本文方法进行圆弧曲梁面内动力学分析是十分精确和可靠的。1状态方程列式lRsszuwqtqnqm考虑任意动力荷载作用下的圆弧曲梁，如图 1所示，其中曲线长度和半径为分别为 和。沿曲梁中曲线建立坐标轴 轴，将 轴逆时针旋转 90得到法向坐标轴 轴。、分别是轴线上的环向位移和径向位移，为截面的转角，、和分别为作用在曲梁轴线上的环向分布荷载、径向分布荷载以及分布弯矩。曲梁的微元及其横截面示意图如图 2 所示。MVNMAIA=bhI=bh3/12bhEGG=E/2(1+)、分别为整个截面上的弯矩、剪力和轴力，弯矩正方向的选取是为了保证弯矩在相应的正转角上做正功；、分别为截面面积和截面惯性矩，并且有，、分别为截面的宽度和高度；、分别是密度、弹性模量和剪切模量，并且有，为泊松比。根据微元的平衡，可以得到：Ns+VR=A2ut2qt,VsNR=A2wt2qn,Ms+V=I2t2qm(1)lRz,ws,uqn(s,t)qm(s,t)qt(s,t)图1平面圆弧曲梁模型Fig.1PlanarcircularlycurvedbeammodelVNV+dVN+dNqn(s,t)qt(s,t)qm(s,t)Au,ttA,ttAw,ttMM+dMbh,E,G,I,A图2曲梁微元及其横截面Fig.2Infinitesimalelementandtypicalcross-sectionofacurvedbeamMVN根据 Timoshenko 曲梁基本理论14，弯矩、剪力 和轴力可以表示为：M=EIs,V=GA(wsuR),N=EA(us+wR)(2)式中，为截面剪切修正系数。将式(1)和式(2)整理成矩阵的形式，即可得到曲梁面内动力学问题的状态方程：vs=Av+q1q2(3)v=wuVNMTTAq1q2式中：为状态向量，上标 为矩阵的转置；系统矩阵、惯性力列向量、荷载列向量的表达式分别为：20工程力学A=01/R11/(GA)001/R0001/(EA)0000001/(EI)00001/R00001/R00000100(4)q1=0 0 0 A2wt2A2ut2I2t2T(5)q2=000qnqtqmT(6)v=wuVNMT本文在这里之所以选择为状态向量：一方面是因为它与边界条件和连续性条件直接相关；另一方面在后续的模态关于质量和旋转惯量的正交关系式证明中，为了利用辛内积的概念，必须将三个广义力和相应的三个广义位移形成一一对应的共轭物理量，即剪力对应于径向位移、轴力对应于环向位移、弯矩对应于转角，并要使得广义力在相应的广义位移上做正功11,13。2自由振动分析2.1状态方程的解对于曲梁面内自由振动问题，可假设：v=w(s)u(s)(s)V(s)N(s)M(s)eit=v(s)eit(7)v(s)vq2=0式中：为固有频率(圆频率)；为 的振幅向量。将式(7)代入式(3)，并且令，得到：dvds=Hv(8)式中：H=01/R11/(GA)001/R0001/(EA)0000001/(EI)A20001/R00A201/R0000I2100(9)H如果曲梁的几何、材料属性沿轴线方向是不变的，即系统矩阵的元素都是常数，则可采用传统传递矩阵法对状态方程式(8)进行求解，得到其解为：v(s)=eHsv(0)(10)eHs式中，为传递矩阵。2.2固有频率及振动模态s=l在式(10)中令，得到：v(l)=Tv(0)(11)T=eHl式中，为曲梁两端状态向量间的传递矩阵。对于两端简支的曲梁，其边界条件可以表示为：w(0)=u(0)=M(0)=w(l)=u(l)=M(l)=0(12)将式(12)代入式(11)，可以得到：T13T14T15T23T24T25T63T64T65(0)V(0)N(0)=TSS(0)V(0)N(0)=000(13)TijTij式中，为矩阵第 行第 列的元素。式(13)是一个齐次线性方程组，有非零解的条件是其系数矩阵行列式的值为 0，即：|TSS|=0(14)s=0式(14)为两端简支曲梁的频率方程，它是一个关于频率 的超越方程，可以采用二分法来搜索它的根15。一旦求得某一阶频率后，将其代回式(13)，即可求得在处的状态向量(包含一个任意常数因子)，再通过式(10)可确定曲梁内任意位置处的状态向量，这样就可以确定该频率对应的振动模态(也包含一个任意常数因子)。对于其他边界条件，同样能求得相应的频率和振动模态。2.3模态关于质量和旋转惯量的正交关系式vivj为了运用模态叠加法计算曲梁在动力荷载作用下的响应，必须建立振动模态关于质量和旋转惯量的正交关系式。本节将采用 SHEN 等11和 LIU等13的方法，即基于状态空间列式和辛内积的概念来进行推导。首先考察对应于两个不同模态的状态向量 和，定义如下辛内积：(vi,vj)wl0vTiJvjds=wl0(wiVj+uiNj+iMj)dswl0(Viwj+Niuj+Mij)ds(15)J=0I33I330I33JT=JJ2=I66I66式中，为单位辛矩阵，为三阶单位矩阵。可以证明单位辛矩阵具有如下性质：，是六阶单位矩阵。wiwj状态向量和对应于两个不同频率的振动工程力学21模态，因此有：dvkds=Hkvk,k=i,j(16)式中：Hk=01/R11/(GA)001/R0001/(EA)0000001/(EI)A2k0001/R00A2k01/R0000I2k100(17)ijvTjJvTiJ对于第 阶模态和第 阶模态，将式(16)两边同时左乘和，并在曲梁的整个弧长上进行积分，然后将两式相减，可以得到：wl0(vTjJdvidsvTiJdvjds)ds=wl0(vTjHivivTiHjvj)ds=(2j2i)wl0(Awiwj+Auiuj+Iij)ds(18)JT=J对式(18)左边进行分部积分并利用，可以得到：wl0(vTjJdvidsvTiJdvjds)ds=(vTjJvi)?l0=(Viwj+Mij+NiujVjwi+Mji+Njui)?l0(19)由式(18)和式(19)可以得到：(2j2i)wl0(Awiwj+Auiuj+Iij)ds=(vTjJvi)|l0=(Viwj+Mij+NiujVjwi+Mji+Njui)|l0(20)考虑实际工程中常见的 3 种边界条件，即简支(S)、固支(C)和自由(F)，它们可以分别表示为：S:w=u=M=0C:w=u=0F:V=M=N=0(21)将这些边界条件代入式(21)，可以得到：wl0(Awiwj+Auiuj+Iij)ds=0,i,j(22)式(22)就是曲梁的两个不同模态关于质量和旋转惯量的正交关系式。3强迫振动分析3.1模态叠加法本节将采用模态叠加法研究曲梁的强迫振动问题，为此，假设：v(s,t)=W(s)(t)(23)W(s)(t)式中，、分别为模态矩阵和广义坐标矢量，可分别定义为：W(s)=v1(s)v2(s)vn(s),(t)=1(t)2(t)n(t)T(24)wk(s)k式中，为第 阶模态。将式(24)的第一式代入式(8)，可以得到：dW(s)ds=AW(s)M(s)(25)M(s)6n式中，是阶矩阵，并且只有第 4 行、第5 行、第 6 行的元素是非零的，即：M(4,k)=A2kwk,M(5,k)=A2kuk,M(6,k)=A2kk,k=1,2,n(26)将式(23)代入式(3)，得到：dW(s)ds(t)=AW(s)(t)+V(s)d2(t)dt2q2(s,t)(27)V(s)6n式中，为阶矩阵，并且只有第 4 行、第5 行、第 6 行的元素是非零的，即：V(4,k)=Awk,V(5,k)=Auk,V(6,k)=Ak,k=1,2,n(28)将式(25)代入式(27)，可以得到：V(s)d2(t)dt2+M(s)(t)=q2(s,t)(29)W(s)TJs将式(29)两边同时乘以，并在曲梁整个弧长上对 进行积分，得到：ni=1wl0(Awkwi+Aukui+Iki)dsd2i(t)dt2+2ii(t)=wl0qn(s,t)wk+qt(s,t)uk+qm(s,t)kds(30)利用模态关于质量和旋转惯量的正交关系式(22)，得到：d2k(t)dt2+2kk(t)=AkBk,k=1,2,n(31)式中：Ak=wl0qn(s,t)wk+qt(s,t)uk+qm(s,t)kds(32)Bk=wl0(Awkwk+Aukuk+Ikk)ds(33)k(t)关于的微分方程的通解可用杜哈梅积分表示为：22工程力学k(t)=k(0)cos(kt)+sin(kt)kdk(0)dt+1kwt0AkBksink(t)d,k=1,2,n(34)k(0)dk(0)dtk(0)=dk(0)dt=0式中，、为曲梁的初始条件，若初始时曲梁处于静止状态，则有。3.2竖向移动集中荷载作用下的动力响应Pm0ctd=l/c考虑曲梁上作用一个沿着梁轴线移动的竖向集中荷载，大小为，方向竖直向上为正，竖直向下为负，移动速率为，为移动荷载作用在梁轴线上的时间，那么该动力荷载可以表示为：qn(s,t)=Pm0(sct)1+y(ct)2,qt(s,t)=Pm0(sct)y(ct)1+y(ct)2(35)()y(ct)s=ct式中：为 DiracDelta 函数；为处曲梁轴线的斜率。将式(35)代入式(32)，得到：Ak=Pm0wk(ct)1+y(ct)2+Pm0y(ct)uk(ct)1+y(ct)2(36)0t tdk(0)=dk(0)dt=0k(t)当时，移动荷载仍然作用在曲梁上，且初始时曲梁处于静止状态，即初始条件为，那么由式(34)可以得到的表达式为：k(t)=Pm0kBkwt0wk(c)+y(c)uk(c)1+y(c)2sink(t)d,0t td(37)ttd当时，移动荷载已经移出曲梁外，此时的k(td)dk(td)dt解是由初始条件和引起的自由振动，即：k(t)=k(td)cos(kt)+sin(kt)kdk(td)dt,ttd(38)这样，可以分时段求得圆弧曲梁在竖向移动荷载作用下的瞬态响应：v(s,t)=kPm0vk(s)kBkwt0wk(c)+y(c)uk(c)1+y(c)2sink(t)d,0t td(39)v(s,t)=kk(td)cos(kt)+sin(kt)kdk(td)dt vk(s),ttd(40)4数值算例4.1自由振动l=1.1781mR=0.75mE=70GPa=0.85G/E=0.3A=4m2I=0.01m4=2777kg/m3l2A/EIl2A/EI考虑四种不同边界条件下(SS：两端简支；CC：两端固支；CS：一端固支一端简支；CF：一端固支一端自由)的等截面圆弧曲梁，其轴线为 1/4 圆弧，曲梁的基本参数如下：，。表 1 给出了 SS 圆弧曲梁的前 10 阶无量纲频率值，并将其与 YANG 等6、EISENBERGER 和 EFRAIM16的无量纲频率值进行对比。从表 1 可以看出，本文方法计算结果与文献 6、文献 16 计算结果吻合得很好，特别是几乎与文献 16 的结果完全相同，从而验证了本文方法在进行自由振动分析时的正确性。表2 还给出了CC、CS 和CF 这三种边界条件下圆弧曲梁的前 10 阶无量纲频率。l2A/EI表1SS 圆弧曲梁前 10 阶无量纲频率值l2A/EITable1Thefirst10dimensionlessfrequenciesoftheSScircularlycurvedbeam频率阶数本文文献6文献16频率阶数本文文献6文献161阶29.28029.30629.2806阶143.618143.679143.6182阶33.30533.24333.3057阶156.666156.629156.6663阶67.12467.12367.1248阶190.477190.596190.4774阶79.97179.95079.9719阶225.361225.349225.3615阶107.851107.844107.85110阶234.524234.809234.5244.2竖向移动集中荷载R=2mb=0.1mh=0.05mE=200GPa=5/6=0.3=7800kg/m3Pm0=1kNc=5m/s考察竖向移动荷载作用下圆弧曲梁的动力响应，曲梁轴线为 1/2 圆弧，其基本参数如下：，。竖向移动荷载为，移动速度为。表 3 给出了采用状态空间法和模态叠加法计算 SS、CC 和 CS 这三种边界条件下圆弧曲梁跨中径向最大位移的收敛性，从表 3可以看出，取前 10 阶模态完全可以满足精度要求。EFTEKHARI9给出的 SS 和 CC边界条件下曲梁跨中径向最大位移分别为0.73mm和0.45mm，与本文的结果吻合得很好，从而验证了运用本文方法进行强迫振动分析时的正确性。图 3 给出了三种边界条件下跨中径向位移和环向位移随时间工程力学23的变化曲线，从图 3 中可以看出，SS 边界条件下跨中位移最大，而 CC 边界条件下跨中位移最小。5结论针对圆弧曲梁面内动力学问题，本文基于所建立的状态空间列式给出了固有频率和振动模态的求解过程。根据 SHEN 等11和 LIU 等13的方法，通过引入辛内积，建立了普通边界条件下圆弧曲梁振动模态关于质量和旋转惯量的正交关系式。基于所建立的正交关系式，运用模态叠加法给出了强迫振动方程的解析解，并具体考察了圆弧曲梁在竖向移动集中荷载作用下的瞬态响应。得到主要结论如下：(1)本文得到的圆弧曲梁自由振动频率、移动荷载下的强迫振动响应与已有结果吻合得很好，证明本文方法是十分精确和可靠的。(2)本文方法很容易扩展到几何、材料属性沿梁轴线方向是变化的曲梁面内动力学分析中，可以为研究复杂结构的振动问题提供一种有效的分析工具。参考文献：叶康生,殷振炜.平面曲梁面内自由振动有限元分析的p型超收敛算法J.工程力学,2019,36(5):2836,52.YEKangsheng,YINZhenwei.Ap-typesuperconvergentrecoverymethodforFEanalysisofin-planefreevibrationofplanarcurvedbeamsJ.EngineeringMechanics,2019,36(5):2836,52.(inChinese)1DENHARTOGJP.XL.Thelowestnaturalfrequencyofcircular arcsJ.The London,Edinburgh,and DublinPhilosophical Magazine and Journal of Science,1928,5(28):400408.2BICKFORDWB,STROMBT.VibrationofplanecurvedbeamsJ.Journal of Sound and Vibration,1975,39(2):135146.3HOWSON W P,JEMAH A K.Exact dynamic stiffnessmethod for planar natural frequencies of curvedTimoshenkobeamsJ.JournalofMechanicalEngineeringScience,1999,213(7):687696.4TSENGTP,HUANGCS,LINCJ.Dynamicstiffnessanalysis for in-plane vibrations of arches with variablecurvatureJ.Journal of Sound and Vibration,1997,207(1):1531.5YANGF,SEDAGHATIR,ESMAILZADEHE.Freein-plane vibration of general curved beams using finiteelementmethodJ.JournalofSoundandVibration,2008,318(4/5):850867.6(参考文献 716 转第 55 页)l2A/EI表2CC、CS 和 CF 三种边界条件下圆弧曲梁前 10 阶无量纲频率值l2A/EITable2Thefirst10dimensionlessfrequenciesofthecircularlycurvedbeamsunderthreedifferentboundaryconditions(CC,CSandCF)频率阶数CCCSCF频率阶数CCCSCF1阶35.07332.2183.6306阶147.261145.802114.5192阶41.14936.75815.5357阶160.291158.232123.5813阶79.58273.45840.2458阶195.774193.194162.2774阶79.99479.99354.4149阶225.365225.363188.3045阶117.748112.96982.70710阶237.933236.264204.282表3竖向移动集中荷载作用下圆弧曲梁跨中径向最大位移计算的收敛性/mmTable3Convergenceofthemaximumradialdeflectionatmid-spanofthecircularlycurvedbeamsunderverticalmovingconcentratedload边界条件模态阶数1阶2阶3阶4阶5阶6阶7阶8阶9阶10阶SS0.00 0.67 0.67 0.71 0.71 0.72 0.72 0.73 0.73 0.73CC0.00 0.41 0.41 0.43 0.43 0.45 0.45 0.45 0.45 0.45CS0.25 0.60 0.62 0.64 0.64 0.65 0.65 0.66 0.66 0.660.0SSCCCS0.01.00.01.21.5时间/sSSCCCS0.90.60.30.01.21.5时间/s0.90.60.31.51.00.50.51.5跨中环向位移/mm跨中径向位移/mm(a)跨中径向位移(b)跨中环向位移0.20.20.40.60.8图3竖向移动集中荷载作用下圆弧曲梁的动力响应Fig.3Dynamicresponsesofthecircularlycurvedbeamsunderverticalmovingconcentratedload24工程力学YANGLu.ThestudyofultimatebearingcapacityoflightaluminumalloymobilehousewithinitialimperfectionsD.Xian:ChanganUniversity,2008.(inChinese)焦隽隽.轻型铝合金活动房屋在地震作用下的动力响应D.西安:长安大学,2008.JIAOJunjun.Thedynamicresponseoflightaluminumalloy mobile building structure D.Xian:ChanganUniversity,2008.(inChinese)33毕磊.野营轻型铝合金房屋在爆炸荷载作用下的动力响应分析D.西安:长安大学,2010.BI Lei.Camping lightweight aluminum housing in theblast loading of dynamic response analysis D.Xian:ChanganUniversity,2010.(inChinese)34MEIMANDV,TORABIANS,SCHAFERBW,etal.Incremental dynamic analysis and seismic performanceevaluation of an aluminum framed building comparedwith steel C/Proceedings of the Eighth InternationalConference on Steel and Aluminium Structures.HongKong,China:The University of Hong Kong,2016,13261340.35FEMAP695,Quantificationofbuildingseismicperformance factors S.Washington,DC:FederalEmergencyManagementAgency,2009.36TANXW,ZHANGZ,HUANGJJ,etal.Second-order37nonlinear analysis of aluminium alloy frame C/Proceedings of the 2017 3rd International Forum onEnergy,Environment Science and Materials(IFEESM2017).Amsterdam:AtlantisPress,2018:124127.刘翔.集成式铝合金结构房屋标准编制及足尺试验研究D.广州:华南理工大学,2018.LIUXiang.Standardpreparationandfull-scalemodelteststudy of integrated aluminum alloy housing D.Guangzhou:South China University of Technology,2018.(inChinese)38WANGZX,WANGYQ,LIBB,etal.ExperimentalandnumericalstudyonseismicbehaviourofaluminiumalloyframesJ.JournalofBuildingEngineering,2022,50:104231.39GB504292007,铝合金结构设计规范S.北京:中国计划出版社,2008.GB 504292007,Code for design of aluminiumstructuresS.Beijing:ChinaPlanningPress,2008.(inChinese)40GB500112010,建筑抗震设计规范S.北京:中国建筑工业出版社,2010.GB500112010,CodeforseismicdesignofbuildingsS.Beijing:ChinaArchitecture&BuildingPress,2010.(inChinese)41(上接第 24 页)袁驷,叶康生,王珂.平面曲梁面内自由振动分析的自适应有限元法J.工程力学,2009,26(增刊2):126132.YUANSi,YEKangsheng,WANGKe.Aself-adaptivefem for free vibration analysis of planar curved beamswithvariablecross-sectionsJ.EngineeringMechanics,2009,26(Suppl2):126132.(inChinese)7吴玉华,楼文娟.水平随机地震激励下弹性圆拱的半解析法J.浙江大学学报(工学版),2010,44(1):156159.WUYuhua,LOUWenjuan.Semi-analyticalsolutionofelastic circular arch under horizontal random seismicexcitationJ.JournalofZhejiangUniversity(Engineering Science),2010,44(1):156 159.(inChinese)8EFTEKHARISA.Differentialquadratureprocedureforin-planevibrationanalysisofvariablethicknesscirculararches traversed by a moving point loadJ.AppliedMathematicalModelling,2016,40(7/8):46404663.9黄伟明,王金昌,徐日庆,等.基于弹性地基曲梁理论的盾构隧道管片分析方法J.浙江大学学报(工学版),2020,54(4):787795.HUANGWeiming,WANGJinchang,XURiqing,etal.StructuralanalysisofshieldtunnelliningusingtheoryofcurvedbeamrestingonelasticfoundationJ.Journalof10ZhejiangUniversity(EngineeringScience),2020,54(4):787795.(inChinese)SHEN X D,CHEN W Q,WU Y F,et al.Dynamicanalysis of partial-interaction composite beamsJ.Composites Science and Technology,2011,71(10):12861294.11DEMIRHANPA,TASKINV.LevysolutionforbendinganalysisoffunctionallygradedsandwichplatesbasedonfourvariableplatetheoryJ.CompositeStructures,2017,177:8095.12LIU X X,WANG Y,WANG G N,et al.DynamicanalysisofRCbeamsexternallybondedwithFRPplatesusing state space methodJ.Engineering Structures,2022,253:113788.13WU J S,CHIANG L K.A new approach for freevibration analysis of arches with effects of sheardeformationandrotaryinertiaconsideredJ.JournalofSoundandVibration,2004,277(1/2):4971.14BAOGJ,XURQ.Dynamicstiffnessmatrixofpartial-interactioncompositebeamsJ.AdvancesinMechanicalEngineering,2015,7(3):17.15EISENBERGERM,EFRAIME.In-planevibrationsofsheardeformablecurvedbeamsJ.InternationalJournalfor Numerical Methods in Engineering,2001,52(11):12211234.16工程力学55