[基础科学]平面桁架的几何非线性有限元分析.pdfVIP

第2 1 卷 第4 期 长 沙 交 通 学 院 学 报 Vol.2 1 No.4 2005 年 12 月 JOURNAL OF CHANGSHA COMMUNICATIONS UNIVERSITY Dec. 2005 : 1000- 77 (2005) 04- 003 - 03 1 2 , ( 1. , 4 10076; 2. , 4 10081) : 平面杆系结构的几何非线性分析方法有多种 以平面桁架为例, 从杆的平衡方程出 发, 直接利用变分方法导出杆单元在大转动、小应变条件下的考虑几何非线性的单元切线刚 度矩阵, 依照此原理编制了相应的有限元程序 算例表明, 推导结果是正确的, 该程序可用于 平面桁架的几何非线性分析 : 平面桁架; 几何非线性; 有限元; 变分方法 : T U311. 1 : A 杆系结构的几何非线性分析一直是结构分析中的一个研究热点, 指导研究理论的不同也就会得出 [1, 2] 不尽相同的结果 就目前而言, 对杆系结构进行考虑几何非线性的结构分析一般有两类分析方法 : 一类是非线性有限元方法, 该方法是通过选择适当的位移模式, 结合非线性的应变- 位移关系, 运用虚 功原理或虚位移原理得出杆端力与杆端位移之间的平衡关系; 另一类是稳定函数法, 此种方法是根据梁 杆的微分方程, 通过引入稳定函数系数来考虑轴力与弯矩的耦合作用及其对结构单元刚度的影响 计 算表明, 上述两种方法均具有较好的精度 本文以平面杆单元为例, 从杆的平衡方程出发, 直接利用变 分方法导出杆单元在大转动、小应变条件下的考虑几何非线性的单元切线刚度矩阵 1 图1 所示为一平面桁架杆元ij , oxy 为结 构坐标系, 杆元初始 t 时刻构形的单元节点坐 t t t t t T 标向量{X}= ( x , y , x , y ) ; 在荷载作用下i i j j 单元 t + $ t 时刻构形的节点坐标向量为 t+ $ t t+ $ t t+ $ t t+ $ t t+ $ t T { X}= ( x ,i y ,i x ,j y ) ; 节点j e T 的位移向量为{U} = ( u , v , u , v ) , 根据几i i j j 何关系可得出: t+ $ t t e { X}= {X}+ {U} ( 1) 设在两种构形中桁架杆单元长度分别设 t t+ $ t 为l 和 l,有 t t 2 t 2 l = l + l x y , t+ $ t t+ $ t 2 t+ $ t 2 (2) 1 l = ( l ) + ( l ) x y 其中, l = x - x , l = y - y ,t t t+ $ tl = l + u - u ,t t+ $ tl = l + v - vt 则变形前后桁架杆单元的轴 x j i y j