有限元法绪论(已排)

《有限元法绪论(已排)》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第2部分有限元分析及应用FiniteElementAnalysisandApplications1第6章　有限元法的基本概念2在工程技术领域内，经常会遇到两类典型的问题。第一类问题，可以归结为有限个已知单元体的组合。例如，材料力学中的连续梁、建筑结构框架和桁架结构。这类问题称为离散系统。如下图所示平面桁架结构，是由6个承受轴向力的“杆单元”组成。6.1工程和科学中典型问题3第二类问题，通常可以建立它们应遵循的基本方程，即微分方程和相应的边界条件。例如弹性力学问题，热传导问题等。由于建立基本方程所研究的对象通常是无限小的单元，这类问题称为连续系统，或场问题。尽管已经建立了连续系统的基本方程，由于

2、边界条件的限制，通常只能得到少数简单问题的精确解答。对于许多实际的工程问题，还无法给出精确的解答。为解决这个困难，工程师们和数学家们提出了许多近似方法。6.1工程和科学中典型问题46.2场问题的一般描述实例：二维热传导（稳态）问题原理：从两个方向传入微元体的热量与微元体内热源产生的热量Q平衡基本方程：边界条件：56.3场问题的求解策略及方法6.3.1求解策略1、直接法：求解基本方程和相应定解条件的解；2、间接法：基于变分原理，构造基本方程及相应定解条件的泛函形式，通过求解泛函的极值来获得原问题的近似解。即将微分形式转化与其等价的泛函变分的积分形式。6.3.2求解方法1、解析或半解析法：2、数值

3、法：A）基于直接法的数值法，如差分法；B）基于间接法的数值法，如等效积分法（如里兹法）、有限元法等。6特点优缺点差分法均匀离散求解域；差分代替微分；解代数方程组要求规则边界，几何形状复杂时精度低等效积分法（加权余量法或泛函变分法）整体场函数用近似函数代替；（近似函数常为含n个待定系数的多项式，）微分方程及定解条件的等效积分转化为某个泛函的变分，--求极值问题，（利用极值条件建立n个代数方程），解代数方程组适合简单问题，复杂问题很难解决有限元法可非均匀离散求解域；分片连续函数近似整体未知场函数；解线性方程组。有限元法的数学基础仍是变分法（同上）。节点可任意配置，边界适应性好；适应任意支撑条件和载

4、荷；计算精度与网格疏密和单元形态有关，精度可控。对裂缝和无限域的分析存在不足数值计算方法分类7先将求解域离散为有限个单元，单元与单元只在节点相互连接；----即原始连续求解域用有限个单元的集合近似代替每个单元选择一个简单的场函数近似表示真实场函数在其上的分布规律，该简单函数可由单元节点上物理量来表示----通常称为插值函数或位移函数基于问题的基本方程，建立单元节点的平衡方程（即单元刚度方程）借助于矩阵表示，把所有单元的刚度方程组合成整体的刚度方程，这是一组以节点物理量为未知量的线形方程组，引入边界条件求解该方程组即可。6.4有限元法基本思想8整体平衡分片近似单元平衡结构离散方程求解问题分析力学

5、模型节点单元位移函数单刚方程总刚方程节点位移有限元法基本思想9节点位移向量表示：节点力向量表示：节点1沿x方向的位移　、其余节点位移全为0时轴向压力为：实例1:(1)求右图离散结构2的点位移10同理，节点2作用于单元1上的力，其大小与之相等，方向相反，x和y方向的分量分别记为：注：表示第e个单元的第j个自由度产生单位位移，而其它自由度上的位移为零时，第i个自由度上所受的力。常称其为单元的刚度系数。实例1:(2)单元分析节点1作用于单元1上的力，在x和y方向的分量分别为：11单元2节点力平衡方程实例1:(2)单元分析同理可求分别作单位位移时相应的刚度系数，考虑到节点的实际受力为和实际位移为　　　

6、　　　，则据各个节点节点力平衡得：12结合前式推导得：实例1:(3)整体分析整体分析：作用于每个节点上的节点力平衡，即13整体矩阵记为：求解上述整体方程，可得问题的节点位移。实例1:(4)引入约束求解将　代入可得整体方程14实例2连续问题例：求等截面直杆在自重作用下的拉伸。图(a)中单位杆长重量为q，杆长为L，截面面积为A，弹性模数为E。15实例2材料力学方法求解直杆拉伸：考虑微段dx，内力N=q(L-x)dx的伸长为：x截面上的位移：根据几何方程求应变，物理方程求应力。这里应变：应力：16iL1iL+图2-3i+1ii-12)LL(q1ii++1、离散化2、外载荷集中到结点上，即把阴影部分的

7、重量作用在结点i上实例2（1）结构离散有限单元法求解直杆拉伸：直接公式法17实例2（2）单元分析3、假设线单元上的位移为线性函数184、以i结点为对象，列力的平衡方程令将位移和内力的关系代入得用结点位移表示的平衡方程，其中i=1，2，…n有n个方程未知数也有n个，解方程组，得出结点位移，进而计算应力。实例2（2）单元分析19假设线单元数为3个的情况，平衡方程有3个：i=1时，i=2时,i=3时，联