梁系结构二阶分析的单元模型.pdfVIP

第42卷第2期 西 安 建 筑 科技 大 学 学 报(自然科学版) V01．42’No．2 2010年4月 Arch．Tech．(NaturalScience Univ．of Edition) Apr．2010 J．Xi’an 梁系结构二阶分析的单元模型 齐 麟 (天津大学建筑工程学院，天津300072) 摘要：以梁单元平衡微分方程为基础，推导了考虑剪切变形影响的空间梁单元位移插值函数，用于模拟结构 中无轴力的杆件．利用Maclaurin级数统一了空间梁柱单元在拉、压两种情况下的不同位移插值函数，与稳定 函数表示的位移插值函数等价，用于模拟结构中有轴力的杆件．推导了包含轴向变形、剪切变形、双向弯曲、扭 转及其耦合效应的二阶单元切线刚度矩阵．从计算精度与总体刚度矩阵的正定性两方面考虑确定了位移插 值函数中级数的展开项数．对承受不同轴压力的悬臂梁梁端位移进行计算，表明所提出的单元模型可以很好 地体现二阶效应的影响．利用不同单元模型对单层柱面网壳进行对比分析，表明所提出的梁系结构二阶分析 单元模型具有较高的计算精度与效率，可以很好地反映单层柱面网壳的几何非线性． 关键词：梁系结构；二阶；位移插值函数；级数；单元切线刚度矩阵 中图分类号：TU393．3 文献标识码：A 文章编号：1006-7930(2010)02-0191-05 梁系结构杆件同时承受横向荷载与轴向荷载，杆件截面会产生附加弯矩，结构分析时应以变形后的 几何关系建立平衡方程，进行二阶分析．二阶分析普遍采用Oran梁柱单元模型[11和空间梁元模型心]． Oran梁柱单元位移插值函数对高阶项无省略，但没有考虑剪切变形影响，将平面单元扩展为空间单元 后某些耦合项被遗失．空间梁元模型忽略了应变函数中位移的高阶项，使计算精度受到影响．文献E33～ Es3提出的SIF(稳定位移插值函数)单元考虑了剪切变形的影响，对高阶项无省略，但受拉与受压采用 不同的表达式，且无法模拟无轴力的杆件．文献E63采用Taylor级数简化稳定函数，得到的位移插值函 数可以用于轴力为零的情况，但由于采用了一些近似替代，轴力接近零时二者计算结果有明显差别．针 对上述问题本文推导了考虑剪切变形影响的空间梁单元位移插值函数，用于模拟无轴力杆件；统一梁柱 单元的位移插值函数，与稳定位移插值函数等价，用于模拟有轴力杆件；推导了考虑轴向变形、剪切变 形、双向弯曲、扭转及其耦合效应的空间梁柱单元二阶切线刚度矩阵，对梁系结构进行二阶分析． 1 空间梁单元的位移插值函数 平面梁单元受力变形如图1所示． 令弯矩与剪力引起的平动位移分别为YM与YQ，由yllM一一 收稿日期：2009—11—13修改稿日期：2010—02一01 基金项目：国家自然科学基金重大研究计划重点支持项目 作者简介：齐麟(1982一)，男，山西太原人，博士研究生，主要从事空间结构建模理论的研究 192 西安建筑科技大学学报(自然科学版) 第42卷 式中：点w=1／(t3／12EI+∥／GA)，其他参数的物理意义见图1． 距离坐标原点为z的任意截面上的剪力Q(z)为 Q(x7)=一Q=一kN『I专(只+Ok)一(文一向)] (3) 根据Timoshenko理论，单元内任意截面转角口(z)为 『茜㈨2—2z)+钟t (4) 推广到空间的情况，空间梁单元任意截面形心处的位移增量△110与杆端位移增量之间△H的关系 为： Nl 0 0 0 0 0 N2 0 0 0 0 0 N。 0 N3y 0 0 0 N姆0Ns，00 0 Nh N0