原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
第三章 平面问题的有限单元解法
对于一个工程问题的弹性体，如果具有某种特殊的形状，并且承受某种特殊的外力，就可以把问题简化为平面问题。本章主要介绍弹性力学平面问题有限元方法的全过程。
§3.1 平面问题三节点三角形单元的有限元格式
由于三角形单元对复杂的工程结构有较强的适应能力，因此容易将一个二维区域离散成有限个三角形单元，如图3－1所示。在边界上以若干直线近似曲线边界，随着单元增多，这种拟合将越精确，本节主要讨论了三结点三角形单元的有限元格式。
图3－1边界的直线近似
§ 3.1.1 单元位移模式及插值函数
如图3－2所示，三结点为：i, j, m。每个结点有两个位移分量（两个自由度），表示为
因此，三节点三角形单元有六个节点自由度。
在有限单元法中，单元的位移函数有各种选取方法，主要采用广义坐标法和自然坐标法。通过下面的具体分析，可以看到他们各自的特点。
因为选择多项式作为近似函数，微积分运算简便，并且随着项数的增加可以逼近任意一段光滑的函数曲线，故可以用多项式作为位移的近似函数。多项式的选取可由低次向高次选取。
三结点三角形单元的位移模式选取一次多项式时，
（3.1.1）
其矩阵形式为：
其中：，
这样，单元中的位移是坐标的线性函数，是待定系数，称之为广义坐标。六个广义坐标由单元的六个结点位移来表示。当取i结点时，则位移为：
同理可得和，则有：
（3.1.2）
即：
解上式，可得到由，和表示的广义坐标。系数行列式为：
(3.1.3)
A为三角形单元的面积。这样广义坐标为：
(3.1.4)
同样可以得到的表达式：
(3.1.5)
(3.1.4)和 (3.1.5) 式中：
(3.1.6)
上式中 表示坐标轮换，如。
将上面求得的广义坐标代入(2.1.1)式中，整理得到：
(3.1.7)
其中
(3.1.8)
在(3.1.7)式中是采用结点的位移来表示单元内任一点的位移，称为单元的插值函数或形函数。它是的一次函数，系数是由(3.1.6)式所示的三个节点坐标值来表示的。
(3.1.9)
(3.1.7)可写成矩阵形式：
(3.1.10)
其中
§ 3.1.2．应变矩阵和应力矩阵
当位移确定后，可以利用几何关系得到应变，进而由本构方程得到应力。对于平面问题，有：
(3.1.11)
称为应变矩阵，为平面问题微分算子，其中
(3.1.12)
由(3.1.8)式知：
(3.1.13)
代入(3.1.12)中，有：
(3.1.14)
则可以写为
(3.1.15)
其中是单元结点坐标的参数，当网格划分后，这些量则确定。因此对于采用线性函数的位移模式，B是常数矩阵。而一旦求出后，可得，也为常数。由此可知三结点三角形单元为常应变单元。在应变梯度较大的部位，应该尽量避免采用此种单元。在必须采用的情况下，则必须增加单元的密度，否则会导致较大的误差。
由物理方程，可以求得单元的应力
(3.1.16)
其中：
(3.1.17)
称之为应力矩阵。
在平面问题中，对于平面应力有：
(3.1.18)
如若是平面应变问题，上式中E应由替换，由替换。
利用(3.1.15) 和(3.1.18)式，可以得到：
文档评论（0）