5第5讲空间杆件结构的有限元法.pptVIP

第5讲 空间杆件结构的有限元法 第一节 局部坐标系下的单元分析 第二节 空间单元坐标变换 第三节空间刚架分析举例 * * 式（2-1） 设i、j、k三点在整体坐标系xyz中的坐标分别为（xi、yi、zi）、(xj、yj、zj)、(xk、yk、zk)，那么如何根据这三个点的坐标值来确定坐标系的关系矩阵t中的九个元素呢？t中的第一行元素较容易确定。如图2-4可得 * 局部坐标系下的单元分析 图2-1 所示为空间刚架中的任一 杆件单元。选取局部坐标系时，去 形心轴为轴，横截面的主轴分 别为坐标系的轴和轴。 轴的方向按右手定则确定。这样， 单元在 平面内的位移与 平面内的位移是彼此独立的。设杆截面面积为A，在平面内的抗弯刚度为，线刚度；在平面内的抗弯刚度为，线刚度；杆件的抗扭刚度为。 空间刚架单元的两端分别与结点i和j相联结。每一个结点有六个结点位移分量和六个结点力分量。在局部坐标系下空间杆件的杆端位移列阵和杆端力列阵分别为 其中为轴向位移，为横向位移，为杆件的扭转角，分别为绕轴和轴弯曲时的转角；为杆件单元的轴力，分别为沿轴和轴作用的剪力，为作用在杆端的力偶矩。这里力偶矩和角位移的指向按照右手定则用双箭头表示；力和线位移的指向用单箭头表示。图2-1中所示的杆端力和杆端位移为正方向。 与平面单元的推导方法一样，首先求出当杆端位移中的一个分量为1，而其余分量均为零时的杆端力。图2-2所示为当单元的i端发生单位位移时，杆端力与杆端位移之间的关系。图中未绘出的杆端力和杆端位移分量，在该情况下数值为零。 依同样方法可以确定当单元j端发生单位位移时，杆端力与杆端位移之间的关系。 当单元的杆端位移分量为任意值时，可以写出空间单元刚度方程，以矩阵表示为 式（2-1）可以简写为 （2-2） 其中单元刚度矩阵为 （ （2-3） 式（2-3）为局部坐标系中的空间单元刚度矩阵。它是12阶方阵，其性质也与平面结构 的相同。 空间单元坐标变换 将局部坐标系下的单元刚度矩阵转换为整体坐标系下的单元刚度矩阵，是通过坐标转换矩阵完成的。首先考虑单元在端点i的三个杆端力分量，在局部坐标系中，它们是；在整体坐标系xyz中，是与之间的关系。设轴与x、y、z轴的夹角分别为（图2-3），则轴在xyz坐标系中的方向余弦为 将杆端力在轴上投影， 可求得杆端力 同理可得 综合上三式 (2-4) 这就是在端点i由整体坐标系中的杆端力推算局部坐标系中杆端力的转换关系式。其中两坐标系的转换矩阵（简称“关系矩阵”）为 （2-5） 参照上述方法，同样可以推出以表示，以表示，以表示的表达式，其转换矩阵也是t。 综合以上分析，整体坐标系中的单元杆端力分量列阵F与局部坐标系中单元杆端力分量列阵之间的关系，可用下时表达 （2-6） 同理，可导出整体坐标系与局部坐标系杆端位移之间的转换关系 （2-7） 在以上两式中 （2-8） 称为“单元坐标转换矩阵”；它是12×12阶矩阵，是一个正交矩阵，故有 （2-9） 在平面结构中，确定了单元的两个结点i和j的坐标，就确定了杆件的位置。在空间结构中，仅仅确定两个端点的坐标还不能完全确定刚架杆件在空间的位置，因为相同的ij杆，其截面形心主轴仍可由不同的方向。为确定刚架杆件在空间的确切位置，还需要在杆轴线外在取一点k，以确定其形心主轴的方向。 取结构的整体坐标系为xyz，单 元局部坐标系为，为杆件截 面形心主轴之一，如图2-4所示。单 元的位置由i、j、k三个点的坐标决 定。这里i为单元起始结点号，j为 单元终点号，由i、j两点可确定 的方向。K点在单元所在的平面 内，但又不在轴上，如果刚架上找 不到合适的点，可用一个假想的点代替。 (2-10) 其中l为杆长，可按下式求得 （2-11） 设i、j、k分别为三个坐标轴方向的单位矢量，轴矢量x可表示为 （2-12） 因为轴的矢量z与平面ijk垂直，所以有 （2-13） 为后面的运算方便，可设 则有 轴的方向余弦为 为后面的运算方便，可设 则有