平面问题有限元.doc

《平面问题有限元.doc》由会员分享，可在线阅读，更多相关《平面问题有限元.doc（8页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、【精品文档】如有侵权，请联系网站删除，仅供学习与交流平面问题有限元.精品文档.CHAPTER02 平面问题的有限元21 引言l 平面应力和平面应变问题的有限元分析l 位移有限元建立的一般方法l 三角形单元、矩形单元l 线性单元、高阶单元l 轴对称单元22 平面问题的三角形单元格式221 平面应力和平面应变平面应力：和， 平面应变： 沿某一方向没有位移，只有和平面应力和平面应变的相互替换以下主要讨论平面应力222 三角形单元任意的一个平面采用三角形单元离散节点编号节点位移：如果节点位移已知，如何给出单元内部任意一点的位移？FEM的思想：单元内部插值223 单元分析若节点位移已知，如何求单元内部任

2、意一点的位移，称为位移模式或位移函数一般取为多项式，多项式选取：由低向高 为待定在节点上应满足节点位移可以求出，它应该是节点坐标和节点位移的函数。同理可以求出最后给出位移的形式为：或：或 由于N是线性函数，因此任意一点的位移就表现为单元节点位移（取决于）的线性组合（取决于N）或就是插值，因此N也叫做插值函数。整个有限元的核心就是构造N和N所具有的性质（也是今后构造N所必须满足的条件）（1） （2） 任意点插值函数的和为1 意味着如果单元三个节点位移相同（具有刚体位移），则内部任意一点也有这个位移。从而有：由于单元受相邻单元的影响，必然产生刚体位移，因此我们在构造N的时候必须使得位移模式包含刚体

3、位移。（3） 单元与单元边界上的位移是唯一的。保证唯一的连续性，否则单元之间出现间隙或重叠。（验证留做作业）有了位移后，应变和应力可以方便的给出注意由于是和的线性函数，因此是常数因此我们得到的常应变单元和常应力单元。作业21 验证上述三角形单元边界上位移的连续性。应变是否连续？224 有限元方程与刚度矩阵采用最小势能原理建立单元的有限元平衡方程如下：上述由三部分组成，一部分为体积力的贡献，另一部分为边界力的贡献，还有一部分为周边单元的作用， 认为是通过节点之间相互作用的。应用最小势能原理有：为刚度矩阵。整个有限元问题最后归结为求解上述方程单元刚度矩阵的的若干性质：（1） 单元的平衡方程： ，为

4、单元节点力 由三部分组成：体积力、边界力、单元相互作用力。（2）对称的矩阵（3）元素的性质：产生单位位移的力（4）的奇异性 单元的变形模式中，刚体位移没有消除（5）主元恒为正（可以直接验证）225 等效节点载荷列阵（1）自重（2）均布侧向力（3）其他的分布力226 有限元平衡方程每个单元， 都按照上述思想列出平衡方程整个结构当然也是处于平衡状态将此式写在一起， 并注意到右边单元之间的作用力互相抵消为零， 于是有：单元刚度阵在总体刚度阵中的位置K的带状形质,半带宽的计算(过去人工编号,现在可以实现自动最优排列)合理的编号使得半带宽尽可能小227 边界条件的施加将有限元平衡方程合成总体方程后，得到

5、其中为总体刚度矩阵，为总体未知节点位移列向量，为总体节点力列向量。由于总体方程是根据每一个节点的平衡得到的，因此中各元素的对应值应该是该节点所作用的等效外载荷。如果该节点没有外载荷的作用，则中的该元素为零，否则应该是各单元中该节点等效外载荷之和。由于到目前为此平衡方程中，位移施加一个刚体位移，并不影响平衡状态，因此需要有具体的位移边界约束，从而使得位移位移。以下讨论位移边界条件的施加方法。假设要求满足的位移边界为（1） 直接代入直接将代入总体平衡方程，将原来的总体自由度为的平衡方程，变为自由度为的平衡方程。代价： 刚度矩阵、节点位移、节点外载荷编号将打乱原有的顺序，不利于程序的编制。（2） 对

6、角元素为1如果对应的位移边界为零，即。此时可以将的元素改为1，而行和列的其他各元素为零，的第个元素也为0，从而实现所要求的边界条件。注意这一方法仅适应于的情况，事实上如果，则在刚度矩阵中必须考虑的影响，从而刚度矩阵中的行和列又不能为零。（3） 对角元素乘以大数其中为大数。此时原来的方程变为：由于很大，上式相当于。这一方法并不改变原有的编号顺序。作业2.2 请说明单元刚度矩阵的近似性原因.作业2.3 证明常应变单元发生刚体位移时,单元内不产生应力.作业2.4 对上图的平面结构由4个三角形单元组成,请对各节点编号,并写出总体刚度矩阵的半带宽.作业2.5 下图平面应力矩形板受均布载荷的作用,板厚为,左边刚性固定,右端为给定水平位移。在图示三角形三节点有限元网格情况下, 依照对角元素乘以大数法，写出总体平衡方程中的等效节点载荷列阵（刚度菊阵中的元素可用符号表示）