差分法下关于非均质变截面弹性直杆的固有频率的研究.docVIP

差分法下关于非均质变截面弹性直杆的固有频率的研究 研 究 生 课 程 论 文 (2013-2014学年第一学期) 差分法下关于非均质变截面弹性直杆的固 有频率的研究 差分法下关于非均质变截面弹性直杆的固有频 率的研究 刘延禄 摘要:采用差分法解非均质变截面弹性直杆的固有频率，并进行相应的研究，研究非均质、变截面弹性直杆的纵向自由振动问题，为其提供一种既能保证较高的精度，计算起来又较为简单的固有频率差分解法。此外，本文还给出了非均质变截面弹性直杆固有频率的计算公式。 关键字:数学物理方程，差分法 1 模型 已知一非均质变截面直杆长1m，其弹性模量服从函数E?E0e??，截面面积服 ???)从函数A?A0e(0.5，单位长度的质量为???0e0.5?，其中，E0?210GPa、 2、?0?78.5kg/m是杆左端对应的参数，??x/L，L是杆的长度。?100cmA0 求杆一端固定一段自由和两端固定情况下的固有频率。 2 纵向振动的微分方程 已知长为L的非均质变截面直杆，其弹性模量E，截面面积A，单位 杆长质量?均是杆位置的函数，其纵向的振动方程 为 ???u??2u ?E(x)A(x)???(x)2?0 (0?x?L) (1) ?x??x??t 式中t为时间，x为杆的纵向坐标，u是杆的位移，E(x)是杆在x处的弹性模量，A(x)是杆在x处的弹性模量，?(x)是杆在x处的单位长度质量。 设振动方程的解为 u(x,t)?U(x)?sin(?t??) (2) 式中U(x)为杆的纵向振型函数，?为纵向振动频率，?为相位常数。 引入无量纲变量? x (0???1) (3) L 将上述(2)式带入(1)式，并考虑(3)式，得到 ?? ??? ??u?22E(?)A(?)?l??(?)U(?)?0 (4) ?????? 设E，A，?是?的任意函数，即设 E(?)?E0F1(?) A(?)?A0F2(?) ?(?)??0F3(?) 将(4)式展开，得 d2U(?)dU(?)2 ?(?)?(?)U(?)?0 (5) kBB212 d?d? 式中: d(?)d(?)?? B1(?)??F1(?)2?F2(?)1?/F1(?)F2(?) d?d??? B1(?)?F3(?)/F1(?)F2(?) k??0?l/E0A0 2 22 3 差分方程和特征方程 采用差分法，将单位杆长分为n等份，网格点为n+1个，每份长1/n。 数值Ui用差分法公式可以得到 n 2 n2 (Ui?1?2Ui?Ui?1)?B1(?i)?(Ui?1?Ui?1)?B2(?i)kUi?0 (6) 2 其中B1(?i)、B2(?i)表示在i点的取值(i=0,1,2,3……，n) 于是可将上式简化为 Ui?1?P1(i)Ui?1?P2(i)Ui (7) 式中: 1 B1(?i)?n P1(i)? 1 n?B1(?i)2 ?B2(?i)k2n P2(i)? 2 n?2B1(?i) 2 2 引入位移列阵 ?Ui?1? ?i??? ???Ui? 传递列阵 1??0 Fi??? (i)(i)P2??P1于是(7)式可以写成如下形式 ??i?1????Fi????i? (8) 经传递运算，到下式 ??1?TT??? ??a?????Fi????b???1112???b? ???T21T22???i?n?1 (9) 式中:??a?和??b?分别为??a?和??a?的两个边界条件下对应的位移列阵 由此可得一下边界条件下对应的位移列阵和频率方程。 (1) 一端固定一端自由非均质变截面弹性直杆边界位移列阵和频率方程 ??0 U(0)?0 ??1 E(?)A(?) ?U(?) ?0 ?? 由上可得 Un?1?Un?1 对应的位移列阵 ?0? ??b????1???? ?U1? ?Un??Un? ??a????n?1??????? ?????Un?1??Un?1? ?Un???1?T11T12??0??? ?? ?????Fi????b????? ????i?n?1 ???T21T22??U1??Un?1? 得到 Un?T12U1 Un?1?T22U1 由(7)式得 Un?1?P1(n)Un?1?P2(n)Un 由于 ?Un?1Un?1 得到 ?(P(n)?1)T 1 22 ?P2(n)T12??U1?0 由于U1?0得对应的频率方程为 D(k)?(P1(n)?1)T22?P2(n)T12?0 (10) 2 (2) 两端固定非均质变截面弹性直杆边界位移列阵和频率方程 ??0 U(0)?0 ??1 U(1)?0