现代机械设计师手册_下册_991-1100.docxVIP

10-96第10篇现代设计方法

作为示例，该平板分成8个平面三角形单元、9个节点。有限元求解具体过程如下：

1）建立平面直角坐标系Oxy，原点为图10.2-21所示节点1处，水平为x轴，垂直为y轴。列写节点坐标见表10.2-10所示。

表10.2-10节点坐标

节点号

1

2

3

4

5

6

7

8

9

x坐标值/m

0

5

10

0

5

10

0

5

10

y坐标值/m

0

0

0

10

10

10

20

20

20

列写单元编码如表10.2-11所示。

表10.2-11单元编码

单元号

1

2

3

4

5

6

7

8

节点1

1

2

1

2

4

5

4

5

节点2

2

3

5

6

5

6

8

9

节点3

5

6

4

5

8

9

7

8

2）计算各单元的单元刚度矩阵并进行扩展。首先计算出所需的系数b1、b2、b3、c1、c2、c3。对于单元1，三个节点对应的整体编码为（i，j，k）=（1，2，5），求得

类似地求得

b2=10，b3=0，c2=-5，c3=5

应变矩阵为

(-0.200.2000)

000

000-0.100.1

■■

■0-0.2-0.1100.10,

弹性力学平面问题的弹性矩阵为

得到单元1的单元刚度矩阵为

■■(2.30770-2.30770.34620-

■■

■0.5769,

对单元1的单元刚度矩阵进行扩展，得到一个18×18阶的方阵ke（x1t）。ke（x1t）只在上述（1，2，5）三个节点对应的元素上有值，其他元素上均为0。

全部8个单元均按上述同样的过程进行计算。

3）对上述8个单元的扩展刚度矩阵进行叠加，得到该结构的整体刚度矩阵为

4）在考虑位移约束条件的情况下列写结构节点位移列阵。在本示例中，节点1、2、3处均为全约束，即这三个节点的x、y方向对应的位移分量为0，即

=（u1，v1，u2，v2，u3，v3，u4，v4，…,u9，v9）T=（0，0，0，0，0，0，u4，v4，…,u9，v9）T

5）考虑结构的外载荷，构造结构载荷列阵。在

第2章有限元设计10-97

本示例中，只在节点9处作用水平和垂直载荷，因此可以得到

=（0，0，…,100，100）

6）根据本节内容引入边界条件，即根据约束情况修正结构有限元方程，特别是消除整体刚度矩阵的奇异性，得到考虑约束条件的、可解的有限元方程。

7）利用线性方程组的数值解法，对上述结构的有限元方程进行求解，得到所有各节点的位移矢量。最后根据需要求解单元应力。

求得节点位移（m）如下：

U=1.0×10-8×

（0000000.24590.15500.2302

0.04720.2190-0.06000.55770.2026

0.55640.08010.58190.0016）T