[发明专利]一种组合梁弯曲振动的固有频率分析方法

权利要求书

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)的一半标示，即两类梁单元在交界面F_i和F_i+1处通过双边V缝周边对焊联接，材料体密度均为ρ，弹性模量均为E；TB-i+1中，梁段宽b、高h、长l_i+1，满足h/l_i+1∈[1/5,+∞]，截面面积A_i+1，相对于x轴的二次惯性矩为I_i+1，截面的剪切修正系数为κ，材料的剪切弹性模量为G；EB-i和EB-i+2中，梁段宽均为b、高均为h、长分别为l_i和l_i+2；所述的l_i+2为l_i＝l_i+2，壁厚均为h₀，满足h/l_i和h/l_i+2∈[0,1/5]，截面面积分别为A_i和A_i+2，所述的A_i+2为A_i+1＝A_i+2，相对于x轴的二次惯性矩分别为I_i和I_i+2；所述的I_i+2为I_i+1＝I_i+2；

步骤2、混合梁单元ETE-B的动力学分析：在Euler-Bernoulli梁单元EB-i的固有圆频率ω满足：

λ4=ρAiEIiω2---(2)]]>

式(2)中λ为EB频率参数，且有振型函数：

Y_i(x)＝e_i,1sinλx+e_i,2cosλx+e_i,3sinhλx+e_i,4coshλx.(3)

式(3)中e_i,1,e_i,2,e_i,3和e_i,4为EB-i梁单元边界条件决定的系数；

对于其他等截面等长Euler-Bernoulli梁单元，只需对应更改梁单元序号i；

由于是同一根组合梁的不同单元，Timoshenko梁单元TB-i+1的固有圆频率与步骤2中EB-i的固有圆频率ω相同，在满足

ω<2Ai+1Eκ2G2ρIi+1(E+κG)---(6)]]>

的前提下，有一组虚根：

r1,2=±r‾·i,---(7)]]>

式(7)中为TB虚频率参数，i为虚数单位，满足

r‾=2κGEIi+1ω(κGρIi+1ω+EρIi+1ω+ρ2Ii+12ω2(κG-E)2+4κ2G2ρAi+1EIi+1)/(2κGEIi+1),---(8)]]>

而还有一组实根

r3,4=±r~,---(9)]]>

式(9)中为TB实频率参数，满足

r~=-2κGEIi+1ω(κGρIi+1ω+EρIi+1ω-ρ2Ii+12ω2(κG-E)2+4κ4G2ρAi+1EIi+1)/(2κGEIi+1).---(10)]]> 1

TB-i+1的弯曲变形振型表达式：

Yi+1(x)=ti+1,1sin(r‾x)+ti+1,2cos(r‾x)+ti+1,3sinh(r~x)+ti+1,4cosh(r~x).---(11)]]>

式中，t_i+1,1,t_i+1,2,t_i+1,3和t_i+1,4为TB-i+1梁单元边界条件决定的系数；

TB-i+1的截面转角振型表达式：

φi+1(x)=-(r‾-ρω2κGr‾)ti+1,2·sin(r‾x)+(r‾-ρω2κGr‾)ti+1,1·cos(r‾x)+(r~+ρω2κGr~)ti+1,4·sinh(r~x)+(r~+ρω2κGr~)ti+1,3·cosh(r~x).---(12)]]>

对于其他等截面等长Timoshenko梁单元，只需对应更改梁单元序号i+1；

步骤3、混合梁单元ETE-B联接条件的数学描述：定义EB系数和TB系数则在交界面F_i＝(EB-i)∩(TB-i+1)处，存在振型系数界面联接向量：

和振型界面联接特征矩阵

类似地，在交界面F_i+1＝(TB-i+1)∩(EB-i+2)处存在振型系数界面联接向量

和振型界面联接特征矩阵

因此，任意单段混合梁单元的振型联接特征矩阵和振型系数联接向量分别为：

[Di+11]8×12=[DiF]4×8[0]4×4[0]4×4[Di+1F]4×8---(27)]]>

和

[Ci+11]12×1=[CiF∪Ci+1F]=EiFTi+1FEi+2F,---(28)]]>

且满足：

Di+11·Ci+11=0.---(29)]]>

左上标为ETE-B的段数，左下标为首段组合梁中TB梁序号，

对于任意连续两段ETE-B，TB梁单元序号为i+1和i+3，共同拥有EB-i+2，振型联接特征矩阵和振型系数联接向量分别为：

和

[Ci+12]20×1=[Ci+11∪Ci+31]=[CiF∪Ci+1F∪Ci+2F∪Ci+3F]=EiFTi+1FEi+2FTi+3FEi+4F,---(31)]]>

且满足：

Di+12·Ci+12=0,---(32)]]>

其中，和分别由和将其各含右下标元素的右下标数值上加2获得，未给出

大小的“0”的维数视矩阵整体而定；

对于n段ETE-B构成的组合梁，TB梁单元序号依次为2,4,6,…,2j,…,2n，振型联接特征矩阵和振型系数联接向量分别为：

和

[C2n](8n+4)×1=[C21∪C41∪...∪C2j-21∪C2j1∪C2j+21∪...∪C2n-21∪C2n1],---(34)]]>

且满足：

D2n·C2n=0,---(35)]]>

其中，左上标“n”为组合梁中ETE-B的段数，左下标“2”为首段混合梁单元中TB梁单元序号；

步骤4、组合梁边界条件的数学描述：对于n段ETE-B构成的组合梁，其经典边界条件下的数学表达如下：

左端固支(CL)

左端简支(PL)

右端自由(FR)

右端简支(PR)

右端固支(CR)

步骤5、组合梁参数化频率方程的建立：两端经典边界条件(分别是式(41)、(42)描述的和式(43)～(45)描述的)的n段ETE-B构成的组合梁的振型特征矩阵D和振型系数向量C分别为：

和

[C](8n+4)×1=[C2n](8n+4)×1,---(47)]]>

且满足：

[D]_{(8n+4)×(8n+4)}·[C]_(8n+4)×1＝0.(48)

令组合梁的特征方程系数矩阵D的行列式(以特征函数f(λ)＝|D|表示)为零，求解参数化频率方程：

f(λ)＝|D|＝0 (49)

即获得整根组合梁的固有圆频率ω，按式f＝ω/2π计算组合梁的各阶固有频率f。