清华大学结构力学求解器简介

《清华大学结构力学求解器简介》由会员分享，可在线阅读，更多相关《清华大学结构力学求解器简介（39页珍藏版）》请在金锄头文库上搜索。

1、 16-6 用求解器求极限荷载 16-6 用求解器求极限荷载 求解器可以计算一般平面结构的极限荷载并能静态或动画显示破坏机构的单向运动模 态。荷载可以是集中荷载或者均布荷载。由于极限荷载和各个杆件刚度无关，因此可以不输 入杆件刚度（当然输入也无妨） 。因此，除了按常规输入结构体系外，还要输入各杆的极限弯 矩和对解答精度控制的误差限，求解器解出的解与精确解相比其误差不超过用户指定的误差 限。由于求解器采用精确求解，通常一个杆件取为一个单元即可。对于低版本的求解器，要 求集中荷载作用在结点上，这样杆件上如有集中荷载作用，就要再分单元以使集中荷载作用 点处成为一个结点。 杆件的极限弯矩的输入同杆件的

2、刚度输入步骤一样，而且在同一个对话框中，这里不再 详述。为了输入误差限的值，依次选择菜单： “命令” 、 “其他参数” 、 “极限荷载” ，从而打开 “误差限”对话框，然后输入数值即可。如果不输入误差限，求解器将使用缺省值 0.05。以 下通过具体的例题作进一步的介绍。 例 16-6-1 试用求解器重新求解例 15-6。单元及结点编码见图 16-6-1。为比较方便，有关 数据采用单位值：m1=l，。 mkN1= u M TITLE,例 16-6-1 结点,1,0,0 结点,2,0,1 结点,3,1,1 结点,4,2,1 结点,5,2,0 单元,1,2,1,1,1,1,1,1 单元生成,3,1,

3、1,1 结点支承,1,6,0,0,0,0 结点支承,5,6,0,0,0,0 结点荷载,2,1,1,0 结点荷载,3,1,1,-90 单元材料性质,1,4,0,0,0,1,-1 单元材料性质,2,3,0,0,0,1.5,-1 极限荷载参数,0.005 END 图 16-6-1 例 16-6-1 图 解 输入的结构和命令文档如图16-6-1 所示。计算的极限荷载为3.5 kN，与精确解 完全一样。求解器中显示的破坏机构模态在 图16-6-2中给出。 例16-6-2 试用求解器求解图16-6-3中 所示的均布荷载作用下门式刚架的极限荷 载和破坏机构模态，有关数据为各杆长 ，。 m1=lmkN10=

4、u M 图 16-6-2 例 16-6-1 的破坏机构运动模态 1 16-6 用求解器求极限荷载 TITLE,例 16-6-2 结点,1,0,0 结点,2,0,1 结点,3,1,1 结点,4,1,0 单元,1,2,1,1,1,1,1,1 单元,2,3,1,1,1,1,1,1 单元,3,4,1,1,1,1,1,1 结点支承,1,6,0,0,0,0 结点支承,4,6,0,0,0,0 单元材料性质,1,3,0,0,0,10,-1 单元荷载,1,3,0.5,0,1,90 单元荷载,2,3,1,0,1,90 单元荷载,3,3,1,0,1,90 极限荷载参数,0.000005 END 图 16-6-3 例

5、 16-6-2 图 解 输入的命令文档如上图所示，计算的极限荷载为116.568 542 kN/m，破坏机构运动模 态在图16-6-4中给出。 图 16-6-4 例 16-6-2 的破坏机构运动模态 习 题 16 - 15 试用求解器求解题 16-1016-13 的极限荷载和破坏机构模态。 16 16 试用求解器求解图示刚架的极限荷载，各杆。 mkN1= u M 7.5 P qF= 题 16-16 图 2 1-6 结构力学求解器简介 1-6 结构力学求解器简介 把繁琐交给求解器，我们留下创造力！ 本书所附的结构力学求解器 （SM Solver for Windows，简称求解器）是一个方便好用

6、 的计算机辅助分析计算软件（见所附光盘） ，其求解内容涵盖了本教材中所涉及的几乎所有问 题，包括：二维平面结构（体系）的几何组成、静定、超静定、位移、内力、影响线、自由 振动、弹性稳定、极限荷载等。对所有这些问题，求解器全部采用精确算法给出精确解答。 有关求解器的安装和使用简介可以参阅软件光盘上的相关文档和联机帮助，这里不作详 细讲解。本书中各个章节中大都会出现如何使用求解器的内容，供读者酌情选学或自学。由 于本书中主要将求解器当作工具来用，因此只在目录中列出相关内容的标题（以便为读者提 供一个循序渐进的向导） ，而具体内容则放到了光盘上，同时其内在的算法和原理将在其他教 材（参见袁驷编著，

