对有限元分析结果判定方法的探讨.pdfVIP

第 19 卷第 6 期 机械设计与研究 Vol . 19 No . 6 　 2003 年 12 月 Machine Design and Research Dec . ,2003 文章编号 (2003) 060 1303 对有限元分析结果判定方法的探讨 吴晓宇 , 　吴清文 , 　杨洪波 , 　何惠阳 ( 中国科学院　长春光学精密机械与物理研究所 ,长春 　130022) 　　摘 　要 :有限元法是一种工程近似求解的方法 ,存在误差 ,所以在用有限元软件进行工程分析计算时 ,很重要 的一点是分析过程中必须审视分析的结果 ,对结果有正确的认识 ,从而指导工程分析的正确进行 。采用有效的方 法对有限元分析结果的准确度进行估测 ,是十分必要的。举例说明了用公式法和模态分析法对有限元分析结果准 确度的判定 ,并根据经验列出了一些相关的判定方法 。 关键词 :有限元结果 ; 模态试验 ; 判定 中图分类号 : TP3 19 . 77 　　　文献标识码 :A 　　用有限元软件进行工程分析计算时 ,很重要的一点是分 析过程中必须审视分析的结果 ,对结果有正确的认识 ,从而 指导工程分析的正确进行 。有限元法是求解微分方程的一 种非常有效的数值分析方法 ,其基本思想是用分片函数去逼 近原函数 ,即把无限自由度问题转化为有限自由度的问题 , 再求解一个线性方程组 ,得到原方程的近似解 。既然是一种 ▲图 1 　圆截面悬臂梁 　　首先用材料力学公式进行计算 : 近似求解 ,就存在人为造成的误差 , 比如象网格划分疏密这 样的因素就会造成的不同程度的误差 ,或者说根据图纸建立 　　x = 1000mm时 ,最大截面计算应力 的有限元模型无法反映实件在加工和装配中产生的误差 ,也 σ = M = Fx = 100 ×1 ×32 = 15. 9 15MPa max W π 3 π 3 z d ×0. 04 会造成有限元的分析结果与实际不符 ,另外误差还与所分析 32 ( ) 的结构 最终表现为刚度矩阵 的性质有密切关系 ,例如当结 　　x = 1300mm时 ,最大截面计算应力 构在不同方向的刚度相差过于悬殊 ,可能使最后的代数方程 M Fx 100 ×1. 3 ×32 σ = = = = 20. 69MPa 组成为病态 ,从而使解的误差很大 ,可能导致结果失败 。当 max Wz π×d3 π×0. 043 32 然 ,还有的误差有可能是由于对软件的使用问题造成的 ,造 　　再用 M SC/ Nastran 软件进行计算 。 成的误差可能比较大 ,甚至远远背离了实际情况 ,这样就失 去了进行工程分析的意义[ 1 ] 。所以,采用有效的方法对工程 有限元分析中 ,当所选取的单元满足有限元解的收敛准 则 ,既完备又协调时 ,理论上讲 ,单元尺寸趋于零 ,则有限元 分析结果的准确程度进行估测 ,并且反过来指导工程分