一种新的集成的非线性杆件单元刚度矩阵的方法-佳构STRAT-通用.PDF

该文发表于：《计算力学学报》，2003 年8 月，20 卷第4 期，490~494 一种新的集成的非线性杆件单元刚度矩阵的方法 谢靖中，李国强 ( 同济大学建筑工程系，上海，200092) 摘 要：对于非线性杆件单元，本文提出一种新的简便有效的集成单元刚度矩阵的算法。该方法直接从结构 力学中的位移法的概念出发，通过解析积分或数值积分求解积分算子，由积分算子线性组合，能快速求解考 虑弯、剪、扭、轴压等各种非线形刚度的杆件单元的刚度矩阵。该方法具有广泛的普适性，能适用于所有多 项式、插值多项式、解析式、离散点描述的变刚度、变截面直杆的单元刚度矩阵集成计算。文中通过求解线 性直杆单元刚度矩阵验证了该方法的正确性。 关键词：非线形杆件单元，变刚度杆单元，刚度矩阵，积分算子 中图分类号：TU313 文献标志码：A 1 概 述 EI(x) 、EA(x) 、GA(x)和GJ(x)等参数已知。 对于非线形杆件单元，在杆件分段截面刚度参 2 单元刚度公式推导 数 EI(x) 、GA(x)、GJ(x)、EA(x)等确定之后，由这 些截面刚度集成单元刚度矩阵时，现有的计算方法 2.1 抗弯、抗剪刚度 大多采用设定的位移插值函数和力场插值函数，由 对于图1 所示两端嵌固变刚度杆单元，抗弯模 虚位移原理通过高斯数值积分，得出单元刚度矩阵 量EI(x) ，剪切模量EA(x) ，没有杆中分布荷载作用。 [1~3] 。 当杆端存在横向位移 或转角 时，则按照力平衡 L L 对于力场插值函数，根据静力平衡条件，一般 的原则，杆中剪力是一恒定值，杆中分布弯矩呈线 能符合实际情况。但对于位移插值函数，由于非线 形变化，假定弯矩M(x)、剪力Q(x)为: 形杆变形受刚度分布影响，很难准确假定，因此有 M(x)=Ax+B Q(x)=A (1) 很大的随意性，这直接影响了计算的精度和可靠 性。 L Y EI(x) 本文直接从结构力学位移法基本概念出发，通 y L 过引入重积分的定积分表示方法，实现了集成刚度 X L 的一般广义积分在杆长度范围内的定积分表示。进 ML 一步通过几个规格化的定积分算子的求解，以及这 些积分算子的线形组合，得出非线形杆件单元的弯 剪刚度矩阵、轴向刚度矩阵、扭转刚度矩阵。这种 M 0 直接积分求解刚度矩阵的方法，没有采用任何假 Q 定，这里称之为积分算子直接刚度法。 本文不涉及非线形杆件的材料本构关系假定， 以及截面刚度参数的算法。本文集中论述由单元截