有限元方法及软件应用有限元平面问题3.doc

有限元方法及软件应用有限元平面问题3

FiniteElementMethod

3．面力的移置

设三角形单元某边界s上受面力q作用，分量为，，则

取ds则

由一般公式：

积分在边界s上

以上三种载荷的等效节点荷载由公式e导出

通常我们称：

为荷载移量的一般公式：

几点说明：

虚功等效静力等效。唯一性

一般

更多节点的单元公式形式不变，但不同

虽然公式e导出但对于面力和体力的计算都是很麻烦和困难的

N为x,y的函数，若p,q再为x,y的函数则更难，且单移分限不好定。

因此，我们将来还要进一步把这个问题解决好。

四．三角形单元的面积坐标（自然坐标，局部坐标）

1．面积坐标的定义：

图示三角形单元 I,j,k中任意一点m，其位置可由xoy坐标系中两个坐标来确定，即m(x,y)

若我们连接,,,则形成了3个小三角形ijm,ikm,jkm.

则有：若m(x,y)确定ijm,ikm,jkm.面积确定。

反之，ijm,ikm,jkm.面积确定m(x,y)确定

（用同底等高的概念解释！！）

因此，三角形单元内任一点可以

我们如何用三角形面积来描述m点的位置呢？

定义：节点I对边为底的三角形面积为;

节点j对边为底的三角形面积为;

节点k对边为底的三角形面积为;

设三角形单元的面积为A

令（2-37）

则三个比值，，称为三角形单元中m点的面积坐标.

2.三角形面积坐标的性质：

这就是三角形单元内m点作用有的等效节点荷载。只要计算出(I=I,j,k)即可。作为特例，考虑三角形单元形心处重力的移置。

形心坐标：

===0===-R

故：重力作用于形心时各节点均担。

体积力的移置

设单元作用有体力

则等效节点荷载为：

=

若为x,y的函数，则把用面积坐标表示（转换）

在常体力的作用下有：

===

即：常体力作用下，总体力均分三节点。

面力的移置。

设三角形单元I,j边上作用有梯形分布的面力q

由面力移置公式得：（可分别由节点合力表示及用节点分力表示）

==（q为合力，非分力）

则==

q为x,y的函数，把q表示面积坐标的函数有q=，在门边上是线性坐标，可利用两点式方程写出。

则：=====

同理：

=

==+

注意到在s上=0

=0

故：==

或：=

此法：1.避免复杂的分离。2.便于编程计算。

特例：若分布荷载为三角形分布。令（或）

则有：=（近端为2，远端为1）

说明：用以上积分的方法求等效节点荷载适用于任意节点的三角形单元，形函数也未必是线性的。

六．三角形单元节点荷载的形成

经过荷载处理后，我们已把非节点的荷载转化为常点荷载。

实际计算的荷载为：

计算荷载=原节点荷载+等效节点荷载

即：（2-49）

等效节点荷载要注意：

同时贡献的问题

用哪个单元计算的问题。

七．计算结果的整理：

有限元计算提供的结果一般为：1。节点位移2。单元应力

节点位移的处理：

一般把节点位移按比例标出，提供出结构常点位移分布规律

连成折线（线性位移函数）

连成光滑曲线（实际变形）

单元应力的处理：

输出的单元应力一般为，，（形心处）（三节点Ⅱ单元为常应力元，无所谓）

变换为单元的主应力。，,（材力）

变换为节点应力的主应力。（）（平均法）

Ex.节点5的应力为：

即：=（x,y,xy）

(,,)(,)

然后标出应力变化曲线。

计算结果的工作量随结构的单元，节点划分增加面增大。

要关注的是：1）位移的变化规律2）应力的最大值及发生地点。

小概念：位移最大的地方，应力未必最大。

八：有限元计算小结：

基本原理：

连续法有限个节点连接，有限大小的单元的组合法。

建立的节点位移为未知数，总刚为系的ｎ阶线性代数方程组。

研究方法

（确定）节点位移单元位移单元应变单元应力

单刚总刚计算荷载（等效节点荷载）

约束处理求解方程整理结果

３．解答特点：

１．假定单元内的位移分布规律，近似离散的数值解。

２．误差主要来源于：结构离散（连续离散），假定位移分布。

３．收敛性：单元缩小划分细密收敛于精确解。

Chap3.平面问题较精密单元的分析（矩形，高阶单元，等参单元）

3-1．问题的提出

在三角形单元中，我们假定位移函数是线性的。即：单元内的位移按线性规律变化。这是最简单，最基本的一个有限单元。而实际结构中在外载荷的作用下位移分布常非按线性变化。

设单元位移曲线为图示的f(x).显然，用线性插值解的精度较差。提示解的精度的方法：

增加单元，节点数（工作量大，费用高）

提高插