基于状态空间法的圆弧曲梁面内动力学分析.pdf
《基于状态空间法的圆弧曲梁面内动力学分析.pdf》由会员分享,可在线阅读,更多相关《基于状态空间法的圆弧曲梁面内动力学分析.pdf(7页珍藏版)》请在咨信网上搜索。
1、文章编号:1000-4750(2023)Suppl-0019-06基于状态空间法的圆弧曲梁面内动力学分析刘兴喜1,杨博1,徐荣桥2(1.浙江理工大学土木工程系,浙江,杭州310018;2.浙江大学土木工程系,浙江,杭州310058)摘要:该文采用状态空间法研究了圆弧曲梁面内动力学特性。通过选择合适的状态变量,建立了相应的状态空间列式,并给出了固有频率和振动模态的求解过程。进一步通过引入辛内积的概念,建立了三种工程常见边界条件(简支、固支和自由)下圆弧曲梁面内振动模态关于质量和旋转惯量的正交关系式。在此基础上,运用模态叠加法给出了非齐次状态方程的解析解,并得到了圆弧曲梁在竖向移动集中荷载作用下的
2、瞬态响应。数值算例结果表明该文方法是十分精确和可靠的。关键词:圆弧曲梁;自由振动;强迫振动;状态空间法;辛内积;模态正交性中图分类号:TB122文献标志码:Adoi:10.6052/j.issn.1000-4750.2022.06.S018IN-PLANEDYNAMICANALYSISOFCIRCULARLYCURVEDBEAMSBASEDONTHESTATESPACEMETHODLIUXing-xi1,YANGBo1,XURong-qiao2(1.DepartmentofCivilEngineering,ZhejiangSci-TechUniversity,Hangzhou,Zhejiang
3、310018,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 andcorrespond
4、ingmodalshapesofcircularlycurvedbeamsarethenobtainedonthebasisofstate-spaceformulae.Byutilizingtheconceptionofsymplecticinnerproduct,themodeorthogonalitywithrespecttothemassandrotaryinertiapropertiesforcurvedbeamsunderthreecommonboundaryconditionsinengineering(simply-supported,clampedandfree)isestab
5、lished.Basedontheorthogonalityrelation,themodesuperpositionmethodisutilizedtoobtainthesolutionoftheinhomogeneousequationforforcedvibrationanddynamicresponseofcurvedbeamssubjectedtoaverticalmovingconstantforce.Thenumericalresultsshowthattheproposedmethodisveryaccurateandreliable.Keywords:circularly c
6、urved beams;free vibration;forced vibration;state space method;sympletic innerproduct;modalorthogonality由于曲梁结构具有造型优美、施工简便及受力性能好等优点,已被广泛应用于土木工程、航空航天、工业机械等领域1。在实际工程中,曲梁不可避免地会受到动力荷载的作用(如地震作用、风荷载、移动荷载等),动力学分析是这类结构计算分析的重要内容。与直梁相比,由于受到曲率的影响,曲梁的动力特性很复杂,目前缺乏对其进行较为精确分析的方法。因此,建立一套精确、可靠的计算方法,高效地获得曲梁动力学问题的精确解,具
7、有十分重要的意义。DENHARTOG2最早使用 Rayleigh-Ritz 法得到 了 圆 弧 拱 面 内 振 动 的 频 率。BICKFORD 和收稿日期:2022-06-27;修改日期:2022-11-25基金项目:国家自然科学基金项目(12272335)通讯作者:徐荣桥(1972),男,浙江人,教授,博士,从事组合桥梁结构、FRP加固分析理论等方面的研究(E-mail:).作者简介:刘兴喜(1989),男,湖北人,讲师,博士,从事 FRP 加固分析理论、结构动力分析理论等方面的研究(E-mail:);杨博(1979),男,宁夏人,教授,博士,主要从事复合材料力学与结构方面的研究(E-ma
8、il:).第40卷增刊Vol.40Suppl工程力学2023 年 6月June2023ENGINEERINGMECHANICS19STROM3运用传递矩阵法对曲梁面内和面外自由振动进行了研究。基于精确的动力刚度法,HOWSON和 JEMAH4、TSENG 等5分析了曲梁的面内自振特性。叶康生和殷振炜1、YANG 等6和袁驷等7采用有限元这一应用范围广泛的数值方法研究了曲梁的面内自由振动。虽然国内外学者对曲梁的自由振动特性做了大量的研究,关于其强迫振动的研究则相对不足。吴玉华和楼文娟8采用虚拟激励法与 Galerkin 法对水平随机地震作用下弹性圆拱的随机响应进行了研究。EFTEKHARI9运用
9、基于微分求积法得到了移动集中荷载作用下曲梁的瞬态响应。但文献 8 和文献 9 都是基于 Euler曲梁理论,没有考虑剪切变形和转动惯量的影响。状态空间法通过选择能量对偶的两类物理量如位移和内力作为状态变量,能将描述结构力学特性的控制方程整理成简洁的矩阵形式10,即状态方程。这种方法具有计算效率高、易于编程等优点,已被广泛地应用于结构的静力和动力分析中1113。然而,在曲梁的强迫振动分析中尚鲜见报道。本文基于 Timoshenko 曲梁理论,通过选择合适的状态变量,建立了圆弧曲梁面内动力学问题的状态空间列式,并给出了固有频率和振动模态的求解过程。采用 SHEN 等11和 LIU 等13的方法,即
10、通过引入辛内积,建立了普通边界条件下圆弧曲梁振动模态关于质量和旋转惯量的正交关系式。基于所建立的正交关系式,运用模态叠加法导出了非齐次状态方程即强迫振动方程的解析解,并得到了圆弧曲梁在竖向移动荷载作用下的瞬态响应。最后,给出了数值算例,结果表明采用本文方法进行圆弧曲梁面内动力学分析是十分精确和可靠的。1状态方程列式lRsszuwqtqnqm考虑任意动力荷载作用下的圆弧曲梁,如图 1所示,其中曲线长度和半径为分别为 和。沿曲梁中曲线建立坐标轴 轴,将 轴逆时针旋转 90得到法向坐标轴 轴。、分别是轴线上的环向位移和径向位移,为截面的转角,、和分别为作用在曲梁轴线上的环向分布荷载、径向分布荷载以及
11、分布弯矩。曲梁的微元及其横截面示意图如图 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+dN
12、qn(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式中:为状态向量,上标 为矩阵的转置;系统矩阵、惯性力列
13、向量、荷载列向量的表达式分别为: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
14、。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.
15、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)为两端简支曲梁的频率方程,它是一个关于频
16、率 的超越方程,可以采用二分法来搜索它的根15。一旦求得某一阶频率后,将其代回式(13),即可求得在处的状态向量(包含一个任意常数因子),再通过式(10)可确定曲梁内任意位置处的状态向量,这样就可以确定该频率对应的振动模态(也包含一个任意常数因子)。对于其他边界条件,同样能求得相应的频率和振动模态。2.3模态关于质量和旋转惯量的正交关系式vivj为了运用模态叠加法计算曲梁在动力荷载作用下的响应,必须建立振动模态关于质量和旋转惯量的正交关系式。本节将采用 SHEN 等11和 LIU等13的方法,即基于状态空间列式和辛内积的概念来进行推导。首先考察对应于两个不同模态的状态向量 和,定义如下辛内积:
17、(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)两边同时左乘和,并
18、在曲梁的整个弧长上进行积分,然后将两式相减,可以得到: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)
19、考虑实际工程中常见的 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(
20、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
21、)代入式(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、)关于的微分方程的通解可用杜哈梅积分表示为: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+
23、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
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 状态 空间 圆弧 曲梁面内 动力学 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。