7、程序结构力学 ）中介绍。 在结构力学的学习中，求解器可以提供多种功用和帮助，有待于读者去熟悉和发现，这 里只列举几点。 1. 结构计算求解的利器 求解器首先是一个计算求解的强有效的工具。犹如计算器的功能是做数值计算，求解器 的功能是做结构计算。任意平面结构输进求解器，鼠标一点，变形图和内力图即出。 2. 体验大结构的受力性能 求解器并非专为课堂上小问题的求解而设计的。有了求解器，即便是结构力学的初学者， 也可以用求解器来方便地求解大结构的各类问题，以增强对大结构受力特性的直观感受和切 实体验。 3. 求解传统方法难以求解的问题 本书中的传统方法并非都适用于精确求解大型实际结构，例如，几何构造分

8、析的三角形 法则只适用于符合该法则的体系，超静定结构分析方法大都假定结构无轴向变形等等；而求 解器则全无这些限制。 4. 演示和作图的工具 求解器可以用静态图形显示结构简图、变形图、内力图，还可以用动画显示机构模态、 振型等动态图形。利用复制到剪贴板的功能，可以将结构简图、变形图、内力图以点阵图或 矢量图的形式粘贴到 MS Word 或 MS PowerPoint 中，并可以方便地进行再编辑。 5. 让求解器当教师答疑 几何构造分析选不出刚片吗？静定桁架的截面法剖不出截面单杆？静定组合结构分析不 知道从哪里开始下手？求解器的智能求解功能可以将求解步骤一步步地演示出来。 6. 定性分析的数值试验

9、平台 将桁架的结点都设为刚结点，看看弯曲应力是不是次应力？对于多层刚架，强梁弱柱和 强柱弱梁的受力特点有什么不同？梁的剪切变形是否小到可以忽略？这些都可以用求解器来 做具体的试验和验证。 1 2-4 在求解器中输入平面结构体系 2-4 在求解器中输入平面结构体系 本节介绍如何在求解器中输入结构体系。 在求解器中定义结构体系要满足数值化和可视化两方面的要求，另外还要有较好的人机 交互功能。如图 2-4-1 所示的几个铰支座，从计算角度上看都是等效的，但在图形显示上却不 同，因此需要补充一些有关的数据信息。另外，对于各类组合结点，在求解器中也要求能够 比较方便地定义。因此，在求解器中如何定义一个结

