第4章 杆梁结构的有限元分析原理.ppt

4.4 梁单元及其坐标变换 一般平面梁单元的描述 纯弯梁单元 i v j v ① i j e i ? e j ? l 由于单元有四个位移分 量，可设梁单元的位移 模式 v ( x ) 为包含 4 个待定 常数的三次多项式： 2 3 0 1 2 3 2 1 2 3 ( ) ( ) 2 3 v x a a x a x a x x a a x a x ? ? ? ? ? ? ? ? 由材料力学知 , 各截面的转角 : x v ? ? ? ? 根据边界条件可以确定待定系数，将其进一步回代，可以得到用节点 位移表示的梁单元位移。 ? ? ( ) v x N q ? ? ? ? ? ? ? 1 2 3 4 ( ) i i j j v v x N N N N N q v ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 式中 2 3 2 3 2 3 2 3 1 2 3 4 2 3 2 2 3 2 1 3 2 , 2 , 3 2 , x x x x x x x x N N x N N l l l l l l l l ? ? ? ? ? ? ? ? ? ? ? 4.3 杆单元及其坐标变换 - 局部坐标 i j x y ? 由于杆单元只有两个节点位 移，故可以设杆单元的位移 模式为之包含两个待定常数 的形式 u ( x )= a 1 + a 2 x 根据有限元法的基本思路，将弹性体离散成有限个单元体的组合，以 结点的位移作为未知量。弹性体内实际的位移分布可以用单元内的位 移分布函数来分块近似地表示。在单元内的位移变化可以假定一个函 数来表示，这个函数称为 单元位移函数 、或 单元位移模式 。 ? 回代得 写成矩阵形式为 1 2 ( ) 1 j i i i j iu i ju j u x a a x u u u x l x x u u l l N u N u ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? i iu ju j u u N N u ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 其中 N i ， N j 是形函数。 ? 根据位移条件有 u (0)= u 0 ， u ( l )= u l ，从而得 ? ? 1 2 , 1, 2 j i i u u a u a i j l ? ? ? ? ? ? 根据几何方程得 ? 根据物理方程得 ? 从而，根据单元分析结果，进行整体分析，求解整体方程组， 进行结果分析 ? ? 2 d 1 d x j i u a u u x l ? ? ? ? ? ? ? 1 1 1 i i iu ju j j u u N N u u l ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? d d x j i u E E u u x l ? ? ? ? ? ? 1 1 i i iu ju j j u u E E N N u u l ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 4.3.2 杆单元的坐标变换 i j x y i j x y 规定： 杆端位移和杆端力取在截面形心上，符号以与单元系坐标正向相 同为正，相反为负。下面讨论整体坐标系下与局部坐标系下的转换关系 式。整体坐标系单元杆端位移和杆端力仍定义在截面形心上，符号以与 坐标正向同向为正反之为负。 局部坐标系 整体坐标系 4.3.2 杆单元的坐标变换 - 平面问题 其中是一个单位正交矩阵，单位正交矩阵的逆即等于其转置。 o R A x y ? x y B u u v 从上图可以得出，整体坐标系逆针旋转 α 角后与单元系相重合。 cos sin cos sin cos sin cos sin i i i j j j i i i j j j u u v u u v X X Y X X Y ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 写成矩阵形式为 e e e e q Tq P TP ? ? 0 0 0 0 c s T c s ? ? ? ? ? ? ? 由于单元的势能是一个标量 ( 能量 ) ，不会因坐标系的不同而改变，因此，可 将节点位移的坐标变换关系代入原来基于局部坐标系的势能表达式中， 整体坐标系下的刚度方程 根据 e e e e q Tq P TP ? ? 得 其中 2 2 2 2 2 2 2 2 - - - - - - - - e e e eT e e e c cs c cs cs s cs s E A k T k T l c cs c cs cs s cs s ? ? ? ?