矩阵论大作业.docxVIP

4彦靖 “矩阵论”课程研究报告 科 目:矩阵理论及其应用 教 师: 名： 学 号: 专 业：机械设计及理论 类 别: 上课时间： 2013年2月至2013年5月 考生成绩： 阅卷评语： 阅卷教师（签名） 利用矩阵论相关知识求解传动轴 固有频率的有限元分析法 摘要：在结构力学中，求解结构自由振动的固有频率是十分重要的内容。本文通 过对某机器传动轴各个单元进行单元刚度矩阵、单元质量矩阵等特性分析，再把 各个单元的特性矩阵组集起来组成结构的总刚度矩阵、总质量矩阵，从而形成结 构的自由振动方程式。最后利用矩阵论相关知识求解自由振动方程式的广义特征 值，并通过广义特征值与固有频率的关系求得传动轴的固有频率。 正文 一、 问题描述 已知某机器传动轴两端固定，其传动轴受扭长度L为1500mm,传动轴的横截 面积是环形，其外径D为50mm,内径d为45mm,弹性模量E为2.1x10sA^/mm2 0利 用矩阵论及有限元分析法求解传动轴的固有频率。 二、 方法简述 1?建立传动轴的有限元分析模型 由于传动轴两端固定，采用平面梁单元分析该传动轴。考虑到本次计算是手 算，为了简化计算，将该传动轴划分为（1）、（2）两个单元，共1、2、3三个结 点。由于该结构中一个结点有两个自由度，故总共有1、2、3、4、5、6六个自由 图1传动轴有限元分析模型 显然，B也是x 显然，B也是x的函数对 （2）式所有项积分得平面梁单元的单元刚度矩阵 (2) (3) 2.平面梁单元的单元刚度矩阵 由《机械结构有限元分析》⑴中形状函数N的构造方法可知，对于该结构的平 面梁单元，它有两个节点，四个自由度，采用自然坐标系，通过构造计算可得单 元的形函数为 N=［（NJ （N） （N） （N）； = ［（1 —3 孑+ 2f） （§-2F + F）（3孑一2孑）（孑一孑）］ （1） 其中为结点编号，v为结点位移，为结点转角，年，I为梁单元的长度。 平面梁单元的单元刚度矩阵 k = E/j BTBdx I 其中E为弹性模量，I为惯性矩。 而B = 驾,通过对形状函数（1）式求两次导数得 d* B = ［H{ B2 B3 B4］ z 6 z 4 /X、 z 6 Z 2 , X、 =（7+127t）（7+6F}（—正一“产）（—7+6正） (4) (4) 12 61 -12 6/ 6/ 4/2 -6/ 2卩 -12 -61 12 -61 6/ 2/2 -6/ 4/2 EI ：=下 3.平面梁单元的单元质量矩阵 由（1）式我们得到了平面梁单元的形状函数，而单元质量矩阵公式为 m = N pNdV = pAl^ NrNd^ 其中p梁单元密度，A为梁单元的横截面积。 将形状函数（1）式代入（5）式，对所有项积分得 _156 22/ 54 -13/ pAl m = 221 4/2 13/ -3Z2 (6) 420 54 13/ 156 -22/ -13/ -312 -22/ 4/2 故加⑴=m ⑵=m o 4.传动轴总刚度矩阵 按照有限元分析法中总体刚度矩阵组集的叠加方法，将两个单元刚度矩阵叠 加后即可得到总刚度矩阵。由于该结构有1、2、3、4、5、6共六个自由度，故总 刚度矩阵是六阶方阵。而￡⑴由1、2、3、4号自由度组成，代2）由3、4、5、6号自 由度组成，故叠加时应自由度对号入座和互叠加。叠加所得总刚度矩阵为 1 2 3 4 5 6 1 ■ 12 6/ -12 6/ 0 0 ■ 2 6/ 4/2 -6/ 2/2 0 0 3 EI -12 -6/ 24 0 -12 61 (7) K_ ： 2/2 8/2 2/2 4 /3 6/ 0 -61 5 0 0 -12 -6/ 12 -6/ 6 0 0 6/ 2/2 -6/ 4/2 5.传动轴总质量矩阵 按照有限元分析法中总质量矩阵组集的叠加方法，将两个单元质量矩阵叠加 后得到总质量矩阵，叠加方式与获取总刚度矩阵叠加方式一样，自由度对号入座 和互叠加。故可得到总质量矩阵为 「156 221 54 -13/ 0 0 _ 22/ 4厂 13/ - 3广 0 0 M = P 54 13/ 312 0 54 -13/ (8) 420 -13/ -引2 0 8/2 1引 -3/2 0 0 54 13/ 156 -22Z 0 0 -13/ 一 3广 -22/ 4/2 _ 6.振幅列阵 因为传动轴总共有三个结点，结点处有1、2、3、4、5、6六个自由度，分别 是儿01，卩2 02,303 故可以得到振幅列阵为 vl 厂V2 叫 7.传动轴自由振动方程式 在有限元分析中，结构自由振动的方程式推导为 Kq = OTMq (10) 由于支撑条件给出了儿=乞=叫=仇3=0,故在总刚度矩阵、总质量矩阵、振幅列 阵中应删掉片0|,叫,彳3所对应的行和列，也就是自