ANASYS有限元计算与材力公式计算结果比较
- 格式：doc
- 大小：820.00 KB
- 文档页数：16
ANASYS有限元计算与材力公式计算结果比较摘要：基于有限元单元法理论，使用ANASYS软件计算悬臂和两端固定两种梁在简单荷载作用下的位移与应力，并与使用材料力学公式计算的结果作比较，分析误差产生的原因，以加深对有限单元法的理解。
关键词：ANASYS；有限元；材料力学ANASYS FEM calculation and build formula results Abstract：Based on the theory of finite element method yuan, calculated using software ANASYS cantilever beam and two fixed ends in a simple load of displacement and stress, and the use of the mechanical formula for the results of comparative analysis of the reasons for the error, to deepen the understanding of the finite element method.Key words: ANASYS; finite element method; material mechanics1.前言有限单元法是当今工程分析中获得最广泛应用的数值计算方法，其分析的基本概念是用较简单的问题代替复杂问题后再求解。
它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适且较简单的近似解，然后推导求解这个域总的满足条件，从而得到问题的解。
这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。
由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。
ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件。
由世界上最大的有限元分析软件公司之一的美国ANSYS开发。
软件主要包括三个部分：前处理模块，分析计算模块和后处理模块。
前处理模块提供了一个强大的实体建模及网格划分工具，可以方便地构造有限元模型，分析计算模块包括结构分析（可进行线性分析、非线性分析和高度非线性分析）、流体动力学分析、电磁场分析、声场分析等，软件提供了100种以上的单元类型，用来模拟工程中的各种结构和材料。
本文基于有限元模拟悬臂梁和两端固定梁受集中力和均布力过程的分析材料应力和位移的方法。
通过比较材料力学计算所得的应力和位移，修正后得出最后结果。
经检验，ANSYS 有限元分析软件得到的应力和位移与材料实际的应力位移吻合得较好，因此，验证了有限元模型和材料模拟的正确性。
2.基本原理和公式变分原理实际是把求解偏微分方程边值问题转换为求解某一泛函的最小值问题。
有限元法是以变分原理和近似插值离散为基础的。
该方法首先利用变分原理把所要求解的边值问题转化为相应的变分问题，也就是泛函极值问题，然后利用对场域的网格剖分离散和在单元上对场函数的插值近似，将变分问题转化为普通多元函数的极值问题，最终归结为一个代数方程组，解之即得待求边值问题的数值解。
主要过程为：1.建立计算模型（结构离散化，划分单元，节点编号等）；2. 计算单元刚度矩阵；3. 形成整体刚度矩阵；4. 计算整体节点力；5. 引入位移边界条件；6. 求解整体平衡方程；7. 回代求解单元应力等对于不同物理性质和数学模型的问题,有限元求解法的基本步骤是相同的,只是具体公式推导和运算求解不同。
材料力学计算公式：（1）均布力作用下的悬臂梁， 2()()2q l x M x -=，()z M x y I σ= （2）集中力作用下两端固支梁，()82Fl Fx M x =- (0<x<2l ) 3()28Fx Fl M x =- (2l <x<l ) ()z M x yI σ=3.ANSYS 有限元计算在结构分析中， ANSYS 有限元分析软件对构件进行模拟，应考虑各种分析因素，选择合适的有限元单元类型，确定线性或非线性的问题等等。
ANSYS 分析GUIF 方式操作主要为三个步骤：1. 创建有限元模型。
2. 施加荷载进行求解。
3. 查看结果。
对于钢筋混凝土结构分析中单元划分通常基于两个层次：梁杆单元和实体单元。
前者着重分析单元力（包括力和弯矩）与位移（包括位移与转角）之间的关系，而后者着重于分析单元应力—应变分布，由于本次计算着重于位移和应力分析，所以采用PLANE42三角形三节点单元，创建二维平面单元。
图1给出了PLANE42单元的几何形状。
现在选取截面分别为1000mm ×300mm ，1000mm ×500mm ，1000mm ×600mm ，跨度为3000mm 的三根梁（悬臂梁和两端固定两种梁，其中悬臂梁施加5 KN/m 均布力，两端固定梁施加5 KN 集中力），即l/h ＝10，6，5，对这六根梁进行ANSYS 模拟，研究悬臂梁在均布力作用下，两端固支梁在集中力作用下的应力和位移，同时考虑跨高比的影响，由于是平面单元的模拟，梁厚度取单位1m ，最后输出应力和位移结果图如图2～图15所示。
4.材料力学数值计算1）悬臂梁 （1）应力每个截面边缘最大应力：()z M x W σ= 26z bh W =所以223()q l x bh σ-=（2）位移因为2()()2q l x M x -=-2''()()2qEIW M x l x =-=- 311'()23q EIW l x C =-⋅-+412()24qEIW l x C x C =-++在x=0 处，0W =,'0W =所以 316q C l = 4224qC l =-434()24624q q qEIW l x l x l =-+-当x=l 时W 最大， 48q W l EI= 所以截面尺寸为1000300⨯的梁：最大应力为 2222335(3000) 1.51000300ql bh σ⨯⨯===⨯ N/mm 2 最大位移为444312530000.75883101000300q W l EI ⨯⨯===⨯⨯⨯⨯ mm 截面尺寸为的1000500⨯梁：最大应力为 2222335(3000)0.541000500ql bh σ⨯⨯===⨯ N/mm 2 最大位移为444312530000.162883101000500q W l EI ⨯⨯===⨯⨯⨯⨯mm截面尺寸为的1000600⨯梁：最大应力为 2222335(3000)0.3751000600ql bh σ⨯⨯===⨯ N/mm 2最大位移为444312530000.0937883101000600q W l EI ⨯⨯===⨯⨯⨯⨯ mm 2）两端固定梁查结构力学可知：梁弯矩图如图：取半边结构（1）应力跨中最大应力z MW σ= 26z bh W =所以234FLbh σ=（2）位移因为()2()82l F x FlM x -=- ''()82Fl FxEIW M x =-=-211'84Fl EIW x Fx C =-+231211612Flx EIW Fx C x C =-++在x=0 处，0W =,'0W = 所以 10C = 20C = 2311612Flx EIW Fx =-当x=2l 时W 最大， 33319216F Fl W l EI Ebh ==所以截面尺寸为1000300⨯的梁：最大应力为234FLbhσ= =0.125 N/mm 2 最大位移为3316Fl W Ebh= =0.0104mm截面尺寸为的1000500⨯梁：最大应力为234FLbhσ= =0.045N/mm 2 最大位移为3316Fl W Ebh= =0.00225mm截面尺寸为的1000600⨯梁：最大应力为234FLbhσ= =0.03125 N/mm 2 最大位移为3316Fl W Ebh = =0.001303mm5.误差分析与比较将ANSYS 分析输出的结果与材料力学计算得的结果列表如下，注：1.应力单位为N/mm 2，位移单位为mm 。
2.跨度为3000mm ， l/h ＝10，6，5。
由上表所列数据，比较ANSYS分析结果和材料力学结果得出的误差可知：（1）最大误差为33.8％，不超过50％。
（2）均布力作用下的悬臂梁比集中力左右下的两端固支梁结果误差要小。
（3）对悬臂梁，应力结果比起位移结果误差偏大，相反，对两端固支梁，应力结果比位移结果误差小。
（4）在容许跨高比范围内，随着跨高比增加，位移误差越来越小。
应力误差趋之于大。
分析误差产生原因：有限元法分析的概念是用较简单的问题代替复杂问题，得到的是一个近似解，非准确解，存在误差是不可避免的。
基于有限元法的理论，用ANSYS软件分析结构，对结果的主要影响因素有：模拟的单元类型选择，节点选择，建模方式，网格划分精度和约束边界条件简化以及荷载施加条件，另外，网格质量，材料参数等都对结果有一定的影响。
在静力分析当中，网格划分越密，应力分析结果准确度越高。
对于简单的结构，经过ANSYS简化分析，得出的结果一般比较准确，与材料力学结果吻合得较好。
而对于较复杂的结构，材料力学计算与实际值相差就较大，ANSYS分析结果的准确性很大程度上就取决于建模和边界条件等是否合理了。
ANSYS分析的结果值误差通常在10％左右，如果不选择其他单元，而只选择在节点上施加荷载或应力集中处的单元，误差值有可能达到50％或以上。
针对本次计算分析，基于ANSYS误差估计，主要有：（1）悬臂梁比两端固定梁结构简单，因此应力与位移结果较准确，与材料力学结果相近。
（2）由于只进行ANSYS一次计算分析，没有重新划分网格进行二次分析比较可能网格划分精度不足，导致应力值误差较大。
（3）位移误差比应力误差小，悬臂梁位移结果与材料力学结果基本吻合，两端固支梁由于跨中作用一集中力，选取结果值的单元为施加荷载或应力集中附近处，所以误差偏离较大，但是也在可理解的合理范围内。
（4）由于选择的单元类型为三节点的三角形单元，此单元计算精度不高，一般结构分析中不推荐使用。
6.结论本文通过比较悬臂梁和两端固支梁的有限元分析和材料力学计算的应力、位移结果，得出以下结论：ANSYS通用有限元分析软件，具有极强的前后处理功能和较高的分析精度，并且使用方便，静力分析建模迅速快捷，简单结构的ANSYS分析结果和材料力学的结果非常接近，精度高。
参考文献【1】王勖成.有限单元法.北京：清华大学出版社，2003【2】李开泰.有限单元方法及其应用.陕西：西安交通大学出版社，1988【3】刘尔烈等.有限单元法及程序设计.天津：天津大学出版社，1999【4】刘涛等.精通ANSYS.北京：清华大学出版社，2002【5】郝文化.ANSYS土木工程应用实例[M].北京：中国水利水电出版社，2005【6】ANSYS中国公司.ANSYS建模及网络划分指南.2002【7】ANSYS结构有限元高级分析方法与范例应用.北京：中国水利水电出版社，2006。
