有限元大作业—蔡进、张宗元、张超、母东、陈仡湑.docVIP

有限元基本理论及工程应用大作业 2012年5月14日 PAGE PAGE 14 指导教师：韩省亮硕1005班 指导教师：韩省亮 硕1005班： 蔡 进 3111003006 母 东 3111003010 张宗元 3111003024 张 超 3111003035 机械硕81班： 陈仡湑2012-5-14 有限元基本理论及工程应用大作业 有限元基本理论及工程应用大作业  1.对一维有限元问题：（1）推导出三点二次单元的形函数，给出相应的图形，并与两点线性单元进行比较；（2）若用三点二次单元对图示一端固定，弹性模量为E，截面面积为A的等直杆轴向拉压问题进行有限元分析，试给出总体有限元方程。（3）给出左端点的位移。 L/4 L/4 L3/4 F 3F 解：（1）三点二次单元、两点线性单元的母体单元如下图所示： 可以推导出各结点的形函数为： 相应形函数的图形如下图所示： 两点线性单元的母体单元如下图所示： 可以推导出各结点的形函数为： 相应形函数的图形如下图所示: （2）如下图建立坐标系，并在不考虑约束条件情况下取等直杆作为一个三点二次单元： 可知结点的坐标为 。 整个结构的结点自由度为 可以构造各结点形函数为： 结合题意如下图所示，对单元进行加载： 可以求得边界力的等效结点力为： 由形函数可以求出几何矩阵[B]： 因此，可以求得单元的刚度矩阵为： 引入约束条件后：因此有限元方程中的表达式如下： 为外载荷在结点的等效结点力，所以： 因此，总体有限元方程为： （3）由（2）得出的有限元方程可以得到如下方程组： 从而可以得到左端点的位移为： 2. 热传导问题的有限元分析。（1）推导对应二维热传导的“能量泛函”（2）用四节点线性形函数表示该“能量泛函”并写出其对应的一般的矩阵形式。 解：（1）二维热传导的能量泛函可推导如下： 易知二维温度场相应的微分方程和边界条件为 ? 能量泛函Ue的变分为（对ε求导后令ε=0 δ = 通过化简可得到： δ 泛函取极值的必要条件是δU - 再有变分法基本预备定理便得， ? 由题意满足边界条件。即能量泛函为 U （2）由题意：四节点等参元的插值函数为 ε= 代入温度值T: T= 其中形状函数为 N 坐标函数为 ξ= x= 通过坐标变换，令 A=?T?ξ ， XC=?x?ξ= YC=?y?ξ= 将他们代入泛函Ge式中有 G 由Ge ? 其中i=1、2、3、4，写成矩阵的形式 ? kk 该式中i=1，2，3，4；j=1，2，3，4。 同理可以得到 ? 其中 zzij= 由 ? 可得 ? 令 [KK] 可得 ? 其中[KK]4×4 3. 使用四节点四边形单元的形函数推导热传导问题的单元刚度矩阵 解：假设四节点四边形单元如下： 温度函数由下式给出： （3.1） 式中，，和是节点温度，形函数参考有限元讲义上的四节点四边形等参数单元，即： （3.2） 对，和有相似的表达式，其中，和由下列方程得到： （3.3） 定义梯度矩阵的方法与定义应力分析问题中的应变矩阵相似： （3.4） 把方程（3.1）带入方程（3.4）有： （3.5） 以紧凑的矩阵表达式表达梯度矩阵与应力分析问题的应变矩阵相似，由下式给出： （3.6） 式中矩阵是把方程（3.2）的三个方程带入方程（3.5）的右边的矩形矩阵中得到的： （3.7） 热流/温度梯度关系为： （3.8） 式中材料性能矩阵为： （3.9） 热传导问题的能量泛函如下：（由参考文献：有限元方法基础教程（第三版），[美]Daryl L. Logan著） （3.10） 式中 ： （3.11） 式中，分别指定热流（指向表面时为正）和对流损失的表面。正如方程（3.11）所示，不能在同意平面上指定和，因为它们不可能同时发生在同一表面上。 把方程（3.6）和方程（3.9）代入方程（3.11），并应用方程（3.10），可以把写为如下的矩阵形式： （3.12） 把方程（3.6）代入方程（3.12），并注意到节点温度与一般的坐标和无关，因此可放到积分外，于是得到： （3.13） 在方程（3.13）中，极小很容易进行，