六面体单元matlab后处理,《有限元基础教程》_【MATLAB算例】基于节点六面体单元的空间块体分析(HexahedralDNode).doc...

《有限元基础教程》_【MATLAB算例】基于节点六面体单元的空间块体分析(HexahedralDNode)

【MATLAB算例】4.8.2(1) 基于8节点面体单元。基于MATLAB平台，用一个空间8节点六面体单元计算各个节点位移、支座反力以及单元的应力。

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

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

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

(1)结构的离散化与编号

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

表4-13 节点的坐标

节点节点坐标/mxyz10.20020.20.80300.80400050.200.660.20.80.6700.80.68000.6节点位移列阵

(4-194)

总的节点载荷列阵

(4-195)

其中，节点外载；支反力为，，，，，，，，，，，；其余节点载荷分量为零。

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

首先在MATLAB环境下，输入弹性模量E和泊松比NU，然后针对题中单元节点坐标，调用函数Hexahedral3D8Node_Stiffness，就可以得到单元的刚度矩阵k1 (24×24)。

>>E=1.0e10;

>>NU=0.25;

>>lx=0.2;

>>ly=0.8;

>>lz=0.6;

>>k1