10、构体系是一个很基本的问题。 图2-4-1 铰支座的不同图形显示 1. 坐标系 求解器中的坐标系遵循以下的统一定义： 对于整个结构体系，规定一个整体坐标系， 用表示，如图 2-4-2 所示。整体坐标系中，),(yx x和方向（即水平和竖直方向）的线位移分别 记为u和， 均规定与整体坐标正方向相同为正； 角位移 y v 规定由x向方向转动为正 。 y 图 2-4-2 整体坐标和局部坐标 对于结构体系的每一个杆件单元，建立一个局部坐标系，用),(yx表示 。在局部坐标系 中，x和y方向的线位移分别记为u和v，相应的力记为 x F和 y F，均规定与局部坐标正方向 相同为正；角位移和相应的力矩M也规定

11、由x向y方向转动为正。 局部坐标系和整体坐标系之间要满足一定的关系。 图 2-4-2 所示为一个典型的杆件单元 e。 为单元规定一个方向，即指定一个始端 1 和一个末端 2。单元方向也可用轴向的箭头来表示： 箭头从始端指向末端。规定局部坐标的原点取在杆端 1，x轴指向杆端 2。y轴的选取应与整 体坐标转向一致，即当x轴向轴方向转y角度使得x轴与x轴同向后，轴应与yy轴同向。 按照这样的定义，角位移和力矩在局部坐标下和整体坐标下方向是一致的，所以没有必要加 以区分。 2. 虚拟刚结点 首先引入的一个概念是虚拟刚结点。 用求解器输入一个结构体系时，首先输入一些结点。这些结点被理解为是虚拟刚结点，

12、即将每一个结点看作是一个具有 3 个自由度（两个平移，一个转动）的小刚体。虚拟刚结点 是杆件之间连接的中介。在随后定义杆件单元时，不是去定义单元杆端和杆端的连接，而是 定义杆端与虚拟刚结点的连接。每个杆端（整体坐标中）有三个位移自由度，为它们各自建 立一个连接码，用 0 和 1 表示：0 表示不连接，1 表示连接。这样，若杆端的连接码为(1,1,1)， 按顺时针或逆时针表示角位移正方向容易引起混乱。如图 2-4-2 中角位移逆时针为正，但若将轴取为向下为正（yx轴方 向不变） ，则角位移变为顺时针为正。 在符号表示上，通常的规则是局部坐标系中的量比整体坐标系中的量在顶上多加一横线。 1 2-4

13、 在求解器中输入平面结构体系 则该杆端与虚拟刚结点为刚结；(1,1,0)为铰结；(1,0,0)为水平链杆连接；(0,1,1)则为定向（水 平滑动）连接，等等。 下面结合图 2-4-32-4-4 中的例子来讨论。 图 2-4-3 刚结点的连接 图 2-4-4 组合结点的连接 图 2-4-3 是刚结点的连接示例，其中图 2-4-3a 中定义了一个虚拟刚结点和杆端的连接码； 各个杆端与虚拟刚结点连接后成为图 2-4-3b 的形式，去除虚拟刚结点后的效果为图 2-4-3c 所 示的刚结点；求解器中显示的是最后的图 2-4-3c。图 2-4-4 是组合结点的连接示例，同理，无 需重复。铰结点是最常见的结

14、点之一，其连接示例在图 2-4-5 中给出。这里，共有四种连接方 式，都等效于图 2-4-5e 中的铰结点，通常采用图 2-4-5a 所示方式即可。值得一提的是，如果 将三个杆件固定住，图 2-4-5bd 中的虚拟刚结点也随之被固定不动，而图 2-4-5a 中的虚拟刚 结点仍然存在一个转动自由度，可以绕结点自由转动。这是一种结点转动机构，在求解器中 会自动将其排除不计 。结点机构实际上也潜存于经典的结构力学之中，如将一个集中力矩加 在铰结点上，便可以理解为加在了结点机构上（犹如加在可自由转动的销钉上） ，是无意义的。 图 2-4-5 铰结点的连接 综上所述，求解器中单元对话框中的“连接方式”是

15、指各杆端与虚拟刚结点的连接方式， 而不是杆件之间的连接方式。这样，各杆件通过虚拟刚结点这一中介再和其他杆件间接地连 极限分析中的结点机构将被视为破坏机构的一种。 2 2-4 在求解器中输入平面结构体系 接。这种处理的好处是可以避免结点的重复编码（如本书中矩阵位移法中所介绍的） ，同时可 以方便地构造各种复杂的组合结点。 另外，在定义位移约束时，结点处的支座约束也是首先加在虚拟刚结点上，再通过虚拟 刚结点施加给其他相关的杆端。 有了虚拟刚结点的概念，其他的具体做法都已比较显而易见，通过具体例题便很容易搞 懂。因此，以下基本上以例题来代替叙述。疑问之处，可以参见求解器的联机帮助。 3. 输入结构体

16、系 欲输入一个结构体系，首先在“编辑器”中打开一个新文件，然后输入命令。在求解器 中输入命令有两种方法： (1)利用“命令”菜单中的子菜单，打开相应的对话框，在对话框中根据提示和选项输入 命令。 (2)在文件中直接键入命令行。 第(1)种方法可以免除用户记忆命令格式，且可以利用对话框中的预览功能随时修正输入 的命令。如果用户对命令格式较熟悉，则直接键入命令更为快捷。本盘在附录中给出了所有 命令的语法规则和格式，需要时可以查阅。 修改命令也有两种方法： 1) 将光标置于要修改的命令行，单击工具栏上的“修改”按钮（或在“命令”菜单下选 “修改”子菜单） ，则求解器将自动打开相应的命令对话框，对话框中的各选项再现了该条命 令中的各个参数，用户可以查看或修改。 2)在文件中直接修改命令行。 求解器中的每一条命令都由一关键词作先导，自版本 v1.5 开始，在保留原来的英文关键 词的同时，允许使用相应的中文关键词，在功能上它们是等效的。英文关键词键入方便，中 文关键词易读易辨。求解器中默认的选择是中文关键词，如果希望用英文关键