论文部分内容阅读
本文主要研究了框架几何非线性分析的若干问题及一些问题解决办法。平面梁单元修正拉格朗日列式分析中,以截面形心位移函数为参数,推导得出了截面转角与形心线横向位移函数和轴向位移函数的准确关系。考虑二阶影响后,得出截面转角不仅与截面形心横线位移导数有关,还与纵向位移导数有关。推导得出了修正拉格朗日列式几何非线性分析的切线刚度矩阵。分析表明,梁单元几何刚度矩阵能否通过刚体检验与截面转角与横向和纵向位移导数之间的关系近似程度有关,通常分析中,仅考虑截面转角与横向位移导数的线性关系,因而得出的几何刚度矩阵不能通过刚体检验。分析了采用同一切线刚度矩阵计算增量节点位移和增量节点力问题,研究表明,采用同一切线刚度矩阵计算增量节点力,与已有节点力的和是欲求构形节点力,这个力是用已知构形坐标描述的,必须采用由增量节点位移确定的近似坐标变换,把这个节点力变换到欲求构形,才得欲求构形单元真实节点力;采用同一刚度矩阵计算单元增量节点力,数值计算结果表明几何刚度矩阵能否通过刚体检验对分析结果并无明显影响。平面梁单元分析中,推导了考虑加载变形弯曲影响的修正拉格朗日列式,计算增量节点力采用了协同转动列式,分析表明对于某些情形考虑加载弯曲变形,会有效减少单元划分数量;但对于某些情形,考虑加载弯曲会导致膜锁。以悬臂曲梁为例,分析了产生膜锁的原因,得出忽略初始弯曲变形二阶项以及高阶项是产生膜锁的主要原因,这些项虽然很小,但系数很大,其乘积可能比有关的线性项还大,从而影响了分析结果,导致误差甚至错误。推导分析了绕定轴旋转、绕从轴旋转、半切向旋转的数学表达式,得出绕定轴转角、绕从轴转角与转角顺序有关;由空间梁单元变形后单元两端的旋转矩阵,精确得出了单元定位矩阵。分析了由于截面转动导致的节点弯矩增量和节点弯矩虚功,得出保守力矩的虚功与转角类型和力矩类型有关。以半切向转角为变量,采用虚功原理推导得出了空间梁单元修正拉格朗日列式的切线刚度矩阵;计算增量节点力采用了协同转动列式,精确考虑增量节点位移确定的刚体旋转,导致计算增量节点力的刚度矩阵是不对称矩阵;若仅考虑增量节点位移二阶项确定的刚体旋转,得到的刚度矩阵为对称矩阵。对于保守弯矩,节点外荷载弯矩的虚功产生外荷载刚度矩阵,这个矩阵会改变几何刚度矩阵一些项。全面搜集整理了空间梁单元几何非线性分析算例,众多的数值计算结果表明,对于空间梁单元采用修正拉格朗日列式计算单元增量节点位移,用协同转动列式计算单元的增量节点力是可行、有效的。以杆件非线性理论为基础,推导得出了杆件稳定分析的平衡方程。以悬臂梁在悬臂端作用弯矩的临界弯矩分析为例,推导得出了半切向弯矩和两类准切向弯矩的临界弯矩;运用有限元软件采用实体单元,分析了悬臂梁在悬臂端作用弯矩的临界荷载,理论分析结果和有限元数值计算结果吻合很好。