积分算子法-单元荷载教案.pdf

该文发表于：《工程力学》，2002 年 12 月, 19 卷第6 期，126~130 非线性杆单元荷载矩阵积分算子直接集成算法 谢靖中 李国强 ( 同济大学建筑工程系, 上海200092) 摘要：在文献[1]的基础上，本文将积分算子直接集成单元刚度矩阵的方法扩展应用于非线性杆单元的 荷载矩阵集成。该方法将分布荷载产生的内力函数分解为由未知参数表示的线形项，和与分布荷载特性相关 的已知函数项的组合，计算积分算子，由积分算子的线形组合得到分布荷载产生的杆件内力函数，从而得到 荷载矩阵。该方法没有采用任何假定，具有广泛的普适性，能用于求解任意刚度分布杆件上作用任意分布荷 载或集中荷载的单元荷载矩阵。对于刚度非连续、荷载分布非连续的情况，该方法避免了高斯数值积分的误 差，具有更高的计算精度， 关键词：非线性杆单元，荷载矩阵 中图分类号：TU313 文献标志码：A [2~4] 对于材料非线性直杆单元，当前集成刚度矩阵和荷载矩阵多采用高斯积分法 。而高 斯积分得到精确解的前提是各参数表达式是连续多项式，且多项式的次数与高斯积分点数对 应。在大多数应用中，这一条件并不能得到满足。如当前大量采用的用于计算钢筋混凝土梁 柱单元的非线性纤维单元，其截面参数本身即为一些离散点，并且在混凝土开裂后局部刚度 突变，如仍采用高斯积分法，其误差非常大。对于纤维单元，采用高斯积分法还有一个限制， 即要求纤维分段必须是等间隔的，分段数不能太多，不能进行两端加密、跨中加密等处理， 限制了纤维单元的应用。 相对于刚度矩阵，荷载矩阵集成的问题更多。对于局部分布荷载、分段分布荷载、集中 力荷载，采用高斯积分即存在误差。如果同时单元刚度分布是离散、突变的，高斯积分的误 差更大。 文献[1]在一般积分定积分表示方法的基础上，提出基于积分算子直接集成非线性、变截 面杆件单元刚度矩阵的算法。本文在该方法的基础上，将积分算子直接集成方法扩展到非线 性杆单元荷载矩阵的计算。 本文方法没有采用任何假定，具有广泛的普适性。不但适用于非线性纤维单元，也适用 于多段变截面单元、刚臂单元，加入适当的积分边界条件，即可用于任意位置塑性铰单元。 对于荷载，可以精确导算任意类型分布荷载、集中荷载，包括抛物线预应力索荷载。 在使用本文方法计算荷载矩阵之前，假定非线性杆单元的截面参数EI(x) 、EA(x) 、GA(x) 和GJ(x)等均已知，同时杆上局部分布荷载p(x)或集中荷载P 的描叙函数已知。 1. 横向分布荷载的导算 杆件单元中间分布荷载在有限元计算中，需要在假定两端嵌固的条件下，求出杆端作用 力。嵌固杆件是超静定的结构，杆端力与刚度分布EI(x)有关。对于图(1a)所示嵌固杆件，可 以等价地分解为作用中间作用分布荷载铰支的静定杆( 图1b)，和两端作用端弯矩的的铰支静 1 定杆( 图1c)两个部分叠加之和。 p(x) p(x) Y Y Y M0 ML