CN104778377B - 一种组合梁弯曲振动的固有频率分析方法 - Google Patents

一种组合梁弯曲振动的固有频率分析方法 Download PDF

Info

Publication number
CN104778377B
CN104778377B CN201510222482.XA CN201510222482A CN104778377B CN 104778377 B CN104778377 B CN 104778377B CN 201510222482 A CN201510222482 A CN 201510222482A CN 104778377 B CN104778377 B CN 104778377B
Authority
CN
China
Prior art keywords
mrow
msub
mmultiscripts
msubsup
msup
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510222482.XA
Other languages
English (en)
Other versions
CN104778377A (zh
Inventor
江海深
董良
赵跃民
刘初升
彭利平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN201510222482.XA priority Critical patent/CN104778377B/zh
Publication of CN104778377A publication Critical patent/CN104778377A/zh
Application granted granted Critical
Publication of CN104778377B publication Critical patent/CN104778377B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)

Abstract

一种组合梁弯曲振动的固有频率分析方法，属于横梁振动固有频率分析方法。根据组合梁几何特征，构建一种基于经典Euler‑Bernoulli梁(EB)与Timoshenko梁(TB)理论的混合梁单元ETE‑B，各ETE‑B中EB梁单元和TB梁单元通过边界条件连续进行联接，建立组合梁的动力学模型、弯曲振动方程，获得经典边界条件下组合梁的参数化频率方程，最后利用一维搜索法确定组合梁的固有频率；所述的一维搜索法在于先确定固有圆频率的可行域，进而在各个可行域内使用二分法获得固有频率数值。优点：1)对组合梁动力学特性的分析涉及的物理意义明确；2)固有频率的获取只与横梁材料、尺寸有关，建立的参数化频率方程具有普适性，无需模型试验或三模建模下的有限元分析。

Description

