几何材料非线性的新梁柱单元及程序编制.pdf

2010正 5 月 郑 州 大 学 学 报 (工 学 版 ) May 2010 第 31卷 第 3期 JournalofZhengzhouUniversity(EngineeringScience) Vol_31 No．3 文章编号 ：1671—6833(2010)03—0065—05 几何材料非线性的新梁柱单元及程序编制 张俊峰 ，郝际平 ，李 天 (1．郑州大学 土木工程学院，河南 郑州 450001；2．西安建筑科技大学 土木工程学院，陕西 西安 710055) 摘 要 ：给 出了基于uL法的三维粱柱单元虚功增量方程 ，详细推导 了考虑翘曲和剪切变形影响的三维 空间粱柱单元的几何非线性切线刚度矩阵，同时对 Orbison截面塑性面进行 了修正，以考虑扭矩和翘 曲 对截面强度 的影响，采用塑性流动法则推导了单元弹塑性切线刚度矩阵．根据面 向对象的程序设计思 想 ，将整个有 限元域划分为 8个基本类 ，在单元基类的基础上派生 了新 的单元类 ，采用 c++语 言编制 了 面向对象的空间钢框架分析程序．算例证明，只需要一到两个单元就可 以准确预测空间钢框架的极限承 栽力和失稳模 态 ． 关键词 ：空 间钢框架 ；翘 曲；剪切变形 ；塑性铰 ；极 限承载力；面向对象程序设计 中图分类号 ：TU313．3 文献标识码 ：A 横向位移和转角采用三次插值 函数 ，考虑了剪切 0 引言 变形和翘曲的影响，推导了空间薄壁梁柱单元 的 结构的几何材料非线性分析一直是结构分析 几何非线性刚度矩阵 ，在此基础上采用修正 的 研究 中的热点 ，在三维钢结构的非线性分析中，通 Orbison屈服面方程 得到了单元 的弹塑性刚度 常采用两种方法建立梁柱单元 的刚度方程 ：梁柱 矩阵． 法和有限单元法．基于梁柱理论的梁柱法 ，可以很 1 空间非线性梁柱单元刚度矩阵 好地模拟平面结构的受力性能 ，但对于空间结构 ， 采用 由平面直接扩展到三维的梁柱单元，不能有 采用的基本假定为构件是等截面的，并且双 效地预测弯扭屈曲．采用有限单元法 ，基于虚功原 轴对称 ；材料是理想弹塑性 ；构件的塑性 出现在杆 理的平衡方程可 以考虑轴向变形与弯曲、轴向变 端 ；大位移小应变 ；不考虑局部屈曲的影响． 形与剪切 、轴 向变形与扭转、双向弯 曲、弯曲与扭 1．1 虚功增量方程 转等耦合作用，能够很好 的模拟空间结构的受力 如图 1所示 ，在几何非线性增量分析中，采用 性能…．在位移插值 函数的选取上 ，有三次插值 三种基本 的构形来描述运动的物体 ：以c。表示初 多项式 I3、五次多项式 以及基于稳定 函数 的 始未变形 的构形 ；C，表示 己知的变形后的某一参 插值多项式 等．基于稳定函数的插值多项式 考构形 ；C：表示变形后 的构形 ，则 Lagrangian构 建立在梁柱方程基础上，可以很好 的反映结构 的 形的增量虚功方程式为 ： 受力变形 ，但 由于拉压不同状态下位移插值 函数 J(Ee…6e +4Ge6e +4Ge6e；)dI，+ 不同，使用起来很不方便 ，LIEW 等 对稳定插值 函数和三次插值多项式 的对 比认为：当轴 向力不 f( 6 +21 叩+21：艿：)d 大时，两种方法的最大误差不超过 5％，这一误差 是工程中所允许的．TEH 也对三次插值函数做 = R一尺