基于8节点六面体单元的空间块体分析.pdf

如图4-23所示的一个空间块体，在右端部受两个集中力F作用，其中的参数为：

105

E110Pa,=0.25,t=0.2m,F=110N。基于MATLAB平台，用一个空间8节点六

面体单元计算各个节点位移、支座反力以及单元的应力。

(a)问题描述(b)有限元分析模型

图4-23右端部受集中力作用的空间块体

解答：对该问题进行有限元分析的过程如下。

（1）结构的离散化与编号

将结构离散为一个8节点六面体单元，节点编号如图4-23(b)所示，节点的几何坐标见

表4-13。

表4-13节点的坐标

节点坐标/m

节点

yz

10.200

20.20.80

300.80

4000

50.200.6

60.20.80.6

700.80.6

8000.6

节点位移列阵

q[uvwuvwuvw]eT（4-194）

111222888

(241)

总的节点载荷列阵

P[PPPPPPPPP]eT（4-195）

x1y1z1x2y2z2x8y8z8

(241)

其中，节点外载PPF110N5；支反力为PR，PR，PR，

z6z7x1x1y1y1z1z1

PRPRPRPRPRPRPRPR

，，，，，，，，

x4x4y4y4z4z4x5x5y5y5z5z5x8x8y8y8

PR；其余节点载荷分量为零。

z8z8

（2）计算单元的刚度矩阵(以国际标准单位)

首先在MATLAB环境下，输入弹性模量E和泊松比NU，然后针对题中单元节点坐标，

调用函数Hexahedral3D8Node_Stiffness，就可以得到单元的刚度矩阵k1(24×24)。

E=1.0e10;

NU=0.25;

l