力矩分配法的方程解法
中国农业大学,北京100083
A SOlVING EQUATION FOR MOMENT DISTRIBUTION METHOD
China Agricultural University, Beijing 100083, China
Received: 2021-07-7 Revised: 2021-11-4
作者简介 About authors
在传统力矩分配法的基础上,将渐进的计算过程转化为方程组的求解,可以更好地理解力矩分配法的分析思路及计算过程,同时提高计算效率和精度。该方法所得到的方程与矩阵位移法所得方程形式相同,但未知量完全不同,方程的推导过程更简洁且易于理解,在教学中可以启发学生对矩阵位移法分析思想的思考。
关键词:
This paper transformed the gradual calculation process in the traditional moment distribution method into the solution of equations. The method can improve the calculation efficiency and accuracy. The equations obtained by this method are similar to the equations obtained by the matrix displacement method, but the unknowns are completely different. The derivation process of the equation is more concise and easier to understand. It can inspire students to think about the analysis of the matrix displacement method in teaching.
Keywords:
本文引用格式
朱晓江, 张琳楠, 秦太验.
ZHU Xiaojiang, ZHANG Linnan, QIN Taiyan.
本文提出了力矩分配法的方程解法,该方法将渐进的计算过程转化为方程组的求解,方程组的系数矩阵为三对角阵,且主对角线元素始终为"1"。该方程组求解简单,适于少结点问题的手算;对于多结点问题,采用相同顺序的结点编号时,应用该方法可以得到与矩阵位移法形式相同的方程,方程未知量为各分配点的分配力矩代数和,方程的推导不需要单元刚度矩阵的集结。力矩分配法的计算精度随计算轮次的增加而提高,最后收敛于精确解,本文方法将渐进计算过程转化为方程组的求解,因而可以直接得到杆端弯矩的精确解。
力矩分配法是为了避免解算力法和位移法典型方程而发展起来的、简捷的、适于手算的渐进计算方法,但力矩分配法出现了迭代次数和精度的新问题,本文提出的方法解决了这两个新问题。虽然本文方法又回到方程的求解,但是一方面该方法所需解算的方程较力法和位移法典型方程更为简单,另一方面,本文方法其实正是力矩分配法向矩阵位移法的过渡,该方法让我们进一步理解,从想摆脱方程求解,到回到方程求解的必要性,在教学中指出这一必要性,有助于学生乐于接受更繁琐的矩阵位移法。此外,对于多结点问题,相比于矩阵位移法,本文方法的方程建立过程更简洁,有助于启发学生的编程思维及相应的计算练习。
1 力矩分配法的方程解法原理
传统的力矩分配法只适用于求解无侧移结构,以李廉锟[1]书中例9-2含有3个分配点的连续梁为例,说明力矩分配法方程解法的原理。
观察图1计算表格,每个椭圆圈内所有数的代数和为"0";$B$左侧截面只有分配力矩,$B$右侧截面则既有分配力矩又有传递力矩,$D$结点与$B$结点情况相反,$C$为中间的分配点,两侧均有分配力矩和传递力矩。由上述观察可知,$B$点两侧所有传递力矩和分配力矩的代数和等于$B$两侧截面固端力矩代数和的相反数,$C$和$D$两结点处也有相同的结论,即矩形框内所有数值的代数和等于各结点两侧截面固端力矩代数和。
图1
设$B$,$C$,$D$三个分配点处各次分配力矩的代数和分别为$M_{B}$,$M_{C}$和$M_{D}$,同时考虑各结点所连杆端分配系数代数和为"1", 则有
该方程组主对角线元素均为"1",求解方程组可得三个分配点处各次分配力矩的代数和${M}_{{B}}$,${M}_{{C}}$和${M}_{{D}}$。分配点杆端弯矩可按式(2)计算
非分配点$A$和$E$的杆端弯矩数值,等于求解方程组得到的相邻分配点分配力矩代数和乘以相应传递系数,再叠加截面相应的固端弯矩数值。
2 力矩分配法的方程解法算例
为方便对比,仍以李廉锟[1]书中例9-2为例,用力矩分配法的方程解法计算各杆端弯矩。
由力矩分配法的方程解法计算过程可知,该方法将传统的渐进计算过程转化为三对角矩阵方程组的求解,简化了计算步骤及计算工作量,并得到精确解。
此例题用位移法求解时的典型方程为
对比方程(3)与图2计算表格中的方程可知,两方程形式一致,但推导过程及未知量完全不同,本文解法给出的方程主对角线元素均为"1",求解更简单。
图2
3 力矩分配法的方程解法一般方程及其解
对有$n$个分配点的多跨梁或单层无侧移刚架,力矩分配法的方程解法即求解式(4)中系数矩阵为三对角阵,且主对角线元素均为"1"的方程组。
其中${C}_{{ij}}$,${\mu}_{{ij}}$及${M}_{{i}}^{{f}}$均为题目的已知条件,求解方程得到$M_{i}$ ($i=1$,2,$\cdots$,$n$),各分配点杆端弯矩表达式为
非分配点杆端弯矩等于固端弯矩加上相邻分配点传来的传递力矩。
对于多分配点连续梁问题及如图3(a)所示的只有一层结点的无侧移刚架问题,上述方法将一般力矩分配法的渐进计算过程转化为三对角矩阵方程组的求解,在不改变力矩分配法的分析思路,保留杆端弯矩作为分析对象的同时,极大简化了计算步骤及计算工作量,并得到精确解。
图3
对于图3(b)所示的两层无侧移刚架,或其他的多层无侧移刚架,应用本文方法仍可快速写出需要求解的方程组,此时的系数矩阵主对角元素仍为"1",带宽取决于各结点的编号顺序。关于本文方法求解多层无侧移刚架及与矩阵位移法的对比还可以进行进一步的讨论。
4 结论
本文提出了力矩分配法的方程解法,该方法适用于所有可用力矩分配法求解的无侧移结构。对于较多分配点的连续梁问题及各分配点连接较多杆件的单层无侧移刚架问题,该方法可简化计算步骤,提高计算效率,并得到问题的精确解。
在教学中,本文提出的方法没有增加新的参数,易于学生的理解、接受和应用;力矩分配法与矩阵位移法都是在位移法的基础上发展起来的,力矩分配法直接以杆端弯矩作为分析对象,采用本文方法可以将问题转化为方程组的求解,由此可以引导学生思考,是否可以在将结点位移作为基本未知量的求解过程中,同样引入适于计算机处理的方程解法,而这一思想正是矩阵位移法基本思路,因此本文方法有助于拓展学生的解题思路,引发学生的深入思考,提高学习积极性。
参考文献
含3个分配点结构的弯矩分配公式法精确解
,
Moment distribution formula method to seek exact solution of structures with three distribution points
一次性分配的力矩分配法
,
Moment distribution method of one-off distribution
/
〈 | 〉 |