matlab patch 六面体,《有限元基础教程》_【MATLAB算例】4.8.2(1) 基于8节点六面体单元的空间块体分析(Hexahedral3D8Node)...

【MATLAB 算例】4.8.2(1) 基于8节点六面体单元的空间块体分析(Hexahedral3D8Node)

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

105110Pa,=0.25,=0.2m,=110N E t F μ=??。基于MATLAB 平台，用一个空间8节点六

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

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

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

解答：对该问题进行有限元分析的过程如下。 (1)结构的离散化与编号

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

表4-13 节点的坐标

节点 节点坐标/m x y z 1 0.2 0 0 2 0.2 0.8 0 3 0 0.8 0 4 0 0 0 5 0.2 0 0.6 6 0.2 0.8 0.6 7 0 0.8 0.6 8

0.6

节点位移列阵

T

1

11

2

2

2

8

88(241)

[]e u v w

u

v

w

u

v

w ?=q (4-194)

总的节点载荷列阵

T 111222888(241)

[]e x y z x y z x y z P P P P P P P P P ?=P (4-195)

其中，节点外载5

67110z z P P F N ==-=-?；支反力为11x x P

R =，11y y P R =，11z z P R =，