BernoulliEuler薄壁梁双向弯曲与扭转耦合的振动分析.pdf

第十三届全国工程建设计算机应用学术会议论文集 广东佛山·2006 Bernoulli—Euler薄壁梁双向弯曲与 扭转耦合的振动分析 黄君毅王 勇 (华南理工大学建筑学院土木工程系) 【摘要】本文对任意截面的Bernoulli．Euler薄壁梁的弯扭耦合的振动微分方程组进行详细的数学推导， 得到其双向弯曲与扭转耦合的振动分析的精确解，并通过算例比较其与有限元软件计算的结果。 【关键词】Bernoulli．Euler薄壁梁；双向弯曲与扭转耦合 1前言 7在已有的梁动力响应问题的研究中【1’2～1，大多数理论假定其弯曲振动和扭转振动互不 耦合，或者仅有一个方向弯曲振动与扭转振动相耦合，这样使得梁的振动计算方便可行， 然而在实际问题中，却有不少梁的截面无对称可言，现有理论对这些梁的动力问题均无法 求解。 有限元法是计算上述振动问题的有效途径，但其缺陷是随着计算结构的自由度增加而 扭耦合的振动微分方程组15]进行详细的数学推导，得出了能精确有效计算双向弯曲与扭转 耦合振动问题的方法。 2动力方程组及其边界条件 ‘ 壁梁的振动微分方程： IEA“’一胆西=一吼 l目∥V：钔一pI螂V7+PAV，+肚z。p=q，一m： … l E1。w：们一∥。《+pal以一p4y，臼=q：+m； 、 “ 式中：E为弹性模量；A为截面面积；Y，、Z，为剪切中心坐标；P为材料密度： ，、L、露、，k分别为截面对Z轴和Y轴的惯性矩、剪切中心扇性惯性矩和抗扭惯性矩； U、v，、W。、秒分别为剪切中心的轴向位移、两个横向位移和转角位移；q。，qy，q：、m”my、 m：为三个方向均布力和局部扭矩；(，)、(·)分别为位移函数对X轴的求导和对时间t的求导。 第十三届全国工程建设计算机应用学术会议论文集 广东佛山·2006 图1 其中第一个方程是独立的，而后面三个方程则是相互耦合的，这三个方程就组成了 Bernoulli．Euler薄壁梁双向弯曲扭转耦合的自由振动微分方程。 各种边界条件如下： (a)自由边界：a2v／孤2=a3v／孤3=0，a2w／孤2=a3w／孤3=0， a2∥孤2=一日：a3∥巩3+GI々a纠孤=0； (2a) (2b) (b)固定边界：1，=av／0x=0，W=a叫孤生0，0=a纠a％=0； (c)简支边界：’，=a2v／融2=0，W=a2w／孤2=0，0=a2∥孤2=0。 (2c) 3特征方程的建立 CO，分离变量并引入过渡变量，则双向弯扭耦合的方程组化为： f(一无+无D4+D2)移一无磊痧=o {(一九+无D4+D2)协+疋弘痧=0，。 (3) I一九Z，哥+七y,Cv+(一A．D2+元D4+D2一五)0=0 式中：帚、访、0为三个方向的振型函数。 ：。 A／，≈，屯=rA／C，五=GI／po矿l，， 变量定义为：丸=譬A／I。，丸=酬Up02，覆=r L2GI 疋=￡Qw+I盔+A或+A乏、)|l‰，九W=ElWs f k，并AD暑Ld／dx 各种边界条件化为：