薄壁杆件截面几何特性计算比拟有限元法.pdfVIP

第19卷第1期 工 程 力 学 V01．19No．1 2002年2月 同婀D旺匹RDfG^怔’cHANICS Feb．2002 文章编号：1000-4750(2002)01-03：8-04 薄壁杆件截面几何特性计算的比拟有限元法 周建春1，刘光栋2，魏琴1 《I．华南理工大学．广末广州510640；2．湖南大学．湖南长沙410082) 擒蔓：针对任意截面形式的开闭口薄壁截面几何特性的计算，提出了～种便于计算机实现的比拟有限元法。 根据能量驻值原理导出的截面单元扇性坐标控制方程，通过采用有限元比拟方法，建立了求解截面扇牲坐标等 几何特性的计算机算法，并用FORTRAN语言编制了相应的计算机程序CPTBGS。本文所提出算法的优点是物 理概念明确且效率报高．特别适用于具有任意截面形式的开闭口薄壁杆件截面几何特性的计算机计算。 关键词：任意截面薄壁杆件：截面几何特性；算法及程序实现：比拟有限元法 中围分类号：TU323文献标识码：A 1 引言 点为i和j，沿杆件轴线方向的长度为单位长度 沿单元壁厚中心线的局部坐标为J。 在薄壁杆件的约束扭转分析理论中，计算截面 杠～‘ 的扇性坐标等一些截面几何性质是最重要和最基 本的计算内容，一般常用薄壁杆件的约束扭转理论 来求解．对于复杂截面形式．如广泛应用于桥粱工 产氐 程中的各种闭口和开闭口组合截面的多室薄壁箱 粱，其计算工作量是相当大的。为此，本文提出了 D 一种便于计算机实现的比拟有限元法。根据能量驻 图1单元￡的扇性坐标 值原理导出的截面单元扇性坐标控制方程，通过采 Fig．1Sectionalcoordinateofelementt 用有限元比拟方法，建立了求解截面扇性坐标等几 设截面的扇性极点为o，则单元￡中任一点j 何特性的计算机算法，并用FORTRAN语言编制了 的扇性坐标为 相应的计算机程序CPTBGS。本文提出的算法的突 ∞(j)=∞。+)(国J一∞．)÷ (I) 出优点是物理概念明确且效率很高．特别适用于具 ‘^ 有任意截面形式的开闭口薄壁杆件截面几何特性 其中∞，和∞，分别为单元t的节点i和J处的扇性坐 的计算机计算。 标．f。为单元长度，s为局部坐标。 2截面扇性坐标控制方程的建立【1l 利用自由扭转时薄壁杆件的翘曲剪应变r。的 表达式，并考虑式(1)有 在计算中假定截面可按一定长度分为一系列 k k：娑+娄：呻，粤+州 等厚度的直线微段，一个直线微段即称为一个单 2面+嚣2呻。百“妒。 一 元，因此每个单元上的扇性坐标都是线性分布的； ：妒’(^一掌)：妒，(^^一等兰) 口J 单