矩阵论大作业.docVIP

“矩阵论”课程研究报告 科 目： 矩阵理论及其应用 教 师： 姓 名： 学 号： 专 业： 机械设计及理论 类 别： 上课时间： 2013年2月至2013年5月 考 生 成 绩： 阅卷评语： 阅卷教师 (签名) 利用矩阵论相关知识求解传动轴 固有频率的有限元分析法 摘要：在结构力学中，求解结构自由振动的固有频率是十分重要的内容。本文通过对某机器传动轴各个单元进行单元刚度矩阵、单元质量矩阵等特性分析，再把各个单元的特性矩阵组集起来组成结构的总刚度矩阵、总质量矩阵，从而形成结构的自由振动方程式。最后利用矩阵论相关知识求解自由振动方程式的广义特征值，并通过广义特征值与固有频率的关系求得传动轴的固有频率。 正文 一、问题描述 已知某机器传动轴两端固定，其传动轴受扭长度L为1500mm，传动轴的横截面积是环形，其外径D为50mm，内径d为45mm，弹性模量E为。利用矩阵论及有限元分析法求解传动轴的固有频率。 二、方法简述 1.建立传动轴的有限元分析模型 由于传动轴两端固定，采用平面梁单元分析该传动轴。考虑到本次计算是手算，为了简化计算，将该传动轴划分为（1）、（2）两个单元，共1、2、3三个结点。由于该结构中一个结点有两个自由度，故总共有1、2、3、4、5、6六个自由度。建立有限元分析模型及各个部分编号如图1所示。 图1 传动轴有限元分析模型 2.平面梁单元的单元刚度矩阵 由《机械结构有限元分析》中形状函数N的构造方法可知，对于该结构的平面梁单元，它有两个节点，四个自由度，采用自然坐标系，通过构造计算可得单元的形函数为 （1） 其中为结点编号，v为结点位移，为结点转角，，为梁单元的长度。 平面梁单元的单元刚度矩阵 （2） 其中E为弹性模量，I为惯性矩。 而，通过对形状函数（1）式求两次导数得 （3） 显然，B也是x的函数对（2）式所有项积分得平面梁单元的单元刚度矩阵 （4） 故。 3.平面梁单元的单元质量矩阵 由（1）式我们得到了平面梁单元的形状函数，而单元质量矩阵公式为 （5） 其中梁单元密度，A为梁单元的横截面积。 将形状函数（1）式代入（5）式，对所有项积分得 （6） 故。 4.传动轴总刚度矩阵 按照有限元分析法中总体刚度矩阵组集的叠加方法，将两个单元刚度矩阵叠加后即可得到总刚度矩阵。由于该结构有1、2、3、4、5、6共六个自由度，故总刚度矩阵是六阶方阵。而由1、2、3、4号自由度组成，由3、4、5、6号自由度组成，故叠加时应自由度对号入座相互叠加。叠加所得总刚度矩阵为 （7） 5.传动轴总质量矩阵 按照有限元分析法中总质量矩阵组集的叠加方法，将两个单元质量矩阵叠加后得到总质量矩阵，叠加方式与获取总刚度矩阵叠加方式一样，自由度对号入座相互叠加。故可得到总质量矩阵为 （8） 6.振幅列阵 因为传动轴总共有三个结点，结点处有1、2、3、4、5、6六个自由度，分别是，故可以得到振幅列阵为 （9） 7.传动轴自由振动方程式 在有限元分析中，结构自由振动的方程式推导为 （10） 由于支撑条件给出了，故在总刚度矩阵、总质量矩阵、振幅列阵中应删掉所对应的行和列，也就是自由度1、2、5、6所对应的行和列。这样，我们就得到了简化的