一种组合梁弯曲振动的固有频率分析方法
技术领域
本发明涉及一种横梁振动固有频率分析方法，特别涉及一种组合梁弯曲振动的固有频率分析方法。
背景技术
当横梁结构跨度很大时，自重作用下，横梁中部挠度最大，一般通过内焊筋板的方式来加大横梁刚度，减少变形。由于筋板与梁体结构材料相似，这种内焊筋板横梁可以看作是一种组合梁。组合梁一般用于重型设备的大跨度承重结构(如起重机桥臂)、大型建筑物的承载结构(如房屋钢结构)和选矿设备承重结构(如大型振动筛承重梁)等。定量分析组合梁弯曲振动的固有频率，可以避免实际中横梁承受近共振区载荷而引起破坏，对指导结构动力学设计、保障设备可靠性和安全生产具有重要意义。
在结构动力学领域，目前对横梁结构固有频率的分析方法主要包括三个途径：(1)基于结构动力学基础理论的理论计算；(2)基于结构试验模型的试验模态分析；(3)基于有限元原理的商业软件分析。对于组合梁，若采用途径(1)，先建立梁体结构的动力学模型，沿用传统结构动力学理论分析其弯曲振动特性时，需同时考虑一维梁结构和二维板结构的耦合作用，微分方程分析过程比较复杂，数值求解困难；若采用途径(2)，需先制作满足相似原理的梁体结构的结构试验模型，确定合适的梁体结构实际边界条件实现方法，同时对试验模态分析所需的传感器、采集卡和分析软件的投入成本高，试验过程操作复杂；若采用途径(3)，需先根据梁体结构几何参数建立其三维模型，导入有限元商业软件中，设定基本的结构材料特性和边界条件后进行分析，一旦结构变化，则重复上述过程，因而比较繁琐。
中国专利CN201410145273.5公开了一种应用改进微分变换法计算欧拉-伯努利梁固有频率的方法，应用改进的微分变换法求解均质欧拉-伯努利梁的自由振动问题，通过迭代以收敛级数的形式得到非线性问题的近似解，得到了四阶固有频率与模态振型等闭式解；中国专利CN201210378363.X公开了一种适于梁弯曲振动分析的传递矩阵计算方法，直接从梁的挠曲线函数出发建立传递矩阵，然后将边界条件代入进行计算，无需改变传递矩阵，也无需考虑梁的变形形式。
发明内容
本发明的目的是要提供一种组合梁弯曲振动的固有频率分析方法，该方法简单，物理意义明确，计算精度高，无需重复性的模型试验或三维建模下的有限元分析。
本发明目的是通过以下技术方案实现的：根据组合梁几何特征，构建一种基于经典Euler-Bernoulli梁(EB)与Timoshenko梁(TB)理论的混合梁单元ETE-B，各ETE-B中EB梁单元和TB梁单元通过边界条件连续进行联接，建立组合梁的动力学模型、弯曲振动方程，获得经典边界条件下组合梁的参数化频率方程，最后利用一维搜索法确定组合梁的固有频率；具体步骤如下：
步骤1、组合梁的动力学模型建立：根据经典梁结构动力学理论，梁截面尺寸相比于跨度而言很小，一般都可以认为是Euler-Bernoulli梁(EB)；然而，当梁体或梁单元的深跨比比1/5大很多时，通常被看作Timoshenko梁(TB)，研究其弯曲振动特征时，需考虑转动惯量和截面剪切变形的影响；混合梁单元包括一段代号为TB-i+1的Timoshenko梁和两段代号分别为EB-i和EB-i+2的Euler-Bernoulli梁，以Timoshenko梁的代号数值的一半标示，即两类梁单元在交界面Fi和Fi+1处通过双边V缝周边对焊联接，材料体密度均为ρ，弹性模量均为E。TB-i+1中，梁段宽b、高h、长li+1，满足h/li+1∈[1/5,+∞]，截面面积Ai+1，相对于x轴的二次惯性矩为Ii+1，截面的剪切修正系数为κ，材料的剪切弹性模量为G。EB-i和EB-i+2中，梁段宽均为b、高均为h、长分别为li和li+2；所述的li+2为li＝li+2，壁厚均为h0，满足h/li和h/li+2∈[0,1/5)，截面面积分别为Ai和Ai+2(Ai+1＝Ai+2)，相对于x轴的二次惯性矩分别为Ii和Ii+2(Ii+1＝Ii+2)；所述的Ii+2为Ii+1＝Ii+2.；
步骤2、混合梁单元ETE-B的动力学分析：在Euler-Bernoulli梁单元EB-i的固有圆频率ω满足：
式(2)中λ为EB频率参数，且有振型函数：
Yi(x)＝ei,1sinλx+ei,2cosλx+ei,3sinhλx+ei,4coshλx. (3)
式(3)中ei,1,ei,2,ei,3和ei,4为EB-i梁单元边界条件决定的系数。
对于其他等截面等长Euler-Bernoulli梁单元，只需对应更改梁单元序号i；
由于是同一根组合梁的不同单元，Timoshenko梁单元TB-i+1的固有圆频率与步骤2中EB-i的固有圆频率ω相同，在满足
的前提下，有一组虚根：
式(7)中为TB虚频率参数，i为虚数单位，满足
而还有一组实根
式(9)中为TB实频率参数，满足
TB-i+1的弯曲变形振型表达式：
式中，ti+1,1,ti+1,2,ti+1,3和ti+1,4为TB-i+1梁单元边界条件决定的系数。
TB-i+1的截面转角振型表达式：
对于其他等截面等长Timoshenko梁单元，只需对应更改梁单元序号i+1。
步骤3、混合梁单元ETE-B联接条件的数学描述：定义EB系数和TB系数则在交界面Fi＝(EB-i)∩(TB-i+1)处，其中∩符号表示两个梁单元联接位置，存在振型系数界面联接向量：
和振型界面联接特征矩阵
类似地，在交界面Fi+1＝(TB-i+1)∩(EB-i+2)处存在振型系数界面联接向量：
和振型界面联接特征矩阵
因此，任意单段混合梁单元的振型联接特征矩阵和振型系数联接向量分别为：
和
且满足：
左上标为ETE-B的段数，左下标为首段组合梁中TB梁序号，
对于任意连续两段ETE-B，TB即梁单元序号为i+1和i+3，共同拥有EB-i+2，振型联接特征矩阵和振型系数联接向量分别为：
和
且满足：
其中，和分别由和将其各含右下标元素的右下标数值上加2获得，未给出大小的“0”的维数视矩阵整体而定。
对于n段ETE-B构成的组合梁，TB梁单元序号依次为2,4,6,…,2j,…,2n，振型联接特征矩阵和振型系数联接向量分别为：
和
且满足：
其中，左上标“n”为组合梁中ETE-B的段数，左下标“2”为首段混合梁单元中TB梁单元序号。
步骤4、组合梁边界条件的数学描述：对于n段ETE-B构成的组合梁，其经典边界条件下的数学表达如下：
左端固支(CL)
左端简支(PL)
右端自由(FR)
右端简支(PR)
右端固支(CR)
步骤5、组合梁参数化频率方程的建立：两端经典边界条件分别是式(41)、(42)描述的和式(43)～(45)描述的的n段ETE-B构成的组合梁的振型特征矩阵D和振型系数向量C分别为：
和
且满足：
[D](8n+4)×(8n+4)·[C](8n+4)×1＝0. (48)
令组合梁的特征方程系数矩阵D的行列式为零，行列式以特征函数f(λ)＝|D|表示，求解参数化频率方程：
f(λ)＝|D|＝0 (49)
即获得整根组合梁的固有圆频率ω，按式f＝ω/2π计算组合梁的各阶固有频率。
所述的混合梁单元ETE-B，将任意两个EB及其间的一个TB看作一个混合梁单元ETE-B，则组合梁结构可认为是多个混合梁单元ETE-B构成的；所有EB等截面、等长、空心，而所有TB等截面、等长、实心；组合梁的首尾两段梁单元都是EB梁。
所述的一维搜索法：以式(2)中EB频率参数λ为自变量，给定其取值范围、递增步长Δλ和初始值λ0，在λ数值递增形成的各个子区间中确定f(λ)变号的子区间[λk,λk+Δλ](即满足f(λk)·f(λk+Δλ)＜0)，作为固有圆频率ω的可行域，在各个可行域内使用二分法，就能得到固有圆频率值，按式f＝ω/2π计算对应固有频率。
有益效果及优点：
1)根据区分经典的Euler-Bernoulli梁(EB)和Timoshenko梁(TB)类型的几何尺寸标准，提出了边界连续条件的EB和TB组成的混合梁单元，用来构造工程中常见的组合梁结构，物理意义明确；
2)固有频率的获取只与横梁材料、尺寸有关，建立的参数化频率方程具有普适性，无需重复性的模型试验或三维建模下的有限元分析。
附图说明
图1是本发明方法的技术流程图。
图2是本发明的实施组合梁固有频率分析的任意连续两段ETE-B的动力学模型。
图3(a)是与本发明的实施组合梁等长的Euler-Bernoulli梁E动力学模型。
图3(b)是本发明的实施固有频率分析的单段ETE-B单元的组合梁ETE-1。
图3(c)是本发明的实施固有频率分析的两段ETE-B单元的组合梁ETE-2。
具体实施方式
下面结合附图对本发明实施例作进一步的说明：
根据组合梁几何特征，构建一种基于经典Euler-Bernoulli梁(EB)与Timoshenko梁(TB)理论的混合梁单元ETE-B，各ETE-B中EB梁单元和TB梁单元通过边界条件连续进行联接，建立组合梁的动力学模型、弯曲振动方程，获得经典边界条件下组合梁的参数化频率方程，最后利用一维搜索法确定组合梁的固有频率；具体步骤如下：
步骤1、组合梁的动力学模型建立：根据经典梁结构动力学理论，梁截面尺寸相比于跨度而言很小，一般都可以认为是Euler-Bernoulli梁(EB)；然而，当梁体或梁单元的深跨比比1/5大很多时，通常被看作Timoshenko梁(TB)，研究其弯曲振动特征时，需考虑转动惯量和截面剪切变形的影响；混合梁单元包括一段代号为TB-i+1的Timoshenko梁和两段代号分别为EB-i和EB-i+2的Euler-Bernoulli梁，以Timoshenko梁的代号数值的一半标示，即两类梁单元在交界面Fi和Fi+1处通过双边V缝周边对焊联接，材料体密度均为ρ，弹性模量均为E；TB-i+1中，梁段宽b、高h、长li+1，满足h/li+1∈[1/5,+∞]，截面面积Ai+1，相对于x轴的二次惯性矩为Ii+1，截面的剪切修正系数为κ，材料的剪切弹性模量为G；EB-i和EB-i+2中，梁段宽均为b、高均为h、长分别为li和li+2(li＝li+2)，所述的li+2为li＝li+2；壁厚均为h0，满足h/li和h/li+2∈[0,1/5)，截面面积分别为Ai和Ai+2(Ai+1＝Ai+2)，相对于x轴的二次惯性矩分别为Ii和Ii+2(Ii+1＝Ii+2)；所述的Ii+2为Ii+1＝Ii+2；
步骤2、混合梁单元ETE-B的动力学分析：对于Euler-Bernoulli梁单元EB-i，按经典梁理论，其弯曲变形yi随时间t和位置坐标x的微分方程如下：
取
其中，λ为EB频率参数，ω为固有圆频率，且有振型函数：
Yi(x)＝ei,1sinλx+ei,2cosλx+ei,3sinhλx+ei,4coshλx (3)
式中，ei,1,ei,2,ei,3和ei,4为EB-i梁单元边界条件决定的系数；
对于其他等截面等长Euler-Bernoulli梁单元，可类似进行上述分析，同时对应更改Euler-Bernoulli梁单元序号i即可；
对于Timoshenko梁单元TB-i+1，设弯曲变形和截面转角分别为yi+1(x,t)和随时间t和位置坐标x的微分方程如下：
设弯曲变形和截面转角振型函数分别为Yi+1(x)和φi+1(x)，固有圆频率为ω，同一根组合梁时，固有圆频率应与式(2)中相同，且满足：
κGAi+1EIi+1·r4+(κGAi+1ρIi+1+ρAi+1EIi+1)ω2r2+ρ2ω4Ai+1Ii+1-κGAi 2 +1ρω2＝0 (5)
式(5)中，r为TB频率参数。
上述方程在
的前提下，有一组虚根：
其中为TB虚频率参数，i为虚数单位，满足
而还有一组实根
其中为TB实频率参数，满足
根据常微分方程解的特性，此时有TB-i+1的弯曲变形振型表达式：
式中，ti+1,1,ti+1,2,ti+1,3和ti+1,4为TB-i+1梁单元边界条件决定的系数；
将式(11)代入式(4)，得到TB-i+1的截面转角振型表达式：
对于其他等截面等长Timoshenko梁单元，可类似进行上述分析，同时对应更改Timoshenko梁单元序号i+1即可；
步骤3、混合梁单元ETE-B联接条件的数学描述：
EB和TB梁单元分别依次建立其局部坐标系内振动微分方程，组合成混合梁单元不考虑焊接对材料属性的影响，焊接联接后，EB和TB梁单元在结构上属于一体，则EB-i在x＝li和TB-i+1对应x＝0处的几何边界条件即位移与转角和力边界条件即交界面剪力V和弯矩M应相同，归结为两类梁单元的联接条件，即有：
Yi(li)＝Yi+1(0) (13)
Yi′(li)＝φi+1(0) (14)
EIiYi″(li)＝EIi+1φ′i+1(0) (15)
EIiYi″′(li)＝κGAi+1[Y′i+1(0)-φi+1(0)] (16)
TB-i+1梁单元在x＝li+1和EB-i+2梁单元对应x＝0处的边界条件相同，有：
Yi+1(li+1)＝Yi+2(0) (17)
φi+1(li+1)＝Y′i+2(0) (18)
EIi+1φ′i+1(li+1)＝EIi+2Y″i+2(0) (19)
κGAi+1[Y′i+1(li+1)-φi+1(li+1)]＝EIi+2Y″′+2(0) (20)
将EB-i振型函数式(3)、TB-i+1弯曲变形振型函数式(11)和截面转角振型函数式(12)及对应的EB-i+2两类振型函数：
Yi+2(x)＝ei+2,1sinλx+ei+2,2cosλx+ei+2,3sinhλx+ei+2,4coshλx (21)
代入到式(13)～(20)中，定义和则在交界面Fi＝(EB-i)∩(TB-i+1)处，∩符号表示两个梁单元联接位置，存在振型系数界面联接向量：
和振型界面联接特征矩阵：
类似地，在交界面Fi+1＝(TB-i+1)∩(EB-i+2)处存在振型系数界面联接向量：
和振型界面联接特征矩阵：
因此，任意单段混合梁单元的振型联接特征矩阵和振型系数联接向量分别为：
和
且满足：
左上标为组合梁中ETE-B的段数，左下标为首段组合梁中TB梁序号；
类似地推导，对于任意连续两段ETE-B，TB梁单元序号为i+1和i+3，共同拥有EB-i+2，振型联接特征矩阵和振型系数联接向量分别为：
和
且满足：
其中，和分别由和将其各元素的右下标+2获得，未给出大小的“0”的维数视矩阵整体而定；
依此类推，对于n段ETE-B构成的组合梁，TB梁单元序号依次为2,4,6,…,2j,…,2n，振型联接特征矩阵和振型系数联接向量分别为：
和
且满足：
其中，左上标“n”为组合梁中ETE-B的段数，左下标“2”为首段混合梁单元中TB梁单元序号；
步骤4、组合梁边界条件的数学描述：
对于n段ETE-B构成的组合梁，其经典边界条件下的数学表达如下；若仅考虑左端
固支或简支情况，则EB-1局部坐标系x＝0处的边界条件为：
左端固支(CL)Y1(0)＝0,Y1′(0)＝0 (36)
左端固支(PL)Y1(0)＝0,Y1″(0)＝0 (37)仅考虑右端固支、简支和自由的情况，EB-2n+1对应局部坐标系中x＝l2n+1处的边界条件为：
右端自由(FR)Y″2n+1(l2n+1)＝0,Y″′2n+1(l2n+1)＝0 (38)
右端简支(PR)Y2n+1(l2n+1)＝0,Y″2n+1(l2n+1)＝0 (39)
右端固支(CR)Y2n+1(l2n+1)＝0,Y′2n+1(l2n+1)＝0 (40)
提取边界条件式(36)～(40)中含λ的三角函数，有以下表达式：
左端固支(CL)
左端简支(PL)
右端自由(FR)
右端简支(PR)
右端固支(CR)
步骤5、组合梁的参数化频率方程的建立：
结合振型联接特征矩阵和振型系数联接向量两端经典边界条件(分别是式(41)、(42)描述的和式(43)～(45)描述的)的n段ETE-B构成的组合梁的振型特征矩阵D和振型系数向量C分别为：
和
且满足：
[D](8n+4)×(8n+4)·[C](8n+4)×1＝0 (48)
令组合梁的特征方程系数矩阵D的行列式为零，其中：行列式以特征函数f(λ)＝|D|表示，求解参数化频率方程：
f(λ)＝|D|＝0 (49)
即获得整根组合梁的固有圆频率ω，按式f＝ω/2π计算组合梁的各阶固有频率；
根据上面的理论推导，利用Matlab编程获取组合梁的各阶固有频率的求解；确定横梁几何参数后，按式(41)～(45)选择需考虑的两端边界条件，并按式(49)建立组合梁的参数化频率方程；考虑到该方程的非线性，采用一维搜索法进行方程根的搜索；由于各个解离散分布，为减少计算工作量，可以先找到解的邻域，再定位精确解，确定求解的基本方法是：以式(2)中EB频率参数λ为自变量，给定其取值范围、递增步长Δλ和初始值λ0，在λ数值递增形成的各个子区间中确定f(λ)变号的子区间[λk,λk+Δλ]，即满足f(λk)·f(λk+Δλ)＜0，其中λk＝λ0+k·Δλ，即为固有圆频率的可行域，这些可行域是离散、等长、短距的，依次在各个可行域内使用二分法，就能得到固有圆频率值，按式f＝ω/2π计算组合梁的各阶固有频率。
所述的混合梁单元ETE-B，将任意两个EB及其间的一个TB看作一个混合梁单元ETE-B，则组合梁结构可认为是多个混合梁单元ETE-B构成的；所有EB等截面、等长、空心，而所有TB等截面、等长、实心；组合梁的首尾两段梁单元都是EB梁。
所述的一维搜索法：以式(2)中EB频率参数λ为自变量，给定其取值范围、递增步长Δλ和初始值λ0，在λ数值递增形成的各个子区间中确定f(λ)变号的子区间[λk,λk+Δλ](即满足f(λk)·f(λk+Δλ)＜0，其中λk＝λ0+k·Δλ)，作为固有圆频率的可行域，在各个可行域内使用二分法，就能得到固有圆频率值，按式f＝ω/2π计算对应固有频率。
具体的：
图2为任意连续两段ETE-B的动力学模型。针对图3(b)、(c)所示的单段(n＝1)、两段(n＝2)混合梁单元的组合梁ETE-1、ETE-2进行数值仿真，讨论相应组合梁的弯曲振动固有频率，同时和等截面等长空心Euler-Bernoulli横梁E(图3a)进行比较，E梁相当于n＝0的组合梁情况。具体的各个梁单元基本几何参数如表1所示。TB梁单元的深跨比为1.5，满足hT/lTi∈[1/5,+∞]，截面面积均AT＝bT×hT，相对于x轴的二次惯性矩矩形截面剪切修正系数κ＝10(v+1)/(11v+12)，泊松比v＝0.28，剪切弹性模量G＝0.5E/(1+v)；EB梁单元最大深跨比为0.056，满足hE/lEi∈[0,1/5)，截面面积AE＝bE×hE-(bE-2h0)×(hE-2h0)，相对于x轴的二次惯性矩分别为梁单元均为不锈钢，密度ρ＝7750kg/m3，弹性模量E＝193GPa；
表1两种组合梁(ETE-1、ETE-2)及等长Euler-Bernoulli梁(E)的几何参数
具体步骤如下：
图3(a)中，对于E，根据表2中经典边界条件下Euler-Bernoulli梁固有频率计算通式，可以得到E梁的三种边界条件(CL-FR、PL-PR和CL-CR)下的固有频率。
表2经典边界条件下Euler-Bernoulli梁固有频率计算通式中λiL的数值
图3(b)中ETE-1，相当于图2中n＝1段的ETE-B构成，EB序号分别为1和3，长度分别l1和l3，且满足l1＝l3，截面尺寸为b×h，侧壁厚度均为h0。TB序号为2，长度为l2，截面尺寸为b×h。将n＝1代入式(33)中，引入式(27)，则ETE-1的振型联接特征矩阵为：
其中，和分别为i＝1下的式(24)和式(26)的数学表达。对应的振型系数联接向量为：
其中，和分别为i＝1下的式(23)和式(25)的数学表达。两端经典边界条件的1段ETE-B构成的组合梁的振型特征矩阵[D]12×12和振型系数向量[C]12×1分别为：
和
则求解参数化频率方程：
f1(λ)＝|[D]12×12|＝0 (54)
即获得ETE-1固有频率。
图3(c)中ETE-2，相当于图2中n＝2段的ETE-B构成，EB序号分别为1、3和5，长度分别l1、l3和l5，且满足l1＝l3＝l5，截面尺寸均为b×h，侧壁厚度均为h0。TB序号为2和4，长度分别为l2和l4，且满足l2＝l4，截面尺寸均为b×h。将n＝2代入式(33)中，引入式(27)，则ETE-2的振型联接特征矩阵为：
其中，未给出的“0”的维数视矩阵整体而定，和分别为i＝1下的式(24)和式(26)的数学表达，和分别为i＝3下的式(24)和式(26)的数学表达。对应的振型系数联接向量为：
其中，和分别为i＝1下的式(23)和式(25)的数学表达，和分别为i＝3下的式(23)和式(25)的数学表达。两端经典边界条件的2段ETE-B构成的组合梁的振型特征矩阵[D]20×20和振型系数向量[C]20×1分别为：
和
则求解参数化频率方程：
f2(λ)＝|[D]20×20|＝0 (59)
即获得ETE-2固有频率。
利用一维搜索法求解参数化频率方程(54)、(59)，以λ为自变量，设定其取值范围[0.1,10]、递增步长Δλ＝0.01和初始值λ0＝0.1，在λ数值递增形成的各个子区间中确定f(λ)变号的子区间[λk,λk+Δλ]，即满足f(λk)·f(λk+Δλ)＜0，其中λk＝λ0+k·Δλ，作为固有圆频率的可行域，依次在各个可行域内使用二分法，分别获取本发明实施例的组合梁ETE-1、ETE-2的前四阶固有频率和等长Euler-Bernoulli梁(E)前四阶固有频率，并通过与有限元软件分析的结果比较进行验证，如表3～5所示。可以发现：各表中利用经典理论和本发明方法获得的E梁前四阶固有频率几乎完全相同，表明本发明原理和模型推导过程的准确性；就本发明方法获得的固有频率值而言，边界条件CL-FR、PL-PR和CL-CR的变化趋势，横梁端部约束自由度增多，故各阶固有频率整体增大，而模态阶数越高，固有频率也增大，与经典的横梁振动理论的性质一致。CL-FR、PL-PR和CL-CR边界条件下本发明方法获得的固有频率值与软件分析结果相对差异绝对值分别在5％、5％和8％以内，主要来源软件分析设置以及TB单元质量因素，实施例中TB单元相较于组合梁整体在轴向长度的占比较大的影响。实际中，一方面组合梁跨度大(如起重机桥臂达20米，振动筛承重梁跨度达4米以上)，另一方面，TB单元类似板结构，长度小，因而相较于组合梁整体在轴向长度的占比很小，TB单元的集中质量因素影响可以忽略，因此，可以认为本方法的分析结果在工程误差允许范围内。
表3 CL-FR边界条件下ETE-1、ETE-2和E梁的前四阶固有频率(Hz)
表4 PL-PR边界条件下ETE-1、ETE-2和E梁的前四阶固有频率(Hz)
表5 CL-CR边界条件下ETE-1、ETE-2和E梁的前四阶固有频率
该实施例说明了本发明一种组合梁弯曲振动的固有频率分析方法的有效性。

