新编力学教程课件作者穆能伶14结构的计算机分析方法.ppt

14 结构的计算机分析方法 本章内容 本章在介绍了矩阵位移法的基本原理后，还将简单介绍杆系结构的计算机分析方法及ANSYS软件的结构静力学分析。 14 结构的计算机分析方法 14.1 矩阵位移法基本方程 14.2 矩阵位移法计算平面刚架示例 14.3 ANSYS软件在结构静力学分析中的应用 由于刚架的结构刚度矩阵是奇异矩阵，无法从所对应的方程中求出结点位移。因此，解题时必须引入位移边界条件。也就是应从上面的单元刚度矩阵中删去与单元始端3的支座约束条件，即相对应的行和列 ，得 第二节 矩阵位移法计算平面刚架示例 再从结点1的变形谐调条件即 ，将单元刚度矩阵上方标注出的位移、、（亦即局部码）分别换成相应的总码 ，同时在其左侧也予以标注，于是得 1 2 3 （d） 第二节 矩阵位移法计算平面刚架示例 这种在编写程序时要进行的编号转换，叫做将单元刚度矩阵的局部码换成总码或简称换码。式（d）即为以总码表示的单元（2）的刚度矩阵。经换码后，单元刚度矩阵中的各元素所表示的即为 第二节 矩阵位移法计算平面刚架示例 同理，也可得到其他两个经换码后单元刚度矩阵的各元素。然后，使各元素在式（c）的结构刚度矩阵中一一就位，也就是所谓的“对号入座”，即得到刚架的结构刚度矩阵K。由此，写出整体结构最后的结构刚度方程即为 (d) 第二节 矩阵位移法计算平面刚架示例 由此可见，结构刚度矩阵是直接利用刚度系数集合而成，所以称为直接刚度法。实际上，根据位移边界条件和变形谐调条件在进行换码时，即引入了单元的定位向量λ（e）。所谓单元的定位向量，就是指将该单元的两端位移分量对应的总码，按顺序排列而形成的列向量。例如，单元（1）的定位向量为 第二节 矩阵位移法计算平面刚架示例 对于单元（2），由于其底端的3个位移分量大小为0，而其顶端的位移分量对应的总码为1、2、3，因此其定位向量为 显然，在中明确地指出了单元λ（2）顶端的位移分量所对应的总码及其底端的约束条件。所以，由，进而在单元刚度矩阵中删去与中的 零元素相对应的行和列，即容易得到式（d）所示的K（2）。这种在集成刚度矩阵时，采用了经边界条件处理的非自由式单元，于是便直接得到了缩减的结构刚度矩阵，故将其称为“先处理法”，这在编制一些专用程序时常常要用到。 第二节 矩阵位移法计算平面刚架示例 如果不是先考虑支座对结点的约束作用，而认为全部结点都有可能发生位移，然后据此对其全部结点位移一一进行编码并形成结构刚度矩阵。最后，再考虑支座约束条件，而对结构刚度矩阵进行修改，故将其称为“后处理法”。 （3）求结构刚度矩阵K的逆矩阵K–1。 （4）由式 未知自由结点位移。 （5）由自由结点位移计算单元的杆端力。 第二节 矩阵位移法计算平面刚架示例 例14-1 试用先处理法建立图14-7a所示刚架结构的结构刚度矩阵及刚度方程。设各杆的E、A、I都相同。 第二节 矩阵位移法计算平面刚架示例 解：（1）对刚架的结点和单元进行编号，选取图所示的结构坐标系和单元坐标系，并进行结点位移分量编号，如图14-7b所示。 (2)建立自由结点的位移列向量（图14-7c）和相应的结点荷载列向量，亦即为 第二节 矩阵位移法计算平面刚架示例 (a) (b) （3）建立相对于结构坐标系的单元刚度矩阵，列出单元定位向量λ（e），并对各单元刚度矩阵按定位向量进行换码。 单元（1）： ，因α1=0，故可由式14-4得 第二节 矩阵位移法计算平面刚架示例 1 2 3 4 5 6 第二节 矩阵位移法计算平面刚架示例 单元（2）： ，因α2=900，故由式14-4得 0 0 0 1 2 3 第二节 矩阵位移法计算平面刚架示例 在上式中，删去与定位向量中零元素相对应的行和列，得 1 2 3 单元（3）： ， ，同样可得 1 2 3 第二节 矩阵位移法计算