【2017年整理】_基于8节点六面体单元的空间块体分析.docVIP

如图4-23所示的一个空间块体，在右端部受两个集中力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 0 0.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=Hexahedral3D8Node_Stiffness(E,NU,lx,0,0,lx,ly,0,0,ly,0,0,0,0,lx,0,lz,lx,ly,lz,0,ly,lz,0,0,lz); （3） 建立整体刚度方程 由于该结构共有8个节点，则总共的自由度数为24，因此，结构总的刚度矩阵为K K(24×24)，先对KK清零，然后调用函数Hexahedral3D8Node_Assembly进行刚度矩阵的组装。由于本题中只用了一个单元，因此总体刚度矩阵KK与单元刚度k1相同，此处不再列出，调用函数的过程如下： KK=zeros(24,24); KK=Hexahedral3D8Node_Assembly(KK,k1,1,2,3,4,5,6,7,8); （4） 边界条件的处理及刚度方程求解 由图4-23(b)可以看出，节点1，4，5和8的3个方向的位移将为零，即。因此，将针对节点2，3，6和7的位移进行求解，这4个节点的位移将对应KK矩阵中的4～9行、4～9列，4～9行、16～21列，16～21行、4～9列，以及16～21行、16～21列，则需从KK(24×24)中提取相应行和列的数据，置给k，然后生成对应的载荷列阵p，再采用高斯消去法进行求解，即MATLAB中的反斜线符号“\”求解。 k=[KK(4:9,4:9),KK(4:9,16:21);KK(16:21,4:9),KK(16:21,16:21)]; p=[0;0;0;0;0;0;0;0;-1e5;0;0;-1e5]; u=k\p u = 1.0e-003 * 0.0223 -0.2769 -0.6728 -0.0223 -0.2769 -0.6728 [将列排成行] -0.0129 0.3108 -0.7774 0.0129 0.3108 -0.7774 [将列排成行] 附件2 一级建造师注册专业对照表（本科） 98年－现在专业名称 93－98年专业名称 93年前专业名称 土木工程 矿井建设 矿井建设 建筑工程 土建结构工程，工业与民用建筑工程，岩土工程，地下工程与隧道工程 城镇建设 城镇建设 交通土建工程 铁道工程，公路与城市道路工程，地下工程与隧道工程，桥梁工程 工业设备安装工程 工业设备安装工程 饭店工程 涉外建筑工程 土木工程 建筑学 建筑学 建筑学，风景园林，室内设计 电子信息 科学与技术 无线电物理学 无线电物理学，物理电子学，无线电波传播与天线 电子学与信息系统 电子学与信息系统，生物医学与信息系统 信息与电子科学 电子科学 与技术 电子材料与无器件 电子材料与元器件，磁性物理与器件 微电子技术 半导体物理与器件 物理电子技术 物理电子技术，电光源 光电子技术 光电子技术，红外技术，光电成像技术 物理电子和光电子技术 计算机 科学与技术 计算机及应用 计算机及应用 计算机软件 计算机软件 计算机科学教育 计算机科学教育 软件工程 计算机器件及设备 计算机科学与技术 采矿工程 采矿工程 采矿工程，露天开采，矿山工程物理 矿物加工 工程 选矿工程 选矿工程 矿物加工工程 勘察技术 与工程 水文地质与工程地质 水文地质与工程地质 应用地球化学 地球化学与勘察 应用地球物理 勘查地球物理，矿场