理论力学求解器.ppt

理论力学求解器 太原理工大学 静力学 平面静定桁架的计算 平面一般力系的计算 平面静定桁架程序设计 基本原理 力学模型与数学模型 数据准备 程序使用 算例 计算 基本原理（桁架） 用节点法求解具有NJ个节点、NG根杆件及3个外约束反力平面静定桁架。因每个节点都受平面汇交力系作用，该节点有两个独立的平衡方程，整个桁架共有2NJ个独立的平衡方程，联立求解这些方程，即可解得全部轴力及外约束反力。 本程序设计只讨论节点数NJ及杆件数NG满足式： 2NJ=NG+3 的情况。 力学模型与数学模型1（桁架） 研究桁架中的任一节点S，称之为主节点，通过杆件k与之相连的节点j称为从节点(图示)。若在主节点S上连有m根杆件，作用l个主动力和n个外约束反力，并设各杆的轴力均为拉力。由平衡方程得 力学模型与数学模型2（桁架2） 上式还可改写为： 力学模型与数学模型3（桁架） Ri——作用在主节点S上的编号为i的外约束反力Ri的大 小； ——外约束反力Ri的方向角； Fi——作用在主节点S上的编号为i的主动力Fi的大小； ——主动力Fi的方向角。 求解上述方程组，即可求得全部轴力和外约束反力。 数据准备（桁架） 基本信息 主动力信息 节点杆件信息 外约束反力信息 数据准备基本信息（桁架） 基本信息： N 主动力个数； NJ 节点个数； NG 杆件个数。 数据准备主动力信息（桁架） 主动力信息包括各力的大小，方向及作用节点的编号。这些信息可用两个数组F(N，2)，和IIF(N)来描述。 数据准备节点杆件信息（桁架） 节点杆件信息包括节点相关信息、辅助信息及节坐标信息。它们可用四个数组，JJ(NJ，NJ)、 JD(NJ—1)、X(NJ)和Y(N了)来描述。 (1)节点相关数组JJ(NJ，NJ) 数组JJ代表一个方阵，用来描述各节点之间的约束关系。其元素按下式计算： 0 (当i、j节点不相关时) i、j节点之间所连杆件的编号k (当i、j节点相关时) 数组JJ是一个主对角线元素均为零的对称矩阵，而且大部分元素都是零。 数据准备节点杆件信息（桁架） (2)辅助数组JD(NJ一1)其元素按下式计算： JD(i)＝JJ阵之上三角阵第i行非零元素的个数 i＝1，2‥‥‥, (NJ一1) (3)节点坐标数组X(NJ)、Y(NJ) 分别用来存放各节点的坐标xi和yi。 数据准备外约束反力信息（桁架） 外约束反力信息包括外约束反力Ri的编号、其作用节点的编号S及该反力的方向角 。它们可用两个数组IR(3，2)和R(3)来描述。 其中， 按右手法则量取，单位是度；外约束反力Ri 的编号应排在各轴力的编号之后，即三个外约束反力的编号分别为NG+i (i＝l，2，3)。 程序使用1（桁架） 在程序中所有输入数据一律按数据文件的形式输入微机。输入数据的顺序如下： 1) 输入 N、NJ、NG。 2) N 次循环输入F(i，1)、F(i，2) 和IIF(I)，即输入各力的大小Fi、方向角 及作用节点的编号S。 3) (N一1) 次循环输入辅助列阵{J D}各元素。 4) 用双重循环输入[JJ]上三角阵中的各非零元素。 5) NJ 次循环输入节点坐标X(i)、y(i)。 6) 三次循环输入IR(i，1)、IR(i，2) 和 R(i)，即输入外约束反力Ri，的编号NC十i 作用节点的编号S 和它的方向角 。 计算结果同时以屏幕显示和输出数据文件两种方式输出。 程序使用2（桁架） 程序运行过程如下 1) 运行程序命令。 2) 系统提示键入输入／输出数据文件名。 3) 键入文件名，文件名必须用单引号括住，二文件名之间用逗号分隔。 算例（桁架） 例1 例2 例1 例1 求如图所示桁架支座反力及各杆轴力。 解：取坐标系xy，给杆件编号后，再接着给外约束反力编号，按图中圆圈內所标数字顺序给节点编号。则有： 例1 例1 外约束反力信息： 例1 输入数据文件： 输入数据文件： 2，6，9 R(1)=50.000 50，270，1