Claims (3)

1.一种组合梁弯曲振动的固有频率分析方法，其特征在于：根据组合梁几何特征，构建一种基于经典Euler-Bernoulli梁(EB)与Timoshenko梁(TB)理论的混合梁单元ETE-B，各ETE-B中EB梁单元和TB梁单元通过边界条件连续进行联接，建立组合梁的动力学模型、弯曲振动方程，获得经典边界条件下组合梁的参数化频率方程，最后利用一维搜索法确定组合梁的固有频率；具体步骤如下：
步骤1、组合梁的动力学模型建立：根据经典梁结构动力学理论，梁截面尺寸相比于跨度而言很小，一般都可以认为是Euler-Bernoulli梁(EB)；然而，当梁体或梁单元的深跨比比1/5大很多时，通常被看作Timoshenko梁(TB)，研究其弯曲振动特征时，需考虑转动惯量和截面剪切变形的影响；混合梁单元包括一段代号为TB-i+1的Timoshenko梁和两段代号分别为EB-i和EB-i+2的Euler-Bernoulli梁，以Timoshenko梁的序号数值(i+1)的一半标示，即两类梁单元在交界面Fi和Fi+1处通过双边V缝周边对焊联接，材料体密度均为ρ，弹性模量均为E；TB-i+1中，梁段宽b、高h、长li+1，满足h/li+1∈[1/5,+∞]，截面面积Ai+1，相对于x轴的二次惯性矩为Ii+1，截面的剪切修正系数为κ，材料的剪切弹性模量为G；EB-i和EB-i+2中，梁段宽均为b、高均为h、长分别为li和li+2；所述的li+2为li＝li+2，壁厚均为h0，满足h/li和h/li+2∈[0,1/5]，截面面积分别为Ai和Ai+2，所述的Ai+2为Ai+1＝Ai+2，相对于x轴的二次惯性矩分别为Ii和Ii+2；所述的Ii+2为Ii+1＝Ii+2；
步骤2、混合梁单元ETE-B的动力学分析：在Euler-Bernoulli梁单元EB-i的固有圆频率ω满足：
<mrow> <msup> <mi>&amp;lambda;</mi> <mn>4</mn> </msup> <mo>=</mo> <mfrac> <mrow> <msub> <mi>&amp;rho;A</mi> <mi>i</mi> </msub> </mrow> <mrow> <msub> <mi>EI</mi> <mi>i</mi> </msub> </mrow> </mfrac> <msup> <mi>&amp;omega;</mi> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
式(2)中λ为EB频率参数，且有振型函数：
Yi(x)＝ei,1sinλx+ei,2cosλx+ei,3sinhλx+ei,4coshλx. (3)
式(3)中ei,1,ei,2,ei,3和ei,4为EB-i梁单元边界条件决定的系数；
对于其他等截面等长Euler-Bernoulli梁单元，只需对应更改梁单元序号i；
由于是同一根组合梁的不同单元，Timoshenko梁单元TB-i+1的固有圆频率与步骤2中EB-i的固有圆频率ω相同，在满足
<mrow> <mi>&amp;omega;</mi> <mo>&lt;</mo> <msqrt> <mfrac> <mrow> <mn>2</mn> <msub> <mi>A</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <msup> <mi>E&amp;kappa;</mi> <mn>2</mn> </msup> <msup> <mi>G</mi> <mn>2</mn> </msup> </mrow> <mrow> <msub> <mi>&amp;rho;I</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>E</mi> <mo>+</mo> <mi>&amp;kappa;</mi> <mi>G</mi> <mo>)</mo> </mrow> </mrow> </mfrac> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
的前提下，有一组虚根：
<mrow> <msub> <mi>r</mi> <mrow> <mn>1</mn> <mo>,</mo> <mn>2</mn> </mrow> </msub> <mo>=</mo> <mo>&amp;PlusMinus;</mo> <mover> <mi>r</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;CenterDot;</mo> <mi>i</mi> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
式(7)中为TB虚频率参数，i为虚数单位，满足
<mrow> <mover> <mi>r</mi> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <msqrt> <mrow> <mn>2</mn> <msub> <mi>&amp;kappa;GEI</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mi>&amp;omega;</mi> <mrow> <mo>(</mo> <msub> <mi>&amp;kappa;G&amp;rho;I</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mi>&amp;omega;</mi> <mo>+</mo> <msub> <mi>E&amp;rho;I</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mi>&amp;omega;</mi> <mo>+</mo> <msqrt> <mrow> <msup> <mi>&amp;rho;</mi> <mn>2</mn> </msup> <msubsup> <mi>I</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>2</mn> </msubsup> <msup> <mi>&amp;omega;</mi> <mn>2</mn> </msup> <msup> <mrow> <mo>(</mo> <mrow> <mi>&amp;kappa;</mi> <mi>G</mi> <mo>-</mo> <mi>E</mi> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mn>4</mn> <msup> <mi>&amp;kappa;</mi> <mn>2</mn> </msup> <msup> <mi>G</mi> <mn>2</mn> </msup> <msub> <mi>&amp;rho;A</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <msub> <mi>EI</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mrow> </msqrt> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>/</mo> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>&amp;kappa;GEI</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
而还有一组实根
<mrow> <msub> <mi>r</mi> <mrow> <mn>3</mn> <mo>,</mo> <mn>4</mn> </mrow> </msub> <mo>=</mo> <mo>&amp;PlusMinus;</mo> <munder> <mi>r</mi> <mo>~</mo> </munder> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
式(9)中为TB实频率参数，满足
<mrow> <munder> <mi>r</mi> <mo>~</mo> </munder> <mo>=</mo> <msqrt> <mrow> <mo>-</mo> <mn>2</mn> <msub> <mi>&amp;kappa;GEI</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mi>&amp;omega;</mi> <mrow> <mo>(</mo> <msub> <mi>&amp;kappa;G&amp;rho;I</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mi>&amp;omega;</mi> <mo>+</mo> <msub> <mi>E&amp;rho;I</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mi>&amp;omega;</mi> <mo>-</mo> <msqrt> <mrow> <msup> <mi>&amp;rho;</mi> <mn>2</mn> </msup> <msubsup> <mi>I</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>2</mn> </msubsup> <msup> <mi>&amp;omega;</mi> <mn>2</mn> </msup> <msup> <mrow> <mo>(</mo> <mrow> <mi>&amp;kappa;</mi> <mi>G</mi> <mo>-</mo> <mi>E</mi> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mn>4</mn> <msup> <mi>&amp;kappa;</mi> <mn>4</mn> </msup> <msup> <mi>G</mi> <mn>2</mn> </msup> <msub> <mi>&amp;rho;A</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <msub> <mi>EI</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mrow> </msqrt> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>/</mo> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>&amp;kappa;GEI</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>.</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow> 1
TB-i+1的弯曲变形振型表达式：
<mrow> <msub> <mi>Y</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>t</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mn>1</mn> </mrow> </msub> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mover> <mi>r</mi> <mo>&amp;OverBar;</mo> </mover> <mi>x</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>t</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> </mrow> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mover> <mi>r</mi> <mo>&amp;OverBar;</mo> </mover> <mi>x</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>t</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mn>3</mn> </mrow> </msub> <mi>sin</mi> <mi>h</mi> <mrow> <mo>(</mo> <munder> <mi>r</mi> <mo>~</mo> </munder> <mi>x</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>t</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mn>4</mn> </mrow> </msub> <mi>cosh</mi> <mrow> <mo>(</mo> <munder> <mi>r</mi> <mo>~</mo> </munder> <mi>x</mi> <mo>)</mo> </mrow> <mo>.</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
式中，ti+1,1,ti+1,2,ti+1,3和ti+1,4为TB-i+1梁单元边界条件决定的系数；
TB-i+1的截面转角振型表达式：
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&amp;phi;</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <mrow> <mo>(</mo> <mover> <mi>r</mi> <mo>&amp;OverBar;</mo> </mover> <mo>-</mo> <mfrac> <mrow> <msup> <mi>&amp;rho;&amp;omega;</mi> <mn>2</mn> </msup> </mrow> <mrow> <mi>&amp;kappa;</mi> <mi>G</mi> <mover> <mi>r</mi> <mo>&amp;OverBar;</mo> </mover> </mrow> </mfrac> <mo>)</mo> </mrow> <msub> <mi>t</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> </mrow> </msub> <mo>&amp;CenterDot;</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mover> <mi>r</mi> <mo>&amp;OverBar;</mo> </mover> <mi>x</mi> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <mover> <mi>r</mi> <mo>&amp;OverBar;</mo> </mover> <mo>-</mo> <mfrac> <mrow> <msup> <mi>&amp;rho;&amp;omega;</mi> <mn>2</mn> </msup> </mrow> <mrow> <mi>&amp;kappa;</mi> <mi>G</mi> <mover> <mi>r</mi> <mo>&amp;OverBar;</mo> </mover> </mrow> </mfrac> <mo>)</mo> </mrow> <msub> <mi>t</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mn>1</mn> </mrow> </msub> <mo>&amp;CenterDot;</mo> <mi>cos</mi> <mrow> <mo>(</mo> <mover> <mi>r</mi> <mo>&amp;OverBar;</mo> </mover> <mi>x</mi> <mo>)</mo> </mrow> <mo>+</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>(</mo> <munder> <mi>r</mi> <mo>~</mo> </munder> <mo>+</mo> <mfrac> <mrow> <msup> <mi>&amp;rho;&amp;omega;</mi> <mn>2</mn> </msup> </mrow> <mrow> <mi>&amp;kappa;</mi> <mi>G</mi> <munder> <mi>r</mi> <mo>~</mo> </munder> </mrow> </mfrac> <mo>)</mo> <msub> <mi>t</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mn>4</mn> </mrow> </msub> <mo>&amp;CenterDot;</mo> <mi>sinh</mi> <mo>(</mo> <munder> <mi>r</mi> <mo>~</mo> </munder> <mi>x</mi> <mo>)</mo> <mo>+</mo> <mo>(</mo> <munder> <mi>r</mi> <mo>~</mo> </munder> <mo>+</mo> <mfrac> <mrow> <msup> <mi>&amp;rho;&amp;omega;</mi> <mn>2</mn> </msup> </mrow> <mrow> <mi>&amp;kappa;</mi> <mi>G</mi> <munder> <mi>r</mi> <mo>~</mo> </munder> </mrow> </mfrac> <mo>)</mo> <msub> <mi>t</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mn>3</mn> </mrow> </msub> <mo>&amp;CenterDot;</mo> <mi>cosh</mi> <mo>(</mo> <munder> <mi>r</mi> <mo>~</mo> </munder> <mi>x</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> <mo>.</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
对于其他等截面等长Timoshenko梁单元，只需对应更改梁单元序号i+1；
步骤3、混合梁单元ETE-B联接条件的数学描述：定义EB系数和TB系数则在交界面Fi＝(EB-i)∩(TB-i+1)处，存在振型系数界面联接向量：
和振型界面联接特征矩阵
类似地，在交界面Fi+1＝(TB-i+1)∩(EB-i+2)处存在振型系数界面联接向量
和振型界面联接特征矩阵
因此，任意单段混合梁单元的振型联接特征矩阵和振型系数联接向量分别为：
<mrow> <msub> <mrow> <mo>&amp;lsqb;</mo> <mmultiscripts> <mi>D</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>1</mn> </mmultiscripts> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mn>8</mn> <mo>&amp;times;</mo> <mn>12</mn> </mrow> </msub> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mrow> <mo>&amp;lsqb;</mo> <msubsup> <mi>D</mi> <mi>i</mi> <mi>F</mi> </msubsup> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mn>4</mn> <mo>&amp;times;</mo> <mn>8</mn> </mrow> </msub> </mtd> <mtd> <msub> <mrow> <mo>&amp;lsqb;</mo> <mn>0</mn> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mn>4</mn> <mo>&amp;times;</mo> <mn>4</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mrow> <mo>&amp;lsqb;</mo> <mn>0</mn> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mn>4</mn> <mo>&amp;times;</mo> <mn>4</mn> </mrow> </msub> </mtd> <mtd> <msub> <mrow> <mo>&amp;lsqb;</mo> <msubsup> <mi>D</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>F</mi> </msubsup> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mn>4</mn> <mo>&amp;times;</mo> <mn>8</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>27</mn> <mo>)</mo> </mrow> </mrow>
和
<mrow> <msub> <mrow> <mo>&amp;lsqb;</mo> <mmultiscripts> <mi>C</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>1</mn> </mmultiscripts> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mn>12</mn> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mo>&amp;lsqb;</mo> <msubsup> <mi>C</mi> <mi>i</mi> <mi>F</mi> </msubsup> <mo>&amp;cup;</mo> <msubsup> <mi>C</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>F</mi> </msubsup> <mo>&amp;rsqb;</mo> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msubsup> <mi>E</mi> <mi>i</mi> <mi>F</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>T</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>F</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>E</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>2</mn> </mrow> <mi>F</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>28</mn> <mo>)</mo> </mrow> </mrow>
且满足：
<mrow> <mmultiscripts> <mi>D</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>1</mn> </mmultiscripts> <mo>&amp;CenterDot;</mo> <mmultiscripts> <mi>C</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>1</mn> </mmultiscripts> <mo>=</mo> <mn>0.</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>29</mn> <mo>)</mo> </mrow> </mrow>
左上标为ETE-B的段数，左下标为首段组合梁中TB梁序号，
对于任意连续两段ETE-B，TB梁单元序号为i+1和i+3，共同拥有EB-i+2，振型联接特征矩阵和振型系数联接向量分别为：
和
<mrow> <msub> <mrow> <mo>&amp;lsqb;</mo> <mmultiscripts> <mi>C</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>2</mn> </mmultiscripts> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mn>20</mn> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mo>&amp;lsqb;</mo> <mmultiscripts> <mi>C</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>1</mn> </mmultiscripts> <mo>&amp;cup;</mo> <mmultiscripts> <mi>C</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>3</mn> </mrow> <mn>1</mn> </mmultiscripts> <mo>&amp;rsqb;</mo> <mo>=</mo> <mo>&amp;lsqb;</mo> <msubsup> <mi>C</mi> <mi>i</mi> <mi>F</mi> </msubsup> <mo>&amp;cup;</mo> <msubsup> <mi>C</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>F</mi> </msubsup> <mo>&amp;cup;</mo> <msubsup> <mi>C</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>2</mn> </mrow> <mi>F</mi> </msubsup> <mo>&amp;cup;</mo> <msubsup> <mi>C</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>3</mn> </mrow> <mi>F</mi> </msubsup> <mo>&amp;rsqb;</mo> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msubsup> <mi>E</mi> <mi>i</mi> <mi>F</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>T</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>F</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>E</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>2</mn> </mrow> <mi>F</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>T</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>3</mn> </mrow> <mi>F</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>E</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>4</mn> </mrow> <mi>F</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>31</mn> <mo>)</mo> </mrow> </mrow>
且满足：
<mrow> <mmultiscripts> <mi>D</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>2</mn> </mmultiscripts> <mo>&amp;CenterDot;</mo> <mmultiscripts> <mi>C</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>2</mn> </mmultiscripts> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>32</mn> <mo>)</mo> </mrow> </mrow>
其中，和分别由和将其各含右下标元素的右下标数值上加2获得，未给出
大小的“0”的维数视矩阵整体而定；
对于n段ETE-B构成的组合梁，TB梁单元序号依次为2,4,6,…,2j,…,2n，振型联接特征矩阵和振型系数联接向量分别为：
和
<mrow> <msub> <mrow> <mo>&amp;lsqb;</mo> <mmultiscripts> <mi>C</mi> <mn>2</mn> <mi>n</mi> </mmultiscripts> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mo>(</mo> <mn>8</mn> <mi>n</mi> <mo>+</mo> <mn>4</mn> <mo>)</mo> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mo>&amp;lsqb;</mo> <mmultiscripts> <mi>C</mi> <mn>2</mn> <mn>1</mn> </mmultiscripts> <mo>&amp;cup;</mo> <mmultiscripts> <mi>C</mi> <mn>4</mn> <mn>1</mn> </mmultiscripts> <mo>&amp;cup;</mo> <mn>...</mn> <mo>&amp;cup;</mo> <mmultiscripts> <mi>C</mi> <mrow> <mn>2</mn> <mi>j</mi> <mo>-</mo> <mn>2</mn> </mrow> <mn>1</mn> </mmultiscripts> <mo>&amp;cup;</mo> <mmultiscripts> <mi>C</mi> <mrow> <mn>2</mn> <mi>j</mi> </mrow> <mn>1</mn> </mmultiscripts> <mo>&amp;cup;</mo> <mmultiscripts> <mi>C</mi> <mrow> <mn>2</mn> <mi>j</mi> <mo>+</mo> <mn>2</mn> </mrow> <mn>1</mn> </mmultiscripts> <mo>&amp;cup;</mo> <mn>...</mn> <mo>&amp;cup;</mo> <mmultiscripts> <mi>C</mi> <mrow> <mn>2</mn> <mi>n</mi> <mo>-</mo> <mn>2</mn> </mrow> <mn>1</mn> </mmultiscripts> <mo>&amp;cup;</mo> <mmultiscripts> <mi>C</mi> <mrow> <mn>2</mn> <mi>n</mi> </mrow> <mn>1</mn> </mmultiscripts> <mo>&amp;rsqb;</mo> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>34</mn> <mo>)</mo> </mrow> </mrow>
且满足：
<mrow> <mmultiscripts> <mi>D</mi> <mn>2</mn> <mi>n</mi> </mmultiscripts> <mo>&amp;CenterDot;</mo> <mmultiscripts> <mi>C</mi> <mn>2</mn> <mi>n</mi> </mmultiscripts> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>35</mn> <mo>)</mo> </mrow> </mrow>
其中，左上标“n”为组合梁中ETE-B的段数，左下标“2”为首段混合梁单元中TB梁单元序号；
步骤4、组合梁边界条件的数学描述：对于n段ETE-B构成的组合梁，其经典边界条件下的数学表达如下：
左端固支(CL)
左端简支(PL)
右端自由(FR)
右端简支(PR)
右端固支(CR)
步骤5、组合梁参数化频率方程的建立：两端经典边界条件(分别是式(41)、(42)描述的和式(43)～(45)描述的)的n段ETE-B构成的组合梁的振型特征矩阵D和振型系数向量C分别为：
和
<mrow> <msub> <mrow> <mo>&amp;lsqb;</mo> <mi>C</mi> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mo>(</mo> <mn>8</mn> <mi>n</mi> <mo>+</mo> <mn>4</mn> <mo>)</mo> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msub> <mrow> <mo>&amp;lsqb;</mo> <mmultiscripts> <mi>C</mi> <mn>2</mn> <mi>n</mi> </mmultiscripts> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mo>(</mo> <mn>8</mn> <mi>n</mi> <mo>+</mo> <mn>4</mn> <mo>)</mo> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>47</mn> <mo>)</mo> </mrow> </mrow>
且满足：
[D](8n+4)×(8n+4)·[C](8n+4)×1＝0. (48)
令组合梁的特征方程系数矩阵D的行列式(以特征函数f(λ)＝|D|表示)为零，求解参数化频率方程：
f(λ)＝|D|＝0 (49)
即获得整根组合梁的固有圆频率ω，按式f＝ω/2π计算组合梁的各阶固有频率f。
2.根据权利要求1所述的一种组合梁弯曲振动的固有频率分析方法，其特征在于：所述的混合梁单元ETE-B，将任意两个EB及其间的一个TB看作一个混合梁单元ETE-B，则组合梁结构可认为是多个混合梁单元ETE-B构成的；所有EB等截面、等长、空心，而所有TB等截面、等长、实心；组合梁的首尾两段梁单元都是EB梁。
3.根据权利要求1所述的一种组合梁弯曲振动的固有频率分析方法，其特征在于：所述的一维搜索法：以式(2)中EB频率参数λ为自变量，给定其取值范围、递增步长Δλ和初始值λ0，在λ数值递增形成的各个子区间中确定f(λ)变号的子区间[λk,λk+Δλ]，即满足f(λk)·f(λk+Δλ)＜0，其中λk＝λ0+k·Δλ，作为固有圆频率ω的可行域，在各个可行域内使用二分法，就能得到固有圆频率值，按式f＝ω/2π计算对应固有频率。
CN201510222482.XA 2015-05-04 2015-05-04 一种组合梁弯曲振动的固有频率分析方法 Active CN104778377B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510222482.XA CN104778377B (zh) 2015-05-04 2015-05-04 一种组合梁弯曲振动的固有频率分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510222482.XA CN104778377B (zh) 2015-05-04 2015-05-04 一种组合梁弯曲振动的固有频率分析方法

