6 平面问题总刚度矩阵的组集.docVIP

平面问题总刚度矩阵的组集 单元刚度矩阵 → 扩展 → 叠加 三角形单元的单元刚度矩阵——单元结点力与单元结点位移之间的关系 将单元刚度矩阵按单元结点号分块、命名 9个子块！ 例： 设系统的结点数为6，e为系统中的一个三角形单元，三个定点的编码分别为1，5，6 扩展：将各子块按下标放在总刚度矩阵对应的位置 叠加扩展后的单元刚度矩阵得原始总刚度矩阵 原始总刚度矩阵：解除所有约束的结构的总刚度矩阵 引入边界条件 69页例题 原始方程 ——12行、12列 (1) 手算：划去方程组的第1、3、7、8、10、12个方程以及总刚余下的行中的第1、3、7、8、10、12列元素，得引入边界条件后的总刚度矩阵 [K]——6行、6列 （2）机解：置大数法 将位移为零的自由度所对应的总刚度矩阵的对角线元素乘以一个大数，载荷对应项置零。 例如，上例中,对应的自由度的序号为3，则将原始总刚中的乘以一个大数，比如，并将赋零（原为），这时方程组的第三个方程成为： 注意，方程左边的各项与第三项相比均可忽略，即上述方程近似等价于： 而不为零，所以这时方程所表示的意义就是 置大数后各矩阵的改变： [K]的维数不变，与受约束自由度对应的对角线元素被改变； {δ}没有改变； {P}的维数不变，与受约束自由度对应的元素均为零。 等效结点载荷—非结点载荷的处理 静力等效原则：转换后的结点载荷与原载荷在任意虚位移上的虚功相等 1．集中力的等效： 例： 图示三角形单元ijk的边界jk上作用有沿x方向的集中力P，力的作用点的坐标为（0.5,0.75），求力P的等效结点载荷 解： (1) 求 (2) 求 解得 (3) 求 解得 形函数的几何意义 在下图中，Δijk位于xoy平面内，ii’, jj’, kk’ 分别垂直于xoy平面 Δi’ij和Δi’ik为直角三角形 Δj’ji和Δj’jk为直角三角形 Δk’kj和Δk’ki为直角三角形 上例：三角形单元ijk的边界jk上作用有沿x方向的集中力P，力的作用点的坐标为（0.5,0.75），计算力P的等效节点载荷 x = 1- 0.75 = 0.25=Nj(0.5,0.75) x = 1- 0.25 = 0.75=Nk(0.5,0.75) 2．分布力的等效（面力、体力） 如果单元上受的载荷是分布载荷{q(x,y)},作用在面积微元dxdy上的力为{q(x,y)}dxdy,他的等效结点载荷为，即 3．组集结点载荷向量 在上例中，如果系统所受载荷如下图所示， 则结点3的等效结点载荷应为： 第i个结点的载荷： 作业 1．证明Δ单元的形函数，，满足 2．已知矩形单元的四个顶点坐标分别为: i(-1,-1), j(1,-1), k(1,1), m(-1,1)。求该单元的4个形函数，，，。（二次项取xy） 3. 求图示等腰直角三角形单元的分布载荷和重力Q（作用在重心）的等效结点载荷。 4．求图示系统4个结点的等效载荷，集中力Q均作用在三角形的重心。