有限元大作业 .doc

PAGE 6 有限元原理及工程应用作业 1 对一维有限元问题 （1）推导出三点二次单元的形函数，给出相应的图形，并与两点线性单元进行比较； （2）若用三点二次单元对图示一端固定，弹性模量为E，截面面积为A的等直杆轴向拉压问题进行有限元分析，试给出总体有限元方程。 （3）给出左端点的位移。L/4 L/4 L3/4 F 3F 解：（1）三点二次单元： -10 -1 0 1 i j m 设： 解得： 所以形函数： 两点线性单元： （2）对一个单元ij 对于1号单元： 对于2号单元： 组装总体刚度矩阵： （3）考虑约束条件：=0， 解得： 3． 使用四节点四边形单元的形函数推导热传导问题的单元刚度矩阵。 四结点等参数单元的插值函数为，采用它作为温度插值函数并代入温度值T，得到 形函数为 （i=1~4） 假设局部坐标系下的单元形状、结点排列和边的定义如下图所示，节点坐标为1：（-1，-1），2：（1，-1），3：（1,1），4：（-1,1）。 因此坐标函数为 （3-1） （3-2） 根据 得到 令，，， 则 得到 （3-3） ds与或之间的变换为 当ds在2号、4号边时，=常数，=0，则；当ds在1号，3号边时，=常数，=0，则。 为便于推导，写成统一的格式为： （3-4） 令 其中为单元内任一点处的温度值，（i=1,2,3,4）为结点处的温度值。 由式（3-2），有 ， ， 将它们代入能量泛函式，经过坐标变换，有 （3-5） 为方便后面求泛函极值，这里先写出对温度的偏导数： 将上式按i=1,2,3,4展开，并写成矩阵形式，有 式中包含的XC，XA，YC，YA，，均与温度T无关，只是的函数，而A和B则是结点温度，，，的函数。 又由于 所以 i=1,2,3,4；j=1,2,3,4 类似的，可以得到 i=1,2,3,4；j=1,2,3,4 回到温度泛函式，对温度求导，有 令 则 其中， 即为单刚阵最后形式。 4.在总体坐标系中有一个任意三角形，设三个节点的坐标分别为：节点1：（x1, y1） 节点2：（x2, y2）节点3：（x3, y3） （1）试推导其总体坐标系中的形函数。 （2）证明（1）给出的形函数与用面积坐标表达的形函数等价。 解： （1）形函数的推导： 单元分析 任取一个一般性的单元，如图1所示。三个结点的编号为i, j, k。结点位移为( ui, vi )、( uj, vj ) 、( uk, vk )　单元结点位移为 假定单元内位移场u, v 是x, y的一次函数 y, v, q y, v, q i j k 图（1） x, u, p (1 (1-1) 为待定常数，在结点处应有 可解出： 其中 当 i, j, k的位置为逆时针排列时，2Δ恒正，且等于三角形单元面积的两倍。将这些结果代入（1-1）有 类似可得到 可以合并成 上式中的、、即为推导出的形函数。 证明（1）中给出的形函数与用面积坐标表达的形函数等价： 面积坐标 单元内任一点 p ： p(x,y) —— 直角坐标表示 —— 面积坐标表示 —— 面积坐标（独立分量只有二个） 由解析几何，有： 同上面一样令： 得： 可见（1）给出的形函数与用面积坐标表达的形函数等价。 5．二跨两端固定等截面梁如图示。在右跨受有向下的均布荷载q。试用有限元方法完成：（1） 简支处的转角及约束反力；（2）梁的剪力图与弯矩图。 解：A、有限元分析过程如下： 结构的离散化和编号。 划分左右两个单元，标记为①(梁单元1-2）②（梁单元2-3），节点为1、2、3，每个节点处有两个自由度，挠度和转角。因此节点总自由度. 对分布外载荷进行等效节点计算，有节点载荷如下： 以下处理均匀分布力： 由相关资料得到对应形函数： 对应形函数： 对应形函数： 对应形函数： 对应形函数： 最终得到节点载荷 2）对各个单元进行描述，梁单元的单元刚度矩阵。 梁单元1-2： 梁单元2-3： 组装总体刚度矩阵 整理有限元方程，应用边界条件： 代入方程即为：将约束自由度赋值为零，划行划列进行