Publications (2)

Publication Number Publication Date
CN104778377A CN104778377A (zh) 2015-07-15
CN104778377B true CN104778377B (zh) 2017-11-10

Family

ID=53619837

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510222482.XA Active CN104778377B (zh) 2015-05-04 2015-05-04 一种组合梁弯曲振动的固有频率分析方法

Country Status (1)

Country Link
CN (1) CN104778377B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105138751A (zh) * 2015-08-13 2015-12-09 中国航空工业集团公司西安飞机设计研究所 一种结构动力学边界参数确定方法
CN105478351B (zh) * 2016-01-14 2017-07-21 河海大学常州校区 一种振动筛处理量实时检测系统及方法
CN106777825B (zh) * 2017-01-24 2020-03-27 湖南科技大学 一种基于谱有限元的矩形板振动模态计算方法
CN108416108B (zh) * 2018-02-05 2021-12-03 东南大学 用于钢-混凝土组合梁栓钉连接件设计的有限元建模方法
CN108509709B (zh) * 2018-03-28 2022-05-17 南京理工大学 双条裂纹fgm简支梁固有振型的数值计算方法
CN111310269B (zh) * 2020-02-25 2021-01-26 东南大学 一种考虑边界位移的固支梁结构非线性动特性分析方法
CN111611693B (zh) * 2020-04-27 2024-03-19 苏州科技大学 一种多段连续梁固有频率的计算方法
CN111797457B (zh) * 2020-07-15 2021-03-23 哈尔滨工业大学 一种钢-混组合连续梁桥自振频率快速估算方法
CN114580825B (zh) * 2021-12-09 2024-02-23 北京交通大学 一种基于数值分析的连接件组合梁分析系统
CN117171864B (zh) * 2023-11-02 2024-02-20 鲁东大学 一种梁结构线性振动预测方法
CN117494486B (zh) * 2024-01-03 2024-04-02 南通泰胜蓝岛海洋工程有限公司 一种局部集中载荷作用下的组合梁精细化应力位移分析方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6598717B1 (en) * 2000-06-09 2003-07-29 The Penn State Research Foundation Active-passive hybrid constrained layer for structural damping augmentation
US20040065152A1 (en) * 2002-09-27 2004-04-08 Hull Andrew J. Inverse method to estimate the properties of a flexural beam and the corresponding boundary parameters
CN102880807A (zh) * 2012-09-25 2013-01-16 中国海洋大学 适于梁弯曲振动分析的传递矩阵计算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6598717B1 (en) * 2000-06-09 2003-07-29 The Penn State Research Foundation Active-passive hybrid constrained layer for structural damping augmentation
US20040065152A1 (en) * 2002-09-27 2004-04-08 Hull Andrew J. Inverse method to estimate the properties of a flexural beam and the corresponding boundary parameters
CN102880807A (zh) * 2012-09-25 2013-01-16 中国海洋大学 适于梁弯曲振动分析的传递矩阵计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《Dynamics of Transversely Vibrating Beams Using Four Engineering Theories》;Seon M.Han 等;《Journal of Sound and Vibration》;19991231;全文 *
《混联弹簧减振系统下自同步直线振动筛的动力学特性》;彭利平 等;《煤炭学报》;20141030;第39卷(第10期);全文 *

