[2018年最新整理]5杆系结构单元.ppt

例5-1 一维拉杆 图示阶梯形直杆，各段长度均为，横截面积分别为3A，2A，A，材料重度为γ，弹性模量E。求结点位移和各段杆中内力。 离散化：将单元划分为3个单元，4个结点。 单元刚度矩阵： 等效结点荷载：按静力等效原则，有: 整体结构平衡方程 例5-4 变截面梁 有一变截面梁，一端固定，另一端铰支。梁长为2l，固支端的截面尽寸为b×1.6h，铰支端的截面尺寸为b×h。梁上作用均布载荷p0。求梁端的约束反力。 离散化 将梁划分成2个单元，3个结点。每个单元 长度为,截面取平均截面。 对号入座，组合整体刚度矩阵 荷载等效结点力向量 引入边界条件 3? 用式（5-77）算A1、A2、A3 和A。 在单元坐标系中，有 3、结构坐标单元刚度矩阵 上式两端左乘[T]T， 注意到式（5-58）、（5-59），有 （5-58） （5-59） [k]——结构坐标单元刚度矩阵。 得 （5-60） 式（5-60）给出了把单元坐标单元刚度矩阵转换为结构坐标单元刚度矩阵的转换式。 引入 5.6 坐标变换矩阵 坐标变换矩阵因单元类型不同而异。 1、平面铰接杆单元（桁架元） 设OXY为结构坐标，oxy为单元坐标。?为从单元 i 端出发的任一矢量。它在结构坐标系中的分量为?X、?Y, 在单元坐标系中的分量为?x、?y。结构坐标系中的分量?X、?Y 在单元坐标x轴上投影的代数和给出?x 。同理， ?X、?Y 在单元坐标y轴上投影的代数和给出?y 。 （5-61） 写成矩阵形式， x y i ? ia+ab X Y ?X ?Y ?x ?y ? b a c d cd-ca 取 ——i节点在单元坐标系中的位移向量 ——i节点在结构坐标系中的位移向量 x对X、Y的方向余弦 y对X、Y的方向余弦 同理可得单元j节点在单元坐标系和结构坐标系中的位移向量： 有 组合上述结果，得平面铰接杆单元的单元坐标单元位移和结构坐标单元位移之间关系： i、j两节点间的位移变换关系互不耦合。 上式可写成 坐标变换矩阵[T]的计算式： （5-62） 2、弯剪平面梁单元 如果在连续梁中使用这类单元，通常可将单元坐标和结构坐标方向取得一致。此时，无须进行坐标变换。 ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ i j (e) X Y i j (e) x y i j x y ?2 ?3 ?5 ?6 l ?1 ?4 于是得到： X Y 由于?1 、 ?2 、?4、 ?5的性质和平面铰接杆相同，因而有相同的[T]矩阵。又因单元坐标系xy平面和结构坐标系XY平面在同一平面上， 因而单元坐标系z轴和结构坐标系的Z轴总有相同指向，所以恒有： 3、轴剪弯平面梁单元（刚架） （5-63） i j x y z 4、面外弯剪扭平面梁单元 此时，xy平面和结构坐标系XY平面仍在同一平面上，因而z轴和结构坐标系的Z轴指向相同，只须取 l ?xi ?xj ?yi ?yj Wi Wj X Y Z 在单元坐标系中，单元每个节点（如i）有3个位移分量： ?xi、?yi和wi 它的变换式和承受轴弯剪平面梁单元中向量 的变换式相同。 并且恒有 i j x y vi ?i vj ?j l ui uj 向量 因此，变换矩阵与轴弯剪平面梁单元相同。 轴弯剪平面梁单元 由此知： 5、空间杆单元 空间杆单元的每个节点有3个相互垂直的线位移分量（u、v、w）。单元自由度为6，如下图。 其中，坐标变换矩阵[T]与轴弯剪平面梁单元相同。 x y z i j 设向量?在单元坐标系和结构坐标系两个坐标系中的分量被表示 为： ?1 ?2 ?3 ?4 ?5 ?6 X Y Z x y z ? ?x、 ?y、 ?z、——向量 ? 在单元坐标轴上的分量 ?X、 ?Y、 ?Z、——向量 ? 在结构坐标轴上的分量 有 （5-65） [?]是坐标系的旋转矩阵，是单元坐标轴x、y、z在结构坐标系XYZ中的方向余弦： ?11、?12、?13——x轴在结构坐标系XYZ中的方向余弦 ?21、?22、?23——y轴在结构坐标系XYZ中的方向余弦 ?31、?32、?33——z轴在结构坐标系XYZ中的方向余弦 （5-64） 容易理解，式（5-64）可代表空间铰接杆中一个节点的节点位移坐标变换。空间杆单元有2个节点，所以坐标变换矩阵一般可表示为： （5-66） 下面讨论[?]矩阵中元素?ij(i=1、2、3，j=1、2、3)。 对于空间杆单元