13杆系结构单元.pptVIP

例13-1 图示阶梯形直杆，各段长度均为，横截面积分别为3A，2A，A，材料重度为γ，弹性模量E。求结点位移和各段杆中内力。 离散化：将单元划分为3个单元，4个结点。 单元刚度矩阵： 等效结点荷载：按静力等效原则，有: 整体结构平衡方程 例13-2 变截面梁 有一变截面梁，一端固定，另一端铰支。梁长为2l，固支端的截面尽寸为b×1.6h，铰支端的截面尺寸为b×h。梁上作用均布载荷p0。求梁端的约束反力。 离散化 将梁划分成2个单元，3个结点。每个单元 长度为,截面取平均截面。 对号入座，组合整体刚度矩阵 荷载等效结点力向量 引入边界条件 3? 用式（13-77）算A1、A2、A3 和A。 13.4 坐标变换 在13.2、13.3节中，单元位移和单元力都是按单元坐标系的坐标轴分量定义的，由此建立的单元刚度矩阵属于单元坐标单元刚度矩阵。 进行系统分析时，需要把单元力按统一的结构坐标轴的分量表示出来，以便建立结点平衡方程。因此，在进行系统分析之前，必须把单元坐标系中的单元力以及单元刚度矩阵都转换到结构坐标系中去。此外，还需要把结构坐标系中的节点位移转换到单元坐标系中去，以计算结构内力。这一转换过程称为坐标变换。(一维杆和弯剪梁单元不需要坐标变换，因两种坐标系统方向一致) 结构坐标符号约定： ???——结构坐标单元位移 ?F?——结构坐标单元力 [k]——结构坐标单元刚度矩阵 1、坐标变换矩阵定义 把单元位移从结构坐标系转换到单元坐标系的变换矩阵定义为坐标变换矩阵，用符号[T]表示。有 单元坐标中的符号约定： ???e——单元坐标单元位移 ?F?e——单元坐标单元力 [k]e——单元坐标单元刚度矩阵 式（13-58）给出了结构坐标单元位移转换为单元坐标单元位移的转换式，同时是坐标变换矩阵[T]的定义式。 2、结构坐标单元力 单元力在单元位移上作的功，不因其坐标系的改变而变。则有 （13-58） 将式（13-58）代入， 对上式两端进行转置，注意到 消去???，得 即得 （13-59） 式（13-59）表明：结构坐标单元力等于单元坐标单元力前乘坐标变换矩阵的转置。 在单元坐标系中，有 3、结构坐标单元刚度矩阵 上式两端左乘[T]T， 注意到式（13-58）、（13-59），有 （13-58） （13-59） [k]——结构坐标单元刚度矩阵。 得 （13-60） 式（13-60）给出了把单元坐标单元刚度矩阵转换为结构坐标单元刚度矩阵的转换式。 引入 13.5 坐标变换矩阵 坐标变换矩阵因单元类型不同而异。 1、平面铰接杆单元（桁架元） 设OXY为结构坐标，oxy为单元坐标。?为从单元 i 端出发的任一矢量。它在结构坐标系中的分量为?X、?Y, 在单元坐标系中的分量为?x、?y。结构坐标系中的分量?X、?Y 在单元坐标x轴上投影的代数和给出?x 。同理， ?X、?Y 在单元坐标y轴上投影的代数和给出?y 。 （13-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]的计算式： （13-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、轴剪弯平面梁单元（刚架） （13-63） i j x y z 4、面外弯剪扭平面梁单元 此时，xy平面和结构坐标系XY平面仍在同一平面上，因而z轴和结构坐标系的Z轴指向相同，只须取 l ?xi ?xj ?yi ?yj Wi Wj X Y Z 在单元坐标系中，单元