Also Published As

Publication number Publication date
CN104778377A (zh) 2015-07-15

Similar Documents

Publication Publication Date Title
CN104778377B (zh) 一种组合梁弯曲振动的固有频率分析方法
CN104573202B (zh) 框架‑核心筒结构体系施工全过程的结构分析方法
CN104036095B (zh) 基于区域分解的耦合高精度复杂外形流场快速算法
CN109141820A (zh) 船体模型扭矩及剪力测量方法
CN114638046B (zh) 一种铁路桥墩数字孪生变截面仿真计算方法
CN102354343B (zh) 基于特征频率法计算超声导波频散关系的方法
CN105181327B (zh) 一种裂纹轮齿啮合刚度计算的方法
CN115081148B (zh) 一种基于势能理论的加筋板等效参数确定方法
CN104679955A (zh) 一种三角形网格加筋圆筒结构有限元参数化建模方法
CN101615219A (zh) 一种模拟输移扩散问题的高精度差分方法
CN104091065A (zh) 一种求解浅水问题模拟间断水流数值的方法
CN107012860A (zh) 一种基于bim技术的柱下单桩插筋定位装置的施工方法
CN104794332B (zh) 一种高层建筑风致响应分析模型的不确定性分析方法
Li et al. Adaptive extended isogeometric upper-bound limit analysis of cracked structures
Saraçoğlu et al. Linear static analysis of perforated plates with round and staggered holes under their self-weights
CN105893699A (zh) 一种含混合不确定性参数的连杆机构鲁棒综合方法
Pargana et al. Fully integrated design and analysis of Tensioned Fabric Structures: Finite elements and case studies
Wang et al. Experimental investigation of special-shaped concrete-filled square steel tube composite columns with steel hoops under axial loads
CN117171864B (zh) 一种梁结构线性振动预测方法
CN106802979A (zh) 基于有限元分析焊接机器人模型简化方法
Lé et al. Mixed dimensional modeling of reinforced structures
CN105005667B (zh) 劈裂注浆复合土体等效弹性模量及泊松比计算方法
Cai et al. A robust algorithm for the generation of integration cells in Numerical Manifold Method
CN112116713A (zh) 一种线性承台类构件高精度钢筋布置方法
CN107134007A (zh) 交通基础设施的三维快速建模方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 221116 Research Institute of China University of Mining and Technology,, Jiangsu

Applicant after: China University of Mining & Technology

Address before: 221116 Research Institute, China University of Mining and Technology, Xuzhou University, Jiangsu, China,

Applicant before: China University of Mining & Technology

CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Jiang Haishen

Inventor after: Dong Liang

Inventor after: Zhao Yuemin

Inventor after: Liu Chusheng

Inventor after: Peng Liping

Inventor before: Zhao Yuemin

Inventor before: Dong Liang

Inventor before: Liu Chusheng

Inventor before: Peng Liping

GR01 Patent grant
GR01 Patent grant