工程结构力学 教学课件 ppt 作者 程选生 第八章 杆系结构的计算机分析 .ppt

按单元定位向量 ①＝（1 2 3 0 0 0）T，将P①中的元素在P中定位，得 P= （6）解基本方程 10-3× = 求得 = （7）求各杆杆端力 单元①：先求F①,然后求 ①， F①=k① ①+ ①= k① ① + TT ① = k ① ① - P ① =10-3× + = ③ = TF ③= （8）根据杆端力绘制内力图,如图8-20示。 FN图(单位 kN) M图(单位 kN.m) FS图(单位 kN) 图 8-20 第七节 桁架的整体分析 如图8-21，局部坐标系下桁架单元的刚度方程已在第二节中给出，其矩阵形式为 实际上，上式是一般单元刚度方程式（8-4）的特殊形式，即将式（8-4）中删去第2、3、5、6行和列后自动得出。 对于斜杆单元，其轴力和轴向位移在整体坐标系中将有沿x轴和y轴的两个分量。因此，如图8-22所示，整体坐标系中杆端力向量和杆端位移向量为 图 8-22 ， 由于第二节推导的一般单元的转换关系具有一般性，而桁架单元是一般单元的特殊情况，只需根据桁架单元的特点，把第二节的坐标转换矩阵T稍加修改，则第二节的其它转换公式均可应用。 为了便于利用以前的坐标转换关系，我们将局部坐标系下的单元刚度方程扩大为四阶形式： （8-46） 式（8-12）是一般单元的杆端力转换式。对于桁架单元，由于不存在杆端弯矩，所以删去弯矩相应的行和列，便得到桁架单元的坐标转换矩阵T如下： 表8-1 两种编码的对应关系 ②=（2 3）T ② ①=（1 2）T ① 对应关系 单元 单元定位向量 （1） 1 （2） 2 （1） 2 （2） 3 局码 总码 （3）注意单元刚度矩阵 和单元贡献矩阵Ke中元素排列方式 在 中，元素按局码排列，或者说，元素按局码对号入座；而在Ke中，元素按总码排列，或者说，元素按总码对号入座 实际上，由单元刚度矩阵 得出单元贡献矩阵Ke的过程，就是换码重排座的过程，具体表8-2 表8-2 元素的换码重排座 改排在 行 列 原排在（i）行（j）列的元素 重排座 （i） （j） 换成新行码 新列码 元素的原行码（i） 新列码（j） 换码 在单元贡献矩阵Ke中 在单元刚度矩阵 中 现按表8-1、2的作法，由k①、k②求K①、K②，详见表8-3： ② ① Ke 单元 表8-3 由k①、k②求K①、K② 由表8-3可以看出，由 求Ke的问题实质上就是 中的元素在Ke中如何定位的问题。定位的规则是 3、单元集成法的实施方案 上面我们所讲的单元集成法之所以分为两步，是为了理解方便。实际工作中，为了使计算程序更为简单，我们经常合二为一，采取边定位边累加的办法，由 直接形成K。 详细地说，按单元集成法，形成K的过程，就是依次将每个 的元素在K中按 定位并累加的过程，其步骤如下： （1）先将K置零，这时K=0； （2）将k①的元素在K中按①定位，这时K=K①； （3）将k②的元素在K中按②定位并进行累加，这时K=K①+K②； （4）按此作法，对所有单元循环一遍，最后得到K= K① 仍以图8-7a所示的连续梁为例，来说明K的集成过程： 将k①集成后，得到阶段性结果： 在此基础上再将k②集成，得到最终结果： 可以看出，由单元集成法所得结果与传统位移法式（8-26）相同。 4、整体刚度矩阵的性质 （1）整体刚度系数的意义 K中的元素Kij称为整体刚度系数。它表示第j个结点位移分量 （其它结点位移分量为零）时所产生的第i个结点力Fi； （2）K是对称矩阵； （3）按本节方法计算连续梁时，K是可逆矩阵。 ② ① 图8-12 如图8-12，由结点力F推算结点位移 时，结构仍是一个几何不变体系，因此，当F指定为任意值时，均可求得 的唯一解，故知K-1是存在的。 （4）K是稀疏矩阵和带状矩阵 对于图8-13的n跨连续梁，根据单元集成法，不难得出其整体刚度方程： n ② ① 图8-13 （8-39） 可以看出，K有许多零元素，故称为稀疏矩阵。又非零元素都分布在以主对角线为中线的倾斜带状区域内，故K为带状矩阵。 例题 8-2 试求图8-14a所示连续梁的整体刚度矩阵K。 ③ ② ① ② ③ ① 图8-14 解：（1）结点位移分量总码 由于支座D的结点位移分量为零，所以该连续