论文部分内容阅读
本文研究结构稳定极限承载力分析的简便方法。根据物理概念,提出了系统(结构)静、动力平衡状态稳定性分析的力素增量方法。建立了结构单元弹塑性增量平衡方程及结构单元弹塑性增量割线刚度矩阵,推导了五种常见结构单元的弹塑性增量割线刚度矩阵的显式表达式,避免了数值积分运算。提出了一致显式增量割线刚度直接迭代算法,提高迭代效率。根据虚位移原理,采用Updated-Lagrange列式及本文提出的力素增量方法建立了结构单元大变位分析的虚功增量平衡方程,揭示了该式左边各项的物理意义。以平面梁单元为例,得到了反映单元大变形抗力增量与荷载增量相互平衡的代数式,将复杂的积分方程改造成为代数方程,简化了计算。根据以下物理事实:单元大变形过程可分解为自然变形阶段和整体刚性运动阶段,单元只在自然变形过程中才会改变节点抗力大小,而刚体运动不会改变单元抗力的大小,只会改变其方向,导出了结构单元大变位分析的增量割线刚度方法并推导了四种结构单元大变位分析的增量割线刚度矩阵。主要的研究内容和取得的成果如下:(1)提出了结构静力平衡状态稳定性分析的力素增量评判准则:当抗力增量大于荷载增量时,平衡状态稳定;小于时,平衡状态不稳定;等于时,结构处于临界失稳状态。建立了系统运动稳定性分析的位移变分法,解决了非线性系统稳定性简便分析以及系统稳定性分布区域的确定问题。(2)推演了杆单元、抛物线索单元、平面梁单元、空间梁单元、板壳单元等的弹塑性刚度矩阵计算表达式,无需数值积分运算,简化了计算过程,提高运算效率。(3)提出了一致显式增量割线刚度直接迭代算法,上一迭代步“校正”阶段形成的结构单元增量割线刚度矩阵可以用于本次迭代步“预测”阶段计算结构单元增量位移,迭代过程简单,提高求解效率。(4)以平面梁单元为例,揭示了结构单元大变形三维虚功增量平衡方程中虚应变能增量的物理意义:分别表示单元抵抗自然变形的抗力虚功,单元刚体变位引起的单元自然变形抗力增量虚功以及单元初始节点力增量虚功,从而建立了单元大变形力素增量平衡方程。(5)建立了结构单元几何非线性分析的增量割线刚度列式、四种常见结构单元增量割线刚度计算表达式,以及提出了这些单元大变形分析的增量割线刚度方法。(6)提出了基于增量割线刚度的柱面弧长算法,求解结构弹塑性大变形非线性平衡方程,追踪结构从开始加载直到顺利跨越极值点全过程的平衡路径。该算法可以有效地避免传统柱面弧长法中由于荷载增量因子选取不当引起求解路径回溯并最终导致迭代过程发散的问题,求解过程简单,可以收敛于正确解。(7)基于本文理论研究成果,编制了结构稳定极限承载力分析的计算程序。分析了文献中选取的11例结构非线性响应计算的经典考题,得到的结果均与解析解或者与其它学者提出的数值解非常吻合,计算了三汊矶湘江大桥在全桥满布荷载工况下的稳定极限承载力,从而验证了本文分析方法的正确性和有效性。