计算力学结课论文平面框架结构有限元分析.docVIP

PAGE PAGE 1平面框架结构有限元分析1 引言1.1 有限元法的重要性有限元方法是求解各种复杂数学物理问题的重要方法，是处理各种复杂工程问题的重要分析手段，也是进行科学研究的重要工具。利用有限元分析可以获取几乎任意复杂工程结构的各种机械性能信息，还可以直接就工程设计进行各种评判，可以就各种工程事故进行技术分析。1.2 有限元法的发展过程有限元的思想可以追溯到20世纪40年代。1941年A.Hrennikoff首次提出用离散元素法求救弹性力学问题，当时仅限于用杆系结构来构造离散模型，但能很好地说明有限元的思想。如果原结构是杆系，这种方法的解是精确的，发展到现在就是大家熟知的矩阵分析法。究其实质这还不能说就是有限元的思想，但结合以后的有限元理论，统称为广义有限单元法。1943年R.Courant在求解扭转问题时为了表征翘曲函数而将截面分成若干三角形区域，在各三角形区域设定一个线性的翘曲函数，这实质上就是有限单元法的基本思想。20世纪50年代因航空工业的需要，美国波音公司的专家首次采用三节点三角形单元，将矩阵位移法用到平面问题上。1960年美国的R.W.Clough教授在一篇题为“平面应力分析的有限单元法”的论文中首先使用“有限单元法”一词，此后这一名称得到广泛承认。20世纪60年代有限元法发展迅速，除力学界外，许多数学家也参与了这一工作，奠定了有限元法的理论基础，搞清了有限元法与变分法之间的关系，发展了各种各样单元模式，扩大了有限元法的应用范围。20世纪70年代以来，有限元法进一步得到蓬勃发展，其应用范围扩展到所有工程领域，成为连续介质问题数值解法中最活跃的分支。由变分法有限元扩展到加权残数法与能量平衡法有限元，由弹性力学平面问题扩展到空间问题、板壳问题，由静力平衡问题扩展到稳定性问题、动力问题和波动问题，由线性问题扩展到非线性问题，分析的对象从弹性材料扩展到塑性、粘弹性、粘塑性和复合材料等，由结构分析扩展到结构优化乃至于设计自动化，从固体力学扩展到流体力学、传热学、电磁学等领域。它使许多复杂的工程分析问题迎刃而解。1.3 有限元法的基本思想有限元法的基本思想是将物体（即连续的求解域）离散成有限个且按一定方式相互联结在一起的单元的组合，来模拟或逼近原来的物体，从而将一个连续的无限自由度问题简化为离散的有限自由度问题求解的一种数值分析方法。物体被离散后，通过对其中各个单元进行单元分析，最终得到对整个物体的分析。网格划分中每一个小的块体称为单元。确定单元形状、单元之间相互联结的点称为节点。单元上节点处的结构内力为节点力，外力（有集中力、分布力等）为节点载荷。2 理论分析如图1所示框架结构，其顶端受均布力作用，用有限元方法分析该结构的位移。结构中各个截面的参数都为：E=3.0×1011Pa,I=6.5×10-7 m4 ，A=6.8×10-4 m2。 图1 框架结构受一均布力作用对该问题进行有限元分析的过程如下：（1）结构的离散化与编号将该结构离散为三个单元，节点位移及单元编号如图2所示。图2 节点位移及单位编号节点位移矩阵为节点外载矩阵为支反力矩阵为其中：Rx3、Rx4、Rθ3为节点3的沿x方向支反力、沿y方向支反力和支反力矩，Rx4、Ry4、Rθ4为节点4的沿x方向支反力、沿y方向支反力和支反力矩，均为代求值。总的节点载荷矩阵为（2）各个单元的描述一般而言，一个框架单元至少需要两个参照系，即一个整体坐标系和一个局部坐标系。整体坐标系（X,Y）作用如下：（1）描述每个连接点（节点）的位置，用角度表示每个单元的方向；（2）按照整体坐标系下各单元的自由度分量施加约束和载荷；（3）表示问题的解。此外，还要用局部坐标系描述单元的轴向效应。局部自由度（指位移和转角）通过转换矩阵与整体自由度联系起来。单元1的局部坐标与整体坐标是一致的，则可以得到单元1的刚度矩阵为单元2和单元3的情况相同，只是节点编号不同而已，其局部坐标系下的单元刚度矩阵为这两个单元轴线的方向余弦为，有坐标转换矩阵为则可以计算出整体坐标下的单元刚度矩阵（单元2和单元3）注意这两个单元所对应的节点位移矩阵分别为对于单元2：对于单元3：（3）建立整体刚度矩阵方程组装整体刚度矩阵并形成整体刚度方程，其中刚度矩阵的装配关系为（4）边界条件的处理及刚度方程的求解该问题的位移边界条件为，处理该边界条件后的刚度方程为=求解后的结果为3 ANSYS分析 ANSYS软件提供了结构静力分析、结构动力学分析、结构非线性分析、动力学分析、热分析、电磁场分析、流体动力学分析、声场分析、压电分析等9种分析类型。本例为平面框架结构，进行结构静力分析，得结果如下：图3 生成的几何模型结果显示图4 施加位移边界条件及载荷的结果显示图5 变形前和未变形的几何形状显示图