计算力学课程设计模板-改 - 副本研讨.docVIP

本科课程设计 计算力学课程设计说明书 专业班级 姓 名 学 号 指导教师 完成日期  目录 1.引言 有限元法是一种高效能、常用的计算方法。有限元法在早期是以变分原理为基础发展起来的，所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类物理场中（这类场与泛函的极值问题有着紧密的联系）。它最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。在河道数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法，从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。不同的组合同样构成不同的有限元计算格式。对于权函数，伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小；在配置法中，先在计算域内选取N个配置点。令近似解在选定的N个配置点上严格满足微分方程，即在配置点上令方程余量为0。插值函数一般由不同次幂的多项式组成，但也有采用三角函数或指数函数组成的乘积表示，但最常用的多项式插值函数。有限元插值函数分为两大类，一类只要求插值多项式本身在插值点取已知值，称为拉格朗日(Lagrange)多项式插值；另一种不仅要求插值多项式本身，还要求它的导数值在插值点取已知值，称为哈密特(Hermite)多项式插值。单元坐标有笛卡尔直角坐标系和无因次自然坐标，有对称和不对称等。常采用的无因次坐标是一种局部坐标系，它的定义取决于单元的几何形状，一维看作长度比，二维看作面积比，三维看作体积比。在二维有限元中，三角形单元应用的最早，近来四边形等参元的应用也越来越广。对于二维三角形和四边形电源单元，常采用的插值函数为有Lagrange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等。 2.基本力学模型 如图所示的一个块体，在右端面上端点受集中力F作用。基于MATLAB平台，计算各个节点位移，支座反力以及单元的应力。取相关的参数E=1.5e10Pa,u=0.25,F=1.5e5N。 对该问题进行有限元分析的过程如下： 图一 (1)结构的离散化与编号 将结构离散为5个4节点四面体单元，单元编号和坐标如上图所示，连接关系见表一，节点的坐标见表二。 表一 单元连接关系 单元号 节点号 1 1 4 2 6 2 1 4 3 7 3 6 7 5 1 4 6 7 8 4 5 1 4 6 7 表二 节点的坐标 节点 节点坐标/m x y z 1 2 3 4 5 6 7 8 0 0 0 0.4 0 0 0 2.0 0 0.4 2.0 0 0 0 1.0 0.4 0 1.0 0 2.0 1.0 0.4 2.0 1.0 节点位移列阵 节点外载列阵 其中 约束的支反力列阵 其中 总的节点载荷列阵 （2）计算各单元的刚度矩阵(以国际标准单位) 首先在MATLAB环境下，输入弹性模量E、泊松比NU，然后针对单元1和单元2，分别5次调用函数Tetrahedron3D4Node_Stiffness，就可以得到单元的刚度矩阵k1(6×6) ～ k5(6×6)。 E=2.5e10; NU=0.25; k1 = Tetrahedron3D4Node_Stiffness(E,NU,0,0,0,0.4,2.0,0,0.4,0,0,0.4,0,1.0); k2 = Tetrahedron3D4Node_Stiffness(E,NU,0,0,0,0.4,2.0,0,0,2.0,0,0,2.0,1.0); k3 = Tetrahedron3D4Node_Stiffness(E,NU,0.4,0,1.0,0,2.0,1.0,0,0,1.0,0,0,0); k4=Tetrahedron3D4Node