一维定长悬索问题的有限元法求解①.doc

一维定长悬索问题的有限元法求解①作者：居加颖 罗阳 江山来源：《科技创新导报》2011年第18期 ????????摘 要:本文讨论一维定长的悬索问题,即长度为L的悬索在均匀载荷力和弹性力作用下产生位移变化,并将其有限元解与真解进行比较,验证有限元解的有效性和精确性。 ????????关键词:有限元方法悬索问题刚度矩阵 ????????中图分类号:U448 文献标识码:A 文章编号:1674-098X(2011)06(c)-0037-01 ????????1 引言 ????????有限元方法[1-3]是求解各种复杂数学物理方程的重要方法,近半个世纪来在工程计算中的作用越来越重要。利用该方法可以获取复杂工程问题的各种信息,直接对工程设计进行具体评价,还能对工程中出现的问题进行具体详细的技术分析。如今有限元方法已经成为科学计算、工程设计、产品优化不可缺少的工具。 ????????2 一维定长悬索问题 ????????许多工程问题可以转化为与之等价的控制方程,通常由其基本方程和平衡方程给出。 ????????2.1 解析解 ????????设模型问题的基本模量为k,作用在悬索上垂直向上的外力为f,悬索沿y轴的正向向上偏移,设荷载f(x)和竖向位移v(x)均为正值,由位移理论知,悬索上的张力T为常数,将y方向上的所有力求和,有(1) ????????位移理论表明和,可得 ????????(2) ????????(2)除以,有(3) ????????(3)中的第一项,有(4) ????????采用相同的求极限过程可知 (5) ????????将(4)和(5)代入(3),得到控制方程 ????????(6) ????????设和为常数,此时方程为常系数线性方程,有指数形式的解。相应的齐次方程解为(7) ????????这里,且(6)的一个特解为(8) ????????联合(7)和(8),并由边界条件,得到问题的解析解 ????????(9) ????????其变分函数为: ????????(10) ????????2.2 有限元的单元刚度矩阵 ????????设横截面面积为常数,变分函数(10)关于形函数的矩阵形式为 ????????(11) ????????将(11)中函数关于取极小值,有 ????????. ????????(12) ????????一般情形下,由于悬索的横截面面积不会发生变化,此处删去面积项、(12)中第一项和第三项,即为 ????????第二项包括了基本模量和形状函数,相应矩阵可以写为 ????????积分并合并同类项,得到单元刚度矩阵 ????????悬索发生小位移的完整的单元刚度矩阵为 ????????. ????????3 数值例子 ????????3.1 有限元解 ????????将长度为L的悬索问题分为11个单元,两端固定且满足,。将已知值代入(10),可得单元刚度矩阵,即。将十一个单元刚度矩阵通过公共节点连接起来,得到12×12阶总刚度矩阵,其主对角线元素为,紧靠主对角线的相邻对角线以上(下)元素为,(其余元素为零),总右端向量元素为。 ????????引入边界条件,可以删去矩阵的第一行、列和第十二行、列,也可以将矩阵的第一行、列和第十二行、列的两个对角元素置为1,而将其中其他元素置为0;同时将右端项矩阵的第一项和第十二项置为0,通过计算机编程求得节点处的位移值。 ????????3.2 比较有限元解与真解 ????????相对误差定义为 | 有限元解-真解|/|真解|,则误差的范数为,各节点解的相对误差≤10-2。通过表1可以看出,有限元解非常好地逼近真解,其误差在允许范围内。 ????????注1．通过总刚度矩阵发现其稀疏性。若总的节点按顺序编号,对于总刚度矩阵而言,主对角线两边多个位置以外的零项没有影响,因此系数矩阵是带状的,只需要存储矩阵中位于上半带内的项。 ????????注2．该问题的单元刚度矩阵是对称的。在大多数的线性问题中,都出现对称的系数矩阵,则在计算的过程中只需计算对角线项和对角线以上的项。 ????????3 结语 ????????本文研究有限元方法在一维定长悬索实际问题背景下的简单应用,通过理论分析和数值试验证明了有限元法的精确性,获得的位移离散值非常有效地逼近了真解。 ????????参考文献 ????????[1] 邹仲康译.J. N. Reddy,有限元法概论[M].长沙:湖南科学技术出版社,1988. ????????[2] 李开泰,黄艾香,黄庆怀.有限元方法及其应用[M].北京:科学出版社,2006. 龙源期刊网